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:
  
  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