Ex06PrimaryGenerator Class Reference

The primary generator. More...

#include <Ex06PrimaryGenerator.h>

Inheritance diagram for Ex06PrimaryGenerator:

Public Member Functions
	Ex06PrimaryGenerator (TVirtualMCStack *stack)
	Ex06PrimaryGenerator (const Ex06PrimaryGenerator &origin, TVirtualMCStack *stack)
	Ex06PrimaryGenerator ()
virtual	~Ex06PrimaryGenerator ()
void	GeneratePrimaries ()
void	SetParticle (Int_t pdg)
void	SetKinEnergy (Double_t kinEnergy)
void	SetDirection (Double_t dirX, Double_t dirY, Double_t dirZ)
void	SetOptPhotonPolar (Double_t angle)
void	SetNofPrimaries (Int_t nofPrimaries)

Private Member Functions
void	GeneratePrimary ()

Private Attributes
TVirtualMCStack *	fStack
	VMC stack.
Int_t	fPdg
	Particle PDG encoding.
Double_t	fKinEnergy
	Particle kinetic energy.
Double_t	fDirX
	Particle direction - x component.
Double_t	fDirY
	Particle direction - y component.
Double_t	fDirZ
	Particle direction - z component.
Double_t	fPolAngle
	Particle polarization angle.
Int_t	fNofPrimaries
	Number of primary particles.

Detailed Description

The primary generator.

Date

16/05/2005

Author

I. Hrivnacova; IPN, Orsay

Definition at line 33 of file Ex06PrimaryGenerator.h.

Constructor & Destructor Documentation

Ex06PrimaryGenerator() [1/3]

Ex06PrimaryGenerator::Ex06PrimaryGenerator

(

TVirtualMCStack *

stack

)

Standard constructor

Parameters

stack

The VMC stack

Definition at line 35 of file Ex06PrimaryGenerator.cxx.

  : TObject(),
    fStack(stack),
    fPdg(kPositron),
    fKinEnergy(500.e-06),
    fDirX(1.),
    fDirY(0.),
    fDirZ(0.),
    fPolAngle(0.),
    fNofPrimaries(1)
{
  /// Standard constructor
  /// \param stack  The VMC stack
}

Ex06PrimaryGenerator() [2/3]

Ex06PrimaryGenerator::Ex06PrimaryGenerator	(	const Ex06PrimaryGenerator &	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 51 of file Ex06PrimaryGenerator.cxx.

  : TObject(origin),
    fStack(stack),
    fPdg(origin.fPdg),
    fKinEnergy(origin.fKinEnergy),
    fDirX(origin.fDirX),
    fDirY(origin.fDirY),
    fDirZ(origin.fDirZ),
    fPolAngle(origin.fPolAngle),
    fNofPrimaries(origin.fNofPrimaries)
{
  /// Copy constructor (for clonig on worker thread in MT mode).
  /// \param origin  The source object (on master).
  /// \param stack   The VMC stack
}

Ex06PrimaryGenerator() [3/3]

Ex06PrimaryGenerator::Ex06PrimaryGenerator

(

)

Default constructor

Definition at line 69 of file Ex06PrimaryGenerator.cxx.

  : TObject(),
    fStack(0),
    fPdg(0),
    fKinEnergy(0.),
    fDirX(0.),
    fDirY(0.),
    fDirZ(0.),
    fPolAngle(0.),
    fNofPrimaries(0)
{
  /// Default constructor
}

~Ex06PrimaryGenerator()

Ex06PrimaryGenerator::~Ex06PrimaryGenerator ( )

virtual

Destructor

Definition at line 84 of file Ex06PrimaryGenerator.cxx.

{
  /// Destructor
}

Member Function Documentation

GeneratePrimaries()

void Ex06PrimaryGenerator::GeneratePrimaries

(

)

Fill the user stack (derived from TVirtualMCStack) with primary particle

Definition at line 158 of file Ex06PrimaryGenerator.cxx.

{
  /// Fill the user stack (derived from TVirtualMCStack) with primary particle
 
  for (Int_t i = 0; i < fNofPrimaries; i++) GeneratePrimary();
}

SetParticle()

void Ex06PrimaryGenerator::SetParticle ( Int_t pdg )

inline

Set particle type

Parameters

pdg	The new particle PDG encoding

Definition at line 73 of file Ex06PrimaryGenerator.h.

{ fPdg = pdg; }

SetKinEnergy()

void Ex06PrimaryGenerator::SetKinEnergy ( Double_t kinEnergy )

inline

Set kinetic energy

Parameters

kinEnergy

The new particle kinetic energy (in GeV)

Definition at line 77 of file Ex06PrimaryGenerator.h.

{
  fKinEnergy = kinEnergy;
}

SetDirection()

void Ex06PrimaryGenerator::SetDirection	(	Double_t	dirX,
		Double_t	dirY,
		Double_t	dirZ )

Set normalized direction

Parameters

dirX	The new direction - x component
dirY	The new direction - y component
dirZ	The new direction - z component

Definition at line 166 of file Ex06PrimaryGenerator.cxx.

{
  /// Set normalized direction
  /// \param dirX  The new direction - x component
  /// \param dirY  The new direction - y component
  /// \param dirZ  The new direction - z component
 
  Double_t norm = TMath::Sqrt(dirX * dirX + dirY * dirY + dirZ * dirZ);
 
  fDirX = dirX / norm;
  fDirY = dirY / norm;
  fDirZ = dirZ / norm;
}

SetOptPhotonPolar()

void Ex06PrimaryGenerator::SetOptPhotonPolar ( Double_t angle )

inline

Set photon polarization

Parameters

angle

The new polarization angle (in degrees)

Definition at line 84 of file Ex06PrimaryGenerator.h.

{
  fPolAngle = angle;
}

SetNofPrimaries()

void Ex06PrimaryGenerator::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 91 of file Ex06PrimaryGenerator.h.

{
  fNofPrimaries = nofPrimaries;
}

GeneratePrimary()

void Ex06PrimaryGenerator::GeneratePrimary ( )

private

Add one primary particle to the user stack (derived from TVirtualMCStack).

Definition at line 95 of file Ex06PrimaryGenerator.cxx.

{
  /// Add one primary particle to the user stack
  /// (derived from TVirtualMCStack).
 
  // Track ID (filled by stack)
  Int_t ntr;
 
  // Option: to be tracked
  Int_t toBeDone = 1;
 
  // Particle type
  Int_t pdg = fPdg;
  TParticlePDG* particlePDG = TDatabasePDG::Instance()->GetParticle(fPdg);
 
  // Position
  Double_t vx = 0.;
  Double_t vy = 0.;
  Double_t vz = 0.;
  Double_t tof = 0.;
 
  // Energy (in GeV)
  Double_t kinEnergy = fKinEnergy;
  Double_t mass = particlePDG->Mass();
  // Double_t mass = 0.51099906*1e-03;
  Double_t e = mass + kinEnergy;
 
  // Particle momentum
  Double_t p0, px, py, pz;
  p0 = sqrt(e * e - mass * mass);
  px = p0 * fDirX;
  py = p0 * fDirY;
  pz = p0 * fDirZ;
 
  // Polarization
  TVector3 polar;
  if (fPdg == 50000050) {
    TVector3 normal(1., 0., 0.);
    TVector3 kphoton = TVector3(fDirX, fDirY, fDirZ);
    TVector3 product = normal.Cross(kphoton);
    Double_t modul2 = product * product;
 
    TVector3 e_perpend(0., 0., 1.);
    if (modul2 > 0.) e_perpend = (1. / sqrt(modul2)) * product;
    TVector3 e_paralle = e_perpend.Cross(kphoton);
 
    polar = TMath::Cos(fPolAngle * TMath::DegToRad()) * e_paralle +
            TMath::Sin(fPolAngle * TMath::DegToRad()) * e_perpend;
  }
  // else
  //  Warning("GeneratePrimary",
  //          "The primary particle is not an opticalphoton");
 
  // Add particle to stack
  fStack->PushTrack(toBeDone, -1, pdg, px, py, pz, e, vx, vy, vz, tof,
    polar.X(), polar.Y(), polar.Z(), kPPrimary, ntr, 1., 0);
}

Member Data Documentation

fStack

TVirtualMCStack* Ex06PrimaryGenerator::fStack

private

VMC stack.

Definition at line 57 of file Ex06PrimaryGenerator.h.

fPdg

Int_t Ex06PrimaryGenerator::fPdg

private

Particle PDG encoding.

Definition at line 58 of file Ex06PrimaryGenerator.h.

fKinEnergy

Double_t Ex06PrimaryGenerator::fKinEnergy

private

Particle kinetic energy.

Definition at line 59 of file Ex06PrimaryGenerator.h.

fDirX

Double_t Ex06PrimaryGenerator::fDirX

private

Particle direction - x component.

Definition at line 60 of file Ex06PrimaryGenerator.h.

fDirY

Double_t Ex06PrimaryGenerator::fDirY

private

Particle direction - y component.

Definition at line 61 of file Ex06PrimaryGenerator.h.

fDirZ

Double_t Ex06PrimaryGenerator::fDirZ

private

Particle direction - z component.

Definition at line 62 of file Ex06PrimaryGenerator.h.

fPolAngle

Double_t Ex06PrimaryGenerator::fPolAngle

private

Particle polarization angle.

Definition at line 63 of file Ex06PrimaryGenerator.h.

fNofPrimaries

Int_t Ex06PrimaryGenerator::fNofPrimaries

private

Number of primary particles.

Definition at line 64 of file Ex06PrimaryGenerator.h.

---

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

• examples/E06/include/Ex06PrimaryGenerator.h
• examples/E06/src/Ex06PrimaryGenerator.cxx