#include <HGCalSD.h>

Inheritance diagram for HGCalSD:

Public Member Functions
	HGCalSD (const std::string &, const HGCalDDDConstants , 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 SensitiveDetectorCatalog &clg, edm::ParameterSet const &p, const SimTrackManager *, float timeSlice=1., bool ignoreTkID=false, const std::string &newcolname="")

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	newCollection (const std::string &name, edm::ParameterSet const &p)

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, const std::string &newcollname="")

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, const std::string &newcollname="")

	~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	processSecondHit (const G4Step , const G4Track ) override

void	update (const BeginOfEvent *) override

void	update (const BeginOfRun *) override

void	update (const ::EndOfEvent *) override

void	update (const BeginOfTrack *trk) override

void	update (const EndOfTrack *trk) 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 (int k=0)

CaloG4Hit *	createNewHit (const G4Step , const G4Track , int k)

virtual void	endEvent ()

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 (int k=0)

double	getResponseWt (const G4Track *, int k=0)

virtual int	getTrackID (const G4Track *)

bool	hitExists (const G4Step *, int k)

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, int k=0)

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 *, int k)

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	calibCell (const uint32_t &id)

bool	isItinFidVolume (const G4ThreeVector &)

uint32_t	setDetUnitId (int, int, int, int, G4ThreeVector &)

Private Attributes
std::vector< double >	angles_

bool	calibCell_

std::vector< int >	calibCellFullHD_

std::vector< int >	calibCellFullLD_

std::vector< int >	calibCellPartHD_

std::vector< int >	calibCellPartLD_

double	calibCellRHD_

double	calibCellRLD_

bool	calibCells_

std::unique_ptr< HGCalCellOffset >	cellOffset_

bool	checkID_

int	cornerMinMask_

const double	cos30deg_

bool	dd4hep_

double	distanceFromEdge_

double	eminHit_

bool	fiducialCut_

double	fraction_

HGCalGeometryMode::GeometryMode	geom_mode_

std::unique_ptr< HGCGuardRing >	guardRing_

double	guardRingOffset_

std::unique_ptr< HGCGuardRingPartial >	guardRingPartial_

const HGCalDDDConstants *	hgcons_

int	levelT1_

int	levelT2_

std::string	missingFile_

std::unique_ptr< HGCMouseBite >	mouseBite_

double	mouseBiteCut_

DetId::Detector	mydet_

std::string	myName_

std::string	nameX_

std::unique_ptr< HGCalNumberingScheme >	numberingScheme_

edm::ParameterSet const &	ps_

bool	rejectMB_

double	slopeMin_

bool	storeAllG4Hits_

const double	tan30deg_

int	useSimWt_

int	verbose_

bool	waferRot_

double	waferSize_

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
std::string	collName_ [2]

CaloG4Hit *	currentHit [2]

CaloHitID	currentID [2]

std::string	detName_

float	edepositEM

float	edepositHAD

double	eminHit

double	energyCut

G4ThreeVector	entranceLocal

G4ThreeVector	entrancePoint

bool	forceSave

std::vector< std::string >	hcn_

float	incidentEnergy

double	kmaxIon

double	kmaxNeutron

double	kmaxProton

int	nHC_

G4ThreeVector	posGlobal

CaloHitID	previousID [2]

bool	suppressHeavy

double	tmaxHit

std::vector< int >	useResMap_

Detailed Description

Definition at line 22 of file HGCalSD.h.

Constructor & Destructor Documentation

◆ HGCalSD()

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

Definition at line 33 of file HGCalSD.cc.

References angles_, cellOffset_, checkID_, cornerMinMask_, dd4hep_, distanceFromEdge_, eminHit_, fiducialCut_, DetId::Forward, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), guardRing_, guardRingPartial_, DetId::HGCalEE, DetId::HGCalHSi, missingFile_, mouseBite_, mydet_, myName_, Skims_PA_cff::name, nameX_, CaloSD::nHC_, numberingScheme_, AlCaHLTBitMon_ParallelJobs::p, rejectMB_, CaloSD::setNumberCheckedHits(), CaloSD::setUseMap(), storeAllG4Hits_, AlCaHLTBitMon_QueryRunRegistry::string, verbose_, and waferRot_.

     : CaloSD(name,
              clg,
              p,
              manager,
              static_cast<float>(p.getParameter<edm::ParameterSet>("HGCSD").getParameter<double>("TimeSliceUnit")),
              p.getParameter<edm::ParameterSet>("HGCSD").getParameter<bool>("IgnoreTrackID"),
              ("Calibration" + name)),
       myName_(name),
       hgcons_(hgc),
       ps_(p),
       slopeMin_(0),
       levelT1_(99),
       levelT2_(99),
       useSimWt_(0),
       calibCells_(false),
       tan30deg_(std::tan(30.0 * CLHEP::deg)),
       cos30deg_(std::cos(30.0 * CLHEP::deg)) {
   numberingScheme_.reset(nullptr);
   guardRing_.reset(nullptr);
   guardRingPartial_.reset(nullptr);
   mouseBite_.reset(nullptr);
   cellOffset_.reset(nullptr);
 
   edm::ParameterSet m_HGC = p.getParameter<edm::ParameterSet>("HGCSD");
   eminHit_ = m_HGC.getParameter<double>("EminHit") * CLHEP::MeV;
   fiducialCut_ = m_HGC.getParameter<bool>("FiducialCut");
   storeAllG4Hits_ = m_HGC.getParameter<bool>("StoreAllG4Hits");
   rejectMB_ = m_HGC.getParameter<bool>("RejectMouseBite");
   waferRot_ = m_HGC.getParameter<bool>("RotatedWafer");
   cornerMinMask_ = m_HGC.getParameter<int>("CornerMinMask");
   nHC_ = m_HGC.getParameter<int>("HitCollection");
   angles_ = m_HGC.getUntrackedParameter<std::vector<double>>("WaferAngles");
   missingFile_ = m_HGC.getUntrackedParameter<std::string>("MissingWaferFile");
   checkID_ = m_HGC.getUntrackedParameter<bool>("CheckID");
   verbose_ = m_HGC.getUntrackedParameter<int>("Verbosity");
   dd4hep_ = p.getParameter<bool>("g4GeometryDD4hepSource");
 
   if (storeAllG4Hits_) {
     setUseMap(false);
     setNumberCheckedHits(0);
   }
 
   //this is defined in the hgcsens.xml
   mydet_ = DetId::Forward;
   nameX_ = "HGCal";
   if (myName_.find("HitsEE") != std::string::npos) {
     mydet_ = DetId::HGCalEE;
     nameX_ = "HGCalEESensitive";
   } else if (myName_.find("HitsHEfront") != std::string::npos) {
     mydet_ = DetId::HGCalHSi;
     nameX_ = "HGCalHESiliconSensitive";
   }
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCSim") << "**************************************************"
                              << "\n"
                              << "*                                                *"
                              << "\n"
                              << "* Constructing a HGCalSD  with name " << name << "\n"
                              << "*                                                *"
                              << "\n"
                              << "**************************************************";
 #endif
   edm::LogVerbatim("HGCSim") << "HGCalSD:: Threshold for storing hits: " << eminHit_ << " for " << nameX_
                              << " detector " << mydet_ << " with " << nHC_ << " hit collections";
   edm::LogVerbatim("HGCSim") << "Flag for storing individual Geant4 Hits " << storeAllG4Hits_;
   edm::LogVerbatim("HGCSim") << "Fiducial volume cut with cut from eta/phi "
                              << "boundary " << fiducialCut_ << " at " << distanceFromEdge_;
   edm::LogVerbatim("HGCSim") << "Reject MosueBite Flag: " << rejectMB_ << " cuts along " << angles_.size()
                              << " axes: " << angles_[0] << ", " << angles_[1];
 }

◆ ~HGCalSD()

HGCalSD::~HGCalSD ( )

overridedefault

Member Function Documentation

◆ calibCell()

bool HGCalSD::calibCell ( const uint32_t & id )

private

Definition at line 377 of file HGCalSD.cc.

References calibCellFullHD_, calibCellFullLD_, calibCellPartHD_, calibCellPartLD_, calibCellRHD_, calibCellRLD_, cellOffset_, HGCalCell::cellPlacementIndex(), spr::find(), RemoveAddSevLevel::flag, fraction_, full, hgcons_, info(), nano_mu_digi_cff::layer, HGCalTypes::layerFrontBack(), HGCalDDDConstants::layerType(), M_PI, EgammaValidation_cff::num, HGCSiliconDetId::unpack(), HGCalTypes::waferFull(), HGCalTypes::waferHD(), HGCalDDDConstants::waferInfo(), HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), and ecaldqm::zside().

Referenced by setDetUnitId().

                                           {
   bool flag(false);
   int type, zside, layer, waferU, waferV, cellU, cellV;
   HGCSiliconDetId(id).unpack(type, zside, layer, waferU, waferV, cellU, cellV);
   HGCalParameters::waferInfo info = hgcons_->waferInfo(layer, waferU, waferV);
   bool hd = HGCalTypes::waferHD(info.type);
   bool full = HGCalTypes::waferFull(info.part);
   int indx = 100 * cellU + cellV;
   if (hd) {
     if (full)
       flag = (std::find(calibCellFullHD_.begin(), calibCellFullHD_.end(), indx) != calibCellFullHD_.end());
     else
       flag = (std::find(calibCellPartHD_.begin(), calibCellPartHD_.end(), indx) != calibCellPartHD_.end());
   } else {
     if (full)
       flag = (std::find(calibCellFullLD_.begin(), calibCellFullLD_.end(), indx) != calibCellFullLD_.end());
     else
       flag = (std::find(calibCellPartLD_.begin(), calibCellPartLD_.end(), indx) != calibCellPartLD_.end());
   }
   if (flag) {
     int32_t place =
         HGCalCell::cellPlacementIndex(zside, HGCalTypes::layerFrontBack(hgcons_->layerType(layer)), info.orient);
     int32_t type = hd ? 0 : 1;
     double num = hd ? (M_PI * calibCellRHD_ * calibCellRHD_) : (M_PI * calibCellRLD_ * calibCellRLD_);
     double bot = cellOffset_->cellAreaUV(cellU, cellV, place, type, true);
     fraction_ = (bot > 0 && bot > num) ? (num / bot) : 1.0;
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCSim") << HGCSiliconDetId(id) << " CalibrationCell flag " << flag << " and fraction " << num
                                << ":" << bot << ":" << fraction_;
 #endif
   } else {
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCSim") << HGCSiliconDetId(id) << " CalibrationCell flag " << flag;
 #endif
   }
   return flag;
 }

◆ filterHit()

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

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 345 of file HGCalSD.cc.

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

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

◆ getEnergyDeposit()

double HGCalSD::getEnergyDeposit ( const G4Step * aStep )

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 110 of file HGCalSD.cc.

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

                                                     {
   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("HGCSim") << "HGCalSD: 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("HGCSim") << "HGCalSD: 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("HGCSim") << "HGCalSD: weights= " << weight_ << ":" << wt1 << ":" << wt2 << " Total weight "
                              << weight_ * wt1 * wt2 << " deStep: " << aStep->GetTotalEnergyDeposit() << ":" << destep;
 #endif
   return destep;
 }

◆ initRun()

void HGCalSD::initRun ( )

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 343 of file HGCalSD.cc.

343 {}

◆ isItinFidVolume()

bool HGCalSD::isItinFidVolume ( const G4ThreeVector & )

private

◆ processSecondHit()

void HGCalSD::processSecondHit	(	const G4Step *	aStep,
		const G4Track *	theTrack
	)

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 349 of file HGCalSD.cc.

References calibCell_, CaloSD::createNewHit(), CaloSD::currentHit, CaloSD::currentID, CaloSD::edepositEM, CaloSD::edepositHAD, fraction_, and CaloSD::hitExists().

                                                                            {
   if (calibCell_) {
     float edEM(edepositEM), edHad(edepositHAD);
     currentID[1] = currentID[0];
     edepositEM *= fraction_;
     edepositHAD *= fraction_;
     if (!hitExists(aStep, 1)) {
       currentHit[1] = createNewHit(aStep, theTrack, 1);
     }
     edepositEM = edEM;
     edepositHAD = edHad;
   }
 }

◆ setDetUnitId() [1/2]

uint32_t HGCalSD::setDetUnitId ( const G4Step * step )

overridevirtual

Implements CaloSD.

Definition at line 141 of file HGCalSD.cc.

References calibCell(), calibCell_, calibCells_, HGCalDDDConstants::cassetteShiftSilicon(), HGCSiliconDetId::cellU(), HGCSiliconDetId::cellV(), checkID_, ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), angle_units::operators::convertRadToDeg(), debug, TauDecayModes::dec, change_name::diff, PVValHelper::dx, PVValHelper::dy, fiducialCut_, fraction_, guardRing_, guardRingPartial_, hgcons_, l1ctLayer2EG_cff::id, HGCalDDDConstants::isValidHex8(), nano_mu_digi_cff::layer, HGCSiliconDetId::layer(), HGCalTypes::layerFrontBack(), HGCalDDDConstants::layerType(), levelT1_, levelT2_, DTRecHitClients_cfi::local, HGCalDDDConstants::locateCell(), mouseBite_, CaloSD::printDetectorLevels(), rejectMB_, mathSSE::sqrt(), useSimWt_, verbose_, HGCalDDDConstants::waferFileInfo(), HGCalWaferIndex::waferIndex(), HGCalDDDConstants::waferTypeRotation(), HGCSiliconDetId::waferU(), HGCSiliconDetId::waferUV(), HGCSiliconDetId::waferV(), geometryCSVtoXML::xx, geometryCSVtoXML::xy, and HGCSiliconDetId::zside().

                                                   {
   const G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
   const G4VTouchable* touch = preStepPoint->GetTouchable();
   fraction_ = 1.0;
   calibCell_ = false;
   int dn = touch->GetHistoryDepth();
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCSim") << "DepthsTop: " << dn << ":" << levelT1_ << ":" << levelT2_;
   printDetectorLevels(touch);
 #endif
   //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(0), moduleLev(-1), cell(-1);
   if (useSimWt_ > 0) {
     layer = touch->GetReplicaNumber(2);
     moduleLev = 1;
   } else if (touch->GetHistoryDepth() > levelT2_) {
     layer = touch->GetReplicaNumber(4);
     cell = touch->GetReplicaNumber(1);
     moduleLev = 3;
   } else {
     layer = touch->GetReplicaNumber(3);
     moduleLev = 2;
   }
   int module = touch->GetReplicaNumber(moduleLev);
   if (verbose_ && (cell == -1))
     edm::LogVerbatim("HGCSim") << "Top " << touch->GetVolume(0)->GetName() << " Module "
                                << touch->GetVolume(moduleLev)->GetName();
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCSim") << "DepthsTop: " << touch->GetHistoryDepth() << ":" << levelT1_ << ":" << levelT2_ << ":"
                              << useSimWt_ << " name " << touch->GetVolume(0)->GetName() << " layer:module:cell "
                              << layer << ":" << moduleLev << ":" << module << ":" << cell;
   printDetectorLevels(touch);
   G4Material* mat = aStep->GetPreStepPoint()->GetMaterial();
   edm::LogVerbatim("HGCSim") << "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_ || fiducialCut_) && id != 0) {
     auto uv = HGCSiliconDetId(id).waferUV();
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCSim") << "ID " << std::hex << id << std::dec << " " << HGCSiliconDetId(id);
 #endif
     G4ThreeVector local = (touch->GetHistory()->GetTransform(dn - moduleLev).TransformPoint(hitPoint));
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCSim") << "Global Point " << hitPoint << " Down0 "
                                << touch->GetHistory()->GetTransform(dn).TransformPoint(hitPoint) << " Down1 "
                                << touch->GetHistory()->GetTransform(dn - 1).TransformPoint(hitPoint) << " Down2 "
                                << touch->GetHistory()->GetTransform(dn - 2).TransformPoint(hitPoint) << " Down3 "
                                << touch->GetHistory()->GetTransform(dn - 3).TransformPoint(hitPoint) << " Local "
                                << local;
 #endif
     if (fiducialCut_) {
       int layertype = hgcons_->layerType(layer);
       int frontBack = HGCalTypes::layerFrontBack(layertype);
       if (guardRing_->exclude(local, iz, frontBack, layer, uv.first, uv.second) ||
           guardRingPartial_->exclude(local, iz, frontBack, layer, uv.first, uv.second)) {
         id = 0;
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HGCSim") << "Rejected by GuardRing cutoff *****";
 #endif
       }
     }
     if ((rejectMB_) && (mouseBite_->exclude(local, iz, layer, uv.first, uv.second))) {
       id = 0;
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HGCSim") << "Rejected by MouseBite cutoff *****";
 #endif
     }
   }
 #ifdef EDM_ML_DEBUG
   if (id != 0)
     edm::LogVerbatim("HGCSim") << HGCSiliconDetId(id);
 #endif
   if ((id != 0) && checkID_) {
     HGCSiliconDetId hid1(id);
     bool cshift = (hgcons_->cassetteShiftSilicon(hid1.zside(), hid1.layer(), hid1.waferU(), hid1.waferV()));
     std::string_view pid = (cshift ? "HGCSim" : "HGCalSim");
     bool debug = (verbose_ > 0) ? true : false;
     auto xy = hgcons_->locateCell(hid1, debug);
     double xx = (hid1.zside() > 0) ? xy.first : -xy.first;
     double dx = xx - (hitPoint.x() / CLHEP::cm);
     double dy = xy.second - (hitPoint.y() / CLHEP::cm);
     double diff = (dx * dx + dy * dy);
     constexpr double tol = 2.0 * 2.0;
     bool valid1 = hgcons_->isValidHex8(hid1.layer(), hid1.waferU(), hid1.waferV(), hid1.cellU(), hid1.cellV(), true);
     if ((diff > tol) || (!valid1))
       pid = "HGCalError";
     auto partn = hgcons_->waferTypeRotation(hid1.layer(), hid1.waferU(), hid1.waferV(), false, false);
     int indx = HGCalWaferIndex::waferIndex(layer, hid1.waferU(), hid1.waferV());
     edm::LogVerbatim(pid) << "CheckID " << HGCSiliconDetId(id) << " Layer:Module:Cell:ModuleLev " << layer << ":"
                           << module << ":" << cell << ":" << moduleLev << " SimWt:history " << useSimWt_ << ":"
                           << touch->GetHistoryDepth() << ":" << levelT1_ << ":" << levelT2_ << " input position: ("
                           << hitPoint.x() / CLHEP::cm << ", " << hitPoint.y() / CLHEP::cm << ":"
                           << convertRadToDeg(std::atan2(hitPoint.y(), hitPoint.x())) << "); position from ID (" << xx
                           << ", " << xy.second << ") distance " << dx << ":" << dy << ":" << std::sqrt(diff)
                           << " Valid " << valid1 << " Wafer type|rotation " << partn.first << ":" << partn.second
                           << " Part:Orient:Cassette " << std::get<1>(hgcons_->waferFileInfo(indx)) << ":"
                           << std::get<2>(hgcons_->waferFileInfo(indx)) << ":"
                           << std::get<3>(hgcons_->waferFileInfo(indx)) << " CassetteShift " << cshift;
     if ((diff > tol) || (!valid1)) {
       printDetectorLevels(touch);
       hgcons_->locateCell(hid1, true);
     }
   }
 
   if ((id != 0) && calibCells_)
     calibCell_ = calibCell(id);
 
   return id;
 }

◆ setDetUnitId() [2/2]

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

private

Definition at line 363 of file HGCalSD.cc.

References cornerMinMask_, hgcons_, l1ctLayer2EG_cff::id, CaloSD::ignoreRejection(), nano_mu_digi_cff::layer, HGCalDDDConstants::maskCell(), numberingScheme_, HGCalDDDConstants::waferHexagon8File(), and weight_.

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

◆ update() [1/6]

void CaloSD::update

overrideprotected

Definition at line 749 of file CaloSD.cc.

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

749 { initRun(); }

HGCalSD::initRun

void initRun() override

Definition: HGCalSD.cc:343

◆ update() [2/6]

void CaloSD::update

overrideprotected

Definition at line 780 of file CaloSD.cc.

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

                                      {
   endEvent();
   for (int k = 0; k < nHC_; ++k) {
     slave[k].get()->ReserveMemory(theHC[k]->entries());
 
     int count(0);
     double eEM(0.0);
     double eHAD(0.0);
     double eEM2(0.0);
     double eHAD2(0.0);
 #ifdef EDM_ML_DEBUG
     int wrong(0);
     double tt(0.0);
     double zloc(0.0);
     double zglob(0.0);
     double ee(0.0);
 #endif
     int hc_entries = theHC[k]->entries();
     for (int i = 0; i < hc_entries; ++i) {
 #ifdef EDM_ML_DEBUG
       if (!saveHit((*theHC[k])[i], k)) {
         ++wrong;
       }
 #else
       saveHit((*theHC[k])[i], k);
 #endif
 
       ++count;
       double x = (*theHC[k])[i]->getEM();
       eEM += x;
       eEM2 += x * x;
       x = (*theHC[k])[i]->getHadr();
       eHAD += x;
       eHAD2 += x * x;
 #ifdef EDM_ML_DEBUG
       tt += (*theHC[k])[i]->getTimeSlice();
       ee += (*theHC[k])[i]->getIncidentEnergy();
       zglob += std::abs((*theHC[k])[i]->getEntry().z());
       zloc += std::abs((*theHC[k])[i]->getEntryLocal().z());
 #endif
     }
 
     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);
 #ifdef EDM_ML_DEBUG
     tt *= norm;
     ee *= norm;
     zglob *= norm;
     zloc *= norm;
     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
     std::vector<std::unique_ptr<CaloG4Hit>>().swap(reusehit[k]);
     if (useMap)
       hitMap[k].erase(hitMap[k].begin(), hitMap[k].end());
   }
   tkMap.erase(tkMap.begin(), tkMap.end());
   boundaryCrossingParentMap_.clear();
 }

◆ update() [3/6]

void CaloSD::update

overrideprotected

Definition at line 759 of file CaloSD.cc.

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

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

◆ update() [4/6]

void CaloSD::update

overrideprotected

Definition at line 751 of file CaloSD.cc.

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

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

◆ update() [5/6]

void CaloSD::update

overrideprotected

Definition at line 1013 of file CaloSD.cc.

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

                                            {
   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
 
     bool clean(false);
     for (int k = 0; k < nHC_; ++k)
       if (theHC[k]->entries() > 0)
         clean = true;
     if (clean)
       cleanHitCollection();
   }
 }

◆ update() [6/6]

void HGCalSD::update ( const BeginOfJob * )

overrideprotectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob *>.

Definition at line 267 of file HGCalSD.cc.

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

                                           {
   if (hgcons_ != nullptr) {
     geom_mode_ = hgcons_->geomMode();
     slopeMin_ = hgcons_->minSlope();
     levelT1_ = hgcons_->levelTop(0);
     levelT2_ = hgcons_->levelTop(1);
     if (dd4hep_) {
       ++levelT1_;
       ++levelT2_;
     }
     useSimWt_ = hgcons_->getParameter()->useSimWt_;
     int useOffset = hgcons_->getParameter()->useOffset_;
     waferSize_ = hgcons_->waferSize(false);
     double mouseBite = hgcons_->mouseBite(false);
     guardRingOffset_ = hgcons_->guardRingOffset(false);
     double sensorSizeOffset = hgcons_->sensorSizeOffset(false);
     if (useOffset > 0) {
       rejectMB_ = true;
       fiducialCut_ = true;
     }
     double mouseBiteNew = (fiducialCut_) ? (mouseBite + guardRingOffset_ + sensorSizeOffset / cos30deg_) : mouseBite;
     mouseBiteCut_ = waferSize_ * tan30deg_ - mouseBiteNew;
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCSim") << "HGCalSD::Initialized with mode " << geom_mode_ << " Slope cut " << slopeMin_
                                << " top Level " << levelT1_ << ":" << levelT2_ << " useSimWt " << useSimWt_ << " wafer "
                                << waferSize_ << ":" << mouseBite << ":" << guardRingOffset_ << ":" << sensorSizeOffset
                                << ":" << mouseBiteNew << ":" << mouseBiteCut_ << " useOffset " << useOffset
                                << " dd4hep " << dd4hep_;
 #endif
 
     numberingScheme_ = std::make_unique<HGCalNumberingScheme>(*hgcons_, mydet_, nameX_, missingFile_);
     if (rejectMB_)
       mouseBite_ = std::make_unique<HGCMouseBite>(*hgcons_, angles_, mouseBiteCut_, waferRot_);
     if (fiducialCut_) {
       guardRing_ = std::make_unique<HGCGuardRing>(*hgcons_);
       guardRingPartial_ = std::make_unique<HGCGuardRingPartial>(*hgcons_);
     }
 
     //Now for calibration cells
     calibCellRHD_ = hgcons_->calibCellRad(true);
     calibCellRLD_ = hgcons_->calibCellRad(false);
     calibCellFullHD_ = hgcons_->calibCells(true, true);
     calibCellPartHD_ = hgcons_->calibCells(true, false);
     calibCellFullLD_ = hgcons_->calibCells(false, true);
     calibCellPartLD_ = hgcons_->calibCells(false, false);
     calibCells_ = ((!calibCellFullHD_.empty()) || (!calibCellPartHD_.empty()));
     calibCells_ |= ((!calibCellFullLD_.empty()) || (!calibCellPartLD_.empty()));
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCSim") << "HGCalSD::Calibration cells initialized with  flag " << calibCells_;
     edm::LogVerbatim("HGCSim") << " Radii " << calibCellRHD_ << ":" << calibCellRLD_;
     std::ostringstream st1;
     for (const auto& cell : calibCellFullHD_)
       st1 << " " << cell;
     edm::LogVerbatim("HGCSim") << calibCellFullHD_.size() << " cells for High Density full wafers: " << st1.str();
     std::ostringstream st2;
     for (const auto& cell : calibCellPartHD_)
       st2 << " " << cell;
     edm::LogVerbatim("HGCSim") << calibCellPartHD_.size() << " cells for High Density partial wafers: " << st2.str();
     std::ostringstream st3;
     for (const auto& cell : calibCellFullLD_)
       st3 << " " << cell;
     edm::LogVerbatim("HGCSim") << calibCellFullLD_.size() << " cells for Low Density full wafers: " << st3.str();
     std::ostringstream st4;
     for (const auto& cell : calibCellPartLD_)
       st4 << " " << cell;
     edm::LogVerbatim("HGCSim") << calibCellPartLD_.size() << " cells for Low Density partial wafers: " << st4.str();
 #endif
   } else {
     throw cms::Exception("Unknown", "HGCalSD") << "Cannot find HGCalDDDConstants for " << nameX_ << "\n";
   }
   if ((nHC_ > 1) && calibCells_)
     newCollection(collName_[1], ps_);
   cellOffset_ = std::make_unique<HGCalCellOffset>(
       waferSize_, hgcons_->getUVMax(0), hgcons_->getUVMax(1), guardRingOffset_, mouseBiteCut_);
 }

Member Data Documentation

◆ angles_

std::vector<double> HGCalSD::angles_

private

Definition at line 68 of file HGCalSD.h.

Referenced by HGCalSD(), and update().

◆ calibCell_

bool HGCalSD::calibCell_

private

Definition at line 70 of file HGCalSD.h.

Referenced by processSecondHit(), and setDetUnitId().

◆ calibCellFullHD_

std::vector<int> HGCalSD::calibCellFullHD_

private

Definition at line 65 of file HGCalSD.h.

Referenced by calibCell(), and update().

◆ calibCellFullLD_

std::vector<int> HGCalSD::calibCellFullLD_

private

Definition at line 66 of file HGCalSD.h.

Referenced by calibCell(), and update().

◆ calibCellPartHD_

std::vector<int> HGCalSD::calibCellPartHD_

private

Definition at line 65 of file HGCalSD.h.

Referenced by calibCell(), and update().

◆ calibCellPartLD_

std::vector<int> HGCalSD::calibCellPartLD_

private

Definition at line 66 of file HGCalSD.h.

Referenced by calibCell(), and update().

◆ calibCellRHD_

double HGCalSD::calibCellRHD_

private

Definition at line 64 of file HGCalSD.h.

Referenced by calibCell(), and update().

◆ calibCellRLD_

double HGCalSD::calibCellRLD_

private

Definition at line 64 of file HGCalSD.h.

Referenced by calibCell(), and update().

◆ calibCells_

bool HGCalSD::calibCells_

private

Definition at line 63 of file HGCalSD.h.

Referenced by setDetUnitId(), and update().

◆ cellOffset_

std::unique_ptr<HGCalCellOffset> HGCalSD::cellOffset_

private

Definition at line 53 of file HGCalSD.h.

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

◆ checkID_

bool HGCalSD::checkID_

private

Definition at line 61 of file HGCalSD.h.

Referenced by HGCalSD(), and setDetUnitId().

◆ cornerMinMask_

int HGCalSD::cornerMinMask_

private

Definition at line 59 of file HGCalSD.h.

Referenced by HGCalSD(), and setDetUnitId().

◆ cos30deg_

const double HGCalSD::cos30deg_

private

Definition at line 67 of file HGCalSD.h.

Referenced by update().

◆ dd4hep_

bool HGCalSD::dd4hep_

private

Definition at line 63 of file HGCalSD.h.

Referenced by HGCalSD(), and update().

◆ distanceFromEdge_

double HGCalSD::distanceFromEdge_

private

Definition at line 57 of file HGCalSD.h.

Referenced by HGCalSD().

◆ eminHit_

double HGCalSD::eminHit_

private

Definition at line 57 of file HGCalSD.h.

Referenced by filterHit(), and HGCalSD().

◆ fiducialCut_

bool HGCalSD::fiducialCut_

private

Definition at line 61 of file HGCalSD.h.

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

◆ fraction_

double HGCalSD::fraction_

private

Definition at line 71 of file HGCalSD.h.

Referenced by calibCell(), processSecondHit(), and setDetUnitId().

◆ geom_mode_

HGCalGeometryMode::GeometryMode HGCalSD::geom_mode_

private

Definition at line 56 of file HGCalSD.h.

Referenced by update().

◆ guardRing_

std::unique_ptr<HGCGuardRing> HGCalSD::guardRing_

private

Definition at line 50 of file HGCalSD.h.

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

◆ guardRingOffset_

double HGCalSD::guardRingOffset_

private

Definition at line 58 of file HGCalSD.h.

Referenced by update().

◆ guardRingPartial_

std::unique_ptr<HGCGuardRingPartial> HGCalSD::guardRingPartial_

private

Definition at line 51 of file HGCalSD.h.

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

◆ hgcons_

const HGCalDDDConstants* HGCalSD::hgcons_

private

Definition at line 47 of file HGCalSD.h.

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

◆ levelT1_

int HGCalSD::levelT1_

private

Definition at line 59 of file HGCalSD.h.

Referenced by setDetUnitId(), and update().

◆ levelT2_

int HGCalSD::levelT2_

private

Definition at line 59 of file HGCalSD.h.

Referenced by setDetUnitId(), and update().

◆ missingFile_

std::string HGCalSD::missingFile_

private

Definition at line 69 of file HGCalSD.h.

Referenced by HGCalSD(), and update().

◆ mouseBite_

std::unique_ptr<HGCMouseBite> HGCalSD::mouseBite_

private

Definition at line 52 of file HGCalSD.h.

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

◆ mouseBiteCut_

double HGCalSD::mouseBiteCut_

private

Definition at line 58 of file HGCalSD.h.

Referenced by update().

◆ mydet_

DetId::Detector HGCalSD::mydet_

private

Definition at line 54 of file HGCalSD.h.

Referenced by HGCalSD(), and update().

◆ myName_

std::string HGCalSD::myName_

private

Definition at line 46 of file HGCalSD.h.

Referenced by HGCalSD().

◆ nameX_

std::string HGCalSD::nameX_

private

Definition at line 55 of file HGCalSD.h.

Referenced by HGCalSD(), and update().

◆ numberingScheme_

std::unique_ptr<HGCalNumberingScheme> HGCalSD::numberingScheme_

private

Definition at line 49 of file HGCalSD.h.

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

◆ ps_

edm::ParameterSet const& HGCalSD::ps_

private

Definition at line 48 of file HGCalSD.h.

Referenced by update().

◆ rejectMB_

bool HGCalSD::rejectMB_

private

Definition at line 61 of file HGCalSD.h.

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

◆ slopeMin_

double HGCalSD::slopeMin_

private

Definition at line 57 of file HGCalSD.h.

Referenced by getEnergyDeposit(), and update().

◆ storeAllG4Hits_

bool HGCalSD::storeAllG4Hits_

private

Definition at line 60 of file HGCalSD.h.

Referenced by HGCalSD().

◆ tan30deg_

const double HGCalSD::tan30deg_

private

Definition at line 67 of file HGCalSD.h.

Referenced by update().

◆ useSimWt_

int HGCalSD::useSimWt_

private

Definition at line 62 of file HGCalSD.h.

Referenced by setDetUnitId(), and update().

◆ verbose_

int HGCalSD::verbose_

private

Definition at line 62 of file HGCalSD.h.

Referenced by HGCalSD(), and setDetUnitId().

◆ waferRot_

bool HGCalSD::waferRot_

private

Definition at line 61 of file HGCalSD.h.

Referenced by HGCalSD(), and update().

◆ waferSize_

double HGCalSD::waferSize_

private

Definition at line 58 of file HGCalSD.h.

Referenced by update().

◆ weight_

double HGCalSD::weight_

private

Definition at line 58 of file HGCalSD.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

◆ HGCalSD()

◆ ~HGCalSD()

Member Function Documentation

◆ calibCell()

◆ filterHit()

◆ getEnergyDeposit()

◆ initRun()

◆ isItinFidVolume()

◆ processSecondHit()

◆ 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_

◆ calibCell_

◆ calibCellFullHD_

◆ calibCellFullLD_

◆ calibCellPartHD_

◆ calibCellPartLD_

◆ calibCellRHD_

◆ calibCellRLD_

◆ calibCells_

◆ cellOffset_

◆ checkID_

◆ cornerMinMask_

◆ cos30deg_

◆ dd4hep_

◆ distanceFromEdge_

◆ eminHit_

◆ fiducialCut_

◆ fraction_

◆ geom_mode_

◆ guardRing_

◆ guardRingOffset_

◆ guardRingPartial_

◆ hgcons_

◆ levelT1_

◆ levelT2_

◆ missingFile_

◆ mouseBite_

◆ mouseBiteCut_

◆ mydet_

◆ myName_

◆ nameX_

◆ numberingScheme_

◆ ps_

◆ rejectMB_

◆ slopeMin_

◆ storeAllG4Hits_

◆ tan30deg_

◆ useSimWt_

◆ verbose_

◆ waferRot_

◆ waferSize_

◆ weight_