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