DreamSD.h

Go to the documentation of this file.

#ifndef SimG4CMS_DreamSD_h
#define SimG4CMS_DreamSD_h

#include "SimG4CMS/Calo/interface/CaloSD.h"
#include "DetectorDescription/Core/interface/DDCompactView.h"
#include "DetectorDescription/DDCMS/interface/DDCompactView.h"

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ParameterSetfwd.h"

#include "G4PhysicsFreeVector.hh"

#include <DD4hep/DD4hepUnits.h>

#include <map>

const int MAXPHOTONS = 500;  // Maximum number of photons we can store

class G4LogicalVolume;

class DreamSD : public CaloSD {
public:
  DreamSD(const std::string &,
          const DDCompactView *,
          const cms::DDCompactView *,
          const SensitiveDetectorCatalog &,
          edm::ParameterSet const &,
          const SimTrackManager *);
  ~DreamSD() override {}

  uint32_t setDetUnitId(const G4Step *) override;

protected:
  double getEnergyDeposit(const G4Step *) override;
  void initRun() override;

private:
  typedef std::pair<double, double> Doubles;
  typedef std::map<G4LogicalVolume *, Doubles> DimensionMap;

  void initMap(const std::string &);
  void fillMap(const std::string &, double, double);
  double curve_LY(const G4Step *, int);
  double crystalLength(G4LogicalVolume *) const;
  double crystalWidth(G4LogicalVolume *) const;

  double cherenkovDeposit_(const G4Step *aStep);
  double getAverageNumberOfPhotons_(const double charge,
                                    const double beta,
                                    const G4Material *aMaterial,
                                    const G4MaterialPropertyVector *rIndex);
  double getPhotonEnergyDeposit_(const G4ParticleMomentum &p, const G4ThreeVector &x, const G4Step *aStep);
  bool setPbWO2MaterialProperties_(G4Material *aMaterial);

  static constexpr double k_ScaleFromDDDToG4 = 1.0;
  static constexpr double k_ScaleFromDD4hepToG4 = 1.0 / dd4hep::mm;

  const DDCompactView *cpvDDD_;
  const cms::DDCompactView *cpvDD4hep_;

  bool useBirk_, doCherenkov_, readBothSide_, dd4hep_;
  double birk1_, birk2_, birk3_;
  double slopeLY_;
  DimensionMap xtalLMap_;  // Store length and width

  int side_;

  std::unique_ptr<G4PhysicsFreeVector> chAngleIntegrals_;
  G4MaterialPropertiesTable *materialPropertiesTable_;

  int nphotons_;
};

#endif  // DreamSD_h