Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/electromagnetic/TestEm17/src/PhysicsList.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/extended/electromagnetic/TestEm17/src/PhysicsList.cc (Version 11.3.0) and /examples/extended/electromagnetic/TestEm17/src/PhysicsList.cc (Version 9.0.p2)


  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 /// \file electromagnetic/TestEm17/src/Physics << 
 27 /// \brief Implementation of the PhysicsList c << 
 28 //                                             << 
 29 //                                                 26 //
                                                   >>  27 // $Id: PhysicsList.cc,v 1.4 2006/10/24 10:06:21 maire Exp $
                                                   >>  28 // GEANT4 tag $Name: geant4-08-03-patch-01 $
 30 //                                                 29 //
 31 //....oooOO0OOooo........oooOO0OOooo........oo     30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 32 //....oooOO0OOooo........oooOO0OOooo........oo     31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 33                                                    32 
 34 #include "PhysicsList.hh"                          33 #include "PhysicsList.hh"
                                                   >>  34 #include "PhysicsListMessenger.hh"
 35                                                    35 
 36 #include "MuNuclearBuilder.hh"                 << 
 37 #include "PhysListEmStandard.hh"                   36 #include "PhysListEmStandard.hh"
 38 #include "PhysicsListMessenger.hh"             <<  37 #include "MuNuclearBuilder.hh"
                                                   >>  38 
                                                   >>  39 #include "G4LossTableManager.hh"
 39                                                    40 
 40 #include "G4BaryonConstructor.hh"              << 
 41 #include "G4BosonConstructor.hh"               << 
 42 #include "G4Electron.hh"                       << 
 43 #include "G4EmParameters.hh"                   << 
 44 #include "G4EmStandardPhysics.hh"              << 
 45 #include "G4EmStandardPhysics_option4.hh"      << 
 46 #include "G4Gamma.hh"                              41 #include "G4Gamma.hh"
 47 #include "G4IonConstructor.hh"                 <<  42 #include "G4Electron.hh"
                                                   >>  43 #include "G4Positron.hh"
                                                   >>  44 
                                                   >>  45 #include "G4BosonConstructor.hh"
 48 #include "G4LeptonConstructor.hh"                  46 #include "G4LeptonConstructor.hh"
 49 #include "G4MesonConstructor.hh"                   47 #include "G4MesonConstructor.hh"
 50 #include "G4Positron.hh"                       <<  48 #include "G4BosonConstructor.hh"
                                                   >>  49 #include "G4BaryonConstructor.hh"
                                                   >>  50 #include "G4IonConstructor.hh"
 51 #include "G4ShortLivedConstructor.hh"              51 #include "G4ShortLivedConstructor.hh"
 52 #include "G4SystemOfUnits.hh"                  <<  52 
 53                                                    53 
 54 //....oooOO0OOooo........oooOO0OOooo........oo     54 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 55                                                    55 
 56 PhysicsList::PhysicsList()                     <<  56 PhysicsList::PhysicsList() : G4VModularPhysicsList()
 57   : G4VModularPhysicsList(), fEmPhysicsList(0) << 
 58 {                                                  57 {
 59   SetVerboseLevel(1);                              58   SetVerboseLevel(1);
 60   fMessenger = new PhysicsListMessenger(this); <<  59   pMessenger = new PhysicsListMessenger(this);
 61                                                    60 
                                                   >>  61   // cuts
                                                   >>  62   currentDefaultCut   = 1.0*mm;
                                                   >>  63   cutForGamma         = currentDefaultCut;
                                                   >>  64   cutForElectron      = currentDefaultCut;
                                                   >>  65   cutForPositron      = currentDefaultCut;
                                                   >>  66   
 62   // EM physics                                    67   // EM physics
 63   fEmName = G4String("emstandard_opt0");       <<  68   emName = G4String("standard");
 64   fEmPhysicsList = new G4EmStandardPhysics();  <<  69   emPhysicsList = new PhysListEmStandard(emName);
 65                                                << 
 66   // extend energy range of PhysicsTables      << 
 67   //                                           << 
 68   G4EmParameters* param = G4EmParameters::Inst << 
 69   param->SetMinEnergy(100 * eV);               << 
 70   param->SetMaxEnergy(1000 * PeV);             << 
 71                                                    70 
 72   fMuNuclPhysicsList = 0;                      <<  71   muNuclPhysicsList = 0;
                                                   >>  72   
                                                   >>  73   // instanciate EnergyLossTable
                                                   >>  74   G4LossTableManager::Instance();
 73 }                                                  75 }
 74                                                    76 
 75 //....oooOO0OOooo........oooOO0OOooo........oo     77 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 76                                                    78 
 77 PhysicsList::~PhysicsList()                        79 PhysicsList::~PhysicsList()
 78 {                                                  80 {
 79   delete fMessenger;                           <<  81   delete pMessenger;
 80 }                                                  82 }
 81                                                    83 
 82 //....oooOO0OOooo........oooOO0OOooo........oo     84 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 83                                                    85 
 84 void PhysicsList::ConstructParticle()              86 void PhysicsList::ConstructParticle()
 85 {                                                  87 {
 86   G4BosonConstructor pBosonConstructor;        <<  88   G4BosonConstructor  pBosonConstructor;
 87   pBosonConstructor.ConstructParticle();           89   pBosonConstructor.ConstructParticle();
 88                                                    90 
 89   G4LeptonConstructor pLeptonConstructor;          91   G4LeptonConstructor pLeptonConstructor;
 90   pLeptonConstructor.ConstructParticle();          92   pLeptonConstructor.ConstructParticle();
 91                                                    93 
 92   G4MesonConstructor pMesonConstructor;            94   G4MesonConstructor pMesonConstructor;
 93   pMesonConstructor.ConstructParticle();           95   pMesonConstructor.ConstructParticle();
 94                                                    96 
 95   G4BaryonConstructor pBaryonConstructor;          97   G4BaryonConstructor pBaryonConstructor;
 96   pBaryonConstructor.ConstructParticle();          98   pBaryonConstructor.ConstructParticle();
 97                                                    99 
 98   G4IonConstructor pIonConstructor;               100   G4IonConstructor pIonConstructor;
 99   pIonConstructor.ConstructParticle();            101   pIonConstructor.ConstructParticle();
100                                                   102 
101   G4ShortLivedConstructor pShortLivedConstruct    103   G4ShortLivedConstructor pShortLivedConstructor;
102   pShortLivedConstructor.ConstructParticle();     104   pShortLivedConstructor.ConstructParticle();
103 }                                                 105 }
104                                                   106 
105 //....oooOO0OOooo........oooOO0OOooo........oo    107 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
106                                                   108 
107 void PhysicsList::ConstructProcess()              109 void PhysicsList::ConstructProcess()
108 {                                                 110 {
109   // transportation                               111   // transportation
110   //                                              112   //
111   AddTransportation();                            113   AddTransportation();
112                                                << 114   
113   // electromagnetic Physics List                 115   // electromagnetic Physics List
114   //                                              116   //
115   fEmPhysicsList->ConstructProcess();          << 117   emPhysicsList->ConstructProcess();
116   if (fMuNuclPhysicsList) fMuNuclPhysicsList-> << 118   if(muNuclPhysicsList) muNuclPhysicsList->ConstructProcess();
117 }                                                 119 }
118                                                   120 
119 //....oooOO0OOooo........oooOO0OOooo........oo    121 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
120                                                   122 
121 void PhysicsList::AddPhysicsList(const G4Strin    123 void PhysicsList::AddPhysicsList(const G4String& name)
122 {                                                 124 {
123   if (verboseLevel > -1) {                     << 125   if (verboseLevel>1) {
124     G4cout << "PhysicsList::AddPhysicsList: <"    126     G4cout << "PhysicsList::AddPhysicsList: <" << name << ">" << G4endl;
125   }                                               127   }
126                                                   128 
127   if (name == fEmName) return;                 << 129   if (name == emName) return;
                                                   >> 130 
                                                   >> 131   if (name == "standard" && name != emName) {
                                                   >> 132 
                                                   >> 133     emName = name;
                                                   >> 134     delete emPhysicsList;
                                                   >> 135     emPhysicsList = new PhysListEmStandard(name);
                                                   >> 136 
                                                   >> 137   } else if (name == "muNucl") {
                                                   >> 138     muNuclPhysicsList = new MuNuclearBuilder(name);
                                                   >> 139 
                                                   >> 140   } else {
128                                                   141 
129   if (name == "emstandard_opt0") {             << 
130     fEmName = name;                            << 
131     delete fEmPhysicsList;                     << 
132     fEmPhysicsList = new G4EmStandardPhysics() << 
133   }                                            << 
134   else if (name == "emstandard_opt4") {        << 
135     fEmName = name;                            << 
136     delete fEmPhysicsList;                     << 
137     fEmPhysicsList = new G4EmStandardPhysics_o << 
138   }                                            << 
139   else if (name == "local") {                  << 
140     fEmName = name;                            << 
141     delete fEmPhysicsList;                     << 
142     fEmPhysicsList = new PhysListEmStandard(na << 
143   }                                            << 
144   else if (name == "muNucl") {                 << 
145     fMuNuclPhysicsList = new MuNuclearBuilder( << 
146   }                                            << 
147   else {                                       << 
148     G4cout << "PhysicsList::AddPhysicsList: <"    142     G4cout << "PhysicsList::AddPhysicsList: <" << name << ">"
149            << " is not defined" << G4endl;     << 143            << " is not defined"
                                                   >> 144            << G4endl;
150   }                                               145   }
151 }                                                 146 }
152                                                   147 
153 //....oooOO0OOooo........oooOO0OOooo........oo    148 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
                                                   >> 149 
                                                   >> 150 void PhysicsList::SetCuts()
                                                   >> 151 {
                                                   >> 152      
                                                   >> 153   if (verboseLevel >0){
                                                   >> 154     G4cout << "PhysicsList::SetCuts:";
                                                   >> 155     G4cout << "CutLength : " << G4BestUnit(defaultCutValue,"Length") << G4endl;
                                                   >> 156   }  
                                                   >> 157 
                                                   >> 158   // set cut values for gamma at first and for e- second and next for e+,
                                                   >> 159   // because some processes for e+/e- need cut values for gamma
                                                   >> 160   SetCutValue(cutForGamma, "gamma");
                                                   >> 161   SetCutValue(cutForElectron, "e-");
                                                   >> 162   SetCutValue(cutForPositron, "e+");   
                                                   >> 163     
                                                   >> 164   if (verboseLevel>0) DumpCutValuesTable();
                                                   >> 165 }
                                                   >> 166 
                                                   >> 167 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
                                                   >> 168 
                                                   >> 169 #include "G4Gamma.hh"
                                                   >> 170 #include "G4Electron.hh"
                                                   >> 171 #include "G4Positron.hh"
                                                   >> 172 
                                                   >> 173 void PhysicsList::SetCutForGamma(G4double cut)
                                                   >> 174 {
                                                   >> 175   cutForGamma = cut;
                                                   >> 176   SetParticleCuts(cutForGamma, G4Gamma::Gamma());
                                                   >> 177 }
                                                   >> 178 
                                                   >> 179 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
                                                   >> 180 
                                                   >> 181 void PhysicsList::SetCutForElectron(G4double cut)
                                                   >> 182 {
                                                   >> 183   cutForElectron = cut;
                                                   >> 184   SetParticleCuts(cutForElectron, G4Electron::Electron());
                                                   >> 185 }
                                                   >> 186 
                                                   >> 187 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
                                                   >> 188 
                                                   >> 189 void PhysicsList::SetCutForPositron(G4double cut)
                                                   >> 190 {
                                                   >> 191   cutForPositron = cut;
                                                   >> 192   SetParticleCuts(cutForPositron, G4Positron::Positron());
                                                   >> 193 }
                                                   >> 194 
                                                   >> 195 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
                                                   >> 196 
154                                                   197