Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/standard/src/G4eBremsstrahlung.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 /processes/electromagnetic/standard/src/G4eBremsstrahlung.cc (Version 11.3.0) and /processes/electromagnetic/standard/src/G4eBremsstrahlung.cc (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 // $Id: G4eBremsstrahlung.cc,v 1.56 2009-02-20 12:06:37 vnivanch Exp $
                                                   >>  27 // GEANT4 tag $Name: not supported by cvs2svn $
 26 //                                                 28 //
 27 // -------------------------------------------     29 // -------------------------------------------------------------------
 28 //                                                 30 //
 29 // GEANT4 Class file                               31 // GEANT4 Class file
 30 //                                                 32 //
 31 //                                                 33 //
 32 // File name:     G4eBremsstrahlung                34 // File name:     G4eBremsstrahlung
 33 //                                                 35 //
 34 // Author:        Michel Maire                     36 // Author:        Michel Maire
 35 //                                                 37 //
 36 // Creation date: 26.06.1996                       38 // Creation date: 26.06.1996
 37 //                                                 39 //
 38 // Modified by Michel Maire, Vladimir Ivanchen <<  40 // Modifications:
 39 //                                                 41 //
                                                   >>  42 // 26-09-96 extension of the total crosssection above 100 GeV, M.Maire
                                                   >>  43 //  1-10-96 new type G4OrderedTable; ComputePartialSumSigma(), M.Maire
                                                   >>  44 // 16-10-96 DoIt() call to the non static GetEnergyCuts(), L.Urban
                                                   >>  45 // 13-12-96 Sign corrected in grejmax and greject
                                                   >>  46 //          error definition of screenvar, L.Urban
                                                   >>  47 // 20-03-97 new energy loss+ionisation+brems scheme, L.Urban
                                                   >>  48 // 07-04-98 remove 'tracking cut' of the diffracted particle, MMa
                                                   >>  49 // 13-08-98 new methods SetBining() PrintInfo()
                                                   >>  50 // 03-03-99 Bug fixed in LPM effect, L.Urban
                                                   >>  51 // 10-02-00 modifications , new e.m. structure, L.Urban
                                                   >>  52 // 07-08-00 new cross section/en.loss parametrisation, LPM flag , L.Urban
                                                   >>  53 // 21-09-00 corrections in the LPM implementation, L.Urban
                                                   >>  54 // 28-05-01 V.Ivanchenko minor changes to provide ANSI -wall compilation
                                                   >>  55 // 09-08-01 new methods Store/Retrieve PhysicsTable (mma)
                                                   >>  56 // 17-09-01 migration of Materials to pure STL (mma)
                                                   >>  57 // 21-09-01 completion of RetrievePhysicsTable() (mma)
                                                   >>  58 // 29-10-01 all static functions no more inlined (mma)
                                                   >>  59 // 08-11-01 particleMass becomes a local variable
                                                   >>  60 // 30-04-02 V.Ivanchenko update to new design
                                                   >>  61 // 23-12-02 Change interface in order to move to cut per region (VI)
                                                   >>  62 // 26-12-02 Secondary production moved to derived classes (VI)
                                                   >>  63 // 23-05-03 Define default integral + BohrFluctuations (V.Ivanchenko)
                                                   >>  64 // 08-08-03 STD substitute standard  (V.Ivanchenko)
                                                   >>  65 // 12-11-03 G4EnergyLossSTD -> G4EnergyLossProcess (V.Ivanchenko)
                                                   >>  66 // 04-11-04 add gamma threshold (V.Ivanchenko)
                                                   >>  67 // 08-11-04 Migration to new interface of Store/Retrieve tables (V.Ivantchenko)
                                                   >>  68 // 08-04-05 Major optimisation of internal interfaces (V.Ivantchenko)
                                                   >>  69 // 22-05-06 Use gammaThreshold from manager (V.Ivantchenko)
                                                   >>  70 // 15-01-07 use SetEmModel() from G4VEnergyLossProcess (mma)
                                                   >>  71 //          use RelEmModel above 1GeV (AS &  VI)  
                                                   >>  72 // 13-11-08 reenable LPM switch (A.Schaelicke)
 40 // -------------------------------------------     73 // -------------------------------------------------------------------
 41 //                                                 74 //
 42 //....oooOO0OOooo........oooOO0OOooo........oo     75 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 43 //....oooOO0OOooo........oooOO0OOooo........oo     76 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 44                                                    77 
 45 #include "G4eBremsstrahlung.hh"                    78 #include "G4eBremsstrahlung.hh"
 46 #include "G4SystemOfUnits.hh"                  << 
 47 #include "G4Gamma.hh"                              79 #include "G4Gamma.hh"
                                                   >>  80 #include "G4eBremsstrahlungModel.hh"
 48 #include "G4SeltzerBergerModel.hh"                 81 #include "G4SeltzerBergerModel.hh"
 49 #include "G4eBremsstrahlungRelModel.hh"            82 #include "G4eBremsstrahlungRelModel.hh"
 50 #include "G4UnitsTable.hh"                         83 #include "G4UnitsTable.hh"
 51 #include "G4LossTableManager.hh"                   84 #include "G4LossTableManager.hh"
 52                                                    85 
 53 #include "G4ProductionCutsTable.hh"                86 #include "G4ProductionCutsTable.hh"
 54 #include "G4MaterialCutsCouple.hh"                 87 #include "G4MaterialCutsCouple.hh"
 55 #include "G4EmParameters.hh"                   << 
 56                                                    88 
 57 //....oooOO0OOooo........oooOO0OOooo........oo     89 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 58                                                    90 
 59 using namespace std;                               91 using namespace std;
 60                                                    92  
 61 G4eBremsstrahlung::G4eBremsstrahlung(const G4S     93 G4eBremsstrahlung::G4eBremsstrahlung(const G4String& name):
 62   G4VEnergyLossProcess(name)                   <<  94   G4VEnergyLossProcess(name), 
                                                   >>  95   isInitialised(false)
 63 {                                                  96 {
 64   SetProcessSubType(fBremsstrahlung);              97   SetProcessSubType(fBremsstrahlung);
 65   SetSecondaryParticle(G4Gamma::Gamma());      << 
 66   SetIonisation(false);                        << 
 67   SetCrossSectionType(fEmTwoPeaks);            << 
 68 }                                                  98 }
 69                                                    99 
 70 //....oooOO0OOooo........oooOO0OOooo........oo    100 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 71                                                   101 
 72 G4eBremsstrahlung::~G4eBremsstrahlung() = defa << 102 G4eBremsstrahlung::~G4eBremsstrahlung()
                                                   >> 103 {}
 73                                                   104 
 74 //....oooOO0OOooo........oooOO0OOooo........oo    105 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 75                                                   106 
 76 G4bool G4eBremsstrahlung::IsApplicable(const G    107 G4bool G4eBremsstrahlung::IsApplicable(const G4ParticleDefinition& p)
 77 {                                                 108 {
 78   return (&p == G4Electron::Electron() || &p =    109   return (&p == G4Electron::Electron() || &p == G4Positron::Positron());
 79 }                                                 110 }
 80                                                   111 
 81 //....oooOO0OOooo........oooOO0OOooo........oo    112 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 82                                                   113 
 83 void                                              114 void 
 84 G4eBremsstrahlung::InitialiseEnergyLossProcess    115 G4eBremsstrahlung::InitialiseEnergyLossProcess(const G4ParticleDefinition*,
 85                  const G4ParticleDefinition*)     116                  const G4ParticleDefinition*)
 86 {                                                 117 {
 87   if(!isInitialised) {                            118   if(!isInitialised) {
 88     G4EmParameters* param = G4EmParameters::In << 119     SetSecondaryParticle(G4Gamma::Gamma());
                                                   >> 120     SetIonisation(false);
 89                                                   121 
 90     G4double emax = param->MaxKinEnergy();     << 122     //    if (!EmModel(1)) { SetEmModel(new G4eBremsstrahlungModel(), 1); }
 91     G4VEmFluctuationModel* fm = nullptr;       << 123     if (!EmModel(1)) { SetEmModel(new G4SeltzerBergerModel(), 1); }
 92                                                << 124     if (!EmModel(2)) { SetEmModel(new G4eBremsstrahlungRelModel(), 2); }
 93     if (nullptr == EmModel(0)) { SetEmModel(ne << 125 
 94     G4double energyLimit = std::min(EmModel(0) << 126     G4double energyLimit = 1*GeV;
 95     EmModel(0)->SetHighEnergyLimit(energyLimit << 127 
 96     EmModel(0)->SetSecondaryThreshold(param->B << 128     EmModel(1)->SetLowEnergyLimit(MinKinEnergy());
 97     AddEmModel(1, EmModel(0), fm);             << 129     EmModel(1)->SetHighEnergyLimit(energyLimit);
 98                                                << 130     EmModel(2)->SetLowEnergyLimit(energyLimit);
 99     if(emax > energyLimit) {                   << 131     EmModel(2)->SetHighEnergyLimit(MaxKinEnergy());
100       if (nullptr == EmModel(1)) {             << 132                 
101   SetEmModel(new G4eBremsstrahlungRelModel()); << 133     G4VEmFluctuationModel* fm = 0;
102       }                                        << 134     AddEmModel(1, EmModel(1), fm);
103       EmModel(1)->SetLowEnergyLimit(energyLimi << 135     AddEmModel(2, EmModel(2), fm);
104       EmModel(1)->SetHighEnergyLimit(emax);    << 
105       EmModel(1)->SetSecondaryThreshold(param- << 
106       AddEmModel(1, EmModel(1), fm);           << 
107     }                                          << 
108     isInitialised = true;                         136     isInitialised = true;
109   }                                               137   }
                                                   >> 138   G4LossTableManager* man = G4LossTableManager::Instance(); 
                                                   >> 139   G4double eth = man->BremsstrahlungTh(); 
                                                   >> 140   EmModel(1)->SetSecondaryThreshold(eth);
                                                   >> 141   EmModel(2)->SetSecondaryThreshold(eth);
                                                   >> 142 
                                                   >> 143   // Only high energy model LMP flag is ON/OFF
                                                   >> 144   EmModel(1)->SetLPMFlag(false);
                                                   >> 145   EmModel(2)->SetLPMFlag(man->LPMFlag());
110 }                                                 146 }
111                                                   147 
112 //....oooOO0OOooo........oooOO0OOooo........oo    148 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
113                                                   149 
114 void G4eBremsstrahlung::StreamProcessInfo(std: << 150 void G4eBremsstrahlung::PrintInfo()
115 {                                                 151 {
116   if(nullptr != EmModel(0)) {                  << 152   if(EmModel(1)) {
117     G4EmParameters* param = G4EmParameters::In << 153     G4LossTableManager* man = G4LossTableManager::Instance(); 
118     G4double eth = param->BremsstrahlungTh();  << 154     G4double eth = man->BremsstrahlungTh(); 
119     out << "      LPM flag: " << param->LPM()  << 155     G4cout << "      LPM flag: " << man->LPMFlag() << " for E > " 
120   << EmModel(0)->HighEnergyLimit()/GeV << " Ge << 156      << EmModel(1)->HighEnergyLimit()/GeV << " GeV";
121     if(eth < DBL_MAX) {                           157     if(eth < DBL_MAX) { 
122       out << ",  VertexHighEnergyTh(GeV)= " << << 158       G4cout << ",  HighEnergyThreshold(GeV)= " << eth/GeV; 
123     }                                             159     }
124     out << G4endl;                             << 160     G4cout << G4endl;
125   }                                               161   }
126 }                                              << 
127                                                << 
128 //....oooOO0OOooo........oooOO0OOooo........oo << 
129                                                << 
130 void G4eBremsstrahlung::ProcessDescription(std << 
131 {                                              << 
132   out << "  Bremsstrahlung";                   << 
133   G4VEnergyLossProcess::ProcessDescription(out << 
134 }                                                 162 }
135                                                   163 
136 //....oooOO0OOooo........oooOO0OOooo........oo    164 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 
137                                                   165