PrimaryGeneratorAction.cc

Go to the documentation of this file.

#include "PrimaryGeneratorAction.hh"

#include "DetectorConstruction.hh"
#include "PrimaryGeneratorMessenger.hh"


#include "G4ios.hh"
#include "PrimaryGeneratorAction.hh"

void PrimaryGeneratorAction::SetParticle(G4int Z, G4int A, G4double E, G4double Q)
{
  G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();
  particleTable->GetIonTable()->CreateAllIon();

  G4ParticleDefinition* particle = particleTable->GetIon(Z, A, E);
  if(particle!=NULL)
    {
      particleGun->SetParticleCharge(Q);                                        
      particleGun->SetParticleDefinition(particle);
      myRunAction->SetParticleZ(Z);
      myRunAction->SetParticleA(A);
    }
}


PrimaryGeneratorAction::PrimaryGeneratorAction(DetectorConstruction* DC, RunAction* run)
  :Detector(DC),myRunAction(run),rndmFlag("on")
{
  G4int n_particle = 1;
  particleGun  = new G4ParticleGun(n_particle);
  
  //create a messenger for this class
  gunMessenger = new PrimaryGeneratorMessenger(this);

  // default particle kinematic
  particleGun->SetParticleMomentumDirection(G4ThreeVector(1.,0.,0.));
  particleGun->SetParticleEnergy(11*GeV);                       // old value: 2.2 GeV

  //default particle gun position.
  G4double position = -(Detector->GetWorldSizeX());
  particleGun->SetParticlePosition(G4ThreeVector(position,0.*cm,0.*cm));
}

PrimaryGeneratorAction::~PrimaryGeneratorAction()
{
  delete particleGun;
  delete gunMessenger;
}

void PrimaryGeneratorAction::GeneratePrimaries(G4Event* anEvent)
{
  //this function is called at the begining of event
  // 
  
  //Always place the gun outside the world.
  G4double x0 = -0.45*(Detector->GetWorldSizeX());
  G4double y0 = 0.*cm, z0 = 0.*cm;
  if (rndmFlag == "on")
     {
       y0 = (Detector->GetBackScintillatorPlateSizeY())*(G4UniformRand()-0.5);
       z0 = (Detector->GetScintillatorBoxSizeZ())*1*(G4UniformRand()-0.5);
     } 

  particleGun->SetParticlePosition(G4ThreeVector(x0,y0,z0));

  particleGun->GeneratePrimaryVertex(anEvent);

    /*if (rndmFlag == "on")
     {y0 = (Detector->GetScintillatorBoxSizeY())*12*(G4UniformRand()-0.5);
      z0 = (Detector->GetScintillatorBoxSizeZ())*1*(G4UniformRand()-0.5);
     } 

  particleGun->SetParticlePosition(G4ThreeVector(x0,y0,z0));

  particleGun->SetParticleTime(1.0*s);
  particleGun->GeneratePrimaryVertex(anEvent);*/
}


G4ParticleGun* PrimaryGeneratorAction::GetParticleGun()
{
        return particleGun;
}

void PrimaryGeneratorAction::PrintTargetSize()
{
  G4cout << "Shooting at area: " << Detector->GetBackScintillatorPlateSizeY()/cm << "cm x " << Detector->GetScintillatorBoxSizeZ()/cm << "cm  => " << (Detector->GetBackScintillatorPlateSizeY())*(Detector->GetScintillatorBoxSizeZ())/cm/cm << "cm^2"<< G4endl;
}