CCalPrimaryGeneratorAction Class Reference

#include <CCalPrimaryGeneratorAction.hh>

Inheritance 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 48 of file CCalPrimaryGeneratorAction.cc.

References python.hepunit::degree, G4ParticleTable::FindParticle(), G4ParticleTable::GetParticleTable(), python.hepunit::GeV, SetMaximumEnergy(), SetMaximumEta(), SetMaximumPhi(), SetMinimumEnergy(), SetMinimumEta(), SetMinimumPhi(), G4ParticleGun::SetParticleDefinition(), G4ParticleGun::SetParticleEnergy(), G4ParticleGun::SetParticleMomentumDirection(), G4VPrimaryGenerator::SetParticlePosition(), SetStepsEta(), SetStepsPhi(), and python.hepunit::TeV.

                                                      : 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);
}

CCalPrimaryGeneratorAction::~CCalPrimaryGeneratorAction

(

)

Definition at line 81 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 88 of file CCalPrimaryGeneratorAction.cc.

References python.hepunit::deg, G4cout, G4endl, G4ParticleGun::GeneratePrimaryVertex(), python.hepunit::GeV, G4ParticleGun::SetParticleEnergy(), G4ParticleGun::SetParticleMomentumDirection(), CLHEP::RandFlat::shoot(), singleRandom, and singleScan.

                                                                   {

  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);
} 

G4double CCalPrimaryGeneratorAction::GetParticleEnergy ( )

inline

Definition at line 91 of file CCalPrimaryGeneratorAction.hh.

References G4ParticleGun::GetParticleEnergy().

Referenced by CCalEndOfEventAction::EndOfEventAction().

{return particleGun->GetParticleEnergy();}

G4ThreeVector CCalPrimaryGeneratorAction::GetParticlePosition ( )

inline

Definition at line 90 of file CCalPrimaryGeneratorAction.hh.

References G4VPrimaryGenerator::GetParticlePosition().

Referenced by CCalEndOfEventAction::EndOfEventAction().

{return particleGun->GetParticlePosition();}

void CCalPrimaryGeneratorAction::SetGunPosition

(

const G4ThreeVector &

pos

)

const

Definition at line 327 of file CCalPrimaryGeneratorAction.cc.

References G4VPrimaryGenerator::SetParticlePosition().

                                                                               {

  particleGun->SetParticlePosition(pos);
}

void CCalPrimaryGeneratorAction::SetMaximumEnergy

(

G4double

p

)

Definition at line 215 of file CCalPrimaryGeneratorAction.cc.

References G4cerr, G4cout, G4endl, and python.hepunit::GeV.

Referenced by CCalPrimaryGeneratorAction(), and CCalPrimaryGeneratorMessenger::SetNewValue().

                                                           {

  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
  }
}

void CCalPrimaryGeneratorAction::SetMaximumEta

(

G4double

p

)

Definition at line 299 of file CCalPrimaryGeneratorAction.cc.

References G4cout, and G4endl.

Referenced by CCalPrimaryGeneratorAction(), and CCalPrimaryGeneratorMessenger::SetNewValue().

                                                        {

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

void CCalPrimaryGeneratorAction::SetMaximumPhi

(

G4double

p

)

Definition at line 251 of file CCalPrimaryGeneratorAction.cc.

References G4cerr, G4cout, G4endl, and python.hepunit::pi.

Referenced by CCalPrimaryGeneratorAction(), and CCalPrimaryGeneratorMessenger::SetNewValue().

                                                        {

  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
  }
}

void CCalPrimaryGeneratorAction::SetMinimumEnergy

(

G4double

p

)

Definition at line 197 of file CCalPrimaryGeneratorAction.cc.

References G4cerr, G4cout, G4endl, and python.hepunit::GeV.

Referenced by CCalPrimaryGeneratorAction(), and CCalPrimaryGeneratorMessenger::SetNewValue().

                                                           {

  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
  }
}

void CCalPrimaryGeneratorAction::SetMinimumEta

(

G4double

p

)

Definition at line 287 of file CCalPrimaryGeneratorAction.cc.

References G4cout, and G4endl.

Referenced by CCalPrimaryGeneratorAction(), and CCalPrimaryGeneratorMessenger::SetNewValue().

                                                        {

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

void CCalPrimaryGeneratorAction::SetMinimumPhi

(

G4double

p

)

Definition at line 233 of file CCalPrimaryGeneratorAction.cc.

References G4cerr, G4cout, G4endl, and python.hepunit::pi.

Referenced by CCalPrimaryGeneratorAction(), and CCalPrimaryGeneratorMessenger::SetNewValue().

                                                        {

  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
  }
}

void CCalPrimaryGeneratorAction::SetRandom

(

G4String

val

)

Definition at line 172 of file CCalPrimaryGeneratorAction.cc.

References singleFixed, and singleRandom.

Referenced by CCalPrimaryGeneratorMessenger::SetNewValue().

                                                       { 

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

void CCalPrimaryGeneratorAction::SetRunNo

(

G4int

val

)

Definition at line 332 of file CCalPrimaryGeneratorAction.cc.

References G4RunManager::GetRunManager(), and G4RunManager::SetRunIDCounter().

Referenced by CCalPrimaryGeneratorMessenger::SetNewValue().

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

void CCalPrimaryGeneratorAction::SetScan

(

G4String

val

)

Definition at line 184 of file CCalPrimaryGeneratorAction.cc.

References singleFixed, and singleScan.

Referenced by CCalPrimaryGeneratorMessenger::SetNewValue().

                                                     { 

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

void CCalPrimaryGeneratorAction::SetStepsEta

(

G4int

val

)

Definition at line 311 of file CCalPrimaryGeneratorAction.cc.

References G4cerr, G4cout, and G4endl.

Referenced by CCalPrimaryGeneratorAction(), and CCalPrimaryGeneratorMessenger::SetNewValue().

                                                     {

  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
  }
}

void CCalPrimaryGeneratorAction::SetStepsPhi

(

G4int

val

)

Definition at line 269 of file CCalPrimaryGeneratorAction.cc.

References G4cerr, G4cout, and G4endl.

Referenced by CCalPrimaryGeneratorAction(), and CCalPrimaryGeneratorMessenger::SetNewValue().

                                                     {

  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
  }
}

void CCalPrimaryGeneratorAction::SetVerboseLevel

(

G4int

val

)

Definition at line 167 of file CCalPrimaryGeneratorAction.cc.

Referenced by CCalPrimaryGeneratorMessenger::SetNewValue().

                                                         {
  verboseLevel = val;
}

---

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

• CCalPrimaryGeneratorAction.hh
• CCalPrimaryGeneratorAction.cc