ModuleDistortionCalculatorProcessor.h

/*
 * ModuleDistortionCalculatorProcessor.h
 *
 *  Created on: June 6, 2013
 *      Author: amunnich
 */

#ifndef MODULEDISTORTIONCALCULATORPROSECCOR_H
#define MODULEDISTORTIONCALCULATORPROSECCOR_H

#include <marlin/Processor.h>
#include <lcio.h>

#include "ModuleDistortionCalculator.h"


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

    virtual Processor*  newProcessor() { return new ModuleDistortionCalculatorProcessor; }
    
    ModuleDistortionCalculatorProcessor();
    
    virtual void init();
    
    virtual void processRunHeader(lcio::LCRunHeader* run);
    
    virtual void processEvent(lcio::LCEvent* evt);
    
    virtual void check(lcio::LCEvent* evt);
    
    virtual void end();
    
    void WriteDistortionFiles(std::string filename,std::map< std::pair<int,int>, std::pair<double,double> > results);
    void WriteRootFile();
    
  protected:
    /* the place for protected and private member data and functions */
    std::string _inputColNameTrackHits; 
    std::string _inputColNameTrackFit; 

  private:
    int _usedEvents;
    int _fillCount;
    std::string _distortionXYOutputFileName;
    std::string _distortionRowOutputFileName;
    std::string _distortionZOutputFileName;
    std::string _rootOutputFileName;
    
    ModuleDistortionCalculator * _myModuleDistortionCalculator;
    
    std::map< std::pair<int,int>, TH1*> _HistogramsXY;
    std::map< std::pair<int,int>, TH1*> _HistogramsRow;
    std::map< std::pair<int,int>, TH1*> _HistogramsZ;
 
    std::map< std::pair<int,int>, std::pair<double,double> >  _ResultsXY;
    std::map< std::pair<int,int>, std::pair<double,double> >  _ResultsRow;
    std::map< std::pair<int,int>, std::pair<double,double> >  _ResultsZ;
 
    const gear::TPCParameters* _GearParameters;
  };
}

#endif //MODULEDISTORTIONCALCULATORPROSECCOR_H