G4LEPTSDistribution Class Reference

#include <G4LEPTSDistribution.hh>

Collaboration diagram for G4LEPTSDistribution:

Public Member Functions
	G4LEPTSDistribution ()
void	ReadFile (G4String fileName)
G4bool	ReadFile (FILE *fp, G4int nData)
G4double	Sample (G4double, G4double)
G4bool	IsFileFound () const

Private Attributes
int	NoBins
G4double	E [NMAX]
G4double	f [NMAX]
G4double	F [NMAX]
G4double	eF [NMAX]
G4bool	bFileFound

Detailed Description

Definition at line 36 of file G4LEPTSDistribution.hh.

Constructor & Destructor Documentation

G4LEPTSDistribution()

G4LEPTSDistribution::G4LEPTSDistribution

(

)

Definition at line 29 of file G4LEPTSDistribution.cc.

{
}

Member Function Documentation

IsFileFound()

G4bool G4LEPTSDistribution::IsFileFound ( ) const

inline

Definition at line 45 of file G4LEPTSDistribution.hh.

                             {
    return bFileFound;
  }

ReadFile() [1/2]

void G4LEPTSDistribution::ReadFile

(

G4String

fileName

)

Definition at line 34 of file G4LEPTSDistribution.cc.

                                                    {

  G4int eB, out, out2;
  float  float_data1,float_data2;
  G4double sum, esum;
  FILE * fp;

  for (eB=0;eB<10000;eB++){
    E[eB]=0.0;
    f[eB]=0.0;
    F[eB]=0.0;
    eF[eB]=0.0;
  }

  if ((fp=fopen(fileName.c_str(), "r"))==NULL){
    //G4cout << "Error reading " << fileName << G4endl;
    NoBins = 0;
    bFileFound = false;
    return;
  }
  else{
    bFileFound = true;
    //    G4cout << "Read Distro (" << fileName << ") " << G4endl;
    out=1;
    eB=1;
    while (out==1){
      out  = fscanf(fp,"%f \n",&float_data1);
      out2 = fscanf(fp,"%f \n",&float_data2);
      if (out==1 && out2==1){
    E[eB]=(G4double)float_data1;
    f[eB]=(G4double)float_data2;
    eB++;
      }
    }

    fclose(fp);
  }

  NoBins=eB-1;  //=1272+1 or 9607+1;

  if( NoBins >= NMAX )
    printf("ERROR !!!!  Eloss NoBins= %d \n", NoBins);

  sum=0.0;
  esum=0.0;
  for (eB=0;eB<=NoBins;eB++) {
    if( f[eB] > 0) {
      sum+=f[eB];
      esum+=E[eB]*f[eB];
    }
    F[eB]=sum;
    eF[eB]=esum;
  }

  //  if( verboseLevel >= 1 ) G4cout << "Norm: " << F[NoBins] << " NoBins: "<< NoBins << G4endl;
  
  for (eB=0;eB<=NoBins;eB++) {
    eF[eB] = eF[eB]/F[eB];
    F[eB] = F[eB]/F[NoBins];
  }
  //for (eB=0;eB<=NoBins;eB++)
  //G4cout << "eff " << E[eB] << " " << f[eB] << " " << F[eB] << "\n";
}

Here is the call graph for this function:

Here is the caller graph for this function:

ReadFile() [2/2]

G4bool G4LEPTSDistribution::ReadFile	(	FILE *	fp,
		G4int	nData
	)

Definition at line 98 of file G4LEPTSDistribution.cc.

{

  G4int eB, out, out2;
  float  float_data1,float_data2;
  G4double sum, esum;

  for (eB=0;eB<10000;eB++){
    E[eB]=0.0;
    f[eB]=0.0;
    F[eB]=0.0;
    eF[eB]=0.0;
  }

  bFileFound = true;
  out=1;
  eB=1;

  for( G4int id = 0; id < nData; id++ ){      
    out  = fscanf(fp,"%f \n",&float_data1);
    out2 = fscanf(fp,"%f \n",&float_data2);
    if (out==1 && out2==1){
      E[eB]=(G4double)float_data1;
      f[eB]=(G4double)float_data2;
      eB++;
    }else{
      return 1;
    }
  }
    
  NoBins=eB-1;  //=1272+1 or 9607+1;

  if( NoBins >= NMAX )
    printf("ERROR !!!!  Eloss NoBins= %d \n", NoBins);

  sum=0.0;
  esum=0.0;
  for (eB=0;eB<=NoBins;eB++) {
    if( f[eB] > 0) {
      sum+=f[eB];
      esum+=E[eB]*f[eB];
    }
    F[eB]=sum;
    eF[eB]=esum;
  }

  //if( verboseLevel >= 1 ) G4cout << "Norm: " << F[NoBins] << " NoBins: "<< NoBins << G4endl;
  
  for (eB=0;eB<=NoBins;eB++) {
    eF[eB] = eF[eB]/F[eB];
    F[eB] = F[eB]/F[NoBins];
  }
  //for (eB=0;eB<=NoBins;eB++)
  //G4cout << "eff " << E[eB] << " " << f[eB] << " " << F[eB] << "\n";

  return 0;
}

Here is the call graph for this function:

Sample()

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

Definition at line 158 of file G4LEPTSDistribution.cc.

                                                                  {
// Sample Energy from Cumulative distr. G4interval [eMin, eMax]

  if( eMin > eMax) return 0.0;

  G4int i,j,k=0, iMin, iMax;

  i=0; j=NoBins;
  while ((j-i)>1) {
    k=(i+j)/2;
    if( E[k] < eMax ) i=k;
    else              j=k;
  }
  iMax = i;

  i=0; j=NoBins;
  while ((j-i)>1) {
    k=(i+j)/2;
    if( E[k] < eMin ) i=k;
    else              j=k;
  }
  iMin = i;

  G4double rnd = F[iMin] + (F[iMax] - F[iMin]) * G4UniformRand();

  i=0; j=NoBins;
  while ((j-i)>1) {
    k=(i+j)/2;
    if( F[k]<rnd) i=k;
    else          j=k;
  }

  G4double Sampled = E[k];

  if(      Sampled < eMin) Sampled = eMin;
  else if( Sampled > eMax) Sampled = eMax;

  return Sampled;
}

Here is the caller graph for this function:

Member Data Documentation

bFileFound

G4bool G4LEPTSDistribution::bFileFound

private

Definition at line 53 of file G4LEPTSDistribution.hh.

E

G4double G4LEPTSDistribution::E[NMAX]

private

Definition at line 52 of file G4LEPTSDistribution.hh.

eF

G4double G4LEPTSDistribution::eF[NMAX]

private

Definition at line 52 of file G4LEPTSDistribution.hh.

F

G4double G4LEPTSDistribution::F[NMAX]

private

Definition at line 52 of file G4LEPTSDistribution.hh.

f

G4double G4LEPTSDistribution::f[NMAX]

private

Definition at line 52 of file G4LEPTSDistribution.hh.

NoBins

int G4LEPTSDistribution::NoBins

private

Definition at line 50 of file G4LEPTSDistribution.hh.

---

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

• G4LEPTSDistribution.hh
• G4LEPTSDistribution.cc