ResolutionPerformanceProcessor.h

#ifndef RESOLUTIONPERFORMANCEPROCESSOR_H
#define RESOLUTIONPERFORMANCEPROCESSOR_H

#include <marlin/Processor.h>
#include "marlin/Global.h"

#include <lcio.h>
#include <string>
#include "TrackFitterBase.h"
#include <EVENT/TrackerHit.h>

//MarlinTPC
#include "ModuleDistortionCalculator.h"

//GEAR
#include "gear/TPCParameters.h"

//ROOT
#include "TH1.h"
#include "TH2.h"
#include "TNtuple.h"


// In case you want AIDA histograms

namespace marlintpc
{

class ResolutionPerformanceProcessor : public marlin::Processor
{
public:

        virtual Processor*  newProcessor() { return new ResolutionPerformanceProcessor; }

        ResolutionPerformanceProcessor();

        virtual void init();

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

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

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

        virtual void end();



protected:
        /* the place for protected and private member data and functions */
        std::string _inputColName; 

private:
        double calcHitPhi(TrackerHit * firstHit, TrackerHit * secondHit);

        double _bfieldScaleFactor; 
        int _xyResolutionBins; 
        float _xyResolutionRange; 
        int _zResolutionBins; 
        float _zResolutionRange; 
        int _chargeBins; 
        int _chargeUpperBound; 
        bool _useHitsWithOverflow; 

        std::map< std::pair<int,int>, TH1*> _zDistanceRow;
        std::map< std::pair<int,int>, TH1*> _xyDistanceRow;
        std::map< std::pair<int,int>, TH1*> _zResidualRow;
        std::map< std::pair<int,int>, TH1*> _xyResidualRow;
        std::map< std::pair<int,int>, TH1*> _phiRow;

        std::map< std::pair<int,int>, TH2*> _zResidualVSChargeRow;

        TH1 * _zDistanceAll;
        TH1 * _xyDistanceAll;
        TH1 * _zResidualAll;
        TH1 * _xyResidualAll;

        TH2 * _zResidualVSChargeAllRows;
        TH2 * _xyResidualVSChargeAllRows;

        TNtuple * _allRowsTuple;
        TNtuple * _singleRowsTuple;
        TNtuple * _bestRowTuple;

        TNtuple * _zResidualVSChargeAllRowsTuple;
        TNtuple * _zResidualVSChargeSingleRowsTuple;

        TrackFitterBase * _trackFitter;
        int _meanMethodZ; 
        int _meanMethodXY; 
        int _bestModule; 
        int _bestRow; 

        ModuleDistortionCalculator * _myModuleDistortionCalculator;
        const gear::TPCParameters* _gearParameters;

        float _maxChargeCut;
        float _minChargeCut;
        int _minBinsRiseTimeCut;
};
} // namespace marlintpc
#endif // RESOLUTIONPERFORMANCEPROCESSOR_H