CutApplicationProcessor.h

#ifndef CutApplicationProcessor_h
#define CutApplicationProcessor_h 1

// C++
#include <string>
#include <iostream>
#include <cmath>
#include <algorithm>

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

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

// Marlin
#include "marlin/Processor.h"

namespace AIDA{
  class IHistogram1D;
  class IHistogram2D;
}

namespace marlintpc{

    class CutApplicationProcessor : public marlin::Processor {
        
    public:
        
        /* This method will be called by the marlin package
         * It returns a processor of the currend type (here CutApplicationProcessor)
         */
        virtual Processor*  newProcessor() { return new CutApplicationProcessor ; }
        
        CutApplicationProcessor() ;
        
        virtual void init() ;
        
        virtual void processRunHeader( lcio::LCRunHeader* run ) ;
        
        virtual void processEvent( lcio::LCEvent * evt ) ; 
        
        static bool compareHit(const std::pair<TrackerHit*, double> & pair1,
                               const std::pair<TrackerHit*, double> & pair2)
          {
            return pair1.second < pair2.second;
          }
        
    protected:
        
        std::string _inputColName ;   
        std::string _outputColName ;  
        double _chipHeight;           
        double _chipLength;           

        int _outputIsTransient;

        // event specific cuts
        int _lowcutNTracks;   
        int _highcutNTracks;  

        // track specific cuts
        int _lowcutNHitsOnTrack;   
        int _highcutNHitsOnTrack;  
        int _lowcutNHitsPerRowOnTrack;   
        int _highcutNHitsPerRowOnTrack;  
        float _lowcutPhi;          
        float _highcutPhi;         
        float _lowcutTanLambda;    
        float _highcutTanLambda;   
        float _lowcutdEdx;         
        float _highcutdEdx;        
        float _lowcutOmega;        
        float _highcutOmega;       
        float _lowcutD0;           
        float _highcutD0;          
        float _lowcutZ0;           
        float _highcutZ0;          
        bool  _cutIfFitFailed;     

        // hit specific cuts
        float _lowcutCharge;       
        float _highcutCharge;      
        double _lowcutX;           
        double _highcutX;          
        double _lowcutY;           
        double _highcutY;          
        double _lowcutZ;           
        double _highcutZ;          
        double _lowCutTrackLength; 
        double _highCutTrackLength;
        double _lowCutXTrack;      
        double _highCutXTrack;     
        double _lowCutYTrack;      
        double _highCutYTrack;     
//private:

    } ;
    
}

#endif