Geant4 Cross Reference |
1 // 2 // ******************************************************************** 3 // * License and Disclaimer * 4 // * * 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. * 10 // * * 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitation of liability. * 17 // * * 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************************************** 25 // 26 27 // GEANT 4 class implementation file --- Copyright CERN 1995 28 29 // History: 30 // 29.02.04 V.Ivanchenko create 31 // 28.07.05, P.Gumplinger add G4ProcessType to constructor 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(const G4String& processName, 45 G4ProcessType type) 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::~G4VTransitionRadiation() 61 { 62 Clear(); 63 theManager->DeRegister(this); 64 } 65 66 void G4VTransitionRadiation::ProcessDescription(std::ostream& out) const 67 { 68 out << "Generic process of transition radiation.\n"; 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::PostStepDoIt(const G4Track& track, 85 const G4Step& step) 86 { 87 // Fill temporary vectors 88 const G4Material* material = track.GetMaterial(); 89 G4double length = step.GetStepLength(); 90 G4ThreeVector direction = track.GetMomentumDirection(); 91 92 if(nSteps == 0) 93 { 94 nSteps = 1; 95 materials.push_back(material); 96 steps.push_back(length); 97 const G4StepPoint* point = step.GetPreStepPoint(); 98 startingPosition = point->GetPosition(); 99 startingDirection = point->GetMomentumDirection(); 100 G4bool valid = true; 101 G4ThreeVector n = G4TransportationManager::GetTransportationManager() 102 ->GetNavigatorForTracking() 103 ->GetLocalExitNormal(&valid); 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 = G4TransportationManager::GetTransportationManager() 122 ->GetNavigatorForTracking() 123 ->GetLocalExitNormal(&valid); 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()->GetRegion() != region || 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(*pParticleChange, materials, steps, normals, 142 startingPosition, track); 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 G4Region* reg) { region = reg; } 159 160 /////////////////////////////////////////////////////////////////////// 161 void G4VTransitionRadiation::SetModel(G4VTRModel* mod) { model = mod; } 162 163 /////////////////////////////////////////////////////////////////////// 164 G4double G4VTransitionRadiation::GetMeanFreePath(const G4Track& track, G4double, 165 G4ForceCondition* condition) 166 { 167 if(nSteps > 0) 168 { 169 *condition = StronglyForced; 170 } 171 else 172 { 173 *condition = NotForced; 174 if(track.GetKineticEnergy() / track.GetDefinition()->GetPDGMass() + 1.0 > 175 gammaMin && 176 track.GetVolume()->GetLogicalVolume()->GetRegion() == region) 177 { 178 *condition = StronglyForced; 179 } 180 } 181 return DBL_MAX; // so TR doesn't limit mean free path 182 } 183