CMS 3D CMS Logo

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

#include <HGCSD.h>

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

Public Member Functions

 HGCSD (const std::string &, const HGCalDDDConstants *, const SensitiveDetectorCatalog &, edm::ParameterSet const &, const SimTrackManager *)
 
uint32_t setDetUnitId (const G4Step *step) override
 
 ~HGCSD () override=default
 
- 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)
 
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
 
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)
 
- 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)
 
 ~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 ()
 
- Public Member Functions inherited from Observer< const BeginOfJob *>
 Observer ()
 
void slotForUpdate (const BeginOfJob * iT)
 
virtual ~Observer ()
 

Protected Member Functions

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

bool isItinFidVolume (const G4ThreeVector &)
 
uint32_t setDetUnitId (ForwardSubdetector &, int, int, int, int, G4ThreeVector &)
 

Private Attributes

std::vector< double > angles_
 
double eminHit_
 
HGCalGeometryMode::GeometryMode geom_mode_
 
const HGCalDDDConstantshgcons_
 
int levelT_
 
std::unique_ptr< HGCMouseBitemouseBite_
 
double mouseBiteCut_
 
ForwardSubdetector myFwdSubdet_
 
std::string nameX_
 
std::unique_ptr< HGCNumberingSchemenumberingScheme_
 
bool rejectMB_
 
double slopeMin_
 
bool storeAllG4Hits_
 
std::vector< double > t_Angle_
 
std::vector< double > t_dEStep1_
 
std::vector< double > t_dEStep2_
 
uint32_t t_EventID_
 
std::vector< int > t_Layer_
 
std::vector< int > t_Parcode_
 
std::vector< double > t_TrackE_
 
TTree * tree_
 
bool waferRot_
 

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
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 22 of file HGCSD.h.

Constructor & Destructor Documentation

◆ HGCSD()

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

Definition at line 32 of file HGCSD.cc.

References angles_, CaloSD::eminHit, eminHit_, ForwardEmpty, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HGCEE, HGCHEB, HGCHEF, mouseBite_, mouseBiteCut_, myFwdSubdet_, Skims_PA_cff::name, nameX_, numberingScheme_, AlCaHLTBitMon_ParallelJobs::p, rejectMB_, CaloSD::setNumberCheckedHits(), CaloSD::setUseMap(), storeAllG4Hits_, funct::tan(), and waferRot_.

37  : CaloSD(name,
38  clg,
39  p,
40  manager,
41  (float)(p.getParameter<edm::ParameterSet>("HGCSD").getParameter<double>("TimeSliceUnit")),
42  p.getParameter<edm::ParameterSet>("HGCSD").getParameter<bool>("IgnoreTrackID")),
43  hgcons_(hgc),
44  slopeMin_(0),
45  levelT_(99),
46  tree_(nullptr) {
47  numberingScheme_.reset(nullptr);
48  mouseBite_.reset(nullptr);
49 
50  edm::ParameterSet m_HGC = p.getParameter<edm::ParameterSet>("HGCSD");
51  eminHit_ = m_HGC.getParameter<double>("EminHit") * CLHEP::MeV;
52  storeAllG4Hits_ = m_HGC.getParameter<bool>("StoreAllG4Hits");
53  rejectMB_ = m_HGC.getParameter<bool>("RejectMouseBite");
54  waferRot_ = m_HGC.getParameter<bool>("RotatedWafer");
55  angles_ = m_HGC.getUntrackedParameter<std::vector<double>>("WaferAngles");
56  double waferSize = m_HGC.getUntrackedParameter<double>("WaferSize") * CLHEP::mm;
57  double mouseBite = m_HGC.getUntrackedParameter<double>("MouseBite") * CLHEP::mm;
58  mouseBiteCut_ = waferSize * tan(30.0 * CLHEP::deg) - mouseBite;
59 
60  if (storeAllG4Hits_) {
61  setUseMap(false);
63  }
64  //this is defined in the hgcsens.xml
65  G4String myName = name;
67  nameX_ = "HGCal";
68  if (myName.find("HitsEE") != std::string::npos) {
70  nameX_ = "HGCalEESensitive";
71  } else if (myName.find("HitsHEfront") != std::string::npos) {
73  nameX_ = "HGCalHESiliconSensitive";
74  } else if (myName.find("HitsHEback") != std::string::npos) {
76  nameX_ = "HGCalHEScintillatorSensitive";
77  }
78 
79 #ifdef EDM_ML_DEBUG
80  edm::LogVerbatim("HGCSim") << "**************************************************"
81  << "\n"
82  << "* *"
83  << "\n"
84  << "* Constructing a HGCSD with name " << name << "\n"
85  << "* *"
86  << "\n"
87  << "**************************************************";
88 #endif
89  edm::LogVerbatim("HGCSim") << "HGCSD:: Threshold for storing hits: " << eminHit << " for " << nameX_ << " subdet "
90  << myFwdSubdet_;
91  edm::LogVerbatim("HGCSim") << "Flag for storing individual Geant4 Hits " << storeAllG4Hits_;
92  edm::LogVerbatim("HGCSim") << "Reject MosueBite Flag: " << rejectMB_ << " Size of wafer " << waferSize
93  << " Mouse Bite " << mouseBite << ":" << mouseBiteCut_ << " along " << angles_.size()
94  << " axes";
95 }
Log< level::Info, true > LogVerbatim
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
void setNumberCheckedHits(int val)
Definition: CaloSD.h:122
int levelT_
Definition: HGCSD.h:54
double mouseBiteCut_
Definition: HGCSD.h:56
std::string nameX_
Definition: HGCSD.h:47
std::vector< double > angles_
Definition: HGCSD.h:57
bool storeAllG4Hits_
Definition: HGCSD.h:55
bool waferRot_
Definition: HGCSD.h:55
void setUseMap(bool val)
Definition: CaloSD.h:111
double eminHit
Definition: CaloSD.h:144
T getUntrackedParameter(std::string const &, T const &) const
double slopeMin_
Definition: HGCSD.h:53
std::unique_ptr< HGCMouseBite > mouseBite_
Definition: HGCSD.h:50
const HGCalDDDConstants * hgcons_
Definition: HGCSD.h:46
Tan< T >::type tan(const T &t)
Definition: Tan.h:22
std::unique_ptr< HGCNumberingScheme > numberingScheme_
Definition: HGCSD.h:49
double eminHit_
Definition: HGCSD.h:51
bool rejectMB_
Definition: HGCSD.h:55
TTree * tree_
Definition: HGCSD.h:59
CaloSD(const std::string &aSDname, const SensitiveDetectorCatalog &clg, edm::ParameterSet const &p, const SimTrackManager *, float timeSlice=1., bool ignoreTkID=false)
Definition: CaloSD.cc:33
ForwardSubdetector myFwdSubdet_
Definition: HGCSD.h:52

◆ ~HGCSD()

HGCSD::~HGCSD ( )
overridedefault

Member Function Documentation

◆ endEvent()

void HGCSD::endEvent ( )
overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 248 of file HGCSD.cc.

References tree_.

248  {
249 #ifdef plotDebug
250  if (tree_)
251  tree_->Fill();
252 #endif
253 }
TTree * tree_
Definition: HGCSD.h:59

◆ filterHit()

bool HGCSD::filterHit ( CaloG4Hit aHit,
double  time 
)
overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 255 of file HGCSD.cc.

References eminHit_, CaloG4Hit::getEnergyDeposit(), protons_cff::time, and CaloSD::tmaxHit.

255  {
256  return ((time <= tmaxHit) && (aHit->getEnergyDeposit() > eminHit_));
257 }
double tmaxHit
Definition: CaloSD.h:144
double eminHit_
Definition: HGCSD.h:51
double getEnergyDeposit() const
Definition: CaloG4Hit.h:79

◆ getEnergyDeposit()

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

Reimplemented from CaloSD.

Definition at line 97 of file HGCSD.cc.

References funct::abs(), angle(), CaloSD::getResponseWt(), phase1PixelTopology::layer, levelT_, alignCSCRings::r, slopeMin_, t_Angle_, t_dEStep1_, t_dEStep2_, t_Layer_, t_Parcode_, t_TrackE_, and z.

97  {
98  double r = aStep->GetPreStepPoint()->GetPosition().perp();
99  double z = std::abs(aStep->GetPreStepPoint()->GetPosition().z());
100 
101 #ifdef EDM_ML_DEBUG
102  G4int parCode = aStep->GetTrack()->GetDefinition()->GetPDGEncoding();
103  G4LogicalVolume* lv = aStep->GetPreStepPoint()->GetPhysicalVolume()->GetLogicalVolume();
104  edm::LogVerbatim("HGCSim") << "HGCSD: Hit from standard path from " << lv->GetName() << " for Track "
105  << aStep->GetTrack()->GetTrackID() << " ("
106  << aStep->GetTrack()->GetDefinition()->GetParticleName() << ":" << parCode << ") R = " << r
107  << " Z = " << z << " slope = " << r / z << ":" << slopeMin_;
108 #endif
109 
110  // Apply fiductial volume
111  if (r < z * slopeMin_) {
112  return 0.0;
113  }
114 
115  double wt1 = getResponseWt(aStep->GetTrack());
116  double wt2 = aStep->GetTrack()->GetWeight();
117  double destep = wt1 * aStep->GetTotalEnergyDeposit();
118  if (wt2 > 0)
119  destep *= wt2;
120 
121 #ifdef plotDebug
122  const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
123  G4double tmptrackE = aStep->GetTrack()->GetKineticEnergy();
124  G4int parCode = aStep->GetTrack()->GetDefinition()->GetPDGEncoding();
125  G4double angle = (aStep->GetTrack()->GetMomentumDirection().theta()) / CLHEP::deg;
126  G4int layer = ((touch->GetHistoryDepth() == levelT_) ? touch->GetReplicaNumber(0) : touch->GetReplicaNumber(2));
127  G4int ilayer = (layer - 1) / 3;
128  if (aStep->GetTotalEnergyDeposit() > 0) {
129  t_Layer_.emplace_back(ilayer);
130  t_Parcode_.emplace_back(parCode);
131  t_dEStep1_.emplace_back(aStep->GetTotalEnergyDeposit());
132  t_dEStep2_.emplace_back(destep);
133  t_TrackE_.emplace_back(tmptrackE);
134  t_Angle_.emplace_back(angle);
135  }
136 #endif
137 
138  return destep;
139 }
std::vector< double > t_TrackE_
Definition: HGCSD.h:62
Log< level::Info, true > LogVerbatim
std::vector< double > t_dEStep1_
Definition: HGCSD.h:62
std::vector< double > t_dEStep2_
Definition: HGCSD.h:62
int levelT_
Definition: HGCSD.h:54
std::vector< int > t_Layer_
Definition: HGCSD.h:61
constexpr std::array< uint8_t, layerIndexSize > layer
double slopeMin_
Definition: HGCSD.h:53
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
std::vector< int > t_Parcode_
Definition: HGCSD.h:61
std::vector< double > t_Angle_
Definition: HGCSD.h:63
double getResponseWt(const G4Track *)
Definition: CaloSD.cc:849
T angle(T x1, T y1, T z1, T x2, T y2, T z2)
Definition: angle.h:11

◆ initEvent()

void HGCSD::initEvent ( const BeginOfEvent g4Event)
overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 235 of file HGCSD.cc.

References t_Angle_, t_dEStep1_, t_dEStep2_, t_EventID_, t_Layer_, t_Parcode_, and t_TrackE_.

235  {
236  const G4Event* evt = (*g4Event)();
237  t_EventID_ = evt->GetEventID();
238 #ifdef plotDebug
239  t_Layer_.clear();
240  t_Parcode_.clear();
241  t_dEStep1_.clear();
242  t_dEStep2_.clear();
243  t_TrackE_.clear();
244  t_Angle_.clear();
245 #endif
246 }
std::vector< double > t_TrackE_
Definition: HGCSD.h:62
std::vector< double > t_dEStep1_
Definition: HGCSD.h:62
std::vector< double > t_dEStep2_
Definition: HGCSD.h:62
std::vector< int > t_Layer_
Definition: HGCSD.h:61
std::vector< int > t_Parcode_
Definition: HGCSD.h:61
std::vector< double > t_Angle_
Definition: HGCSD.h:63
uint32_t t_EventID_
Definition: HGCSD.h:60

◆ initRun()

void HGCSD::initRun ( )
overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 219 of file HGCSD.cc.

References t_Angle_, t_dEStep1_, t_dEStep2_, t_EventID_, t_Layer_, t_Parcode_, t_TrackE_, compare::tfile, and tree_.

219  {
220 #ifdef plotDebug
222  if (tfile.isAvailable()) {
223  tree_ = tfile->make<TTree>("TreeHGCSD", "TreeHGCSD");
224  tree_->Branch("EventID", &t_EventID_);
225  tree_->Branch("Layer", &t_Layer_);
226  tree_->Branch("ParticleCode", &t_Parcode_);
227  tree_->Branch("dEStepOriginal", &t_dEStep1_);
228  tree_->Branch("dEStepWeighted", &t_dEStep2_);
229  tree_->Branch("TrackEnergy", &t_TrackE_);
230  tree_->Branch("ThetaAngle", &t_Angle_);
231  }
232 #endif
233 }
std::vector< double > t_TrackE_
Definition: HGCSD.h:62
std::vector< double > t_dEStep1_
Definition: HGCSD.h:62
std::vector< double > t_dEStep2_
Definition: HGCSD.h:62
std::vector< int > t_Layer_
Definition: HGCSD.h:61
Definition: tfile.py:1
std::vector< int > t_Parcode_
Definition: HGCSD.h:61
std::vector< double > t_Angle_
Definition: HGCSD.h:63
TTree * tree_
Definition: HGCSD.h:59
uint32_t t_EventID_
Definition: HGCSD.h:60

◆ isItinFidVolume()

bool HGCSD::isItinFidVolume ( const G4ThreeVector &  )
inlineprivate

Definition at line 44 of file HGCSD.h.

44 { return true; }

◆ setDetUnitId() [1/2]

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

Implements CaloSD.

Definition at line 141 of file HGCSD.cc.

References TauDecayModes::dec, triggerObjects_cff::id, phase1PixelTopology::layer, levelT_, callgraph::module, mouseBite_, myFwdSubdet_, rejectMB_, HGCalTestNumbering::unpackHexagonIndex(), and z.

141  {
142  const G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
143  const G4VTouchable* touch = preStepPoint->GetTouchable();
144 
145  //determine the exact position in global coordinates in the mass geometry
146  G4ThreeVector hitPoint = preStepPoint->GetPosition();
147  float globalZ = touch->GetTranslation(0).z();
148  int iz(globalZ > 0 ? 1 : -1);
149 
150  //convert to local coordinates (=local to the current volume):
151  G4ThreeVector localpos = touch->GetHistory()->GetTopTransform().TransformPoint(hitPoint);
152 
153  //get the det unit id with
155 
156  int layer, module, cell;
157  if (touch->GetHistoryDepth() == levelT_) {
158  layer = touch->GetReplicaNumber(0);
159  module = -1;
160  cell = -1;
161 #ifdef EDM_ML_DEBUG
162  edm::LogVerbatim("HGCSim") << "Depths: " << touch->GetHistoryDepth() << " name " << touch->GetVolume(0)->GetName()
163  << " layer:module:cell " << layer << ":" << module << ":" << cell;
164 #endif
165  } else {
166  layer = touch->GetReplicaNumber(2);
167  module = touch->GetReplicaNumber(1);
168  cell = touch->GetReplicaNumber(0);
169  }
170 #ifdef EDM_ML_DEBUG
171  const G4Material* mat = aStep->GetPreStepPoint()->GetMaterial();
172  edm::LogVerbatim("HGCSim") << "Depths: " << touch->GetHistoryDepth() << " name " << touch->GetVolume(0)->GetName()
173  << ":" << touch->GetReplicaNumber(0) << " " << touch->GetVolume(1)->GetName() << ":"
174  << touch->GetReplicaNumber(1) << " " << touch->GetVolume(2)->GetName() << ":"
175  << touch->GetReplicaNumber(2) << " layer:module:cell " << layer << ":" << module << ":"
176  << cell << " Material " << mat->GetName() << ":" << mat->GetRadlen();
177 //for (int k = 0; k< touch->GetHistoryDepth(); ++k)
178 // edm::LogVerbatim("HGCSim") << "Level [" << k << "] " << touch->GetVolume(k)->GetName() << ":" << touch->GetReplicaNumber(k);
179 #endif
180  // The following statement should be examined later before elimination
181  // VI: this is likely a check if media is vacuum - not needed
182  if (aStep->GetPreStepPoint()->GetMaterial()->GetRadlen() > 100000.)
183  return 0;
184 
185  uint32_t id = setDetUnitId(subdet, layer, module, cell, iz, localpos);
186  if (rejectMB_ && id != 0) {
187  int det, z, lay, wafer, type, ic;
188  HGCalTestNumbering::unpackHexagonIndex(id, det, z, lay, wafer, type, ic);
189 #ifdef EDM_ML_DEBUG
190  edm::LogVerbatim("HGCSim") << "ID " << std::hex << id << std::dec << " Input " << subdet << ":" << layer << ":"
191  << module << ":" << cell << ":" << iz << localpos.x() << ":" << localpos.y()
192  << " Decode " << det << ":" << z << ":" << lay << ":" << wafer << ":" << type << ":"
193  << ic;
194 #endif
195  if (mouseBite_->exclude(hitPoint, z, wafer, 0))
196  id = 0;
197  }
198  return id;
199 }
Log< level::Info, true > LogVerbatim
int levelT_
Definition: HGCSD.h:54
ForwardSubdetector
constexpr std::array< uint8_t, layerIndexSize > layer
std::unique_ptr< HGCMouseBite > mouseBite_
Definition: HGCSD.h:50
bool rejectMB_
Definition: HGCSD.h:55
uint32_t setDetUnitId(const G4Step *step) override
Definition: HGCSD.cc:141
ForwardSubdetector myFwdSubdet_
Definition: HGCSD.h:52
static void unpackHexagonIndex(const uint32_t &idx, int &subdet, int &z, int &lay, int &wafer, int &celltyp, int &cell)

◆ setDetUnitId() [2/2]

uint32_t HGCSD::setDetUnitId ( ForwardSubdetector subdet,
int  layer,
int  module,
int  cell,
int  iz,
G4ThreeVector &  pos 
)
private

Definition at line 259 of file HGCSD.cc.

References triggerObjects_cff::id, phase1PixelTopology::layer, callgraph::module, and numberingScheme_.

259  {
260  uint32_t id = numberingScheme_ ? numberingScheme_->getUnitID(subdet, layer, module, cell, iz, pos) : 0;
261  return id;
262 }
constexpr std::array< uint8_t, layerIndexSize > layer
std::unique_ptr< HGCNumberingScheme > numberingScheme_
Definition: HGCSD.h:49

◆ update() [1/6]

void CaloSD::update
overrideprotected

Definition at line 715 of file CaloSD.cc.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), and MatrixUtil.Steps::overwrite().

715  {
716  endEvent();
717  slave.get()->ReserveMemory(theHC->entries());
718 
719  int count(0);
720  int wrong(0);
721  double eEM(0.0);
722  double eHAD(0.0);
723  double eEM2(0.0);
724  double eHAD2(0.0);
725  double tt(0.0);
726  double zloc(0.0);
727  double zglob(0.0);
728  double ee(0.0);
729  int hc_entries = theHC->entries();
730  for (int i = 0; i < hc_entries; ++i) {
731  if (!saveHit((*theHC)[i])) {
732  ++wrong;
733  }
734  ++count;
735  double x = (*theHC)[i]->getEM();
736  eEM += x;
737  eEM2 += x * x;
738  x = (*theHC)[i]->getHadr();
739  eHAD += x;
740  eHAD2 += x * x;
741  tt += (*theHC)[i]->getTimeSlice();
742  ee += (*theHC)[i]->getIncidentEnergy();
743  zglob += std::abs((*theHC)[i]->getEntry().z());
744  zloc += std::abs((*theHC)[i]->getEntryLocal().z());
745  }
746 
747  double norm = (count > 0) ? 1.0 / count : 0.0;
748  eEM *= norm;
749  eEM2 *= norm;
750  eHAD *= norm;
751  eHAD2 *= norm;
752  eEM2 = std::sqrt(eEM2 - eEM * eEM);
753  eHAD2 = std::sqrt(eHAD2 - eHAD * eHAD);
754  tt *= norm;
755  ee *= norm;
756  zglob *= norm;
757  zloc *= norm;
758 
759 #ifdef EDM_ML_DEBUG
760  edm::LogVerbatim("CaloSim") << "CaloSD: " << GetName() << " store " << count << " hits; " << wrong
761  << " track IDs not given properly and " << totalHits - count
762  << " hits not passing cuts\n EmeanEM= " << eEM << " ErmsEM= " << eEM2
763  << "\n EmeanHAD= " << eHAD << " ErmsHAD= " << eHAD2 << " TimeMean= " << tt
764  << " E0mean= " << ee << " Zglob= " << zglob << " Zloc= " << zloc << " ";
765 #endif
766  tkMap.erase(tkMap.begin(), tkMap.end());
767  std::vector<std::unique_ptr<CaloG4Hit>>().swap(reusehit);
768  if (useMap)
769  hitMap.erase(hitMap.begin(), hitMap.end());
771 }
Int_t getEntry(TBranch *branch, EntryNumber entryNumber)
Definition: RootTree.cc:527
Log< level::Info, true > LogVerbatim
int totalHits
Definition: CaloSD.h:177
bool useMap
Definition: CaloSD.h:171
void swap(Association< C > &lhs, Association< C > &rhs)
Definition: Association.h:117
std::unordered_map< unsigned int, unsigned int > boundaryCrossingParentMap_
Definition: CaloSD.h:189
std::unique_ptr< CaloSlaveSD > slave
Definition: CaloSD.h:163
T sqrt(T t)
Definition: SSEVec.h:19
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
std::map< CaloHitID, CaloG4Hit * > hitMap
Definition: CaloSD.h:187
CaloG4HitCollection * theHC
Definition: CaloSD.h:166
std::map< int, TrackWithHistory * > tkMap
Definition: CaloSD.h:188
std::vector< std::unique_ptr< CaloG4Hit > > reusehit
Definition: CaloSD.h:190
void endEvent() override
Definition: HGCSD.cc:248
bool saveHit(CaloG4Hit *)
Definition: CaloSD.cc:870

◆ update() [2/6]

void CaloSD::update
overrideprotected

Definition at line 933 of file CaloSD.cc.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), and MatrixUtil.Steps::overwrite().

933  {
934  int primary = -1;
935  TrackInformation* trkInfo = cmsTrackInformation((*trk)());
936  if (trkInfo->isPrimary())
937  primary = (*trk)()->GetTrackID();
938 
939 #ifdef EDM_ML_DEBUG
940  edm::LogVerbatim("CaloSim") << "New track: isPrimary " << trkInfo->isPrimary() << " primary ID = " << primary
941  << " primary ancestor ID " << primAncestor;
942 #endif
943 
944  // update the information if a different primary track ID
945 
946  if (primary > 0 && primary != primAncestor) {
947  primAncestor = primary;
948 
949  // clean the hits information
950 
951  if (theHC->entries() > 0)
953  }
954 }
Log< level::Info, true > LogVerbatim
bool isPrimary() const
int primAncestor
Definition: CaloSD.h:175
void cleanHitCollection()
Definition: CaloSD.cc:956
TrackInformation * cmsTrackInformation(const G4Track *aTrack)
CaloG4HitCollection * theHC
Definition: CaloSD.h:166

◆ update() [3/6]

void CaloSD::update
overrideprotected

Definition at line 679 of file CaloSD.cc.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), and MatrixUtil.Steps::overwrite().

679 { initRun(); }
void initRun() override
Definition: HGCSD.cc:219

◆ update() [4/6]

void CaloSD::update
overrideprotected

Definition at line 681 of file CaloSD.cc.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), and MatrixUtil.Steps::overwrite().

681  {
682 #ifdef EDM_ML_DEBUG
683  edm::LogVerbatim("CaloSim") << "CaloSD: Dispatched BeginOfEvent for " << GetName() << " !";
684 #endif
685  clearHits();
686  initEvent(g4Event);
687 }
Log< level::Info, true > LogVerbatim
void initEvent(const BeginOfEvent *) override
Definition: HGCSD.cc:235
void clearHits() override
Definition: CaloSD.cc:773

◆ update() [5/6]

void CaloSD::update
overrideprotected

Definition at line 689 of file CaloSD.cc.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), and MatrixUtil.Steps::overwrite().

689  {
690  int id = (*trk)()->GetTrackID();
691  TrackInformation* trkI = cmsTrackInformation((*trk)());
692  int lastTrackID = -1;
693  if (trkI)
694  lastTrackID = trkI->getIDonCaloSurface();
695  if (id == lastTrackID) {
696  const TrackContainer* trksForThisEvent = m_trackManager->trackContainer();
697  if (trksForThisEvent != nullptr) {
698  int it = (int)(trksForThisEvent->size()) - 1;
699  if (it >= 0) {
700  TrackWithHistory* trkH = (*trksForThisEvent)[it];
701  if (trkH->trackID() == (unsigned int)(id))
702  tkMap[id] = trkH;
703 #ifdef EDM_ML_DEBUG
704  edm::LogVerbatim("CaloSim") << "CaloSD: get track " << it << " from Container of size "
705  << trksForThisEvent->size() << " with ID " << trkH->trackID();
706  } else {
707  edm::LogVerbatim("CaloSim") << "CaloSD: get track " << it << " from Container of size "
708  << trksForThisEvent->size() << " with no ID";
709 #endif
710  }
711  }
712  }
713 }
Log< level::Info, true > LogVerbatim
std::vector< TrackWithHistory * > TrackContainer
Definition: TrackContainer.h:8
TrackInformation * cmsTrackInformation(const G4Track *aTrack)
int getIDonCaloSurface() const
const TrackContainer * trackContainer() const
std::map< int, TrackWithHistory * > tkMap
Definition: CaloSD.h:188
const SimTrackManager * m_trackManager
Definition: CaloSD.h:161
unsigned int trackID() const

◆ update() [6/6]

void HGCSD::update ( const BeginOfJob )
overrideprotectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob *>.

Definition at line 201 of file HGCSD.cc.

References angles_, Exception, geom_mode_, HGCalDDDConstants::geomMode(), hgcons_, levelT_, HGCalDDDConstants::levelTop(), HGCalDDDConstants::minSlope(), mouseBite_, mouseBiteCut_, nameX_, numberingScheme_, rejectMB_, slopeMin_, and waferRot_.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

201  {
202  if (hgcons_ != nullptr) {
205  levelT_ = hgcons_->levelTop();
206  numberingScheme_ = std::make_unique<HGCNumberingScheme>(*hgcons_, nameX_);
207  if (rejectMB_)
208  mouseBite_ = std::make_unique<HGCMouseBite>(*hgcons_, angles_, mouseBiteCut_, waferRot_);
209  } else {
210  edm::LogError("HGCSim") << "HGCSD : Cannot find HGCalDDDConstants for " << nameX_;
211  throw cms::Exception("Unknown", "HGCSD") << "Cannot find HGCalDDDConstants for " << nameX_ << "\n";
212  }
213 #ifdef EDM_ML_DEBUG
214  edm::LogVerbatim("HGCSim") << "HGCSD::Initialized with mode " << geom_mode_ << " Slope cut " << slopeMin_
215  << " top Level " << levelT_;
216 #endif
217 }
Log< level::Info, true > LogVerbatim
int levelT_
Definition: HGCSD.h:54
double mouseBiteCut_
Definition: HGCSD.h:56
std::string nameX_
Definition: HGCSD.h:47
std::vector< double > angles_
Definition: HGCSD.h:57
Log< level::Error, false > LogError
HGCalGeometryMode::GeometryMode geomMode() const
bool waferRot_
Definition: HGCSD.h:55
HGCalGeometryMode::GeometryMode geom_mode_
Definition: HGCSD.h:48
double slopeMin_
Definition: HGCSD.h:53
std::unique_ptr< HGCMouseBite > mouseBite_
Definition: HGCSD.h:50
const HGCalDDDConstants * hgcons_
Definition: HGCSD.h:46
std::unique_ptr< HGCNumberingScheme > numberingScheme_
Definition: HGCSD.h:49
bool rejectMB_
Definition: HGCSD.h:55
double minSlope() const
int levelTop(int ind=0) const

Member Data Documentation

◆ angles_

std::vector<double> HGCSD::angles_
private

Definition at line 57 of file HGCSD.h.

Referenced by HGCSD(), and update().

◆ eminHit_

double HGCSD::eminHit_
private

Definition at line 51 of file HGCSD.h.

Referenced by filterHit(), and HGCSD().

◆ geom_mode_

HGCalGeometryMode::GeometryMode HGCSD::geom_mode_
private

Definition at line 48 of file HGCSD.h.

Referenced by update().

◆ hgcons_

const HGCalDDDConstants* HGCSD::hgcons_
private

Definition at line 46 of file HGCSD.h.

Referenced by update().

◆ levelT_

int HGCSD::levelT_
private

Definition at line 54 of file HGCSD.h.

Referenced by getEnergyDeposit(), setDetUnitId(), and update().

◆ mouseBite_

std::unique_ptr<HGCMouseBite> HGCSD::mouseBite_
private

Definition at line 50 of file HGCSD.h.

Referenced by HGCSD(), setDetUnitId(), and update().

◆ mouseBiteCut_

double HGCSD::mouseBiteCut_
private

Definition at line 56 of file HGCSD.h.

Referenced by HGCSD(), and update().

◆ myFwdSubdet_

ForwardSubdetector HGCSD::myFwdSubdet_
private

Definition at line 52 of file HGCSD.h.

Referenced by HGCSD(), and setDetUnitId().

◆ nameX_

std::string HGCSD::nameX_
private

Definition at line 47 of file HGCSD.h.

Referenced by HGCSD(), and update().

◆ numberingScheme_

std::unique_ptr<HGCNumberingScheme> HGCSD::numberingScheme_
private

Definition at line 49 of file HGCSD.h.

Referenced by HGCSD(), setDetUnitId(), and update().

◆ rejectMB_

bool HGCSD::rejectMB_
private

Definition at line 55 of file HGCSD.h.

Referenced by HGCSD(), setDetUnitId(), and update().

◆ slopeMin_

double HGCSD::slopeMin_
private

Definition at line 53 of file HGCSD.h.

Referenced by getEnergyDeposit(), and update().

◆ storeAllG4Hits_

bool HGCSD::storeAllG4Hits_
private

Definition at line 55 of file HGCSD.h.

Referenced by HGCSD().

◆ t_Angle_

std::vector<double> HGCSD::t_Angle_
private

Definition at line 63 of file HGCSD.h.

Referenced by getEnergyDeposit(), initEvent(), and initRun().

◆ t_dEStep1_

std::vector<double> HGCSD::t_dEStep1_
private

Definition at line 62 of file HGCSD.h.

Referenced by getEnergyDeposit(), initEvent(), and initRun().

◆ t_dEStep2_

std::vector<double> HGCSD::t_dEStep2_
private

Definition at line 62 of file HGCSD.h.

Referenced by getEnergyDeposit(), initEvent(), and initRun().

◆ t_EventID_

uint32_t HGCSD::t_EventID_
private

Definition at line 60 of file HGCSD.h.

Referenced by initEvent(), and initRun().

◆ t_Layer_

std::vector<int> HGCSD::t_Layer_
private

Definition at line 61 of file HGCSD.h.

Referenced by getEnergyDeposit(), initEvent(), and initRun().

◆ t_Parcode_

std::vector<int> HGCSD::t_Parcode_
private

Definition at line 61 of file HGCSD.h.

Referenced by getEnergyDeposit(), initEvent(), and initRun().

◆ t_TrackE_

std::vector<double> HGCSD::t_TrackE_
private

Definition at line 62 of file HGCSD.h.

Referenced by getEnergyDeposit(), initEvent(), and initRun().

◆ tree_

TTree* HGCSD::tree_
private

Definition at line 59 of file HGCSD.h.

Referenced by endEvent(), and initRun().

◆ waferRot_

bool HGCSD::waferRot_
private

Definition at line 55 of file HGCSD.h.

Referenced by HGCSD(), and update().