ECalSD.h

Go to the documentation of this file.

#ifndef SimG4CMS_ECalSD_h
#define SimG4CMS_ECalSD_h
// File: ECalSD.h
// Description: Stores hits of Electromagnetic calorimeters in appropriate
//              container
// Use in your sensitive detector builder:
//    ECalSD* ecalSD = new ECalSD(SDname, new CaloNumberingScheme());
//#define plotDebug

#include "SimG4CMS/Calo/interface/CaloSD.h"
#include "SimG4CMS/Calo/interface/EnergyResolutionVsLumi.h"
#include "Geometry/EcalCommonData/interface/EcalNumberingScheme.h"
#include "DetectorDescription/Core/interface/DDsvalues.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "G4String.hh"
#ifdef plotDebug
#include <TH2F.h>
#endif
#include <string>
#include <map>

class EcalBaseNumber;
class G4LogicalVolume;
class EnergyResolutionVsLumi;

class ECalSD : public CaloSD {
public:
  ECalSD(const std::string &,
         const DDCompactView &,
         const SensitiveDetectorCatalog &,
         edm::ParameterSet const &p,
         const SimTrackManager *);
  ~ECalSD() override;
  uint32_t setDetUnitId(const G4Step *) override;
  void setNumberingScheme(EcalNumberingScheme *);

protected:
  double getEnergyDeposit(const G4Step *) override;
  int getTrackID(const G4Track *) override;
  uint16_t getDepth(const G4Step *) override;

private:
  void initMap(const G4String &, const DDCompactView &);
  uint16_t getRadiationLength(const G4StepPoint *hitPoint, const G4LogicalVolume *lv);
  uint16_t getLayerIDForTimeSim();
  double curve_LY(const G4LogicalVolume *);

  void getBaseNumber(const G4Step *);
  double getBirkL3(const G4Step *);

  std::vector<double> getDDDArray(const std::string &, const DDsvalues_type &);
  std::vector<std::string> getStringArray(const std::string &, const DDsvalues_type &);

  // initialised before run
  bool isEB;
  bool isEE;
  EcalNumberingScheme *numberingScheme;
  bool useWeight, storeTrack, storeRL, storeLayerTimeSim;
  bool useBirk, useBirkL3;
  double birk1, birk2, birk3, birkSlope, birkCut;
  double slopeLY, scaleRL;
  std::string crystalMat, depth1Name, depth2Name;
  std::map<const G4LogicalVolume *, double> xtalLMap;
  std::vector<const G4LogicalVolume *> useDepth1, useDepth2, noWeight;
  EcalBaseNumber theBaseNumber;
  EnergyResolutionVsLumi ageing;
  bool ageingWithSlopeLY;

  // run time cache
  G4ThreeVector currentLocalPoint;
  double crystalLength;
  double crystalDepth;
  uint16_t depth;

#ifdef plotDebug
  TH2F *g2L_[4];
#endif
};

#endif  // ECalSD_h