Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/hadronic/cross_sections/src/G4PiData.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/hadronic/cross_sections/src/G4PiData.cc (Version 11.3.0) and /processes/hadronic/cross_sections/src/G4PiData.cc (Version 1.0)


  1 //                                                  1 
  2 // *******************************************    
  3 // * License and Disclaimer                       
  4 // *                                              
  5 // * The  Geant4 software  is  copyright of th    
  6 // * the Geant4 Collaboration.  It is provided    
  7 // * conditions of the Geant4 Software License    
  8 // * LICENSE and available at  http://cern.ch/    
  9 // * include a list of copyright holders.         
 10 // *                                              
 11 // * Neither the authors of this software syst    
 12 // * institutes,nor the agencies providing fin    
 13 // * work  make  any representation or  warran    
 14 // * regarding  this  software system or assum    
 15 // * use.  Please see the license in the file     
 16 // * for the full disclaimer and the limitatio    
 17 // *                                              
 18 // * This  code  implementation is the result     
 19 // * technical work of the GEANT4 collaboratio    
 20 // * By using,  copying,  modifying or  distri    
 21 // * any work based  on the software)  you  ag    
 22 // * use  in  resulting  scientific  publicati    
 23 // * acceptance of all terms of the Geant4 Sof    
 24 // *******************************************    
 25 //                                                
 26 //                                                
 27 // -------------------------------------------    
 28 // by J.P Wellisch, Sun Sep 15 2002.              
 29                                                   
 30 #include "G4PiData.hh"                            
 31 #include "G4SystemOfUnits.hh"                     
 32                                                   
 33 //////////////////////////////////////////////    
 34                                                   
 35 G4PiData::G4PiData(const G4double * aT, const     
 36                    const G4double * anE, G4int    
 37 {                                                 
 38   for(G4int i = 0; i < nP; ++i )                  
 39   {                                               
 40     std::pair<G4double, G4double> x;              
 41     x.first=aT[i]*millibarn;                      
 42     x.second=aIn[i]*millibarn;                    
 43     std::pair<G4double, std::pair<G4double, G4    
 44     aP.first=anE[i]*GeV;                          
 45     aP.second=x;                                  
 46     push_back(aP);                                
 47     if (i == 0) {                                 
 48       fMinE = aP.first;                           
 49       fTot0 = x.first;                            
 50       fInel0 = x.second;                          
 51     } else if (i+1 == nP) {                       
 52       fMaxE = aP.first;                           
 53       fTot1 = x.first;                            
 54       fInel1 = x.second;                          
 55     }                                             
 56   }                                               
 57 }                                                 
 58                                                   
 59 //////////////////////////////////////////////    
 60                                                   
 61 G4double G4PiData::ReactionXSection(G4double k    
 62 {                                                 
 63   if (kineticEnergy <= fMinE) { return fInel0;    
 64   if (kineticEnergy >= fMaxE) { return fInel1;    
 65                                                   
 66   G4PiData::iterator it=begin();                  
 67   while(it!=end()&&kineticEnergy>(*it).first)     
 68   if(it==end())                                   
 69   {                                               
 70     G4ExceptionDescription ed;                    
 71     ed << "This cross section is applied for E    
 72        << " outside allowed energy interval" <    
 73     G4Exception("G4PiData::ReactionXSection",     
 74   }                                               
 75   if(it==begin()) it++;                           
 76   G4double x1,x2,e1,e2;                           
 77   e1=(*(it-1)).first;                             
 78   x1=(*(it-1)).second.second;                     
 79   e2=(*(it)).first;                               
 80   x2=(*(it)).second.second;                       
 81   return std::max(0., x1 + (kineticEnergy-e1)*    
 82 }                                                 
 83                                                   
 84 //////////////////////////////////////////////    
 85                                                   
 86 G4double G4PiData::ElasticXSection(G4double ki    
 87 {                                                 
 88   if (kineticEnergy <= fMinE) { return fTot0 -    
 89   if (kineticEnergy >= fMaxE) { return fTot1 -    
 90                                                   
 91   G4PiData::iterator it=begin();                  
 92   while(it!=end()&&kineticEnergy>(*it).first)     
 93   if(it==end())                                   
 94   {                                               
 95     G4ExceptionDescription ed;                    
 96     ed << "This cross section is applied for E    
 97        << " outside allowed energy interval" <    
 98     G4Exception("G4PiData::ElasticXSection", "    
 99   }                                               
100   if(it==begin()) it++;                           
101   G4double x1,x2,e1,e2;                           
102   e1=(*(it-1)).first;                             
103   x1=(*(it-1)).second.first - (*(it-1)).second    
104   e2=(*(it)).first;                               
105   x2=(*(it)).second.first - (*(it)).second.sec    
106   return std::max(0., x1 + (kineticEnergy-e1)*    
107 }                                                 
108                                                   
109 //////////////////////////////////////////////    
110                                                   
111 G4double G4PiData::TotalXSection(G4double kine    
112 {                                                 
113   if (kineticEnergy <= fMinE) { return fTot0;     
114   if (kineticEnergy >= fMaxE) { return fTot1;     
115                                                   
116   G4PiData::iterator it=begin();                  
117   while(it!=end()&&kineticEnergy>(*it).first)     
118   if(it==end())                                   
119   {                                               
120     G4ExceptionDescription ed;                    
121     ed << "This cross section is applied for E    
122        << " outside allowed energy interval" <    
123     G4Exception("G4PiData::TotalXSection", "ha    
124   }                                               
125   if(it==begin()) it++;                           
126   G4double x1,x2,e1,e2;                           
127   e1=(*(it-1)).first;                             
128   x1=(*(it-1)).second.first;                      
129   e2=(*(it)).first;                               
130   x2=(*(it)).second.first;                        
131   return std::max(0., x1 + (kineticEnergy-e1)*    
132 }                                                 
133