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 10.1)


  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     26 /// \file electromagnetic/TestEm17/src/PhysicsList.cc
 27 /// \brief Implementation of the PhysicsList c     27 /// \brief Implementation of the PhysicsList class
 28 //                                                 28 //
 29 //                                                 29 //
                                                   >>  30 // $Id: PhysicsList.cc 85311 2014-10-27 14:23:25Z gcosmo $
 30 //                                                 31 //
 31 //....oooOO0OOooo........oooOO0OOooo........oo     32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 32 //....oooOO0OOooo........oooOO0OOooo........oo     33 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 33                                                    34 
 34 #include "PhysicsList.hh"                          35 #include "PhysicsList.hh"
                                                   >>  36 #include "PhysicsListMessenger.hh"
 35                                                    37 
 36 #include "MuNuclearBuilder.hh"                 <<  38 #include "G4EmStandardPhysics.hh"
 37 #include "PhysListEmStandard.hh"                   39 #include "PhysListEmStandard.hh"
 38 #include "PhysicsListMessenger.hh"             <<  40 #include "MuNuclearBuilder.hh"
 39                                                    41 
 40 #include "G4BaryonConstructor.hh"              << 
 41 #include "G4BosonConstructor.hh"               << 
 42 #include "G4Electron.hh"                       << 
 43 #include "G4EmParameters.hh"                       42 #include "G4EmParameters.hh"
 44 #include "G4EmStandardPhysics.hh"              <<  43 
 45 #include "G4EmStandardPhysics_option4.hh"      << 
 46 #include "G4Gamma.hh"                              44 #include "G4Gamma.hh"
 47 #include "G4IonConstructor.hh"                 <<  45 #include "G4Electron.hh"
                                                   >>  46 #include "G4Positron.hh"
                                                   >>  47 
                                                   >>  48 #include "G4BosonConstructor.hh"
 48 #include "G4LeptonConstructor.hh"                  49 #include "G4LeptonConstructor.hh"
 49 #include "G4MesonConstructor.hh"                   50 #include "G4MesonConstructor.hh"
 50 #include "G4Positron.hh"                       <<  51 #include "G4BosonConstructor.hh"
                                                   >>  52 #include "G4BaryonConstructor.hh"
                                                   >>  53 #include "G4IonConstructor.hh"
 51 #include "G4ShortLivedConstructor.hh"              54 #include "G4ShortLivedConstructor.hh"
                                                   >>  55 
 52 #include "G4SystemOfUnits.hh"                      56 #include "G4SystemOfUnits.hh"
 53                                                    57 
 54 //....oooOO0OOooo........oooOO0OOooo........oo     58 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 55                                                    59 
 56 PhysicsList::PhysicsList()                     <<  60 PhysicsList::PhysicsList() : G4VModularPhysicsList(),
 57   : G4VModularPhysicsList(), fEmPhysicsList(0) <<  61   fEmPhysicsList(0),
                                                   >>  62   fMuNuclPhysicsList(0),
                                                   >>  63   fMessenger(0)
 58 {                                                  64 {
 59   SetVerboseLevel(1);                              65   SetVerboseLevel(1);
 60   fMessenger = new PhysicsListMessenger(this);     66   fMessenger = new PhysicsListMessenger(this);
 61                                                    67 
                                                   >>  68   //extend energy range of PhysicsTables
                                                   >>  69   //
                                                   >>  70   G4EmParameters* param = G4EmParameters::Instance();
                                                   >>  71   param->SetMinEnergy(100*eV);  
                                                   >>  72   param->SetMaxEnergy(1000*PeV);
                                                   >>  73 
 62   // EM physics                                    74   // EM physics
 63   fEmName = G4String("emstandard_opt0");           75   fEmName = G4String("emstandard_opt0");
 64   fEmPhysicsList = new G4EmStandardPhysics();      76   fEmPhysicsList = new G4EmStandardPhysics();
 65                                                    77 
 66   // extend energy range of PhysicsTables      << 
 67   //                                           << 
 68   G4EmParameters* param = G4EmParameters::Inst << 
 69   param->SetMinEnergy(100 * eV);               << 
 70   param->SetMaxEnergy(1000 * PeV);             << 
 71                                                << 
 72   fMuNuclPhysicsList = 0;                          78   fMuNuclPhysicsList = 0;
 73 }                                                  79 }
 74                                                    80 
 75 //....oooOO0OOooo........oooOO0OOooo........oo     81 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 76                                                    82 
 77 PhysicsList::~PhysicsList()                        83 PhysicsList::~PhysicsList()
 78 {                                                  84 {
 79   delete fMessenger;                               85   delete fMessenger;
 80 }                                                  86 }
 81                                                    87 
 82 //....oooOO0OOooo........oooOO0OOooo........oo     88 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 83                                                    89 
 84 void PhysicsList::ConstructParticle()              90 void PhysicsList::ConstructParticle()
 85 {                                                  91 {
 86   G4BosonConstructor pBosonConstructor;        <<  92   G4BosonConstructor  pBosonConstructor;
 87   pBosonConstructor.ConstructParticle();           93   pBosonConstructor.ConstructParticle();
 88                                                    94 
 89   G4LeptonConstructor pLeptonConstructor;          95   G4LeptonConstructor pLeptonConstructor;
 90   pLeptonConstructor.ConstructParticle();          96   pLeptonConstructor.ConstructParticle();
 91                                                    97 
 92   G4MesonConstructor pMesonConstructor;            98   G4MesonConstructor pMesonConstructor;
 93   pMesonConstructor.ConstructParticle();           99   pMesonConstructor.ConstructParticle();
 94                                                   100 
 95   G4BaryonConstructor pBaryonConstructor;         101   G4BaryonConstructor pBaryonConstructor;
 96   pBaryonConstructor.ConstructParticle();         102   pBaryonConstructor.ConstructParticle();
 97                                                   103 
 98   G4IonConstructor pIonConstructor;               104   G4IonConstructor pIonConstructor;
 99   pIonConstructor.ConstructParticle();            105   pIonConstructor.ConstructParticle();
100                                                   106 
101   G4ShortLivedConstructor pShortLivedConstruct    107   G4ShortLivedConstructor pShortLivedConstructor;
102   pShortLivedConstructor.ConstructParticle();     108   pShortLivedConstructor.ConstructParticle();
103 }                                                 109 }
104                                                   110 
105 //....oooOO0OOooo........oooOO0OOooo........oo    111 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
106                                                   112 
107 void PhysicsList::ConstructProcess()              113 void PhysicsList::ConstructProcess()
108 {                                                 114 {
109   // transportation                               115   // transportation
110   //                                              116   //
111   AddTransportation();                            117   AddTransportation();
112                                                << 118   
113   // electromagnetic Physics List                 119   // electromagnetic Physics List
114   //                                              120   //
115   fEmPhysicsList->ConstructProcess();             121   fEmPhysicsList->ConstructProcess();
116   if (fMuNuclPhysicsList) fMuNuclPhysicsList-> << 122   if(fMuNuclPhysicsList) fMuNuclPhysicsList->ConstructProcess();
117 }                                                 123 }
118                                                   124 
119 //....oooOO0OOooo........oooOO0OOooo........oo    125 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
120                                                   126 
121 void PhysicsList::AddPhysicsList(const G4Strin    127 void PhysicsList::AddPhysicsList(const G4String& name)
122 {                                                 128 {
123   if (verboseLevel > -1) {                     << 129   if (verboseLevel>1) {
124     G4cout << "PhysicsList::AddPhysicsList: <"    130     G4cout << "PhysicsList::AddPhysicsList: <" << name << ">" << G4endl;
125   }                                               131   }
126                                                   132 
127   if (name == fEmName) return;                    133   if (name == fEmName) return;
128                                                   134 
129   if (name == "emstandard_opt0") {                135   if (name == "emstandard_opt0") {
                                                   >> 136 
130     fEmName = name;                               137     fEmName = name;
131     delete fEmPhysicsList;                        138     delete fEmPhysicsList;
132     fEmPhysicsList = new G4EmStandardPhysics() << 139     fEmPhysicsList = new PhysListEmStandard(name);
133   }                                            << 140 
134   else if (name == "emstandard_opt4") {        << 141   } else if (name == "local") {
135     fEmName = name;                            << 142 
136     delete fEmPhysicsList;                     << 
137     fEmPhysicsList = new G4EmStandardPhysics_o << 
138   }                                            << 
139   else if (name == "local") {                  << 
140     fEmName = name;                               143     fEmName = name;
141     delete fEmPhysicsList;                        144     delete fEmPhysicsList;
142     fEmPhysicsList = new PhysListEmStandard(na    145     fEmPhysicsList = new PhysListEmStandard(name);
143   }                                            << 146 
144   else if (name == "muNucl") {                 << 147   } else if (name == "muNucl") {
145     fMuNuclPhysicsList = new MuNuclearBuilder(    148     fMuNuclPhysicsList = new MuNuclearBuilder(name);
146   }                                            << 149 
147   else {                                       << 150   } else {
                                                   >> 151 
148     G4cout << "PhysicsList::AddPhysicsList: <"    152     G4cout << "PhysicsList::AddPhysicsList: <" << name << ">"
149            << " is not defined" << G4endl;     << 153            << " is not defined"
                                                   >> 154            << G4endl;
150   }                                               155   }
151 }                                                 156 }
152                                                   157 
153 //....oooOO0OOooo........oooOO0OOooo........oo    158 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
                                                   >> 159 
                                                   >> 160 
154                                                   161