G4LEPTSElossDistr Class Reference

#include <G4LEPTSElossDistr.hh>

Public Member Functions
	G4LEPTSElossDistr (const G4String &)
void	ReadFile ()
G4double	Sample (G4double, G4double)
G4bool	IsFileFound () const

Detailed Description

Definition at line 40 of file G4LEPTSElossDistr.hh.

Constructor & Destructor Documentation

G4LEPTSElossDistr()

G4LEPTSElossDistr::G4LEPTSElossDistr

(

const G4String &

fname

)

Definition at line 37 of file G4LEPTSElossDistr.cc.

{
  fileName = fname;
  ReadFile();
}

Member Function Documentation

IsFileFound()

G4bool G4LEPTSElossDistr::IsFileFound ( ) const

inline

Definition at line 48 of file G4LEPTSElossDistr.hh.

{ return bFileFound; }

ReadFile()

void G4LEPTSElossDistr::ReadFile

(

)

Definition at line 43 of file G4LEPTSElossDistr.cc.

{
  theNDistributions = 0;
 
  FILE * fp;
 
  if ((fp=fopen(fileName.c_str(), "r"))==nullptr)
  {
    NoBins = 0;
    bFileFound = false;
    return;
  } 
 
  bFileFound = true;
 
  G4int nEnergies;
  G4int nAngles;
  G4int nData;
  (void) fscanf(fp,"%i \n",&nEnergies);
  for( G4int ie = 0; ie < nEnergies; ++ie )
  {
    G4float energySep; 
    (void) fscanf(fp,"%f \n",&energySep);
    (void) fscanf(fp,"%i \n",&nAngles);
    for( G4int ia = 0; ia < nAngles; ++ia )
    {
      G4float angleSep; 
      (void) fscanf(fp,"%f \n",&angleSep);
      auto dist = new G4LEPTSDistribution();
      ++theNDistributions;
      std::map<G4double, G4LEPTSDistribution *> angleDist;
      angleDist[angleSep] = dist;
      theDistributions[energySep] = std::move(angleDist);
      
      (void) fscanf(fp,"%i \n",&nData);
      if( dist->ReadFile( fp, nData ) )
      {
        std::ostringstream message;
        message << "End of file found while reading file: "
                << fileName;
        G4Exception("G4LEPTSElossDistr::ReadFile()", "ReadError",
                    FatalException, message); 
      }
    }
  }
  
  fclose(fp);
}

Referenced by G4LEPTSElossDistr().

Sample()

G4double G4LEPTSElossDistr::Sample	(	G4double	eMin,
		G4double	eMax )

Definition at line 92 of file G4LEPTSElossDistr.cc.

{
  // Sample Energy from Cumulative distr. G4interval [eMin, eMax]
 
  if( eMin > eMax) return 0.0;
 
  // Get the distribution to do the sampling
  G4LEPTSDistribution* distr  = nullptr;
  if( theNDistributions == 1 )
  {
    distr = (*( (*(theDistributions.cbegin())).second ).cbegin()).second;
  }
  else
  {
    for( auto itedd = theDistributions.cbegin(); itedd != theDistributions.cend(); ++itedd ){
      G4double energySep = (*itedd).first;
      if( eMax < energySep )
      {
        const auto& dist1 = (*itedd).second;
        for( auto ited = dist1.cbegin(); ited != dist1.cend(); ++ited )
        {
          G4double angleSep = (*ited).first;
      if( 1 < angleSep )
          {
            distr = (*ited).second;
            break;
          }
        }
        break;
      }
    }
  }
 
  return (nullptr != distr) ? distr->Sample(eMin, eMax) : 0.0;
}

---

The documentation for this class was generated from the following files:

• G4LEPTSElossDistr.hh
• G4LEPTSElossDistr.cc