#include <ML2PrimaryGenerationAction.hh>

Inheritance diagram for CML2PrimaryGenerationAction:

Collaboration diagram for CML2PrimaryGenerationAction:

Public Member Functions
	CML2PrimaryGenerationAction (void)

	~CML2PrimaryGenerationAction (void)

void	design (G4double accTargetZPosition)

void	GeneratePrimaries (G4Event *anEvent)

void	inizialize (SPrimaryParticle *primaryParticleData)

void	setNRecycling (G4int val)

void	setNLoopsPhSpParticles (G4int val)

void	setNMaxParticlesInRamPhaseSpace (G4int val)

void	setGunMeanEnergy (G4double val)

void	setGunStdEnergy (G4double val)

void	setGunRadious (G4double val)

void	setCalculatedPhaseSpaceFileIN (G4String val)

void	setSourceTypeName (G4String val)

void	setRotation (G4RotationMatrix *val)

G4int	getNrecycling ()

G4int	getSourceTypeName ()

Public Member Functions inherited from G4VUserPrimaryGeneratorAction
	G4VUserPrimaryGeneratorAction ()

virtual	~G4VUserPrimaryGeneratorAction ()

Static Public Member Functions
static CML2PrimaryGenerationAction *	GetInstance (void)

Private Member Functions
void	setGunRandom ()

void	setGunCalculatedPhaseSpace ()

void	GenerateFromRandom ()

void	GenerateFromCalculatedPhaseSpace ()

void	fillParticlesContainer ()

void	applySourceRotation ()

Private Attributes
G4int	nRecycling

G4int	nLoopsPhSpParticles

G4int	nMaxParticlesInRamPhaseSpace

G4int	idParticleSource

G4double	GunMeanEnegy

G4double	GunStdEnegy

G4double	GunRadious

G4String	calculatedPhaseSpaceFileIN

CML2PrimaryGenerationActionMessenger *	PrimaryGenerationActionMessenger

G4double	accTargetZPosition

G4ThreeVector	dir

G4ThreeVector	pos

G4double	ek

G4RotationMatrix *	rm

G4Timer	myTime

G4double	sinTheta

G4double	cosTheta

G4double	phi

G4double	ro

G4double	alfa

G4ParticleGun *	particleGun

G4ParticleDefinition *	gamma

G4ParticleDefinition *	electron

G4ParticleDefinition *	positron

SPrimaryParticle *	primaryParticleData

Sparticle *	particles

Sparticle *	particle

int	nParticle

int	nPhSpParticles

int	nRandomParticles

int	idCurrentParticleSource

G4String	sourceTypeName

Static Private Attributes
static CML2PrimaryGenerationAction *	instance = 0

Detailed Description

Definition at line 65 of file ML2PrimaryGenerationAction.hh.

Constructor & Destructor Documentation

◆ CML2PrimaryGenerationAction()

CML2PrimaryGenerationAction::CML2PrimaryGenerationAction ( void )

Definition at line 47 of file ML2PrimaryGenerationAction.cc.

 :particleGun(0),gamma(0),electron(0),positron(0),primaryParticleData(0),particles(0)
 {
 }

Here is the caller graph for this function:

◆ ~CML2PrimaryGenerationAction()

CML2PrimaryGenerationAction::~CML2PrimaryGenerationAction ( void )

Definition at line 108 of file ML2PrimaryGenerationAction.cc.

 {
     delete particleGun;
     delete [] particles;
     delete particles;
 }

Member Function Documentation

◆ applySourceRotation()

void CML2PrimaryGenerationAction::applySourceRotation ( )

private

Definition at line 191 of file ML2PrimaryGenerationAction.cc.

 {
     pos=*rm*pos;
     dir=*rm*dir;
 }

Here is the caller graph for this function:

◆ design()

void CML2PrimaryGenerationAction::design ( G4double accTargetZPosition )

Definition at line 80 of file ML2PrimaryGenerationAction.cc.

 {
     accTargetZPosition=aTZ;
     switch (idParticleSource)
     {
     case id_randomTarget:
         setGunRandom();
         break;
     case id_phaseSpace:
         setGunCalculatedPhaseSpace();
         break;
     }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ fillParticlesContainer()

void CML2PrimaryGenerationAction::fillParticlesContainer ( )

private

Definition at line 196 of file ML2PrimaryGenerationAction.cc.

 {
     static int currentFilePosition=0;
     static int currentFileSize=0;
     int startDataFilePosition;
     std::ifstream in;
     in.open(calculatedPhaseSpaceFileIN, std::ios::in);
     if (in)
     {
         std::cout <<"ERROR phase space file: "<< calculatedPhaseSpaceFileIN << " NOT found. Run abort "<< G4endl;
         G4RunManager::GetRunManager()->AbortRun(true);
     }
 
     static bool bFirstTime=true;
     if (bFirstTime)
     {
         in.seekg(-1,std::ios::end);
         currentFileSize=in.tellg();
         in.seekg(0,std::ios::beg);
         bFirstTime=false;
     }
 
     char a[1000];
     in.getline(a,1000,'\n');
     in.getline(a,1000,'\n');
     startDataFilePosition=in.tellg();
     if (currentFilePosition>0)
     {in.seekg(currentFilePosition, std::ios::beg);}
     int i;
     G4double x,y,z;
     G4int d;
 
     static bool checkFileRewind=false;
     static bool bRewindTheFile=false;
     static int nPhSpFileRewind=0;
 
     for (i=0;i<nMaxParticlesInRamPhaseSpace;i++)
     {
         if (bRewindTheFile) // to read the phase space file again to fill the container
         {
             in.close();
             in.open(calculatedPhaseSpaceFileIN, std::ios::in);
             in.seekg(startDataFilePosition, std::ios::beg);
             checkFileRewind=true;
             bRewindTheFile=false;
             std::cout<<"\n################\nI have reached the end of the phase space file "<<++nPhSpFileRewind <<" times, I rewind the file\n" << G4endl;
             std::cout <<"loaded " <<i <<"/"<< nMaxParticlesInRamPhaseSpace<<" particles" << G4endl;
         }
         in >> d; 
         in >> x; in >>y; in >> z; 
 /*          std::cout <<"x:" <<x << G4endl;
             std::cout <<"y:" <<y << G4endl;
             std::cout <<"z:" <<z << G4endl;*/
         particles[i].pos.set(x,y,z-accTargetZPosition);
         in >> x; in >>y; in >> z; 
         particles[i].dir.set(x,y,z);
         in >> x; 
         particles[i].kinEnergy=x;
         in >> d; 
         particles[i].partPDGE=d;
         in >> d; in >> d; 
         if (in.eof())   {bRewindTheFile=true;}
         if (checkFileRewind)    {checkFileRewind=false;}
     }
     std::cout <<"loaded " <<i <<"/"<< nMaxParticlesInRamPhaseSpace<<" particles" << G4endl;
     currentFilePosition=in.tellg(); // to remind the actual position in the phase space file
     if (currentFilePosition>=currentFileSize) // to read the phase space file again
     {currentFilePosition=startDataFilePosition;} 
     in.close();
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ GenerateFromCalculatedPhaseSpace()

void CML2PrimaryGenerationAction::GenerateFromCalculatedPhaseSpace ( )

private

Definition at line 162 of file ML2PrimaryGenerationAction.cc.

 {
     static bool bFirstTime=true;
     if (bFirstTime)
     {bFirstTime=false;fillParticlesContainer();}
     if (nParticle==nMaxParticlesInRamPhaseSpace) // once all the particles stored in RAM hae been processed a new set is loaded
     {
         fillParticlesContainer();
         nParticle=0;
     }
 
     pos=particles[nParticle].pos;
     dir=particles[nParticle].dir;
     ek=particles[nParticle].kinEnergy;
     switch (particles[nParticle].partPDGE)
     {
     case -11:
         particleGun->SetParticleDefinition(positron);
         break;
     case 11:
         particleGun->SetParticleDefinition(electron);
         break;
     case 22:
         particleGun->SetParticleDefinition(gamma);
         break;
     }
     nPhSpParticles++;
     nParticle++;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ GenerateFromRandom()

void CML2PrimaryGenerationAction::GenerateFromRandom ( )

private

Definition at line 146 of file ML2PrimaryGenerationAction.cc.

 {
     sinTheta=RandGauss::shoot(0., 0.003);
     cosTheta=std::sqrt(1 - sinTheta*sinTheta);
     phi=twopi*G4UniformRand();
     dir.set(sinTheta*std::cos(phi), sinTheta*std::sin(phi), cosTheta);
 
     ro=G4UniformRand()*GunRadious; 
     alfa=G4UniformRand()*twopi;
 
     pos.setX(ro*std::sin(alfa));
     pos.setY(ro*std::cos(alfa));
     pos.setZ(-(accTargetZPosition +5.)*mm); // the primary electrons are generated 5 mm before the target
     ek=RandGauss::shoot(GunMeanEnegy, GunStdEnegy);
     nRandomParticles++;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ GeneratePrimaries()

void CML2PrimaryGenerationAction::GeneratePrimaries ( G4Event * anEvent )

virtual

Implements G4VUserPrimaryGeneratorAction.

Definition at line 114 of file ML2PrimaryGenerationAction.cc.

 {
     static int currentRecycle=nRecycling;
     static G4ThreeVector pos0, dir0;
     if (currentRecycle==nRecycling)
     {
         currentRecycle=0;
         switch (idCurrentParticleSource)
         {
         case id_randomTarget:
                 GenerateFromRandom();
             break;
         case id_phaseSpace:
                 GenerateFromCalculatedPhaseSpace();
             break;
         }
         pos0=pos;
         dir0=dir;
     }
     currentRecycle++;
     pos=pos0;
     dir=dir0;
     applySourceRotation(); // to follow the accelerator rotation
 
     primaryParticleData->partPDGE=particleGun->GetParticleDefinition()->GetPDGEncoding();
     primaryParticleData->nPrimaryParticle++;
 
     particleGun->SetParticleEnergy(ek*MeV);
     particleGun->SetParticlePosition(pos*mm);
     particleGun->SetParticleMomentumDirection((G4ParticleMomentum)dir);
     particleGun->GeneratePrimaryVertex(anEvent);
 }

Here is the call graph for this function:

◆ GetInstance()

CML2PrimaryGenerationAction * CML2PrimaryGenerationAction::GetInstance ( void )

static

Definition at line 53 of file ML2PrimaryGenerationAction.cc.

 {
   if (instance == 0)
     {
       instance = new CML2PrimaryGenerationAction();
     }
   return instance;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ getNrecycling()

G4int CML2PrimaryGenerationAction::getNrecycling ( )

inline

Definition at line 96 of file ML2PrimaryGenerationAction.hh.

96 {return nRecycling;};

CML2PrimaryGenerationAction::nRecycling

G4int nRecycling

Definition: ML2PrimaryGenerationAction.hh:109

Here is the caller graph for this function:

◆ getSourceTypeName()

G4int CML2PrimaryGenerationAction::getSourceTypeName ( )

inline

Definition at line 97 of file ML2PrimaryGenerationAction.hh.

97 {return idParticleSource;};

CML2PrimaryGenerationAction::idParticleSource

G4int idParticleSource

Definition: ML2PrimaryGenerationAction.hh:109

Here is the call graph for this function:

◆ inizialize()

void CML2PrimaryGenerationAction::inizialize ( SPrimaryParticle * primaryParticleData )

Definition at line 61 of file ML2PrimaryGenerationAction.cc.

 {
     rm=new G4RotationMatrix();
     PrimaryGenerationActionMessenger=new CML2PrimaryGenerationActionMessenger(this);
     particle=new Sparticle;
     nParticle=nPhSpParticles=nRandomParticles=0;
 
     G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();
 
     gamma=particleTable->FindParticle("gamma");
     electron=particleTable->FindParticle("e-");
     positron=particleTable->FindParticle("e+");
     particleGun=new G4ParticleGun();
 
     primaryParticleData=pData;
     primaryParticleData->nPrimaryParticle=0;
     primaryParticleData->partPDGE=0;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ setCalculatedPhaseSpaceFileIN()

void CML2PrimaryGenerationAction::setCalculatedPhaseSpaceFileIN ( G4String val )

inline

Definition at line 82 of file ML2PrimaryGenerationAction.hh.

82 {calculatedPhaseSpaceFileIN=val;}

CML2PrimaryGenerationAction::calculatedPhaseSpaceFileIN

G4String calculatedPhaseSpaceFileIN

Definition: ML2PrimaryGenerationAction.hh:111

Here is the caller graph for this function:

◆ setGunCalculatedPhaseSpace()

void CML2PrimaryGenerationAction::setGunCalculatedPhaseSpace ( )

private

Definition at line 101 of file ML2PrimaryGenerationAction.cc.

 {
     particles=new Sparticle[nMaxParticlesInRamPhaseSpace];
     particleGun->SetNumberOfParticles(1);
     idCurrentParticleSource=idParticleSource;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ setGunMeanEnergy()

void CML2PrimaryGenerationAction::setGunMeanEnergy ( G4double val )

inline

Definition at line 79 of file ML2PrimaryGenerationAction.hh.

79 {GunMeanEnegy=val;}

CML2PrimaryGenerationAction::GunMeanEnegy

G4double GunMeanEnegy

Definition: ML2PrimaryGenerationAction.hh:110

Here is the caller graph for this function:

◆ setGunRadious()

void CML2PrimaryGenerationAction::setGunRadious ( G4double val )

inline

Definition at line 81 of file ML2PrimaryGenerationAction.hh.

81 {GunRadious=val;}

CML2PrimaryGenerationAction::GunRadious

G4double GunRadious

Definition: ML2PrimaryGenerationAction.hh:110

Here is the caller graph for this function:

◆ setGunRandom()

void CML2PrimaryGenerationAction::setGunRandom ( )

private

Definition at line 94 of file ML2PrimaryGenerationAction.cc.

 {
     particleGun->SetParticleDefinition(electron);
     particleGun->SetNumberOfParticles(1);
     idCurrentParticleSource=idParticleSource;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ setGunStdEnergy()

void CML2PrimaryGenerationAction::setGunStdEnergy ( G4double val )

inline

Definition at line 80 of file ML2PrimaryGenerationAction.hh.

80 {GunStdEnegy=val;}

CML2PrimaryGenerationAction::GunStdEnegy

G4double GunStdEnegy

Definition: ML2PrimaryGenerationAction.hh:110

Here is the caller graph for this function:

◆ setNLoopsPhSpParticles()

void CML2PrimaryGenerationAction::setNLoopsPhSpParticles ( G4int val )

inline

Definition at line 76 of file ML2PrimaryGenerationAction.hh.

76 {nLoopsPhSpParticles=val;}

CML2PrimaryGenerationAction::nLoopsPhSpParticles

G4int nLoopsPhSpParticles

Definition: ML2PrimaryGenerationAction.hh:109

◆ setNMaxParticlesInRamPhaseSpace()

void CML2PrimaryGenerationAction::setNMaxParticlesInRamPhaseSpace ( G4int val )

inline

Definition at line 77 of file ML2PrimaryGenerationAction.hh.

77 {nMaxParticlesInRamPhaseSpace=val;std::cout<<"Current nMaxParticlesInRamPhaseSpace: " << nMaxParticlesInRamPhaseSpace<< G4endl;}

CML2PrimaryGenerationAction::nMaxParticlesInRamPhaseSpace

G4int nMaxParticlesInRamPhaseSpace

Definition: ML2PrimaryGenerationAction.hh:109

G4endl

#define G4endl

Definition: G4ios.hh:61

Here is the caller graph for this function:

◆ setNRecycling()

void CML2PrimaryGenerationAction::setNRecycling ( G4int val )

inline

Definition at line 75 of file ML2PrimaryGenerationAction.hh.

75 {nRecycling=val;}

CML2PrimaryGenerationAction::nRecycling

G4int nRecycling

Definition: ML2PrimaryGenerationAction.hh:109

Here is the caller graph for this function:

◆ setRotation()

void CML2PrimaryGenerationAction::setRotation ( G4RotationMatrix * val )

inline

Definition at line 95 of file ML2PrimaryGenerationAction.hh.

95 {rm=val;};

CML2PrimaryGenerationAction::rm

G4RotationMatrix * rm

Definition: ML2PrimaryGenerationAction.hh:118

Here is the caller graph for this function:

◆ setSourceTypeName()

void CML2PrimaryGenerationAction::setSourceTypeName ( G4String val )

inline

Definition at line 83 of file ML2PrimaryGenerationAction.hh.

     {
         sourceTypeName=val;
         if (sourceTypeName=="randomTarget")
         {
             idParticleSource=id_randomTarget;
         }
         else if (sourceTypeName=="phaseSpace")
         {
             idParticleSource=id_phaseSpace;
         }
     }