#include <HFNoseSD.h>

Inheritance diagram for HFNoseSD:

Public Member Functions
	HFNoseSD (const std::string &, const edm::EventSetup &, const SensitiveDetectorCatalog &, edm::ParameterSet const &, const SimTrackManager *)

uint32_t	setDetUnitId (const G4Step *step) override

	~HFNoseSD () 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 ()

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 ::EndOfEvent *) override

void	update (const BeginOfEvent *) override

void	update (const BeginOfJob *) override
	This routine will be called when the appropriate signal arrives. More...

void	update (const BeginOfRun *) override

void	update (const BeginOfTrack *trk) override

void	update (const EndOfTrack *trk) override

Protected Member Functions inherited from CaloSD
bool	checkHit ()

CaloG4Hit *	createNewHit (const G4Step , const G4Track )

virtual void	endEvent ()

virtual double	EnergyCorrected (const G4Step &step, const G4Track *)

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 *)

void	hitBookkeepingFineCalo (const G4Step step, const G4Track currentTrack, CaloG4Hit *hit)

bool	hitExists (const G4Step *)

void	ignoreRejection ()

virtual void	initEvent (const BeginOfEvent *)

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)

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
TrackInformation *	cmsTrackInformation (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 HGCalDDDConstants *	hgcons_

int	levelT1_

int	levelT2_

std::unique_ptr< HGCMouseBite >	mouseBite_

double	mouseBiteCut_

std::string	nameX_

std::unique_ptr< HFNoseNumberingScheme >	numberingScheme_

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 }

Static Protected Member Functions inherited from CaloSD
static std::string	printableDecayChain (const std::vector< unsigned int > &decayChain)

Protected Attributes inherited from CaloSD
CaloG4Hit *	currentHit

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 20 of file HFNoseSD.h.

Constructor & Destructor Documentation

◆ HFNoseSD()

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

Definition at line 30 of file HFNoseSD.cc.

     : CaloSD(name,
              es,
              clg,
              p,
              manager,
              (float)(p.getParameter<edm::ParameterSet>("HFNoseSD").getParameter<double>("TimeSliceUnit")),
              p.getParameter<edm::ParameterSet>("HFNoseSD").getParameter<bool>("IgnoreTrackID")),
       hgcons_(nullptr),
       slopeMin_(0),
       levelT1_(99),
       levelT2_(99),
       tan30deg_(std::tan(30.0 * CLHEP::deg)) {
   numberingScheme_.reset(nullptr);
   mouseBite_.reset(nullptr);
  
   edm::ParameterSet m_HFN = p.getParameter<edm::ParameterSet>("HFNoseSD");
   eminHit_ = m_HFN.getParameter<double>("EminHit") * CLHEP::MeV;
   fiducialCut_ = m_HFN.getParameter<bool>("FiducialCut");
   distanceFromEdge_ = m_HFN.getParameter<double>("DistanceFromEdge");
   storeAllG4Hits_ = m_HFN.getParameter<bool>("StoreAllG4Hits");
   rejectMB_ = m_HFN.getParameter<bool>("RejectMouseBite");
   waferRot_ = m_HFN.getParameter<bool>("RotatedWafer");
   cornerMinMask_ = m_HFN.getParameter<int>("CornerMinMask");
   angles_ = m_HFN.getUntrackedParameter<std::vector<double>>("WaferAngles");
  
   nameX_ = ((name.find("HFNoseHits") != std::string::npos) ? "HGCalHFNoseSensitive" : "HFNoseSensitive");
  
   if (storeAllG4Hits_) {
     setUseMap(false);
     setNumberCheckedHits(0);
   }
  
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HFNSim") << "**************************************************"
                              << "\n"
                              << "*                                                *"
                              << "\n"
                              << "* Constructing a HFNoseSD  with name " << name << "\n"
                              << "*                                                *"
                              << "\n"
                              << "**************************************************";
 #endif
   edm::LogVerbatim("HFNSim") << "HFNoseSD:: Threshold for storing hits: " << eminHit_ << " for " << name;
   edm::LogVerbatim("HFNSim") << "Flag for storing individual Geant4 Hits " << storeAllG4Hits_;
   edm::LogVerbatim("HFNSim") << "Fiducial volume cut with cut from eta/phi "
                              << "boundary " << fiducialCut_ << " at " << distanceFromEdge_;
   edm::LogVerbatim("HFNSim") << "Reject MosueBite Flag: " << rejectMB_ << " cuts along " << angles_.size()
                              << " axes: " << angles_[0] << ", " << angles_[1];
 }

References angles_, cornerMinMask_, distanceFromEdge_, eminHit_, fiducialCut_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), MeV, mouseBite_, Skims_PA_cff::name, nameX_, numberingScheme_, AlCaHLTBitMon_ParallelJobs::p, rejectMB_, CaloSD::setNumberCheckedHits(), CaloSD::setUseMap(), storeAllG4Hits_, and waferRot_.

◆ ~HFNoseSD()

HFNoseSD::~HFNoseSD ( )

overridedefault

Member Function Documentation

◆ filterHit()

bool HFNoseSD::filterHit	(	CaloG4Hit *	aHit,
		double	time
	)

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 205 of file HFNoseSD.cc.

                                                      {
   return ((time <= tmaxHit) && (aHit->getEnergyDeposit() > eminHit_));
 }

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

◆ getEnergyDeposit()

double HFNoseSD::getEnergyDeposit ( const G4Step * aStep )

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 85 of file HFNoseSD.cc.

                                                      {
   double r = aStep->GetPreStepPoint()->GetPosition().perp();
   double z = std::abs(aStep->GetPreStepPoint()->GetPosition().z());
 #ifdef EDM_ML_DEBUG
   G4int parCode = aStep->GetTrack()->GetDefinition()->GetPDGEncoding();
   G4String parName = aStep->GetTrack()->GetDefinition()->GetParticleName();
   G4LogicalVolume* lv = aStep->GetPreStepPoint()->GetPhysicalVolume()->GetLogicalVolume();
   edm::LogVerbatim("HFNSim") << "HFNoseSD: Hit from standard path from " << lv->GetName() << " for Track "
                              << aStep->GetTrack()->GetTrackID() << " (" << parCode << ":" << parName << ") R = " << r
                              << " Z = " << z << " slope = " << r / z << ":" << slopeMin_;
 #endif
   // Apply fiducial cut
   if (r < z * slopeMin_) {
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HFNSim") << "HFNoseSD: Fiducial Volume cut";
 #endif
     return 0.0;
   }
  
   double wt1 = getResponseWt(aStep->GetTrack());
   double wt2 = aStep->GetTrack()->GetWeight();
   double destep = weight_ * wt1 * (aStep->GetTotalEnergyDeposit());
   if (wt2 > 0)
     destep *= wt2;
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HFNSim") << "HFNoseSD: weights= " << weight_ << ":" << wt1 << ":" << wt2 << " Total weight "
                              << weight_ * wt1 * wt2 << " deStep: " << aStep->GetTotalEnergyDeposit() << ":" << destep;
 #endif
   return destep;
 }

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

◆ initRun()

void HFNoseSD::initRun ( )

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 203 of file HFNoseSD.cc.

203 {}

◆ isItinFidVolume()

bool HFNoseSD::isItinFidVolume ( const G4ThreeVector & pos )

private

Definition at line 218 of file HFNoseSD.cc.

                                                        {
   if (fiducialCut_) {
     return (hgcons_->distFromEdgeHex(pos.x(), pos.y(), pos.z()) > distanceFromEdge_);
   } else {
     return true;
   }
 }

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

◆ setDetUnitId() [1/2]

uint32_t HFNoseSD::setDetUnitId ( const G4Step * step )

overridevirtual

Implements CaloSD.

Definition at line 116 of file HFNoseSD.cc.

                                                    {
   const G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
   const G4VTouchable* touch = preStepPoint->GetTouchable();
  
   //determine the exact position in global coordinates in the mass geometry
   G4ThreeVector hitPoint = preStepPoint->GetPosition();
   float globalZ = touch->GetTranslation(0).z();
   int iz(globalZ > 0 ? 1 : -1);
  
   int layer, module, cell;
   if ((touch->GetHistoryDepth() == levelT1_) || (touch->GetHistoryDepth() == levelT2_)) {
     layer = touch->GetReplicaNumber(0);
     module = -1;
     cell = -1;
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HFNSim") << "DepthsTop: " << touch->GetHistoryDepth() << ":" << levelT1_ << ":" << levelT2_
                                << " name " << touch->GetVolume(0)->GetName() << " layer:module:cell " << layer << ":"
                                << module << ":" << cell;
 #endif
   } else {
     layer = touch->GetReplicaNumber(3);
     module = touch->GetReplicaNumber(2);
     cell = touch->GetReplicaNumber(1);
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HFNSim") << "DepthsInside: " << touch->GetHistoryDepth() << " name "
                                << touch->GetVolume(0)->GetName() << " layer:module:cell " << layer << ":" << module
                                << ":" << cell;
 #endif
   }
 #ifdef EDM_ML_DEBUG
   G4Material* mat = aStep->GetPreStepPoint()->GetMaterial();
   edm::LogVerbatim("HFNSim") << "Depths: " << touch->GetHistoryDepth() << " name " << touch->GetVolume(0)->GetName()
                              << ":" << touch->GetReplicaNumber(0) << "   " << touch->GetVolume(1)->GetName() << ":"
                              << touch->GetReplicaNumber(1) << "   " << touch->GetVolume(2)->GetName() << ":"
                              << touch->GetReplicaNumber(2) << "   " << touch->GetVolume(3)->GetName() << ":"
                              << touch->GetReplicaNumber(3) << "   " << touch->GetVolume(4)->GetName() << ":"
                              << touch->GetReplicaNumber(4) << "   "
                              << " layer:module:cell " << layer << ":" << module << ":" << cell << " Material "
                              << mat->GetName() << ":" << mat->GetRadlen();
 #endif
   // The following statement should be examined later before elimination
   if (aStep->GetPreStepPoint()->GetMaterial()->GetRadlen() > 100000.)
     return 0;
  
   uint32_t id = setDetUnitId(layer, module, cell, iz, hitPoint);
   if (rejectMB_ && id != 0) {
     auto uv = HFNoseDetId(id).waferUV();
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HFNSim") << "ID " << std::hex << id << std::dec << " " << HFNoseDetId(id);
 #endif
     if (mouseBite_->exclude(hitPoint, iz, uv.first, uv.second)) {
       id = 0;
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HFNSim") << "Rejected by mousebite cutoff *****";
 #endif
     }
   }
   return id;
 }

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

◆ setDetUnitId() [2/2]

uint32_t HFNoseSD::setDetUnitId	(	int	layer,
		int	module,
		int	cell,
		int	iz,
		G4ThreeVector &	pos
	)

private

Definition at line 209 of file HFNoseSD.cc.

                                                                                            {
   uint32_t id = numberingScheme_ ? numberingScheme_->getUnitID(layer, module, cell, iz, pos, weight_) : 0;
   if (cornerMinMask_ > 2) {
     if (hgcons_->maskCell(DetId(id), cornerMinMask_))
       id = 0;
   }
   return id;
 }

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

◆ update() [1/6]

void CaloSD::update

overrideprotected

Definition at line 743 of file CaloSD.cc.

                                      {
   endEvent();
   slave.get()->ReserveMemory(theHC->entries());
  
   int count(0);
   int wrong(0);
   double eEM(0.0);
   double eHAD(0.0);
   double eEM2(0.0);
   double eHAD2(0.0);
   double tt(0.0);
   double zloc(0.0);
   double zglob(0.0);
   double ee(0.0);
   int hc_entries = theHC->entries();
   for (int i = 0; i < hc_entries; ++i) {
     if (!saveHit((*theHC)[i])) {
       ++wrong;
     }
     ++count;
     double x = (*theHC)[i]->getEM();
     eEM += x;
     eEM2 += x * x;
     x = (*theHC)[i]->getHadr();
     eHAD += x;
     eHAD2 += x * x;
     tt += (*theHC)[i]->getTimeSlice();
     ee += (*theHC)[i]->getIncidentEnergy();
     zglob += std::abs((*theHC)[i]->getEntry().z());
     zloc += std::abs((*theHC)[i]->getEntryLocal().z());
   }
  
   double norm = (count > 0) ? 1.0 / count : 0.0;
   eEM *= norm;
   eEM2 *= norm;
   eHAD *= norm;
   eHAD2 *= norm;
   eEM2 = std::sqrt(eEM2 - eEM * eEM);
   eHAD2 = std::sqrt(eHAD2 - eHAD * eHAD);
   tt *= norm;
   ee *= norm;
   zglob *= norm;
   zloc *= norm;
  
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CaloSim") << "CaloSD: " << GetName() << " store " << count << " hits; " << wrong
                               << " track IDs not given properly and " << totalHits - count
                               << " hits not passing cuts\n EmeanEM= " << eEM << " ErmsEM= " << eEM2
                               << "\n EmeanHAD= " << eHAD << " ErmsHAD= " << eHAD2 << " TimeMean= " << tt
                               << " E0mean= " << ee << " Zglob= " << zglob << " Zloc= " << zloc << " ";
 #endif
   tkMap.erase(tkMap.begin(), tkMap.end());
   std::vector<std::unique_ptr<CaloG4Hit>>().swap(reusehit);
   if (useMap)
     hitMap.erase(hitMap.begin(), hitMap.end());
 }

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

◆ update() [2/6]

void CaloSD::update

overrideprotected

Definition at line 709 of file CaloSD.cc.

                                                {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CaloSim") << "CaloSD: Dispatched BeginOfEvent for " << GetName() << " !";
 #endif
   clearHits();
   initEvent(g4Event);
 }

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

◆ update() [3/6]

void HFNoseSD::update ( const BeginOfJob * )

overrideprotectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob * >.

Definition at line 176 of file HFNoseSD.cc.

                                            {
   const edm::EventSetup* es = (*job)();
   edm::ESHandle<HGCalDDDConstants> hdc;
   es->get<IdealGeometryRecord>().get(nameX_, hdc);
   if (hdc.isValid()) {
     hgcons_ = hdc.product();
     geom_mode_ = hgcons_->geomMode();
     slopeMin_ = hgcons_->minSlope();
     levelT1_ = hgcons_->levelTop(0);
     levelT2_ = hgcons_->levelTop(1);
     double waferSize = hgcons_->waferSize(false);
     double mouseBite = hgcons_->mouseBite(false);
     mouseBiteCut_ = waferSize * tan30deg_ - mouseBite;
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HFNSim") << "HFNoseSD::Initialized with mode " << geom_mode_ << " Slope cut " << slopeMin_
                                << " top Level " << levelT1_ << ":" << levelT2_ << " wafer " << waferSize << ":"
                                << mouseBite;
 #endif
  
     numberingScheme_ = std::make_unique<HFNoseNumberingScheme>(*hgcons_);
     if (rejectMB_)
       mouseBite_ = std::make_unique<HGCMouseBite>(*hgcons_, angles_, mouseBiteCut_, waferRot_);
   } else {
     throw cms::Exception("Unknown", "HFNoseSD") << "Cannot find HGCalDDDConstants for " << nameX_ << "\n";
   }
 }

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

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

◆ update() [4/6]

void CaloSD::update

overrideprotected

Definition at line 707 of file CaloSD.cc.

707 { initRun(); }

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

◆ update() [5/6]

void CaloSD::update

overrideprotected

Definition at line 965 of file CaloSD.cc.

                                            {
   int primary = -1;
   TrackInformation* trkInfo = cmsTrackInformation((*trk)());
   if (trkInfo->isPrimary())
     primary = (*trk)()->GetTrackID();
  
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CaloSim") << "New track: isPrimary " << trkInfo->isPrimary() << " primary ID = " << primary
                               << " primary ancestor ID " << primAncestor;
 #endif
  
   // update the information if a different primary track ID
  
   if (primary > 0 && primary != primAncestor) {
     primAncestor = primary;
  
     // clean the hits information
  
     if (theHC->entries() > 0)
       cleanHitCollection();
   }
 }

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

◆ update() [6/6]

void CaloSD::update

overrideprotected

Definition at line 717 of file CaloSD.cc.

                                          {
   int id = (*trk)()->GetTrackID();
   TrackInformation* trkI = cmsTrackInformation((*trk)());
   int lastTrackID = -1;
   if (trkI)
     lastTrackID = trkI->getIDonCaloSurface();
   if (id == lastTrackID) {
     const TrackContainer* trksForThisEvent = m_trackManager->trackContainer();
     if (trksForThisEvent != nullptr) {
       int it = (int)(trksForThisEvent->size()) - 1;
       if (it >= 0) {
         TrackWithHistory* trkH = (*trksForThisEvent)[it];
         if (trkH->trackID() == (unsigned int)(id))
           tkMap[id] = trkH;
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("CaloSim") << "CaloSD: get track " << it << " from Container of size "
                                     << trksForThisEvent->size() << " with ID " << trkH->trackID();
       } else {
         edm::LogVerbatim("CaloSim") << "CaloSD: get track " << it << " from Container of size "
                                     << trksForThisEvent->size() << " with no ID";
 #endif
       }
     }
   }
 }

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

Member Data Documentation

◆ angles_

std::vector<double> HFNoseSD::angles_

private

Definition at line 53 of file HFNoseSD.h.

Referenced by HFNoseSD(), and update().

◆ cornerMinMask_

int HFNoseSD::cornerMinMask_

private

Definition at line 49 of file HFNoseSD.h.

Referenced by HFNoseSD(), and setDetUnitId().

◆ distanceFromEdge_

double HFNoseSD::distanceFromEdge_

private

Definition at line 48 of file HFNoseSD.h.

Referenced by HFNoseSD(), and isItinFidVolume().

◆ eminHit_

double HFNoseSD::eminHit_

private

Definition at line 47 of file HFNoseSD.h.

Referenced by filterHit(), and HFNoseSD().

◆ fiducialCut_

bool HFNoseSD::fiducialCut_

private

Definition at line 51 of file HFNoseSD.h.

Referenced by HFNoseSD(), and isItinFidVolume().

◆ geom_mode_

HGCalGeometryMode::GeometryMode HFNoseSD::geom_mode_

private

Definition at line 46 of file HFNoseSD.h.

Referenced by update().

◆ hgcons_

const HGCalDDDConstants* HFNoseSD::hgcons_

private

Definition at line 42 of file HFNoseSD.h.

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

◆ levelT1_

int HFNoseSD::levelT1_

private

Definition at line 49 of file HFNoseSD.h.

Referenced by setDetUnitId(), and update().

◆ levelT2_

int HFNoseSD::levelT2_

private

Definition at line 49 of file HFNoseSD.h.

Referenced by setDetUnitId(), and update().

◆ mouseBite_

std::unique_ptr<HGCMouseBite> HFNoseSD::mouseBite_

private

Definition at line 44 of file HFNoseSD.h.

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

◆ mouseBiteCut_

double HFNoseSD::mouseBiteCut_

private

Definition at line 48 of file HFNoseSD.h.

Referenced by update().

◆ nameX_

std::string HFNoseSD::nameX_

private

Definition at line 45 of file HFNoseSD.h.

Referenced by HFNoseSD(), and update().

◆ numberingScheme_

std::unique_ptr<HFNoseNumberingScheme> HFNoseSD::numberingScheme_

private

Definition at line 43 of file HFNoseSD.h.

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

◆ rejectMB_

bool HFNoseSD::rejectMB_

private

Definition at line 51 of file HFNoseSD.h.

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

◆ slopeMin_

double HFNoseSD::slopeMin_

private

Definition at line 47 of file HFNoseSD.h.

Referenced by getEnergyDeposit(), and update().

◆ storeAllG4Hits_

bool HFNoseSD::storeAllG4Hits_

private

Definition at line 50 of file HFNoseSD.h.

Referenced by HFNoseSD().

◆ tan30deg_

const double HFNoseSD::tan30deg_

private

Definition at line 52 of file HFNoseSD.h.

Referenced by update().

◆ waferRot_

bool HFNoseSD::waferRot_

private

Definition at line 51 of file HFNoseSD.h.

Referenced by HFNoseSD(), and update().

◆ weight_

double HFNoseSD::weight_

private

Definition at line 47 of file HFNoseSD.h.

Referenced by getEnergyDeposit(), and setDetUnitId().

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HFNoseSD()

◆ ~HFNoseSD()

Member Function Documentation

◆ filterHit()

◆ getEnergyDeposit()

◆ initRun()

◆ isItinFidVolume()

◆ setDetUnitId() [1/2]

◆ setDetUnitId() [2/2]

◆ update() [1/6]

◆ update() [2/6]

◆ update() [3/6]

◆ update() [4/6]

◆ update() [5/6]

◆ update() [6/6]

Member Data Documentation

◆ angles_

◆ cornerMinMask_

◆ distanceFromEdge_

◆ eminHit_

◆ fiducialCut_

◆ geom_mode_

◆ hgcons_

◆ levelT1_

◆ levelT2_

◆ mouseBite_

◆ mouseBiteCut_

◆ nameX_

◆ numberingScheme_

◆ rejectMB_

◆ slopeMin_

◆ storeAllG4Hits_

◆ tan30deg_

◆ waferRot_

◆ weight_