Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/highenergy/src/G4eeToPGammaModel.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/highenergy/src/G4eeToPGammaModel.cc (Version 11.3.0) and /processes/electromagnetic/highenergy/src/G4eeToPGammaModel.cc (Version 1.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 header file                       
 30 //                                                
 31 //                                                
 32 // File name:     G4eeToPGammaModel               
 33 //                                                
 34 // Author:        Vladimir Ivanchenko             
 35 //                                                
 36 // Creation date: 25.10.2003                      
 37 //                                                
 38 // Modifications:                                 
 39 //                                                
 40 //                                                
 41 // -------------------------------------------    
 42 //                                                
 43                                                   
 44                                                   
 45 //....oooOO0OOooo........oooOO0OOooo........oo    
 46 //....oooOO0OOooo........oooOO0OOooo........oo    
 47                                                   
 48 #include "G4eeToPGammaModel.hh"                   
 49 #include "Randomize.hh"                           
 50 #include "G4PhysicalConstants.hh"                 
 51 #include "G4SystemOfUnits.hh"                     
 52 #include "G4PionZero.hh"                          
 53 #include "G4Eta.hh"                               
 54 #include "G4Gamma.hh"                             
 55 #include "G4DynamicParticle.hh"                   
 56 #include "G4PhysicsVector.hh"                     
 57 #include "G4PhysicsLinearVector.hh"               
 58 #include "G4eeCrossSections.hh"                   
 59                                                   
 60 //....oooOO0OOooo........oooOO0OOooo........oo    
 61                                                   
 62 using namespace std;                              
 63                                                   
 64 G4eeToPGammaModel::G4eeToPGammaModel(G4eeCross    
 65              const G4String& npart,               
 66              G4double maxkinEnergy,               
 67              G4double binWidth)                   
 68 : G4Vee2hadrons(cr,                               
 69     npart=="pi0" ? 782.62*MeV:1019.46*MeV,        
 70     maxkinEnergy,                                 
 71     binWidth)                                     
 72 {                                                 
 73   G4cout << "####G4eeToPGammaModel & particle:    
 74    << "####" << G4endl;                           
 75                                                   
 76   pi0 = G4PionZero::PionZero();                   
 77   if(npart == "pi0") {                            
 78     massR = 782.62*MeV;                           
 79     particle = pi0;                               
 80   } else {                                        
 81     massR = 1019.46*MeV;                          
 82     particle = G4Eta::Eta();                      
 83   }                                               
 84   massP = particle->GetPDGMass();                 
 85 }                                                 
 86                                                   
 87 //....oooOO0OOooo........oooOO0OOooo........oo    
 88                                                   
 89 G4eeToPGammaModel::~G4eeToPGammaModel()           
 90 {}                                                
 91                                                   
 92 //....oooOO0OOooo........oooOO0OOooo........oo    
 93                                                   
 94 G4double G4eeToPGammaModel::PeakEnergy() const    
 95 {                                                 
 96   return massR;                                   
 97 }                                                 
 98                                                   
 99 //....oooOO0OOooo........oooOO0OOooo........oo    
100                                                   
101 G4double G4eeToPGammaModel::ComputeCrossSectio    
102 {                                                 
103   G4double xs;                                    
104   if(particle == pi0) xs = cross->CrossSection    
105   else                xs = cross->CrossSection    
106   return xs;                                      
107 }                                                 
108                                                   
109 //....oooOO0OOooo........oooOO0OOooo........oo    
110                                                   
111 void G4eeToPGammaModel::SampleSecondaries(std:    
112       G4double e, const G4ThreeVector& directi    
113 {                                                 
114   G4double egam = 0.5*e*(1.0 - massP*massP/(ma    
115   G4double tkin = e - egam - massP;               
116   if(tkin < 0.0) tkin = 0.0;                      
117   G4double cost;                                  
118   do {                                            
119     cost = 2.0*G4UniformRand() - 1.0;             
120     // Loop checking, 07-Aug-2015, Vladimir Iv    
121   } while( 2.0*G4UniformRand() > 1.0 + cost*co    
122                                                   
123   G4double sint = sqrt(1.0 - cost*cost);          
124   G4double phi  = twopi * G4UniformRand();        
125                                                   
126   G4ThreeVector dir(sint*cos(phi),sint*sin(phi    
127   dir.rotateUz(direction);                        
128                                                   
129   // create G4DynamicParticle objects             
130   G4DynamicParticle* p1 =                         
131      new G4DynamicParticle(particle,dir,tkin);    
132   G4DynamicParticle* p2 =                         
133      new G4DynamicParticle(G4Gamma::Gamma(),-d    
134   newp->push_back(p1);                            
135   newp->push_back(p2);                            
136 }                                                 
137                                                   
138 //....oooOO0OOooo........oooOO0OOooo........oo    
139                                                   
140