NewLDCCaloDigi.h

#ifndef DIGITIZER_NEWLDCCALODIGI_H
#define DIGITIZER_NEWLDCCALODIGI_H 1

#include "marlin/Processor.h"
#include <IMPL/CalorimeterHitImpl.h>
#include "lcio.h"
#include <string>
#include <vector>

using namespace lcio ;
using namespace marlin ;

const int MAX_LAYERS = 200;
const int MAX_STAVES =  16;

class NewLDCCaloDigi : public Processor {
  
 public:
  
  virtual Processor*  newProcessor() { return new NewLDCCaloDigi ; }
  
  
  NewLDCCaloDigi() ;
  
  virtual void init() ;
  
  virtual void processRunHeader( LCRunHeader* run ) ;
  
  virtual void processEvent( LCEvent * evt ) ; 
   
  virtual void check( LCEvent * evt ) ; 
  
  virtual void end() ;

  virtual void fillECALGaps() ;
  
  
 protected:

  int _nRun ;
  int _nEvt ;
  
  std::vector<std::string> _ecalCollections;
  std::vector<std::string> _hcalCollections;


  std::string _outputEcalCollection0;
  std::string _outputEcalCollection1;
  std::string _outputEcalCollection2;
  std::string _outputHcalCollection0;
  std::string _outputHcalCollection1;
  std::string _outputHcalCollection2;
  std::vector<std::string> _outputEcalCollections;
  std::vector<std::string> _outputHcalCollections;
  std::string _outputRelCollection;

  float _thresholdEcal;
  std::vector<float> _thresholdHcal;

  int _digitalEcal;
  int _digitalHcal;

  std::vector<float> _calibrCoeffEcal;
  std::vector<float> _calibrCoeffHcal;

  std::vector<int> _ecalLayers;
  std::vector<int> _hcalLayers;

  int _ecalGapCorrection;
  float _ecalGapCorrectionFactor;
  float _ecalModuleGapCorrectionFactor;
  float _ecalEndcapCorrectionFactor;
  float _hcalEndcapCorrectionFactor;


  std::string _cellIDLayerString ;
  std::string _cellIDModuleString ;
  std::string _cellIDStaveString ;


  std::vector<CalorimeterHitImpl*> _calHitsByStaveLayer[MAX_STAVES][MAX_LAYERS];
  std::vector<int> _calHitsByStaveLayerModule[MAX_STAVES][MAX_LAYERS];

  float _zOfEcalEndcap;
  float _barrelPixelSizeT[MAX_LAYERS];
  float _barrelPixelSizeZ[MAX_LAYERS];
  float _endcapPixelSizeX[MAX_LAYERS];
  float _endcapPixelSizeY[MAX_LAYERS];
  float _barrelStaveDir[MAX_STAVES][2];


} ;

#endif