Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/brachytherapy/src/BrachyPhysicsList.cc

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 /examples/advanced/brachytherapy/src/BrachyPhysicsList.cc (Version 11.3.0) and /examples/advanced/brachytherapy/src/BrachyPhysicsList.cc (Version 9.4)


  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 Author: Susanna Guatelli                       <<  27 // Code developed by:
 28 */                                             <<  28 // S. Agostinelli, F. Foppiano, S. Garelli , M. Tropeano, S.Guatelli
                                                   >>  29 //
                                                   >>  30 // Code review: MGP, 5 November 2006 (still to be completed)
 29 //                                                 31 //
 30 //    **********************************           32 //    **********************************
 31 //    *                                *           33 //    *                                *
 32 //    *     BrachyPhysicsList.cc       *           34 //    *     BrachyPhysicsList.cc       *
 33 //    *                                *           35 //    *                                *
 34 //    **********************************           36 //    **********************************
 35 //                                                 37 //
                                                   >>  38 // $Id: BrachyPhysicsList.cc,v 1.15 2010/11/24 15:56:19 cirrone Exp $
                                                   >>  39 // GEANT4 tag $Name: geant4-09-04 $
                                                   >>  40 //
 36 #include "BrachyPhysicsList.hh"                    41 #include "BrachyPhysicsList.hh"
 37 #include "BrachyPhysicsListMessenger.hh"       <<  42 
 38 #include "G4EmStandardPhysics_option4.hh"      << 
 39 #include "G4EmLivermorePhysics.hh"             << 
 40 #include "G4EmStandardPhysics.hh"              << 
 41 #include "G4EmStandardPhysics_option1.hh"      << 
 42 #include "G4EmStandardPhysics_option2.hh"      << 
 43 #include "G4EmStandardPhysics_option3.hh"      << 
 44 #include "G4EmStandardPhysics_option4.hh"      << 
 45 #include "G4DecayPhysics.hh"                   << 
 46 #include "G4RadioactiveDecayPhysics.hh"        << 
 47 #include "G4EmPenelopePhysics.hh"              << 
 48 #include "G4VPhysicsConstructor.hh"            << 
 49 #include "G4ParticleDefinition.hh"                 43 #include "G4ParticleDefinition.hh"
 50 #include "G4ProductionCutsTable.hh"                44 #include "G4ProductionCutsTable.hh"
 51 #include "G4ProcessManager.hh"                     45 #include "G4ProcessManager.hh"
 52 #include "G4ParticleTypes.hh"                      46 #include "G4ParticleTypes.hh"
 53 #include "G4ios.hh"                            <<  47 #include "G4UnitsTable.hh"
 54 #include "G4StepLimiter.hh"                    <<  48 #include "G4ios.hh"              
 55 #include "G4ParticleDefinition.hh"             <<  49 // gamma
 56 #include "globals.hh"                          <<  50 #include "G4PhotoElectricEffect.hh"
 57 #include "G4SystemOfUnits.hh"                  <<  51 #include "G4LivermorePhotoElectricModel.hh"
 58 #include "G4UAtomicDeexcitation.hh"            <<  52 
 59 #include "G4LossTableManager.hh"               <<  53 #include "G4ComptonScattering.hh"
 60                                                <<  54 #include "G4LivermoreComptonModel.hh"
 61 BrachyPhysicsList::BrachyPhysicsList():  G4VMo <<  55 
 62 {                                              <<  56 #include "G4GammaConversion.hh"
 63 SetVerboseLevel(1);                            <<  57 #include "G4LivermoreGammaConversionModel.hh"
 64                                                <<  58 
 65 G4ProductionCutsTable::GetProductionCutsTable( <<  59 #include "G4RayleighScattering.hh" 
 66 SetDefaultCutValue(0.05 *mm);                  <<  60 #include "G4LivermoreRayleighModel.hh"
 67 DumpCutValuesTable();                          <<  61 
 68                                                <<  62 // e-
 69 // EM physics: default                         <<  63 #include "G4eMultipleScattering.hh"
 70 fEmPhysicsList = new G4EmLivermorePhysics();   <<  64 #include "G4eIonisation.hh"
 71 fEmName="emlivermore";                         <<  65 #include "G4LivermoreIonisationModel.hh"
 72                                                <<  66 #include "G4eBremsstrahlung.hh"
 73 // Add Decay                                   <<  67 #include "G4LivermoreBremsstrahlungModel.hh"
 74 fDecPhysicsList = new G4DecayPhysics();        <<  68 
 75 fRadDecayPhysicsList = new G4RadioactiveDecayP <<  69 // e+
 76 fMessenger = new BrachyPhysicsListMessenger(th <<  70 #include "G4eIonisation.hh" 
                                                   >>  71 #include "G4eBremsstrahlung.hh" 
                                                   >>  72 #include "G4eplusAnnihilation.hh"
                                                   >>  73 
                                                   >>  74 BrachyPhysicsList::BrachyPhysicsList():  G4VUserPhysicsList()
                                                   >>  75 {
                                                   >>  76   SetVerboseLevel(1);
 77 }                                                  77 }
 78                                                    78 
 79 BrachyPhysicsList::~BrachyPhysicsList()            79 BrachyPhysicsList::~BrachyPhysicsList()
 80 {                                              <<  80 {
 81 delete fMessenger;                             << 
 82 delete fDecPhysicsList;                        << 
 83 delete fRadDecayPhysicsList;                   << 
 84 delete fEmPhysicsList;                         << 
 85 }                                                  81 }
 86                                                    82 
 87 void BrachyPhysicsList::ConstructParticle()        83 void BrachyPhysicsList::ConstructParticle()
 88 {                                                  84 {
 89 fDecPhysicsList -> ConstructParticle();        <<  85   // In this method, static member functions should be called
                                                   >>  86   // for all particles which you want to use.
                                                   >>  87   // This ensures that objects of these particle types will be
                                                   >>  88   // created in the program. 
                                                   >>  89 
                                                   >>  90   ConstructBosons();
                                                   >>  91   ConstructLeptons();
                                                   >>  92 
 90 }                                                  93 }
 91                                                    94 
 92 void BrachyPhysicsList::ConstructProcess()     <<  95 void BrachyPhysicsList::ConstructBosons()
                                                   >>  96 { 
                                                   >>  97   // photons
                                                   >>  98   G4Gamma::GammaDefinition();
                                                   >>  99 }
                                                   >> 100 
                                                   >> 101 void BrachyPhysicsList::ConstructLeptons()
 93 {                                                 102 {
 94 AddTransportation();                           << 103   // leptons
 95 fEmPhysicsList -> ConstructProcess();          << 104   G4Electron::ElectronDefinition();
                                                   >> 105   G4Positron::PositronDefinition();
                                                   >> 106 }
 96                                                   107 
 97 // decay physics list                          << 108 void BrachyPhysicsList::ConstructProcess()
 98 fDecPhysicsList -> ConstructProcess();         << 109 {
 99 fRadDecayPhysicsList -> ConstructProcess();    << 110   AddTransportation();
100                                                << 111   ConstructEM();
101 // Deexcitation                                << 
102 // Both Fluorescence and Auger e- emission act << 
103 G4VAtomDeexcitation* de = new G4UAtomicDeexcit << 
104 G4LossTableManager::Instance()->SetAtomDeexcit << 
105 de -> SetFluo(true);                           << 
106 de -> SetAuger(true);                          << 
107                                                << 
108 // To model full Auger cascade include in the  << 
109 // the following UI commands:                  << 
110 // process/em/augerCascade true                << 
111 // process/em/deexcitationIgnoreCut true       << 
112 }                                                 112 }
113                                                   113 
114 void BrachyPhysicsList::AddPhysicsList(const G << 114 void BrachyPhysicsList::ConstructEM()
115 {                                                 115 {
116                                                << 116   theParticleIterator->reset();
117   if (name == fEmName) return;                 << 117 
                                                   >> 118   while( (*theParticleIterator)() ){
118                                                   119 
119   if (name == "emstandard_opt0"){              << 120     G4ParticleDefinition* particle = theParticleIterator->value();
120     fEmName = name;                            << 121     G4ProcessManager* pmanager = particle->GetProcessManager();
121     delete fEmPhysicsList;                     << 122     G4String particleName = particle->GetParticleName();
122     fEmPhysicsList = new G4EmStandardPhysics() << 123     
123     G4cout << "PhysicsList::AddPhysicsList: <" << 124     // Processes
124                                                << 125     
125   } else if (name == "emstandard_opt1"){       << 126     if (particleName == "gamma") {
126     fEmName = name;                            << 127       
127     delete fEmPhysicsList;                     << 128       // Photon     
128     fEmPhysicsList = new G4EmStandardPhysics_o << 129       G4RayleighScattering* theRayleigh = new G4RayleighScattering();
129     G4cout << "PhysicsList::AddPhysicsList: <" << 130       theRayleigh->SetModel(new G4LivermoreRayleighModel());  //not strictly necessary
130   } else if (name == "emstandard_opt2"){       << 131       pmanager->AddDiscreteProcess(theRayleigh);
131     fEmName = name;                            << 132 
132     delete fEmPhysicsList;                     << 133       G4PhotoElectricEffect* thePhotoElectricEffect = new G4PhotoElectricEffect();
133     fEmPhysicsList = new G4EmStandardPhysics_o << 134       thePhotoElectricEffect->SetModel(new G4LivermorePhotoElectricModel());
134     G4cout << "PhysicsList::AddPhysicsList: <" << 135       pmanager->AddDiscreteProcess(thePhotoElectricEffect);
135   } else if (name == "emstandard_opt3"){       << 136   
136     fEmName = name;                            << 137       G4ComptonScattering* theComptonScattering = new G4ComptonScattering();
137     delete fEmPhysicsList;                     << 138       theComptonScattering->SetModel(new G4LivermoreComptonModel());
138     fEmPhysicsList = new G4EmStandardPhysics_o << 139       pmanager->AddDiscreteProcess(theComptonScattering);
139     G4cout << "PhysicsList::AddPhysicsList: <" << 140   
140   } else if (name == "emstandard_opt4"){       << 141       G4GammaConversion* theGammaConversion = new G4GammaConversion();
141     fEmName = name;                            << 142       theGammaConversion->SetModel(new G4LivermoreGammaConversionModel());
142     delete fEmPhysicsList;                     << 143       pmanager->AddDiscreteProcess(theGammaConversion);
143     fEmPhysicsList = new G4EmStandardPhysics_o << 144       
144     G4cout << "PhysicsList::AddPhysicsList: <" << 145     } else if (particleName == "e-") {
145   } else if (name == "empenelope"){            << 146       // Electron
146     fEmName = name;                            << 147       
147     delete fEmPhysicsList;                     << 148       G4eMultipleScattering* msc = new G4eMultipleScattering();
148     fEmPhysicsList = new G4EmPenelopePhysics() << 149       msc->SetStepLimitType(fUseDistanceToBoundary);
149     G4cout << "PhysicsList::AddPhysicsList: <" << 150       pmanager->AddProcess(msc,-1, 1, 1);
150   } else if (name == "emlivermore"){           << 151 
151     fEmName = name;                            << 152       // Ionisation
152     delete fEmPhysicsList;                     << 153       G4eIonisation* eIonisation = new G4eIonisation();
153     fEmPhysicsList = new G4EmLivermorePhysics( << 154       eIonisation->SetEmModel(new G4LivermoreIonisationModel());
154     G4cout << "PhysicsList::AddPhysicsList: <" << 155       eIonisation->SetStepFunction(0.2, 100*um); //improved precision in tracking  
155   } else {                                     << 156       pmanager->AddProcess(eIonisation,-1, 2, 2);
156                                                << 157   
157     G4cout << "PhysicsList::AddPhysicsList: <" << 158       // Bremsstrahlung
158            << " is not defined"                << 159       G4eBremsstrahlung* eBremsstrahlung = new G4eBremsstrahlung();
159            << G4endl;                          << 160       eBremsstrahlung->SetEmModel(new G4LivermoreBremsstrahlungModel());
160   }                                            << 161       pmanager->AddProcess(eBremsstrahlung, -1,-3, 3);
161   G4cout << "PhysicsList::AddPhysicsList: <" < << 162       
162            << " is activated"                  << 163     } else if (particleName == "e+") {
163            << G4endl;                          << 164       // Positron      
                                                   >> 165       G4eMultipleScattering* msc = new G4eMultipleScattering();
                                                   >> 166       msc->SetStepLimitType(fUseDistanceToBoundary);
                                                   >> 167       pmanager->AddProcess(msc,-1, 1, 1);
                                                   >> 168   
                                                   >> 169       // Ionisation
                                                   >> 170       G4eIonisation* eIonisation = new G4eIonisation();
                                                   >> 171       eIonisation->SetStepFunction(0.2, 100*um); //     
                                                   >> 172       pmanager->AddProcess(eIonisation,                 -1, 2, 2);
                                                   >> 173 
                                                   >> 174       //Bremsstrahlung (use default, no low-energy available)
                                                   >> 175       pmanager->AddProcess(new G4eBremsstrahlung(), -1,-1, 3);
                                                   >> 176 
                                                   >> 177       //Annihilation
                                                   >> 178       pmanager->AddProcess(new G4eplusAnnihilation(),0,-1, 4);     
                                                   >> 179       
                                                   >> 180     }
                                                   >> 181   }  
164 }                                                 182 }
165                                                   183 
                                                   >> 184 void BrachyPhysicsList::SetCuts()
                                                   >> 185 {
                                                   >> 186   // The production threshold is fixed to 0.1 mm for all the particles
                                                   >> 187   // Secondary particles with a range bigger than 0.1 mm 
                                                   >> 188   // are generated; otherwise their energy is considered deposited locally
                                                   >> 189 
                                                   >> 190   defaultCutValue = 0.1 * mm;
                                                   >> 191 
                                                   >> 192   const G4double cutForGamma = defaultCutValue;
                                                   >> 193   const G4double cutForElectron = defaultCutValue;
                                                   >> 194   const G4double cutForPositron = defaultCutValue;
                                                   >> 195 
                                                   >> 196   SetCutValue(cutForGamma, "gamma");
                                                   >> 197   SetCutValue(cutForElectron, "e-");
                                                   >> 198   SetCutValue(cutForPositron, "e+");
                                                   >> 199 
                                                   >> 200   // Set the secondary production cut lower than 990. eV
                                                   >> 201   // Very important for high precision of lowenergy processes at low energies
                                                   >> 202  
                                                   >> 203   G4double lowLimit = 250. * eV;
                                                   >> 204   G4double highLimit = 100. * GeV;
                                                   >> 205   G4ProductionCutsTable::GetProductionCutsTable()->SetEnergyRange(lowLimit, highLimit);
                                                   >> 206   
                                                   >> 207   if (verboseLevel>0) DumpCutValuesTable();
                                                   >> 208 }
166                                                   209