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 ()

Detailed Description

Definition at line 66 of file CCalPrimaryGeneratorAction.hh.

Constructor & Destructor Documentation

CCalPrimaryGeneratorAction::CCalPrimaryGeneratorAction

(

)

Definition at line 49 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

(

)

Definition at line 82 of file CCalPrimaryGeneratorAction.cc.

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

Member Function Documentation

void CCalPrimaryGeneratorAction::GeneratePrimaries ( G4Event *  anEvent )

virtual

Implements G4VUserPrimaryGeneratorAction.

Definition at line 89 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(G4Exp(-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(G4Exp(-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:

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:

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:

void CCalPrimaryGeneratorAction::SetGunPosition

(

const G4ThreeVector &

pos

)

const

Definition at line 328 of file CCalPrimaryGeneratorAction.cc.

                                                                               {

  particleGun->SetParticlePosition(pos);
}

Here is the call graph for this function:

void CCalPrimaryGeneratorAction::SetMaximumEnergy

(

G4double

p

)

Definition at line 216 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:

void CCalPrimaryGeneratorAction::SetMaximumEta

(

G4double

p

)

Definition at line 300 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:

void CCalPrimaryGeneratorAction::SetMaximumPhi

(

G4double

p

)

Definition at line 252 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:

void CCalPrimaryGeneratorAction::SetMinimumEnergy

(

G4double

p

)

Definition at line 198 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:

void CCalPrimaryGeneratorAction::SetMinimumEta

(

G4double

p

)

Definition at line 288 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:

void CCalPrimaryGeneratorAction::SetMinimumPhi

(

G4double

p

)

Definition at line 234 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:

void CCalPrimaryGeneratorAction::SetRandom

(

G4String

val

)

Definition at line 173 of file CCalPrimaryGeneratorAction.cc.

                                                       { 

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

Here is the caller graph for this function:

void CCalPrimaryGeneratorAction::SetRunNo

(

G4int

val

)

Definition at line 333 of file CCalPrimaryGeneratorAction.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

void CCalPrimaryGeneratorAction::SetScan

(

G4String

val

)

Definition at line 185 of file CCalPrimaryGeneratorAction.cc.

                                                     { 

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

Here is the caller graph for this function:

void CCalPrimaryGeneratorAction::SetStepsEta

(

G4int

val

)

Definition at line 312 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:

void CCalPrimaryGeneratorAction::SetStepsPhi

(

G4int

val

)

Definition at line 270 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:

void CCalPrimaryGeneratorAction::SetVerboseLevel

(

G4int

val

)

Definition at line 168 of file CCalPrimaryGeneratorAction.cc.

                                                         {
  verboseLevel = val;
}

Here is the caller graph for this function:

---

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

• source/geant4.10.03.p02/examples/advanced/composite_calorimeter/include/CCalPrimaryGeneratorAction.hh
• source/geant4.10.03.p02/examples/advanced/composite_calorimeter/src/CCalPrimaryGeneratorAction.cc