LArSoft  v09_90_00
Liquid Argon Software toolkit - https://larsoft.org/
larg4::SimEnergyDepositSD Class Reference

#include "SimEnergyDepositSD.h"

Inheritance diagram for larg4::SimEnergyDepositSD:

Public Member Functions

 SimEnergyDepositSD (G4String)
 
 ~SimEnergyDepositSD ()
 
void Initialize (G4HCofThisEvent *)
 
G4bool ProcessHits (G4Step *, G4TouchableHistory *)
 
const sim::SimEnergyDepositCollectionGetHits () const
 

Private Attributes

sim::SimEnergyDepositCollection hitCollection
 

Detailed Description

Definition at line 19 of file SimEnergyDepositSD.h.

Constructor & Destructor Documentation

larg4::SimEnergyDepositSD::SimEnergyDepositSD ( G4String  name)

Definition at line 32 of file SimEnergyDepositSD.cc.

References hitCollection.

32  : G4VSensitiveDetector(name)
33  {
34  hitCollection.clear();
35  }
sim::SimEnergyDepositCollection hitCollection
larg4::SimEnergyDepositSD::~SimEnergyDepositSD ( )

Definition at line 39 of file SimEnergyDepositSD.cc.

39 {}

Member Function Documentation

const sim::SimEnergyDepositCollection& larg4::SimEnergyDepositSD::GetHits ( ) const
inline

Definition at line 25 of file SimEnergyDepositSD.h.

References hitCollection.

Referenced by larg4::LArG4DetectorService::doFillEventWithArtHits().

25 { return hitCollection; }
sim::SimEnergyDepositCollection hitCollection
void larg4::SimEnergyDepositSD::Initialize ( G4HCofThisEvent HCE)

Definition at line 43 of file SimEnergyDepositSD.cc.

References hitCollection.

44  {
45  hitCollection.clear();
46  }
sim::SimEnergyDepositCollection hitCollection
G4bool larg4::SimEnergyDepositSD::ProcessHits ( G4Step *  aStep,
G4TouchableHistory *   
)

Definition at line 49 of file SimEnergyDepositSD.cc.

References edep, util::end(), and hitCollection.

50  {
51  G4double edep = aStep->GetTotalEnergyDeposit() / CLHEP::MeV;
52 
53  if (edep == 0.) return false;
54  //std::cout << "7777777777777777: "<< aStep->GetTotalEnergyDeposit()/CLHEP::MeV << " " << aStep->GetTotalEnergyDeposit() <<std::endl;
55  const int electronsperMeV = 10000;
56  int nrelec = (int)round(edep * electronsperMeV);
57  if (aStep->GetTrack()->GetDynamicParticle()->GetCharge() == 0) return false;
58  G4int photons = 0;
59  G4SteppingManager* fpSteppingManager =
60  G4EventManager::GetEventManager()->GetTrackingManager()->GetSteppingManager();
61  G4StepStatus stepStatus = fpSteppingManager->GetfStepStatus();
62  if (stepStatus != fAtRestDoItProc) {
63  G4ProcessVector* procPost = fpSteppingManager->GetfPostStepDoItVector();
64  size_t MAXofPostStepLoops = fpSteppingManager->GetMAXofPostStepLoops();
65  for (size_t i3 = 0; i3 < MAXofPostStepLoops; i3++) {
66  if (!(*procPost)[i3]) continue;
67  /*
68  if ((*procPost)[i3]->GetProcessName() == "Cerenkov") {
69  G4Cerenkov* proc =(G4Cerenkov*) (*procPost)[i3];
70  photons+=proc->GetNumPhotons();
71  }
72  */
73  if ((*procPost)[i3]->GetProcessName() == "Scintillation") {
74  G4Scintillation* proc1 = (G4Scintillation*)(*procPost)[i3];
75  photons += proc1->GetNumPhotons();
76  }
77  }
78  }
79  geo::Point_t start = geo::Point_t(aStep->GetPreStepPoint()->GetPosition().x() / CLHEP::cm,
80  aStep->GetPreStepPoint()->GetPosition().y() / CLHEP::cm,
81  aStep->GetPreStepPoint()->GetPosition().z() / CLHEP::cm);
82  geo::Point_t end = geo::Point_t(aStep->GetPostStepPoint()->GetPosition().x() / CLHEP::cm,
83  aStep->GetPostStepPoint()->GetPosition().y() / CLHEP::cm,
84  aStep->GetPostStepPoint()->GetPosition().z() / CLHEP::cm);
85  sim::SimEnergyDeposit newHit =
86  sim::SimEnergyDeposit(photons,
87  nrelec,
88  1.0,
89  edep,
90  start,
91  end,
92  aStep->GetPreStepPoint()->GetGlobalTime() / CLHEP::ns,
93  aStep->GetPostStepPoint()->GetGlobalTime() / CLHEP::ns,
94  aStep->GetTrack()->GetTrackID(),
95  aStep->GetTrack()->GetParticleDefinition()->GetPDGEncoding(),
96  aStep->GetTrack()->GetTrackID() //original track id
97  );
98  hitCollection.push_back(newHit);
99  return true;
100  } // end ProcessHits
sim::SimEnergyDepositCollection hitCollection
decltype(auto) constexpr end(T &&obj)
ADL-aware version of std::end.
Definition: StdUtils.h:77
Double_t edep
Definition: macro.C:13
ROOT::Math::PositionVector3D< ROOT::Math::Cartesian3D< double >, ROOT::Math::GlobalCoordinateSystemTag > Point_t
Type for representation of position in physical 3D space.
Definition: geo_vectors.h:180
Energy deposition in the active material.

Member Data Documentation

sim::SimEnergyDepositCollection larg4::SimEnergyDepositSD::hitCollection
private

Definition at line 28 of file SimEnergyDepositSD.h.

Referenced by GetHits(), Initialize(), ProcessHits(), and SimEnergyDepositSD().


The documentation for this class was generated from the following files: