Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/lowenergy/src/G4teoCrossSection.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 ]

  1 //
  2 // ********************************************************************
  3 // * License and Disclaimer                                           *
  4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.                             *
 10 // *                                                                  *
 11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                                                  *
 18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // ********************************************************************
 25 //
 26 //
 27 //         
 28 //
 29 // History:
 30 // -----------
 31 //  21 Apr 2009   ALF   1st implementation
 32 //  29 Apr 2009   ALF Updated Desing for Integration
 33 //  15 Mar 2011   ALF introduced the usage of G4AtomicShellEnumerator
 34 //  20 Oct 2011   ALF updated to take into account ECPSSR Form Factor
 35 //  09 Mar 2012   LP  update methods
 36 //  09 Mar 2012   ALF  update for M-shells Simulation
 37 //  10 Nov 2021   S. Guatelli & S. Bakr Added ECPSSR form factor documented 
 38 //                in Bakr et al, NIM B, vol. 436, pp:285-291, 2018 and
 39 //                called here ECPSSR_ANSTO
 40 
 41 #include "globals.hh"
 42 #include "G4teoCrossSection.hh"
 43 #include "G4Proton.hh"
 44 #include "G4ecpssrBaseKxsModel.hh"
 45 #include "G4ecpssrBaseLixsModel.hh"
 46 #include "G4ecpssrFormFactorKxsModel.hh"
 47 #include "G4ecpssrFormFactorLixsModel.hh"
 48 #include "G4ecpssrFormFactorMixsModel.hh"
 49 #include "G4ANSTOecpssrKxsModel.hh"
 50 #include "G4ANSTOecpssrLixsModel.hh"
 51 #include "G4ANSTOecpssrMixsModel.hh"
 52 
 53 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 54 G4teoCrossSection::G4teoCrossSection(const G4String& nam)
 55   :G4VhShellCrossSection(nam),totalCS(0.0)
 56 { 
 57   ecpssrShellMi = nullptr;
 58   if (nam == "ECPSSR_Analytical") 
 59     {
 60       ecpssrShellK  = new G4ecpssrBaseKxsModel();  
 61       ecpssrShellLi = new G4ecpssrBaseLixsModel();      
 62     }
 63   else if (nam == "ECPSSR_FormFactor")
 64     {
 65       ecpssrShellK  = new G4ecpssrFormFactorKxsModel();  
 66       ecpssrShellLi = new G4ecpssrFormFactorLixsModel();
 67       ecpssrShellMi = new G4ecpssrFormFactorMixsModel();
 68     }
 69   else if (nam == "ECPSSR_ANSTO")
 70     {
 71       ecpssrShellK  = new G4ANSTOecpssrKxsModel();  
 72       ecpssrShellLi = new G4ANSTOecpssrLixsModel();
 73       ecpssrShellMi = new G4ANSTOecpssrMixsModel();
 74     }   
 75   else 
 76     { 
 77       G4cout << "G4teoCrossSection::G4teoCrossSection: ERROR " 
 78        << " in cross section name ECPSSR_Analytical is used"
 79        << G4endl;
 80       ecpssrShellK  = new G4ecpssrBaseKxsModel();  
 81       ecpssrShellLi = new G4ecpssrBaseLixsModel();      
 82     }
 83 }
 84 
 85 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 86 
 87 G4teoCrossSection::~G4teoCrossSection()
 88 { 
 89   delete ecpssrShellK;
 90   delete ecpssrShellLi;
 91   delete ecpssrShellMi;
 92 }
 93 
 94 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 95 
 96 std::vector<G4double> G4teoCrossSection::GetCrossSection(G4int Z,
 97                G4double incidentEnergy,
 98                G4double mass,
 99                G4double,
100                const G4Material*)
101 {
102   std::vector<G4double> crossSections;
103 
104   crossSections.push_back( ecpssrShellK->CalculateCrossSection(Z, mass, incidentEnergy) );
105   
106   crossSections.push_back( ecpssrShellLi->CalculateL1CrossSection(Z, mass, incidentEnergy) );
107   crossSections.push_back( ecpssrShellLi->CalculateL2CrossSection(Z, mass, incidentEnergy) );
108   crossSections.push_back( ecpssrShellLi->CalculateL3CrossSection(Z, mass, incidentEnergy) );
109 
110   if (ecpssrShellMi) {
111     crossSections.push_back( ecpssrShellMi->CalculateM1CrossSection(Z, mass, incidentEnergy) );
112     crossSections.push_back( ecpssrShellMi->CalculateM2CrossSection(Z, mass, incidentEnergy) );
113     crossSections.push_back( ecpssrShellMi->CalculateM3CrossSection(Z, mass, incidentEnergy) );
114     crossSections.push_back( ecpssrShellMi->CalculateM4CrossSection(Z, mass, incidentEnergy) );
115     crossSections.push_back( ecpssrShellMi->CalculateM5CrossSection(Z, mass, incidentEnergy) );
116   }
117   return crossSections;
118 }
119 
120 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
121 
122 G4double G4teoCrossSection::CrossSection(G4int Z, G4AtomicShellEnumerator shell,
123            G4double incidentEnergy,
124            G4double mass,
125            const G4Material*)
126 {
127   G4double res = 0.0;
128   if(shell > 3 && !ecpssrShellMi) {
129     return res; 
130   } 
131   else if(shell > 8) { 
132     return res;
133   } 
134   else if(fKShell  == shell) 
135     { 
136       res = ecpssrShellK->CalculateCrossSection(Z, mass, incidentEnergy);
137     } 
138   else if(fL1Shell == shell) 
139     { 
140       res = ecpssrShellLi->CalculateL1CrossSection(Z, mass, incidentEnergy);
141     }   
142   else if(fL2Shell == shell) 
143     { 
144       res = ecpssrShellLi->CalculateL2CrossSection(Z, mass, incidentEnergy);
145     }   
146   else if(fL3Shell == shell) 
147     { 
148       res = ecpssrShellLi->CalculateL3CrossSection(Z, mass, incidentEnergy);
149     }
150   else if(fM1Shell == shell) 
151     { 
152       res = ecpssrShellMi->CalculateM1CrossSection(Z, mass, incidentEnergy);
153     }
154   else if(fM2Shell == shell) 
155     { 
156       res = ecpssrShellMi->CalculateM2CrossSection(Z, mass, incidentEnergy);
157     }
158   else if(fM3Shell == shell) 
159     { 
160       res = ecpssrShellMi->CalculateM3CrossSection(Z, mass, incidentEnergy);
161     }
162   else if(fM4Shell == shell) 
163     { 
164       res = ecpssrShellMi->CalculateM4CrossSection(Z, mass, incidentEnergy);
165     }
166   else if(fM5Shell == shell) 
167     { 
168       res = ecpssrShellMi->CalculateM5CrossSection(Z, mass, incidentEnergy);
169     }
170   return res;
171 }
172 
173 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
174 
175 std::vector<G4double> G4teoCrossSection::Probabilities(G4int Z,
176                    G4double incidentEnergy,
177                    G4double mass,
178                    G4double deltaEnergy,
179                    const G4Material*) 
180 {
181   std::vector<G4double> crossSections = 
182     GetCrossSection(Z, incidentEnergy, mass, deltaEnergy);
183 
184   for (size_t i=0; i<crossSections.size(); ++i ) {    
185     if (totalCS) {
186       crossSections[i] = crossSections[i]/totalCS;
187     }    
188   }
189   return crossSections;
190 }
191 
192 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
193 
194 void G4teoCrossSection::SetTotalCS(G4double val){
195   totalCS = val;
196 }
197 
198 
199 
200