AnalysisCutSelectorProcessor.h

#ifndef ANALYSISCUTSELECTORPROCESSOR_H
#define ANALYSISCUTSELECTORPROCESSOR_H

// C++
#include <string>

// LCIO
#include "lcio.h"
#include <EVENT/TrackerHit.h>
#include <EVENT/LCCollection.h>
#include <EVENT/Track.h>
#include <IMPL/LCCollectionVec.h>
#include <IMPL/TrackImpl.h>
#include <IMPL/LCFlagImpl.h>

// Marlin
#include "marlin/Processor.h"
#include "gear/TPCParameters.h"
//ROOT
#include "TH1F.h"
#include "TNtuple.h"

// Aida
#ifdef MARLIN_USE_AIDA
//AIDA
#include <marlin/AIDAProcessor.h>
#endif

namespace marlintpc{

    class AnalysisCutSelectorProcessor : public marlin::Processor {

    public:

        /* This method will be called by the marlin package
         * It returns a processor of the current type (here AnalysisCutSelectorProcessor)
         */
        virtual Processor*  newProcessor() { return new AnalysisCutSelectorProcessor ; }

        AnalysisCutSelectorProcessor() ;

        virtual void init() ;

        virtual void processRunHeader( lcio::LCRunHeader* run ) ;

        virtual void processEvent( lcio::LCEvent * evt ) ;

        virtual void end();

        protected:

        std::string _inputColName ;   
        std::string _outputColName ;  

        int _outputIsTransient;

        // event specific cuts
        int _lowcutNTracks;   
        int _highcutNTracks;  

        // track specific cuts
        int _lowcutNHitsOnTrack;   
        int _highcutNHitsOnTrack;  
        float _lowcutPhi;          
        float _highcutPhi;         
        float _lowcutChiSq;          
        float _highcutChiSq;         
        float _lowcutTanLambda;    
        float _highcutTanLambda;   
        float _lowcutdEdx;         
        float _highcutdEdx;        
        float _lowcutOmega;        
        float _highcutOmega;       
        float _lowcutD0;           
        float _highcutD0;          
        float _lowcutZ0;           
        float _highcutZ0;          
        double _runZOverride;       
        double _runZErrorOverride;  
        bool  _cutIfFitFailed;     
        int _bestModule; 
        int _bestRow; 

        float _zPosHistogramMax; 

    private:
        TH1F * _cutFlowEvents;
        TH1F * _cutFlowTracks;
        TH1F * _hitZPosition;
        TH1F * _hitTime;
        TNtuple* _runZPosTuple;
        TNtuple* _runHitTimeTuple;
        // angle between the hit and the track
        TH1F * _hitPhi;
        TNtuple* _runHitPhiTuple;
        TH1F * _hitPhiBestRow;
        TNtuple* _runHitPhiBestRowTuple;

    const gear::TPCParameters* _gearParameters;

    } ;


}

#endif // ANALYSISCUTSELECTORPROCESSOR_H