CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
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 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
< HGCNumberingScheme
numberingScheme_
 
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 ( 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_, mergeVDriftHistosByStation::name, nameX_, numberingScheme_, 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 getUntrackedParameter(std::string const &, T const &) const
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
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
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
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 ( )
overridedefault

Member Function Documentation

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
bool HGCSD::filterHit ( CaloG4Hit aHit,
double  time 
)
overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 255 of file HGCSD.cc.

References eminHit_, CaloG4Hit::getEnergyDeposit(), 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
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
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
void HGCSD::initRun ( )
overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 219 of file HGCSD.cc.

References edm::Service< T >::isAvailable(), TFileService::make(), 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
T * make(const Args &...args) const
make new ROOT object
Definition: TFileService.h:64
bool isAvailable() const
Definition: Service.h:40
std::vector< int > t_Parcode_
Definition: HGCSD.h:61
std::vector< double > t_Angle_
Definition: HGCSD.h:63
TTree * tree_
Definition: HGCSD.h:59
tuple tfile
Definition: compare.py:324
uint32_t t_EventID_
Definition: HGCSD.h:60
bool HGCSD::isItinFidVolume ( const G4ThreeVector &  )
inlineprivate

Definition at line 44 of file HGCSD.h.

44 { return true; }
uint32_t HGCSD::setDetUnitId ( const G4Step *  step)
overridevirtual

Implements CaloSD.

Definition at line 141 of file HGCSD.cc.

References TauDecayModes::dec, gpuClustering::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
uint16_t *__restrict__ id
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)
tuple module
Definition: callgraph.py:69
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 gpuClustering::id, and numberingScheme_.

259  {
260  uint32_t id = numberingScheme_ ? numberingScheme_->getUnitID(subdet, layer, module, cell, iz, pos) : 0;
261  return id;
262 }
uint16_t *__restrict__ id
constexpr std::array< uint8_t, layerIndexSize > layer
std::unique_ptr< HGCNumberingScheme > numberingScheme_
Definition: HGCSD.h:49
tuple module
Definition: callgraph.py:69
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
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
double minSlope() const
HGCalGeometryMode::GeometryMode geomMode() const
std::unique_ptr< HGCNumberingScheme > numberingScheme_
Definition: HGCSD.h:49
bool rejectMB_
Definition: HGCSD.h:55
int levelTop(int ind=0) const

Member Data Documentation

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

Definition at line 57 of file HGCSD.h.

Referenced by HGCSD(), and update().

double HGCSD::eminHit_
private

Definition at line 51 of file HGCSD.h.

Referenced by filterHit(), and HGCSD().

HGCalGeometryMode::GeometryMode HGCSD::geom_mode_
private

Definition at line 48 of file HGCSD.h.

Referenced by update().

const HGCalDDDConstants* HGCSD::hgcons_
private

Definition at line 46 of file HGCSD.h.

Referenced by update().

int HGCSD::levelT_
private

Definition at line 54 of file HGCSD.h.

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

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

Definition at line 50 of file HGCSD.h.

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

double HGCSD::mouseBiteCut_
private

Definition at line 56 of file HGCSD.h.

Referenced by HGCSD(), and update().

ForwardSubdetector HGCSD::myFwdSubdet_
private

Definition at line 52 of file HGCSD.h.

Referenced by HGCSD(), and setDetUnitId().

std::string HGCSD::nameX_
private

Definition at line 47 of file HGCSD.h.

Referenced by HGCSD(), and update().

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

Definition at line 49 of file HGCSD.h.

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

bool HGCSD::rejectMB_
private

Definition at line 55 of file HGCSD.h.

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

double HGCSD::slopeMin_
private

Definition at line 53 of file HGCSD.h.

Referenced by getEnergyDeposit(), and update().

bool HGCSD::storeAllG4Hits_
private

Definition at line 55 of file HGCSD.h.

Referenced by HGCSD().

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

Definition at line 63 of file HGCSD.h.

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

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

Definition at line 62 of file HGCSD.h.

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

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

Definition at line 62 of file HGCSD.h.

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

uint32_t HGCSD::t_EventID_
private

Definition at line 60 of file HGCSD.h.

Referenced by initEvent(), and initRun().

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

Definition at line 61 of file HGCSD.h.

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

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

Definition at line 61 of file HGCSD.h.

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

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

Definition at line 62 of file HGCSD.h.

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

TTree* HGCSD::tree_
private

Definition at line 59 of file HGCSD.h.

Referenced by endEvent(), and initRun().

bool HGCSD::waferRot_
private

Definition at line 55 of file HGCSD.h.

Referenced by HGCSD(), and update().