FP420SD Class Reference

#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	(	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::getParameter(), edm::ParameterSet::getUntrackedParameter(), numberingScheme, 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 ( )

override

Definition at line 82 of file FP420SD.cc.

References numberingScheme, and slave.

                  {
  delete slave;
  delete numberingScheme;
}

Member Function Documentation

void FP420SD::clear ( void  )

override

Definition at line 402 of file FP420SD.cc.

{}

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();
}

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);
}

void FP420SD::DrawAll ( )

override

Definition at line 404 of file FP420SD.cc.

{}

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();
}

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();
  }
}

double FP420SD::getEnergyDeposit ( G4Step *  step )

virtual

Definition at line 87 of file FP420SD.cc.

Referenced by GetStepInfo().

{ return aStep->GetTotalEnergyDeposit(); }

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, 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();
}

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;
  }
}

void FP420SD::Initialize ( G4HCofThisEvent *  HCE )

override

Definition at line 89 of file FP420SD.cc.

References ecaldqm::collectionName, hcID, LogDebug, mergeVDriftHistosByStation::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;

}

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();
}

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;
}

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();
}

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));
}

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;
}

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;
}

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);
}

void FP420SD::Summarize ( )

private

Definition at line 400 of file FP420SD.cc.

Referenced by EndOfEvent().

{}

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, and gammaPDG.

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();
}

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();
}

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

{}

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;
}

Member Data Documentation

FP420G4Hit* FP420SD::currentHit

private

Definition at line 97 of file FP420SD.h.

Referenced by CreateNewHit(), HitExists(), and UpdateHit().

G4VPhysicalVolume* FP420SD::currentPV

private

Definition at line 99 of file FP420SD.h.

Referenced by CreateNewHit(), and GetStepInfo().

float FP420SD::edeposit

private

Definition at line 109 of file FP420SD.h.

Referenced by GetStepInfo(), and ProcessHits().

float FP420SD::edepositEM

protected

Definition at line 133 of file FP420SD.h.

Referenced by GetStepInfo(), and UpdateHit().

float FP420SD::edepositHAD

protected

Definition at line 133 of file FP420SD.h.

Referenced by GetStepInfo(), and UpdateHit().

float FP420SD::Eloss

private

Definition at line 117 of file FP420SD.h.

Referenced by CreateNewHit(), GetStepInfo(), and UpdateHit().

G4int FP420SD::emPDG

protected

Definition at line 134 of file FP420SD.h.

Referenced by GetStepInfo(), and update().

G4ThreeVector FP420SD::entrancePoint

private

Definition at line 86 of file FP420SD.h.

Referenced by ResetForNewPrimary().

G4int FP420SD::epPDG

protected

Definition at line 135 of file FP420SD.h.

Referenced by GetStepInfo(), and update().

int FP420SD::eventno

private

Definition at line 130 of file FP420SD.h.

Referenced by update().

G4ThreeVector FP420SD::exitPoint

private

Definition at line 86 of file FP420SD.h.

Referenced by ResetForNewPrimary().

G4int FP420SD::gammaPDG

protected

Definition at line 136 of file FP420SD.h.

Referenced by GetStepInfo(), and update().

G4int FP420SD::hcID

private

Definition at line 92 of file FP420SD.h.

Referenced by Initialize().

G4ThreeVector FP420SD::hitPoint

private

Definition at line 111 of file FP420SD.h.

Referenced by CreateNewHit(), GetStepInfo(), and ResetForNewPrimary().

G4ThreeVector FP420SD::hitPointExit

private

Definition at line 112 of file FP420SD.h.

Referenced by GetStepInfo(), and ResetForNewPrimary().

G4ThreeVector FP420SD::hitPointLocal

private

Definition at line 113 of file FP420SD.h.

Referenced by CreateNewHit(), and GetStepInfo().

G4ThreeVector FP420SD::hitPointLocalExit

private

Definition at line 114 of file FP420SD.h.

Referenced by CreateNewHit(), and GetStepInfo().

float FP420SD::incidentEnergy

private

Definition at line 90 of file FP420SD.h.

Referenced by CreateNewHit(), and ResetForNewPrimary().

FP420NumberingScheme* FP420SD::numberingScheme

private

Definition at line 84 of file FP420SD.h.

Referenced by FP420SD(), setDetUnitId(), and ~FP420SD().

float FP420SD::Pabs

private

Definition at line 115 of file FP420SD.h.

Referenced by CreateNewHit(), and GetStepInfo().

int FP420SD::ParentId

private

Definition at line 123 of file FP420SD.h.

Referenced by CreateNewHit(), and GetStepInfo().

short FP420SD::ParticleType

private

Definition at line 118 of file FP420SD.h.

Referenced by CreateNewHit(), and GetStepInfo().

float FP420SD::PhiAtEntry

private

Definition at line 121 of file FP420SD.h.

Referenced by CreateNewHit(), and GetStepInfo().

G4StepPoint* FP420SD::postStepPoint

private

Definition at line 108 of file FP420SD.h.

Referenced by GetStepInfo(), SetToLocalExit(), and UpdateHit().

G4StepPoint* FP420SD::preStepPoint

private

Definition at line 107 of file FP420SD.h.

Referenced by GetStepInfo(), ResetForNewPrimary(), and SetToLocal().

uint32_t FP420SD::previousUnitID

private

Definition at line 101 of file FP420SD.h.

Referenced by CreateNewHit(), HitExists(), and UpdateHit().

unsigned int FP420SD::primaryID

private

Definition at line 103 of file FP420SD.h.

Referenced by CreateNewHit(), GetStepInfo(), HitExists(), and UpdateHit().

unsigned int FP420SD::primID

private

Definition at line 103 of file FP420SD.h.

Referenced by CreateNewHit(), HitExists(), Initialize(), and UpdateHit().

TrackingSlaveSD* FP420SD::slave

private

Definition at line 83 of file FP420SD.h.

Referenced by clearHits(), EndOfEvent(), fillHits(), FP420SD(), and ~FP420SD().

G4ThreeVector FP420SD::theEntryPoint

private

Definition at line 87 of file FP420SD.h.

Referenced by SetToLocal().

G4ThreeVector FP420SD::theExitPoint

private

Definition at line 88 of file FP420SD.h.

Referenced by SetToLocalExit().

FP420G4HitCollection* FP420SD::theHC

private

Definition at line 93 of file FP420SD.h.

Referenced by EndOfEvent(), HitExists(), Initialize(), PrintAll(), ProcessHits(), and StoreHit().

const SimTrackManager* FP420SD::theManager

private

Definition at line 94 of file FP420SD.h.

float FP420SD::ThetaAtEntry

private

Definition at line 120 of file FP420SD.h.

Referenced by CreateNewHit(), and GetStepInfo().

G4Track* FP420SD::theTrack

private

Definition at line 98 of file FP420SD.h.

Referenced by CreateNewHit(), and GetStepInfo().

float FP420SD::Tof

private

Definition at line 116 of file FP420SD.h.

Referenced by CreateNewHit(), and GetStepInfo().

G4int FP420SD::tsID

private

Definition at line 96 of file FP420SD.h.

Referenced by CreateNewHit(), HitExists(), Initialize(), and UpdateHit().

G4double FP420SD::tSlice

private

Definition at line 105 of file FP420SD.h.

Referenced by CreateNewHit(), and GetStepInfo().

G4int FP420SD::tSliceID

private

Definition at line 102 of file FP420SD.h.

Referenced by CreateNewHit(), GetStepInfo(), HitExists(), and UpdateHit().

uint32_t FP420SD::unitID

private

Definition at line 101 of file FP420SD.h.

Referenced by CreateNewHit(), GetStepInfo(), HitExists(), and UpdateHit().

float FP420SD::Vx

private

Definition at line 124 of file FP420SD.h.

Referenced by CreateNewHit(), and GetStepInfo().

float FP420SD::Vy

private

Definition at line 124 of file FP420SD.h.

Referenced by CreateNewHit(), and GetStepInfo().

float FP420SD::Vz

private

Definition at line 124 of file FP420SD.h.

Referenced by CreateNewHit(), and GetStepInfo().

float FP420SD::X

private

Definition at line 125 of file FP420SD.h.

Referenced by svgfig.Curve.Sample::__repr__(), CreateNewHit(), and GetStepInfo().

float FP420SD::Y

private

Definition at line 125 of file FP420SD.h.

Referenced by CreateNewHit(), and GetStepInfo().

float FP420SD::Z

private

Definition at line 125 of file FP420SD.h.

Referenced by CreateNewHit(), and GetStepInfo().