ThrustReconstruction.h

#ifndef ThrustReconstruction_h
#define ThrustReconstruction_h 1
#include <vector>
#include "marlin/Processor.h"
#include "lcio.h"
#include <iostream>
#include <string>
#include <IMPL/ReconstructedParticleImpl.h>

#include <CLHEP/Vector/ThreeVector.h>
#include <CLHEP/Random/RanluxEngine.h>

namespace CLHEP{}    // declare namespace CLHEP for backward compatibility
using namespace CLHEP ;


using namespace lcio ;
using namespace marlin ;



class ThrustReconstruction : public Processor {
  
public:
  
  ThrustReconstruction(const ThrustReconstruction&) = delete;
  ThrustReconstruction& operator=(const ThrustReconstruction&) = delete;

  virtual Processor* newProcessor() { return new ThrustReconstruction;}
  
  
  ThrustReconstruction() ;
  
  virtual void init() ;
  
  virtual void processRunHeader( LCRunHeader* run ) ;
  virtual void modifyRunHeader( LCRunHeader* /*run*/ ) {}
  
  virtual void processEvent( LCEvent * evt ) ; 
  virtual void modifyEvent( LCEvent * /*evt*/ ) {}
  
 
  
  
  virtual void end() ;
  
protected:
  int TassoThrust();
  int JetsetThrust();
  double sign(double a,double b);
  double min(double a,double b);

  std::string  _inputCollectionName{};
  int _typeOfThrustFinder{};

  float _principleThrustValue{};
  float _majorThrustValue{};
  float _minorThrustValue{};
  Hep3Vector _principleThrustAxis{};
  Hep3Vector _majorThrustAxis{};
  Hep3Vector _minorThrustAxis{};
  float _min{},_max{};
  LCCollection* _inParVec{};
  std::vector<Hep3Vector> _partMom{};
  std::string filename{};
  RanluxEngine myrnd{};
} ;

#endif