MokkaCaloDigi.h

#ifndef MokkaCaloDigi1_h
#define MokkaCaloDigi1_h 1

#include "marlin/Processor.h"
#include "EVENT/SimCalorimeterHit.h"
#include "IMPL/CalorimeterHitImpl.h"
#include "IMPL/LCCollectionVec.h"
#include "lcio.h"
#include <string>


using namespace lcio ;
using namespace marlin ;

struct MyHit {
  CalorimeterHitImpl * hit{};
  std::vector<SimCalorimeterHit*> simHits{};
};


class MokkaCaloDigi : public Processor {
  
 public:
  
  MokkaCaloDigi(const MokkaCaloDigi&) = delete;
  MokkaCaloDigi& operator=(const MokkaCaloDigi&) = delete;

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


  MyHit * ProcessHitInBarrel(SimCalorimeterHit * hit);
  MyHit * ProcessHitInEndcap(SimCalorimeterHit * hit);

 protected:

  std::vector<std::string>  _ecalCollections{};
  std::vector<std::string>  _hcalCollections{};
  std::string _newCollNameHCAL{};
  std::string _newCollNameECAL{};
  std::string _relationCollName{};
  int _nRun{};
  int _nEvt{};

  float _thresholdEcal{};
  float _thresholdHcal{};

  int _digitalEcal{};
  int _digitalHcal{};

  std::vector<int> _ecalLayers{};
  std::vector<int> _hcalLayers{};
  std::vector<float> _calibrCoeffEcal{};
  std::vector<float> _calibrCoeffHcal{};
  float * _endBarrelChamberLength{};
  float * _barrelLateralWidth{};
  float * _barrelOffsetMaxX{};
  float * _endBarrelOffsetMaxZ{};
  float _regularBarrelOffsetMaxZ{};
  float _lateralPlateThickness{};
  float _modulesGap{};
  float _innerHcalRadius{};
  int _numberOfHcalLayers{};
  int _nStaves{};
  int _nModules{};
  int _cellScaleX{};
  int _cellScaleZ{};
  float _newCellSizeX{};
  float _newCellSizeZ{};
  float _hcalLayerThickness{};
  float _hcalAbsorberThickness{};
  float _hcalSensitiveThickness{};
  float _virtualCellSizeX{};
  float _virtualCellSizeZ{};
  float _regularBarrelModuleLength{};
  float _regularBarrelChamberLength{};
  float _deltaPhi{};
  std::vector< std::vector<MyHit*> > _calorimeterHitVec{};
  LCCollectionVec * _relationCollection{};
  int _startIEndcap{};
  int _startJEndcap{};
  float _startXEndcap{};
  float _startZEndcap{};



} ;

#endif