G4LEPTSElossDistr Class Reference

#include <G4LEPTSElossDistr.hh>

Public Member Functions
	G4LEPTSElossDistr (std::string)
void	ReadFile ()
G4double	Sample (G4double, G4double)
G4bool	IsFileFound () const

Detailed Description

Definition at line 42 of file G4LEPTSElossDistr.hh.

Constructor & Destructor Documentation

G4LEPTSElossDistr()

G4LEPTSElossDistr::G4LEPTSElossDistr

(

std::string

file

)

Definition at line 30 of file G4LEPTSElossDistr.cc.

                                                   {
  fileName = file;
 
  ReadFile();
}

Member Function Documentation

IsFileFound()

G4bool G4LEPTSElossDistr::IsFileFound ( ) const

inline

Definition at line 50 of file G4LEPTSElossDistr.hh.

                             {
    return bFileFound;
  }

ReadFile()

void G4LEPTSElossDistr::ReadFile

(

)

Definition at line 37 of file G4LEPTSElossDistr.cc.

{
  theNDistributions = 0;
 
  FILE * fp;
 
  if ((fp=fopen(fileName.c_str(), "r"))==nullptr){
    //G4cout << "Error reading " << fileName << G4endl;
    NoBins = 0;
    bFileFound = false;
    return;
  } 
 
  bFileFound = true;
  //  G4cout << "Read Eloss Distro (" << fileName << ") " << G4endl;
  G4int nEnergies;
  G4int nAngles;
  G4int nData;
  (void) fscanf(fp,"%i \n",&nEnergies);
  for( G4int ie = 0; ie < nEnergies; ie++ ){
    float energySep; 
    (void) fscanf(fp,"%f \n",&energySep);
    (void) fscanf(fp,"%i \n",&nAngles);
    for( G4int ia = 0; ia < nAngles; ia++ ){
      float angleSep; 
      (void) fscanf(fp,"%f \n",&angleSep);
      auto  dist = new G4LEPTSDistribution();
      theNDistributions ++;
      mddist angleDist;
      angleDist[angleSep] = dist;
      theDistributions[energySep] = angleDist;
      
      (void) fscanf(fp,"%i \n",&nData);
      if( dist->ReadFile( fp, nData ) ) {
    G4Exception("G4LEPTSElossDistr",
          "",
            FatalException,
          ("End of file found while reading file"+ fileName).c_str());  
      }
    }
  }
  
  fclose(fp);
 
}

Referenced by G4LEPTSElossDistr().

Sample()

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

Definition at line 85 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.begin())).second ).begin()).second;
  } else{
    mdmddist::const_iterator itedd;
    for( itedd = theDistributions.begin(); itedd != theDistributions.end(); itedd++ ){
      G4double energySep = (*itedd).first;
      if( eMax < energySep ) {
    //tt    if( eMax <= energySep ) { 
    mddist dist1 = (*itedd).second;
    mddist::const_iterator ited;
    for( ited = dist1.begin(); ited != dist1.end(); ited++ ){
      G4double angleSep = (*ited).first;
      if( 1 < angleSep ) {
        distr = (*ited).second;
        break;
      }
    }
    break;
      }
    }
  }
 
  //  G4cout << " LEPTSElossDistr::Sample(  " << distr << " NDIST " << theNDistributions << G4endl; //GDEB
  return distr->Sample(eMin, eMax);
}

---

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

• G4LEPTSElossDistr.hh
• G4LEPTSElossDistr.cc