#include <SimG4CMS/Forward/interface/TotemSD.h>

Inheritance diagram for TotemSD:

Public Member Functions
void	clearHits () override

void	EndOfEvent (G4HCofThisEvent *eventHC) override

void	fillHits (edm::PSimHitContainer &, const std::string &) override

void	Initialize (G4HCofThisEvent *HCE) override

void	PrintAll () override

bool	ProcessHits (G4Step , G4TouchableHistory ) override

uint32_t	setDetUnitId (const G4Step *) override

	TotemSD (const std::string &, const SensitiveDetectorCatalog &, edm::ParameterSet const &, const SimTrackManager *)

	~TotemSD () override

Public Member Functions inherited from SensitiveTkDetector
	SensitiveTkDetector (const std::string &iname, const SensitiveDetectorCatalog &clg)

Public Member Functions inherited from SensitiveDetector
void	EndOfEvent (G4HCofThisEvent *eventHC) override

const std::vector< std::string > &	getNames () const

void	Initialize (G4HCofThisEvent *eventHC) override

bool	isCaloSD () const

	SensitiveDetector (const std::string &iname, const SensitiveDetectorCatalog &, bool calo, const std::string &newcollname="")

	~SensitiveDetector () override

Public Member Functions inherited from Observer< const BeginOfEvent *>
	Observer ()

void	slotForUpdate (const BeginOfEvent * iT)

virtual	~Observer ()

Protected Member Functions
void	update (const BeginOfEvent *) override
	This routine will be called when the appropriate signal arrives. More...

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

Private Member Functions
void	createNewHit ()

void	createNewHitEvo ()

void	getStepInfo (const G4Step *aStep)

bool	hitExists ()

G4ThreeVector	posizioEvo (const G4ThreeVector &, double, double, double, double, int &)

void	resetForNewPrimary ()

G4ThreeVector	setToLocal (const G4ThreeVector &globalPoint)

void	storeHit (TotemG4Hit *)

void	updateHit ()

Private Attributes
TotemG4Hit *	currentHit

G4VPhysicalVolume *	currentPV

float	edeposit

float	Eloss

G4ThreeVector	entrancePoint

G4int	hcID

G4ThreeVector	hitPoint

float	incidentEnergy

TotemVDetectorOrganization *	numberingScheme

float	Pabs

int	ParentId

short	ParticleType

float	PhiAtEntry

G4ThreeVector	Posizio

const G4StepPoint *	postStepPoint

const G4StepPoint *	preStepPoint

uint32_t	previousUnitID

int	primaryID

G4int	primID

TrackingSlaveSD *	slave

TotemG4HitCollection *	theHC

const SimTrackManager *	theManager

float	ThetaAtEntry

G4Track *	theTrack

float	Tof

int	tsID

double	tSlice

int	tSliceID

uint32_t	unitID

float	Vx

float	Vy

float	Vz

Additional Inherited Members
Protected Types inherited from SensitiveDetector
enum	coordinates { WorldCoordinates, LocalCoordinates }

Detailed Description

Description: Stores hits of Totem in appropriate container

Usage: Used in sensitive detector builder

Definition at line 46 of file TotemSD.h.

Constructor & Destructor Documentation

◆ TotemSD()

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

Definition at line 47 of file TotemSD.cc.

References edm::ParameterSet::getUntrackedParameter(), Skims_PA_cff::name, numberingScheme, AlCaHLTBitMon_ParallelJobs::p, and slave.

     : SensitiveTkDetector(name, clg),
       numberingScheme(nullptr),
       hcID(-1),
       theHC(nullptr),
       theManager(manager),
       currentHit(nullptr),
       theTrack(nullptr),
       currentPV(nullptr),
       unitID(0),
       previousUnitID(0),
       preStepPoint(nullptr),
       postStepPoint(nullptr) {
   //Parameters
   edm::ParameterSet m_p = p.getParameter<edm::ParameterSet>("TotemSD");
   int verbn = m_p.getUntrackedParameter<int>("Verbosity");
 
   SetVerboseLevel(verbn);
 
   slave = new TrackingSlaveSD(name);
 
   if (name == "TotemHitsT1") {
     numberingScheme = dynamic_cast<TotemVDetectorOrganization*>(new TotemT1NumberingScheme(1));
   } else if (name == "TotemHitsT2Si") {
     numberingScheme = dynamic_cast<TotemVDetectorOrganization*>(new TotemT2NumberingSchemeGem(3));
   } else if (name == "TotemHitsT2Gem") {
     numberingScheme = dynamic_cast<TotemVDetectorOrganization*>(new TotemT2NumberingSchemeGem(4));
   } else if (name == "TotemHitsRP") {
     numberingScheme = dynamic_cast<TotemVDetectorOrganization*>(new TotemRPNumberingScheme(3));
   } else {
     edm::LogWarning("ForwardSim") << "TotemSD: ReadoutName not supported\n";
   }
 }

◆ ~TotemSD()

TotemSD::~TotemSD ( )

override

Definition at line 84 of file TotemSD.cc.

References numberingScheme, and slave.

                   {
   delete slave;
   delete numberingScheme;
 }

Member Function Documentation

◆ clearHits()

void TotemSD::clearHits ( )

overridevirtual

Implements SensitiveDetector.

Definition at line 164 of file TotemSD.cc.

References TrackingSlaveSD::Initialize(), and slave.

Referenced by update().

164 { slave->Initialize(); }

TrackingSlaveSD::Initialize

virtual void Initialize()

Definition: TrackingSlaveSD.cc:17

TotemSD::slave

TrackingSlaveSD * slave

Definition: TotemSD.h:76

◆ createNewHit()

void TotemSD::createNewHit ( )

private

Definition at line 249 of file TotemSD.cc.

References currentHit, currentPV, Eloss, incidentEnergy, Pabs, ParentId, ParticleType, PhiAtEntry, Posizio, primaryID, TotemG4Hit::setEnergyLoss(), TotemG4Hit::setEntry(), TotemG4Hit::setIncidentEnergy(), TotemG4Hit::setPabs(), TotemG4Hit::setParentId(), TotemG4Hit::setParticleType(), TotemG4Hit::setPhiAtEntry(), TotemG4Hit::setThetaAtEntry(), TotemG4Hit::setTimeSlice(), TotemG4Hit::setTof(), TotemG4Hit::setTrackID(), TotemG4Hit::setUnitID(), TotemG4Hit::setVx(), TotemG4Hit::setVy(), TotemG4Hit::setVz(), storeHit(), ThetaAtEntry, theTrack, Tof, tSlice, tSliceID, unitID, updateHit(), Vx, Vy, and Vz.

Referenced by ProcessHits().

                            {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("ForwardSim") << "TotemSD CreateNewHit for PV " << currentPV->GetName()
                                  << " PVid = " << currentPV->GetCopyNo() << " Unit " << unitID << "\n primary "
                                  << primaryID << " time slice " << tSliceID << " For Track  " << theTrack->GetTrackID()
                                  << " which is a " << theTrack->GetDefinition()->GetParticleName();
 
   if (theTrack->GetTrackID() == 1) {
     edm::LogVerbatim("ForwardSim") << " of energy " << theTrack->GetTotalEnergy();
   } else {
     edm::LogVerbatim("ForwardSim") << " daughter of part. " << theTrack->GetParentID();
   }
 
   if (theTrack->GetCreatorProcess() != nullptr)
     edm::LogVerbatim("ForwardSim") << theTrack->GetCreatorProcess()->GetProcessName();
   else
     edm::LogVerbatim("ForwardSim") << "NO process";
 #endif
 
   currentHit = new TotemG4Hit;
   currentHit->setTrackID(primaryID);
   currentHit->setTimeSlice(tSlice);
   currentHit->setUnitID(unitID);
   currentHit->setIncidentEnergy(incidentEnergy);
 
   currentHit->setPabs(Pabs);
   currentHit->setTof(Tof);
   currentHit->setEnergyLoss(Eloss);
   currentHit->setParticleType(ParticleType);
   currentHit->setThetaAtEntry(ThetaAtEntry);
   currentHit->setPhiAtEntry(PhiAtEntry);
 
   currentHit->setEntry(Posizio.x(), Posizio.y(), Posizio.z());
 
   currentHit->setParentId(ParentId);
   currentHit->setVx(Vx);
   currentHit->setVy(Vy);
   currentHit->setVz(Vz);
 
   updateHit();
 
   storeHit(currentHit);
 }

◆ createNewHitEvo()

void TotemSD::createNewHitEvo ( )

private

Definition at line 293 of file TotemSD.cc.

References currentHit, Eloss, incidentEnergy, Pabs, ParentId, ParticleType, PhiAtEntry, Posizio, posizioEvo(), primaryID, TotemG4Hit::setEnergyLoss(), TotemG4Hit::setEntry(), TotemG4Hit::setIncidentEnergy(), TotemG4Hit::setPabs(), TotemG4Hit::setParentId(), TotemG4Hit::setParticleType(), TotemG4Hit::setPhiAtEntry(), TotemG4Hit::setThetaAtEntry(), TotemG4Hit::setTimeSlice(), TotemG4Hit::setTof(), TotemG4Hit::setTrackID(), TotemG4Hit::setUnitID(), TotemG4Hit::setVx(), TotemG4Hit::setVy(), TotemG4Hit::setVz(), storeHit(), ThetaAtEntry, Tof, tSlice, unitID, updateHit(), Vx, Vy, and Vz.

Referenced by ProcessHits().

                               {
   // edm::LogVerbatim("ForwardSim") << "INSIDE CREATE NEW HIT EVO ";
 
   currentHit = new TotemG4Hit;
   currentHit->setTrackID(primaryID);
   currentHit->setTimeSlice(tSlice);
   currentHit->setUnitID(unitID);
   currentHit->setIncidentEnergy(incidentEnergy);
 
   currentHit->setPabs(Pabs);
   currentHit->setTof(Tof);
   currentHit->setEnergyLoss(Eloss);
   currentHit->setParticleType(ParticleType);
   currentHit->setThetaAtEntry(ThetaAtEntry);
   currentHit->setPhiAtEntry(PhiAtEntry);
 
   //  edm::LogVerbatim("ForwardSim") << Posizio.x() << " " << Posizio.y() << " " << Posizio.z();
 
   currentHit->setParentId(ParentId);
   currentHit->setVx(Vx);
   currentHit->setVy(Vy);
   currentHit->setVz(Vz);
 
   G4ThreeVector _PosizioEvo;
   int flagAcc = 0;
   _PosizioEvo = posizioEvo(Posizio, Vx, Vy, Vz, Pabs, flagAcc);
 
   if (flagAcc == 1) {
     currentHit->setEntry(_PosizioEvo.x(), _PosizioEvo.y(), _PosizioEvo.z());
 
     updateHit();
 
     storeHit(currentHit);
   }
   // edm::LogVerbatim("ForwardSim") << "STORED HIT IN: " << unitID;
 }

◆ EndOfEvent()

void TotemSD::EndOfEvent ( G4HCofThisEvent * eventHC )

override

Definition at line 119 of file TotemSD.cc.

References TotemG4Hit::getEnergyLoss(), TotemG4Hit::getEntry(), TotemG4Hit::getPabs(), TotemG4Hit::getParticleType(), TotemG4Hit::getPhiAtEntry(), TotemG4Hit::getThetaAtEntry(), TotemG4Hit::getTof(), TotemG4Hit::getTrackID(), TotemG4Hit::getUnitID(), dqmiolumiharvest::j, TrackingSlaveSD::processHits(), slave, and theHC.

                                          {
   // here we loop over transient hits and make them persistent
   int thehc_entries = theHC->entries();
   for (int j = 0; j < thehc_entries && j < 15000; j++) {
     TotemG4Hit* aHit = (*theHC)[j];
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("ForwardSim") << "HIT NUMERO " << j << "unit ID = " << aHit->getUnitID()
                                    << "\n               enrty z " << aHit->getEntry().z() << "\n               theta   "
                                    << aHit->getThetaAtEntry() << "\n";
 #endif
     Local3DPoint theExitPoint(0, 0, 0);
     Local3DPoint Entrata(aHit->getEntry().x(), aHit->getEntry().y(), aHit->getEntry().z());
     slave->processHits(PSimHit(Entrata,
                                theExitPoint,
                                aHit->getPabs(),
                                aHit->getTof(),
                                aHit->getEnergyLoss(),
                                aHit->getParticleType(),
                                aHit->getUnitID(),
                                aHit->getTrackID(),
                                aHit->getThetaAtEntry(),
                                aHit->getPhiAtEntry()));
   }
 }

◆ fillHits()

void TotemSD::fillHits	(	edm::PSimHitContainer &	cc,
		const std::string &	hname
	)

overridevirtual

Implements SensitiveTkDetector.

Definition at line 151 of file TotemSD.cc.

References gpuPixelDoublets::cc, TrackingSlaveSD::hits(), TrackingSlaveSD::name(), and slave.

                                                                       {
   if (slave->name() == hname) {
     cc = slave->hits();
   }
 }

◆ getStepInfo()

void TotemSD::getStepInfo ( const G4Step * aStep )

private

Definition at line 173 of file TotemSD.cc.

References currentPV, edeposit, Eloss, hitPoint, createfilelist::int, Skims_PA_cff::name, Pabs, ParentId, PbPb_ZMuSkimMuonDPG_cff::particleType, ParticleType, PhiAtEntry, Posizio, postStepPoint, preStepPoint, primaryID, setDetUnitId(), ThetaAtEntry, theTrack, Tof, tSlice, tSliceID, unitID, Vx, Vy, and Vz.

Referenced by ProcessHits().

                                              {
   preStepPoint = aStep->GetPreStepPoint();
   postStepPoint = aStep->GetPostStepPoint();
   theTrack = aStep->GetTrack();
   hitPoint = preStepPoint->GetPosition();
   currentPV = preStepPoint->GetPhysicalVolume();
 
   // double weight = 1;
   G4String name = currentPV->GetName();
   name.assign(name, 0, 4);
   G4String particleType = theTrack->GetDefinition()->GetParticleName();
   edeposit = aStep->GetTotalEnergyDeposit();
 
   tSlice = (postStepPoint->GetGlobalTime()) / CLHEP::nanosecond;
   tSliceID = (int)tSlice;
   unitID = setDetUnitId(aStep);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("ForwardSim") << "UNITa " << unitID;
 #endif
   primaryID = theTrack->GetTrackID();
 
   Posizio = hitPoint;
   Pabs = aStep->GetPreStepPoint()->GetMomentum().mag() / CLHEP::GeV;
   Tof = aStep->GetPostStepPoint()->GetGlobalTime() / CLHEP::nanosecond;
 
   Eloss = aStep->GetTotalEnergyDeposit() / CLHEP::GeV;
   ParticleType = theTrack->GetDefinition()->GetPDGEncoding();
 
   ThetaAtEntry = aStep->GetPreStepPoint()->GetPosition().theta() / CLHEP::deg;
   PhiAtEntry = aStep->GetPreStepPoint()->GetPosition().phi() / CLHEP::deg;
 
   ParentId = theTrack->GetParentID();
   Vx = theTrack->GetVertexPosition().x();
   Vy = theTrack->GetVertexPosition().y();
   Vz = theTrack->GetVertexPosition().z();
 }

◆ hitExists()

bool TotemSD::hitExists ( )

private

Definition at line 210 of file TotemSD.cc.

References currentHit, newFWLiteAna::found, TotemG4Hit::getTimeSliceID(), TotemG4Hit::getTrackID(), TotemG4Hit::getUnitID(), dqmiolumiharvest::j, previousUnitID, primaryID, primID, resetForNewPrimary(), theHC, tsID, tSliceID, unitID, and updateHit().

Referenced by ProcessHits().

                         {
   if (primaryID < 1) {
     edm::LogWarning("ForwardSim") << "***** TotemSD error: primaryID = " << primaryID << " maybe detector name changed";
   }
 
   // Update if in the same detector, time-slice and for same track
   //  if (primaryID == primID && tSliceID == tsID && unitID==previousUnitID) {
   if (tSliceID == tsID && unitID == previousUnitID) {
     updateHit();
     return true;
   }
 
   // Reset entry point for new primary
   if (primaryID != primID)
     resetForNewPrimary();
 
   //look in the HitContainer whether a hit with the same primID, unitID,
   //tSliceID already exists:
 
   bool found = false;
   int thehc_entries = theHC->entries();
   for (int j = 0; j < thehc_entries && !found; j++) {
     TotemG4Hit* aPreviousHit = (*theHC)[j];
     if (aPreviousHit->getTrackID() == primaryID && aPreviousHit->getTimeSliceID() == tSliceID &&
         aPreviousHit->getUnitID() == unitID) {
       currentHit = aPreviousHit;
       found = true;
       break;
     }
   }
 
   if (found) {
     updateHit();
     return true;
   } else {
     return false;
   }
 }

◆ Initialize()

void TotemSD::Initialize ( G4HCofThisEvent * HCE )

override

Definition at line 105 of file TotemSD.cc.

References bysipixelclustmulteventfilter_cfi::collectionName, hcID, primID, theHC, and tsID.

                                              {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("ForwardSim") << "TotemSD : Initialize called for " << GetName();
 #endif
 
   theHC = new TotemG4HitCollection(GetName(), collectionName[0]);
   if (hcID < 0)
     hcID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]);
   HCE->AddHitsCollection(hcID, theHC);
 
   tsID = -2;
   primID = -2;
 }

◆ posizioEvo()

G4ThreeVector TotemSD::posizioEvo	(	const G4ThreeVector &	Pos,
		double	vx,
		double	vy,
		double	vz,
		double	pabs,
		int &	accettanza
	)

private

Definition at line 330 of file TotemSD.cc.

References funct::abs(), dqmiolumiharvest::j, mathSSE::sqrt(), btvMC_cff::vx, btvMC_cff::vy, and btvMC_cff::vz.

Referenced by createNewHitEvo().

                                                                                              {
   accettanza = 0;
   //Pos.xyz() in mm
   G4ThreeVector PosEvo;
   double ThetaX = std::atan((Pos.x() - vx) / (Pos.z() - vz));
   double ThetaY = std::atan((Pos.y() - vy) / (Pos.z() - vz));
   double X_at_0 = (vx - ((Pos.x() - vx) / (Pos.z() - vz)) * vz) / 1000.;
   double Y_at_0 = (vy - ((Pos.y() - vy) / (Pos.z() - vz)) * vz) / 1000.;
 
   //  double temp_evo_X;
   //  double temp_evo_Y;
   //  double temp_evo_Z;
   //  temp_evo_Z = std::abs(Pos.z())/Pos.z()*220000.;
 
   //csi=-dp/d
   double csi = std::abs((7000. - pabs) / 7000.);
 
   // all in m
   const int no_rp = 4;
   double x_par[no_rp + 1];
   double y_par[no_rp + 1];
   //rp z position
   double rp[no_rp] = {141., 149., 198., 220.};
   //{lx0,mlx} for each rp; Lx=lx0+mlx*csi
   double leffx[][2] = {{122.5429, -46.9312}, {125.4194, -49.1849}, {152.6, -81.157}, {98.8914, -131.8390}};
   //{ly0,mly} for each rp; Ly=ly0+mly*csi
   double leffy[][2] = {{124.2314, -55.4852}, {127.7825, -57.4503}, {179.455, -76.274}, {273.0931, -40.4626}};
   //{vx0,mvx0} for each rp; vx=vx0+mvx*csi
   double avx[][2] = {{0.515483, -1.0123}, {0.494122, -1.0534}, {0.2217, -1.483}, {0.004633, -1.0719}};
   //{vy0,mvy0} for each rp; vy=vy0+mvy*csi
   double avy[][2] = {{0.371418, -1.6327}, {0.349035, -1.6955}, {0.0815, -2.59}, {0.007592, -4.0841}};
   //{D0,md,a,b} for each rp; D=D0+(md+a*thetax)*csi+b*thetax
   double ddx[][4] = {{-0.082336, -0.092513, 112.3436, -82.5029},
                      {-0.086927, -0.097670, 114.9513, -82.9835},
                      {-0.092117, -0.0915, 180.6236, -82.443},
                      {-0.050470, 0.058837, 208.1106, 20.8198}};
   // {10sigma_x+0.5mm,10sigma_y+0.5mm}
   double detlim[][2] = {{0, 0}, {0.0028, 0.0021}, {0, 0}, {0.0008, 0.0013}};
   //{rmax,dmax}
   double pipelim[][2] = {{0.026, 0.026}, {0.04, 0.04}, {0.0226, 0.0177}, {0.04, 0.04}};
 
   for (int j = 0; j < no_rp; j++) {
     //y=Ly*thetay+vy*y0
     //x=Lx*thetax+vx*x0-csi*D
     y_par[j] = ThetaY * (leffy[j][0] + leffy[j][1] * csi) + (avy[j][0] + avy[j][1] * csi) * Y_at_0;
     x_par[j] = ThetaX * (leffx[j][0] + leffx[j][1] * csi) + (avx[j][0] + avx[j][1] * csi) * X_at_0 -
                csi * (ddx[j][0] + (ddx[j][1] + ddx[j][2] * ThetaX) * csi + ddx[j][3] * ThetaX);
   }
 
   //pass TAN@141
   if (std::abs(y_par[0]) < pipelim[0][1] && sqrt((y_par[0] * y_par[0]) + (x_par[0] * x_par[0])) < pipelim[0][0]) {
     //pass 149
     if ((sqrt((y_par[1] * y_par[1]) + (x_par[1] * x_par[1])) < pipelim[1][0]) &&
         (std::abs(y_par[1]) > detlim[1][1] || x_par[1] > detlim[1][0])) {
       accettanza = 1;
     }
   }
 
   //pass TAN@141
   if (std::abs(y_par[0]) < pipelim[0][1] && sqrt((y_par[0]) * (y_par[0]) + (x_par[0]) * (x_par[0])) < pipelim[0][0]) {
     //pass Q5@198
     if (std::abs(y_par[2]) < pipelim[2][1] && sqrt((y_par[2] * y_par[2]) + (x_par[2] * x_par[2])) < pipelim[2][0]) {
       //pass 220
       if ((sqrt((y_par[3] * y_par[3]) + (x_par[3] * x_par[3])) < pipelim[3][0]) &&
           (std::abs(y_par[3]) > detlim[3][1] || x_par[3] > detlim[3][0])) {
         accettanza = 1;
 
         PosEvo.setX(1000 * x_par[3]);
         PosEvo.setY(1000 * y_par[3]);
         PosEvo.setZ(1000 * rp[3]);
         if (Pos.z() < vz)
           PosEvo.setZ(-1000 * rp[3]);
       }
     }
   }
   /*
   edm::LogVerbatim("ForwardSim") << "\n"
              << "ACCETTANZA: "<<accettanza << "\n" 
              << "CSI: "<< csi << "\n"
              << "Theta_X: " << ThetaX << "\n"
              << "Theta_Y: " << ThetaY << "\n"
              << "X_at_0: "<< X_at_0 << "\n"
              << "Y_at_0: "<< Y_at_0 << "\n" 
              << "x_par[3]: "<< x_par[3] << "\n"
              << "y_par[3]: "<< y_par[3] << "\n"
              << "pos " << Pos.x() << " " << Pos.y() << " " 
              << Pos.z() << "\n" << "V "<< vx << " " << vy << " "
              << vz << "\n"
 */
   // --------------
   return PosEvo;
 }

◆ PrintAll()

void TotemSD::PrintAll ( )

override

Definition at line 144 of file TotemSD.cc.

References theHC.

                        {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("ForwardSim") << "TotemSD: Collection " << theHC->GetName();
 #endif
   theHC->PrintAllHits();
 }

◆ ProcessHits()

bool TotemSD::ProcessHits	(	G4Step *	aStep,
		G4TouchableHistory *
	)

overridevirtual

Implements SensitiveDetector.

Definition at line 89 of file TotemSD.cc.

References createNewHit(), createNewHitEvo(), edeposit, getStepInfo(), hitExists(), ParentId, ParticleType, and unitID.

                                                             {
   getStepInfo(aStep);
   if (!hitExists() && edeposit > 0.) {
     createNewHit();
     return true;
   }
   if (!hitExists() && (((unitID == 1111 || unitID == 2222) && ParentId == 0 && ParticleType == 2212))) {
     createNewHitEvo();
   }
   return true;
 }

◆ resetForNewPrimary()

void TotemSD::resetForNewPrimary ( )

private

Definition at line 452 of file TotemSD.cc.

References entrancePoint, hitPoint, incidentEnergy, preStepPoint, and setToLocal().

Referenced by hitExists().

                                  {
   entrancePoint = setToLocal(hitPoint);
   incidentEnergy = preStepPoint->GetKineticEnergy();
 }

◆ setDetUnitId()

uint32_t TotemSD::setDetUnitId ( const G4Step * aStep )

overridevirtual

Implements SensitiveDetector.

Definition at line 101 of file TotemSD.cc.

References TotemVDetectorOrganization::getUnitID(), and numberingScheme.

Referenced by getStepInfo().

                                                   {
   return (numberingScheme == nullptr ? 0 : numberingScheme->getUnitID(aStep));
 }

◆ setToLocal()

G4ThreeVector TotemSD::setToLocal ( const G4ThreeVector & globalPoint )

private

Definition at line 166 of file TotemSD.cc.

References preStepPoint.

Referenced by resetForNewPrimary().

                                                              {
   G4ThreeVector localPoint;
   const G4VTouchable* touch = preStepPoint->GetTouchable();
   localPoint = touch->GetHistory()->GetTopTransform().TransformPoint(global);
   return localPoint;
 }

◆ storeHit()

void TotemSD::storeHit ( TotemG4Hit * hit )

private

Definition at line 441 of file TotemSD.cc.

References primID, and theHC.

Referenced by createNewHit(), and createNewHitEvo().

                                       {
   if (primID < 0)
     return;
   if (hit == nullptr) {
     edm::LogWarning("ForwardSim") << "TotemSD: hit to be stored is NULL !!";
     return;
   }
 
   theHC->insert(hit);
 }

◆ update()

void TotemSD::update ( const BeginOfEvent * )

overrideprotectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent *>.

Definition at line 157 of file TotemSD.cc.

References clearHits().

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

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

◆ updateHit()

void TotemSD::updateHit ( )

private

Definition at line 424 of file TotemSD.cc.

References currentHit, Eloss, postStepPoint, previousUnitID, primaryID, primID, TotemG4Hit::setEnergyLoss(), tsID, tSliceID, and unitID.

Referenced by createNewHit(), createNewHitEvo(), and hitExists().

                         {
   //
   if (Eloss > 0.) {
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("ForwardSim") << "G4TotemT1SD updateHit: add eloss " << Eloss << "\nCurrentHit=" << currentHit
                                    << ", PostStepPoint=" << postStepPoint->GetPosition();
 #endif
 
     currentHit->setEnergyLoss(Eloss);
   }
 
   // buffer for next steps:
   tsID = tSliceID;
   primID = primaryID;
   previousUnitID = unitID;
 }