MollerPrimaryGenAction.hh

Go to the documentation of this file.

#ifndef MollerPrimaryGenAction_h
#define MollerPrimaryGenAction_h 1

#include "G4VUserPrimaryGeneratorAction.hh"
#include "globals.hh"
#include "Randomize.hh"

#include "Rtypes.h"

class G4ParticleGun;
class G4Event;
class MollerPrimaryGeneratorMessenger;


typedef struct {
  double theta;
  double y;
  double xs;
}cms_cross_section_t;

typedef struct {
  double theta;
  double Eprime;
  double xs;
}double_differential_t;

class MollerPrimaryGenAction : public G4VUserPrimaryGeneratorAction
{
 
public:
  MollerPrimaryGenAction();
  ~MollerPrimaryGenAction();

public:
  void GeneratePrimaries(G4Event* anEvent);
  void GeneratePrimaries_moller_dustin(G4Event* anEvent);
  void GeneratePrimaries_elasticep_dustin(G4Event* anevent);
  void GeneratePrimaries_inelasticep_dustin(G4Event* anevent);
  void GeneratePrimaries_phasespace_mollers(G4Event* anevent);
  void GeneratePrimaries_phasespace_motts(G4Event* anevent);
  void GeneratePrimaries_pion(G4Event* anevent);
  void SetGenerator(G4int val) { GeneratorNumber = val;}
  G4int GetGenerator() {return GeneratorNumber;}
  double EnergNumInt(double btt, double a0, double b0);
  double RadProfile(double eloss, double btt);

  void setSeedValue (G4bool value) { fSetSeedConst = value; 
    G4cout<<"Setting fSetSeedConst"<<G4endl;};

  Double_t wiser_sigma(Double_t ebeam, Double_t pf, Double_t thf, Double_t rl, Int_t type);
  static Double_t wiser_all_fit( Double_t *x, Double_t *par );


private:
  G4ParticleGun*                    particleGun;
  MollerPrimaryGeneratorMessenger*  gunMessenger;   //messenger of this class
  G4int                             GeneratorNumber; // Number of generator to use
  G4int                             eventnumber;
  static const G4int nbin = 1000;
  static const G4int nbinIneEP = 100;
  cms_cross_section_t XStable[nbin];
  double_differential_t DDtable[nbinIneEP][nbinIneEP];
  const G4double ebeam;
  const G4double mp;
  const G4double bremcut; //lower limit of bremsstrahlung (1 keV)
  const G4double me;


  G4bool fSetSeedConst;
};

#endif