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

(

std::string

file

)

Definition at line 30 of file G4LEPTSElossDistr.cc.

                                                   {
  fileName = file;

  ReadFile();
}

Here is the call graph for this function:

Member Function Documentation

G4bool G4LEPTSElossDistr::IsFileFound ( ) const

inline

Definition at line 50 of file G4LEPTSElossDistr.hh.

                             {
    return bFileFound;
  }

void G4LEPTSElossDistr::ReadFile

(

)

Definition at line 37 of file G4LEPTSElossDistr.cc.

{
  theNDistributions = 0;

  FILE * fp;

  if ((fp=fopen(fileName.c_str(), "r"))==NULL){
    //G4cout << "Error reading " << fileName << G4endl;
    NoBins = 0;
    bFileFound = false;
    return;
  } 

  bFileFound = true;
  //  G4cout << "Read Eloss Distro (" << fileName << ") " << G4endl;
  G4int nEnergies;
  G4int nAngles;
  G4int nData;
  fscanf(fp,"%i \n",&nEnergies);
  for( G4int ie = 0; ie < nEnergies; ie++ ){
    float energySep; 
    fscanf(fp,"%f \n",&energySep);
    fscanf(fp,"%i \n",&nAngles);
    for( G4int ia = 0; ia < nAngles; ia++ ){
      float angleSep; 
      fscanf(fp,"%f \n",&angleSep);
      G4LEPTSDistribution* dist = new G4LEPTSDistribution();
      theNDistributions ++;
      mddist angleDist;
      angleDist[angleSep] = dist;
      theDistributions[energySep] = angleDist;
      
      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);

}

Here is the call graph for this function:

Here is the caller graph for this function:

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  = 0;
  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);
}

Here is the call graph for this function:

---

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

• geant4.10.03.p01/source/processes/electromagnetic/dna/models/include/G4LEPTSElossDistr.hh
• geant4.10.03.p01/source/processes/electromagnetic/dna/models/src/G4LEPTSElossDistr.cc