Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/dna/models/src/G4LEPTSElossDistr.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/dna/models/src/G4LEPTSElossDistr.cc (Version 11.3.0) and /processes/electromagnetic/dna/models/src/G4LEPTSElossDistr.cc (Version 10.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 // G4LEPTSElossDistr                              
 27 //                                                
 28 // Author: Pedro Arce (CIEMAT), 2014              
 29 // -------------------------------------------    
 30                                                   
 31 #include "G4LEPTSElossDistr.hh"                   
 32 #include "G4LEPTSDistribution.hh"                 
 33                                                   
 34 #include <stdio.h>                                
 35 #include <iostream>                               
 36                                                   
 37 G4LEPTSElossDistr::G4LEPTSElossDistr(const G4S    
 38 {                                                 
 39   fileName = fname;                               
 40   ReadFile();                                     
 41 }                                                 
 42                                                   
 43 void G4LEPTSElossDistr::ReadFile()                
 44 {                                                 
 45   theNDistributions = 0;                          
 46                                                   
 47   FILE * fp;                                      
 48                                                   
 49   if ((fp=fopen(fileName.c_str(), "r"))==nullp    
 50   {                                               
 51     NoBins = 0;                                   
 52     bFileFound = false;                           
 53     return;                                       
 54   }                                               
 55                                                   
 56   bFileFound = true;                              
 57                                                   
 58   G4int nEnergies;                                
 59   G4int nAngles;                                  
 60   G4int nData;                                    
 61   (void) fscanf(fp,"%i \n",&nEnergies);           
 62   for( G4int ie = 0; ie < nEnergies; ++ie )       
 63   {                                               
 64     G4float energySep;                            
 65     (void) fscanf(fp,"%f \n",&energySep);         
 66     (void) fscanf(fp,"%i \n",&nAngles);           
 67     for( G4int ia = 0; ia < nAngles; ++ia )       
 68     {                                             
 69       G4float angleSep;                           
 70       (void) fscanf(fp,"%f \n",&angleSep);        
 71       auto dist = new G4LEPTSDistribution();      
 72       ++theNDistributions;                        
 73       std::map<G4double, G4LEPTSDistribution *    
 74       angleDist[angleSep] = dist;                 
 75       theDistributions[energySep] = std::move(    
 76                                                   
 77       (void) fscanf(fp,"%i \n",&nData);           
 78       if( dist->ReadFile( fp, nData ) )           
 79       {                                           
 80         std::ostringstream message;               
 81         message << "End of file found while re    
 82                 << fileName;                      
 83         G4Exception("G4LEPTSElossDistr::ReadFi    
 84                     FatalException, message);     
 85       }                                           
 86     }                                             
 87   }                                               
 88                                                   
 89   fclose(fp);                                     
 90 }                                                 
 91                                                   
 92 G4double G4LEPTSElossDistr::Sample( G4double e    
 93 {                                                 
 94   // Sample Energy from Cumulative distr. G4in    
 95                                                   
 96   if( eMin > eMax) return 0.0;                    
 97                                                   
 98   // Get the distribution to do the sampling      
 99   G4LEPTSDistribution* distr  = nullptr;          
100   if( theNDistributions == 1 )                    
101   {                                               
102     distr = (*( (*(theDistributions.cbegin()))    
103   }                                               
104   else                                            
105   {                                               
106     for( auto itedd = theDistributions.cbegin(    
107       G4double energySep = (*itedd).first;        
108       if( eMax < energySep )                      
109       {                                           
110         const auto& dist1 = (*itedd).second;      
111         for( auto ited = dist1.cbegin(); ited     
112         {                                         
113           G4double angleSep = (*ited).first;      
114     if( 1 < angleSep )                            
115           {                                       
116             distr = (*ited).second;               
117             break;                                
118           }                                       
119         }                                         
120         break;                                    
121       }                                           
122     }                                             
123   }                                               
124                                                   
125   return (nullptr != distr) ? distr->Sample(eM    
126 }                                                 
127