CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
List of all members | Public Member Functions | Static Public Attributes | Private Member Functions | Private Attributes | Static Private Attributes
TotemRPSD Class Reference

#include <TotemRPSD.h>

Inheritance diagram for TotemRPSD:
SensitiveTkDetector Observer< const BeginOfEvent * > Observer< const EndOfEvent * > SensitiveDetector

Public Member Functions

void clear () override
 
void DrawAll () override
 
void EndOfEvent (G4HCofThisEvent *eventHC) override
 
void fillHits (edm::PSimHitContainer &, const std::string &) override
 
void Initialize (G4HCofThisEvent *HCE) override
 
void PrintAll () override
 
void printHitInfo ()
 
 TotemRPSD (const std::string &, const SensitiveDetectorCatalog &, edm::ParameterSet const &, const SimTrackManager *)
 
 ~TotemRPSD () 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 ()
 
- Public Member Functions inherited from Observer< const EndOfEvent * >
 Observer ()
 
void slotForUpdate (const EndOfEvent *iT)
 
virtual ~Observer ()
 

Static Public Attributes

static constexpr double rp_garage_position_ = 40.0
 

Private Member Functions

void clearHits () override
 
void createNewHit ()
 
G4bool hitExists ()
 
bool isPrimary (const G4Track *track)
 
bool ProcessHits (G4Step *step, G4TouchableHistory *tHistory) override
 
void resetForNewPrimary ()
 
uint32_t setDetUnitId (const G4Step *step) override
 
void setNumberingScheme (TotemRPVDetectorOrganization *scheme)
 
G4ThreeVector setToLocal (const G4ThreeVector &globalPoint)
 
void stepInfo (const G4Step *aStep)
 
void storeHit (TotemRPG4Hit *)
 
void summarize ()
 
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

TotemRPG4HitcurrentHit_
 
G4VPhysicalVolume * currentPV_
 
double Eloss_
 
G4ThreeVector entrancePoint_
 
int eventno_
 
G4ThreeVector exitPoint_
 
G4int hcID_
 
G4ThreeVector hitPoint_
 
double incidentEnergy_
 
G4String name_
 
std::unique_ptr
< TotemRPVDetectorOrganization
numberingScheme_
 
double Pabs_
 
int ParentId_
 
short ParticleType_
 
double PhiAtEntry_
 
G4StepPoint * postStepPoint_
 
G4StepPoint * preStepPoint_
 
G4int primaryID_
 
std::unique_ptr< TrackingSlaveSDslave_
 
TotemRPG4HitCollectiontheHC_
 
G4ThreeVector theLocalEntryPoint_
 
G4ThreeVector theLocalExitPoint_
 
double thePx_
 
double thePy_
 
double thePz_
 
double ThetaAtEntry_
 
G4Track * theTrack_
 
double Tof_
 
G4double tSlice_
 
G4int tSliceID_
 
unsigned int unitID_
 
int verbosity_
 
double Vx_
 
double Vy_
 
double Vz_
 

Static Private Attributes

static constexpr unsigned int maxTotemHits_ = 15000
 

Additional Inherited Members

- Protected Types inherited from SensitiveDetector
enum  coordinates { WorldCoordinates, LocalCoordinates }
 
- Protected Member Functions inherited from SensitiveDetector
TrackInformationcmsTrackInformation (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 30 of file TotemRPSD.h.

Constructor & Destructor Documentation

TotemRPSD::TotemRPSD ( const std::string &  name_,
const SensitiveDetectorCatalog clg,
edm::ParameterSet const &  p,
const SimTrackManager manager 
)

Definition at line 33 of file TotemRPSD.cc.

References ecaldqm::collectionName, edm::ParameterSet::getParameter(), name_, numberingScheme_, slave_, and verbosity_.

37  : SensitiveTkDetector(name_, clg),
38  numberingScheme_(nullptr),
39  hcID_(-1),
40  theHC_(nullptr),
41  currentHit_(nullptr),
42  theTrack_(nullptr),
43  currentPV_(nullptr),
44  unitID_(0),
45  preStepPoint_(nullptr),
46  postStepPoint_(nullptr),
47  eventno_(0) {
48  collectionName.insert(name_);
49 
50  edm::ParameterSet m_Anal = p.getParameter<edm::ParameterSet>("TotemRPSD");
51  verbosity_ = m_Anal.getParameter<int>("Verbosity");
52 
53  slave_ = std::make_unique<TrackingSlaveSD>(name_);
54 
55  if (name_ == "TotemHitsRP") {
56  numberingScheme_ = std::make_unique<PPSStripNumberingScheme>(3);
57  } else {
58  edm::LogWarning("TotemRP") << "TotemRPSD: ReadoutName not supported\n";
59  }
60 
61  edm::LogInfo("TotemRP") << "TotemRPSD: Instantiation completed";
62 }
G4Track * theTrack_
Definition: TotemRPSD.h:86
G4VPhysicalVolume * currentPV_
Definition: TotemRPSD.h:87
std::unique_ptr< TotemRPVDetectorOrganization > numberingScheme_
Definition: TotemRPSD.h:57
int eventno_
Definition: TotemRPSD.h:113
TotemRPG4HitCollection * theHC_
Definition: TotemRPSD.h:83
std::string const collectionName[nCollections]
Definition: Collections.h:47
G4StepPoint * postStepPoint_
Definition: TotemRPSD.h:93
int verbosity_
Definition: TotemRPSD.h:60
unsigned int unitID_
Definition: TotemRPSD.h:88
G4StepPoint * preStepPoint_
Definition: TotemRPSD.h:92
Log< level::Info, false > LogInfo
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
G4int hcID_
Definition: TotemRPSD.h:82
Log< level::Warning, false > LogWarning
SensitiveTkDetector(const std::string &iname, const SensitiveDetectorCatalog &clg)
std::unique_ptr< TrackingSlaveSD > slave_
Definition: TotemRPSD.h:56
TotemRPG4Hit * currentHit_
Definition: TotemRPSD.h:85
G4String name_
Definition: TotemRPSD.h:81
TotemRPSD::~TotemRPSD ( )
override

Definition at line 64 of file TotemRPSD.cc.

64 {}

Member Function Documentation

void TotemRPSD::clear ( void  )
override

Definition at line 259 of file TotemRPSD.cc.

259 {}
void TotemRPSD::clearHits ( )
overrideprivatevirtual

Implements SensitiveDetector.

Definition at line 286 of file TotemRPSD.cc.

References slave_.

Referenced by update().

286 { slave_->Initialize(); }
std::unique_ptr< TrackingSlaveSD > slave_
Definition: TotemRPSD.h:56
void TotemRPSD::createNewHit ( )
private

Definition at line 189 of file TotemRPSD.cc.

References currentHit_, Eloss_, exitPoint_, hitPoint_, incidentEnergy_, Pabs_, ParentId_, ParticleType_, PhiAtEntry_, primaryID_, rp_garage_position_, TotemRPG4Hit::setEnergyLoss(), TotemRPG4Hit::setEntry(), TotemRPG4Hit::setExit(), TotemRPG4Hit::setIncidentEnergy(), TotemRPG4Hit::setLocalEntry(), TotemRPG4Hit::setLocalExit(), TotemRPG4Hit::setP(), TotemRPG4Hit::setParentId(), TotemRPG4Hit::setParticleType(), TotemRPG4Hit::setPhiAtEntry(), TotemRPG4Hit::setPx(), TotemRPG4Hit::setPy(), TotemRPG4Hit::setPz(), TotemRPG4Hit::setThetaAtEntry(), TotemRPG4Hit::setTimeSlice(), TotemRPG4Hit::setTof(), TotemRPG4Hit::setTrackID(), TotemRPG4Hit::setUnitID(), TotemRPG4Hit::setVx(), TotemRPG4Hit::setVy(), TotemRPG4Hit::setVz(), storeHit(), theLocalEntryPoint_, theLocalExitPoint_, thePx_, thePy_, thePz_, ThetaAtEntry_, Tof_, tSlice_, unitID_, Vx_, Vy_, and Vz_.

Referenced by ProcessHits().

189  {
190  // Protect against creating hits in detectors not inserted
191  double outrangeX = hitPoint_.x();
192  double outrangeY = hitPoint_.y();
193  if (fabs(outrangeX) > rp_garage_position_)
194  return;
195  if (fabs(outrangeY) > rp_garage_position_)
196  return;
197  // end protection
198 
204 
211 
216 
221 
225 
227 }
void setEntry(G4ThreeVector xyz)
Definition: TotemRPG4Hit.cc:76
void setTof(double e)
void setP(double e)
double Tof_
Definition: TotemRPSD.h:101
void setIncidentEnergy(double e)
Definition: TotemRPG4Hit.cc:87
void storeHit(TotemRPG4Hit *)
Definition: TotemRPSD.cc:180
void setVx(double p)
void setPy(double p)
double Pabs_
Definition: TotemRPSD.h:99
G4double tSlice_
Definition: TotemRPSD.h:90
static constexpr double rp_garage_position_
Definition: TotemRPSD.h:44
void setLocalEntry(const G4ThreeVector &theLocalEntryPoint)
Definition: TotemRPG4Hit.cc:82
void setExit(G4ThreeVector xyz)
Definition: TotemRPG4Hit.cc:79
G4ThreeVector theLocalEntryPoint_
Definition: TotemRPSD.h:96
int ParentId_
Definition: TotemRPSD.h:108
short ParticleType_
Definition: TotemRPSD.h:103
unsigned int unitID_
Definition: TotemRPSD.h:88
G4int primaryID_
Definition: TotemRPSD.h:89
void setLocalExit(const G4ThreeVector &theLocalExitPoint)
Definition: TotemRPG4Hit.cc:84
void setThetaAtEntry(double t)
void setUnitID(unsigned int i)
Definition: TotemRPG4Hit.cc:93
double PhiAtEntry_
Definition: TotemRPSD.h:106
double Vy_
Definition: TotemRPSD.h:109
G4ThreeVector theLocalExitPoint_
Definition: TotemRPSD.h:97
double ThetaAtEntry_
Definition: TotemRPSD.h:105
double Vz_
Definition: TotemRPSD.h:109
void setPhiAtEntry(double f)
double thePz_
Definition: TotemRPSD.h:100
void setVz(double p)
double thePx_
Definition: TotemRPSD.h:100
double incidentEnergy_
Definition: TotemRPSD.h:79
void setParentId(int p)
void setTrackID(int i)
Definition: TotemRPG4Hit.cc:90
G4ThreeVector exitPoint_
Definition: TotemRPSD.h:95
void setPz(double p)
void setParticleType(short i)
double Eloss_
Definition: TotemRPSD.h:102
void setEnergyLoss(double e)
void setVy(double p)
double thePy_
Definition: TotemRPSD.h:100
void setPx(double p)
G4ThreeVector hitPoint_
Definition: TotemRPSD.h:94
TotemRPG4Hit * currentHit_
Definition: TotemRPSD.h:85
double Vx_
Definition: TotemRPSD.h:109
void setTimeSlice(double d)
Definition: TotemRPG4Hit.cc:96
void TotemRPSD::DrawAll ( )
override

Definition at line 261 of file TotemRPSD.cc.

261 {}
void TotemRPSD::EndOfEvent ( G4HCofThisEvent *  eventHC)
override

Definition at line 236 of file TotemRPSD.cc.

References TotemRPG4Hit::energyLoss(), mps_splice::entry, beamvalidation::exit(), dqmiolumiharvest::j, TotemRPG4Hit::localEntry(), TotemRPG4Hit::localExit(), maxTotemHits_, TotemRPG4Hit::p(), TotemRPG4Hit::particleType(), TotemRPG4Hit::phiAtEntry(), slave_, summarize(), theHC_, TotemRPG4Hit::thetaAtEntry(), TotemRPG4Hit::tof(), TotemRPG4Hit::trackID(), and TotemRPG4Hit::unitID().

236  {
237  // here we loop over transient hits and make them persistent
238  for (unsigned int j = 0; j < (unsigned int)theHC_->entries() && j < maxTotemHits_; j++) {
239  TotemRPG4Hit* aHit = (*theHC_)[j];
240 
241  Local3DPoint entry(aHit->localEntry().x(), aHit->localEntry().y(), aHit->localEntry().z());
242  Local3DPoint exit(aHit->localExit().x(), aHit->localExit().y(), aHit->localExit().z());
243  slave_->processHits(PSimHit(entry,
244  exit,
245  aHit->p(),
246  aHit->tof(),
247  aHit->energyLoss(),
248  aHit->particleType(),
249  aHit->unitID(),
250  aHit->trackID(),
251  aHit->thetaAtEntry(),
252  aHit->phiAtEntry()));
253  }
254  summarize();
255 }
G4ThreeVector localEntry() const
Definition: TotemRPG4Hit.cc:81
unsigned int trackID() const
Definition: TotemRPG4Hit.cc:89
int particleType() const
int unitID() const
Definition: TotemRPG4Hit.cc:92
TotemRPG4HitCollection * theHC_
Definition: TotemRPSD.h:83
double p() const
Definition: TotemRPG4Hit.cc:99
double phiAtEntry() const
double energyLoss() const
double thetaAtEntry() const
static constexpr unsigned int maxTotemHits_
Definition: TotemRPSD.h:47
G4ThreeVector localExit() const
Definition: TotemRPG4Hit.cc:83
double tof() const
list entry
Definition: mps_splice.py:68
void summarize()
Definition: TotemRPSD.cc:257
std::unique_ptr< TrackingSlaveSD > slave_
Definition: TotemRPSD.h:56
void TotemRPSD::fillHits ( edm::PSimHitContainer c,
const std::string &  n 
)
overridevirtual

Implements SensitiveTkDetector.

Definition at line 268 of file TotemRPSD.cc.

References dqmiodumpmetadata::n, and slave_.

268  {
269  if (slave_->name() == n) {
270  c = slave_->hits();
271  }
272 }
const edm::EventSetup & c
std::unique_ptr< TrackingSlaveSD > slave_
Definition: TotemRPSD.h:56
G4bool TotemRPSD::hitExists ( )
private
void TotemRPSD::Initialize ( G4HCofThisEvent *  HCE)
override

Definition at line 66 of file TotemRPSD.cc.

References ecaldqm::collectionName, hcID_, LogDebug, name_, and theHC_.

66  {
67  LogDebug("TotemRP") << "TotemRPSD : Initialize called for " << name_;
68 
69  theHC_ = new TotemRPG4HitCollection(GetName(), collectionName[0]);
70  G4SDManager::GetSDMpointer()->AddNewCollection(name_, collectionName[0]);
71 
72  if (hcID_ < 0)
73  hcID_ = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]);
74  HCE->AddHitsCollection(hcID_, theHC_);
75 }
G4THitsCollection< TotemRPG4Hit > TotemRPG4HitCollection
TotemRPG4HitCollection * theHC_
Definition: TotemRPSD.h:83
std::string const collectionName[nCollections]
Definition: Collections.h:47
G4int hcID_
Definition: TotemRPSD.h:82
G4String name_
Definition: TotemRPSD.h:81
#define LogDebug(id)
bool TotemRPSD::isPrimary ( const G4Track *  track)
private

Definition at line 288 of file TotemRPSD.cc.

References info(), and TrackInformation::isPrimary().

Referenced by stepInfo().

288  {
289  TrackInformation* info = dynamic_cast<TrackInformation*>(track->GetUserInformation());
290  return info && info->isPrimary();
291 }
static const TGPicture * info(bool iBackgroundIsBlack)
bool isPrimary() const
void TotemRPSD::PrintAll ( )
override

Definition at line 263 of file TotemRPSD.cc.

References LogDebug, and theHC_.

263  {
264  LogDebug("TotemRP") << "TotemRPSD: Collection " << theHC_->GetName() << std::endl;
265  theHC_->PrintAllHits();
266 }
TotemRPG4HitCollection * theHC_
Definition: TotemRPSD.h:83
#define LogDebug(id)
void TotemRPSD::printHitInfo ( )

Definition at line 77 of file TotemRPSD.cc.

References currentPV_, Eloss_, LogDebug, postStepPoint_, preStepPoint_, primaryID_, gpuVertexFinder::printf(), theTrack_, tSliceID_, and unitID_.

77  {
78  LogDebug("TotemRP") << theTrack_->GetDefinition()->GetParticleName() << " TotemRPSD CreateNewHit for"
79  << " PV " << currentPV_->GetName() << " PVid = " << currentPV_->GetCopyNo() << " Unit "
80  << unitID_;
81  LogDebug("TotemRP") << " primary " << primaryID_ << " time slice " << tSliceID_ << " of energy "
82  << theTrack_->GetTotalEnergy() << " Eloss_ " << Eloss_ << " positions ";
83  printf("(%10f,%10f,%10f)",
84  preStepPoint_->GetPosition().x(),
85  preStepPoint_->GetPosition().y(),
86  preStepPoint_->GetPosition().z());
87  printf("(%10f,%10f,%10f)",
88  postStepPoint_->GetPosition().x(),
89  postStepPoint_->GetPosition().y(),
90  postStepPoint_->GetPosition().z());
91  LogDebug("TotemRP") << " positions "
92  << "(" << postStepPoint_->GetPosition().x() << "," << postStepPoint_->GetPosition().y() << ","
93  << postStepPoint_->GetPosition().z() << ")"
94  << " For Track " << theTrack_->GetTrackID() << " which is a "
95  << theTrack_->GetDefinition()->GetParticleName();
96 
97  if (theTrack_->GetTrackID() == 1) {
98  LogDebug("TotemRP") << " primary particle ";
99  } else {
100  LogDebug("TotemRP") << " daughter of part. " << theTrack_->GetParentID();
101  }
102 
103  LogDebug("TotemRP") << " and created by ";
104 
105  if (theTrack_->GetCreatorProcess() != nullptr)
106  LogDebug("TotemRP") << theTrack_->GetCreatorProcess()->GetProcessName();
107  else
108  LogDebug("TotemRP") << "NO process";
109 
110  LogDebug("TotemRP") << std::endl;
111 }
G4Track * theTrack_
Definition: TotemRPSD.h:86
G4int tSliceID_
Definition: TotemRPSD.h:89
G4VPhysicalVolume * currentPV_
Definition: TotemRPSD.h:87
printf("params %d %f %f %f\n", minT, eps, errmax, chi2max)
G4StepPoint * postStepPoint_
Definition: TotemRPSD.h:93
unsigned int unitID_
Definition: TotemRPSD.h:88
G4int primaryID_
Definition: TotemRPSD.h:89
G4StepPoint * preStepPoint_
Definition: TotemRPSD.h:92
double Eloss_
Definition: TotemRPSD.h:102
#define LogDebug(id)
bool TotemRPSD::ProcessHits ( G4Step *  step,
G4TouchableHistory *  tHistory 
)
overrideprivatevirtual

Implements SensitiveDetector.

Definition at line 113 of file TotemRPSD.cc.

References createNewHit(), and stepInfo().

113  {
114  if (aStep == nullptr) {
115  return true;
116  } else {
117  stepInfo(aStep);
118 
119  createNewHit();
120  return true;
121  }
122 }
void createNewHit()
Definition: TotemRPSD.cc:189
void stepInfo(const G4Step *aStep)
Definition: TotemRPSD.cc:124
void TotemRPSD::resetForNewPrimary ( )
private
uint32_t TotemRPSD::setDetUnitId ( const G4Step *  step)
overrideprivatevirtual

Implements SensitiveDetector.

Definition at line 176 of file TotemRPSD.cc.

References numberingScheme_.

Referenced by stepInfo().

176  {
177  return (numberingScheme_ == nullptr ? 0 : numberingScheme_->unitID(aStep));
178 }
std::unique_ptr< TotemRPVDetectorOrganization > numberingScheme_
Definition: TotemRPSD.h:57
void TotemRPSD::setNumberingScheme ( TotemRPVDetectorOrganization scheme)
private

Definition at line 273 of file TotemRPSD.cc.

References LogDebug, and numberingScheme_.

273  {
274  if (scheme) {
275  LogDebug("TotemRP") << "TotemRPSD: updates numbering scheme for " << GetName();
276  numberingScheme_.reset(scheme);
277  }
278 }
std::unique_ptr< TotemRPVDetectorOrganization > numberingScheme_
Definition: TotemRPSD.h:57
#define LogDebug(id)
G4ThreeVector TotemRPSD::setToLocal ( const G4ThreeVector &  globalPoint)
private

Definition at line 229 of file TotemRPSD.cc.

References preStepPoint_.

Referenced by stepInfo().

229  {
230  G4ThreeVector localPoint;
231  const G4VTouchable* touch = preStepPoint_->GetTouchable();
232  localPoint = touch->GetHistory()->GetTopTransform().TransformPoint(global);
233  return localPoint;
234 }
G4StepPoint * preStepPoint_
Definition: TotemRPSD.h:92
void TotemRPSD::stepInfo ( const G4Step *  aStep)
private

Definition at line 124 of file TotemRPSD.cc.

References SensitiveDetector::ConvertToLocal3DPoint(), currentPV_, Eloss_, exitPoint_, GeV, hitPoint_, isPrimary(), LogDebug, name_, Pabs_, ParentId_, HLT_FULL_cff::particleType, ParticleType_, PV3DBase< T, PVType, FrameType >::phi(), PhiAtEntry_, postStepPoint_, preStepPoint_, primaryID_, setDetUnitId(), setToLocal(), theLocalEntryPoint_, theLocalExitPoint_, thePx_, thePy_, thePz_, PV3DBase< T, PVType, FrameType >::theta(), ThetaAtEntry_, theTrack_, Tof_, tSlice_, tSliceID_, unitID_, verbosity_, Vx_, Vy_, and Vz_.

Referenced by ProcessHits().

124  {
125  preStepPoint_ = aStep->GetPreStepPoint();
126  postStepPoint_ = aStep->GetPostStepPoint();
127  theTrack_ = aStep->GetTrack();
128  hitPoint_ = preStepPoint_->GetPosition();
129  exitPoint_ = postStepPoint_->GetPosition();
130  currentPV_ = preStepPoint_->GetPhysicalVolume();
133 
134  G4String name_ = currentPV_->GetName();
135  name_.assign(name_, 0, 4);
136  G4String particleType = theTrack_->GetDefinition()->GetParticleName();
137  tSlice_ = (postStepPoint_->GetGlobalTime()) / nanosecond;
138  tSliceID_ = (int)tSlice_;
139  unitID_ = setDetUnitId(aStep);
140 
141  if (verbosity_)
142  LogDebug("TotemRP") << "UNIT " << unitID_ << std::endl;
143 
144  primaryID_ = theTrack_->GetTrackID();
145 
146  Pabs_ = (aStep->GetPreStepPoint()->GetMomentum().mag()) / GeV;
147  thePx_ = (aStep->GetPreStepPoint()->GetMomentum().x()) / GeV;
148  thePy_ = (aStep->GetPreStepPoint()->GetMomentum().y()) / GeV;
149  thePz_ = (aStep->GetPreStepPoint()->GetMomentum().z()) / GeV;
150 
151  Tof_ = aStep->GetPostStepPoint()->GetGlobalTime() / nanosecond;
152  Eloss_ = aStep->GetTotalEnergyDeposit() / GeV;
153  ParticleType_ = theTrack_->GetDefinition()->GetPDGEncoding();
154 
155  //corrected phi and theta treatment
156  G4ThreeVector gmd = aStep->GetPreStepPoint()->GetMomentumDirection();
157  // convert it to local frame
158  G4ThreeVector lmd = ((G4TouchableHistory*)(aStep->GetPreStepPoint()->GetTouchable()))
159  ->GetHistory()
160  ->GetTopTransform()
161  .TransformAxis(gmd);
163  ThetaAtEntry_ = lnmd.theta();
164  PhiAtEntry_ = lnmd.phi();
165 
166  if (isPrimary(theTrack_))
167  ParentId_ = 0;
168  else
169  ParentId_ = theTrack_->GetParentID();
170 
171  Vx_ = theTrack_->GetVertexPosition().x() / mm;
172  Vy_ = theTrack_->GetVertexPosition().y() / mm;
173  Vz_ = theTrack_->GetVertexPosition().z() / mm;
174 }
const double GeV
Definition: MathUtil.h:16
G4Track * theTrack_
Definition: TotemRPSD.h:86
Local3DPoint ConvertToLocal3DPoint(const G4ThreeVector &point) const
double Tof_
Definition: TotemRPSD.h:101
G4ThreeVector setToLocal(const G4ThreeVector &globalPoint)
Definition: TotemRPSD.cc:229
G4int tSliceID_
Definition: TotemRPSD.h:89
G4VPhysicalVolume * currentPV_
Definition: TotemRPSD.h:87
Geom::Phi< T > phi() const
Definition: PV3DBase.h:66
double Pabs_
Definition: TotemRPSD.h:99
G4double tSlice_
Definition: TotemRPSD.h:90
Geom::Theta< T > theta() const
Definition: PV3DBase.h:72
G4ThreeVector theLocalEntryPoint_
Definition: TotemRPSD.h:96
int ParentId_
Definition: TotemRPSD.h:108
uint32_t setDetUnitId(const G4Step *step) override
Definition: TotemRPSD.cc:176
short ParticleType_
Definition: TotemRPSD.h:103
G4StepPoint * postStepPoint_
Definition: TotemRPSD.h:93
int verbosity_
Definition: TotemRPSD.h:60
unsigned int unitID_
Definition: TotemRPSD.h:88
G4int primaryID_
Definition: TotemRPSD.h:89
G4StepPoint * preStepPoint_
Definition: TotemRPSD.h:92
double PhiAtEntry_
Definition: TotemRPSD.h:106
double Vy_
Definition: TotemRPSD.h:109
G4ThreeVector theLocalExitPoint_
Definition: TotemRPSD.h:97
double ThetaAtEntry_
Definition: TotemRPSD.h:105
double Vz_
Definition: TotemRPSD.h:109
double thePz_
Definition: TotemRPSD.h:100
double thePx_
Definition: TotemRPSD.h:100
G4ThreeVector exitPoint_
Definition: TotemRPSD.h:95
bool isPrimary(const G4Track *track)
Definition: TotemRPSD.cc:288
double Eloss_
Definition: TotemRPSD.h:102
double thePy_
Definition: TotemRPSD.h:100
G4ThreeVector hitPoint_
Definition: TotemRPSD.h:94
double Vx_
Definition: TotemRPSD.h:109
G4String name_
Definition: TotemRPSD.h:81
#define LogDebug(id)
void TotemRPSD::storeHit ( TotemRPG4Hit hit)
private

Definition at line 180 of file TotemRPSD.cc.

References LogDebug, theHC_, and verbosity_.

Referenced by createNewHit().

180  {
181  if (hit == nullptr) {
182  if (verbosity_)
183  LogDebug("TotemRP") << "TotemRPSD: hit to be stored is NULL !!" << std::endl;
184  return;
185  }
186  theHC_->insert(hit);
187 }
TotemRPG4HitCollection * theHC_
Definition: TotemRPSD.h:83
int verbosity_
Definition: TotemRPSD.h:60
#define LogDebug(id)
void TotemRPSD::summarize ( )
private

Definition at line 257 of file TotemRPSD.cc.

Referenced by EndOfEvent().

257 {}
void TotemRPSD::update ( const BeginOfEvent )
overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 279 of file TotemRPSD.cc.

References clearHits(), and eventno_.

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

279  {
280  clearHits();
281  eventno_ = (*i)()->GetEventID();
282 }
int eventno_
Definition: TotemRPSD.h:113
void clearHits() override
Definition: TotemRPSD.cc:286
void TotemRPSD::update ( const ::EndOfEvent )
overrideprivate
void TotemRPSD::updateHit ( )
private

Member Data Documentation

TotemRPG4Hit* TotemRPSD::currentHit_
private

Definition at line 85 of file TotemRPSD.h.

Referenced by createNewHit().

G4VPhysicalVolume* TotemRPSD::currentPV_
private

Definition at line 87 of file TotemRPSD.h.

Referenced by printHitInfo(), and stepInfo().

double TotemRPSD::Eloss_
private

Definition at line 102 of file TotemRPSD.h.

Referenced by createNewHit(), printHitInfo(), and stepInfo().

G4ThreeVector TotemRPSD::entrancePoint_
private

Definition at line 78 of file TotemRPSD.h.

int TotemRPSD::eventno_
private

Definition at line 113 of file TotemRPSD.h.

Referenced by update().

G4ThreeVector TotemRPSD::exitPoint_
private

Definition at line 95 of file TotemRPSD.h.

Referenced by createNewHit(), and stepInfo().

G4int TotemRPSD::hcID_
private

Definition at line 82 of file TotemRPSD.h.

Referenced by Initialize().

G4ThreeVector TotemRPSD::hitPoint_
private

Definition at line 94 of file TotemRPSD.h.

Referenced by createNewHit(), and stepInfo().

double TotemRPSD::incidentEnergy_
private

Definition at line 79 of file TotemRPSD.h.

Referenced by createNewHit().

constexpr unsigned int TotemRPSD::maxTotemHits_ = 15000
staticprivate

Definition at line 47 of file TotemRPSD.h.

Referenced by EndOfEvent().

G4String TotemRPSD::name_
private

Definition at line 81 of file TotemRPSD.h.

Referenced by Initialize(), stepInfo(), and TotemRPSD().

std::unique_ptr<TotemRPVDetectorOrganization> TotemRPSD::numberingScheme_
private

Definition at line 57 of file TotemRPSD.h.

Referenced by setDetUnitId(), setNumberingScheme(), and TotemRPSD().

double TotemRPSD::Pabs_
private

Definition at line 99 of file TotemRPSD.h.

Referenced by createNewHit(), and stepInfo().

int TotemRPSD::ParentId_
private

Definition at line 108 of file TotemRPSD.h.

Referenced by createNewHit(), and stepInfo().

short TotemRPSD::ParticleType_
private

Definition at line 103 of file TotemRPSD.h.

Referenced by createNewHit(), and stepInfo().

double TotemRPSD::PhiAtEntry_
private

Definition at line 106 of file TotemRPSD.h.

Referenced by createNewHit(), and stepInfo().

G4StepPoint* TotemRPSD::postStepPoint_
private

Definition at line 93 of file TotemRPSD.h.

Referenced by printHitInfo(), and stepInfo().

G4StepPoint* TotemRPSD::preStepPoint_
private

Definition at line 92 of file TotemRPSD.h.

Referenced by printHitInfo(), setToLocal(), and stepInfo().

G4int TotemRPSD::primaryID_
private

Definition at line 89 of file TotemRPSD.h.

Referenced by createNewHit(), printHitInfo(), and stepInfo().

constexpr double TotemRPSD::rp_garage_position_ = 40.0
static

Definition at line 44 of file TotemRPSD.h.

Referenced by createNewHit().

std::unique_ptr<TrackingSlaveSD> TotemRPSD::slave_
private

Definition at line 56 of file TotemRPSD.h.

Referenced by clearHits(), EndOfEvent(), fillHits(), and TotemRPSD().

TotemRPG4HitCollection* TotemRPSD::theHC_
private

Definition at line 83 of file TotemRPSD.h.

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

G4ThreeVector TotemRPSD::theLocalEntryPoint_
private

Definition at line 96 of file TotemRPSD.h.

Referenced by createNewHit(), and stepInfo().

G4ThreeVector TotemRPSD::theLocalExitPoint_
private

Definition at line 97 of file TotemRPSD.h.

Referenced by createNewHit(), and stepInfo().

double TotemRPSD::thePx_
private

Definition at line 100 of file TotemRPSD.h.

Referenced by createNewHit(), and stepInfo().

double TotemRPSD::thePy_
private

Definition at line 100 of file TotemRPSD.h.

Referenced by createNewHit(), and stepInfo().

double TotemRPSD::thePz_
private

Definition at line 100 of file TotemRPSD.h.

Referenced by createNewHit(), and stepInfo().

double TotemRPSD::ThetaAtEntry_
private

Definition at line 105 of file TotemRPSD.h.

Referenced by createNewHit(), and stepInfo().

G4Track* TotemRPSD::theTrack_
private

Definition at line 86 of file TotemRPSD.h.

Referenced by printHitInfo(), and stepInfo().

double TotemRPSD::Tof_
private

Definition at line 101 of file TotemRPSD.h.

Referenced by createNewHit(), and stepInfo().

G4double TotemRPSD::tSlice_
private

Definition at line 90 of file TotemRPSD.h.

Referenced by createNewHit(), and stepInfo().

G4int TotemRPSD::tSliceID_
private

Definition at line 89 of file TotemRPSD.h.

Referenced by printHitInfo(), and stepInfo().

unsigned int TotemRPSD::unitID_
private

Definition at line 88 of file TotemRPSD.h.

Referenced by createNewHit(), printHitInfo(), and stepInfo().

int TotemRPSD::verbosity_
private

Definition at line 60 of file TotemRPSD.h.

Referenced by stepInfo(), storeHit(), and TotemRPSD().

double TotemRPSD::Vx_
private

Definition at line 109 of file TotemRPSD.h.

Referenced by createNewHit(), and stepInfo().

double TotemRPSD::Vy_
private

Definition at line 109 of file TotemRPSD.h.

Referenced by createNewHit(), and stepInfo().

double TotemRPSD::Vz_
private

Definition at line 109 of file TotemRPSD.h.

Referenced by createNewHit(), and stepInfo().