Generator.h

Go to the documentation of this file.

//! \file
//!
//! Header file for Generator class used for defining the OLYMPUS
//! event generator.
//!
//! This defines the Generator class and the member routines which
//! generate the primary events.
//!
//! \author D.K. Hasell
//! \version 1.0
//! \date 2010-10-14
//!
//! \ingroup event

// *+****1****+****2****+****3****+****4****+****5****+****6****+****7****+****

#ifndef GENERATOR_H
#define GENERATOR_H 

// Include the GEANT4 header files referenced here.

#include "G4VUserPrimaryGeneratorAction.hh"
#include "G4ThreeVector.hh"

// Declare the classes referenced here.

class G4Event;
class G4ParticleGun;

class ep_Kinematics;

//! Generator class.

class Generator : public G4VUserPrimaryGeneratorAction {

public:

     //! Constructor.

     Generator();

     //! Destructor.

     ~Generator();

     // Member functions.

     //! Generate the requested primary particles.

     void GeneratePrimaries( G4Event * event );

     //! Return a location in the target following a triangular distribution.

     G4ThreeVector Target_Distribution();

     //! Set/Get beam kinetic energy.

     G4double setEbeam( G4double );
     G4double getEbeam();

     //! Set/Get beam charge.

     G4int setBeamCharge( G4int );
     G4int getBeamCharge();

     //! Set/Get ep Theta min.

     G4double setepThetaMin( G4double );
     G4double getepThetaMin();

     //! Set/Get ep Theta max.

     G4double setepThetaMax( G4double );
     G4double getepThetaMax();

     //! Set/Get ep Phi min.

     G4double setepPhiMin( G4double );
     G4double getepPhiMin();

     //! Set/Get ep Phi max.

     G4double setepPhiMax( G4double );
     G4double getepPhiMax();

     //! Print

     void Print();

private:

     G4ParticleGun * particleGun;

     G4double Ebeam;
     G4int BeamCharge;

     G4bool ep_Elastic;

     ep_Kinematics * epKin;

     G4double epTheta_Min;
     G4double epTheta_Max;
     G4double cos_epThMin;
     G4double cos_epThMax;
     G4double epPhi_Min;
     G4double epPhi_Max;

};

//! Class for calculating ep elastic kinematics.

class ep_Kinematics {

public:

     //! Constructors.

     ep_Kinematics( G4double Ebeam );

     //! Destructor.

     ~ep_Kinematics();

     //! Return kinematic variables as a function of electron scattering angle.

     void for_e_angle( G4double theta3, G4double & E3, G4double & E4,
                       G4double & Q2, G4double & theta4 );

     //! Print

     void Print();

     // Data members.

private:

     G4double E1;
     G4double p1;
     G4double beta;
     G4double E3cmbyG;

};

//! Calculate the ep elastic cross section.

G4double epCrossSection( const G4double & E, const G4double & theta3,
                         const G4double & Ep, const G4double & Q2 );

#endif