CMS 3D CMS Logo

List of all members | Public Member Functions | Protected Attributes | Private Member Functions | Private Attributes
BHMSD Class Reference

#include <BHMSD.h>

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

Public Member Functions

 BHMSD (std::string, const DDCompactView &, const SensitiveDetectorCatalog &, edm::ParameterSet const &, const SimTrackManager *)
 
virtual void clear ()
 
virtual void DrawAll ()
 
virtual void EndOfEvent (G4HCofThisEvent *eventHC)
 
void fillHits (edm::PSimHitContainer &, std::string use)
 
virtual double getEnergyDeposit (G4Step *step)
 
std::vector< std::string > getNames ()
 
virtual void Initialize (G4HCofThisEvent *HCE)
 
virtual void PrintAll ()
 
virtual bool ProcessHits (G4Step *, G4TouchableHistory *)
 
virtual uint32_t setDetUnitId (G4Step *)
 
virtual ~BHMSD ()
 
- Public Member Functions inherited from SensitiveTkDetector
 SensitiveTkDetector (std::string &iname, const DDCompactView &cpv, const SensitiveDetectorCatalog &clg, edm::ParameterSet const &p)
 
- Public Member Functions inherited from SensitiveDetector
virtual void AssignSD (const std::string &vname)
 
Local3DPoint ConvertToLocal3DPoint (const G4ThreeVector &point)
 
Local3DPoint FinalStepPosition (G4Step *s, coordinates)
 
Local3DPoint InitialStepPosition (G4Step *s, coordinates)
 
std::string nameOfSD ()
 
void NaNTrap (G4Step *step)
 
void Register ()
 
 SensitiveDetector (std::string &iname, const DDCompactView &cpv, const SensitiveDetectorCatalog &, edm::ParameterSet const &p)
 
virtual ~SensitiveDetector ()
 
- 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

virtual void clearHits ()
 
void CreateNewHit ()
 
void GetStepInfo (G4Step *aStep)
 
G4bool HitExists ()
 
void ResetForNewPrimary ()
 
G4ThreeVector SetToLocal (const G4ThreeVector &global)
 
G4ThreeVector SetToLocalExit (const G4ThreeVector &globalPoint)
 
void StoreHit (BscG4Hit *)
 
void Summarize ()
 
void update (const BeginOfRun *)
 This routine will be called when the appropriate signal arrives. More...
 
void update (const BeginOfEvent *)
 This routine will be called when the appropriate signal arrives. More...
 
void update (const ::EndOfEvent *)
 
void UpdateHit ()
 

Private Attributes

BscG4HitcurrentHit
 
G4VPhysicalVolume * currentPV
 
float edeposit
 
float Eloss
 
G4ThreeVector entrancePoint
 
int eventno
 
G4ThreeVector exitPoint
 
G4int hcID
 
G4ThreeVector hitPoint
 
G4ThreeVector hitPointExit
 
G4ThreeVector hitPointLocal
 
G4ThreeVector hitPointLocalExit
 
float incidentEnergy
 
std::string name
 
BHMNumberingSchemenumberingScheme
 
float Pabs
 
int ParentId
 
short ParticleType
 
float PhiAtEntry
 
G4StepPoint * postStepPoint
 
G4StepPoint * preStepPoint
 
uint32_t previousUnitID
 
G4int primaryID
 
G4int primID
 
TrackingSlaveSDslave
 
G4ThreeVector theEntryPoint
 
G4ThreeVector theExitPoint
 
BscG4HitCollectiontheHC
 
const SimTrackManagertheManager
 
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

- Public Types inherited from SensitiveDetector
enum  coordinates { WorldCoordinates, LocalCoordinates }
 
- 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 35 of file BHMSD.h.

Constructor & Destructor Documentation

BHMSD::BHMSD ( std::string  name,
const DDCompactView cpv,
const SensitiveDetectorCatalog clg,
edm::ParameterSet const &  p,
const SimTrackManager manager 
)

Definition at line 30 of file BHMSD.cc.

References SensitiveDetector::AssignSD(), ecaldqm::collectionName, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), LogDebug, SensitiveDetectorCatalog::logicalNames(), numberingScheme, SensitiveDetector::Register(), and slave.

32  :
34  hcID(-1), theHC(0), theManager(manager), currentHit(0), theTrack(0),
36  postStepPoint(0), eventno(0){
37 
38  //Add BHM Sentitive Detector Name
39  collectionName.insert(name);
40 
41 
42  //Parameters
43  edm::ParameterSet m_p = p.getParameter<edm::ParameterSet>("BHMSD");
44  int verbn = m_p.getUntrackedParameter<int>("Verbosity");
45  //int verbn = 1;
46 
47  SetVerboseLevel(verbn);
48  LogDebug("BHMSim")
49  << "*******************************************************\n"
50  << "* *\n"
51  << "* Constructing a BHMSD with name " << name << "\n"
52  << "* *\n"
53  << "*******************************************************";
54 
55 
56  slave = new TrackingSlaveSD(name);
57 
58  //
59  // attach detectors (LogicalVolumes)
60  //
61  std::vector<std::string> lvNames = clg.logicalNames(name);
62 
63  this->Register();
64 
65  for (std::vector<std::string>::iterator it=lvNames.begin();
66  it !=lvNames.end(); it++) {
67  this->AssignSD(*it);
68  edm::LogInfo("BHMSim") << "BHMSD : Assigns SD to LV " << (*it);
69  }
70 
71  if (verbn > 0) {
72  edm::LogInfo("BHMSim") << "name = " <<name <<" and new BHMNumberingScheme";
73  }
75 
76  edm::LogInfo("BHMSim") << "BHMSD: Instantiation completed";
77 }
#define LogDebug(id)
T getUntrackedParameter(std::string const &, T const &) const
std::string name
Definition: BHMSD.h:92
BHMNumberingScheme * numberingScheme
Definition: BHMSD.h:84
G4StepPoint * preStepPoint
Definition: BHMSD.h:105
int eventno
Definition: BHMSD.h:122
const std::vector< std::string > & logicalNames(const std::string &readoutName) const
virtual void AssignSD(const std::string &vname)
TrackingSlaveSD * slave
Definition: BHMSD.h:83
std::string const collectionName[nCollections]
Definition: Collections.h:45
uint32_t previousUnitID
Definition: BHMSD.h:101
uint32_t unitID
Definition: BHMSD.h:101
G4Track * theTrack
Definition: BHMSD.h:99
const SimTrackManager * theManager
Definition: BHMSD.h:95
G4StepPoint * postStepPoint
Definition: BHMSD.h:106
SensitiveTkDetector(std::string &iname, const DDCompactView &cpv, const SensitiveDetectorCatalog &clg, edm::ParameterSet const &p)
BscG4Hit * currentHit
Definition: BHMSD.h:98
BscG4HitCollection * theHC
Definition: BHMSD.h:94
G4VPhysicalVolume * currentPV
Definition: BHMSD.h:100
G4int hcID
Definition: BHMSD.h:93
BHMSD::~BHMSD ( )
virtual

Definition at line 80 of file BHMSD.cc.

References numberingScheme, and slave.

80  {
81 
82  if (slave) delete slave;
83  if (numberingScheme) delete numberingScheme;
84 }
BHMNumberingScheme * numberingScheme
Definition: BHMSD.h:84
TrackingSlaveSD * slave
Definition: BHMSD.h:83

Member Function Documentation

void BHMSD::clear ( void  )
virtual

Definition at line 360 of file BHMSD.cc.

360  {
361 }
void BHMSD::clearHits ( )
privatevirtual

Implements SensitiveDetector.

Definition at line 398 of file BHMSD.cc.

References TrackingSlaveSD::Initialize(), and slave.

Referenced by update().

398  {
399  slave->Initialize();
400 }
TrackingSlaveSD * slave
Definition: BHMSD.h:83
virtual void Initialize()
void BHMSD::CreateNewHit ( )
private

Definition at line 234 of file BHMSD.cc.

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

Referenced by ProcessHits().

234  {
235 
236 #ifdef debug
237  LogDebug("BHMSim") << "BHMSD CreateNewHit for"
238  << " PV " << currentPV->GetName()
239  << " PVid = " << currentPV->GetCopyNo()
240  << " Unit " << unitID <<std::endl;
241  LogDebug("BHMSim") << " primary " << primaryID
242  << " time slice " << tSliceID
243  << " For Track " << theTrack->GetTrackID()
244  << " which is a " << theTrack->GetDefinition()->GetParticleName();
245 
246  if (theTrack->GetTrackID()==1) {
247  LogDebug("BHMSim") << " of energy " << theTrack->GetTotalEnergy();
248  } else {
249  LogDebug("BHMSim") << " daughter of part. " << theTrack->GetParentID();
250  }
251 
252  LogDebug("BHMSim") << " and created by " ;
253  if (theTrack->GetCreatorProcess()!=NULL)
254  LogDebug("BHMSim") << theTrack->GetCreatorProcess()->GetProcessName() ;
255  else
256  LogDebug("BHMSim") << "NO process";
257  LogDebug("BHMSim") << std::endl;
258 #endif
259 
260  currentHit = new BscG4Hit;
265 
272 
274 
277 
279  currentHit->setVx(Vx);
280  currentHit->setVy(Vy);
281  currentHit->setVz(Vz);
282 
283  currentHit->setX(X);
284  currentHit->setY(Y);
285  currentHit->setZ(Z);
286 
287  UpdateHit();
288 
290 }
#define LogDebug(id)
float Z
Definition: BHMSD.h:120
void setTof(float e)
Definition: BSCG4Hit.cc:154
float incidentEnergy
Definition: BHMSD.h:89
G4double tSlice
Definition: BHMSD.h:103
void setEntry(const G4ThreeVector &xyz)
Definition: BSCG4Hit.cc:116
float PhiAtEntry
Definition: BHMSD.h:116
void setVx(float p)
Definition: BSCG4Hit.cc:177
float Eloss
Definition: BHMSD.h:114
void setTrackID(int i)
Definition: BSCG4Hit.cc:134
void setEnergyLoss(float e)
Definition: BSCG4Hit.cc:155
void setParentId(int p)
Definition: BSCG4Hit.cc:174
G4ThreeVector hitPoint
Definition: BHMSD.h:109
#define NULL
Definition: scimark2.h:8
void setUnitID(unsigned int i)
Definition: BSCG4Hit.cc:137
G4int tSliceID
Definition: BHMSD.h:102
G4ThreeVector hitPointLocal
Definition: BHMSD.h:111
void setVy(float p)
Definition: BSCG4Hit.cc:180
void setY(float t)
Definition: BSCG4Hit.cc:168
float X
Definition: BHMSD.h:120
void setThetaAtEntry(float t)
Definition: BSCG4Hit.cc:161
void setX(float t)
Definition: BSCG4Hit.cc:165
void setPabs(float e)
Definition: BSCG4Hit.cc:153
void setEntryLocalP(const G4ThreeVector &xyz1)
Definition: BSCG4Hit.cc:119
float Vx
Definition: BHMSD.h:119
float Vz
Definition: BHMSD.h:119
float Y
Definition: BHMSD.h:120
void setPhiAtEntry(float f)
Definition: BSCG4Hit.cc:162
uint32_t unitID
Definition: BHMSD.h:101
float Vy
Definition: BHMSD.h:119
float Tof
Definition: BHMSD.h:114
G4Track * theTrack
Definition: BHMSD.h:99
void setTimeSlice(double d)
Definition: BSCG4Hit.cc:140
void setParticleType(short i)
Definition: BSCG4Hit.cc:156
void setVz(float p)
Definition: BSCG4Hit.cc:183
short ParticleType
Definition: BHMSD.h:115
G4int primaryID
Definition: BHMSD.h:102
BscG4Hit * currentHit
Definition: BHMSD.h:98
int ParentId
Definition: BHMSD.h:118
void UpdateHit()
Definition: BHMSD.cc:293
float ThetaAtEntry
Definition: BHMSD.h:116
void setZ(float t)
Definition: BSCG4Hit.cc:171
G4ThreeVector hitPointLocalExit
Definition: BHMSD.h:112
G4VPhysicalVolume * currentPV
Definition: BHMSD.h:100
void setExitLocalP(const G4ThreeVector &xyz1)
Definition: BSCG4Hit.cc:122
float Pabs
Definition: BHMSD.h:114
void setIncidentEnergy(double e)
Definition: BSCG4Hit.cc:131
void StoreHit(BscG4Hit *)
Definition: BHMSD.cc:222
void BHMSD::DrawAll ( )
virtual

Definition at line 364 of file BHMSD.cc.

364  {
365 }
void BHMSD::EndOfEvent ( G4HCofThisEvent *  eventHC)
virtual

Reimplemented from SensitiveDetector.

Definition at line 329 of file BHMSD.cc.

References BscG4Hit::getEnergyLoss(), BscG4Hit::getEntry(), BscG4Hit::getPabs(), BscG4Hit::getParticleType(), BscG4Hit::getPhiAtEntry(), BscG4Hit::getThetaAtEntry(), BscG4Hit::getTof(), BscG4Hit::getTrackID(), BscG4Hit::getUnitID(), LogDebug, TrackingSlaveSD::processHits(), slave, Summarize(), and theHC.

329  {
330 
331  // here we loop over transient hits and make them persistent
332  for (int j=0; j<theHC->entries(); j++) {
333  BscG4Hit* aHit = (*theHC)[j];
334  LogDebug("BHMSim") << "hit number" << j << "unit ID = "<<aHit->getUnitID()<< "\n";
335  LogDebug("BHMSim") << "entry z " << aHit->getEntry().z()<< "\n";
336  LogDebug("BHMSim") << "entr theta " << aHit->getThetaAtEntry()<< "\n";
337 
338  Local3DPoint locExitPoint(0,0,0);
339  Local3DPoint locEntryPoint(aHit->getEntry().x(),
340  aHit->getEntry().y(),
341  aHit->getEntry().z());
342  slave->processHits(PSimHit(locEntryPoint,locExitPoint,
343  aHit->getPabs(),
344  aHit->getTof(),
345  aHit->getEnergyLoss(),
346  aHit->getParticleType(),
347  aHit->getUnitID(),
348  aHit->getTrackID(),
349  aHit->getThetaAtEntry(),
350  aHit->getPhiAtEntry()));
351  }
352  Summarize();
353 }
#define LogDebug(id)
float getPhiAtEntry() const
Definition: BSCG4Hit.cc:159
int getParticleType() const
Definition: BSCG4Hit.cc:151
G4int getTrackID() const
Definition: BSCG4Hit.cc:133
float getTof() const
Definition: BSCG4Hit.cc:149
float getPabs() const
Definition: BSCG4Hit.cc:148
unsigned int getUnitID() const
Definition: BSCG4Hit.cc:136
TrackingSlaveSD * slave
Definition: BHMSD.h:83
void Summarize()
Definition: BHMSD.cc:356
float getThetaAtEntry() const
Definition: BSCG4Hit.cc:158
BscG4HitCollection * theHC
Definition: BHMSD.h:94
G4ThreeVector getEntry() const
Definition: BSCG4Hit.cc:115
virtual bool processHits(const PSimHit &)
float getEnergyLoss() const
Definition: BSCG4Hit.cc:150
void BHMSD::fillHits ( edm::PSimHitContainer c,
std::string  use 
)
virtual

Implements SensitiveTkDetector.

Definition at line 374 of file BHMSD.cc.

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

374  {
375  if (slave->name() == n) c=slave->hits();
376 }
std::string name() const
TrackingSlaveSD * slave
Definition: BHMSD.h:83
std::vector< PSimHit > & hits()
double BHMSD::getEnergyDeposit ( G4Step *  step)
virtual

Definition at line 86 of file BHMSD.cc.

Referenced by GetStepInfo().

86  {
87  return aStep->GetTotalEnergyDeposit();
88 }
std::vector< std::string > BHMSD::getNames ( )
virtual

Reimplemented from SensitiveDetector.

Definition at line 402 of file BHMSD.cc.

References TrackingSlaveSD::name(), slave, and groupFilesInBlocks::temp.

402  {
403  std::vector<std::string> temp;
404  temp.push_back(slave->name());
405  return temp;
406 }
std::string name() const
TrackingSlaveSD * slave
Definition: BHMSD.h:83
void BHMSD::GetStepInfo ( G4Step *  aStep)
private

Definition at line 122 of file BHMSD.cc.

References currentPV, edeposit, edepositEM, edepositHAD, Eloss, emPDG, epPDG, gammaPDG, getEnergyDeposit(), GeV, 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().

122  {
123 
124  preStepPoint = aStep->GetPreStepPoint();
125  postStepPoint= aStep->GetPostStepPoint();
126  theTrack = aStep->GetTrack();
127  hitPoint = preStepPoint->GetPosition();
128  currentPV = preStepPoint->GetPhysicalVolume();
129  hitPointExit = postStepPoint->GetPosition();
130 
131  hitPointLocal = preStepPoint->GetTouchable()->GetHistory()->GetTopTransform().TransformPoint(hitPoint);
132  hitPointLocalExit = preStepPoint->GetTouchable()->GetHistory()->GetTopTransform().TransformPoint(hitPointExit);
133 
134 
135  G4int particleCode = theTrack->GetDefinition()->GetPDGEncoding();
136  LogDebug("BHMSim") << " BHMSD :particleType = " << theTrack->GetDefinition()->GetParticleName() <<std::endl;
137  if (particleCode == emPDG ||
138  particleCode == epPDG ||
139  particleCode == gammaPDG ) {
140  edepositEM = getEnergyDeposit(aStep); edepositHAD = 0.;
141  } else {
142  edepositEM = 0.; edepositHAD = getEnergyDeposit(aStep);
143  }
144  edeposit = aStep->GetTotalEnergyDeposit();
145  tSlice = (postStepPoint->GetGlobalTime() )/CLHEP::nanosecond;
146  tSliceID = (int) tSlice;
147  unitID = setDetUnitId(aStep);
148 #ifdef debug
149  LogDebug("BHMSim") << "unitID=" << unitID <<std::endl;
150 #endif
151  primaryID = theTrack->GetTrackID();
152  // Position = hitPoint;
153  Pabs = aStep->GetPreStepPoint()->GetMomentum().mag()/CLHEP::GeV;
154  Tof = aStep->GetPostStepPoint()->GetGlobalTime()/CLHEP::nanosecond;
155  Eloss = aStep->GetTotalEnergyDeposit()/CLHEP::GeV;
156  ParticleType = theTrack->GetDefinition()->GetPDGEncoding();
157  ThetaAtEntry = aStep->GetPreStepPoint()->GetPosition().theta()/CLHEP::deg;
158  PhiAtEntry = aStep->GetPreStepPoint()->GetPosition().phi()/CLHEP::deg;
159 
160  ParentId = theTrack->GetParentID();
161  Vx = theTrack->GetVertexPosition().x();
162  Vy = theTrack->GetVertexPosition().y();
163  Vz = theTrack->GetVertexPosition().z();
164  X = hitPoint.x();
165  Y = hitPoint.y();
166  Z = hitPoint.z();
167 }
#define LogDebug(id)
float Z
Definition: BHMSD.h:120
G4double tSlice
Definition: BHMSD.h:103
G4ThreeVector hitPointExit
Definition: BHMSD.h:110
const double GeV
Definition: MathUtil.h:16
float PhiAtEntry
Definition: BHMSD.h:116
float Eloss
Definition: BHMSD.h:114
G4StepPoint * preStepPoint
Definition: BHMSD.h:105
G4ThreeVector hitPoint
Definition: BHMSD.h:109
G4int tSliceID
Definition: BHMSD.h:102
G4ThreeVector hitPointLocal
Definition: BHMSD.h:111
virtual double getEnergyDeposit(G4Step *step)
Definition: BHMSD.cc:86
float X
Definition: BHMSD.h:120
float edepositEM
Definition: BHMSD.h:126
float edepositHAD
Definition: BHMSD.h:126
float Vx
Definition: BHMSD.h:119
float Vz
Definition: BHMSD.h:119
float Y
Definition: BHMSD.h:120
uint32_t unitID
Definition: BHMSD.h:101
float Vy
Definition: BHMSD.h:119
float Tof
Definition: BHMSD.h:114
G4Track * theTrack
Definition: BHMSD.h:99
float edeposit
Definition: BHMSD.h:107
G4StepPoint * postStepPoint
Definition: BHMSD.h:106
short ParticleType
Definition: BHMSD.h:115
G4int primaryID
Definition: BHMSD.h:102
G4int gammaPDG
Definition: BHMSD.h:129
int ParentId
Definition: BHMSD.h:118
float ThetaAtEntry
Definition: BHMSD.h:116
G4ThreeVector hitPointLocalExit
Definition: BHMSD.h:112
G4VPhysicalVolume * currentPV
Definition: BHMSD.h:100
float Pabs
Definition: BHMSD.h:114
G4int epPDG
Definition: BHMSD.h:128
virtual uint32_t setDetUnitId(G4Step *)
Definition: BHMSD.cc:169
G4int emPDG
Definition: BHMSD.h:127
G4bool BHMSD::HitExists ( )
private

Definition at line 174 of file BHMSD.cc.

References currentHit, runEdmFileComparison::found, BscG4Hit::getTimeSliceID(), BscG4Hit::getTrackID(), BscG4Hit::getUnitID(), previousUnitID, primaryID, primID, ResetForNewPrimary(), theHC, tsID, tSliceID, unitID, and UpdateHit().

Referenced by ProcessHits().

174  {
175  if (primaryID<1) {
176  edm::LogWarning("BHMSim") << "***** BHMSD error: primaryID = "
177  << primaryID
178  << " maybe detector name changed";
179  }
180 
181  // Update if in the same detector, time-slice and for same track
182  if (tSliceID == tsID && unitID==previousUnitID) {
183  UpdateHit();
184  return true;
185  }
186  // Reset entry point for new primary
187  if (primaryID != primID)
189 
190  //look in the HitContainer whether a hit with the same primID, unitID,
191  //tSliceID already exists:
192 
193  G4bool found = false;
194 
195  for (int j=0; j<theHC->entries()&&!found; j++) {
196  BscG4Hit* aPreviousHit = (*theHC)[j];
197  if (aPreviousHit->getTrackID() == primaryID &&
198  aPreviousHit->getTimeSliceID() == tSliceID &&
199  aPreviousHit->getUnitID() == unitID ) {
200  currentHit = aPreviousHit;
201  found = true;
202  }
203  }
204 
205  if (found) {
206  UpdateHit();
207  return true;
208  } else {
209  return false;
210  }
211 }
void ResetForNewPrimary()
Definition: BHMSD.cc:214
G4int getTrackID() const
Definition: BSCG4Hit.cc:133
G4int tSliceID
Definition: BHMSD.h:102
unsigned int getUnitID() const
Definition: BSCG4Hit.cc:136
int getTimeSliceID() const
Definition: BSCG4Hit.cc:141
uint32_t previousUnitID
Definition: BHMSD.h:101
G4int tsID
Definition: BHMSD.h:97
uint32_t unitID
Definition: BHMSD.h:101
G4int primaryID
Definition: BHMSD.h:102
G4int primID
Definition: BHMSD.h:90
BscG4Hit * currentHit
Definition: BHMSD.h:98
BscG4HitCollection * theHC
Definition: BHMSD.h:94
void UpdateHit()
Definition: BHMSD.cc:293
void BHMSD::Initialize ( G4HCofThisEvent *  HCE)
virtual

Reimplemented from SensitiveDetector.

Definition at line 90 of file BHMSD.cc.

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

90  {
91 #ifdef debug
92  LogDebug("BHMSim") << "BHMSD : Initialize called for " << name << std::endl;
93 #endif
94 
96  if (hcID<0)
97  hcID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]);
98  HCE->AddHitsCollection(hcID, theHC);
99 
100  tsID = -2;
101  primID = -2;
102 }
#define LogDebug(id)
std::string name
Definition: BHMSD.h:92
std::string const collectionName[nCollections]
Definition: Collections.h:45
G4int tsID
Definition: BHMSD.h:97
G4int primID
Definition: BHMSD.h:90
BscG4HitCollection * theHC
Definition: BHMSD.h:94
G4int hcID
Definition: BHMSD.h:93
G4THitsCollection< BscG4Hit > BscG4HitCollection
void BHMSD::PrintAll ( )
virtual

Definition at line 368 of file BHMSD.cc.

References LogDebug, and theHC.

368  {
369  LogDebug("BHMSim") << "BHMSD: Collection " << theHC->GetName() << "\n";
370  theHC->PrintAllHits();
371 }
#define LogDebug(id)
BscG4HitCollection * theHC
Definition: BHMSD.h:94
bool BHMSD::ProcessHits ( G4Step *  aStep,
G4TouchableHistory *   
)
virtual

Implements SensitiveDetector.

Definition at line 105 of file BHMSD.cc.

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

105  {
106 
107  if (aStep == NULL) {
108  return true;
109  } else {
110  GetStepInfo(aStep);
111 #ifdef debug
112  LogDebug("BHMSim") << "BHMSD : number of hits = " << theHC->entries() << std::endl;
113 #endif
114  if (HitExists() == false && edeposit>0. ){
115  CreateNewHit();
116  return true;
117  }
118  }
119  return true;
120 }
#define LogDebug(id)
#define NULL
Definition: scimark2.h:8
void CreateNewHit()
Definition: BHMSD.cc:234
float edeposit
Definition: BHMSD.h:107
BscG4HitCollection * theHC
Definition: BHMSD.h:94
void GetStepInfo(G4Step *aStep)
Definition: BHMSD.cc:122
G4bool HitExists()
Definition: BHMSD.cc:174
void BHMSD::ResetForNewPrimary ( )
private

Definition at line 214 of file BHMSD.cc.

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

Referenced by HitExists().

214  {
215 
218  incidentEnergy = preStepPoint->GetKineticEnergy();
219 }
float incidentEnergy
Definition: BHMSD.h:89
G4ThreeVector hitPointExit
Definition: BHMSD.h:110
G4StepPoint * preStepPoint
Definition: BHMSD.h:105
G4ThreeVector hitPoint
Definition: BHMSD.h:109
G4ThreeVector SetToLocalExit(const G4ThreeVector &globalPoint)
Definition: BHMSD.cc:321
G4ThreeVector entrancePoint
Definition: BHMSD.h:86
G4ThreeVector SetToLocal(const G4ThreeVector &global)
Definition: BHMSD.cc:313
G4ThreeVector exitPoint
Definition: BHMSD.h:86
uint32_t BHMSD::setDetUnitId ( G4Step *  aStep)
virtual

Implements SensitiveDetector.

Definition at line 169 of file BHMSD.cc.

References BHMNumberingScheme::getUnitID(), and numberingScheme.

Referenced by GetStepInfo().

169  {
170  return (numberingScheme == 0 ? 0 : numberingScheme->getUnitID(aStep));
171 }
BHMNumberingScheme * numberingScheme
Definition: BHMSD.h:84
virtual unsigned int getUnitID(const G4Step *aStep) const
G4ThreeVector BHMSD::SetToLocal ( const G4ThreeVector &  global)
private

Definition at line 313 of file BHMSD.cc.

References preStepPoint, and theEntryPoint.

Referenced by ResetForNewPrimary().

313  {
314 
315  const G4VTouchable* touch= preStepPoint->GetTouchable();
316  theEntryPoint = touch->GetHistory()->GetTopTransform().TransformPoint(global);
317  return theEntryPoint;
318 }
G4ThreeVector theEntryPoint
Definition: BHMSD.h:87
G4StepPoint * preStepPoint
Definition: BHMSD.h:105
G4ThreeVector BHMSD::SetToLocalExit ( const G4ThreeVector &  globalPoint)
private

Definition at line 321 of file BHMSD.cc.

References postStepPoint, and theExitPoint.

Referenced by ResetForNewPrimary().

321  {
322 
323  const G4VTouchable* touch= postStepPoint->GetTouchable();
324  theExitPoint = touch->GetHistory()->GetTopTransform().TransformPoint(globalPoint);
325  return theExitPoint;
326 }
G4ThreeVector theExitPoint
Definition: BHMSD.h:87
G4StepPoint * postStepPoint
Definition: BHMSD.h:106
void BHMSD::StoreHit ( BscG4Hit hit)
private

Definition at line 222 of file BHMSD.cc.

References primID, and theHC.

Referenced by CreateNewHit().

222  {
223 
224  if (primID<0) return;
225  if (hit == 0 ) {
226  edm::LogWarning("BHMSim") << "BHMSD: hit to be stored is NULL !!";
227  return;
228  }
229 
230  theHC->insert( hit );
231 }
G4int primID
Definition: BHMSD.h:90
BscG4HitCollection * theHC
Definition: BHMSD.h:94
void BHMSD::Summarize ( )
private

Definition at line 356 of file BHMSD.cc.

Referenced by EndOfEvent().

356  {
357 }
void BHMSD::update ( const BeginOfRun )
privatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 385 of file BHMSD.cc.

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

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

385  {
386 
387  G4ParticleTable * theParticleTable = G4ParticleTable::GetParticleTable();
388  G4String particleName;
389  emPDG = theParticleTable->FindParticle(particleName="e-")->GetPDGEncoding();
390  epPDG = theParticleTable->FindParticle(particleName="e+")->GetPDGEncoding();
391  gammaPDG = theParticleTable->FindParticle(particleName="gamma")->GetPDGEncoding();
392 
393 }
G4int gammaPDG
Definition: BHMSD.h:129
G4int epPDG
Definition: BHMSD.h:128
G4int emPDG
Definition: BHMSD.h:127
void BHMSD::update ( const BeginOfEvent )
privatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 378 of file BHMSD.cc.

References clearHits(), eventno, and LogDebug.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

378  {
379  LogDebug("BHMSim") << " Dispatched BeginOfEvent for " << GetName()
380  << " !" ;
381  clearHits();
382  eventno = (*i)()->GetEventID();
383 }
#define LogDebug(id)
int eventno
Definition: BHMSD.h:122
virtual void clearHits()
Definition: BHMSD.cc:398
void BHMSD::update ( const ::EndOfEvent )
private
void BHMSD::UpdateHit ( )
private

Definition at line 293 of file BHMSD.cc.

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

Referenced by CreateNewHit(), and HitExists().

293  {
294 
295  if (Eloss > 0.) {
297 
298 #ifdef debug
299  LogDebug("BHMSim") << "updateHit: add eloss " << Eloss <<std::endl;
300  LogDebug("BHMSim") << "CurrentHit=" << currentHit
301  << ", PostStepPoint=" << postStepPoint->GetPosition();
302 #endif
304  }
305 
306  // buffer for next steps:
307  tsID = tSliceID;
308  primID = primaryID;
310 }
#define LogDebug(id)
float Eloss
Definition: BHMSD.h:114
void setEnergyLoss(float e)
Definition: BSCG4Hit.cc:155
G4int tSliceID
Definition: BHMSD.h:102
float edepositEM
Definition: BHMSD.h:126
float edepositHAD
Definition: BHMSD.h:126
uint32_t previousUnitID
Definition: BHMSD.h:101
G4int tsID
Definition: BHMSD.h:97
uint32_t unitID
Definition: BHMSD.h:101
void addEnergyDeposit(double em, double hd)
Definition: BSCG4Hit.cc:143
G4StepPoint * postStepPoint
Definition: BHMSD.h:106
G4int primaryID
Definition: BHMSD.h:102
G4int primID
Definition: BHMSD.h:90
BscG4Hit * currentHit
Definition: BHMSD.h:98

Member Data Documentation

BscG4Hit* BHMSD::currentHit
private

Definition at line 98 of file BHMSD.h.

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

G4VPhysicalVolume* BHMSD::currentPV
private

Definition at line 100 of file BHMSD.h.

Referenced by CreateNewHit(), and GetStepInfo().

float BHMSD::edeposit
private

Definition at line 107 of file BHMSD.h.

Referenced by GetStepInfo(), and ProcessHits().

float BHMSD::edepositEM
protected

Definition at line 126 of file BHMSD.h.

Referenced by GetStepInfo(), and UpdateHit().

float BHMSD::edepositHAD
protected

Definition at line 126 of file BHMSD.h.

Referenced by GetStepInfo(), and UpdateHit().

float BHMSD::Eloss
private

Definition at line 114 of file BHMSD.h.

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

G4int BHMSD::emPDG
protected

Definition at line 127 of file BHMSD.h.

Referenced by GetStepInfo(), and update().

G4ThreeVector BHMSD::entrancePoint
private

Definition at line 86 of file BHMSD.h.

Referenced by ResetForNewPrimary().

G4int BHMSD::epPDG
protected

Definition at line 128 of file BHMSD.h.

Referenced by GetStepInfo(), and update().

int BHMSD::eventno
private

Definition at line 122 of file BHMSD.h.

Referenced by update().

G4ThreeVector BHMSD::exitPoint
private

Definition at line 86 of file BHMSD.h.

Referenced by ResetForNewPrimary().

G4int BHMSD::gammaPDG
protected

Definition at line 129 of file BHMSD.h.

Referenced by GetStepInfo(), and update().

G4int BHMSD::hcID
private

Definition at line 93 of file BHMSD.h.

Referenced by Initialize().

G4ThreeVector BHMSD::hitPoint
private

Definition at line 109 of file BHMSD.h.

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

G4ThreeVector BHMSD::hitPointExit
private

Definition at line 110 of file BHMSD.h.

Referenced by GetStepInfo(), and ResetForNewPrimary().

G4ThreeVector BHMSD::hitPointLocal
private

Definition at line 111 of file BHMSD.h.

Referenced by CreateNewHit(), and GetStepInfo().

G4ThreeVector BHMSD::hitPointLocalExit
private

Definition at line 112 of file BHMSD.h.

Referenced by CreateNewHit(), and GetStepInfo().

float BHMSD::incidentEnergy
private

Definition at line 89 of file BHMSD.h.

Referenced by CreateNewHit(), and ResetForNewPrimary().

std::string BHMSD::name
private
BHMNumberingScheme* BHMSD::numberingScheme
private

Definition at line 84 of file BHMSD.h.

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

float BHMSD::Pabs
private

Definition at line 114 of file BHMSD.h.

Referenced by CreateNewHit(), and GetStepInfo().

int BHMSD::ParentId
private

Definition at line 118 of file BHMSD.h.

Referenced by CreateNewHit(), and GetStepInfo().

short BHMSD::ParticleType
private

Definition at line 115 of file BHMSD.h.

Referenced by CreateNewHit(), and GetStepInfo().

float BHMSD::PhiAtEntry
private

Definition at line 116 of file BHMSD.h.

Referenced by CreateNewHit(), and GetStepInfo().

G4StepPoint* BHMSD::postStepPoint
private

Definition at line 106 of file BHMSD.h.

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

G4StepPoint* BHMSD::preStepPoint
private

Definition at line 105 of file BHMSD.h.

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

uint32_t BHMSD::previousUnitID
private

Definition at line 101 of file BHMSD.h.

Referenced by HitExists(), and UpdateHit().

G4int BHMSD::primaryID
private

Definition at line 102 of file BHMSD.h.

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

G4int BHMSD::primID
private

Definition at line 90 of file BHMSD.h.

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

TrackingSlaveSD* BHMSD::slave
private

Definition at line 83 of file BHMSD.h.

Referenced by BHMSD(), clearHits(), EndOfEvent(), fillHits(), getNames(), and ~BHMSD().

G4ThreeVector BHMSD::theEntryPoint
private

Definition at line 87 of file BHMSD.h.

Referenced by SetToLocal().

G4ThreeVector BHMSD::theExitPoint
private

Definition at line 87 of file BHMSD.h.

Referenced by SetToLocalExit().

BscG4HitCollection* BHMSD::theHC
private

Definition at line 94 of file BHMSD.h.

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

const SimTrackManager* BHMSD::theManager
private

Definition at line 95 of file BHMSD.h.

float BHMSD::ThetaAtEntry
private

Definition at line 116 of file BHMSD.h.

Referenced by CreateNewHit(), and GetStepInfo().

G4Track* BHMSD::theTrack
private

Definition at line 99 of file BHMSD.h.

Referenced by CreateNewHit(), and GetStepInfo().

float BHMSD::Tof
private

Definition at line 114 of file BHMSD.h.

Referenced by CreateNewHit(), and GetStepInfo().

G4int BHMSD::tsID
private

Definition at line 97 of file BHMSD.h.

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

G4double BHMSD::tSlice
private

Definition at line 103 of file BHMSD.h.

Referenced by CreateNewHit(), and GetStepInfo().

G4int BHMSD::tSliceID
private

Definition at line 102 of file BHMSD.h.

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

uint32_t BHMSD::unitID
private

Definition at line 101 of file BHMSD.h.

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

float BHMSD::Vx
private

Definition at line 119 of file BHMSD.h.

Referenced by CreateNewHit(), and GetStepInfo().

float BHMSD::Vy
private

Definition at line 119 of file BHMSD.h.

Referenced by CreateNewHit(), and GetStepInfo().

float BHMSD::Vz
private

Definition at line 119 of file BHMSD.h.

Referenced by CreateNewHit(), and GetStepInfo().

float BHMSD::X
private

Definition at line 120 of file BHMSD.h.

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

float BHMSD::Y
private

Definition at line 120 of file BHMSD.h.

Referenced by CreateNewHit(), and GetStepInfo().

float BHMSD::Z
private

Definition at line 120 of file BHMSD.h.

Referenced by CreateNewHit(), and GetStepInfo().