DreamSD.h

Go to the documentation of this file.

#ifndef SimG4CMS_DreamSD_h
#define SimG4CMS_DreamSD_h

#include "SimG4CMS/Calo/interface/CaloSD.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "G4String.hh"
#include "G4PhysicsOrderedFreeVector.hh"

#include <map>

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

class G4LogicalVolume;

class TTree;

class DreamSD : public CaloSD {

public:    

  DreamSD(const std::string&, const DDCompactView &, const SensitiveDetectorCatalog &,
      edm::ParameterSet const &, const SimTrackManager*);
  ~DreamSD() override {}
  //  bool   ProcessHits(G4Step * step,G4TouchableHistory * tHistory) override;
  uint32_t setDetUnitId(const G4Step*) override;

protected:

  // G4bool getStepInfo(G4Step* aStep) override;
  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&, const DDCompactView &);
  double         curve_LY(const G4Step*, int); 
  const double   crystalLength(G4LogicalVolume*) const;
  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 );

  bool         useBirk, doCherenkov_, readBothSide_;
  double       birk1, birk2, birk3;
  double       slopeLY;
  DimensionMap xtalLMap; // Store length and width

  int          side;

  std::unique_ptr<G4PhysicsOrderedFreeVector> chAngleIntegrals_;
  G4MaterialPropertiesTable* materialPropertiesTable;
  // Histogramming
  TTree* ntuple_;
  int nphotons_;
  float px_[MAXPHOTONS],py_[MAXPHOTONS],pz_[MAXPHOTONS];
  float x_[MAXPHOTONS],y_[MAXPHOTONS],z_[MAXPHOTONS];

};

#endif // DreamSD_h