InterTrackDistanceProcessor.h

/*
 * InterTrackDistanceProcessor.h
 *
 *  Created on: October 2nd, 2013
 *      Author: amunnich
 */

#ifndef INTERTRACKDISTANCEPROCESSOR_H
#define INTERTRACKDISTANCEPROCESSOR_H

#include <marlin/Processor.h>
#include <lcio.h>
#include "EVENT/Track.h"

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

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

#include "TH1F.h"
#include "TH2.h"
#include "TNtuple.h"


namespace marlintpc
{
  class InterTrackDistanceProcessor : public marlin::Processor
  {
  public:

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

    InterTrackDistanceProcessor();

    virtual void init();

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

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

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

    virtual void end();

    std::vector< std::pair<int,double> > Get_RPhiDistVec(){return _RPhiDistVec;};
    std::vector< std::pair<int,double> > Get_ZDistVec(){return _ZDistVec;};
    void SetInputColName(std::string name){_inputColName=name;};
    std::pair<double,double> Get_PhiIntersection(lcio::Track *track, int module, int row){return phiIntersection(track, module, row);};

  protected:
    std::string _inputColName; 


  private:
    void DoubleTrackProductionEfficiency(std::vector<lcio::Track *> trackVec);
    std::pair<double,double> phiIntersection(lcio::Track *track, int module, int row);
    double zIntersection(lcio::Track *track, int module, int row);
        const gear::TPCParameters* _GearParameters;

        TH2F * _histoInterTrackrphiDistance;
        TH2F * _histoInterTrackZDistance;
        TH2F * _histoInterTrackDistance;
        TH1F * _histoDoubleTrackEfficiency;

        int    _NumRowsWithMergedTracks;
        float _MaxDistBetweenMergedTracks;
        int    _NumRowsWithSeparatedTracks;
        float _MinDistBetweenSeparatedTracks;

        float _Max_Z_DistBetweenMergedTracks;
        float _Min_Z_DistBetweenSeparatedTracks;

        int _eff_counter;
        int _all_counter;

    float _histogramMinDistance; 
    float _histogramMaxDistance; 
    IntVec _moduleList; 

    std::vector< std::pair<int,double> > _RPhiDistVec;
    std::vector< std::pair<int,double> > _ZDistVec;

  };
}

#endif //INTERTRACKDISTANCEPROCESSOR_H