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 // 27 /// \file Par02FastSimModelTracker.cc 28 /// \brief Implementation of the Par02FastSimM 29 30 #include "Par02FastSimModelTracker.hh" 31 32 #include "Par02EventInformation.hh" 33 #include "Par02Output.hh" 34 #include "Par02PrimaryParticleInformation.hh" 35 #include "Par02Smearer.hh" 36 37 #include "G4AnalysisManager.hh" 38 #include "G4Electron.hh" 39 #include "G4Event.hh" 40 #include "G4FieldTrack.hh" 41 #include "G4FieldTrackUpdator.hh" 42 #include "G4Gamma.hh" 43 #include "G4PathFinder.hh" 44 #include "G4Positron.hh" 45 #include "G4RunManager.hh" 46 #include "G4SystemOfUnits.hh" 47 #include "G4Track.hh" 48 #include "Randomize.hh" 49 50 //....oooOO0OOooo........oooOO0OOooo........oo 51 52 Par02FastSimModelTracker::Par02FastSimModelTra 53 G4String aModelName, G4Region* aEnvelope, Pa 54 : G4VFastSimulationModel(aModelName, aEnvelo 55 fCalculateParametrisation(), 56 fParametrisation(aType) 57 {} 58 59 //....oooOO0OOooo........oooOO0OOooo........oo 60 61 Par02FastSimModelTracker::Par02FastSimModelTra 62 : G4VFastSimulationModel(aModelName, aEnvelo 63 fCalculateParametrisation(), 64 fParametrisation(Par02DetectorParametrisat 65 {} 66 67 //....oooOO0OOooo........oooOO0OOooo........oo 68 69 Par02FastSimModelTracker::Par02FastSimModelTra 70 : G4VFastSimulationModel(aModelName), 71 fCalculateParametrisation(), 72 fParametrisation(Par02DetectorParametrisat 73 {} 74 75 //....oooOO0OOooo........oooOO0OOooo........oo 76 77 Par02FastSimModelTracker::~Par02FastSimModelTr 78 79 //....oooOO0OOooo........oooOO0OOooo........oo 80 81 G4bool Par02FastSimModelTracker::IsApplicable( 82 { 83 return aParticleType.GetPDGCharge() != 0; / 84 } 85 86 //....oooOO0OOooo........oooOO0OOooo........oo 87 88 G4bool Par02FastSimModelTracker::ModelTrigger( 89 { 90 return true; // No kinematical restrictions 91 } 92 93 //....oooOO0OOooo........oooOO0OOooo........oo 94 95 void Par02FastSimModelTracker::DoIt(const G4Fa 96 { 97 G4cout << " ________Tracker model triggered 98 99 // Calculate the final position (at the oute 100 // of the particle with the momentum at the 101 102 G4Track track = *aFastTrack.GetPrimaryTrack( 103 G4FieldTrack aFieldTrack('0'); 104 G4FieldTrackUpdator::Update(&aFieldTrack, &t 105 106 G4double retSafety = -1.0; 107 ELimited retStepLimited; 108 G4FieldTrack endTrack('a'); 109 G4double currentMinimumStep = 10.0 * m; // 110 // 111 G4PathFinder* fPathFinder = G4PathFinder::Ge 112 /*G4double lengthAlongCurve = */ 113 fPathFinder->ComputeStep(aFieldTrack, curren 114 aFastTrack.GetPrima 115 retStepLimited, end 116 117 // Place the particle at the tracking detect 118 // (at the place it would reach without the 119 aFastStep.ProposePrimaryTrackFinalPosition(e 120 121 // Consider only primary tracks (do nothing 122 G4ThreeVector Porg = aFastTrack.GetPrimaryTr 123 if (!aFastTrack.GetPrimaryTrack()->GetParent 124 auto info = (Par02EventInformation*)G4Even 125 if (info->GetDoSmearing()) { 126 // Smearing according to the tracking de 127 G4double res = fCalculateParametrisation 128 Par02DetectorParametrisation::eTRACKER 129 G4double eff = fCalculateParametrisation 130 Par02DetectorParametrisation::eTRACKER 131 G4ThreeVector Psm; 132 Psm = Par02Smearer::Instance()->SmearMom 133 Par02Output::Instance()->FillHistogram(0 134 // Setting the values of Psm, res and ef 135 ((Par02PrimaryParticleInformation*)(cons 136 aF 137 138 139 -> 140 ->SetTrackerMomentum(Psm); 141 ((Par02PrimaryParticleInformation*)(cons 142 aF 143 144 145 -> 146 ->SetTrackerResolution(res); 147 ((Par02PrimaryParticleInformation*)(cons 148 aF 149 150 151 -> 152 ->SetTrackerEfficiency(eff); 153 } 154 else { 155 // No smearing: simply setting the value 156 ((Par02PrimaryParticleInformation*)(cons 157 aF 158 159 160 -> 161 ->SetTrackerMomentum(Porg); 162 } 163 } 164 } 165 166 //....oooOO0OOooo........oooOO0OOooo........oo 167