CMS 3D CMS Logo

List of all members | Public Member Functions | Protected Member Functions | Private Member Functions | Private Attributes
EcalTBH4BeamSD Class Reference

#include <EcalTBH4BeamSD.h>

Inheritance diagram for EcalTBH4BeamSD:
CaloSD SensitiveCaloDetector Observer< const BeginOfRun *> Observer< const BeginOfEvent *> Observer< const BeginOfTrack *> Observer< const EndOfTrack *> Observer< const EndOfEvent *> SensitiveDetector

Public Member Functions

 EcalTBH4BeamSD (const std::string &, const SensitiveDetectorCatalog &, edm::ParameterSet const &, const SimTrackManager *)
 
uint32_t setDetUnitId (const G4Step *step) override
 
void setNumberingScheme (EcalNumberingScheme *scheme)
 
 ~EcalTBH4BeamSD () override
 
- Public Member Functions inherited from CaloSD
 CaloSD (const std::string &aSDname, const SensitiveDetectorCatalog &clg, edm::ParameterSet const &p, const SimTrackManager *, float timeSlice=1., bool ignoreTkID=false, const std::string &newcolname="")
 
void clear () override
 
void clearHits () override
 
void DrawAll () override
 
void EndOfEvent (G4HCofThisEvent *eventHC) override
 
void fillHits (edm::PCaloHitContainer &, const std::string &) override
 
void Initialize (G4HCofThisEvent *HCE) override
 
bool isItFineCalo (const G4VTouchable *touch)
 
void newCollection (const std::string &name, edm::ParameterSet const &p)
 
void PrintAll () override
 
G4bool ProcessHits (G4Step *step, G4TouchableHistory *) override
 
bool ProcessHits (G4GFlashSpot *aSpot, G4TouchableHistory *) override
 
void reset () override
 
 ~CaloSD () override
 
- Public Member Functions inherited from SensitiveCaloDetector
 SensitiveCaloDetector (const std::string &iname, const SensitiveDetectorCatalog &clg, const std::string &newcollname="")
 
- Public Member Functions inherited from SensitiveDetector
void EndOfEvent (G4HCofThisEvent *eventHC) override
 
const std::vector< std::string > & getNames () const
 
void Initialize (G4HCofThisEvent *eventHC) override
 
bool isCaloSD () const
 
 SensitiveDetector (const std::string &iname, const SensitiveDetectorCatalog &, bool calo, const std::string &newcollname="")
 
 ~SensitiveDetector () override
 
- Public Member Functions inherited from Observer< const BeginOfRun *>
 Observer ()
 
void slotForUpdate (const BeginOfRun * iT)
 
virtual ~Observer ()
 
- Public Member Functions inherited from Observer< const BeginOfEvent *>
 Observer ()
 
void slotForUpdate (const BeginOfEvent * iT)
 
virtual ~Observer ()
 
- Public Member Functions inherited from Observer< const BeginOfTrack *>
 Observer ()
 
void slotForUpdate (const BeginOfTrack * iT)
 
virtual ~Observer ()
 
- Public Member Functions inherited from Observer< const EndOfTrack *>
 Observer ()
 
void slotForUpdate (const EndOfTrack * iT)
 
virtual ~Observer ()
 
- Public Member Functions inherited from Observer< const EndOfEvent *>
 Observer ()
 
void slotForUpdate (const EndOfEvent * iT)
 
virtual ~Observer ()
 

Protected Member Functions

double getEnergyDeposit (const G4Step *) override
 
- Protected Member Functions inherited from CaloSD
bool checkHit (int k=0)
 
CaloG4HitcreateNewHit (const G4Step *, const G4Track *, int k)
 
virtual void endEvent ()
 
virtual double EnergyCorrected (const G4Step &step, const G4Track *)
 
virtual bool filterHit (CaloG4Hit *, double)
 
unsigned int findBoundaryCrossingParent (const G4Track *track, bool markParentAsSaveable=true)
 
double getAttenuation (const G4Step *aStep, double birk1, double birk2, double birk3) const
 
virtual uint16_t getDepth (const G4Step *)
 
virtual bool getFromLibrary (const G4Step *step)
 
int getNumberOfHits (int k=0)
 
double getResponseWt (const G4Track *, int k=0)
 
virtual int getTrackID (const G4Track *)
 
bool hitExists (const G4Step *, int k)
 
void ignoreRejection ()
 
virtual void initEvent (const BeginOfEvent *)
 
virtual void initRun ()
 
void printDetectorLevels (const G4VTouchable *) const
 
void processHit (const G4Step *step)
 
virtual void processSecondHit (const G4Step *, const G4Track *)
 
void resetForNewPrimary (const G4Step *)
 
void setNumberCheckedHits (int val, int k=0)
 
void setParameterized (bool val)
 
G4ThreeVector setToGlobal (const G4ThreeVector &, const G4VTouchable *) const
 
G4ThreeVector setToLocal (const G4ThreeVector &, const G4VTouchable *) const
 
virtual int setTrackID (const G4Step *)
 
void setUseMap (bool val)
 
std::string shortreprID (const CaloHitID &ID)
 
std::string shortreprID (const CaloG4Hit *hit)
 
void update (const BeginOfRun *) override
 This routine will be called when the appropriate signal arrives. More...
 
void update (const BeginOfEvent *) override
 This routine will be called when the appropriate signal arrives. More...
 
void update (const BeginOfTrack *trk) override
 This routine will be called when the appropriate signal arrives. More...
 
void update (const EndOfTrack *trk) override
 This routine will be called when the appropriate signal arrives. More...
 
void update (const ::EndOfEvent *) override
 
void updateHit (CaloG4Hit *, int k)
 
- Protected Member Functions inherited from SensitiveDetector
TrackInformationcmsTrackInformation (const G4Track *aTrack)
 
Local3DPoint ConvertToLocal3DPoint (const G4ThreeVector &point) const
 
Local3DPoint FinalStepPosition (const G4Step *step, coordinates) const
 
Local3DPoint InitialStepPosition (const G4Step *step, coordinates) const
 
Local3DPoint LocalPostStepPosition (const G4Step *step) const
 
Local3DPoint LocalPreStepPosition (const G4Step *step) const
 
void NaNTrap (const G4Step *step) const
 
void setNames (const std::vector< std::string > &)
 
- Protected Member Functions inherited from Observer< const EndOfEvent *>
virtual void update (const EndOfEvent *)=0
 This routine will be called when the appropriate signal arrives. More...
 

Private Member Functions

void getBaseNumber (const G4Step *aStep)
 

Private Attributes

double birk1
 
double birk2
 
double birk3
 
EcalNumberingSchemenumberingScheme
 
EcalBaseNumber theBaseNumber
 
bool useBirk
 
bool useWeight
 

Additional Inherited Members

- Protected Types inherited from SensitiveDetector
enum  coordinates { WorldCoordinates, LocalCoordinates }
 
- Static Protected Member Functions inherited from CaloSD
static std::string printableDecayChain (const std::vector< unsigned int > &decayChain)
 
- Protected Attributes inherited from CaloSD
std::string collName_ [2]
 
CaloG4HitcurrentHit [2]
 
CaloHitID currentID [2]
 
std::string detName_
 
float edepositEM
 
float edepositHAD
 
double eminHit
 
double energyCut
 
G4ThreeVector entranceLocal
 
G4ThreeVector entrancePoint
 
bool forceSave
 
std::vector< std::string > hcn_
 
float incidentEnergy
 
double kmaxIon
 
double kmaxNeutron
 
double kmaxProton
 
int nHC_
 
G4ThreeVector posGlobal
 
CaloHitID previousID [2]
 
bool suppressHeavy
 
double tmaxHit
 
std::vector< int > useResMap_
 

Detailed Description

Definition at line 21 of file EcalTBH4BeamSD.h.

Constructor & Destructor Documentation

◆ EcalTBH4BeamSD()

EcalTBH4BeamSD::EcalTBH4BeamSD ( const std::string &  name,
const SensitiveDetectorCatalog clg,
edm::ParameterSet const &  p,
const SimTrackManager manager 
)

Definition at line 22 of file EcalTBH4BeamSD.cc.

References birk1, birk2, birk3, g, edm::ParameterSet::getParameter(), mergeVDriftHistosByStation::name, AlCaHLTBitMon_ParallelJobs::p, generator_cfi::scheme, setNumberingScheme(), and useBirk.

26  : CaloSD(name, clg, p, manager), numberingScheme(nullptr) {
27  edm::ParameterSet m_EcalTBH4BeamSD = p.getParameter<edm::ParameterSet>("EcalTBH4BeamSD");
28  useBirk = m_EcalTBH4BeamSD.getParameter<bool>("UseBirkLaw");
29  birk1 = m_EcalTBH4BeamSD.getParameter<double>("BirkC1") * (CLHEP::g / (CLHEP::MeV * CLHEP::cm2));
30  birk2 = m_EcalTBH4BeamSD.getParameter<double>("BirkC2");
31  birk3 = m_EcalTBH4BeamSD.getParameter<double>("BirkC3");
32 
33  EcalNumberingScheme *scheme = nullptr;
34  if (name == "EcalTBH4BeamHits") {
36  } else {
37  edm::LogWarning("EcalTBSim") << "EcalTBH4BeamSD: ReadoutName not supported\n";
38  }
39 
40  if (scheme)
42  edm::LogInfo("EcalTBSim") << "Constructing a EcalTBH4BeamSD with name " << GetName();
43  edm::LogInfo("EcalTBSim") << "EcalTBH4BeamSD:: Use of Birks law is set to " << useBirk
44  << " with three constants kB = " << birk1 << ", C1 = " << birk2
45  << ", C2 = " << birk3;
46 }
T getParameter(std::string const &) const
Definition: ParameterSet.h:307
The Signals That Services Can Subscribe To This is based on ActivityRegistry and is current per Services can connect to the signals distributed by the ActivityRegistry in order to monitor the activity of the application Each possible callback has some defined which we here list in angle e g
Definition: Activities.doc:4
void setNumberingScheme(EcalNumberingScheme *scheme)
EcalNumberingScheme * numberingScheme
Log< level::Info, false > LogInfo
CaloSD(const std::string &aSDname, const SensitiveDetectorCatalog &clg, edm::ParameterSet const &p, const SimTrackManager *, float timeSlice=1., bool ignoreTkID=false, const std::string &newcolname="")
Definition: CaloSD.cc:34
Log< level::Warning, false > LogWarning

◆ ~EcalTBH4BeamSD()

EcalTBH4BeamSD::~EcalTBH4BeamSD ( )
override

Definition at line 48 of file EcalTBH4BeamSD.cc.

References numberingScheme.

48  {
49  if (numberingScheme)
50  delete numberingScheme;
51 }
EcalNumberingScheme * numberingScheme

Member Function Documentation

◆ getBaseNumber()

void EcalTBH4BeamSD::getBaseNumber ( const G4Step *  aStep)
private

Definition at line 79 of file EcalTBH4BeamSD.cc.

References EcalBaseNumber::addLevel(), EcalBaseNumber::getCapacity(), cuy::ii, LogDebug, EcalBaseNumber::reset(), EcalBaseNumber::setSize(), and theBaseNumber.

Referenced by setDetUnitId().

79  {
81  const G4VTouchable *touch = aStep->GetPreStepPoint()->GetTouchable();
82  int theSize = touch->GetHistoryDepth() + 1;
83  if (theBaseNumber.getCapacity() < theSize)
84  theBaseNumber.setSize(theSize);
85  // Get name and copy numbers
86  if (theSize > 1) {
87  for (int ii = 0; ii < theSize; ii++) {
88  theBaseNumber.addLevel(touch->GetVolume(ii)->GetName(), touch->GetReplicaNumber(ii));
89  LogDebug("EcalTBSim") << "EcalTBH4BeamSD::getBaseNumber(): Adding level " << ii << ": "
90  << touch->GetVolume(ii)->GetName() << "[" << touch->GetReplicaNumber(ii) << "]";
91  }
92  }
93 }
void addLevel(const std::string &name, const int &copyNumber)
ii
Definition: cuy.py:589
EcalBaseNumber theBaseNumber
void setSize(const int &size)
#define LogDebug(id)

◆ getEnergyDeposit()

double EcalTBH4BeamSD::getEnergyDeposit ( const G4Step *  aStep)
overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 53 of file EcalTBH4BeamSD.cc.

References birk1, birk2, birk3, CaloSD::getAttenuation(), LogDebug, useBirk, and mps_merge::weight.

53  {
54  // take into account light collection curve for crystals
55  double weight = 1.;
56  if (useBirk)
57  weight *= getAttenuation(aStep, birk1, birk2, birk3);
58  double edep = aStep->GetTotalEnergyDeposit() * weight;
59  LogDebug("EcalTBSim") << "EcalTBH4BeamSD:: " << aStep->GetPreStepPoint()->GetPhysicalVolume()->GetName()
60  << " Light Collection Efficiency " << weight << " Weighted Energy Deposit " << edep / CLHEP::MeV
61  << " MeV";
62  return edep;
63 }
Definition: weight.py:1
double getAttenuation(const G4Step *aStep, double birk1, double birk2, double birk3) const
Definition: CaloSD.cc:728
#define LogDebug(id)

◆ setDetUnitId()

uint32_t EcalTBH4BeamSD::setDetUnitId ( const G4Step *  step)
overridevirtual

Implements CaloSD.

Definition at line 65 of file EcalTBH4BeamSD.cc.

References getBaseNumber(), EcalNumberingScheme::getUnitID(), numberingScheme, and theBaseNumber.

65  {
66  getBaseNumber(aStep);
67  return (numberingScheme == nullptr ? 0 : numberingScheme->getUnitID(theBaseNumber));
68 }
EcalNumberingScheme * numberingScheme
EcalBaseNumber theBaseNumber
void getBaseNumber(const G4Step *aStep)
virtual uint32_t getUnitID(const EcalBaseNumber &baseNumber) const =0

◆ setNumberingScheme()

void EcalTBH4BeamSD::setNumberingScheme ( EcalNumberingScheme scheme)

Definition at line 70 of file EcalTBH4BeamSD.cc.

References numberingScheme, and generator_cfi::scheme.

Referenced by EcalTBH4BeamSD().

70  {
71  if (scheme != nullptr) {
72  edm::LogInfo("EcalTBSim") << "EcalTBH4BeamSD: updates numbering scheme for " << GetName() << "\n";
73  if (numberingScheme)
74  delete numberingScheme;
76  }
77 }
EcalNumberingScheme * numberingScheme
Log< level::Info, false > LogInfo

Member Data Documentation

◆ birk1

double EcalTBH4BeamSD::birk1
private

Definition at line 39 of file EcalTBH4BeamSD.h.

Referenced by EcalTBH4BeamSD(), and getEnergyDeposit().

◆ birk2

double EcalTBH4BeamSD::birk2
private

Definition at line 39 of file EcalTBH4BeamSD.h.

Referenced by EcalTBH4BeamSD(), and getEnergyDeposit().

◆ birk3

double EcalTBH4BeamSD::birk3
private

Definition at line 39 of file EcalTBH4BeamSD.h.

Referenced by EcalTBH4BeamSD(), and getEnergyDeposit().

◆ numberingScheme

EcalNumberingScheme* EcalTBH4BeamSD::numberingScheme
private

Definition at line 36 of file EcalTBH4BeamSD.h.

Referenced by setDetUnitId(), setNumberingScheme(), and ~EcalTBH4BeamSD().

◆ theBaseNumber

EcalBaseNumber EcalTBH4BeamSD::theBaseNumber
private

Definition at line 40 of file EcalTBH4BeamSD.h.

Referenced by getBaseNumber(), and setDetUnitId().

◆ useBirk

bool EcalTBH4BeamSD::useBirk
private

Definition at line 38 of file EcalTBH4BeamSD.h.

Referenced by EcalTBH4BeamSD(), and getEnergyDeposit().

◆ useWeight

bool EcalTBH4BeamSD::useWeight
private

Definition at line 37 of file EcalTBH4BeamSD.h.