Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/lowenergy/src/G4eIonisationCrossSectionHandler.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/lowenergy/src/G4eIonisationCrossSectionHandler.cc (Version 11.3.0) and /processes/electromagnetic/lowenergy/src/G4eIonisationCrossSectionHandler.cc (Version 4.0)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                    <<   3 // * DISCLAIMER                                                       *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th <<   5 // * The following disclaimer summarizes all the specific disclaimers *
  6 // * the Geant4 Collaboration.  It is provided <<   6 // * of contributors to this software. The specific disclaimers,which *
  7 // * conditions of the Geant4 Software License <<   7 // * govern, are listed with their locations in:                      *
  8 // * LICENSE and available at  http://cern.ch/ <<   8 // *   http://cern.ch/geant4/license                                  *
  9 // * include a list of copyright holders.      << 
 10 // *                                                9 // *                                                                  *
 11 // * Neither the authors of this software syst     10 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     11 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     12 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     13 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file  <<  14 // * use.                                                             *
 16 // * for the full disclaimer and the limitatio << 
 17 // *                                               15 // *                                                                  *
 18 // * This  code  implementation is the result  <<  16 // * This  code  implementation is the  intellectual property  of the *
 19 // * technical work of the GEANT4 collaboratio <<  17 // * GEANT4 collaboration.                                            *
 20 // * By using,  copying,  modifying or  distri <<  18 // * By copying,  distributing  or modifying the Program (or any work *
 21 // * any work based  on the software)  you  ag <<  19 // * based  on  the Program)  you indicate  your  acceptance of  this *
 22 // * use  in  resulting  scientific  publicati <<  20 // * statement, and all its terms.                                    *
 23 // * acceptance of all terms of the Geant4 Sof << 
 24 // *******************************************     21 // ********************************************************************
 25 //                                                 22 //
                                                   >>  23 // $Id: G4eIonisationCrossSectionHandler.cc,v 1.6 2001/11/29 19:01:37 vnivanch Exp $
                                                   >>  24 // GEANT4 tag $Name: geant4-04-00 $
 26 //                                                 25 //
 27 // -------------------------------------------     26 // -------------------------------------------------------------------
 28 //                                                 27 //
 29 // GEANT4 Class file                               28 // GEANT4 Class file
 30 //                                                 29 //
 31 //                                                 30 //
 32 // File name:     G4eIonisationCrossSectionHan     31 // File name:     G4eIonisationCrossSectionHandler
 33 //                                                 32 //
 34 // Author:        V.Ivanchenko (Vladimir.Ivanc     33 // Author:        V.Ivanchenko (Vladimir.Ivanchenko@cern.ch)
 35 //                                                 34 // 
 36 // Creation date: 25 Sept 2001                     35 // Creation date: 25 Sept 2001
 37 //                                                 36 //
 38 // Modifications:                                  37 // Modifications: 
 39 // 10 Oct 2001  M.G. Pia     Revision to impro <<  38 // 10 Oct 2001  M.G. Pia        Revision to improve code quality and consistency with design
 40 // 19 Jul 2002   VI          Create composite  << 
 41 // 21 Jan 2003  V.Ivanchenko Cut per region    << 
 42 // 28 Jan 2009  L.Pandola    Added public meth << 
 43 //                           G4LowEnergyIonisa << 
 44 // 15 Jul 2009   Nicolas A. Karakatsanis       << 
 45 //                                             << 
 46 //                           - BuildCrossSecti << 
 47 //                             logarithmic val << 
 48 //                             It retrieves th << 
 49 //                             respective log  << 
 50 //                             The EM data set << 
 51 //                             These initializ << 
 52 //                             operations      << 
 53 //                                             << 
 54 //                                             << 
 55 //                                                 39 //
 56 // -------------------------------------------     40 // -------------------------------------------------------------------
 57                                                    41 
 58 #include "G4eIonisationCrossSectionHandler.hh"     42 #include "G4eIonisationCrossSectionHandler.hh"
 59 #include "G4SystemOfUnits.hh"                  << 
 60 #include "G4VEnergySpectrum.hh"                    43 #include "G4VEnergySpectrum.hh"
 61 #include "G4DataVector.hh"                         44 #include "G4DataVector.hh"
 62 #include "G4CompositeEMDataSet.hh"                 45 #include "G4CompositeEMDataSet.hh"
 63 #include "G4VDataSetAlgorithm.hh"                  46 #include "G4VDataSetAlgorithm.hh"
 64 #include "G4LinLogLogInterpolation.hh"         << 
 65 #include "G4SemiLogInterpolation.hh"               47 #include "G4SemiLogInterpolation.hh"
 66 #include "G4VEMDataSet.hh"                         48 #include "G4VEMDataSet.hh"
 67 #include "G4EMDataSet.hh"                          49 #include "G4EMDataSet.hh"
 68 #include "G4Material.hh"                           50 #include "G4Material.hh"
 69 #include "G4ProductionCutsTable.hh"            <<  51 #include "G4MaterialTable.hh"
                                                   >>  52 
 70                                                    53 
 71 //....oooOO0OOooo........oooOO0OOooo........oo <<  54 G4eIonisationCrossSectionHandler::G4eIonisationCrossSectionHandler(
 72  G4eIonisationCrossSectionHandler::G4eIonisati <<  55     const G4VEnergySpectrum* spec, G4VDataSetAlgorithm* alg, 
 73     const G4VEnergySpectrum* spec, G4VDataSetA << 
 74     G4double emin, G4double emax, G4int nbin)      56     G4double emin, G4double emax, G4int nbin)
 75  :  G4VCrossSectionHandler(),                      57  :  G4VCrossSectionHandler(),
 76     theParam(spec),verbose(0)                  <<  58     theParam(spec)
 77 {                                                  59 {
 78   G4VCrossSectionHandler::Initialise(alg, emin     60   G4VCrossSectionHandler::Initialise(alg, emin, emax, nbin);
 79   interp = new G4LinLogLogInterpolation();     <<  61   interp = new G4SemiLogInterpolation();
 80 }                                                  62 }
 81                                                    63 
 82 //....oooOO0OOooo........oooOO0OOooo........oo << 
 83                                                    64 
 84 G4eIonisationCrossSectionHandler::~G4eIonisati <<  65 G4eIonisationCrossSectionHandler::~G4eIonisationCrossSectionHandler() 
 85 {                                                  66 {
 86   delete interp;                                   67   delete interp;
 87 }                                                  68 }
 88                                                    69 
 89 //....oooOO0OOooo........oooOO0OOooo........oo << 
 90                                                    70 
 91 std::vector<G4VEMDataSet*>* G4eIonisationCross <<  71 G4std::vector<G4VEMDataSet*>* G4eIonisationCrossSectionHandler::BuildCrossSectionsForMaterials(
 92                         const G4DataVector& en <<  72                              const G4DataVector& energyVector, 
 93                         const G4DataVector* en <<  73                              const G4DataVector* energyCuts)
 94 {                                                  74 {
 95   std::vector<G4VEMDataSet*>* set = new std::v <<  75 
                                                   >>  76   G4std::vector<G4VEMDataSet*>* set = new G4std::vector<G4VEMDataSet*>;
 96                                                    77 
 97   G4DataVector* energies;                          78   G4DataVector* energies;
 98   G4DataVector* cs;                                79   G4DataVector* cs;
                                                   >>  80   G4int nOfBins = energyVector.size();
 99                                                    81 
100   G4DataVector* log_energies;                  <<  82   const G4MaterialTable* materialTable = G4Material::GetMaterialTable();
101   G4DataVector* log_cs;                        <<  83   if (materialTable == 0)
102                                                <<  84     G4Exception("G4VCrossSectionHandler::G4VCrossSectionHandler - no MaterialTable found)");
103   std::size_t nOfBins = energyVector.size();   << 
104                                                    85 
105   const G4ProductionCutsTable* theCoupleTable= <<  86   G4int nMaterials = G4Material::GetNumberOfMaterials();
106         G4ProductionCutsTable::GetProductionCu << 
107   G4int numOfCouples = (G4int)theCoupleTable-> << 
108                                                    87 
109   for (G4int mLocal=0; mLocal<numOfCouples; ++ <<  88   for (G4int m=0; m<nMaterials; m++) {
110                                                    89 
111     const G4MaterialCutsCouple* couple = theCo <<  90     const G4Material* material = (*materialTable)[m];
112     const G4Material* material= couple->GetMat << 
113     const G4ElementVector* elementVector = mat     91     const G4ElementVector* elementVector = material->GetElementVector();
                                                   >>  92     //const G4double* nAtomsPerVolume = material->GetVecNbOfAtomsPerVolume();
114     const G4double* nAtomsPerVolume = material     93     const G4double* nAtomsPerVolume = material->GetAtomicNumDensityVector();
115     G4int nElements = (G4int)material->GetNumb <<  94     G4int nElements = material->GetNumberOfElements();
116                                                << 
117     if(verbose > 0)                            << 
118       {                                        << 
119   G4cout << "eIonisation CS for " << mLocal << << 
120          << material->GetName()                << 
121          << "  eEl= " << nElements << G4endl;  << 
122       }                                        << 
123                                                << 
124     G4double tcut  = (*energyCuts)[mLocal];    << 
125                                                    95 
                                                   >>  96     G4double tcut  = (*energyCuts)[m];
                                                   >>  97     
126     G4VDataSetAlgorithm* algo = interp->Clone(     98     G4VDataSetAlgorithm* algo = interp->Clone();
127     G4VEMDataSet* setForMat = new G4CompositeE     99     G4VEMDataSet* setForMat = new G4CompositeEMDataSet(algo,1.,1.);
128                                                   100 
129     for (G4int i=0; i<nElements; ++i) {        << 101     for (G4int i=0; i<nElements; i++) {
130                                                << 102  
131       G4int Z = (G4int) (*elementVector)[i]->G    103       G4int Z = (G4int) (*elementVector)[i]->GetZ();
132       G4int nShells = NumberOfComponents(Z);      104       G4int nShells = NumberOfComponents(Z);
133                                                << 
134       energies = new G4DataVector;                105       energies = new G4DataVector;
135       cs       = new G4DataVector;                106       cs       = new G4DataVector;
136                                                << 
137       log_energies = new G4DataVector;         << 
138       log_cs       = new G4DataVector;         << 
139                                                << 
140       G4double density = nAtomsPerVolume[i];      107       G4double density = nAtomsPerVolume[i];
141                                                   108 
142       for (std::size_t bin=0; bin<nOfBins; ++b << 109       for (G4int bin=0; bin<nOfBins; bin++) {
143                                                   110 
144         G4double e = energyVector[bin];           111         G4double e = energyVector[bin];
145         energies->push_back(e);                   112         energies->push_back(e);
146         log_energies->push_back(std::log10(e)) << 
147         G4double value = 0.0;                     113         G4double value = 0.0;
148         G4double log_value = -300;             << 
149                                                   114 
150         if(e > tcut) {                            115         if(e > tcut) {
151           for (G4int n=0; n<nShells; n++) {       116           for (G4int n=0; n<nShells; n++) {
152             G4double cross = FindValue(Z, e, n    117             G4double cross = FindValue(Z, e, n);
153             G4double p = theParam->Probability << 118             G4double p = theParam->Probability(Z, tcut, e, e, n);  
154             value += cross * p * density;         119             value += cross * p * density;
155                                                << 120       /*
156       if(verbose>0 && mLocal == 0 && e>=1. &&  << 121             G4cout << "G4eIonisationCrossSectionHandler: e= " << e
157       {                                        << 122                    << " n= " << n
158               G4cout << "G4eIonCrossSH: e(MeV) << 123                    << " cross= " << cross
159                      << " n= " << n            << 124                    << " p= " << p
160                      << " cross= " << cross    << 125                    << " value= " << value
161                      << " p= " << p            << 126                    << G4endl;
162                      << " value= " << value    << 127       */
163                      << " tcut(MeV)= " << tcut << 
164                      << " rho= " << density    << 
165                      << " Z= " << Z            << 
166                      << G4endl;                << 
167         }                                      << 
168     }                                             128     }
169           if (value == 0.) value = 1e-300;     << 
170           log_value = std::log10(value);       << 
171   }                                               129   }
172         cs->push_back(value);                     130         cs->push_back(value);
173         log_cs->push_back(log_value);          << 131       } 
174       }                                        << 132       G4VDataSetAlgorithm* algo = interp->Clone();
175       G4VDataSetAlgorithm* algoLocal = interp- << 133       G4VEMDataSet* elSet = new G4EMDataSet(i,energies,cs,algo,1.,1.);
176       G4VEMDataSet* elSet = new G4EMDataSet(i, << 
177                                                << 
178       setForMat->AddComponent(elSet);             134       setForMat->AddComponent(elSet);
179     }                                             135     }
180     set->push_back(setForMat);                    136     set->push_back(setForMat);
181   }                                               137   }
182                                                << 138  
183   return set;                                     139   return set;
184 }                                              << 
185                                                << 
186 //....oooOO0OOooo........oooOO0OOooo........oo << 
187                                                << 
188 G4double G4eIonisationCrossSectionHandler::Get << 
189                        G4double cutEnergy,     << 
190                        G4int Z)                << 
191 {                                              << 
192   G4int nShells = NumberOfComponents(Z);       << 
193   G4double value = 0.;                         << 
194   if(energy > cutEnergy)                       << 
195     {                                          << 
196       for (G4int n=0; n<nShells; ++n) {        << 
197   G4double cross = FindValue(Z, energy, n);    << 
198   G4double p = theParam->Probability(Z, cutEne << 
199   value += cross * p;                          << 
200       }                                        << 
201     }                                          << 
202   return value;                                << 
203 }                                                 140 }
204                                                   141