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 ]

  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 // 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 G4double * aIn, 
 36                    const G4double * anE, G4int nP)
 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, G4double > > aP;
 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 kineticEnergy)
 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) {++it;}  /* Loop checking, 08.01.2016, W. Pokorski */
 68   if(it==end()) 
 69   {
 70     G4ExceptionDescription ed;
 71     ed << "This cross section is applied for E(MeV)= " << kineticEnergy
 72        << " outside allowed energy interval" << G4endl; 
 73     G4Exception("G4PiData::ReactionXSection", "had001", FatalException, ed);
 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)*(x2-x1)/(e2-e1));
 82 }
 83 
 84 ////////////////////////////////////////////////////////////////////////////
 85 
 86 G4double G4PiData::ElasticXSection(G4double kineticEnergy)
 87 {
 88   if (kineticEnergy <= fMinE) { return fTot0 - fInel0; }
 89   if (kineticEnergy >= fMaxE) { return fTot1 - fInel1; }
 90 
 91   G4PiData::iterator it=begin();
 92   while(it!=end()&&kineticEnergy>(*it).first) {it++;}  /* Loop checking, 08.01.2016, W. Pokorski */
 93   if(it==end()) 
 94   {
 95     G4ExceptionDescription ed;
 96     ed << "This cross section is applied for E(MeV)= " << kineticEnergy
 97        << " outside allowed energy interval" << G4endl; 
 98     G4Exception("G4PiData::ElasticXSection", "had001", FatalException, ed);
 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.second;
104   e2=(*(it)).first;
105   x2=(*(it)).second.first - (*(it)).second.second;
106   return std::max(0., x1 + (kineticEnergy-e1)*(x2-x1)/(e2-e1));
107 }
108 
109 ////////////////////////////////////////////////////////////////////////////
110 
111 G4double G4PiData::TotalXSection(G4double kineticEnergy)
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) {it++;}  /* Loop checking, 08.01.2016, W. Pokorski */
118   if(it==end()) 
119   {
120     G4ExceptionDescription ed;
121     ed << "This cross section is applied for E(MeV)= " << kineticEnergy
122        << " outside allowed energy interval" << G4endl; 
123     G4Exception("G4PiData::TotalXSection", "had001", FatalException, ed);
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)*(x2-x1)/(e2-e1));
132 }
133