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 // GEANT 4 class implementation file --- Copyr 28 29 // History: 30 // 29.02.04 V.Ivanchenko create 31 // 28.07.05, P.Gumplinger add G4ProcessType to 32 33 #include "G4VTransitionRadiation.hh" 34 35 #include "G4EmProcessSubType.hh" 36 #include "G4LossTableManager.hh" 37 #include "G4Material.hh" 38 #include "G4ParticleDefinition.hh" 39 #include "G4Region.hh" 40 #include "G4TransportationManager.hh" 41 #include "G4VTRModel.hh" 42 43 ////////////////////////////////////////////// 44 G4VTransitionRadiation::G4VTransitionRadiation 45 46 : G4VDiscreteProcess(processName, type) 47 , region(nullptr) 48 , model(nullptr) 49 , gammaMin(100.) 50 , cosDThetaMax(std::cos(0.1)) 51 , nSteps(0) 52 { 53 SetProcessSubType(fTransitionRadiation); 54 Clear(); 55 theManager = G4LossTableManager::Instance(); 56 theManager->Register(this); 57 } 58 59 ////////////////////////////////////////////// 60 G4VTransitionRadiation::~G4VTransitionRadiatio 61 { 62 Clear(); 63 theManager->DeRegister(this); 64 } 65 66 void G4VTransitionRadiation::ProcessDescriptio 67 { 68 out << "Generic process of transition radiat 69 70 if(model) 71 model->PrintInfo(); 72 } 73 74 ////////////////////////////////////////////// 75 void G4VTransitionRadiation::Clear() 76 { 77 materials.clear(); 78 steps.clear(); 79 normals.clear(); 80 nSteps = 0; 81 } 82 83 ////////////////////////////////////////////// 84 G4VParticleChange* G4VTransitionRadiation::Pos 85 86 { 87 // Fill temporary vectors 88 const G4Material* material = track.GetMateri 89 G4double length = step.GetStepLen 90 G4ThreeVector direction = track.GetMoment 91 92 if(nSteps == 0) 93 { 94 nSteps = 1; 95 materials.push_back(material); 96 steps.push_back(length); 97 const G4StepPoint* point = step.GetPreStep 98 startingPosition = point->GetPosit 99 startingDirection = point->GetMomen 100 G4bool valid = true; 101 G4ThreeVector n = G4TransportationManager: 102 ->GetNavigatorForTrack 103 ->GetLocalExitNormal(& 104 if(valid) 105 normals.push_back(n); 106 else 107 normals.push_back(direction); 108 } 109 else 110 { 111 if(material == materials[nSteps - 1]) 112 { 113 steps[nSteps - 1] += length; 114 } 115 else 116 { 117 ++nSteps; 118 materials.push_back(material); 119 steps.push_back(length); 120 G4bool valid = true; 121 G4ThreeVector n = G4TransportationManage 122 ->GetNavigatorForTra 123 ->GetLocalExitNormal 124 if(valid) 125 normals.push_back(n); 126 else 127 normals.push_back(direction); 128 } 129 } 130 131 // Check PostStepPoint condition 132 if(track.GetTrackStatus() == fStopAndKill || 133 track.GetVolume()->GetLogicalVolume()->Ge 134 startingDirection.x() * direction.x() + 135 startingDirection.y() * direction.y() 136 startingDirection.z() * direction.z() 137 cosDThetaMax) 138 { 139 if(model) 140 { 141 model->GenerateSecondaries(*pParticleCha 142 startingPosit 143 } 144 Clear(); 145 } 146 147 return pParticleChange; 148 } 149 150 ////////////////////////////////////////////// 151 G4bool G4VTransitionRadiation::IsApplicable( 152 const G4ParticleDefinition& aParticle) 153 { 154 return (aParticle.GetPDGCharge() != 0.0); 155 } 156 157 ////////////////////////////////////////////// 158 void G4VTransitionRadiation::SetRegion(const G 159 160 ////////////////////////////////////////////// 161 void G4VTransitionRadiation::SetModel(G4VTRMod 162 163 ////////////////////////////////////////////// 164 G4double G4VTransitionRadiation::GetMeanFreePa 165 166 { 167 if(nSteps > 0) 168 { 169 *condition = StronglyForced; 170 } 171 else 172 { 173 *condition = NotForced; 174 if(track.GetKineticEnergy() / track.GetDef 175 gammaMin && 176 track.GetVolume()->GetLogicalVolume()-> 177 { 178 *condition = StronglyForced; 179 } 180 } 181 return DBL_MAX; // so TR doesn't limit mean 182 } 183