#include <PndBoxGenerator.h>

Inheritance diagram for PndBoxGenerator:

Public Member Functions
	PndBoxGenerator ()

	PndBoxGenerator (Int_t pdgid, Int_t mult=1)

virtual	~PndBoxGenerator ()

void	SetPDGType (Int_t pdg)

void	SetMultiplicity (Int_t mult)

void	SetPRange (Double32_t pmin=0, Double32_t pmax=10)

void	SetPtRange (Double32_t ptmin=0, Double32_t ptmax=10)

void	SetPhiRange (Double32_t phimin=0, Double32_t phimax=360)

void	SetEtaRange (Double32_t etamin=-5, Double32_t etamax=7)

void	SetYRange (Double32_t ymin=-5, Double32_t ymax=7)

void	SetThetaRange (Double32_t thetamin=0, Double32_t thetamax=90)

void	SetCosTheta ()

void	SetInverseP ()

void	SetXYZ (Double32_t x=0, Double32_t y=0, Double32_t z=0)

void	SetBoxXYZ (Double32_t x1=0, Double32_t y1=0, Double32_t x2=0, Double32_t y2=0, Double32_t z=0)

void	SetDebug (Bool_t debug=0)

Bool_t	Init ()

virtual Bool_t	ReadEvent (FairPrimaryGenerator *primGen)

Private Member Functions
	ClassDef (PndBoxGenerator, 1)

Private Attributes
Int_t	fPDGType

Int_t	fMult

Double32_t	fPDGMass

Double32_t	fPtMin

Double32_t	fPtMax

Double32_t	fPhiMin

Double32_t	fPhiMax

Double32_t	fEtaMin

Double32_t	fEtaMax

Double32_t	fYMin

Double32_t	fYMax

Double32_t	fPMin

Double32_t	fPMax

Double32_t	fThetaMin

Double32_t	fThetaMax

Double32_t	fX

Double32_t	fY

Double32_t	fZ

Double32_t	fX1

Double32_t	fY1

Double32_t	fX2

Double32_t	fY2

Bool_t	fEtaRangeIsSet

Bool_t	fYRangeIsSet

Bool_t	fThetaRangeIsSet

Bool_t	fCosThetaIsSet

Bool_t	fInversePIsSet

Bool_t	fPtRangeIsSet

Bool_t	fPRangeIsSet

Bool_t	fPointVtxIsSet

Bool_t	fBoxVtxIsSet

Bool_t	fDebug

Detailed Description

Definition at line 25 of file PndBoxGenerator.h.

Constructor & Destructor Documentation

PndBoxGenerator::PndBoxGenerator ( )

Default constructor.

Definition at line 15 of file PndBoxGenerator.cxx.

                                  :
   fPDGType(0),fMult(0),fPDGMass(0),fPtMin(0),fPtMax(0),
   fPhiMin(0),fPhiMax(0),fEtaMin(0),fEtaMax(0),fYMin(0),fYMax(0),
   fPMin(0),fPMax(0),fThetaMin(0),fThetaMax(0),fX(0),fY(0),fZ(0),
   fX1(0),fY1(0),fX2(0),fY2(0),
   fEtaRangeIsSet(0),fYRangeIsSet(0),fThetaRangeIsSet(0),
   fCosThetaIsSet(0),fInversePIsSet(0),fPtRangeIsSet(0),fPRangeIsSet(0),
   fPointVtxIsSet(0), fBoxVtxIsSet(0)
 {
   // Default constructor
 }

PndBoxGenerator::PndBoxGenerator	(	Int_t	pdgid,
		Int_t	mult = `1`
	)

Constructor with PDG-ID, multiplicity

Parameters

pdgid	Particle type (PDG encoding)
mult	Multiplicity (default is 1)

Definition at line 28 of file PndBoxGenerator.cxx.

References SetPhiRange().

                                                         :
   fPDGType(pdgid),fMult(mult),fPDGMass(0),fPtMin(0),fPtMax(0),
   fPhiMin(0),fPhiMax(0), fEtaMin(0),fEtaMax(0),fYMin(0),fYMax(0),
   fPMin(0),fPMax(0),fThetaMin(0),fThetaMax(0),fX(0),fY(0),fZ(0),
   fX1(0),fY1(0),fX2(0),fY2(0),
   fEtaRangeIsSet(0),fYRangeIsSet(0),fThetaRangeIsSet(0),
   fCosThetaIsSet(0),fInversePIsSet(0),fPtRangeIsSet(0),fPRangeIsSet(0),
   fPointVtxIsSet(0), fBoxVtxIsSet(0)
 
 {
   // Constructor. Set default kinematics limits
   SetPhiRange  ();
 }

virtual PndBoxGenerator::~PndBoxGenerator ( )

inlinevirtual

Destructor

Definition at line 39 of file PndBoxGenerator.h.

39 {};

Member Function Documentation

PndBoxGenerator::ClassDef	(	PndBoxGenerator	,
		1
	)

private

Bool_t PndBoxGenerator::Init ( )

Initializer

Definition at line 42 of file PndBoxGenerator.cxx.

References fBoxVtxIsSet, fEtaRangeIsSet, fInversePIsSet, fPDGMass, fPDGType, fPhiMax, fPhiMin, fPMax, fPMin, fPointVtxIsSet, fPRangeIsSet, fPtMax, fPtMin, fPtRangeIsSet, fThetaRangeIsSet, fYRangeIsSet, and particle.

 {
   // Initialize generator
 
   if (fPhiMax-fPhiMin>360)
     Fatal("Init()","PndBoxGenerator: phi range is too wide: %f<phi<%f",
           fPhiMin,fPhiMax);
   if (fPRangeIsSet && fPtRangeIsSet)
     Fatal("Init()","PndBoxGenerator: Cannot set P and Pt ranges simultaneously");
   if (fPRangeIsSet && fYRangeIsSet)
     Fatal("Init()","PndBoxGenerator: Cannot set P and Y ranges simultaneously");
   if ( (fThetaRangeIsSet && fYRangeIsSet) ||
        (fThetaRangeIsSet && fEtaRangeIsSet) ||
        (fYRangeIsSet     && fEtaRangeIsSet) )
     Fatal("Init()","PndBoxGenerator: Cannot set Y, Theta or Eta ranges simultaneously");
   if (fPointVtxIsSet && fBoxVtxIsSet)
     Fatal("Init()","PndBoxGenerator: Cannot set point and box vertices simultaneously");
   if (fInversePIsSet && ( ( fPRangeIsSet && ( fPMin==0 || fPMax==0 ) ) || ( fPRangeIsSet && ( fPtMin==0 || fPtMax==0 ) ) ) )
     Fatal("Init()","PndBoxGenerator: Cannot use P == 0 as limit for inverse momentum distribution");
   
   // Check for particle type
   TDatabasePDG* pdgBase = TDatabasePDG::Instance();
   TParticlePDG *particle = pdgBase->GetParticle(fPDGType);
   if (! particle) Fatal("PndBoxGenerator","PDG code %d not defined.",fPDGType);
   fPDGMass = particle->Mass();
   return kTRUE;
 }

Bool_t PndBoxGenerator::ReadEvent ( FairPrimaryGenerator * primGen )

virtual

Creates an event with given type and multiplicity.

Parameters

primGen pointer to the FairPrimaryGenerator

Definition at line 71 of file PndBoxGenerator.cxx.

References acos(), CAMath::Cos(), cos(), eta, fBoxVtxIsSet, fCosThetaIsSet, fDebug, fEtaMax, fEtaMin, fEtaRangeIsSet, fInversePIsSet, fMult, fPDGMass, fPDGType, fPhiMax, fPhiMin, fPMax, fPMin, fPRangeIsSet, fPtMax, fPtMin, fPtRangeIsSet, fThetaMax, fThetaMin, fThetaRangeIsSet, fX, fX1, fX2, fY, fY1, fY2, fYMax, fYMin, fYRangeIsSet, fZ, phi, printf(), pt(), pz, CAMath::Sin(), CAMath::Sqrt(), CAMath::Tan(), theta, and y.

 {
   // Generate one event: produce primary particles emitted from one vertex.
   // Primary particles are distributed uniformly along
   // those kinematics variables which were limitted by setters.
   // if SetCosTheta() function is used, the distribution will be uniform in cos(theta)
   
   Double32_t pabs, phi, pt, theta=0, eta, y, mt, px, py, pz, pinv=0;
 
   // Generate particles
   for (Int_t k = 0; k < fMult; k++) {
     phi = gRandom->Uniform(fPhiMin,fPhiMax) * TMath::DegToRad();
 
     if      (fPRangeIsSet ) {
       if (fInversePIsSet)
         {
           pinv = gRandom->Uniform(1./fPMax, 1./fPMin);
           pabs = 1./pinv;
         }
       else
         {
           pabs = gRandom->Uniform(fPMin,fPMax);
         }
     }
     else if (fPtRangeIsSet) {
       if (fInversePIsSet)
         {
           pinv = gRandom->Uniform(1./fPtMax, 1./fPtMin);
           pt = 1./pinv;
         }
       else
         {
           pt   = gRandom->Uniform(fPtMin,fPtMax);
         }
     }
     
     if      (fThetaRangeIsSet) {
       if (fCosThetaIsSet)
         theta = acos(gRandom->Uniform(cos(fThetaMin* TMath::DegToRad()),cos(fThetaMax* TMath::DegToRad())));
       else
         theta = gRandom->Uniform(fThetaMin,fThetaMax) * TMath::DegToRad();
     }
     else if (fEtaRangeIsSet) {
       eta   = gRandom->Uniform(fEtaMin,fEtaMax);
       theta = 2*TMath::ATan(TMath::Exp(-eta));
     }
     else if (fYRangeIsSet) {
       y     = gRandom->Uniform(fYMin,fYMax);
       mt = TMath::Sqrt(fPDGMass*fPDGMass + pt*pt);
       pz = mt * TMath::SinH(y);
     }
     
     if (fThetaRangeIsSet || fEtaRangeIsSet) {
       if      (fPRangeIsSet ) {
         pz = pabs*TMath::Cos(theta);
         pt = pabs*TMath::Sin(theta);
       }
       else if (fPtRangeIsSet) 
         pz = pt/TMath::Tan(theta);
     }
 
     px = pt*TMath::Cos(phi);
     py = pt*TMath::Sin(phi);
 
     if (fBoxVtxIsSet) {
       fX = gRandom->Uniform(fX1,fX2);
       fY = gRandom->Uniform(fY1,fY2);
     }
 
     if (fDebug)
       printf("BoxGen: kf=%d, p=(%.2f, %.2f, %.2f) GeV, x=(%.1f, %.1f, %.1f) cm\n",
              fPDGType, px, py, pz, fX, fY, fZ);
 
     primGen->AddTrack(fPDGType, px, py, pz, fX, fY, fZ);
   }
   return kTRUE;
 
 }