LDCCaloDigi.h

#ifndef DIGITIZER_LDCCALODIGI_H
#define DIGITIZER_LDCCALODIGI_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 LDCCaloDigi : public Processor {
  
 public:
  
  LDCCaloDigi(const LDCCaloDigi&) = delete;
  LDCCaloDigi& operator=(const LDCCaloDigi&) = delete;

  virtual Processor*  newProcessor() { return new LDCCaloDigi ; }
  
  
  LDCCaloDigi() ;
  
  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 _outputEcalCollection{};
  std::string _outputHcalCollection{};
  std::string _outputRelCollection{};

  float _thresholdEcal{};
  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{};

  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