CCalPrimaryGeneratorAction Class Reference

#include <CCalPrimaryGeneratorAction.hh>

Inheritance diagram for CCalPrimaryGeneratorAction:

Collaboration diagram for CCalPrimaryGeneratorAction:

Public Member Functions
	CCalPrimaryGeneratorAction ()
	~CCalPrimaryGeneratorAction ()
void	GeneratePrimaries (G4Event *anEvent)
void	SetVerboseLevel (G4int val)
void	SetRandom (G4String val)
void	SetScan (G4String val)
void	SetMinimumEnergy (G4double p)
void	SetMaximumEnergy (G4double p)
void	SetMinimumPhi (G4double p)
void	SetMaximumPhi (G4double p)
void	SetStepsPhi (G4int val)
void	SetMinimumEta (G4double p)
void	SetMaximumEta (G4double p)
void	SetStepsEta (G4int val)
void	SetGunPosition (const G4ThreeVector &pos) const
void	SetRunNo (G4int val)
G4ThreeVector	GetParticlePosition ()
G4double	GetParticleEnergy ()
Public Member Functions inherited from G4VUserPrimaryGeneratorAction
	G4VUserPrimaryGeneratorAction ()
virtual	~G4VUserPrimaryGeneratorAction ()

Private Member Functions
void	initialize ()
void	print (G4int val)

Private Attributes
CCalPrimaryGeneratorMessenger *	gunMessenger
G4ParticleGun *	particleGun
generatorInputType	generatorInput
G4int	verboseLevel
G4int	n_particle
G4String	particleName
G4double	particleEnergy
G4ThreeVector	particlePosition
G4ThreeVector	particleDir
G4double	energyMin
G4double	energyMax
G4double	etaMin
G4double	etaMax
G4double	phiMin
G4double	phiMax
G4int	etaSteps
G4int	phiSteps
G4int	isInitialized
G4double	etaValue
G4double	phiValue
G4int	scanSteps

Detailed Description

Definition at line 66 of file CCalPrimaryGeneratorAction.hh.

Constructor & Destructor Documentation

CCalPrimaryGeneratorAction()

CCalPrimaryGeneratorAction::CCalPrimaryGeneratorAction

(

)

Definition at line 48 of file CCalPrimaryGeneratorAction.cc.

                                                      : particleGun(0),
  generatorInput(singleFixed),  verboseLevel(0), n_particle(1), 
  particleName("pi-"), particleEnergy(100*GeV), particlePosition(0.,0.,0.),
  particleDir(1.,1.,0.1), isInitialized(0), scanSteps(0) {
  
  //Initialise the messenger
  gunMessenger = new CCalPrimaryGeneratorMessenger(this);
    
  //Default settings:
  SetMinimumEnergy(1.*GeV);
  SetMaximumEnergy(1.*TeV);
  SetMinimumEta(0.);
  SetMaximumEta(3.5);
  SetMinimumPhi(0.*degree);
  SetMaximumPhi(360.*degree);
  SetStepsPhi(1);
  SetStepsEta(1);
    
  // Number of particles per gun
  particleGun = new G4ParticleGun(n_particle);
    
  // Getting particle definition
  G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();
  G4ParticleDefinition* particle = particleTable->FindParticle(particleName);
  particleGun->SetParticleDefinition(particle);
    
  // Setting particle properties
  particleGun->SetParticleEnergy(particleEnergy);
  particleGun->SetParticleMomentumDirection(particleDir);
  particleGun->SetParticlePosition(particlePosition);
  print(0);
}

Here is the call graph for this function:

~CCalPrimaryGeneratorAction()

CCalPrimaryGeneratorAction::~CCalPrimaryGeneratorAction

(

)

Definition at line 81 of file CCalPrimaryGeneratorAction.cc.

                                                        {
  if (gunMessenger)
    delete gunMessenger;
  if (particleGun)
    delete particleGun;
}

Member Function Documentation

GeneratePrimaries()

void CCalPrimaryGeneratorAction::GeneratePrimaries ( G4Event *  anEvent )

virtual

Implements G4VUserPrimaryGeneratorAction.

Definition at line 88 of file CCalPrimaryGeneratorAction.cc.

                                                                   {

  if (isInitialized == 0) initialize();

  if (generatorInput == singleRandom) {
    particleEnergy = CLHEP::RandFlat::shoot(energyMin,energyMax);
    particleGun->SetParticleEnergy(particleEnergy);

    G4double eta = CLHEP::RandFlat::shoot(etaMin,etaMax);
    G4double phi = CLHEP::RandFlat::shoot(phiMin,phiMax);
    G4double theta = std::atan(std::exp(-eta))*2.;
    G4double randomX = std::sin(theta)*std::cos(phi);
    G4double randomY = std::sin(theta)*std::sin(phi);
    G4double randomZ = std::cos(theta);
  
    particleDir = G4ThreeVector(randomX,randomY,randomZ);
    particleGun->SetParticleMomentumDirection(particleDir);
#ifdef debug
    if (verboseLevel >= 2 ) {
      G4cout << "Energy " << particleEnergy/GeV << " GeV; Theta " 
         << theta/deg << " degree; Phi " << phi/deg << " degree" << G4endl;
      G4cout << "Shooting in " << particleDir << " direction "<< G4endl;
    }
#endif
  } else if (generatorInput == singleScan) {
    G4double scanEtaStep, scanPhiStep;
    if (scanSteps == 0) { 
      scanPhiStep = (phiMax - phiMin) / phiSteps;
      phiValue = phiMin - scanPhiStep; //first step is going to change scanCurrentPhiValue
      etaValue = etaMin;
    }

    scanEtaStep = (etaMax - etaMin) / etaSteps;
    scanPhiStep = (phiMax - phiMin) / phiSteps;
#ifdef debug
    if (verboseLevel > 2 ) {
      G4cout << " scanEtaStep " << scanEtaStep << " # of Steps " << etaSteps 
         << G4endl;
      G4cout << " scanPhiStep " << scanPhiStep << " # of Steps " << phiSteps 
         << G4endl;
    }
#endif

    //----- First scan in phi, then in eta
    if (phiMax - phiValue < 1.E-6 * scanPhiStep) { // !only <= 1.E6 steps allowed
      if (etaMax - etaValue < 1.E-6 * scanEtaStep) { // !only <= 1.E6 steps allowed
    G4cout << " Scan completed!" << G4endl;
    return;
      } else {
    etaValue += scanEtaStep; 
    phiValue  = phiMin;
      }
    } else {
      phiValue += scanPhiStep;
    }    
    G4double theta = std::atan(std::exp(-etaValue))*2.;

    G4double scanX = std::sin(theta)*std::cos(phiValue);
    G4double scanY = std::sin(theta)*std::sin(phiValue);
    G4double scanZ = std::cos(theta);
    if (verboseLevel >= 2 ) {
      G4cout << "Scan eta " << etaValue << " Phi " << phiValue/deg
         << " theta " << theta/deg << G4endl;
    }
    particleDir = G4ThreeVector(scanX,scanY,scanZ);
    particleGun->SetParticleMomentumDirection(particleDir);
#ifdef debug
    if (verboseLevel > 2 ) {
      G4cout  << "Shooting in " << particleDir << " direction "<< G4endl;
    }
#endif
    scanSteps++;
  }
  
  // Generate GEANT4 primary vertex
  particleGun->GeneratePrimaryVertex(anEvent);
} 

Here is the call graph for this function:

GetParticleEnergy()

G4double CCalPrimaryGeneratorAction::GetParticleEnergy ( )

inline

Definition at line 91 of file CCalPrimaryGeneratorAction.hh.

{return particleGun->GetParticleEnergy();}

Here is the call graph for this function:

Here is the caller graph for this function:

GetParticlePosition()

G4ThreeVector CCalPrimaryGeneratorAction::GetParticlePosition ( )

inline

Definition at line 90 of file CCalPrimaryGeneratorAction.hh.

{return particleGun->GetParticlePosition();}

Here is the call graph for this function:

Here is the caller graph for this function:

initialize()

void CCalPrimaryGeneratorAction::initialize ( )

private

Definition at line 336 of file CCalPrimaryGeneratorAction.cc.

                                           {

  isInitialized = 1;

  print (1);
}

Here is the call graph for this function:

Here is the caller graph for this function:

print()

void CCalPrimaryGeneratorAction::print ( G4int  val )

private

Definition at line 344 of file CCalPrimaryGeneratorAction.cc.

                                               {

#ifdef debug
  if (verboseLevel >= val) {

    if (generatorInput == singleRandom) {
      G4cout << G4endl
             << "**********************************************************************" << G4endl
         << "*                                                                    *" << G4endl  
             << "* CCalPrimaryGeneratorAction DEFAULT Random Energy/Direction setting:*" << G4endl
         << "*                                                                    *" << G4endl  
         << "*                                                                    *" << G4endl  
         << "*   Energy in    [ "<< energyMin/GeV   << " - " << energyMax/GeV   << "] (GeV) "<< G4endl
         << "*   Phi angle in [ "<< phiMin          << " - " << phiMax          << "] (rad) "<< G4endl
         << "*                [ "<< phiMin/degree   << " - " << phiMax/degree   << "] (deg) "<< G4endl 
         << "*   Eta in       [ "<< etaMin          << " - " << etaMax          << "]       "<< G4endl
         << "*                                                                    *" << G4endl  
         << "*                                                                    *" << G4endl  
             << "**********************************************************************" << G4endl;
    } else if (generatorInput == singleScan) {
      G4cout << G4endl
         << "**********************************************************************" << G4endl
         << "*                                                                    *" << G4endl  
             << "* CCalPrimaryGeneratorAction DEFAULT Scan Direction settings :       *" << G4endl
         << "*                                                                    *" << G4endl  
         << "*                                                                    *" << G4endl  
         << "*   Phi angle in [ " << phiMin/degree   << " - " << phiMax/degree << "] (deg) " << G4endl
         << "*   Eta in       [ " << etaMin          << " - " << etaMax        << "]       " << G4endl
         << "*   Steps along eta " << etaSteps << " and along phi " << phiSteps << G4endl
         << "*                                                                    *" << G4endl  
         << "*                                                                    *" << G4endl  
             << "**********************************************************************" << G4endl;
    } else  if (generatorInput == singleFixed) {
      G4cout << G4endl
         << "*******************************************************************" << G4endl
         << "*                                                                 *" << G4endl  
         << "* CCalPrimaryGeneratorAction: Current settings :                  *" << G4endl
         << "*                                                                 *" << G4endl  
         << "* " << particleGun->GetNumberOfParticles() 
         << "  " << particleGun->GetParticleDefinition()->GetParticleName() 
         << " of " << particleGun->GetParticleEnergy()/GeV << " GeV" << G4endl
         << "* will be shot from " << particleGun->GetParticlePosition() << G4endl;
      G4cout << "* in direction " << particleGun->GetParticleMomentumDirection() << G4endl;
      G4cout  << "*                                                                 *" << G4endl  
          << "*                                                                 *" << G4endl  
          << "*******************************************************************" << G4endl;
    }
  }
#endif
 
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetGunPosition()

void CCalPrimaryGeneratorAction::SetGunPosition

(

const G4ThreeVector &

pos

)

const

Definition at line 327 of file CCalPrimaryGeneratorAction.cc.

                                                                               {

  particleGun->SetParticlePosition(pos);
}

Here is the call graph for this function:

SetMaximumEnergy()

void CCalPrimaryGeneratorAction::SetMaximumEnergy

(

G4double

p

)

Definition at line 215 of file CCalPrimaryGeneratorAction.cc.

                                                           {

  if (p <= 0.) {
    G4cerr << "CCalPrimaryGeneratorAction::SetMaximumEnergy: value " << p/GeV 
       << "GeV is out of bounds, it will not be used" << G4endl;
    G4cerr << " Should be  >0.  Please check" << G4endl; 
  } else {
    energyMax = p;
#ifdef debug
    if (verboseLevel >= 1 ) {
      G4cout << " CCalPrimaryGeneratorAction: setting max. value of energy to "
         << energyMax/GeV << " GeV " << G4endl;
    }
#endif
  }
}

Here is the caller graph for this function:

SetMaximumEta()

void CCalPrimaryGeneratorAction::SetMaximumEta

(

G4double

p

)

Definition at line 299 of file CCalPrimaryGeneratorAction.cc.

                                                        {

  etaMax = p;
#ifdef debug
  if (verboseLevel >= 1 ) {
    G4cout << " CCalPrimaryGeneratorAction: setting max. value of eta to "
       << etaMax << G4endl;
  }
#endif
}

Here is the caller graph for this function:

SetMaximumPhi()

void CCalPrimaryGeneratorAction::SetMaximumPhi

(

G4double

p

)

Definition at line 251 of file CCalPrimaryGeneratorAction.cc.

                                                        {

  if (std::fabs(p)>2.*pi) {
    G4cerr << "CCalPrimaryGeneratorAction::SetMaximumPhi: setting value quite "
       << "large" << G4endl;
    G4cerr << " Should be given in radians - Please check" << G4endl;
  } else {
    phiMax = std::fabs(p);
#ifdef debug
    if (verboseLevel >= 1 ) {
      G4cout << " CCalPrimaryGeneratorAction: setting max. value of phi to "
         << phiMax << G4endl;
    }
#endif
  }
}

Here is the caller graph for this function:

SetMinimumEnergy()

void CCalPrimaryGeneratorAction::SetMinimumEnergy

(

G4double

p

)

Definition at line 197 of file CCalPrimaryGeneratorAction.cc.

                                                           {

  if (p <= 0.) {
    G4cerr << "CCalPrimaryGeneratorAction::SetMinimumEnergy: value " << p/GeV 
       << "GeV is out of bounds, it will not be used" << G4endl;
    G4cerr << " Should be  >0.  Please check" << G4endl; 
  } else {
    energyMin = p;
#ifdef debug
    if (verboseLevel >= 1 ) {
      G4cout << " CCalPrimaryGeneratorAction: setting min. value of energy to "
         << energyMin/GeV << " GeV " << G4endl;
    }
#endif
  }
}

Here is the caller graph for this function:

SetMinimumEta()

void CCalPrimaryGeneratorAction::SetMinimumEta

(

G4double

p

)

Definition at line 287 of file CCalPrimaryGeneratorAction.cc.

                                                        {

  etaMin = p;
#ifdef debug
  if (verboseLevel >= 1 ) {
    G4cout << " CCalPrimaryGeneratorAction: setting min. value of eta to "
       << etaMin << G4endl;
  }
#endif
}

Here is the caller graph for this function:

SetMinimumPhi()

void CCalPrimaryGeneratorAction::SetMinimumPhi

(

G4double

p

)

Definition at line 233 of file CCalPrimaryGeneratorAction.cc.

                                                        {

  if (std::fabs(p)>2.*pi) {
    G4cerr << "CCalPrimaryGeneratorAction::SetMinimumPhi: setting value quite "
       << "large" << G4endl;
    G4cerr << " Should be given in radians - Please check" << G4endl;
  } else {
    phiMin = std::fabs(p);
#ifdef debug
    if (verboseLevel >= 1 ) {
      G4cout << " CCalPrimaryGeneratorAction: setting min. value of phi to "
         << phiMin << G4endl;
    }
#endif
  }
}

Here is the caller graph for this function:

SetRandom()

void CCalPrimaryGeneratorAction::SetRandom

(

G4String

val

)

Definition at line 172 of file CCalPrimaryGeneratorAction.cc.

                                                       { 

  if (val=="on") {
    generatorInput = singleRandom;
    print (1);
  } else {
    generatorInput = singleFixed;
    print (1);
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetRunNo()

void CCalPrimaryGeneratorAction::SetRunNo

(

G4int

val

)

Definition at line 332 of file CCalPrimaryGeneratorAction.cc.

                                                  {
  G4RunManager::GetRunManager()->SetRunIDCounter( val );
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetScan()

void CCalPrimaryGeneratorAction::SetScan

(

G4String

val

)

Definition at line 184 of file CCalPrimaryGeneratorAction.cc.

                                                     { 

  if (val=="on") {
    generatorInput = singleScan;
    scanSteps = 0;
    print (1);
  } else {
    generatorInput = singleFixed;
    print (1);
  }  
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetStepsEta()

void CCalPrimaryGeneratorAction::SetStepsEta

(

G4int

val

)

Definition at line 311 of file CCalPrimaryGeneratorAction.cc.

                                                     {

  if (val <= 0) {
    G4cerr<<"CCalPrimaryGeneratorAction::SetStepsEta: value " << val << " is out of bounds, it will not be used"<<G4endl;
    G4cerr<<" Should be  > 0  Please check"<<G4endl; 
  } else {
    etaSteps = val;
#ifdef debug
    if (verboseLevel >= 1 ) {
      G4cout << " CCalPrimaryGeneratorAction: setting no. of steps in eta to "
         << etaSteps << G4endl;
    }
#endif
  }
}

Here is the caller graph for this function:

SetStepsPhi()

void CCalPrimaryGeneratorAction::SetStepsPhi

(

G4int

val

)

Definition at line 269 of file CCalPrimaryGeneratorAction.cc.

                                                     {

  if (val <= 0) {
    G4cerr << "CCalPrimaryGeneratorAction::SetStepsPhi: value " << val 
       << " is out of bounds, it will not be used" << G4endl;
    G4cerr << " Should be  > 0  Please check" << G4endl; 
  } else {
    phiSteps = val;
#ifdef debug
    if (verboseLevel >= 1 ) {
      G4cout << " CCalPrimaryGeneratorAction: setting no. of steps in phi to "
         << phiSteps << G4endl;
    }
#endif
  }
}

Here is the caller graph for this function:

SetVerboseLevel()

void CCalPrimaryGeneratorAction::SetVerboseLevel

(

G4int

val

)

Definition at line 167 of file CCalPrimaryGeneratorAction.cc.

                                                         {
  verboseLevel = val;
}

Here is the caller graph for this function:

Member Data Documentation

energyMax

G4double CCalPrimaryGeneratorAction::energyMax

private

Definition at line 105 of file CCalPrimaryGeneratorAction.hh.

energyMin

G4double CCalPrimaryGeneratorAction::energyMin

private

Definition at line 105 of file CCalPrimaryGeneratorAction.hh.

etaMax

G4double CCalPrimaryGeneratorAction::etaMax

private

Definition at line 106 of file CCalPrimaryGeneratorAction.hh.

etaMin

G4double CCalPrimaryGeneratorAction::etaMin

private

Definition at line 106 of file CCalPrimaryGeneratorAction.hh.

etaSteps

G4int CCalPrimaryGeneratorAction::etaSteps

private

Definition at line 108 of file CCalPrimaryGeneratorAction.hh.

etaValue

G4double CCalPrimaryGeneratorAction::etaValue

private

Definition at line 111 of file CCalPrimaryGeneratorAction.hh.

generatorInput

generatorInputType CCalPrimaryGeneratorAction::generatorInput

private

Definition at line 96 of file CCalPrimaryGeneratorAction.hh.

gunMessenger

CCalPrimaryGeneratorMessenger* CCalPrimaryGeneratorAction::gunMessenger

private

Definition at line 94 of file CCalPrimaryGeneratorAction.hh.

isInitialized

G4int CCalPrimaryGeneratorAction::isInitialized

private

Definition at line 110 of file CCalPrimaryGeneratorAction.hh.

n_particle

G4int CCalPrimaryGeneratorAction::n_particle

private

Definition at line 99 of file CCalPrimaryGeneratorAction.hh.

particleDir

G4ThreeVector CCalPrimaryGeneratorAction::particleDir

private

Definition at line 103 of file CCalPrimaryGeneratorAction.hh.

particleEnergy

G4double CCalPrimaryGeneratorAction::particleEnergy

private

Definition at line 101 of file CCalPrimaryGeneratorAction.hh.

particleGun

G4ParticleGun* CCalPrimaryGeneratorAction::particleGun

private

Definition at line 95 of file CCalPrimaryGeneratorAction.hh.

particleName

G4String CCalPrimaryGeneratorAction::particleName

private

Definition at line 100 of file CCalPrimaryGeneratorAction.hh.

particlePosition

G4ThreeVector CCalPrimaryGeneratorAction::particlePosition

private

Definition at line 102 of file CCalPrimaryGeneratorAction.hh.

phiMax

G4double CCalPrimaryGeneratorAction::phiMax

private

Definition at line 107 of file CCalPrimaryGeneratorAction.hh.

phiMin

G4double CCalPrimaryGeneratorAction::phiMin

private

Definition at line 107 of file CCalPrimaryGeneratorAction.hh.

phiSteps

G4int CCalPrimaryGeneratorAction::phiSteps

private

Definition at line 108 of file CCalPrimaryGeneratorAction.hh.

phiValue

G4double CCalPrimaryGeneratorAction::phiValue

private

Definition at line 111 of file CCalPrimaryGeneratorAction.hh.

scanSteps

G4int CCalPrimaryGeneratorAction::scanSteps

private

Definition at line 112 of file CCalPrimaryGeneratorAction.hh.

verboseLevel

G4int CCalPrimaryGeneratorAction::verboseLevel

private

Definition at line 98 of file CCalPrimaryGeneratorAction.hh.

---

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

• CCalPrimaryGeneratorAction.hh
• CCalPrimaryGeneratorAction.cc