All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
SimEnergyDepositSD.cc
Go to the documentation of this file.
1 //
2 // __ __ __ __ __
3 // ____ ______/ /_____ _/ // / / /_/ /__
4 // / __ `/ ___/ __/ __ `/ // /_/ __/ //_/
5 // / /_/ / / / /_/ /_/ /__ __/ /_/ ,<
6 // \__,_/_/ \__/\__, / /_/ \__/_/|_|
7 // /____/
8 //
9 // artg4tk: art based Geant 4 Toolkit
10 //
11 //=============================================================================
12 // SimEnergyDeposit.cc: Class representing a sensitive tracking detector
13 // Author: Hans Wenzel (Fermilab)
14 //=============================================================================
15 #include "larg4/Services/SimEnergyDepositSD.h"
16 #include "Geant4/G4HCofThisEvent.hh"
17 #include "Geant4/G4Step.hh"
18 #include "Geant4/G4ThreeVector.hh"
19 #include "Geant4/G4SDManager.hh"
20 #include "Geant4/G4ios.hh"
21 #include "Geant4/G4VVisManager.hh"
22 #include "Geant4/G4Event.hh"
23 #include "Geant4/G4EventManager.hh"
24 #include "Geant4/G4VSolid.hh"
25 #include "Geant4/G4Cerenkov.hh"
26 #include "Geant4/G4Scintillation.hh"
27 #include "Geant4/G4SteppingManager.hh"
28 
29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
30 namespace larg4 {
31 
32  SimEnergyDepositSD::SimEnergyDepositSD(G4String name)
33 : G4VSensitiveDetector(name) {
34  hitCollection.clear();
35 }
36 
37 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
38 
39  SimEnergyDepositSD::~SimEnergyDepositSD() {
40 
41 }
42 
43 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
44 
45  void SimEnergyDepositSD::Initialize(G4HCofThisEvent* HCE) {
46  hitCollection.clear();
47  }
48  //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
49 
50  G4bool SimEnergyDepositSD::ProcessHits(G4Step* aStep, G4TouchableHistory*) {
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 = G4EventManager::GetEventManager()
60  ->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(
80  aStep->GetPreStepPoint()->GetPosition().x()/CLHEP::cm,
81  aStep->GetPreStepPoint()->GetPosition().y()/CLHEP::cm,
82  aStep->GetPreStepPoint()->GetPosition().z()/CLHEP::cm);
83  geo::Point_t end = geo::Point_t(
84  aStep->GetPostStepPoint()->GetPosition().x()/CLHEP::cm,
85  aStep->GetPostStepPoint()->GetPosition().y()/CLHEP::cm,
86  aStep->GetPostStepPoint()->GetPosition().z()/CLHEP::cm);
87  sim::SimEnergyDeposit newHit = sim::SimEnergyDeposit(photons,
88  nrelec,
89  1.0,
90  edep,
91  start,
92  end,
93  aStep->GetPreStepPoint()->GetGlobalTime() / CLHEP::ns,
94  aStep->GetPostStepPoint()->GetGlobalTime() /CLHEP::ns,
95  aStep->GetTrack()->GetTrackID(),
96  aStep->GetTrack()->GetParticleDefinition()->GetPDGEncoding() );
97  hitCollection.push_back(newHit);
98  return true;
99  }// end ProcessHits
100 } // end namespace larg4
larg4::SimEnergyDepositSD::hitCollection
sim::SimEnergyDepositCollection hitCollection
Definition: SimEnergyDepositSD.h:27
larg4::SimEnergyDepositSD::Initialize
void Initialize(G4HCofThisEvent *)
Definition: SimEnergyDepositSD.cc:45
MeV
util::quantities::megaelectronvolt MeV
Definition: TriggerEfficiencyPlots_module.cc:89
icarus::ns::util::end
auto end(FixedBins< T, C > const &) noexcept
Definition: FixedBins.h:585
sim::SimEnergyDeposit
Energy deposition in the active material.
Definition: SimEnergyDeposit.h:42
larg4::SimEnergyDepositSD::ProcessHits
G4bool ProcessHits(G4Step *, G4TouchableHistory *)
Definition: SimEnergyDepositSD.cc:50
name
then echo fcl name
Definition: sbndpoms_genfclwithrunnumber_maker.sh:220
geo::Point_t
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:184
Generated by   doxygen 1.8.5