Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/utils/src/G4EmSaturation.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/utils/src/G4EmSaturation.cc (Version 11.3.0) and /processes/electromagnetic/utils/src/G4EmSaturation.cc (Version 9.3.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: G4EmSaturation.cc,v 1.10 2009/09/25 09:16:40 vnivanch Exp $
                                                   >>  27 // GEANT4 tag $Name: geant4-09-03-patch-01 $
 26 //                                                 28 //
 27 // -------------------------------------------     29 // -------------------------------------------------------------------
 28 //                                                 30 //
 29 // GEANT4 Class file                               31 // GEANT4 Class file
 30 //                                                 32 //
 31 //                                                 33 //
 32 // File name:     G4EmSaturation                   34 // File name:     G4EmSaturation
 33 //                                                 35 //
 34 // Author:        Vladimir Ivanchenko              36 // Author:        Vladimir Ivanchenko
 35 //                                                 37 //
 36 // Creation date: 18.02.2008                       38 // Creation date: 18.02.2008
 37 //                                                 39 //
 38 // Modifications:                                  40 // Modifications:
 39 //                                                 41 //
 40 // -------------------------------------------     42 // -------------------------------------------------------------
 41                                                    43 
 42 //....oooOO0OOooo........oooOO0OOooo........oo     44 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 43 //....oooOO0OOooo........oooOO0OOooo........oo     45 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 44                                                    46 
 45 #include "G4EmSaturation.hh"                       47 #include "G4EmSaturation.hh"
 46 #include "G4PhysicalConstants.hh"              << 
 47 #include "G4SystemOfUnits.hh"                  << 
 48 #include "G4LossTableManager.hh"                   48 #include "G4LossTableManager.hh"
 49 #include "G4NistManager.hh"                        49 #include "G4NistManager.hh"
 50 #include "G4Material.hh"                           50 #include "G4Material.hh"
 51 #include "G4MaterialCutsCouple.hh"                 51 #include "G4MaterialCutsCouple.hh"
 52 #include "G4ParticleTable.hh"                  <<  52 #include "G4Electron.hh"
                                                   >>  53 #include "G4Proton.hh"
 53                                                    54 
 54 //....oooOO0OOooo........oooOO0OOooo........oo     55 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 55                                                    56 
 56 std::size_t G4EmSaturation::nMaterials = 0;    <<  57 G4EmSaturation::G4EmSaturation()
 57 std::vector<G4double> G4EmSaturation::massFact << 
 58 std::vector<G4double> G4EmSaturation::effCharg << 
 59 std::vector<G4double> G4EmSaturation::g4MatDat << 
 60 std::vector<G4String> G4EmSaturation::g4MatNam << 
 61                                                << 
 62 G4EmSaturation::G4EmSaturation(G4int verb)     << 
 63 {                                                  58 {
 64   verbose = verb;                              <<  59   verbose = 1;
 65   nWarnings = nG4Birks = 0;                    <<  60   manager = 0;
 66                                                <<  61   curMaterial = 0;
 67   electron = nullptr;                          <<  62   curBirks = 0.0;
 68   proton   = nullptr;                          <<  63   curRatio = 1.0;
 69   nist     = G4NistManager::Instance();        <<  64   curChargeSq = 1.0;
 70   InitialiseG4Saturation();                    <<  65   nMaterials = 0;
                                                   >>  66   electron = 0;
                                                   >>  67   Initialise(); 
 71 }                                                  68 }
 72                                                    69 
 73 //....oooOO0OOooo........oooOO0OOooo........oo     70 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 74                                                    71 
 75 G4EmSaturation::~G4EmSaturation() = default;   <<  72 G4EmSaturation::~G4EmSaturation()
                                                   >>  73 {}
 76                                                    74 
 77 //....oooOO0OOooo........oooOO0OOooo........oo     75 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 78                                                    76 
 79 G4double G4EmSaturation::VisibleEnergyDepositi     77 G4double G4EmSaturation::VisibleEnergyDeposition(
 80                                       const G4     78                                       const G4ParticleDefinition* p, 
 81                                       const G4 <<  79               const G4MaterialCutsCouple* couple, 
 82                                       G4double <<  80               G4double length,
 83                                       G4double <<  81               G4double edep,
 84                                       G4double <<  82               G4double niel)
 85 {                                                  83 {
 86   // no energy deposition                      <<  84   if(edep <= 0.0) return 0.0;
 87   if(edep <= 0.0) { return 0.0; }              << 
 88                                                << 
 89   // zero step length may happens only if step << 
 90   // is applied, in that case saturation shoul << 
 91   if(length <= 0.0) { return edep; }           << 
 92                                                    85 
 93   G4double evis = edep;                            86   G4double evis = edep;
 94   G4double bfactor = couple->GetMaterial()->Ge <<  87   G4double bfactor = FindBirksCoefficient(couple->GetMaterial());
 95                                                    88 
 96   if(bfactor > 0.0) {                              89   if(bfactor > 0.0) { 
 97                                                    90 
                                                   >>  91     G4int pdgCode = p->GetPDGEncoding();
 98     // atomic relaxations for gamma incident       92     // atomic relaxations for gamma incident
 99     if(22 ==  p->GetPDGEncoding()) {           <<  93     if(22 == pdgCode) {
100       //G4cout << "%% gamma edep= " << edep/ke <<  94       evis /= (1.0 + bfactor*edep/manager->GetRange(electron,edep,couple));
101       evis /= (1.0 + bfactor*edep/             << 
102         G4LossTableManager::Instance()->GetRan << 
103                                                    95 
104       // energy loss                               96       // energy loss
105     } else {                                       97     } else {
106                                                    98 
107       // protections                               99       // protections
108       G4double nloss = std::max(niel, 0.0);    << 100       G4double nloss = niel;
                                                   >> 101       if(nloss < 0.0) nloss = 0.0;
109       G4double eloss = edep - nloss;              102       G4double eloss = edep - nloss;
110                                                   103 
111       // neutrons and neutral hadrons          << 104       // neutrons
112       if(0.0 == p->GetPDGCharge() || eloss < 0 << 105       if(2112 == pdgCode || eloss < 0.0 || length <= 0.0) {
113         nloss = edep;                          << 106   nloss = edep;
114         eloss = 0.0;                              107         eloss = 0.0;
115       } else {                                 << 
116                                                << 
117   // continues energy loss                     << 
118   eloss /= (1.0 + bfactor*eloss/length);       << 
119       }                                           108       }
                                                   >> 109 
                                                   >> 110       // continues energy loss
                                                   >> 111       if(eloss > 0.0) eloss /= (1.0 + bfactor*eloss/length);
                                                   >> 112  
120       // non-ionizing energy loss                 113       // non-ionizing energy loss
121       if(nloss > 0.0) {                           114       if(nloss > 0.0) {
122         std::size_t idx = couple->GetMaterial( << 115         if(!proton) {proton = G4Proton::Proton();}
123         G4double escaled = nloss*massFactors[i << 116         G4double escaled = nloss*curRatio;
124   /*                                           << 117         G4double s = manager->GetRange(proton,escaled,couple)/curChargeSq; 
125         G4cout << "%% p edep= " << nloss/keV < << 118   nloss /= (1.0 + bfactor*nloss/s);
126                << escaled << " MeV  in " << co << 
127                << "  " << p->GetParticleName() << 
128                << G4endl;                      << 
129   G4cout << proton->GetParticleName() << G4end << 
130   */                                           << 
131         G4double range = G4LossTableManager::I << 
132           ->GetRange(proton,escaled,couple)/ef << 
133         nloss /= (1.0 + bfactor*nloss/range);  << 
134       }                                           119       }
                                                   >> 120 
135       evis = eloss + nloss;                       121       evis = eloss + nloss;
136     }                                             122     }
137   }                                               123   }
                                                   >> 124   
138   return evis;                                    125   return evis;
139 }                                                 126 }
140                                                   127 
141 //....oooOO0OOooo........oooOO0OOooo........oo    128 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
142                                                   129 
143 void G4EmSaturation::InitialiseG4Saturation()  << 
144 {                                              << 
145   if(nMaterials == G4Material::GetNumberOfMate << 
146   nMaterials = G4Material::GetNumberOfMaterial << 
147   massFactors.resize(nMaterials, 1.0);         << 
148   effCharges.resize(nMaterials, 1.0);          << 
149                                                << 
150   if(0 == nG4Birks) {  InitialiseG4materials() << 
151                                                << 
152   for(std::size_t i=0; i<nMaterials; ++i) {    << 
153     InitialiseBirksCoefficient((*G4Material::G << 
154   }                                            << 
155   if(verbose > 0) { DumpBirksCoefficients(); } << 
156 }                                              << 
157                                                << 
158 //....oooOO0OOooo........oooOO0OOooo........oo << 
159                                                << 
160 G4double G4EmSaturation::FindG4BirksCoefficien    130 G4double G4EmSaturation::FindG4BirksCoefficient(const G4Material* mat)
161 {                                                 131 {
162   if(0 == nG4Birks) {  InitialiseG4materials() << 
163                                                << 
164   G4String name = mat->GetName();                 132   G4String name = mat->GetName();
165   // is this material in the vector?              133   // is this material in the vector?
166                                                   134   
167   for(G4int j=0; j<nG4Birks; ++j) {            << 135   for(G4int j=0; j<nG4Birks; j++) {
168     if(name == g4MatNames[j]) {                   136     if(name == g4MatNames[j]) {
169       if(verbose > 0)                             137       if(verbose > 0) 
170         G4cout << "### G4EmSaturation::FindG4B << 138   G4cout << "### G4EmSaturation::FindG4BirksCoefficient for "
171                << name << " is " << g4MatData[ << 139          << name << " is " << g4MatData[j]*MeV/mm << " mm/MeV "
172                << G4endl;                      << 140          << G4endl;
173       return g4MatData[j];                        141       return g4MatData[j];
174     }                                             142     }
175   }                                               143   }
176   return 0.0;                                  << 144   return FindBirksCoefficient(mat);
177 }                                                 145 }
178                                                   146 
179 //....oooOO0OOooo........oooOO0OOooo........oo    147 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
180                                                   148 
181 void G4EmSaturation::InitialiseBirksCoefficien << 149 G4double G4EmSaturation::FindBirksCoefficient(const G4Material* mat)
182 {                                                 150 {
183   // electron and proton should exist in any c << 151   // electron should exist in any case
184   if(nullptr == electron) {                    << 152   if(!manager) {
185     electron = G4ParticleTable::GetParticleTab << 153     manager = G4LossTableManager::Instance();
186     proton = G4ParticleTable::GetParticleTable << 154     nist    = G4NistManager::Instance();
187     if(nullptr == electron) {                  << 155     electron= G4Electron::Electron();
188       G4Exception("G4EmSaturation::InitialiseB << 156     proton  = 0;
189       FatalException, "electron should exist") << 
190     }                                          << 
191   }                                               157   }
192                                                   158 
193   G4double curBirks = mat->GetIonisation()->Ge << 159   if(mat == curMaterial) return curBirks;
                                                   >> 160 
                                                   >> 161   curMaterial = mat;
                                                   >> 162   curBirks = 0.0;
                                                   >> 163   curRatio = 1.0;
                                                   >> 164   curChargeSq = 1.0;
                                                   >> 165 
                                                   >> 166   // seach in the run-time list
                                                   >> 167   for(G4int i=0; i<nMaterials; i++) {
                                                   >> 168     if(mat == matPointers[i]) {
                                                   >> 169       curBirks = mat->GetIonisation()->GetBirksConstant();
                                                   >> 170       curRatio = massFactors[i];
                                                   >> 171       curChargeSq = effCharges[i];
                                                   >> 172       return curBirks;
                                                   >> 173     }
                                                   >> 174   }
194                                                   175 
195   G4String name = mat->GetName();                 176   G4String name = mat->GetName();
                                                   >> 177   curBirks = mat->GetIonisation()->GetBirksConstant();
196                                                   178 
197   // material has no Birks coeffitient defined    179   // material has no Birks coeffitient defined
198   // seach in the Geant4 list                     180   // seach in the Geant4 list
199   if(curBirks == 0.0) {                           181   if(curBirks == 0.0) {
200     for(G4int j=0; j<nG4Birks; ++j) {          << 182     for(G4int j=0; j<nG4Birks; j++) {
201       if(name == g4MatNames[j]) {                 183       if(name == g4MatNames[j]) {
202         mat->GetIonisation()->SetBirksConstant << 184   mat->GetIonisation()->SetBirksConstant(g4MatData[j]);
203         curBirks = g4MatData[j];                  185         curBirks = g4MatData[j];
204         break;                                    186         break;
205       }                                           187       }
206     }                                             188     }
207   }                                               189   }
208                                                   190 
209   if(curBirks == 0.0) { return; }              << 191   if(curBirks == 0.0 && verbose > 0) {
                                                   >> 192       G4cout << "### G4EmSaturation::FindBirksCoefficient fails "
                                                   >> 193   " for material " << name << G4endl;
                                                   >> 194   }
210                                                   195 
211   // compute mean mass ratio                      196   // compute mean mass ratio
212   G4double curRatio = 0.0;                     << 197   curRatio = 0.0;
213   G4double curChargeSq = 0.0;                  << 198   curChargeSq = 0.0;
214   G4double norm = 0.0;                            199   G4double norm = 0.0;
215   const G4ElementVector* theElementVector = ma    200   const G4ElementVector* theElementVector = mat->GetElementVector();
216   const G4double* theAtomNumDensityVector = ma    201   const G4double* theAtomNumDensityVector = mat->GetVecNbOfAtomsPerVolume();
217   std::size_t nelm = mat->GetNumberOfElements( << 202   size_t nelm = mat->GetNumberOfElements();
218   for (std::size_t i=0; i<nelm; ++i) {         << 203   for (size_t i=0; i<nelm; i++) {
219     const G4Element* elm = (*theElementVector)    204     const G4Element* elm = (*theElementVector)[i];
220     G4int Z = elm->GetZasInt();                << 205     G4double Z = elm->GetZ();
221     G4double w = theAtomNumDensityVector[i];   << 206     G4double w = Z*Z*theAtomNumDensityVector[i];
222     curRatio += w/nist->GetAtomicMassAmu(Z);   << 207     curRatio += w/nist->GetAtomicMassAmu(G4int(Z));
223     curChargeSq += (Z*Z)*w;                    << 208     curChargeSq = Z*Z*w;
224     norm += w;                                    209     norm += w;
225   }                                               210   }
226   if ( norm > 0.0) { norm = 1.0/norm; }        << 211   curRatio *= proton_mass_c2/norm;
227   curRatio *= (CLHEP::proton_mass_c2*norm);    << 212   curChargeSq /= norm;
228   curChargeSq *= norm;                         << 
229                                                   213 
230   // store results                                214   // store results
231   std::size_t idx = mat->GetIndex();           << 215   matPointers.push_back(mat);
232   massFactors[idx] = curRatio;                 << 216   matNames.push_back(name);
233   effCharges[idx] = curChargeSq;               << 217   massFactors.push_back(curRatio);
                                                   >> 218   effCharges.push_back(curChargeSq);
                                                   >> 219   nMaterials++;
                                                   >> 220   if(curBirks > 0.0 && verbose > 0) {
                                                   >> 221     G4cout << "### G4EmSaturation::FindBirksCoefficient Birks coefficient for "
                                                   >> 222      << name << "  " << curBirks*MeV/mm << " mm/MeV" << G4endl;
                                                   >> 223   }
                                                   >> 224   return curBirks;
234 }                                                 225 }
235                                                   226 
236 //....oooOO0OOooo........oooOO0OOooo........oo    227 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
237                                                   228 
238 void G4EmSaturation::DumpBirksCoefficients()      229 void G4EmSaturation::DumpBirksCoefficients()
239 {                                                 230 {
240   G4cout << "### Birks coefficients used in ru << 231   if(nMaterials > 0) {
241   const G4MaterialTable* mtable = G4Material:: << 232     G4cout << "### Birks coeffitients used in run time" << G4endl;
242   for(std::size_t i=0; i<nMaterials; ++i) {    << 233     for(G4int i=0; i<nMaterials; i++) {
243     const G4Material* mat = (*mtable)[i];      << 234       G4double br = matPointers[i]->GetIonisation()->GetBirksConstant();
244     G4double br = mat->GetIonisation()->GetBir << 235       G4cout << "   " << matNames[i] << "     " 
245     if(br > 0.0) {                             << 
246       G4cout << "   " << mat->GetName() << "   << 
247        << br*MeV/mm << " mm/MeV" << "     "       236        << br*MeV/mm << " mm/MeV" << "     "
248        << br*mat->GetDensity()*MeV*cm2/g       << 237        << br*matPointers[i]->GetDensity()*MeV*cm2/g 
249        << " g/cm^2/MeV  massFactor=  " << mass << 238        << " g/cm^2/MeV" 
250        << " effCharge= " << effCharges[i] << G << 239        << G4endl;
251     }                                             240     }
252   }                                               241   }
253 }                                                 242 }
254                                                   243 
255 //....oooOO0OOooo........oooOO0OOooo........oo    244 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
256                                                   245 
257 void G4EmSaturation::DumpG4BirksCoefficients()    246 void G4EmSaturation::DumpG4BirksCoefficients()
258 {                                                 247 {
259   if(nG4Birks > 0) {                              248   if(nG4Birks > 0) {
260     G4cout << "### Birks coefficients for Gean << 249     G4cout << "### Birks coeffitients for Geant4 materials" << G4endl;
261     for(G4int i=0; i<nG4Birks; ++i) {          << 250     for(G4int i=0; i<nG4Birks; i++) {
262       G4cout << "   " << g4MatNames[i] << "       251       G4cout << "   " << g4MatNames[i] << "   " 
263              << g4MatData[i]*MeV/mm << " mm/Me << 252        << g4MatData[i]*MeV/mm << " mm/MeV" << G4endl;
264     }                                             253     }
265   }                                               254   }
266 }                                                 255 }
267                                                   256 
268 //....oooOO0OOooo........oooOO0OOooo........oo    257 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
269                                                   258 
270 void G4EmSaturation::InitialiseG4materials()   << 259 void G4EmSaturation::Initialise()
271 {                                                 260 {
272   nG4Birks = 4;                                << 
273   g4MatData.reserve(nG4Birks);                 << 
274                                                << 
275   // M.Hirschberg et al., IEEE Trans. Nuc. Sci    261   // M.Hirschberg et al., IEEE Trans. Nuc. Sci. 39 (1992) 511
276   // SCSN-38 kB = 0.00842 g/cm^2/MeV; rho = 1.    262   // SCSN-38 kB = 0.00842 g/cm^2/MeV; rho = 1.06 g/cm^3
277   g4MatNames.push_back("G4_POLYSTYRENE");         263   g4MatNames.push_back("G4_POLYSTYRENE");
278   g4MatData.push_back(0.07943*mm/MeV);            264   g4MatData.push_back(0.07943*mm/MeV);
279                                                   265 
280   // C.Fabjan (private communication)             266   // C.Fabjan (private communication)
281   // kB = 0.006 g/cm^2/MeV; rho = 7.13 g/cm^3     267   // kB = 0.006 g/cm^2/MeV; rho = 7.13 g/cm^3
282   g4MatNames.push_back("G4_BGO");                 268   g4MatNames.push_back("G4_BGO");
283   g4MatData.push_back(0.008415*mm/MeV);           269   g4MatData.push_back(0.008415*mm/MeV);
284                                                   270 
285   // A.Ribon analysis of publications             271   // A.Ribon analysis of publications
286   // Scallettar et al., Phys. Rev. A25 (1982)     272   // Scallettar et al., Phys. Rev. A25 (1982) 2419.
287   // NIM A 523 (2004) 275.                        273   // NIM A 523 (2004) 275. 
288   // kB = 0.022 g/cm^2/MeV; rho = 1.396 g/cm^3    274   // kB = 0.022 g/cm^2/MeV; rho = 1.396 g/cm^3; 
289   // ATLAS Efield = 10 kV/cm provide the stron << 275   // ATLAS Efield = 10 kV/cm provide the strongest effect 
290   // kB = 0.1576*mm/MeV                        << 
291   // A. Kiryunin and P.Strizenec "Geant4 hadro << 
292   // working group meeting " kB = 0.041/9.13 g << 
293   g4MatNames.push_back("G4_lAr");                 276   g4MatNames.push_back("G4_lAr");
294   g4MatData.push_back(0.032*mm/MeV);           << 277   g4MatData.push_back(0.1576*mm/MeV);
295                                                   278 
296   //G4_BARIUM_FLUORIDE                            279   //G4_BARIUM_FLUORIDE
297   //G4_CESIUM_IODIDE                              280   //G4_CESIUM_IODIDE
298   //G4_GEL_PHOTO_EMULSION                         281   //G4_GEL_PHOTO_EMULSION
299   //G4_PHOTO_EMULSION                             282   //G4_PHOTO_EMULSION
300   //G4_PLASTIC_SC_VINYLTOLUENE                    283   //G4_PLASTIC_SC_VINYLTOLUENE
301   //G4_SODIUM_IODIDE                              284   //G4_SODIUM_IODIDE
302   //G4_STILBENE                                   285   //G4_STILBENE
303   //G4_lAr                                        286   //G4_lAr
304                                                << 287   //G4_PbWO4
305   //G4_PbWO4 - CMS value                       << 
306   g4MatNames.push_back("G4_PbWO4");            << 
307   g4MatData.push_back(0.0333333*mm/MeV);       << 
308                                                << 
309   //G4_Lucite                                     288   //G4_Lucite
310                                                   289 
                                                   >> 290   nG4Birks = g4MatData.size();
311 }                                                 291 }
312                                                   292 
313 //....oooOO0OOooo........oooOO0OOooo........oo    293 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
314                                                   294