B5PrimaryGeneratorAction Class Reference

#include <B5PrimaryGeneratorAction.hh>

Inheritance diagram for B5PrimaryGeneratorAction:

Collaboration diagram for B5PrimaryGeneratorAction:

Public Member Functions
	B5PrimaryGeneratorAction ()
virtual	~B5PrimaryGeneratorAction ()
virtual void	GeneratePrimaries (G4Event *)
void	SetMomentum (G4double val)
G4double	GetMomentum () const
void	SetSigmaMomentum (G4double val)
G4double	GetSigmaMomentum () const
void	SetSigmaAngle (G4double val)
G4double	GetSigmaAngle () const
void	SetRandomize (G4bool val)
G4bool	GetRandomize () const
Public Member Functions inherited from G4VUserPrimaryGeneratorAction
	G4VUserPrimaryGeneratorAction ()
virtual	~G4VUserPrimaryGeneratorAction ()

Private Member Functions
void	DefineCommands ()

Private Attributes
G4ParticleGun *	fParticleGun
G4GenericMessenger *	fMessenger
G4ParticleDefinition *	fPositron
G4ParticleDefinition *	fMuon
G4ParticleDefinition *	fPion
G4ParticleDefinition *	fKaon
G4ParticleDefinition *	fProton
G4double	fMomentum
G4double	fSigmaMomentum
G4double	fSigmaAngle
G4bool	fRandomizePrimary

Detailed Description

Primary generator

A single particle is generated. User can select

• the initial momentum and angle
• the momentum and angle spreads
• random selection of a particle type from proton, kaon+, pi+, muon+, e+

Definition at line 51 of file B5PrimaryGeneratorAction.hh.

Constructor & Destructor Documentation

B5PrimaryGeneratorAction()

B5PrimaryGeneratorAction::B5PrimaryGeneratorAction

(

)

Definition at line 43 of file B5PrimaryGeneratorAction.cc.

: G4VUserPrimaryGeneratorAction(),     
  fParticleGun(0), fMessenger(0), 
  fPositron(0), fMuon(0), fPion(0), fKaon(0), fProton(0),
  fMomentum(1000.*MeV),
  fSigmaMomentum(50.*MeV),
  fSigmaAngle(2.*deg),
  fRandomizePrimary(true)
{
    G4int n_particle = 1;
    fParticleGun  = new G4ParticleGun(n_particle);
    
    G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();
    G4String particleName;
    fPositron = particleTable->FindParticle(particleName="e+");
    fMuon = particleTable->FindParticle(particleName="mu+");
    fPion = particleTable->FindParticle(particleName="pi+");
    fKaon = particleTable->FindParticle(particleName="kaon+");
    fProton = particleTable->FindParticle(particleName="proton");
    
    // default particle kinematics
    fParticleGun->SetParticlePosition(G4ThreeVector(0.,0.,-8.*m));
    fParticleGun->SetParticleDefinition(fPositron);
    
    // define commands for this class
    DefineCommands();
}

Here is the call graph for this function:

~B5PrimaryGeneratorAction()

B5PrimaryGeneratorAction::~B5PrimaryGeneratorAction ( )

virtual

Definition at line 73 of file B5PrimaryGeneratorAction.cc.

{
    delete fParticleGun;
    delete fMessenger;
}

Member Function Documentation

DefineCommands()

void B5PrimaryGeneratorAction::DefineCommands ( )

private

Definition at line 127 of file B5PrimaryGeneratorAction.cc.

{
    // Define /B5/generator command directory using generic messenger class
    fMessenger 
      = new G4GenericMessenger(this, 
                               "/B5/generator/", 
                               "Primary generator control");
              
    // momentum command
    G4GenericMessenger::Command& momentumCmd
      = fMessenger->DeclarePropertyWithUnit("momentum", "GeV", fMomentum, 
                                    "Mean momentum of primaries.");
    momentumCmd.SetParameterName("p", true);
    momentumCmd.SetRange("p>=0.");                                
    momentumCmd.SetDefaultValue("1.");
    // ok
    //momentumCmd.SetParameterName("p", true);
    //momentumCmd.SetRange("p>=0.");                                
    
    // sigmaMomentum command
    G4GenericMessenger::Command& sigmaMomentumCmd
      = fMessenger->DeclarePropertyWithUnit("sigmaMomentum",
          "MeV", fSigmaMomentum, "Sigma momentum of primaries.");
    sigmaMomentumCmd.SetParameterName("sp", true);
    sigmaMomentumCmd.SetRange("sp>=0.");                                
    sigmaMomentumCmd.SetDefaultValue("50.");

    // sigmaAngle command
    G4GenericMessenger::Command& sigmaAngleCmd
      = fMessenger->DeclarePropertyWithUnit("sigmaAngle", "deg", fSigmaAngle, 
                                    "Sigma angle divergence of primaries.");
    sigmaAngleCmd.SetParameterName("t", true);
    sigmaAngleCmd.SetRange("t>=0.");                                
    sigmaAngleCmd.SetDefaultValue("2.");

    // randomizePrimary command
    G4GenericMessenger::Command& randomCmd
      = fMessenger->DeclareProperty("randomizePrimary", fRandomizePrimary);
    G4String guidance
       = "Boolean flag for randomizing primary particle types.\n";   
    guidance
       += "In case this flag is false, you can select the primary particle\n";
    guidance += "  with /gun/particle command.";                               
    randomCmd.SetGuidance(guidance);
    randomCmd.SetParameterName("flg", true);
    randomCmd.SetDefaultValue("true");
}

Here is the call graph for this function:

Here is the caller graph for this function:

GeneratePrimaries()

void B5PrimaryGeneratorAction::GeneratePrimaries ( G4Event *  event )

virtual

Implements G4VUserPrimaryGeneratorAction.

Definition at line 81 of file B5PrimaryGeneratorAction.cc.

{
    G4ParticleDefinition* particle;
    
    if (fRandomizePrimary)
    {
        G4int i = (int)(5.*G4UniformRand());
        switch(i)
        {
            case 0:
                particle = fPositron;
                break;
            case 1:
                particle = fMuon;
                break;
            case 2:
                particle = fPion;
                break;
            case 3:
                particle = fKaon;
                break;
            default:
                particle = fProton;
                break;
        }
        fParticleGun->SetParticleDefinition(particle);
    }
    else
    {
        particle = fParticleGun->GetParticleDefinition();
    }
    
    G4double pp = fMomentum + (G4UniformRand()-0.5)*fSigmaMomentum;
    G4double mass = particle->GetPDGMass();
    G4double Ekin = std::sqrt(pp*pp+mass*mass)-mass;
    fParticleGun->SetParticleEnergy(Ekin);
    
    G4double angle = (G4UniformRand()-0.5)*fSigmaAngle;
    fParticleGun->SetParticleMomentumDirection(G4ThreeVector(std::sin(angle),0.,
                                                             std::cos(angle)));
    
    fParticleGun->GeneratePrimaryVertex(event);
}

Here is the call graph for this function:

GetMomentum()

G4double B5PrimaryGeneratorAction::GetMomentum ( ) const

inline

Definition at line 60 of file B5PrimaryGeneratorAction.hh.

{ return fMomentum; }

GetRandomize()

G4bool B5PrimaryGeneratorAction::GetRandomize ( ) const

inline

Definition at line 69 of file B5PrimaryGeneratorAction.hh.

{ return fRandomizePrimary; }

Here is the call graph for this function:

GetSigmaAngle()

G4double B5PrimaryGeneratorAction::GetSigmaAngle ( ) const

inline

Definition at line 66 of file B5PrimaryGeneratorAction.hh.

{ return fSigmaAngle; }

GetSigmaMomentum()

G4double B5PrimaryGeneratorAction::GetSigmaMomentum ( ) const

inline

Definition at line 63 of file B5PrimaryGeneratorAction.hh.

{ return fSigmaMomentum; }

SetMomentum()

void B5PrimaryGeneratorAction::SetMomentum ( G4double  val )

inline

Definition at line 59 of file B5PrimaryGeneratorAction.hh.

{ fMomentum = val; }

SetRandomize()

void B5PrimaryGeneratorAction::SetRandomize ( G4bool  val )

inline

Definition at line 68 of file B5PrimaryGeneratorAction.hh.

{ fRandomizePrimary = val; }

SetSigmaAngle()

void B5PrimaryGeneratorAction::SetSigmaAngle ( G4double  val )

inline

Definition at line 65 of file B5PrimaryGeneratorAction.hh.

{ fSigmaAngle = val; }

SetSigmaMomentum()

void B5PrimaryGeneratorAction::SetSigmaMomentum ( G4double  val )

inline

Definition at line 62 of file B5PrimaryGeneratorAction.hh.

{ fSigmaMomentum = val; }

Member Data Documentation

fKaon

G4ParticleDefinition* B5PrimaryGeneratorAction::fKaon

private

Definition at line 79 of file B5PrimaryGeneratorAction.hh.

fMessenger

G4GenericMessenger* B5PrimaryGeneratorAction::fMessenger

private

Definition at line 75 of file B5PrimaryGeneratorAction.hh.

fMomentum

G4double B5PrimaryGeneratorAction::fMomentum

private

Definition at line 81 of file B5PrimaryGeneratorAction.hh.

fMuon

G4ParticleDefinition* B5PrimaryGeneratorAction::fMuon

private

Definition at line 77 of file B5PrimaryGeneratorAction.hh.

fParticleGun

G4ParticleGun* B5PrimaryGeneratorAction::fParticleGun

private

Definition at line 74 of file B5PrimaryGeneratorAction.hh.

fPion

G4ParticleDefinition* B5PrimaryGeneratorAction::fPion

private

Definition at line 78 of file B5PrimaryGeneratorAction.hh.

fPositron

G4ParticleDefinition* B5PrimaryGeneratorAction::fPositron

private

Definition at line 76 of file B5PrimaryGeneratorAction.hh.

fProton

G4ParticleDefinition* B5PrimaryGeneratorAction::fProton

private

Definition at line 80 of file B5PrimaryGeneratorAction.hh.

fRandomizePrimary

G4bool B5PrimaryGeneratorAction::fRandomizePrimary

private

Definition at line 84 of file B5PrimaryGeneratorAction.hh.

fSigmaAngle

G4double B5PrimaryGeneratorAction::fSigmaAngle

private

Definition at line 83 of file B5PrimaryGeneratorAction.hh.

fSigmaMomentum

G4double B5PrimaryGeneratorAction::fSigmaMomentum

private

Definition at line 82 of file B5PrimaryGeneratorAction.hh.

---

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

• B5PrimaryGeneratorAction.hh
• B5PrimaryGeneratorAction.cc