Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/decay/include/G4UnknownDecay.hh

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /processes/decay/include/G4UnknownDecay.hh (Version 11.3.0) and /processes/decay/include/G4UnknownDecay.hh (Version 9.5.p1)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                         3 // * License and Disclaimer                                           *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th      5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided      6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License      7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/      8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.           9 // * include a list of copyright holders.                             *
 10 // *                                               10 // *                                                                  *
 11 // * Neither the authors of this software syst     11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file      15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitatio     16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                               17 // *                                                                  *
 18 // * This  code  implementation is the result      18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboratio     19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distri     20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  ag     21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publicati     22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Sof     23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // *******************************************     24 // ********************************************************************
 25 //                                                 25 //
 26 //                                                 26 //
                                                   >>  27 // $Id: G4UnknownDecay.hh,v 1.3 2006-06-29 19:30:56 gunter Exp $
                                                   >>  28 // GEANT4 tag $Name: not supported by cvs2svn $
 27 //                                                 29 //
 28 //                                                 30 //
 29 // -------------------------------------------     31 // ------------------------------------------------------------
 30 //      GEANT 4 class header file                  32 //      GEANT 4 class header file
 31 //                                                 33 //
 32 //                                                 34 //
 33 #ifndef G4UnknownDecay_h                           35 #ifndef G4UnknownDecay_h
 34 #define G4UnknownDecay_h 1                         36 #define G4UnknownDecay_h 1
 35                                                    37 
 36 #include <CLHEP/Units/PhysicalConstants.h>     << 
 37                                                << 
 38 #include "G4ios.hh"                                38 #include "G4ios.hh"
 39 #include "globals.hh"                              39 #include "globals.hh"
 40 #include "G4VDiscreteProcess.hh"                   40 #include "G4VDiscreteProcess.hh"
 41 #include "G4ParticleChangeForDecay.hh"             41 #include "G4ParticleChangeForDecay.hh"
 42                                                    42 
 43 class G4UnknownDecay : public G4VDiscreteProce     43 class G4UnknownDecay : public G4VDiscreteProcess 
 44 {                                                  44 {
 45  // Class Description                              45  // Class Description
 46   //  This class is a decay process for "unkno     46   //  This class is a decay process for "unknown" particle.
 47                                                    47 
 48   public:                                          48   public:
 49     //  Constructors                               49     //  Constructors 
 50     G4UnknownDecay(const G4String& processName     50     G4UnknownDecay(const G4String& processName ="UnknownDecay");
 51                                                    51 
 52     //  Destructor                                 52     //  Destructor
 53     virtual ~G4UnknownDecay();                     53     virtual ~G4UnknownDecay();
 54                                                    54 
 55   private:                                         55   private:
 56     //  copy constructor                           56     //  copy constructor
 57       G4UnknownDecay(const G4UnknownDecay &rig     57       G4UnknownDecay(const G4UnknownDecay &right);
 58                                                    58 
 59     //  Assignment Operation (generated)           59     //  Assignment Operation (generated)
 60       G4UnknownDecay & operator=(const G4Unkno     60       G4UnknownDecay & operator=(const G4UnknownDecay &right);
 61                                                    61 
 62   public: //With Description                       62   public: //With Description
 63                                                    63   
 64      virtual G4VParticleChange *PostStepDoIt(      64      virtual G4VParticleChange *PostStepDoIt(
 65            const G4Track& aTrack,                  65            const G4Track& aTrack,
 66                              const G4Step& aSt     66                              const G4Step& aStep
 67                             ) override;        <<  67                             );
 68                                                    68 
 69      virtual void BuildPhysicsTable(const G4Pa <<  69      virtual void BuildPhysicsTable(const G4ParticleDefinition&); 
 70      // In G4UnknownDecay, thePhysicsTable sto     70      // In G4UnknownDecay, thePhysicsTable stores values of
 71     //    beta * std::sqrt( 1 - beta*beta)         71     //    beta * std::sqrt( 1 - beta*beta) 
 72     //  as a function of normalized kinetic en     72     //  as a function of normalized kinetic enregy (=Ekin/mass),
 73     //  becasuse this table is universal for a     73     //  becasuse this table is universal for all particle types,
 74                                                    74 
 75                                                    75 
 76     virtual G4bool IsApplicable(const G4Partic <<  76     virtual G4bool IsApplicable(const G4ParticleDefinition&);
 77     // returns "true" if the decay process can     77     // returns "true" if the decay process can be applied to
 78     // the particle type.                          78     // the particle type. 
 79                                                    79  
 80     virtual void ProcessDescription(std::ostre << 
 81     //                                         << 
 82                                                << 
 83   protected: // With Description                   80   protected: // With Description
 84     virtual G4VParticleChange* DecayIt(            81     virtual G4VParticleChange* DecayIt(
 85            const G4Track& aTrack,                  82            const G4Track& aTrack,
 86            const G4Step&  aStep                    83            const G4Step&  aStep
 87           ) ;                                  <<  84           );
 88     // The DecayIt() method returns by pointer     85     // The DecayIt() method returns by pointer a particle-change object,
 89     // which has information of daughter parti     86     // which has information of daughter particles.
 90                                                    87 
 91   public:                                          88   public:
 92                                                    89 
 93     virtual G4double PostStepGetPhysicalIntera     90     virtual G4double PostStepGetPhysicalInteractionLength(
 94                              const G4Track& tr     91                              const G4Track& track,
 95                              G4double   previo     92                              G4double   previousStepSize,
 96                              G4ForceCondition*     93                              G4ForceCondition* condition
 97                             ) override;        <<  94                             );
 98                                                    95 
 99                                                    96 
100   protected: // With Description                   97   protected: // With Description
101     // GetMeanFreePath returns ctau*beta*gamma     98     // GetMeanFreePath returns ctau*beta*gamma for decay in flight 
102     virtual G4double GetMeanFreePath(const G4T     99     virtual G4double GetMeanFreePath(const G4Track& aTrack,
103                               G4double   previ    100                               G4double   previousStepSize,
104                               G4ForceCondition    101                               G4ForceCondition* condition
105                              ) override;       << 102                              );
                                                   >> 103 
                                                   >> 104   public:
                                                   >> 105      void  SetVerboseLevel(G4int value);
                                                   >> 106      G4int GetVerboseLevel() const;
106                                                   107 
107   private:                                        108   private:
108      G4int verboseLevel;                          109      G4int verboseLevel;
109      // controle flag for output message          110      // controle flag for output message
110      //  0: Silent                                111      //  0: Silent
111      //  1: Warning message                       112      //  1: Warning message
112      //  2: More                                  113      //  2: More
113                                                   114 
114   private:                                        115   private:
115     // HighestValue.                              116     // HighestValue.
116     const G4double HighestValue;                  117     const G4double HighestValue;
117                                                   118  
118     // ParticleChange for decay process           119     // ParticleChange for decay process
119     G4ParticleChangeForDecay fParticleChangeFo    120     G4ParticleChangeForDecay fParticleChangeForDecay;
120                                                   121     
121 };                                                122 };
122                                                   123 
123 inline                                            124 inline 
124  G4double G4UnknownDecay::PostStepGetPhysicalI    125  G4double G4UnknownDecay::PostStepGetPhysicalInteractionLength(
125                              const G4Track& tr    126                              const G4Track& track,
126                              G4double   /*prev    127                              G4double   /*previousStepSize*/,
127                              G4ForceCondition*    128                              G4ForceCondition* condition
128                             )                     129                             )
129 {                                                 130 {
130   // pre-assigned UnknownDecay time               131   // pre-assigned UnknownDecay time
131   G4double pTime = track.GetDynamicParticle()-    132   G4double pTime = track.GetDynamicParticle()->GetPreAssignedDecayProperTime();
132                                                   133 
133   if (pTime < 0.) pTime = DBL_MIN;                134   if (pTime < 0.) pTime = DBL_MIN;
134                                                   135 
135   // condition is set to "Not Forced"             136   // condition is set to "Not Forced"
136   *condition = NotForced;                         137   *condition = NotForced;
137                                                   138   
138   // reminder proper time                         139   // reminder proper time
139   G4double remainder = pTime - track.GetProper    140   G4double remainder = pTime - track.GetProperTime();
140   if (remainder <= 0.0) remainder = DBL_MIN;      141   if (remainder <= 0.0) remainder = DBL_MIN;
141                                                   142   
142   // use pre-assigned Decay time to determine     143   // use pre-assigned Decay time to determine PIL
143   //return GetMeanFreePath(track, previousStep    144   //return GetMeanFreePath(track, previousStepSize, condition);
144   return remainder*CLHEP::c_light;             << 145   return remainder*c_light;
145                                                   146 
146 }                                                 147 }
                                                   >> 148 
                                                   >> 149 inline
                                                   >> 150  void  G4UnknownDecay::SetVerboseLevel(G4int value){ verboseLevel = value; }
                                                   >> 151 
                                                   >> 152 inline
                                                   >> 153  G4int G4UnknownDecay::GetVerboseLevel() const { return verboseLevel; }
147                                                   154 
148 inline                                            155 inline  
149   G4VParticleChange* G4UnknownDecay::PostStepD    156   G4VParticleChange* G4UnknownDecay::PostStepDoIt(
150            const G4Track& aTrack,                 157            const G4Track& aTrack,
151            const G4Step&  aStep                   158            const G4Step&  aStep
152           )                                       159           )
153 {                                                 160 {
154   return DecayIt(aTrack, aStep);                  161   return DecayIt(aTrack, aStep);
155 }                                                 162 }
156                                                   163 
157 #endif                                            164 #endif
158                                                   165 
159                                                   166 
160                                                   167 
161                                                   168 
162                                                   169 
163                                                   170 
164                                                   171 
165                                                   172 
166                                                   173 
167                                                   174 
168                                                   175