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