Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/error_propagation/src/G4ErrorPhysicsList.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 /error_propagation/src/G4ErrorPhysicsList.cc (Version 11.3.0) and /error_propagation/src/G4ErrorPhysicsList.cc (Version 9.2.p4)


  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 // $Id: G4ErrorPhysicsList.cc,v 1.2 2007/05/29 14:41:35 gcosmo Exp $
                                                   >>  27 // GEANT4 tag $Name: geant4-09-02-patch-04 $
 26 //                                                 28 //
 27 // -------------------------------------------     29 // ------------------------------------------------------------
 28 //      GEANT 4 class implementation file      <<  30 //      GEANT 4 class implementation file 
 29 // -------------------------------------------     31 // ------------------------------------------------------------
 30                                                    32 
 31 #include "globals.hh"                              33 #include "globals.hh"
 32                                                    34 
 33 #include "G4PhysicalConstants.hh"              << 
 34 #include "G4SystemOfUnits.hh"                  << 
 35                                                << 
 36 #include "G4ErrorPhysicsList.hh"                   35 #include "G4ErrorPhysicsList.hh"
 37 #include "G4ComptonScattering.hh"                  36 #include "G4ComptonScattering.hh"
 38 #include "G4GammaConversion.hh"                    37 #include "G4GammaConversion.hh"
 39 #include "G4PhotoElectricEffect.hh"                38 #include "G4PhotoElectricEffect.hh"
 40                                                <<  39  
 41 #include "G4eIonisation.hh"                        40 #include "G4eIonisation.hh"
 42 #include "G4eBremsstrahlung.hh"                    41 #include "G4eBremsstrahlung.hh"
 43 #include "G4eplusAnnihilation.hh"                  42 #include "G4eplusAnnihilation.hh"
 44                                                    43 
 45 #include "G4MuIonisation.hh"                       44 #include "G4MuIonisation.hh"
 46 #include "G4MuBremsstrahlung.hh"                   45 #include "G4MuBremsstrahlung.hh"
 47 #include "G4MuPairProduction.hh"                   46 #include "G4MuPairProduction.hh"
 48                                                    47 
 49 #include "G4hIonisation.hh"                        48 #include "G4hIonisation.hh"
 50                                                    49 
 51 #include "G4MuIonisation.hh"                       50 #include "G4MuIonisation.hh"
 52 #include "G4MuBremsstrahlung.hh"                   51 #include "G4MuBremsstrahlung.hh"
 53 #include "G4MuPairProduction.hh"                   52 #include "G4MuPairProduction.hh"
 54                                                    53 
 55 #include "G4hIonisation.hh"                        54 #include "G4hIonisation.hh"
 56                                                    55 
 57 #include "G4ParticleDefinition.hh"                 56 #include "G4ParticleDefinition.hh"
 58 #include "G4ProcessManager.hh"                     57 #include "G4ProcessManager.hh"
 59 #include "G4ProcessVector.hh"                      58 #include "G4ProcessVector.hh"
 60 #include "G4ParticleTypes.hh"                      59 #include "G4ParticleTypes.hh"
 61 #include "G4ParticleTable.hh"                      60 #include "G4ParticleTable.hh"
 62 #include "G4Material.hh"                           61 #include "G4Material.hh"
 63 #include "G4MaterialTable.hh"                      62 #include "G4MaterialTable.hh"
 64 #include "G4ios.hh"                                63 #include "G4ios.hh"
 65 #include "G4PhysicsTable.hh"                       64 #include "G4PhysicsTable.hh"
 66 #include "G4Transportation.hh"                     65 #include "G4Transportation.hh"
 67                                                    66 
 68 #include "G4ErrorEnergyLoss.hh"                    67 #include "G4ErrorEnergyLoss.hh"
 69                                                    68 
 70 //--------------------------------------------     69 //------------------------------------------------------------------------
 71 G4ErrorPhysicsList::G4ErrorPhysicsList()       <<  70 G4ErrorPhysicsList::G4ErrorPhysicsList():  G4VUserPhysicsList()
 72   : G4VUserPhysicsList()                       << 
 73 {                                                  71 {
 74   defaultCutValue =                            <<  72   defaultCutValue = 1.0E+9*cm; // set big step so that AlongStep computes all the energy
 75     1.0E+9 * cm;  // set big step so that Alon << 
 76 }                                                  73 }
 77                                                    74 
                                                   >>  75 
 78 //--------------------------------------------     76 //------------------------------------------------------------------------
 79 G4ErrorPhysicsList::~G4ErrorPhysicsList() {}   <<  77 G4ErrorPhysicsList::~G4ErrorPhysicsList()
                                                   >>  78 {
                                                   >>  79 }
                                                   >>  80 
 80                                                    81 
 81 //--------------------------------------------     82 //------------------------------------------------------------------------
 82 void G4ErrorPhysicsList::ConstructParticle()       83 void G4ErrorPhysicsList::ConstructParticle()
 83 {                                                  84 {
 84   // In this method, static member functions s <<  85 // In this method, static member functions should be called
 85   // for all particles which you want to use.      86   // for all particles which you want to use.
 86   // This ensures that objects of these partic     87   // This ensures that objects of these particle types will be
 87   // created in the program.                   <<  88   // created in the program. 
 88   //  gamma                                        89   //  gamma
 89   G4Gamma::GammaDefinition();                  <<  90   G4Gamma::GammaDefinition(); 
 90   //  e+/-                                         91   //  e+/-
 91   G4Electron::ElectronDefinition();                92   G4Electron::ElectronDefinition();
 92   G4Positron::PositronDefinition();                93   G4Positron::PositronDefinition();
 93   // mu+/-                                         94   // mu+/-
 94   G4MuonPlus::MuonPlusDefinition();                95   G4MuonPlus::MuonPlusDefinition();
 95   G4MuonMinus::MuonMinusDefinition();              96   G4MuonMinus::MuonMinusDefinition();
 96                                                    97 
 97   // pi+/-                                     << 
 98   G4PionPlus::PionPlusDefinition();            << 
 99   G4PionMinus::PionMinusDefinition();          << 
100                                                << 
101   // proton                                    << 
102   G4Proton::ProtonDefinition();                << 
103 }                                                  98 }
104                                                    99 
                                                   >> 100 
105 //--------------------------------------------    101 //------------------------------------------------------------------------
106 void G4ErrorPhysicsList::ConstructProcess()       102 void G4ErrorPhysicsList::ConstructProcess()
107 {                                                 103 {
108   G4Transportation* theTransportationProcess = << 104   G4Transportation* theTransportationProcess= new G4Transportation();
109                                                   105 
110 #ifdef G4VERBOSE                                  106 #ifdef G4VERBOSE
111   if(verboseLevel >= 4)                        << 107     if (verboseLevel >= 4){
112   {                                            << 108       G4cout << "G4VUserPhysicsList::ConstructProcess()  "<< G4endl;
113     G4cout << "G4VUserPhysicsList::ConstructPr << 109     }
114   }                                            << 
115 #endif                                            110 #endif
116                                                   111 
117   // loop over all particles in G4ParticleTabl    112   // loop over all particles in G4ParticleTable
118   auto myParticleIterator = GetParticleIterato << 113   theParticleIterator->reset();
119   myParticleIterator->reset();                 << 114   while( (*theParticleIterator)() ){
120   while((*myParticleIterator)())               << 115     G4ParticleDefinition* particle = theParticleIterator->value();
121   {  // Loop checking, 06.08.2015, G.Cosmo     << 116     G4ProcessManager* pmanager = particle->GetProcessManager();
122     G4ParticleDefinition* particle = myParticl << 117     if (!particle->IsShortLived()) {
123     G4ProcessManager* pmanager     = particle- << 118       G4cout << particle << "G4ErrorPhysicsList:: particle process manager " << particle->GetParticleName() << " = " << particle->GetProcessManager() << G4endl;
124     if(!particle->IsShortLived())              << 
125     {                                          << 
126       G4cout << particle << "G4ErrorPhysicsLis << 
127              << particle->GetParticleName() << << 
128              << particle->GetProcessManager()  << 
129       // Add transportation process for all pa    119       // Add transportation process for all particles other than  "shortlived"
130       if(pmanager == 0)                        << 120       if ( pmanager == 0) {
131       {                                        << 
132         // Error !! no process manager            121         // Error !! no process manager
133         G4String particleName = particle->GetP    122         G4String particleName = particle->GetParticleName();
134         G4Exception("G4ErrorPhysicsList::Const << 123         G4Exception("G4ErrorPhysicsList::ConstructProcess","No process manager",
135                     "No process manager", RunM << 124                     RunMustBeAborted, particleName );
136       }                                        << 125       } else {
137       else                                     << 
138       {                                        << 
139         // add transportation with ordering =     126         // add transportation with ordering = ( -1, "first", "first" )
140         pmanager->AddProcess(theTransportation << 127         pmanager ->AddProcess(theTransportationProcess);
141         pmanager->SetProcessOrderingToFirst(th << 128         pmanager ->SetProcessOrderingToFirst(theTransportationProcess, idxAlongStep);
142                                             id << 129         pmanager ->SetProcessOrderingToFirst(theTransportationProcess, idxPostStep);
143         pmanager->SetProcessOrderingToFirst(th << 
144                                             id << 
145       }                                           130       }
146     }                                          << 131     } else {
147     else                                       << 
148     {                                          << 
149       // shortlived particle case                 132       // shortlived particle case
150     }                                             133     }
151   }                                               134   }
152                                                   135 
153   ConstructEM();                                  136   ConstructEM();
154 }                                                 137 }
155                                                   138 
                                                   >> 139 
156 //--------------------------------------------    140 //------------------------------------------------------------------------
157 #include "G4eBremsstrahlung.hh"                   141 #include "G4eBremsstrahlung.hh"
158 #include "G4eIonisation.hh"                       142 #include "G4eIonisation.hh"
159                                                   143 
160 #include "G4eIonisation.hh"                       144 #include "G4eIonisation.hh"
161                                                   145 
162 #include "G4MuBremsstrahlung.hh"                  146 #include "G4MuBremsstrahlung.hh"
163 #include "G4MuIonisation.hh"                      147 #include "G4MuIonisation.hh"
164 #include "G4MuPairProduction.hh"                  148 #include "G4MuPairProduction.hh"
165                                                   149 
166 #include "G4PhysicsTable.hh"                      150 #include "G4PhysicsTable.hh"
167                                                   151 
                                                   >> 152 #include "G4VEnergyLoss.hh"
                                                   >> 153 
168 #include "G4MuIonisation.hh"                      154 #include "G4MuIonisation.hh"
169                                                   155 
170 #include "G4ErrorStepLengthLimitProcess.hh"       156 #include "G4ErrorStepLengthLimitProcess.hh"
171 #include "G4ErrorMagFieldLimitProcess.hh"         157 #include "G4ErrorMagFieldLimitProcess.hh"
172 #include "G4ErrorMessenger.hh"                    158 #include "G4ErrorMessenger.hh"
173                                                   159 
174 void G4ErrorPhysicsList::ConstructEM()            160 void G4ErrorPhysicsList::ConstructEM()
175 {                                                 161 {
                                                   >> 162 
176   G4ErrorEnergyLoss* eLossProcess = new G4Erro    163   G4ErrorEnergyLoss* eLossProcess = new G4ErrorEnergyLoss;
177   G4ErrorStepLengthLimitProcess* stepLengthLim << 164   G4ErrorStepLengthLimitProcess* stepLengthLimitProcess = new G4ErrorStepLengthLimitProcess;
178     new G4ErrorStepLengthLimitProcess;         << 165   G4ErrorMagFieldLimitProcess* magFieldLimitProcess = new G4ErrorMagFieldLimitProcess;
179   G4ErrorMagFieldLimitProcess* magFieldLimitPr << 166   new G4ErrorMessenger( stepLengthLimitProcess, magFieldLimitProcess, eLossProcess );
180     new G4ErrorMagFieldLimitProcess;           << 167 
181   new G4ErrorMessenger(stepLengthLimitProcess, << 168   theParticleIterator->reset();
182                        eLossProcess);          << 169   while( (*theParticleIterator)() ){
183                                                << 170     G4ParticleDefinition* particle = theParticleIterator->value();
184   auto myParticleIterator = GetParticleIterato << 171     G4ProcessManager* pmanager = particle->GetProcessManager();
185   myParticleIterator->reset();                 << 172     G4String particleName = particle->GetParticleName();
186   while((*myParticleIterator)())               << 173      
187   {  // Loop checking, 06.08.2015, G.Cosmo     << 174     if (particleName == "gamma") {
188     G4ParticleDefinition* particle = myParticl << 175     // gamma
189     G4ProcessManager* pmanager     = particle- << 
190     G4String particleName          = particle- << 
191                                                << 
192     if(particleName == "gamma")                << 
193     {                                          << 
194       // gamma                                 << 
195       pmanager->AddDiscreteProcess(new G4Gamma    176       pmanager->AddDiscreteProcess(new G4GammaConversion());
196       pmanager->AddDiscreteProcess(new G4Compt << 177       pmanager->AddDiscreteProcess(new G4ComptonScattering());      
197       pmanager->AddDiscreteProcess(new G4Photo    178       pmanager->AddDiscreteProcess(new G4PhotoElectricEffect());
198                                                   179 
199       //    } else if (particleName == "e-" ||    180       //    } else if (particleName == "e-" || particleName == "e+"
200       //               || particleName == "mu+    181       //               || particleName == "mu+" || particleName == "mu-" ) {
201     }                                          << 182     }else if (!particle->IsShortLived() && particle->GetPDGCharge() != 0 ) {
202     else if(!particle->IsShortLived() && parti << 183   
203     {                                          << 184       pmanager->AddContinuousProcess(eLossProcess,1);
204       pmanager->AddContinuousProcess(eLossProc << 185       pmanager->AddDiscreteProcess( stepLengthLimitProcess, 2 ); 
205       pmanager->AddDiscreteProcess(stepLengthL << 186       pmanager->AddDiscreteProcess( magFieldLimitProcess, 3 );     
206       pmanager->AddDiscreteProcess(magFieldLim << 187       
207                                                << 188     } else if ((!particle->IsShortLived()) &&
208       /*     } else if ((!particle->IsShortLiv << 189                (particle->GetPDGCharge() != 0.0) && 
209                (particle->GetPDGCharge() != 0. << 
210                (particle->GetParticleName() !=    190                (particle->GetParticleName() != "chargedgeantino")) {
211      // all others charged particles except ge    191      // all others charged particles except geantino
212       //   G4VProcess* aMultipleScattering = n    192       //   G4VProcess* aMultipleScattering = new G4MultipleScattering();
213      G4VProcess* anIonisation        = new G4h << 193      G4VProcess* anIonisation        = new G4hIonisation();     
214      ////G4VProcess*  theUserCuts = new G4User    194      ////G4VProcess*  theUserCuts = new G4UserSpecialCuts();
215                                                << 195      
216      //                                           196      //
217      // add processes                             197      // add processes
218      pmanager->AddProcess(anIonisation);          198      pmanager->AddProcess(anIonisation);
219      //   pmanager->AddProcess(aMultipleScatte << 199      //   pmanager->AddProcess(aMultipleScattering);    
220      ////pmanager->AddProcess(theUserCuts);       200      ////pmanager->AddProcess(theUserCuts);
221                                                << 201      
222      //                                           202      //
223      // set ordering for AlongStepDoIt            203      // set ordering for AlongStepDoIt
224      //   pmanager->SetProcessOrdering(aMultip    204      //   pmanager->SetProcessOrdering(aMultipleScattering, idxAlongStep,1);
225      pmanager->SetProcessOrdering(anIonisation    205      pmanager->SetProcessOrdering(anIonisation, idxAlongStep,1);
226                                                << 206      
227      //                                           207      //
228      // set ordering for PostStepDoIt             208      // set ordering for PostStepDoIt
229      //   pmanager->SetProcessOrdering(aMultip    209      //   pmanager->SetProcessOrdering(aMultipleScattering, idxPostStep,1);
230      pmanager->SetProcessOrdering(anIonisation    210      pmanager->SetProcessOrdering(anIonisation, idxPostStep,1);
231      ////pmanager->SetProcessOrdering(theUserC    211      ////pmanager->SetProcessOrdering(theUserCuts,     idxPostStep,2);
232      */                                        << 
233     }                                             212     }
234   }                                               213   }
235 }                                                 214 }
236                                                   215 
                                                   >> 216 
237 //--------------------------------------------    217 //------------------------------------------------------------------------
238 void G4ErrorPhysicsList::SetCuts()                218 void G4ErrorPhysicsList::SetCuts()
239 {                                                 219 {
240   //  " G4VUserPhysicsList::SetCutsWithDefault << 220   //  " G4VUserPhysicsList::SetCutsWithDefault" method sets 
241   //   the default cut value or all particle t << 221   //   the default cut value or all particle types 
242   SetCutsWithDefault();                        << 222   SetCutsWithDefault(); 
243   // if (verboseLevel>0)                       << 223   // if (verboseLevel>0) 
244   //  DumpCutValuesTable();                       224   //  DumpCutValuesTable();
245 }                                                 225 }
                                                   >> 226 
246                                                   227