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 11.2.2)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                         3 // * License and Disclaimer                                           *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th      5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided      6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License      7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/      8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.           9 // * include a list of copyright holders.                             *
 10 // *                                               10 // *                                                                  *
 11 // * Neither the authors of this software syst     11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file      15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitatio     16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                               17 // *                                                                  *
 18 // * This  code  implementation is the result      18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboratio     19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distri     20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  ag     21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publicati     22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Sof     23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // *******************************************     24 // ********************************************************************
 25 //                                                 25 //
 26 // G4LEPTSElossDistr                           << 
 27 //                                             << 
 28 // Author: Pedro Arce (CIEMAT), 2014           << 
 29 // ------------------------------------------- << 
 30                                                << 
 31 #include "G4LEPTSElossDistr.hh"                    26 #include "G4LEPTSElossDistr.hh"
 32 #include "G4LEPTSDistribution.hh"                  27 #include "G4LEPTSDistribution.hh"
 33                                                    28 
 34 #include <stdio.h>                             << 
 35 #include <iostream>                            << 
 36                                                    29 
 37 G4LEPTSElossDistr::G4LEPTSElossDistr(const G4S <<  30 G4LEPTSElossDistr::G4LEPTSElossDistr(std::string file) {
 38 {                                              <<  31   fileName = file;
 39   fileName = fname;                            <<  32 
 40   ReadFile();                                      33   ReadFile();
 41 }                                                  34 }
 42                                                    35 
                                                   >>  36 
 43 void G4LEPTSElossDistr::ReadFile()                 37 void G4LEPTSElossDistr::ReadFile() 
 44 {                                                  38 {
 45   theNDistributions = 0;                           39   theNDistributions = 0;
 46                                                    40 
 47   FILE * fp;                                       41   FILE * fp;
 48                                                    42 
 49   if ((fp=fopen(fileName.c_str(), "r"))==nullp <<  43   if ((fp=fopen(fileName.c_str(), "r"))==nullptr){
 50   {                                            <<  44     //G4cout << "Error reading " << fileName << G4endl;
 51     NoBins = 0;                                    45     NoBins = 0;
 52     bFileFound = false;                            46     bFileFound = false;
 53     return;                                        47     return;
 54   }                                                48   } 
 55                                                    49 
 56   bFileFound = true;                               50   bFileFound = true;
 57                                                <<  51   //  G4cout << "Read Eloss Distro (" << fileName << ") " << G4endl;
 58   G4int nEnergies;                                 52   G4int nEnergies;
 59   G4int nAngles;                                   53   G4int nAngles;
 60   G4int nData;                                     54   G4int nData;
 61   (void) fscanf(fp,"%i \n",&nEnergies);            55   (void) fscanf(fp,"%i \n",&nEnergies);
 62   for( G4int ie = 0; ie < nEnergies; ++ie )    <<  56   for( G4int ie = 0; ie < nEnergies; ie++ ){
 63   {                                            <<  57     float energySep; 
 64     G4float energySep;                         << 
 65     (void) fscanf(fp,"%f \n",&energySep);          58     (void) fscanf(fp,"%f \n",&energySep);
 66     (void) fscanf(fp,"%i \n",&nAngles);            59     (void) fscanf(fp,"%i \n",&nAngles);
 67     for( G4int ia = 0; ia < nAngles; ++ia )    <<  60     for( G4int ia = 0; ia < nAngles; ia++ ){
 68     {                                          <<  61       float angleSep; 
 69       G4float angleSep;                        << 
 70       (void) fscanf(fp,"%f \n",&angleSep);         62       (void) fscanf(fp,"%f \n",&angleSep);
 71       auto dist = new G4LEPTSDistribution();   <<  63       auto  dist = new G4LEPTSDistribution();
 72       ++theNDistributions;                     <<  64       theNDistributions ++;
 73       std::map<G4double, G4LEPTSDistribution * <<  65       mddist angleDist;
 74       angleDist[angleSep] = dist;                  66       angleDist[angleSep] = dist;
 75       theDistributions[energySep] = std::move( <<  67       theDistributions[energySep] = angleDist;
 76                                                    68       
 77       (void) fscanf(fp,"%i \n",&nData);            69       (void) fscanf(fp,"%i \n",&nData);
 78       if( dist->ReadFile( fp, nData ) )        <<  70       if( dist->ReadFile( fp, nData ) ) {
 79       {                                        <<  71   G4Exception("G4LEPTSElossDistr",
 80         std::ostringstream message;            <<  72       "",
 81         message << "End of file found while re <<  73         FatalException,
 82                 << fileName;                   <<  74       ("End of file found while reading file"+ fileName).c_str());  
 83         G4Exception("G4LEPTSElossDistr::ReadFi << 
 84                     FatalException, message);  << 
 85       }                                            75       }
 86     }                                              76     }
 87   }                                                77   }
 88                                                    78   
 89   fclose(fp);                                      79   fclose(fp);
                                                   >>  80 
 90 }                                                  81 }
 91                                                    82 
 92 G4double G4LEPTSElossDistr::Sample( G4double e <<  83 
                                                   >>  84 
                                                   >>  85 G4double G4LEPTSElossDistr::Sample( G4double eMin, G4double eMax) 
 93 {                                                  86 {
 94   // Sample Energy from Cumulative distr. G4in <<  87 // Sample Energy from Cumulative distr. G4interval [eMin, eMax]
 95                                                    88 
 96   if( eMin > eMax) return 0.0;                     89   if( eMin > eMax) return 0.0;
 97                                                    90 
 98   // Get the distribution to do the sampling   <<  91   //Get the distribution to do the sampling
 99   G4LEPTSDistribution* distr  = nullptr;           92   G4LEPTSDistribution* distr  = nullptr;
100   if( theNDistributions == 1 )                 <<  93   if( theNDistributions == 1 ){
101   {                                            <<  94     distr = (*( (*(theDistributions.begin())).second ).begin()).second;
102     distr = (*( (*(theDistributions.cbegin())) <<  95   } else{
103   }                                            <<  96     mdmddist::const_iterator itedd;
104   else                                         <<  97     for( itedd = theDistributions.begin(); itedd != theDistributions.end(); itedd++ ){
105   {                                            << 
106     for( auto itedd = theDistributions.cbegin( << 
107       G4double energySep = (*itedd).first;         98       G4double energySep = (*itedd).first;
108       if( eMax < energySep )                   <<  99       if( eMax < energySep ) {
109       {                                        << 100   //tt  if( eMax <= energySep ) { 
110         const auto& dist1 = (*itedd).second;   << 101   mddist dist1 = (*itedd).second;
111         for( auto ited = dist1.cbegin(); ited  << 102   mddist::const_iterator ited;
112         {                                      << 103   for( ited = dist1.begin(); ited != dist1.end(); ited++ ){
113           G4double angleSep = (*ited).first;   << 104     G4double angleSep = (*ited).first;
114     if( 1 < angleSep )                         << 105     if( 1 < angleSep ) {
115           {                                    << 106       distr = (*ited).second;
116             distr = (*ited).second;            << 107       break;
117             break;                             << 108     }
118           }                                    << 109   }
119         }                                      << 110   break;
120         break;                                 << 
121       }                                           111       }
122     }                                             112     }
123   }                                               113   }
124                                                   114 
125   return (nullptr != distr) ? distr->Sample(eM << 115   //  G4cout << " LEPTSElossDistr::Sample(  " << distr << " NDIST " << theNDistributions << G4endl; //GDEB
                                                   >> 116   return distr->Sample(eMin, eMax);
126 }                                                 117 }
127                                                   118