TotemSD Class Reference

#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)
	~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	(	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::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) {
  //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 ( )

override

Definition at line 84 of file TotemSD.cc.

References numberingScheme, and slave.

                  {
  delete slave;
  delete numberingScheme;
}

Member Function Documentation

void TotemSD::clearHits ( )

overridevirtual

Implements SensitiveDetector.

Definition at line 164 of file TotemSD.cc.

References TrackingSlaveSD::Initialize(), and slave.

Referenced by update().

{ slave->Initialize(); }

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

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

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

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

overridevirtual

Implements SensitiveTkDetector.

Definition at line 151 of file TotemSD.cc.

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

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

void TotemSD::getStepInfo ( const G4Step *  aStep )

private

Definition at line 173 of file TotemSD.cc.

References currentPV, edeposit, Eloss, GeV, hitPoint, mergeVDriftHistosByStation::name, Pabs, ParentId, ParticleType, HLT_FULL_cff::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()) / 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();
}

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

void TotemSD::Initialize ( G4HCofThisEvent *  HCE )

override

Definition at line 105 of file TotemSD.cc.

References ecaldqm::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;
}

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, and mathSSE::sqrt().

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

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

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

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

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

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

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

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

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

Member Data Documentation

TotemG4Hit* TotemSD::currentHit

private

Definition at line 93 of file TotemSD.h.

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

G4VPhysicalVolume* TotemSD::currentPV

private

Definition at line 95 of file TotemSD.h.

Referenced by createNewHit(), and getStepInfo().

float TotemSD::edeposit

private

Definition at line 102 of file TotemSD.h.

Referenced by getStepInfo(), and ProcessHits().

float TotemSD::Eloss

private

Definition at line 108 of file TotemSD.h.

Referenced by createNewHit(), createNewHitEvo(), getStepInfo(), and updateHit().

G4ThreeVector TotemSD::entrancePoint

private

Definition at line 84 of file TotemSD.h.

Referenced by resetForNewPrimary().

G4int TotemSD::hcID

private

Definition at line 88 of file TotemSD.h.

Referenced by Initialize().

G4ThreeVector TotemSD::hitPoint

private

Definition at line 103 of file TotemSD.h.

Referenced by getStepInfo(), and resetForNewPrimary().

float TotemSD::incidentEnergy

private

Definition at line 85 of file TotemSD.h.

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

TotemVDetectorOrganization* TotemSD::numberingScheme

private

Definition at line 77 of file TotemSD.h.

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

float TotemSD::Pabs

private

Definition at line 106 of file TotemSD.h.

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

int TotemSD::ParentId

private

Definition at line 114 of file TotemSD.h.

Referenced by createNewHit(), createNewHitEvo(), getStepInfo(), and ProcessHits().

short TotemSD::ParticleType

private

Definition at line 109 of file TotemSD.h.

Referenced by createNewHit(), createNewHitEvo(), getStepInfo(), and ProcessHits().

float TotemSD::PhiAtEntry

private

Definition at line 112 of file TotemSD.h.

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

G4ThreeVector TotemSD::Posizio

private

Definition at line 105 of file TotemSD.h.

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

const G4StepPoint* TotemSD::postStepPoint

private

Definition at line 101 of file TotemSD.h.

Referenced by getStepInfo(), and updateHit().

const G4StepPoint* TotemSD::preStepPoint

private

Definition at line 100 of file TotemSD.h.

Referenced by getStepInfo(), resetForNewPrimary(), and setToLocal().

uint32_t TotemSD::previousUnitID

private

Definition at line 96 of file TotemSD.h.

Referenced by hitExists(), and updateHit().

int TotemSD::primaryID

private

Definition at line 97 of file TotemSD.h.

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

G4int TotemSD::primID

private

Definition at line 86 of file TotemSD.h.

Referenced by hitExists(), Initialize(), storeHit(), and updateHit().

TrackingSlaveSD* TotemSD::slave

private

Definition at line 76 of file TotemSD.h.

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

TotemG4HitCollection* TotemSD::theHC

private

Definition at line 89 of file TotemSD.h.

Referenced by EndOfEvent(), hitExists(), Initialize(), PrintAll(), and storeHit().

const SimTrackManager* TotemSD::theManager

private

Definition at line 90 of file TotemSD.h.

float TotemSD::ThetaAtEntry

private

Definition at line 111 of file TotemSD.h.

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

G4Track* TotemSD::theTrack

private

Definition at line 94 of file TotemSD.h.

Referenced by createNewHit(), and getStepInfo().

float TotemSD::Tof

private

Definition at line 107 of file TotemSD.h.

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

int TotemSD::tsID

private

Definition at line 92 of file TotemSD.h.

Referenced by hitExists(), Initialize(), and updateHit().

double TotemSD::tSlice

private

Definition at line 98 of file TotemSD.h.

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

int TotemSD::tSliceID

private

Definition at line 97 of file TotemSD.h.

Referenced by createNewHit(), getStepInfo(), hitExists(), and updateHit().

uint32_t TotemSD::unitID

private

Definition at line 96 of file TotemSD.h.

Referenced by createNewHit(), createNewHitEvo(), getStepInfo(), hitExists(), ProcessHits(), and updateHit().

float TotemSD::Vx

private

Definition at line 115 of file TotemSD.h.

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

float TotemSD::Vy

private

Definition at line 115 of file TotemSD.h.

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

float TotemSD::Vz

private

Definition at line 115 of file TotemSD.h.

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