RPCutProcessor.h

#ifndef RPCutProcessor_h
#define RPCutProcessor_h 1

#include <EVENT/ReconstructedParticle.h>
#include <UTIL/LCRelationNavigator.h>

#include <gear/VXDParameters.h>
#include <gear/GearParameters.h>
#ifdef MCFAIL_DIAGNOSTICS
#include "TH1.h"
#include "TH2.h"
#include "TFile.h"
#endif

#include "marlin/Processor.h"
#include "lcio.h"
#include <string>
using namespace lcio ;
using namespace marlin ;



class RPCutProcessor : public Processor {
  
 public:
  //The usual Marlin processor methods
  virtual Processor*  newProcessor() { return new RPCutProcessor ; }
  RPCutProcessor() ;
  virtual void init() ;
  virtual void processRunHeader( LCRunHeader* run ) ;
  virtual void processEvent( LCEvent * evt ) ; 
  virtual void check( LCEvent * evt ) ; 
  virtual void end() ;
  protected:

  bool  _Chi2OverDOFFail(ReconstructedParticle* RPTrack);
  bool  _D0Fail(ReconstructedParticle* RPTrack);
  bool  _D0ErrFail(ReconstructedParticle* RPTrack);
  bool  _Z0Fail(ReconstructedParticle* RPTrack);
  bool  _Z0ErrFail(ReconstructedParticle* RPTrack);
  bool  _PTFail(ReconstructedParticle* RPTrack);
  bool  _DetectorHitsFail(ReconstructedParticle* RPTrack,std::map<std::string,int> SubdetectorIndex);
  bool  _MCPIDFail( lcio::ReconstructedParticle* RPTrack, UTIL::LCRelationNavigator* pMCRelationNavigator );
  bool  _BadParametersFail(lcio::ReconstructedParticle* RPTrack);
  bool  _MCVertexFail(lcio::ReconstructedParticle* RPTrack, UTIL::LCRelationNavigator* pMCRelationNavigator );
  
  std::string _InRCPColName ;
  std::string _TrackColName;
  std::string _OutRCPColName ;
  bool _WriteNewCollection;
  bool _SubParticleLists;
  
  bool _Chi2OverDOFEnable;
  bool _Chi2OverDOFCutLowerThan;
  float _Chi2OverDOFCutValue;
  
  bool _D0Enable;
  bool _D0CutLowerThan;
  float _D0CutValue;
  
  bool _D0ErrEnable;
  bool _D0ErrCutLowerThan;
  float _D0ErrCutValue;
  
  bool _Z0Enable;
  bool _Z0CutLowerThan;
  float _Z0CutValue;

  bool _Z0ErrEnable;
  bool _Z0ErrCutLowerThan;
  float _Z0ErrCutValue;
  
  bool _PTEnable;
  bool _PTCutLowerThan;
  float _PTCutValue;

  bool _MonteCarloPIDEnable;
  std::vector<int> _MonteCarloPIDsToCut;
  std::string _MonteCarloRelationColName;
  
  bool _BadParametersEnable;
  
  bool _DetectorHitsEnable;
  std::vector<std::string> _DetectorHitsBoundaryDetectorNames;
  std::vector<int> _DetectorHitsBoundaryCuts;
  std::vector<std::string> _DetectorHitsRegion1DetectorNames;
  std::vector<int> _DetectorHitsRegion1Cuts;
  std::vector<std::string> _DetectorHitsRegion2DetectorNames;
  std::vector<int> _DetectorHitsRegion2Cuts;
  std::vector<std::string> _DetectorNames;
  
  bool _MCVertexEnable;
  const gear::VXDParameters* _VxdPar;
  double _BeamPipeInnerRadius;
  double _BeamPipeOuterRadius;
  double _BeamPipeHalfZ;
  double _CutDistance;

#ifdef MCFAIL_DIAGNOSTICS
  TH2F *_diaghist_bpmat_xy;
  TH2F *_diaghist_vxmat_xy;
  TH2F *_diaghist_nomat_xy;
  TH2F *_diaghist_bpmat_rz;
  TH2F *_diaghist_vxmat_rz;
  TH2F *_diaghist_nomat_rz;
  TH1F *_diaghist_dist;
  TH1F *_diaghist_dist_vxmat;
  TH1F *_diaghist_dist_nomat;
#endif
  
  int _nRun ;
  int _nEvt ;
} ;

#endif