A01PrimaryGenerator Class Reference

The primary generator. More...

#include <A01PrimaryGenerator.h>

Inheritance diagram for A01PrimaryGenerator:

Public Member Functions
	A01PrimaryGenerator (TVirtualMCStack *stack)
	A01PrimaryGenerator (const A01PrimaryGenerator &origin, TVirtualMCStack *stack)
	A01PrimaryGenerator ()
virtual	~A01PrimaryGenerator ()
virtual void	GeneratePrimaries ()
void	SetNofPrimaries (Int_t nofPrimaries)
void	SetMomentum (Double_t val)
void	SetSigmaMomentum (Double_t val)
void	SetSigmaAngle (Double_t val)
void	SetRandomize (Bool_t val)
Double_t	GetMomentum () const
Double_t	GetSigmaMomentum () const
Double_t	GetSigmaAngle () const
Bool_t	GetRandomize () const

Private Attributes
TVirtualMCStack *	fStack
	VMC stack.
Int_t	fNofPrimaries
	Number of primary particles.
Int_t	fDefaultParticle
	Default particle PDG.
Double_t	fMomentum
	Default particle momentum.
Double_t	fSigmaMomentum
	The sigma of particle momentum.
Double_t	fSigmaAngle
	The sigma of particle direction.
Bool_t	fRandomizePrimary
	Option to randomize primary particle type.

Detailed Description

The primary generator.

Date

12/05/2012

Author

I. Hrivnacova; IPN, Orsay

Definition at line 33 of file A01PrimaryGenerator.h.

Constructor & Destructor Documentation

A01PrimaryGenerator() [1/3]

A01PrimaryGenerator::A01PrimaryGenerator

(

TVirtualMCStack *

stack

)

Standard constructor

Parameters

stack

The VMC stack

Definition at line 35 of file A01PrimaryGenerator.cxx.

  : TObject(),
    fStack(stack),
    fNofPrimaries(1),
    fDefaultParticle(kMuonPlus),
    fMomentum(1.),           // 1 GeV
    fSigmaMomentum(50.e-03), // 50 MeV;
    fSigmaAngle(2.),         // 2 deg
    fRandomizePrimary(true)
{
  /// Standard constructor
  /// \param stack  The VMC stack
}

A01PrimaryGenerator() [2/3]

A01PrimaryGenerator::A01PrimaryGenerator	(	const A01PrimaryGenerator &	origin,
		TVirtualMCStack *	stack )

Copy constructor (for clonig on worker thread in MT mode).

Parameters

origin	The source object (on master).
stack	The VMC stack

Definition at line 50 of file A01PrimaryGenerator.cxx.

  : TObject(origin),
    fStack(stack),
    fNofPrimaries(origin.fNofPrimaries),
    fDefaultParticle(origin.fDefaultParticle),
    fMomentum(origin.fMomentum),
    fSigmaMomentum(origin.fSigmaMomentum),
    fSigmaAngle(origin.fSigmaAngle),
    fRandomizePrimary(origin.fRandomizePrimary)
{
  /// Copy constructor (for clonig on worker thread in MT mode).
  /// \param origin    The source object (on master).
  /// \param stack  The VMC stack
}

A01PrimaryGenerator() [3/3]

A01PrimaryGenerator::A01PrimaryGenerator

(

)

Default constructor

Definition at line 67 of file A01PrimaryGenerator.cxx.

  : TObject(),
    fStack(0),
    fNofPrimaries(0),
    fDefaultParticle(kMuonPlus),
    fMomentum(1.),           // 1 GeV
    fSigmaMomentum(50.e-03), // 50 MeV;
    fSigmaAngle(2.),         // 2 deg
    fRandomizePrimary(true)
{
  /// Default constructor
}

~A01PrimaryGenerator()

A01PrimaryGenerator::~A01PrimaryGenerator ( )

virtual

Destructor

Definition at line 81 of file A01PrimaryGenerator.cxx.

{
  /// Destructor
}

Member Function Documentation

GeneratePrimaries()

void A01PrimaryGenerator::GeneratePrimaries ( )

virtual

Fill the user stack (derived from TVirtualMCStack) with primary particles. All primaries in one event have the same properties.

Definition at line 91 of file A01PrimaryGenerator.cxx.

{
  /// Fill the user stack (derived from TVirtualMCStack) with primary particles.
  /// All primaries in one event have the same properties.
 
  // Track ID (filled by stack)
  Int_t ntr;
 
  // Option: to be tracked
  Int_t toBeDone = 1;
 
  Int_t pdg = fDefaultParticle;
  if (fRandomizePrimary) {
    // Int_t i = (Int_t)(2.*gMC->GetRandom()->Rndm());
    // Int_t i = (Int_t)(2.*0.3);
    static Int_t counter = 0;
    Int_t i = (counter++) % 5;
    switch (i) {
      case 0:
        pdg = kPositron;
        break;
      case 1:
        pdg = kMuonPlus;
        break;
      case 2:
        pdg = kPiPlus;
        break;
      case 3:
        pdg = kKPlus;
        break;
      case 4:
        pdg = kProton;
        break;
      default:
        pdg = kPositron;
        break;
    }
  }
 
  // Polarization
  Double_t polx = 0.;
  Double_t poly = 0.;
  Double_t polz = 0.;
 
  // Position
  Double_t vx = 0;
  Double_t vy = 0.;
  Double_t vz = -800.0; //- 8.*m;
  Double_t tof = 0.;
 
  // Particle momentum
  TDatabasePDG* databasePDG = TDatabasePDG::Instance();
  TParticlePDG* particlePDG = databasePDG->GetParticle(pdg);
  Double_t mass = particlePDG->Mass();
  // Double_t pp = fMomentum + (gMC->GetRandom()->Rndm()-0.5)*fSigmaMomentum;
  Double_t pp = fMomentum + (0.3 - 0.5) * fSigmaMomentum;
  Double_t e = TMath::Sqrt(pp * pp + mass * mass);
 
  // Double_t angle = (gMC->GetRandom()->Rndm()-0.5)*fSigmaAngle;
  Double_t angle = (0.3 - 0.5) * fSigmaAngle;
  Double_t px, py, pz;
  px = pp * TMath::Sin(angle * TMath::DegToRad());
  py = 0.;
  pz = pp * TMath::Cos(angle * TMath::DegToRad());
 
  for (Int_t i = 0; i < fNofPrimaries; i++) {
    // Add particle to stack
    fStack->PushTrack(toBeDone, -1, pdg, px, py, pz, e, vx, vy, vz, tof, polx,
      poly, polz, kPPrimary, ntr, 1., 0);
  }
}

SetNofPrimaries()

void A01PrimaryGenerator::SetNofPrimaries ( Int_t nofPrimaries )

inline

Set the number of particles to be generated

Parameters

nofPrimaries

The number of particles to be generated

Definition at line 78 of file A01PrimaryGenerator.h.

{
  fNofPrimaries = nofPrimaries;
}

SetMomentum()

void A01PrimaryGenerator::SetMomentum ( Double_t val )

inline

Definition at line 47 of file A01PrimaryGenerator.h.

{ fMomentum = val; }

SetSigmaMomentum()

void A01PrimaryGenerator::SetSigmaMomentum ( Double_t val )

inline

Definition at line 48 of file A01PrimaryGenerator.h.

{ fSigmaMomentum = val; }

SetSigmaAngle()

void A01PrimaryGenerator::SetSigmaAngle ( Double_t val )

inline

Definition at line 49 of file A01PrimaryGenerator.h.

{ fSigmaAngle = val; }

SetRandomize()

void A01PrimaryGenerator::SetRandomize ( Bool_t val )

inline

Definition at line 50 of file A01PrimaryGenerator.h.

{ fRandomizePrimary = val; }

GetMomentum()

Double_t A01PrimaryGenerator::GetMomentum ( ) const

inline

Definition at line 53 of file A01PrimaryGenerator.h.

{ return fMomentum; }

GetSigmaMomentum()

Double_t A01PrimaryGenerator::GetSigmaMomentum ( ) const

inline

Definition at line 54 of file A01PrimaryGenerator.h.

{ return fSigmaMomentum; }

GetSigmaAngle()

Double_t A01PrimaryGenerator::GetSigmaAngle ( ) const

inline

Definition at line 55 of file A01PrimaryGenerator.h.

{ return fSigmaAngle; }

GetRandomize()

Bool_t A01PrimaryGenerator::GetRandomize ( ) const

inline

Definition at line 56 of file A01PrimaryGenerator.h.

{ return fRandomizePrimary; }

Member Data Documentation

fStack

TVirtualMCStack* A01PrimaryGenerator::fStack

private

VMC stack.

Definition at line 63 of file A01PrimaryGenerator.h.

fNofPrimaries

Int_t A01PrimaryGenerator::fNofPrimaries

private

Number of primary particles.

Definition at line 64 of file A01PrimaryGenerator.h.

fDefaultParticle

Int_t A01PrimaryGenerator::fDefaultParticle

private

Default particle PDG.

Definition at line 65 of file A01PrimaryGenerator.h.

fMomentum

Double_t A01PrimaryGenerator::fMomentum

private

Default particle momentum.

Definition at line 66 of file A01PrimaryGenerator.h.

fSigmaMomentum

Double_t A01PrimaryGenerator::fSigmaMomentum

private

The sigma of particle momentum.

Definition at line 67 of file A01PrimaryGenerator.h.

fSigmaAngle

Double_t A01PrimaryGenerator::fSigmaAngle

private

The sigma of particle direction.

Definition at line 68 of file A01PrimaryGenerator.h.

fRandomizePrimary

Bool_t A01PrimaryGenerator::fRandomizePrimary

private

Option to randomize primary particle type.

Definition at line 69 of file A01PrimaryGenerator.h.

---

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

• examples/A01/include/A01PrimaryGenerator.h
• examples/A01/src/A01PrimaryGenerator.cxx