CMS 3D CMS Logo

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

#include <HGCalSD.h>

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

Public Member Functions

 HGCalSD (const std::string &, const DDCompactView &, const SensitiveDetectorCatalog &, edm::ParameterSet const &, const SimTrackManager *)
 
uint32_t setDetUnitId (const G4Step *step) override
 
 ~HGCalSD () override=default
 
- Public Member Functions inherited from CaloSD
 CaloSD (const std::string &aSDname, const DDCompactView &cpv, 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
 
void PrintAll () override
 
G4bool ProcessHits (G4Step *step, G4TouchableHistory *) override
 
bool ProcessHits (G4GFlashSpot *aSpot, G4TouchableHistory *) override
 
 ~CaloSD () override
 
- Public Member Functions inherited from SensitiveCaloDetector
 SensitiveCaloDetector (const std::string &iname, const DDCompactView &cpv, 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 DDCompactView &cpv, 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 ()
 
- Public Member Functions inherited from Observer< const BeginOfJob * >
 Observer ()
 
void slotForUpdate (const BeginOfJob * iT)
 
virtual ~Observer ()
 

Protected Member Functions

bool filterHit (CaloG4Hit *, double) override
 
double getEnergyDeposit (const G4Step *) 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 *)
 
virtual void endEvent ()
 
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 *)
 
virtual void initEvent (const BeginOfEvent *)
 
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 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 (int, int, int, int, G4ThreeVector &)
 

Private Attributes

std::vector< double > angles_
 
int cornerMinMask_
 
double distanceFromEdge_
 
double eminHit_
 
bool fiducialCut_
 
HGCalGeometryMode::GeometryMode geom_mode_
 
const HGCalDDDConstantshgcons_
 
int levelT1_
 
int levelT2_
 
std::unique_ptr< HGCMouseBitemouseBite_
 
double mouseBiteCut_
 
DetId::Detector mydet_
 
std::string nameX_
 
std::unique_ptr< HGCalNumberingSchemenumberingScheme_
 
bool rejectMB_
 
double slopeMin_
 
bool storeAllG4Hits_
 
const double tan30deg_
 
bool waferRot_
 
double weight_
 

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 21 of file HGCalSD.h.

Constructor & Destructor Documentation

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

Definition at line 30 of file HGCalSD.cc.

References angles_, cornerMinMask_, distanceFromEdge_, eminHit_, fiducialCut_, DetId::Forward, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), DetId::HGCalEE, DetId::HGCalHSi, MeV, mouseBite_, mydet_, dataset::name, nameX_, numberingScheme_, rejectMB_, CaloSD::setNumberCheckedHits(), CaloSD::setUseMap(), storeAllG4Hits_, and waferRot_.

32  :
33  CaloSD(name, cpv, clg, p, manager,
34  (float)(p.getParameter<edm::ParameterSet>("HGCSD").getParameter<double>("TimeSliceUnit")),
35  p.getParameter<edm::ParameterSet>("HGCSD").getParameter<bool>("IgnoreTrackID")),
36  hgcons_(nullptr), slopeMin_(0), levelT1_(99), levelT2_(99),
37  tan30deg_(std::tan(30.0*CLHEP::deg)) {
38 
39  numberingScheme_.reset(nullptr); mouseBite_.reset(nullptr);
40 
41  edm::ParameterSet m_HGC = p.getParameter<edm::ParameterSet>("HGCSD");
42  eminHit_ = m_HGC.getParameter<double>("EminHit")*CLHEP::MeV;
43  fiducialCut_ = m_HGC.getParameter<bool>("FiducialCut");
44  distanceFromEdge_= m_HGC.getParameter<double>("DistanceFromEdge");
45  storeAllG4Hits_ = m_HGC.getParameter<bool>("StoreAllG4Hits");
46  rejectMB_ = m_HGC.getParameter<bool>("RejectMouseBite");
47  waferRot_ = m_HGC.getParameter<bool>("RotatedWafer");
48  cornerMinMask_ = m_HGC.getParameter<int>("CornerMinMask");
49  angles_ = m_HGC.getUntrackedParameter<std::vector<double>>("WaferAngles");
50 
51  if(storeAllG4Hits_) {
52  setUseMap(false);
54  }
55 
56  //this is defined in the hgcsens.xml
57  G4String myName = name;
59  nameX_ = "HGCal";
60  if (myName.find("HitsEE")!=std::string::npos) {
62  nameX_ = "HGCalEESensitive";
63  } else if (myName.find("HitsHEfront")!=std::string::npos) {
65  nameX_ = "HGCalHESiliconSensitive";
66  }
67 
68 #ifdef EDM_ML_DEBUG
69  edm::LogVerbatim("HGCSim")<< "**************************************************"
70  << "\n"
71  << "* *"
72  << "\n"
73  << "* Constructing a HGCalSD with name " << name << "\n"
74  << "* *"
75  << "\n"
76  << "**************************************************";
77 #endif
78  edm::LogVerbatim("HGCSim") << "HGCalSD:: Threshold for storing hits: "
79  << eminHit_ << " for " << nameX_ << " detector "
80  << mydet_;
81  edm::LogVerbatim("HGCSim") << "Flag for storing individual Geant4 Hits "
82  << storeAllG4Hits_;
83  edm::LogVerbatim("HGCSim") << "Fiducial volume cut with cut from eta/phi "
84  << "boundary " << fiducialCut_ << " at "
86  edm::LogVerbatim("HGCSim") << "Reject MosueBite Flag: " << rejectMB_
87  << " cuts along " << angles_.size() << " axes: "
88  << angles_[0] << ", " << angles_[1];
89 }
T getParameter(std::string const &) const
T getUntrackedParameter(std::string const &, T const &) const
double distanceFromEdge_
Definition: HGCalSD.h:51
int levelT2_
Definition: HGCalSD.h:53
void setNumberCheckedHits(int val)
Definition: CaloSD.h:111
double slopeMin_
Definition: HGCalSD.h:51
std::unique_ptr< HGCalNumberingScheme > numberingScheme_
Definition: HGCalSD.h:46
std::string nameX_
Definition: HGCalSD.h:49
void setUseMap(bool val)
Definition: CaloSD.h:100
const double MeV
bool rejectMB_
Definition: HGCalSD.h:55
DetId::Detector mydet_
Definition: HGCalSD.h:48
bool storeAllG4Hits_
Definition: HGCalSD.h:54
std::vector< double > angles_
Definition: HGCalSD.h:57
const HGCalDDDConstants * hgcons_
Definition: HGCalSD.h:45
Tan< T >::type tan(const T &t)
Definition: Tan.h:22
const double tan30deg_
Definition: HGCalSD.h:56
bool waferRot_
Definition: HGCalSD.h:55
CaloSD(const std::string &aSDname, const DDCompactView &cpv, const SensitiveDetectorCatalog &clg, edm::ParameterSet const &p, const SimTrackManager *, float timeSlice=1., bool ignoreTkID=false)
Definition: CaloSD.cc:28
int cornerMinMask_
Definition: HGCalSD.h:53
bool fiducialCut_
Definition: HGCalSD.h:55
int levelT1_
Definition: HGCalSD.h:53
std::unique_ptr< HGCMouseBite > mouseBite_
Definition: HGCalSD.h:47
double eminHit_
Definition: HGCalSD.h:51
HGCalSD::~HGCalSD ( )
overridedefault

Member Function Documentation

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

Reimplemented from CaloSD.

Definition at line 236 of file HGCalSD.cc.

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

236  {
237  return ((time <= tmaxHit) && (aHit->getEnergyDeposit() > eminHit_));
238 }
double tmaxHit
Definition: CaloSD.h:134
double eminHit_
Definition: HGCalSD.h:51
double getEnergyDeposit() const
Definition: CaloG4Hit.h:81
double HGCalSD::getEnergyDeposit ( const G4Step *  aStep)
overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 91 of file HGCalSD.cc.

References funct::abs(), CaloSD::getResponseWt(), alignCSCRings::r, slopeMin_, weight_, and z.

91  {
92 
93  double r = aStep->GetPreStepPoint()->GetPosition().perp();
94  double z = std::abs(aStep->GetPreStepPoint()->GetPosition().z());
95 #ifdef EDM_ML_DEBUG
96  G4int parCode = aStep->GetTrack()->GetDefinition()->GetPDGEncoding();
97  G4String parName = aStep->GetTrack()->GetDefinition()->GetParticleName();
98  G4LogicalVolume* lv =
99  aStep->GetPreStepPoint()->GetPhysicalVolume()->GetLogicalVolume();
100  edm::LogVerbatim("HGCSim") << "HGCalSD: Hit from standard path from "
101  << lv->GetName() << " for Track "
102  << aStep->GetTrack()->GetTrackID() << " ("
103  << parCode << ":" << parName << ") R = " << r
104  << " Z = " << z << " slope = " << r/z << ":"
105  << slopeMin_;
106 #endif
107  // Apply fiducial cut
108  if (r < z*slopeMin_) {
109 #ifdef EDM_ML_DEBUG
110  edm::LogVerbatim("HGCSim") << "HGCalSD: Fiducial Volume cut";
111 #endif
112  return 0.0;
113  }
114 
115  double wt1 = getResponseWt(aStep->GetTrack());
116  double wt2 = aStep->GetTrack()->GetWeight();
117  double destep = weight_*wt1*(aStep->GetTotalEnergyDeposit());
118  if (wt2 > 0) destep *= wt2;
119 #ifdef EDM_ML_DEBUG
120  edm::LogVerbatim("HGCSim") << "HGCalSD: weights= " << weight_ << ":"
121  << wt1 << ":" << wt2 << " Total weight "
122  << weight_*wt1*wt2 << " deStep: "
123  << aStep->GetTotalEnergyDeposit()
124  << ":" <<destep;
125 #endif
126  return destep;
127 }
double slopeMin_
Definition: HGCalSD.h:51
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
double weight_
Definition: HGCalSD.h:52
double getResponseWt(const G4Track *)
Definition: CaloSD.cc:652
void HGCalSD::initRun ( )
overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 233 of file HGCalSD.cc.

233  {
234 }
bool HGCalSD::isItinFidVolume ( const G4ThreeVector &  pos)
private

Definition at line 250 of file HGCalSD.cc.

References distanceFromEdge_, HGCalDDDConstants::distFromEdgeHex(), fiducialCut_, and hgcons_.

250  {
251  if (fiducialCut_) {
252  return (hgcons_->distFromEdgeHex(pos.x(),pos.y(),pos.z()) > distanceFromEdge_);
253  } else {
254  return true;
255  }
256 }
double distanceFromEdge_
Definition: HGCalSD.h:51
const HGCalDDDConstants * hgcons_
Definition: HGCalSD.h:45
bool fiducialCut_
Definition: HGCalSD.h:55
double distFromEdgeHex(double x, double y, double z) const
uint32_t HGCalSD::setDetUnitId ( const G4Step *  step)
overridevirtual

Implements CaloSD.

Definition at line 129 of file HGCalSD.cc.

References TauDecayModes::dec, triggerObjects_cff::id, levelT1_, levelT2_, mouseBite_, rejectMB_, and HGCSiliconDetId::waferUV().

129  {
130 
131  const G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
132  const G4VTouchable* touch = preStepPoint->GetTouchable();
133 
134  //determine the exact position in global coordinates in the mass geometry
135  G4ThreeVector hitPoint = preStepPoint->GetPosition();
136  float globalZ = touch->GetTranslation(0).z();
137  int iz( globalZ>0 ? 1 : -1);
138 
139  int layer, module, cell;
140  if ((touch->GetHistoryDepth() == levelT1_) ||
141  (touch->GetHistoryDepth() == levelT2_)) {
142  layer = touch->GetReplicaNumber(0);
143  module = -1;
144  cell = -1;
145 #ifdef EDM_ML_DEBUG
146  edm::LogVerbatim("HGCSim") << "DepthsTop: " << touch->GetHistoryDepth()
147  << ":" << levelT1_ << ":" << levelT2_
148  << " name " << touch->GetVolume(0)->GetName()
149  << " layer:module:cell " << layer << ":"
150  << module << ":" << cell;
151 #endif
152  } else {
153  layer = touch->GetReplicaNumber(3);
154  module = touch->GetReplicaNumber(2);
155  cell = touch->GetReplicaNumber(1);
156 #ifdef EDM_ML_DEBUG
157  edm::LogVerbatim("HGCSim") << "DepthsInside: " << touch->GetHistoryDepth()
158  << " name " << touch->GetVolume(0)->GetName()
159  << " layer:module:cell " << layer << ":"
160  << module << ":" << cell;
161 #endif
162  }
163 #ifdef EDM_ML_DEBUG
164  G4Material* mat = aStep->GetPreStepPoint()->GetMaterial();
165  edm::LogVerbatim("HGCSim") << "Depths: " << touch->GetHistoryDepth()
166  << " name " << touch->GetVolume(0)->GetName()
167  << ":" << touch->GetReplicaNumber(0) << " "
168  << touch->GetVolume(1)->GetName()
169  << ":" << touch->GetReplicaNumber(1) << " "
170  << touch->GetVolume(2)->GetName()
171  << ":" << touch->GetReplicaNumber(2) << " "
172  << touch->GetVolume(3)->GetName()
173  << ":" << touch->GetReplicaNumber(3) << " "
174  << touch->GetVolume(4)->GetName()
175  << ":" << touch->GetReplicaNumber(4) << " "
176  << " layer:module:cell " << layer << ":"
177  << module << ":" << cell << " Material "
178  << mat->GetName() << ":" << mat->GetRadlen();
179 #endif
180  // The following statement should be examined later before elimination
181  if (aStep->GetPreStepPoint()->GetMaterial()->GetRadlen() > 100000.) return 0;
182 
183  uint32_t id = setDetUnitId (layer, module, cell, iz, hitPoint);
184  if (rejectMB_ && id != 0) {
185  auto uv = HGCSiliconDetId(id).waferUV();
186 #ifdef EDM_ML_DEBUG
187  edm::LogVerbatim("HGCSim") << "ID " << std::hex << id << std::dec
188  << " " << HGCSiliconDetId(id);
189 #endif
190  if (mouseBite_->exclude(hitPoint, iz, uv.first, uv.second)) {
191  id = 0;
192 #ifdef EDM_ML_DEBUG
193  edm::LogVerbatim("HGCSim") << "Rejected by mousebite cutoff *****";
194 #endif
195  }
196  }
197 #ifdef EDM_ML_DEBUG
198  if (id != 0) edm::LogVerbatim("HGCSim") << HGCSiliconDetId(id);
199 #endif
200  return id;
201 }
int levelT2_
Definition: HGCalSD.h:53
bool rejectMB_
Definition: HGCalSD.h:55
uint32_t setDetUnitId(const G4Step *step) override
Definition: HGCalSD.cc:129
std::pair< int, int > waferUV() const
int levelT1_
Definition: HGCalSD.h:53
std::unique_ptr< HGCMouseBite > mouseBite_
Definition: HGCalSD.h:47
Definition: vlib.h:208
uint32_t HGCalSD::setDetUnitId ( int  layer,
int  module,
int  cell,
int  iz,
G4ThreeVector &  pos 
)
private

Definition at line 240 of file HGCalSD.cc.

References cornerMinMask_, hgcons_, triggerObjects_cff::id, HGCalDDDConstants::maskCell(), numberingScheme_, and weight_.

241  {
242  uint32_t id = numberingScheme_ ?
243  numberingScheme_->getUnitID(layer, module, cell, iz, pos, weight_) : 0;
244  if (cornerMinMask_ > 2) {
245  if (hgcons_->maskCell(DetId(id), cornerMinMask_)) id = 0;
246  }
247  return id;
248 }
std::unique_ptr< HGCalNumberingScheme > numberingScheme_
Definition: HGCalSD.h:46
const HGCalDDDConstants * hgcons_
Definition: HGCalSD.h:45
Definition: DetId.h:18
int cornerMinMask_
Definition: HGCalSD.h:53
bool maskCell(const DetId &id, int corners) const
double weight_
Definition: HGCalSD.h:52
Definition: vlib.h:208
void HGCalSD::update ( const BeginOfJob )
overrideprotectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob * >.

Definition at line 203 of file HGCalSD.cc.

References angles_, Exception, geom_mode_, HGCalDDDConstants::geomMode(), edm::EventSetup::get(), hgcons_, levelT1_, levelT2_, HGCalDDDConstants::levelTop(), HGCalDDDConstants::minSlope(), HGCalDDDConstants::mouseBite(), mouseBite_, mouseBiteCut_, mydet_, nameX_, numberingScheme_, rejectMB_, slopeMin_, tan30deg_, waferRot_, and HGCalDDDConstants::waferSize().

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

203  {
204 
205  const edm::EventSetup* es = (*job)();
207  es->get<IdealGeometryRecord>().get(nameX_,hdc);
208  if (hdc.isValid()) {
209  hgcons_ = hdc.product();
212  levelT1_ = hgcons_->levelTop(0);
213  levelT2_ = hgcons_->levelTop(1);
214  double waferSize = hgcons_->waferSize(false);
215  double mouseBite = hgcons_->mouseBite(false);
216  mouseBiteCut_ = waferSize*tan30deg_ - mouseBite;
217 #ifdef EDM_ML_DEBUG
218  edm::LogVerbatim("HGCSim") << "HGCalSD::Initialized with mode "
219  << geom_mode_ << " Slope cut " << slopeMin_
220  << " top Level " << levelT1_ << ":" << levelT2_
221  << " wafer " << waferSize << ":" << mouseBite;
222 #endif
223 
225  if (rejectMB_)
227  waferRot_));
228  } else {
229  throw cms::Exception("Unknown", "HGCalSD") << "Cannot find HGCalDDDConstants for " << nameX_ << "\n";
230  }
231 }
int levelT2_
Definition: HGCalSD.h:53
double slopeMin_
Definition: HGCalSD.h:51
std::unique_ptr< HGCalNumberingScheme > numberingScheme_
Definition: HGCalSD.h:46
std::string nameX_
Definition: HGCalSD.h:49
double mouseBiteCut_
Definition: HGCalSD.h:52
bool rejectMB_
Definition: HGCalSD.h:55
DetId::Detector mydet_
Definition: HGCalSD.h:48
double mouseBite(bool reco) const
std::vector< double > angles_
Definition: HGCalSD.h:57
const HGCalDDDConstants * hgcons_
Definition: HGCalSD.h:45
double minSlope() const
const double tan30deg_
Definition: HGCalSD.h:56
HGCalGeometryMode::GeometryMode geomMode() const
double waferSize(bool reco) const
bool waferRot_
Definition: HGCalSD.h:55
HGCalGeometryMode::GeometryMode geom_mode_
Definition: HGCalSD.h:50
T get() const
Definition: EventSetup.h:68
int levelT1_
Definition: HGCalSD.h:53
std::unique_ptr< HGCMouseBite > mouseBite_
Definition: HGCalSD.h:47
int levelTop(int ind=0) const

Member Data Documentation

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

Definition at line 57 of file HGCalSD.h.

Referenced by HGCalSD(), and update().

int HGCalSD::cornerMinMask_
private

Definition at line 53 of file HGCalSD.h.

Referenced by HGCalSD(), and setDetUnitId().

double HGCalSD::distanceFromEdge_
private

Definition at line 51 of file HGCalSD.h.

Referenced by HGCalSD(), and isItinFidVolume().

double HGCalSD::eminHit_
private

Definition at line 51 of file HGCalSD.h.

Referenced by filterHit(), and HGCalSD().

bool HGCalSD::fiducialCut_
private

Definition at line 55 of file HGCalSD.h.

Referenced by HGCalSD(), and isItinFidVolume().

HGCalGeometryMode::GeometryMode HGCalSD::geom_mode_
private

Definition at line 50 of file HGCalSD.h.

Referenced by update().

const HGCalDDDConstants* HGCalSD::hgcons_
private

Definition at line 45 of file HGCalSD.h.

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

int HGCalSD::levelT1_
private

Definition at line 53 of file HGCalSD.h.

Referenced by setDetUnitId(), and update().

int HGCalSD::levelT2_
private

Definition at line 53 of file HGCalSD.h.

Referenced by setDetUnitId(), and update().

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

Definition at line 47 of file HGCalSD.h.

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

double HGCalSD::mouseBiteCut_
private

Definition at line 52 of file HGCalSD.h.

Referenced by update().

DetId::Detector HGCalSD::mydet_
private

Definition at line 48 of file HGCalSD.h.

Referenced by HGCalSD(), and update().

std::string HGCalSD::nameX_
private

Definition at line 49 of file HGCalSD.h.

Referenced by HGCalSD(), and update().

std::unique_ptr<HGCalNumberingScheme> HGCalSD::numberingScheme_
private

Definition at line 46 of file HGCalSD.h.

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

bool HGCalSD::rejectMB_
private

Definition at line 55 of file HGCalSD.h.

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

double HGCalSD::slopeMin_
private

Definition at line 51 of file HGCalSD.h.

Referenced by getEnergyDeposit(), and update().

bool HGCalSD::storeAllG4Hits_
private

Definition at line 54 of file HGCalSD.h.

Referenced by HGCalSD().

const double HGCalSD::tan30deg_
private

Definition at line 56 of file HGCalSD.h.

Referenced by update().

bool HGCalSD::waferRot_
private

Definition at line 55 of file HGCalSD.h.

Referenced by HGCalSD(), and update().

double HGCalSD::weight_
private

Definition at line 52 of file HGCalSD.h.

Referenced by getEnergyDeposit(), and setDetUnitId().