Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer << 3 // * DISCLAIMER * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th << 5 // * The following disclaimer summarizes all the specific disclaimers * 6 // * the Geant4 Collaboration. It is provided << 6 // * of contributors to this software. The specific disclaimers,which * 7 // * conditions of the Geant4 Software License << 7 // * govern, are listed with their locations in: * 8 // * LICENSE and available at http://cern.ch/ << 8 // * http://cern.ch/geant4/license * 9 // * include a list of copyright holders. << 10 // * 9 // * * 11 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file << 14 // * use. * 16 // * for the full disclaimer and the limitatio << 17 // * 15 // * * 18 // * This code implementation is the result << 16 // * This code implementation is the intellectual property of the * 19 // * technical work of the GEANT4 collaboratio << 17 // * GEANT4 collaboration. * 20 // * By using, copying, modifying or distri << 18 // * By copying, distributing or modifying the Program (or any work * 21 // * any work based on the software) you ag << 19 // * based on the Program) you indicate your acceptance of this * 22 // * use in resulting scientific publicati << 20 // * statement, and all its terms. * 23 // * acceptance of all terms of the Geant4 Sof << 24 // ******************************************* 21 // ******************************************************************** 25 // 22 // 26 // 23 // >> 24 // $Id: G4Decay.hh,v 1.8 2002/02/12 02:14:03 kurasige Exp $ >> 25 // GEANT4 tag $Name: geant4-04-01 $ 27 // 26 // 28 // 27 // 29 // ------------------------------------------- 28 // ------------------------------------------------------------ 30 // GEANT 4 class header file 29 // GEANT 4 class header file 31 // 30 // 32 // History: first implementation, based o 31 // History: first implementation, based on object model of 33 // 7 July 1996 H.Kurashige 32 // 7 July 1996 H.Kurashige 34 // ------------------------------------------- 33 // ------------------------------------------------------------ 35 // New Physics scheme 18 Jan. 1997 34 // New Physics scheme 18 Jan. 1997 H.Kurahige 36 // ------------------------------------------- 35 // ------------------------------------------------------------ 37 // modified 4 Feb. 1997 36 // modified 4 Feb. 1997 H.Kurahige 38 // modified 8 Sep. 1997 37 // modified 8 Sep. 1997 H.Kurahige 39 // remove BuildPhysicsTable() 27 Nov. 1997 38 // remove BuildPhysicsTable() 27 Nov. 1997 H.Kurashige 40 // modified for new ParticleChange 12 Mar. 1 39 // modified for new ParticleChange 12 Mar. 1998 H.Kurashige 41 // added aPhysicsTable 2 Aug. 1998 40 // added aPhysicsTable 2 Aug. 1998 H.Kurashige 42 // PreAssignedDecayTime 18 Jan. 2001 41 // PreAssignedDecayTime 18 Jan. 2001 H.Kurashige 43 // Add External Decayer 23 Feb. 2001 42 // Add External Decayer 23 Feb. 2001 H.Kurashige 44 // Remove PhysicsTable 12 Feb. 2002 43 // Remove PhysicsTable 12 Feb. 2002 H.Kurashige 45 // Fixed bug in PostStepGPIL << 44 // 46 // in case of stopping during AlongStepDoIt << 47 // Add GetRemainderLifeTime 10 Aug/2004 H.K << 48 // Add DaughterPolarization 23 July 2008 << 49 << 50 << 51 #ifndef G4Decay_h 45 #ifndef G4Decay_h 52 #define G4Decay_h 1 46 #define G4Decay_h 1 53 47 54 #include "G4ios.hh" 48 #include "G4ios.hh" 55 #include "globals.hh" 49 #include "globals.hh" 56 #include "G4VRestDiscreteProcess.hh" 50 #include "G4VRestDiscreteProcess.hh" 57 #include "G4ParticleChangeForDecay.hh" 51 #include "G4ParticleChangeForDecay.hh" 58 #include "G4DecayProcessType.hh" << 59 << 60 class G4VExtDecayer; 52 class G4VExtDecayer; 61 53 62 class G4Decay : public G4VRestDiscreteProcess 54 class G4Decay : public G4VRestDiscreteProcess 63 { 55 { 64 // Class Description 56 // Class Description 65 // This class is a decay process 57 // This class is a decay process 66 58 67 public: 59 public: 68 // Constructors 60 // Constructors 69 G4Decay(const G4String& processName ="Deca 61 G4Decay(const G4String& processName ="Decay"); 70 62 71 // Destructor 63 // Destructor 72 virtual ~G4Decay(); 64 virtual ~G4Decay(); 73 65 74 private: 66 private: 75 // copy constructor 67 // copy constructor 76 G4Decay(const G4Decay &right); 68 G4Decay(const G4Decay &right); 77 69 78 // Assignment Operation (generated) 70 // Assignment Operation (generated) 79 G4Decay & operator=(const G4Decay &right 71 G4Decay & operator=(const G4Decay &right); 80 72 81 public: //With Description 73 public: //With Description 82 // G4Decay Process has both 74 // G4Decay Process has both 83 // PostStepDoIt (for decay in flight) 75 // PostStepDoIt (for decay in flight) 84 // and 76 // and 85 // AtRestDoIt (for decay at rest) 77 // AtRestDoIt (for decay at rest) 86 78 87 virtual G4VParticleChange *PostStepDoIt( 79 virtual G4VParticleChange *PostStepDoIt( 88 const G4Track& aTrack, 80 const G4Track& aTrack, 89 const G4Step& aSt 81 const G4Step& aStep 90 ) override; << 82 ); 91 83 92 virtual G4VParticleChange* AtRestDoIt( 84 virtual G4VParticleChange* AtRestDoIt( 93 const G4Track& aTrack, 85 const G4Track& aTrack, 94 const G4Step& aStep 86 const G4Step& aStep 95 ) override; << 87 ); 96 88 97 virtual void BuildPhysicsTable(const G4Pa << 89 virtual void BuildPhysicsTable(const G4ParticleDefinition&); 98 // In G4Decay, thePhysicsTable stores val 90 // In G4Decay, thePhysicsTable stores values of 99 // beta * std::sqrt( 1 - beta*beta) << 91 // beta * sqrt( 1 - beta*beta) 100 // as a function of normalized kinetic en 92 // as a function of normalized kinetic enregy (=Ekin/mass), 101 // becasuse this table is universal for a 93 // becasuse this table is universal for all particle types, 102 94 103 95 104 virtual G4bool IsApplicable(const G4Partic << 96 virtual G4bool IsApplicable(const G4ParticleDefinition&); 105 // returns "true" if the decay process can 97 // returns "true" if the decay process can be applied to 106 // the particle type. 98 // the particle type. 107 99 108 protected: // With Description 100 protected: // With Description 109 virtual G4VParticleChange* DecayIt( 101 virtual G4VParticleChange* DecayIt( 110 const G4Track& aTrack, 102 const G4Track& aTrack, 111 const G4Step& aStep 103 const G4Step& aStep 112 ); 104 ); 113 // The DecayIt() method returns by pointer 105 // The DecayIt() method returns by pointer a particle-change object, 114 // which has information of daughter parti 106 // which has information of daughter particles. 115 107 116 // Set daughter polarization << 108 public: 117 // NO OPERATION in the base class of G4De << 118 virtual void DaughterPolarization(const G4 << 119 G4DecayProducts* products); << 120 << 121 public: << 122 virtual G4double AtRestGetPhysicalInteract 109 virtual G4double AtRestGetPhysicalInteractionLength( 123 const G4Track& tr 110 const G4Track& track, 124 G4ForceCondition* 111 G4ForceCondition* condition 125 ) override; << 112 ); 126 113 127 virtual G4double PostStepGetPhysicalIntera 114 virtual G4double PostStepGetPhysicalInteractionLength( 128 const G4Track& tr 115 const G4Track& track, 129 G4double previo 116 G4double previousStepSize, 130 G4ForceCondition* 117 G4ForceCondition* condition 131 ) override; << 118 ); >> 119 132 120 133 protected: // With Description 121 protected: // With Description 134 // GetMeanFreePath returns ctau*beta*gamma 122 // GetMeanFreePath returns ctau*beta*gamma for decay in flight 135 // GetMeanLifeTime returns ctau for decay 123 // GetMeanLifeTime returns ctau for decay at rest 136 virtual G4double GetMeanFreePath(const G4T 124 virtual G4double GetMeanFreePath(const G4Track& aTrack, 137 G4double previ 125 G4double previousStepSize, 138 G4ForceCondition 126 G4ForceCondition* condition 139 ) override; << 127 ); 140 128 141 virtual G4double GetMeanLifeTime(const G4T 129 virtual G4double GetMeanLifeTime(const G4Track& aTrack, 142 G4ForceCondition 130 G4ForceCondition* condition 143 ) override; << 131 ); 144 << 145 public: //With Description << 146 virtual void StartTracking(G4Track*) over << 147 virtual void EndTracking() override; << 148 // inform Start/End of tracking for each << 149 132 150 public: //With Description << 133 public: //With Description 151 void SetExtDecayer(G4VExtDecayer*); 134 void SetExtDecayer(G4VExtDecayer*); 152 const G4VExtDecayer* GetExtDecayer() cons 135 const G4VExtDecayer* GetExtDecayer() const; 153 // Set/Get External Decayer 136 // Set/Get External Decayer 154 << 137 155 G4double GetRemainderLifeTime() const; << 138 public: 156 //Get Remainder of life time at rest decay << 139 void SetVerboseLevel(G4int value); 157 << 140 G4int GetVerboseLevel() const; 158 virtual void ProcessDescription(std::ostre << 159 // << 160 141 161 protected: << 142 private: 162 G4int verboseLevel; 143 G4int verboseLevel; 163 // controle flag for output message 144 // controle flag for output message 164 // 0: Silent 145 // 0: Silent 165 // 1: Warning message 146 // 1: Warning message 166 // 2: More 147 // 2: More 167 148 168 protected: << 149 private: 169 // HighestValue. 150 // HighestValue. 170 const G4double HighestValue; 151 const G4double HighestValue; 171 152 172 // Remainder of life time at rest 153 // Remainder of life time at rest 173 G4double fRemainderLifeTim 154 G4double fRemainderLifeTime; 174 << 155 175 // ParticleChange for decay process 156 // ParticleChange for decay process 176 G4ParticleChangeForDecay fParticleChangeFo 157 G4ParticleChangeForDecay fParticleChangeForDecay; 177 158 178 // External Decayer 159 // External Decayer 179 G4VExtDecayer* pExtDecayer; 160 G4VExtDecayer* pExtDecayer; 180 }; 161 }; 181 162 >> 163 inline G4double G4Decay::PostStepGetPhysicalInteractionLength( >> 164 const G4Track& track, >> 165 G4double previousStepSize, >> 166 G4ForceCondition* condition >> 167 ) >> 168 { >> 169 // pre-assigned Decay time >> 170 G4double pTime = track.GetDynamicParticle()->GetPreAssignedDecayProperTime(); >> 171 >> 172 if (pTime < 0.) { >> 173 // normal case >> 174 return G4VRestDiscreteProcess::PostStepGetPhysicalInteractionLength(track, previousStepSize, condition); >> 175 } >> 176 >> 177 // condition is set to "Not Forced" >> 178 *condition = NotForced; >> 179 >> 180 // reminder proper time >> 181 fRemainderLifeTime = pTime - track.GetProperTime(); >> 182 if (fRemainderLifeTime <= 0.0) fRemainderLifeTime = DBL_MIN; >> 183 >> 184 // use pre-assigned Decay time to determine PIL >> 185 G4double tau = track.GetDefinition()->GetPDGLifeTime(); >> 186 return (fRemainderLifeTime/tau)*GetMeanFreePath(track, previousStepSize, condition); >> 187 >> 188 } >> 189 inline >> 190 G4double G4Decay::AtRestGetPhysicalInteractionLength( >> 191 const G4Track& track, >> 192 G4ForceCondition* condition >> 193 ) >> 194 { >> 195 G4double pTime = track.GetDynamicParticle()->GetPreAssignedDecayProperTime(); >> 196 if (pTime >= 0.) { >> 197 fRemainderLifeTime = pTime - track.GetProperTime(); >> 198 if (fRemainderLifeTime <= 0.0) fRemainderLifeTime = DBL_MIN; >> 199 } else { >> 200 fRemainderLifeTime = >> 201 G4VRestDiscreteProcess::AtRestGetPhysicalInteractionLength(track, condition ); >> 202 } >> 203 return fRemainderLifeTime; >> 204 } >> 205 >> 206 inline >> 207 void G4Decay::SetVerboseLevel(G4int value){ verboseLevel = value; } >> 208 >> 209 inline >> 210 G4int G4Decay::GetVerboseLevel() const { return verboseLevel; } >> 211 182 inline 212 inline 183 G4VParticleChange* G4Decay::AtRestDoIt( 213 G4VParticleChange* G4Decay::AtRestDoIt( 184 const G4Track& aTrack, 214 const G4Track& aTrack, 185 const G4Step& aStep 215 const G4Step& aStep 186 ) 216 ) 187 { 217 { 188 return DecayIt(aTrack, aStep); 218 return DecayIt(aTrack, aStep); 189 } 219 } 190 220 >> 221 inline >> 222 G4VParticleChange* G4Decay::PostStepDoIt( >> 223 const G4Track& aTrack, >> 224 const G4Step& aStep >> 225 ) >> 226 { >> 227 return DecayIt(aTrack, aStep); >> 228 } 191 229 192 inline 230 inline 193 const G4VExtDecayer* G4Decay::GetExtDecayer() << 231 void G4Decay::SetExtDecayer(G4VExtDecayer* val) 194 { 232 { 195 return pExtDecayer; << 233 pExtDecayer = val; 196 } 234 } 197 235 198 inline 236 inline 199 G4double G4Decay::GetRemainderLifeTime() cons << 237 const G4VExtDecayer* G4Decay::GetExtDecayer() const 200 { 238 { 201 return fRemainderLifeTime; << 239 return pExtDecayer; 202 } 240 } 203 241 204 #endif 242 #endif 205 243 206 244 207 245 208 246 209 247 210 248 211 249 212 250 213 251 214 252 215 253