Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/hadronic/models/coherent_elastic/src/G4ChipsElasticModel.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/models/coherent_elastic/src/G4ChipsElasticModel.cc (Version 11.3.0) and /processes/hadronic/models/coherent_elastic/src/G4ChipsElasticModel.cc (Version 8.1)


  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 //                                                
 29 // Geant4 Class : G4ChipsElasticModel             
 30 //                                                
 31 // Author : V.Ivanchenko 29 June 2009             
 32 //                                                
 33 // Modified:                                      
 34 // 13.01.10: M.Kosov: Use G4Q(Pr/Neut)ElasticC    
 35 //                                                
 36 //--------------------------------------------    
 37 // CHIPS model of hadron elastic scattering       
 38 //                                                
 39                                                   
 40 #include "G4ChipsElasticModel.hh"                 
 41 #include "G4ParticleDefinition.hh"                
 42 #include <iostream>                               
 43                                                   
 44 #include "G4CrossSectionDataSetRegistry.hh"       
 45                                                   
 46 G4ChipsElasticModel::G4ChipsElasticModel() : G    
 47 {                                                 
 48     pxsManager    = (G4ChipsProtonElasticXS*)G    
 49     nxsManager    = (G4ChipsNeutronElasticXS*)    
 50     PBARxsManager = (G4ChipsAntiBaryonElasticX    
 51     PIPxsManager  = (G4ChipsPionPlusElasticXS*    
 52     PIMxsManager  = (G4ChipsPionMinusElasticXS    
 53     KPxsManager   = (G4ChipsKaonPlusElasticXS*    
 54     KMxsManager   = (G4ChipsKaonMinusElasticXS    
 55   //Description();                                
 56 }                                                 
 57                                                   
 58 G4ChipsElasticModel::~G4ChipsElasticModel()       
 59 {}                                                
 60                                                   
 61 void G4ChipsElasticModel::ModelDescription(std    
 62 {                                                 
 63                                                   
 64     outFile << "The G4ChipsElasticModel model     
 65             << "scattering using the parameter    
 66             << "of M. Kossov\n";                  
 67                                                   
 68 }                                                 
 69                                                   
 70                                                   
 71 G4double                                          
 72 G4ChipsElasticModel::SampleInvariantT(const G4    
 73          G4double plab, G4int Z, G4int A)         
 74 {                                                 
 75   G4int N = A - Z;                                
 76   if(Z == 1 && N == 2)      { N = 1; }            
 77   else if(Z == 2 && N == 1) { N = 2; }            
 78   G4int projPDG = p->GetPDGEncoding();            
 79   G4double cs = 0.;                               
 80   if     (projPDG==2212) { cs = pxsManager->Ge    
 81   else if(projPDG==2112) { cs = nxsManager->Ge    
 82   else if(projPDG==-2212){ cs = PBARxsManager-    
 83   else if(projPDG== 211) { cs = PIPxsManager->    
 84   else if(projPDG==-211) { cs = PIMxsManager->    
 85   else if(projPDG== 321) { cs = KPxsManager->G    
 86   else if(projPDG==-321) { cs = KMxsManager->G    
 87                                                   
 88   G4double t = 0.0;                               
 89   if(cs > 0.0)                                    
 90   {                                               
 91     if     (projPDG== 2212) { t = pxsManager->    
 92     else if(projPDG== 2112) { t = nxsManager->    
 93     else if(projPDG==-2212) { t = PBARxsManage    
 94     else if(projPDG==  211) { t = PIPxsManager    
 95     else if(projPDG== -211) { t = PIMxsManager    
 96     else if(projPDG==  321) { t = KPxsManager-    
 97     else if(projPDG== -321) { t = KMxsManager-    
 98   }                                               
 99   else  { t = G4HadronElastic::SampleInvariant    
100   return t;                                       
101 }                                                 
102                                                   
103