MollerPrimaryGenAction_inelasticep.cc File Reference

#include "G4ParticleGun.hh"
#include "G4Electron.hh"
#include "G4Event.hh"
#include "globals.hh"
#include "MollerPrimaryGenAction.hh"
#include "RootAnalysis.hh"

Include dependency graph for MollerPrimaryGenAction_inelasticep.cc:

Go to the source code of this file.

Defines
#define	r0   2.818
#define	pi   3.14159
#define	me   0.511E-3
#define	Euler   0.5772157
#define	d2r   pi/180.
#define	r2d   180/pi
#define	gevfm   0.1973
#define	pmu   2.79
#define	alpha   0.007299
Functions
double	SigT (double k)

---

Define Documentation

#define alpha   0.007299

Definition at line 37 of file MollerPrimaryGenAction_inelasticep.cc.

#define d2r   pi/180.

Definition at line 32 of file MollerPrimaryGenAction_inelasticep.cc.

#define Euler   0.5772157

Definition at line 31 of file MollerPrimaryGenAction_inelasticep.cc.

#define gevfm   0.1973

Definition at line 35 of file MollerPrimaryGenAction_inelasticep.cc.

#define me   0.511E-3

Definition at line 30 of file MollerPrimaryGenAction_inelasticep.cc.

#define pi   3.14159

Definition at line 29 of file MollerPrimaryGenAction_inelasticep.cc.

#define pmu   2.79

Definition at line 36 of file MollerPrimaryGenAction_inelasticep.cc.

#define r0   2.818

Definition at line 28 of file MollerPrimaryGenAction_inelasticep.cc.

#define r2d   180/pi

Definition at line 33 of file MollerPrimaryGenAction_inelasticep.cc.

---

Function Documentation

double SigT

(

double

k

)

Definition at line 314 of file MollerPrimaryGenAction_inelasticep.cc.

Referenced by MollerPrimaryGenAction::GeneratePrimaries_inelasticep_dustin().

                      {
  double sigt ; 
  sigt=0 ;
  /* histogram-type fit to world data from PDG over Delta, S11, tc.*/
  /* switch to linear interpolation
   * rather than block
  if(k > 0.175) sigt=0.08 ;
  if(k > 0.196) sigt=0.12 ;
  if(k > 0.214) sigt=0.17 ;
  if(k > 0.233) sigt=0.21 ;
  if(k > 0.251) sigt=0.32 ;
  if(k > 0.251) sigt=0.32 ;
  if(k > 0.251) sigt=0.32 ;
  if(k > 0.271) sigt=0.40 ;
  if(k > 0.290) sigt=0.50 ;
  if(k > 0.310) sigt=0.53 ;
  if(k > 0.330) sigt=0.50 ;
  if(k > 0.350) sigt=0.45 ;
  if(k > 0.370) sigt=0.41 ;
  if(k > 0.390) sigt=0.32 ;
  if(k > 0.410) sigt=0.24 ;
  if(k > 0.430) sigt=0.22 ;
  if(k > 0.450) sigt=0.18 ;
  if(k > 0.470) sigt=0.17 ;
  if(k > 0.520) sigt=0.18 ;
  if(k > 0.550) sigt=0.21 ;
  if(k > 0.600) sigt=0.22 ;
  if(k > 0.650) sigt=0.24 ;
  if(k > 0.700) sigt=0.27 ;
  if(k > 0.750) sigt=0.25 ;
  if(k > 0.800) sigt=0.25 ;
  if(k > 0.750) sigt=0.25 ;
  if(k > 0.800) sigt=0.22 ;
  if(k > 0.900) sigt=0.21 ;
  if(k > 0.950) sigt=0.215;
  if(k > 1.000) sigt=0.225;
  if(k > 1.050) sigt=0.215;
  if(k > 1.100) sigt=0.185;
  if(k > 1.150) sigt=0.165;
  if(k > 1.200) sigt=0.15 ;
  if(k > 1.300) sigt=0.16 ;
  if(k > 1.400) sigt=0.15 ;
  */
  /* use fit to sigt (in mb) from Caldwell */

  static const int nbin = 37;


  double kval[nbin] = { 0.175, 0.196, 0.214, 0.233, 0.251, 0.251, 0.251, 0.271, 0.290, 0.310, 0.330, 0.350, 0.370, 0.390, 0.410, 0.430, 0.450, 0.470, 0.520, 0.550, 0.600, 0.650, 0.700, 0.750, 0.800, 0.750, 0.800, 0.900, 0.950, 1.000, 1.050, 1.100, 1.150, 1.200, 1.300, 1.400, 1.5 };

  double tval[nbin] = { 0.08, 0.12, 0.17, 0.21, 0.32, 0.32, 0.32, 0.40, 0.50, 0.53, 0.50, 0.45, 0.41, 0.32, 0.24, 0.22, 0.18, 0.17, 0.18, 0.21, 0.22, 0.24, 0.27, 0.25, 0.25, 0.25, 0.22, 0.21, 0.215, 0.225, 0.215, 0.185, 0.165, 0.15, 0.16, 0.15, 0.15};


  int idx = 0;

  if(k > 1.5){
      sigt=0.0987 + 0.0649/sqrt(k) ;
  } else { 

      while( !( kval[idx] < k && k <= kval[idx+1] ) && idx < nbin-1 ){ idx++; }

      // Didn't find a good value
      if( idx == nbin-1 ) { 
    assert( k < 0.175 );
    sigt = 0.0; 
      }
      else {
    assert( 0 <= idx && idx < nbin-1 );

    double intval = (k - kval[idx])/(kval[idx+1]-kval[idx]);

    assert( 0.0 <= intval && intval < 1.0 );

    sigt = tval[idx]*(1.0-intval) + tval[idx+1]*intval;
      }
  }

  return sigt ;
}      

Here is the caller graph for this function: