LEPTrackingProcessor.h

#ifndef LEPTrackingProcessor_h
#define LEPTrackingProcessor_h 1

#include "marlin/Processor.h"
#include "lcio.h"
#include <string>
#include <EVENT/TrackerHit.h>
#include <IMPL/LCCollectionVec.h>
#include <EVENT/SimTrackerHit.h>
#include "TFile.h"
#include "TH1.h"
#include "TH2.h"


using namespace lcio ;
using namespace marlin ;



class LEPTrackingProcessor : public Processor {
  
 public:
  
  virtual Processor*  newProcessor() { return new LEPTrackingProcessor ; }
  
  
  LEPTrackingProcessor() ;
  
  virtual void init() ;
  
  virtual void processRunHeader( LCRunHeader* run ) ;
  
  virtual void processEvent( LCEvent * evt ) ; 
  
  
  virtual void check( LCEvent * evt ) ; 
  
  
  virtual void end() ;
  
  void selectTPCHits(LCCollection* tpcTHcol, LCCollectionVec* remainingCol, LCCollectionVec* cutCol, int binHeight, int binWidth, float zmin=-10000., float zmax=+10000.);

  void selectTPCHits(LCCollection* tpcTHcol, LCCollection* usedCol, float zmin = -10000., float zmax = +10000.);

  void FillTPCHitBanks();

 protected:

  std::string _colNameTPC ;
  std::string _colNameVTX ;
  std::string _colNameSIT ;
  std::string _colNameTPCTracks ;
  std::string _colNameTracks ;
  std::string _colNameMCTPCTracksRel ;
  std::string _colNameMCTracksRel ;
  std::string _usedColName;
  std::string _droppedColName;

  LCCollectionVec* _tpcTrackVec ;
  LCCollectionVec* _tpclcRelVec ;
  LCCollectionVec* _usedCol ;
  LCCollectionVec* _droppedCol ;
  LCCollectionVec* _usedColForOutput ;
  LCCollectionVec* _droppedColForOutput ;

  int _nRun ;
  int _nEvt ;

  int _binHeight ;
  int _nSlicesInZ ;
  int _binWidth ;
  int _multiplicityCut ;
  int _AlwaysRunCurlKiller ;
  std::vector<EVENT::TrackerHit*> _goodHits;
  std::vector<EVENT::TrackerHit*> _savedGoodHits;
  int   _histograms;
  TH1F* fTPCR;
  TH1F* fTPCRRaw;
  TH1F* fTPCZ;
  TH1F* fTPCZRaw;
  TH2F* fTPCRZ;
  TH2F* fTPCXYRaw;
  TH2F* fTPCXY;
  TH2F* fTPCRZRaw;


} ;

#endif