CMS 3D CMS Logo

List of all members | Public Member Functions | Static Public Member Functions | Protected Member Functions | Private Member Functions | Private Attributes
HGCalTB16SD01 Class Reference
Inheritance diagram for HGCalTB16SD01:
CaloSD SensitiveCaloDetector Observer< const BeginOfRun * > Observer< const BeginOfEvent * > Observer< const BeginOfTrack * > Observer< const EndOfTrack * > Observer< const EndOfEvent * > SensitiveDetector

Public Member Functions

 HGCalTB16SD01 (const std::string &, const edm::EventSetup &, const SensitiveDetectorCatalog &, edm::ParameterSet const &, const SimTrackManager *)
 
uint32_t setDetUnitId (const G4Step *step) override
 
 ~HGCalTB16SD01 () override=default
 
- Public Member Functions inherited from CaloSD
 CaloSD (const std::string &aSDname, const edm::EventSetup &es, const SensitiveDetectorCatalog &clg, edm::ParameterSet const &p, const SimTrackManager *, float timeSlice=1., bool ignoreTkID=false)
 
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 PrintAll () override
 
bool ProcessHits (G4GFlashSpot *aSpot, G4TouchableHistory *) override
 
G4bool ProcessHits (G4Step *step, G4TouchableHistory *) override
 
void reset () override
 
 ~CaloSD () override
 
- Public Member Functions inherited from SensitiveCaloDetector
 SensitiveCaloDetector (const std::string &iname, const edm::EventSetup &es, const SensitiveDetectorCatalog &clg, edm::ParameterSet const &p)
 
- 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 edm::EventSetup &es, const SensitiveDetectorCatalog &, edm::ParameterSet const &p, bool calo)
 
 ~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 ()
 

Static Public Member Functions

static uint32_t packIndex (int det, int lay, int x, int y)
 
static void unpackIndex (const uint32_t &idx, int &det, int &lay, int &x, int &y)
 

Protected Member Functions

double getEnergyDeposit (const G4Step *) override
 
- Protected Member Functions inherited from CaloSD
bool checkHit ()
 
CaloG4HitcreateNewHit (const G4Step *, const G4Track *)
 
virtual void endEvent ()
 
virtual double EnergyCorrected (const G4Step &step, const G4Track *)
 
virtual bool filterHit (CaloG4Hit *, double)
 
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 ()
 
double getResponseWt (const G4Track *)
 
virtual int getTrackID (const G4Track *)
 
bool hitExists (const G4Step *)
 
void ignoreRejection ()
 
virtual void initEvent (const BeginOfEvent *)
 
virtual void initRun ()
 
void processHit (const G4Step *step)
 
void resetForNewPrimary (const G4Step *)
 
void setNumberCheckedHits (int val)
 
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)
 
void update (const ::EndOfEvent *) override
 
void update (const BeginOfEvent *) override
 This routine will be called when the appropriate signal arrives. More...
 
void update (const BeginOfRun *) 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 updateHit (CaloG4Hit *)
 
- 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 initialize (const G4StepPoint *point)
 

Private Attributes

double birk1_
 
double birk2_
 
double birk3_
 
bool initialize_
 
std::string matName_
 
G4Material * matScin_
 
bool useBirk_
 

Additional Inherited Members

- Protected Types inherited from SensitiveDetector
enum  coordinates { WorldCoordinates, LocalCoordinates }
 
- Protected Attributes inherited from CaloSD
CaloG4HitcurrentHit
 
CaloHitID currentID
 
float edepositEM
 
float edepositHAD
 
double eminHit
 
double energyCut
 
G4ThreeVector entranceLocal
 
G4ThreeVector entrancePoint
 
bool forceSave
 
float incidentEnergy
 
double kmaxIon
 
double kmaxNeutron
 
double kmaxProton
 
G4ThreeVector posGlobal
 
CaloHitID previousID
 
bool suppressHeavy
 
double tmaxHit
 

Detailed Description

Definition at line 15 of file HGCalTB16SD01.cc.

Constructor & Destructor Documentation

◆ HGCalTB16SD01()

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

Definition at line 40 of file HGCalTB16SD01.cc.

45  : CaloSD(name, es, clg, p, manager), initialize_(true) {
46  // Values from NIM 80 (1970) 239-244: as implemented in Geant3
47  edm::ParameterSet m_HC = p.getParameter<edm::ParameterSet>("HGCalTestBeamSD");
48  matName_ = m_HC.getParameter<std::string>("Material");
49  useBirk_ = m_HC.getParameter<bool>("UseBirkLaw");
50  birk1_ = m_HC.getParameter<double>("BirkC1") * (g / (MeV * cm2));
51  birk2_ = m_HC.getParameter<double>("BirkC2");
52  birk3_ = m_HC.getParameter<double>("BirkC3");
53  matScin_ = nullptr;
54 
55  edm::LogVerbatim("HGCSim") << "HGCalTB16SD01:: Use of Birks law is set to " << useBirk_ << " for " << matName_
56  << " with three constants kB = " << birk1_ << ", C1 = " << birk2_ << ", C2 = " << birk3_;
57 }

References birk1_, birk2_, birk3_, g, edm::ParameterSet::getParameter(), matName_, matScin_, MeV, AlCaHLTBitMon_ParallelJobs::p, AlCaHLTBitMon_QueryRunRegistry::string, and useBirk_.

◆ ~HGCalTB16SD01()

HGCalTB16SD01::~HGCalTB16SD01 ( )
overridedefault

Member Function Documentation

◆ getEnergyDeposit()

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

Reimplemented from CaloSD.

Definition at line 59 of file HGCalTB16SD01.cc.

59  {
60  auto const point = aStep->GetPreStepPoint();
61  if (initialize_)
63  double destep = aStep->GetTotalEnergyDeposit();
64  double wt2 = aStep->GetTrack()->GetWeight();
65  double weight = (wt2 > 0.0) ? wt2 : 1.0;
66  if (useBirk_ && matScin_ == point->GetMaterial()) {
68  }
69 #ifdef EDM_ML_DEBUG
70  edm::LogVerbatim("HGCSim") << "HGCalTB16SD01: Detector " << point->GetTouchable()->GetVolume()->GetName() << " with "
71  << point->GetMaterial()->GetName() << " weight " << weight << ":" << wt2;
72 #endif
73  return weight * destep;
74 }

References birk1_, birk2_, birk3_, CaloSD::getAttenuation(), initialize(), initialize_, matScin_, point, and useBirk_.

◆ initialize()

void HGCalTB16SD01::initialize ( const G4StepPoint *  point)
private

Definition at line 122 of file HGCalTB16SD01.cc.

122  {
123  if (matName_ == point->GetMaterial()->GetName()) {
124  matScin_ = point->GetMaterial();
125  initialize_ = false;
126  }
127 #ifdef EDM_ML_DEBUG
128  edm::LogVerbatim("HGCSim") << "HGCalTB16SD01: Material pointer for " << matName_
129  << " is initialized to : " << matScin_;
130 #endif
131 }

References initialize_, matName_, matScin_, and point.

Referenced by getEnergyDeposit().

◆ packIndex()

uint32_t HGCalTB16SD01::packIndex ( int  det,
int  lay,
int  x,
int  y 
)
static

Definition at line 86 of file HGCalTB16SD01.cc.

86  {
87  int ix = 0, ixx = x;
88  if (x < 0) {
89  ix = 1;
90  ixx = -x;
91  }
92  int iy = 0, iyy = y;
93  if (y < 0) {
94  iy = 1;
95  iyy = -y;
96  }
97  uint32_t idx = (det & 15) << 28; // bits 28-31
98  idx += (lay & 127) << 21; // bits 21-27
99  idx += (iy & 1) << 19; // bit 19
100  idx += (iyy & 511) << 10; // bits 10-18
101  idx += (ix & 1) << 9; // bit 9
102  idx += (ixx & 511); // bits 0-8
103 
104 #ifdef EDM_ML_DEBUG
105  edm::LogVerbatim("HGCSim") << "HGCalTB16SD01: Detector " << det << " Layer " << lay << " x " << x << " " << ix << " "
106  << ixx << " y " << y << " " << iy << " " << iyy << " ID " << std::hex << idx << std::dec;
107 #endif
108  return idx;
109 }

References TauDecayModes::dec, heavyIonCSV_trainingSettings::idx, ixx, iyy, x, and y.

Referenced by setDetUnitId().

◆ setDetUnitId()

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

Implements CaloSD.

Definition at line 76 of file HGCalTB16SD01.cc.

76  {
77  const G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
78  const G4VTouchable* touch = preStepPoint->GetTouchable();
79 
80  int det(1), x(0), y(0);
81  int lay = (touch->GetReplicaNumber(0));
82 
83  return packIndex(det, lay, x, y);
84 }

References packIndex(), x, and y.

◆ unpackIndex()

void HGCalTB16SD01::unpackIndex ( const uint32_t &  idx,
int &  det,
int &  lay,
int &  x,
int &  y 
)
static

Definition at line 111 of file HGCalTB16SD01.cc.

111  {
112  det = (idx >> 28) & 15;
113  lay = (idx >> 21) & 127;
114  y = (idx >> 10) & 511;
115  if (((idx >> 19) & 1) == 1)
116  y = -y;
117  x = (idx)&511;
118  if (((idx >> 9) & 1) == 1)
119  x = -x;
120 }

References heavyIonCSV_trainingSettings::idx, x, and y.

Member Data Documentation

◆ birk1_

double HGCalTB16SD01::birk1_
private

Definition at line 35 of file HGCalTB16SD01.cc.

Referenced by getEnergyDeposit(), and HGCalTB16SD01().

◆ birk2_

double HGCalTB16SD01::birk2_
private

Definition at line 35 of file HGCalTB16SD01.cc.

Referenced by getEnergyDeposit(), and HGCalTB16SD01().

◆ birk3_

double HGCalTB16SD01::birk3_
private

Definition at line 35 of file HGCalTB16SD01.cc.

Referenced by getEnergyDeposit(), and HGCalTB16SD01().

◆ initialize_

bool HGCalTB16SD01::initialize_
private

Definition at line 36 of file HGCalTB16SD01.cc.

Referenced by getEnergyDeposit(), and initialize().

◆ matName_

std::string HGCalTB16SD01::matName_
private

Definition at line 33 of file HGCalTB16SD01.cc.

Referenced by HGCalTB16SD01(), and initialize().

◆ matScin_

G4Material* HGCalTB16SD01::matScin_
private

Definition at line 37 of file HGCalTB16SD01.cc.

Referenced by getEnergyDeposit(), HGCalTB16SD01(), and initialize().

◆ useBirk_

bool HGCalTB16SD01::useBirk_
private

Definition at line 34 of file HGCalTB16SD01.cc.

Referenced by getEnergyDeposit(), and HGCalTB16SD01().

DDAxes::y
AlCaHLTBitMon_ParallelJobs.p
p
Definition: AlCaHLTBitMon_ParallelJobs.py:153
CaloSD::getAttenuation
double getAttenuation(const G4Step *aStep, double birk1, double birk2, double birk3) const
Definition: CaloSD.cc:517
MeV
const double MeV
HGCalTB16SD01::matName_
std::string matName_
Definition: HGCalTB16SD01.cc:33
DDAxes::x
iyy
int iyy[18][41][3]
Definition: EcalDAQHandler.cc:317
heavyIonCSV_trainingSettings.idx
idx
Definition: heavyIonCSV_trainingSettings.py:5
HGCalTB16SD01::packIndex
static uint32_t packIndex(int det, int lay, int x, int y)
Definition: HGCalTB16SD01.cc:86
HGCalTB16SD01::matScin_
G4Material * matScin_
Definition: HGCalTB16SD01.cc:37
HGCalTB16SD01::useBirk_
bool useBirk_
Definition: HGCalTB16SD01.cc:34
HGCalTB16SD01::initialize
void initialize(const G4StepPoint *point)
Definition: HGCalTB16SD01.cc:122
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
HGCalTB16SD01::birk2_
double birk2_
Definition: HGCalTB16SD01.cc:35
edm::ParameterSet
Definition: ParameterSet.h:47
CaloSD::CaloSD
CaloSD(const std::string &aSDname, const edm::EventSetup &es, const SensitiveDetectorCatalog &clg, edm::ParameterSet const &p, const SimTrackManager *, float timeSlice=1., bool ignoreTkID=false)
Definition: CaloSD.cc:33
HGCalTB16SD01::birk1_
double birk1_
Definition: HGCalTB16SD01.cc:35
ixx
int ixx[18][41][3]
Definition: EcalDAQHandler.cc:226
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
Skims_PA_cff.name
name
Definition: Skims_PA_cff.py:17
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
HGCalTB16SD01::birk3_
double birk3_
Definition: HGCalTB16SD01.cc:35
HGCalTB16SD01::initialize_
bool initialize_
Definition: HGCalTB16SD01.cc:36
point
*vegas h *****************************************************used in the default bin number in original ***version of VEGAS is ***a higher bin number might help to derive a more precise ***grade subtle point
Definition: invegas.h:5
TauDecayModes.dec
dec
Definition: TauDecayModes.py:143
weight
Definition: weight.py:1
g
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