| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 #ifdef USE_INFERENCE
27 # include "Par04MLFastSimModel.hh"
28
29 # include "Par04InferenceSetup.hh" // for Pa
30
31 # include "G4Electron.hh" // for G4Electron
32 # include "G4FastHit.hh" // for G4FastHit
33 # include "G4FastSimHitMaker.hh" // for G4Fa
34 # include "G4Gamma.hh" // for G4Gamma
35 # include "G4Positron.hh" // for G4Positron
36
37 # include <G4FastStep.hh> // for G4FastStep
38 # include <G4FastTrack.hh> // for G4FastTrac
39 # include <G4Track.hh> // for G4Track
40 # include <G4VFastSimulationModel.hh> // for
41 # include <stddef.h> // for size_t
42 class G4ParticleDefinition;
43 class G4Region;
44
45 //....oooOO0OOooo........oooOO0OOooo........oo
46
47 Par04MLFastSimModel::Par04MLFastSimModel(G4Str
48 : G4VFastSimulationModel(aModelName, aEnvelo
49 fInference(new Par04InferenceSetup),
50 fHitMaker(new G4FastSimHitMaker),
51 fParallelHitMaker(new G4FastSimHitMaker)
52 {
53 fParallelHitMaker->SetNameOfWorldWithSD("par
54 }
55
56 //....oooOO0OOooo........oooOO0OOooo........oo
57
58 Par04MLFastSimModel::Par04MLFastSimModel(G4Str
59 : G4VFastSimulationModel(aModelName),
60 fInference(new Par04InferenceSetup),
61 fHitMaker(new G4FastSimHitMaker),
62 fParallelHitMaker(new G4FastSimHitMaker)
63 {
64 fParallelHitMaker->SetNameOfWorldWithSD("par
65 }
66
67 //....oooOO0OOooo........oooOO0OOooo........oo
68
69 Par04MLFastSimModel::~Par04MLFastSimModel() {}
70
71 //....oooOO0OOooo........oooOO0OOooo........oo
72
73 G4bool Par04MLFastSimModel::IsApplicable(const
74 {
75 return &aParticleType == G4Electron::Electro
76 || &aParticleType == G4Positron::Posi
77 || &aParticleType == G4Gamma::GammaDe
78 }
79
80 //....oooOO0OOooo........oooOO0OOooo........oo
81
82 G4bool Par04MLFastSimModel::ModelTrigger(const
83 {
84 return fInference->IfTrigger(aFastTrack.GetP
85 }
86
87 //....oooOO0OOooo........oooOO0OOooo........oo
88
89 void Par04MLFastSimModel::DoIt(const G4FastTra
90 {
91 // remove particle from further processing b
92 aFastStep.KillPrimaryTrack();
93 aFastStep.SetPrimaryTrackPathLength(0.0);
94 G4double energy = aFastTrack.GetPrimaryTrack
95 aFastStep.SetTotalEnergyDeposited(energy);
96 G4ThreeVector position = aFastTrack.GetPrima
97 G4ThreeVector direction = aFastTrack.GetPrim
98
99 // calculate the incident angle
100 G4float angle = direction.theta();
101
102 // calculate how to deposit energy within th
103 // get it from inference model
104 fInference->GetEnergies(fEnergies, energy, a
105 fInference->GetPositions(fPositions, positio
106
107 // deposit energy in the detector using calc
108 // and positions
109 for (size_t iHit = 0; iHit < fPositions.size
110 if (fEnergies[iHit] > 0.0005) {
111 // Place hit in the physical readout of
112 fParallelHitMaker->make(G4FastHit(fPosit
113 // Place hit in the virtual grid, for va
114 fHitMaker->make(G4FastHit(fPositions[iHi
115 }
116 }
117 }
118 #endif
119