#include <FP420SD.h>

Inheritance diagram for FP420SD:

Public Member Functions
void	clear () override

void	clearHits () override

void	DrawAll () override

void	EndOfEvent (G4HCofThisEvent *eventHC) override

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

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

virtual double	getEnergyDeposit (G4Step *step)

void	Initialize (G4HCofThisEvent *HCE) override

void	PrintAll () override

bool	ProcessHits (G4Step , G4TouchableHistory ) override

uint32_t	setDetUnitId (const G4Step *) override

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

	~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 EndOfEvent *>
	Observer ()

void	slotForUpdate (const EndOfEvent * iT)

virtual	~Observer ()

Protected Attributes
float	edepositEM

float	edepositHAD

G4int	emPDG

G4int	epPDG

G4int	gammaPDG

Private Member Functions
void	CreateNewHit ()

void	GetStepInfo (G4Step *aStep)

G4bool	HitExists ()

void	ResetForNewPrimary ()

G4ThreeVector	SetToLocal (const G4ThreeVector &global)

G4ThreeVector	SetToLocalExit (const G4ThreeVector &globalPoint)

void	StoreHit (FP420G4Hit *)

void	Summarize ()

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

void	UpdateHit ()

Private Attributes
FP420G4Hit *	currentHit

G4VPhysicalVolume *	currentPV

float	edeposit

float	Eloss

G4ThreeVector	entrancePoint

int	eventno

G4ThreeVector	exitPoint

G4int	hcID

G4ThreeVector	hitPoint

G4ThreeVector	hitPointExit

G4ThreeVector	hitPointLocal

G4ThreeVector	hitPointLocalExit

float	incidentEnergy

FP420NumberingScheme *	numberingScheme

float	Pabs

int	ParentId

short	ParticleType

float	PhiAtEntry

G4StepPoint *	postStepPoint

G4StepPoint *	preStepPoint

uint32_t	previousUnitID

unsigned int	primaryID

unsigned int	primID

TrackingSlaveSD *	slave

G4ThreeVector	theEntryPoint

G4ThreeVector	theExitPoint

FP420G4HitCollection *	theHC

const SimTrackManager *	theManager

float	ThetaAtEntry

G4Track *	theTrack

float	Tof

G4int	tsID

G4double	tSlice

G4int	tSliceID

uint32_t	unitID

float	Vx

float	Vy

float	Vz

float	X

float	Y

float	Z

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

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...

Detailed Description

Definition at line 44 of file FP420SD.h.

Constructor & Destructor Documentation

◆ FP420SD()

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

Definition at line 46 of file FP420SD.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),
       eventno(0) {
   //Parameters
   edm::ParameterSet m_p = p.getParameter<edm::ParameterSet>("FP420SD");
   int verbn = m_p.getUntrackedParameter<int>("Verbosity");
   //int verbn = 1;
 
   SetVerboseLevel(verbn);
 
   slave = new TrackingSlaveSD(name);
 
   if (name == "FP420SI") {
     if (verbn > 0) {
       edm::LogInfo("FP420Sim") << "name = FP420SI and  new FP420NumberingSchem";
     }
     numberingScheme = new FP420NumberingScheme();
   } else {
     edm::LogWarning("FP420Sim") << "FP420SD: ReadoutName not supported\n";
   }
 }

◆ ~FP420SD()

FP420SD::~FP420SD ( )

override

Definition at line 82 of file FP420SD.cc.

References numberingScheme, and slave.

                   {
   delete slave;
   delete numberingScheme;
 }

Member Function Documentation

◆ clear()

void FP420SD::clear ( void )

override

Definition at line 402 of file FP420SD.cc.

402 {}

◆ clearHits()

void FP420SD::clearHits ( )

overridevirtual

Implements SensitiveDetector.

Definition at line 441 of file FP420SD.cc.

References TrackingSlaveSD::Initialize(), and slave.

Referenced by update().

                         {
   //AZ:
   slave->Initialize();
 }

◆ CreateNewHit()

void FP420SD::CreateNewHit ( )

private

Definition at line 233 of file FP420SD.cc.

References currentHit, currentPV, Eloss, hitPoint, hitPointLocal, hitPointLocalExit, incidentEnergy, LogDebug, NULL, Pabs, ParentId, ParticleType, PhiAtEntry, previousUnitID, primaryID, primID, FP420G4Hit::setEnergyLoss(), FP420G4Hit::setEntry(), FP420G4Hit::setEntryLocalP(), FP420G4Hit::setExitLocalP(), FP420G4Hit::setIncidentEnergy(), FP420G4Hit::setPabs(), FP420G4Hit::setParentId(), FP420G4Hit::setParticleType(), FP420G4Hit::setPhiAtEntry(), FP420G4Hit::setThetaAtEntry(), FP420G4Hit::setTimeSlice(), FP420G4Hit::setTof(), FP420G4Hit::setTrackID(), FP420G4Hit::setUnitID(), FP420G4Hit::setVx(), FP420G4Hit::setVy(), FP420G4Hit::setVz(), FP420G4Hit::setX(), FP420G4Hit::setY(), FP420G4Hit::setZ(), StoreHit(), ThetaAtEntry, theTrack, Tof, tsID, tSlice, tSliceID, unitID, Vx, Vy, Vz, X, Y, and Z.

Referenced by ProcessHits().

                            {
 #ifdef debug
   //       << " MVid = " << currentPV->GetMother()->GetCopyNo()
   LogDebug("FP420Sim") << "FP420SD CreateNewHit for"
                        << " PV " << currentPV->GetName() << " PVid = " << currentPV->GetCopyNo() << " Unit " << unitID
                        << std::endl;
   LogDebug("FP420Sim") << " primary " << primaryID << " time slice " << tSliceID << " For Track  "
                        << theTrack->GetTrackID() << " which is a " << theTrack->GetDefinition()->GetParticleName();
 
   if (theTrack->GetTrackID() == 1) {
     LogDebug("FP420Sim") << " of energy " << theTrack->GetTotalEnergy();
   } else {
     LogDebug("FP420Sim") << " daughter of part. " << theTrack->GetParentID();
   }
 
   LogDebug("FP420Sim") << " and created by ";
   if (theTrack->GetCreatorProcess() != NULL)
     LogDebug("FP420Sim") << theTrack->GetCreatorProcess()->GetProcessName();
   else
     LogDebug("FP420Sim") << "NO process";
   LogDebug("FP420Sim") << std::endl;
 #endif
 
   currentHit = new FP420G4Hit;
   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(entrancePoint);
   currentHit->setEntry(hitPoint);
 
   currentHit->setEntryLocalP(hitPointLocal);
   currentHit->setExitLocalP(hitPointLocalExit);
 
   currentHit->setParentId(ParentId);
   currentHit->setVx(Vx);
   currentHit->setVy(Vy);
   currentHit->setVz(Vz);
 
   currentHit->setX(X);
   currentHit->setY(Y);
   currentHit->setZ(Z);
   //AZ:12.10.2007
   //  UpdateHit();
   // buffer for next steps:
   tsID = tSliceID;
   primID = primaryID;
   previousUnitID = unitID;
 
   StoreHit(currentHit);
 }

◆ DrawAll()

void FP420SD::DrawAll ( )

override

Definition at line 404 of file FP420SD.cc.

404 {}

◆ EndOfEvent()

void FP420SD::EndOfEvent ( G4HCofThisEvent * eventHC )

override

Definition at line 326 of file FP420SD.cc.

References FP420G4Hit::getEnergyLoss(), FP420G4Hit::getEntry(), FP420G4Hit::getEntryLocalP(), FP420G4Hit::getExitLocalP(), FP420G4Hit::getPabs(), FP420G4Hit::getParticleType(), FP420G4Hit::getPhiAtEntry(), FP420G4Hit::getThetaAtEntry(), FP420G4Hit::getTof(), FP420G4Hit::getTrackID(), FP420G4Hit::getUnitID(), dqmiolumiharvest::j, LogDebug, nhits, TrackingSlaveSD::processHits(), slave, Summarize(), and theHC.

                                          {
   // here we loop over transient hits and make them persistent
 
   //  if(theHC->entries() > 100){
   //    LogDebug("FP420Sim") << "FP420SD: warning!!! Number of hits exceed 100 and =" << theHC->entries() << "\n";
   //  }
   //  for (int j=0; j<theHC->entries() && j<100; j++) {
   int nhitsHPS240 = 0;
   int nhitsFP420 = 0;
   int nhits = theHC->entries();
   for (int j = 0; j < nhits; j++) {
     FP420G4Hit* aHit = (*theHC)[j];
     if ((fabs(aHit->getTof()) > 780. && fabs(aHit->getTof()) < 840.))
       ++nhitsHPS240;
     if ((fabs(aHit->getTof()) > 1380. && fabs(aHit->getTof()) < 1450.))
       ++nhitsFP420;
     //    if(fabs(aHit->getTof()) < 1700.) {
     if ((fabs(aHit->getTof()) > 780. && fabs(aHit->getTof()) < 840. && nhitsHPS240 < 200.) ||
         (fabs(aHit->getTof()) > 1380. && fabs(aHit->getTof()) < 1450. && nhitsFP420 < 200.)) {
 #ifdef ddebug
       //    LogDebug("FP420SD") << " FP420Hit " << j << " " << *aHit << std::endl;
       LogDebug("FP420Sim") << "hit number" << j << "unit ID = " << aHit->getUnitID() << "\n";
       LogDebug("FP420Sim") << "entry z " << aHit->getEntry().z() << "\n";
       LogDebug("FP420Sim") << "entr theta " << aHit->getThetaAtEntry() << "\n";
 #endif
 
       //    Local3DPoint locExitPoint(0,0,0);
       //  Local3DPoint locEntryPoint(aHit->getEntry().x(),
       //     aHit->getEntry().y(),
       //     aHit->getEntry().z());
       Local3DPoint locExitPoint(aHit->getExitLocalP().x(), aHit->getExitLocalP().y(), aHit->getExitLocalP().z());
       Local3DPoint locEntryPoint(aHit->getEntryLocalP().x(), aHit->getEntryLocalP().y(), aHit->getEntryLocalP().z());
       // implicit conversion (slicing) to PSimHit!!!
       slave->processHits(PSimHit(locEntryPoint,
                                  locExitPoint,             //entryPoint(), exitPoint()  Local3DPoint
                                  aHit->getPabs(),          // pabs()  float
                                  aHit->getTof(),           // tof() float
                                  aHit->getEnergyLoss(),    // energyLoss() float
                                  aHit->getParticleType(),  // particleType()   int
                                  aHit->getUnitID(),        // detUnitId() unsigned int
                                  aHit->getTrackID(),       // trackId() unsigned int
                                  aHit->getThetaAtEntry(),  //  thetaAtEntry()   float
                                  aHit->getPhiAtEntry()));  //  phiAtEntry()   float
 
       //PSimHit( const Local3DPoint& entry, const Local3DPoint& exit,
       //       float pabs, float tof, float eloss, int particleType,
       //       unsigned int detId, unsigned int trackId,
       //       float theta, float phi, unsigned short processType=0) :
 
       //  LocalVector direction = hit.exitPoint() - hit.entryPoint();
       //hit.energyLoss
 
       /*      
           aHit->getEM(),               -
           aHit->getHadr(),             -
           aHit->getIncidentEnergy(),   -
           aHit->getTimeSlice(),       -
           aHit->getEntry(),     -
           aHit->getParentId(),     
           aHit->getEntryLocalP(),  -
           aHit->getExitLocalP(),   -
           aHit->getX(),    -
           aHit->getY(),   -
           aHit->getZ(),   -
           aHit->getVx(),  -
           aHit->getVy(),  -
           aHit->getVz()));  -
       */
     }  //if Tof<1600. if nhits<100
   }    // for loop on hits
 
   Summarize();
 }

◆ fillHits()

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

overridevirtual

Implements SensitiveTkDetector.

Definition at line 419 of file FP420SD.cc.

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

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

◆ getEnergyDeposit()

double FP420SD::getEnergyDeposit ( G4Step * step )

virtual

Definition at line 87 of file FP420SD.cc.

Referenced by GetStepInfo().

87 { return aStep->GetTotalEnergyDeposit(); }

◆ GetStepInfo()

void FP420SD::GetStepInfo ( G4Step * aStep )

private

Definition at line 126 of file FP420SD.cc.

References currentPV, edeposit, edepositEM, edepositHAD, Eloss, emPDG, epPDG, gammaPDG, getEnergyDeposit(), hitPoint, hitPointExit, hitPointLocal, hitPointLocalExit, createfilelist::int, LogDebug, Pabs, ParentId, ParticleType, PhiAtEntry, postStepPoint, preStepPoint, primaryID, setDetUnitId(), ThetaAtEntry, theTrack, Tof, tSlice, tSliceID, unitID, Vx, Vy, Vz, X, Y, and Z.

Referenced by ProcessHits().

                                        {
   preStepPoint = aStep->GetPreStepPoint();
   postStepPoint = aStep->GetPostStepPoint();
   theTrack = aStep->GetTrack();
   hitPoint = preStepPoint->GetPosition();
   currentPV = preStepPoint->GetPhysicalVolume();
   hitPointExit = postStepPoint->GetPosition();
 
   hitPointLocal = preStepPoint->GetTouchable()->GetHistory()->GetTopTransform().TransformPoint(hitPoint);
   hitPointLocalExit = preStepPoint->GetTouchable()->GetHistory()->GetTopTransform().TransformPoint(hitPointExit);
 
   G4int particleCode = theTrack->GetDefinition()->GetPDGEncoding();
   if (particleCode == emPDG || particleCode == epPDG || particleCode == gammaPDG) {
     edepositEM = getEnergyDeposit(aStep);
     edepositHAD = 0.;
   } else {
     edepositEM = 0.;
     edepositHAD = getEnergyDeposit(aStep);
   }
   edeposit = aStep->GetTotalEnergyDeposit();
   tSlice = (postStepPoint->GetGlobalTime()) / nanosecond;
   tSliceID = (int)tSlice;
   unitID = setDetUnitId(aStep);
 #ifdef debug
   LogDebug("FP420Sim") << "unitID=" << unitID << std::endl;
 #endif
   primaryID = theTrack->GetTrackID();
   //  Position     = hitPoint;
   Pabs = aStep->GetPreStepPoint()->GetMomentum().mag() / GeV;
   //Tof          = 1400. + aStep->GetPostStepPoint()->GetGlobalTime()/nanosecond;
   Tof = aStep->GetPostStepPoint()->GetGlobalTime() / nanosecond;
   Eloss = aStep->GetTotalEnergyDeposit() / GeV;
   ParticleType = theTrack->GetDefinition()->GetPDGEncoding();
   ThetaAtEntry = aStep->GetPreStepPoint()->GetPosition().theta() / deg;
   PhiAtEntry = aStep->GetPreStepPoint()->GetPosition().phi() / deg;
 
   ParentId = theTrack->GetParentID();
   Vx = theTrack->GetVertexPosition().x();
   Vy = theTrack->GetVertexPosition().y();
   Vz = theTrack->GetVertexPosition().z();
   X = hitPoint.x();
   Y = hitPoint.y();
   Z = hitPoint.z();
 }

◆ HitExists()

G4bool FP420SD::HitExists ( )

private

Definition at line 175 of file FP420SD.cc.

References currentHit, newFWLiteAna::found, FP420G4Hit::getTimeSliceID(), FP420G4Hit::getTrackID(), FP420G4Hit::getUnitID(), dqmiolumiharvest::j, nhits, previousUnitID, primaryID, primID, ResetForNewPrimary(), theHC, tsID, tSliceID, unitID, and UpdateHit().

Referenced by ProcessHits().

                           {
   if (primaryID < 1) {
     edm::LogWarning("FP420Sim") << "***** FP420SD 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) {
     //AZ:
     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:
 
   G4bool found = false;
 
   //    LogDebug("FP420Sim") << "FP420SD: HCollection=  " << theHC->entries()    <<std::endl;
   int nhits = theHC->entries();
   for (int j = 0; j < nhits && !found; j++) {
     FP420G4Hit* aPreviousHit = (*theHC)[j];
     if (aPreviousHit->getTrackID() == primaryID && aPreviousHit->getTimeSliceID() == tSliceID &&
         aPreviousHit->getUnitID() == unitID) {
       //AZ:
       currentHit = aPreviousHit;
       found = true;
     }
   }
 
   if (found) {
     //AZ:
     UpdateHit();
     return true;
   } else {
     return false;
   }
 }

◆ Initialize()

void FP420SD::Initialize ( G4HCofThisEvent * HCE )

override

Definition at line 89 of file FP420SD.cc.

References bysipixelclustmulteventfilter_cfi::collectionName, hcID, LogDebug, Skims_PA_cff::name, primID, theHC, and tsID.

                                              {
 #ifdef debug
   LogDebug("FP420Sim") << "FP420SD : Initialize called for " << name << std::endl;
 #endif
 
   theHC = new FP420G4HitCollection(GetName(), collectionName[0]);
   if (hcID < 0)
     hcID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]);
   HCE->AddHitsCollection(hcID, theHC);
 
   tsID = -2;
   //  primID = -2;
   primID = 0;
 
 }

◆ PrintAll()

void FP420SD::PrintAll ( )

override

Definition at line 406 of file FP420SD.cc.

References LogDebug, and theHC.

                        {
   LogDebug("FP420Sim") << "FP420SD: Collection " << theHC->GetName() << "\n";
   theHC->PrintAllHits();
 }

◆ ProcessHits()

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

overridevirtual

Implements SensitiveDetector.

Definition at line 106 of file FP420SD.cc.

References CreateNewHit(), edeposit, GetStepInfo(), HitExists(), LogDebug, and theHC.

                                                             {
   if (aStep == nullptr) {
     return true;
   } else {
     GetStepInfo(aStep);
     //   LogDebug("FP420Sim") << edeposit <<std::endl;
 
     //AZ
 #ifdef debug
     LogDebug("FP420Sim") << "FP420SD :  number of hits = " << theHC->entries() << std::endl;
 #endif
 
     if (HitExists() == false && edeposit > 0. && theHC->entries() < 200) {
       CreateNewHit();
       return true;
     }
   }
   return true;
 }

◆ ResetForNewPrimary()

void FP420SD::ResetForNewPrimary ( )

private

Definition at line 217 of file FP420SD.cc.

References entrancePoint, exitPoint, hitPoint, hitPointExit, incidentEnergy, preStepPoint, SetToLocal(), and SetToLocalExit().

Referenced by HitExists().

                                  {
   entrancePoint = SetToLocal(hitPoint);
   exitPoint = SetToLocalExit(hitPointExit);
   incidentEnergy = preStepPoint->GetKineticEnergy();
 }

◆ setDetUnitId()

uint32_t FP420SD::setDetUnitId ( const G4Step * aStep )

overridevirtual

Implements SensitiveDetector.

Definition at line 171 of file FP420SD.cc.

References FP420NumberingScheme::getUnitID(), and numberingScheme.

Referenced by GetStepInfo().

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

◆ SetToLocal()

G4ThreeVector FP420SD::SetToLocal ( const G4ThreeVector & global )

private

Definition at line 314 of file FP420SD.cc.

References preStepPoint, and theEntryPoint.

Referenced by ResetForNewPrimary().

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

◆ SetToLocalExit()

G4ThreeVector FP420SD::SetToLocalExit ( const G4ThreeVector & globalPoint )

private

Definition at line 320 of file FP420SD.cc.

References postStepPoint, and theExitPoint.

Referenced by ResetForNewPrimary().

                                                                       {
   const G4VTouchable* touch = postStepPoint->GetTouchable();
   theExitPoint = touch->GetHistory()->GetTopTransform().TransformPoint(globalPoint);
   return theExitPoint;
 }

◆ StoreHit()

void FP420SD::StoreHit ( FP420G4Hit * hit )

private

Definition at line 223 of file FP420SD.cc.

References theHC.

Referenced by CreateNewHit().

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

◆ Summarize()

void FP420SD::Summarize ( )

private

Definition at line 400 of file FP420SD.cc.

Referenced by EndOfEvent().

400 {}

◆ update() [1/3]

void FP420SD::update ( const BeginOfRun * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun *>.

Definition at line 431 of file FP420SD.cc.

References emPDG, epPDG, gammaPDG, and HiggsValidation_cfi::particleName.

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

                                       {
   G4ParticleTable* theParticleTable = G4ParticleTable::GetParticleTable();
   G4String particleName;
   emPDG = theParticleTable->FindParticle(particleName = "e-")->GetPDGEncoding();
   epPDG = theParticleTable->FindParticle(particleName = "e+")->GetPDGEncoding();
   gammaPDG = theParticleTable->FindParticle(particleName = "gamma")->GetPDGEncoding();
 }

◆ update() [2/3]

void FP420SD::update ( const BeginOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent *>.

Definition at line 425 of file FP420SD.cc.

References clearHits(), eventno, and LogDebug.

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

                                           {
   LogDebug("ForwardSim") << " Dispatched BeginOfEvent for " << GetName() << " !";
   clearHits();
   eventno = (*i)()->GetEventID();
 }

◆ update() [3/3]

void FP420SD::update ( const ::EndOfEvent * )

overrideprivate

Definition at line 439 of file FP420SD.cc.

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

439 {}

◆ UpdateHit()

void FP420SD::UpdateHit ( )

private

Definition at line 293 of file FP420SD.cc.

References FP420G4Hit::addEnergyDeposit(), FP420G4Hit::addEnergyLoss(), currentHit, edepositEM, edepositHAD, Eloss, LogDebug, postStepPoint, previousUnitID, primaryID, primID, tsID, tSliceID, and unitID.

Referenced by HitExists().

                         {
   //
   if (Eloss > 0.) {
     currentHit->addEnergyDeposit(edepositEM, edepositHAD);
 
 #ifdef debug
     LogDebug("FP420Sim") << "updateHit: add eloss " << Eloss << std::endl;
     LogDebug("FP420Sim") << "CurrentHit=" << currentHit << ", PostStepPoint=" << postStepPoint->GetPosition()
                          << std::endl;
 #endif
     //AZ
     //   currentHit->setEnergyLoss(Eloss);
     currentHit->addEnergyLoss(Eloss);
   }
 
   // buffer for next steps:
   tsID = tSliceID;
   primID = primaryID;
   previousUnitID = unitID;
 }