PndVolGenerator Class Reference

#include <PndVolGenerator.h>

Inheritance diagram for PndVolGenerator:

Public Member Functions
	PndVolGenerator ()
	PndVolGenerator (Int_t pdgid, Int_t mult=1)
virtual	~PndVolGenerator ()
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	SetVolXYZ (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)
void	SetVolTarg (Int_t quad=2, Int_t abs=5)
void	SetVolTgFc (Int_t quad=2, TF1 *func=0)
Bool_t	Init ()
virtual Bool_t	ReadEvent (FairPrimaryGenerator *primGen)

Private Member Functions
	ClassDef (PndVolGenerator, 3)

Private Attributes
Int_t	fPDGType
Int_t	fMult
Int_t	fQuad
Int_t	fAbs
TF1 *	fFunc
TGeoHMatrix *	fCurrentTransMat
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	fVolVtxIsSet
Bool_t	fDebug
Bool_t	fExt

Detailed Description

Definition at line 30 of file PndVolGenerator.h.

Constructor & Destructor Documentation

PndVolGenerator::PndVolGenerator

(

)

Default constructor.

Definition at line 25 of file PndVolGenerator.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),fPtRangeIsSet(0),fPRangeIsSet(0),fInversePIsSet(0),
  fAbs(0),fQuad(0),fFunc(0),fExt(0)
{
  //fGeoH = PndGeoHandling::Instance();
  // Default constructor
}

PndVolGenerator::PndVolGenerator	(	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 39 of file PndVolGenerator.cxx.

References SetPhiRange().

                                                        :
  fPDGType(pdgid),fMult(mult),fPDGMass(0),
  fEtaMin(0),fEtaMax(0),fYMin(0),fYMax(0),
  fPMin(0),fPMax(0),fX(0),fY(0),fZ(0),
  fX1(0),fY1(0),fX2(0),fY2(0),
  fEtaRangeIsSet(0),fYRangeIsSet(0),fThetaRangeIsSet(0),
  fCosThetaIsSet(0),fPtRangeIsSet(0),fPRangeIsSet(0),fInversePIsSet(0),
  fAbs(0),fQuad(0),fFunc(0),fExt(0)

{
  // Constructor. Set default kinematics limits
  SetPhiRange  ();
}

virtual PndVolGenerator::~PndVolGenerator ( )

inlinevirtual

Destructor

Definition at line 44 of file PndVolGenerator.h.

{};

Member Function Documentation

PndVolGenerator::ClassDef ( PndVolGenerator  , 3    )

private

Bool_t PndVolGenerator::Init

(

)

Initializer

Definition at line 53 of file PndVolGenerator.cxx.

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

{
  // Initialize generator

  if (fPhiMax-fPhiMin>360)
    Fatal("Init()","PndVolGenerator: phi range is too wide: %f<phi<%f",
          fPhiMin,fPhiMax);
  if (fPRangeIsSet && fPtRangeIsSet)
    Fatal("Init()","PndVolGenerator: Cannot set P and Pt ranges simultaneously");
  if (fPRangeIsSet && fYRangeIsSet)
    Fatal("Init()","PndVolGenerator: Cannot set P and Y ranges simultaneously");
  if ( (fThetaRangeIsSet && fYRangeIsSet) ||
       (fThetaRangeIsSet && fEtaRangeIsSet) ||
       (fYRangeIsSet     && fEtaRangeIsSet) )
    Fatal("Init()","PndVolGenerator: Cannot set Y, Theta or Eta ranges simultaneously");
  if (fPointVtxIsSet && fVolVtxIsSet)
    Fatal("Init()","PndVolGenerator: Cannot set point and box vertices simultaneously");
  if (fInversePIsSet && ( ( fPRangeIsSet && ( fPMin==0 || fPMax==0 ) ) || ( fPRangeIsSet && ( fPtMin==0 || fPtMax==0 ) ) ) )
    Fatal("Init()","PndVolGenerator: Cannot use P == 0 as limit for inverse momentum distribution");
  
  // Check for particle type
  //FairGeoLoader* geoLoad = new FairGeoLoader("TGeo","FairGeoLoader");
  //FairGeoInterface *geoFace = geoLoad->getGeoInterface();

 
 std::cout<<" gGeoManager 1 "<<gGeoManager<<" "<<std::endl;

  TDatabasePDG* pdgBase = TDatabasePDG::Instance();
  TParticlePDG *particle = pdgBase->GetParticle(fPDGType);
  if (! particle) Fatal("PndVolGenerator","PDG code %d not defined.",fPDGType);
  fPDGMass = particle->Mass();
}

Bool_t PndVolGenerator::ReadEvent ( FairPrimaryGenerator *  primGen )

virtual

Creates an event with given type and multiplicity.

Parameters

primGen

pointer to the FairPrimaryGenerator

Definition at line 87 of file PndVolGenerator.cxx.

References acos(), CAMath::Cos(), cos(), Double_t, eta, fAbs, fCosThetaIsSet, fCurrentTransMat, fDebug, fEtaMax, fEtaMin, fEtaRangeIsSet, fExt, fFunc, fInversePIsSet, fMult, fPDGMass, fPDGType, fPhiMax, fPhiMin, fPMax, fPMin, fPRangeIsSet, fPtMax, fPtMin, fPtRangeIsSet, fQuad, fThetaMax, fThetaMin, fThetaRangeIsSet, fVolVtxIsSet, fX, fX1, fX2, fY, fY1, fY2, fYMax, fYMin, fYRangeIsSet, fZ, gGeoManager, CAMath::Nint(), phi, printf(), pt(), pz, CAMath::Sin(), CAMath::Sqrt(), CAMath::Tan(), theta, top, 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;

  // TGeoVolume *top ;
  //if(gGeoManager)top = (TGeoVolume*)gGeoManager->GetTopVolume();

  TGeoVolume *top ;
  if(gGeoManager)top = (TGeoVolume*)gGeoManager->GetTopVolume();
  gGeoManager->GetCache()->BuildIdArray();
  std::cout<<" gGeoManager 2 "<<gGeoManager<<std::endl;

  //if(gGeoManager)std::cout<<" gGeoManager 2 "<<gGeoManager<<" "<<top<<std::endl;
  //std::cout<<" first volume "<<gGeoManager->GetCurrentVolume()->GetName()<<std::endl;

  if(top->GetName()!=gGeoManager->GetCurrentVolume()->GetName()){
    gGeoManager->CdTop();
  //  std::cout<<" first volume Top "<<gGeoManager->GetCurrentVolume()->GetName()<<" top "<<top->GetName()<<std::endl;
  }
  gGeoManager->CdDown(0); //complex
  //std::cout<<"  volume 1 "<<gGeoManager->GetPath()<<" "
        //   <<gGeoManager->GetCurrentVolume()->GetNdaughters()<<std::endl;

  Int_t Daught;
  Daught = gGeoManager->GetCurrentVolume()->GetNdaughters();
  
  gGeoManager->CdDown(fQuad); // stgQuad0 , 1 -> stgQuad2 , CdUp(), CdDown(2)

  //std::cout<<"  volume 2 "<<gGeoManager->GetPath()<<" "
        //   <<gGeoManager->GetCurrentVolume()->GetNdaughters()<<std::endl;

  if(fExt) {gGeoManager->CdDown(fAbs); //stglAb0 ,0 -> stglSi0_0 ,CdUp(), CdDown(1)
  }else {fAbs = TMath::Nint(fFunc->GetRandom());
          std::cout<<" random fabs "<<fAbs<<std::endl;
          gGeoManager->CdDown(fAbs);}
 
 fCurrentTransMat = gGeoManager->GetCurrentMatrix();
 Double_t posLab[3], posSens[3];
 
 posSens[0]=0.;  posSens[1]=0.;  posSens[2]=0.;
 fCurrentTransMat->LocalToMaster(posSens,posLab);
  
 fX = posLab[0];
 fY = posLab[1];
 fZ = posLab[2];
 
 std::cout<<" vertex "<<fX<<" "<<fY<<" "<<fZ<<std::endl;
 
  
  // 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 (fVolVtxIsSet) {
      fX = gRandom->Uniform(fX1,fX2);
      fY = gRandom->Uniform(fY1,fY2);
    }

    if (fDebug)
      printf("VolGen: 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;

}

void PndVolGenerator::SetCosTheta ( )

inline

Definition at line 69 of file PndVolGenerator.h.

References fCosThetaIsSet.

    {fCosThetaIsSet=kTRUE;};

void PndVolGenerator::SetDebug ( Bool_t  debug = 0 )

inline

Definition at line 81 of file PndVolGenerator.h.

References fDebug.

{fDebug = debug;}

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

inline

Definition at line 60 of file PndVolGenerator.h.

References fEtaMax, fEtaMin, and fEtaRangeIsSet.

    {fEtaMin=etamin; fEtaMax=etamax; fEtaRangeIsSet=kTRUE;};

void PndVolGenerator::SetInverseP ( )

inline

Definition at line 72 of file PndVolGenerator.h.

References fInversePIsSet.

    {fInversePIsSet=kTRUE;};

void PndVolGenerator::SetMultiplicity ( Int_t  mult )

inline

Definition at line 49 of file PndVolGenerator.h.

References fMult, and mult.

{fMult    = mult; };

void PndVolGenerator::SetPDGType ( Int_t  pdg )

inline

Modifiers

Definition at line 47 of file PndVolGenerator.h.

References fPDGType.

{fPDGType = pdg;  };

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

inline

Definition at line 57 of file PndVolGenerator.h.

References fPhiMax, and fPhiMin.

Referenced by PndVolGenerator().

    {fPhiMin=phimin; fPhiMax=phimax;};

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

inline

Definition at line 51 of file PndVolGenerator.h.

References fPMax, fPMin, and fPRangeIsSet.

Referenced by sim_Gamma().

    {fPMin=pmin; fPMax=pmax; fPRangeIsSet=kTRUE;}

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

inline

Definition at line 54 of file PndVolGenerator.h.

References fPtMax, fPtMin, and fPtRangeIsSet.

    {fPtMin=ptmin; fPtMax=ptmax; fPtRangeIsSet=kTRUE;};

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

inline

Definition at line 66 of file PndVolGenerator.h.

References fThetaMax, fThetaMin, and fThetaRangeIsSet.

    {fThetaMin=thetamin; fThetaMax=thetamax; fThetaRangeIsSet=kTRUE;};

void PndVolGenerator::SetVolTarg ( Int_t  quad = 2, Int_t  abs = 5  )

inline

Definition at line 83 of file PndVolGenerator.h.

References fAbs, fExt, and fQuad.

  {fQuad=quad; fAbs=abs;fExt=kTRUE;};        // Absorber Volume Setter

void PndVolGenerator::SetVolTgFc ( Int_t  quad = 2, TF1 *  func = 0  )

inline

Definition at line 86 of file PndVolGenerator.h.

References fExt, fFunc, fQuad, and func.

Referenced by sim_Gamma().

  {fQuad=quad; fFunc=func;fExt=kFALSE;};        // Absorber Volume Setter

void PndVolGenerator::SetVolXYZ ( Double32_t  x1 = 0, Double32_t  y1 = 0, Double32_t  x2 = 0, Double32_t  y2 = 0, Double32_t  z = 0  )

inline

Definition at line 78 of file PndVolGenerator.h.

References fVolVtxIsSet, fX1, fX2, fY1, fY2, fZ, and z.

                                                                                                      {
    fX1=x1; fY1=y1; fX2=x2; fY2=y2; fZ=z; fVolVtxIsSet=kTRUE;}

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

inline

Definition at line 75 of file PndVolGenerator.h.

References fPointVtxIsSet, fX, fY, fZ, x, y, and z.

                                                                 {
    fX=x; fY=y; fZ=z; fPointVtxIsSet=kTRUE;}

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

inline

Definition at line 63 of file PndVolGenerator.h.

References fYMax, fYMin, and fYRangeIsSet.

    {fYMin=ymin; fYMax=ymax; fYRangeIsSet=kTRUE;};

Member Data Documentation

Int_t PndVolGenerator::fAbs

private

Definition at line 101 of file PndVolGenerator.h.

Referenced by ReadEvent(), and SetVolTarg().

Bool_t PndVolGenerator::fCosThetaIsSet

private

Definition at line 119 of file PndVolGenerator.h.

Referenced by ReadEvent(), and SetCosTheta().

TGeoHMatrix* PndVolGenerator::fCurrentTransMat

private

Definition at line 103 of file PndVolGenerator.h.

Referenced by ReadEvent().

Bool_t PndVolGenerator::fDebug

private

Definition at line 125 of file PndVolGenerator.h.

Referenced by ReadEvent(), and SetDebug().

Double32_t PndVolGenerator::fEtaMax

private

Definition at line 108 of file PndVolGenerator.h.

Referenced by ReadEvent(), and SetEtaRange().

Double32_t PndVolGenerator::fEtaMin

private

Definition at line 108 of file PndVolGenerator.h.

Referenced by ReadEvent(), and SetEtaRange().

Bool_t PndVolGenerator::fEtaRangeIsSet

private

Definition at line 116 of file PndVolGenerator.h.

Referenced by Init(), ReadEvent(), and SetEtaRange().

Bool_t PndVolGenerator::fExt

private

Definition at line 126 of file PndVolGenerator.h.

Referenced by ReadEvent(), SetVolTarg(), and SetVolTgFc().

TF1* PndVolGenerator::fFunc

private

Definition at line 102 of file PndVolGenerator.h.

Referenced by ReadEvent(), and SetVolTgFc().

Bool_t PndVolGenerator::fInversePIsSet

private

Definition at line 120 of file PndVolGenerator.h.

Referenced by Init(), ReadEvent(), and SetInverseP().

Int_t PndVolGenerator::fMult

private

Definition at line 99 of file PndVolGenerator.h.

Referenced by ReadEvent(), and SetMultiplicity().

Double32_t PndVolGenerator::fPDGMass

private

Definition at line 105 of file PndVolGenerator.h.

Referenced by Init(), and ReadEvent().

Int_t PndVolGenerator::fPDGType

private

Definition at line 98 of file PndVolGenerator.h.

Referenced by Init(), ReadEvent(), and SetPDGType().

Double32_t PndVolGenerator::fPhiMax

private

Definition at line 107 of file PndVolGenerator.h.

Referenced by Init(), ReadEvent(), and SetPhiRange().

Double32_t PndVolGenerator::fPhiMin

private

Definition at line 107 of file PndVolGenerator.h.

Referenced by Init(), ReadEvent(), and SetPhiRange().

Double32_t PndVolGenerator::fPMax

private

Definition at line 110 of file PndVolGenerator.h.

Referenced by Init(), ReadEvent(), and SetPRange().

Double32_t PndVolGenerator::fPMin

private

Definition at line 110 of file PndVolGenerator.h.

Referenced by Init(), ReadEvent(), and SetPRange().

Bool_t PndVolGenerator::fPointVtxIsSet

private

Definition at line 123 of file PndVolGenerator.h.

Referenced by Init(), and SetXYZ().

Bool_t PndVolGenerator::fPRangeIsSet

private

Definition at line 122 of file PndVolGenerator.h.

Referenced by Init(), ReadEvent(), and SetPRange().

Double32_t PndVolGenerator::fPtMax

private

Definition at line 106 of file PndVolGenerator.h.

Referenced by Init(), ReadEvent(), and SetPtRange().

Double32_t PndVolGenerator::fPtMin

private

Definition at line 106 of file PndVolGenerator.h.

Referenced by Init(), ReadEvent(), and SetPtRange().

Bool_t PndVolGenerator::fPtRangeIsSet

private

Definition at line 121 of file PndVolGenerator.h.

Referenced by Init(), ReadEvent(), and SetPtRange().

Int_t PndVolGenerator::fQuad

private

Definition at line 100 of file PndVolGenerator.h.

Referenced by ReadEvent(), SetVolTarg(), and SetVolTgFc().

Double32_t PndVolGenerator::fThetaMax

private

Definition at line 111 of file PndVolGenerator.h.

Referenced by ReadEvent(), and SetThetaRange().

Double32_t PndVolGenerator::fThetaMin

private

Definition at line 111 of file PndVolGenerator.h.

Referenced by ReadEvent(), and SetThetaRange().

Bool_t PndVolGenerator::fThetaRangeIsSet

private

Definition at line 118 of file PndVolGenerator.h.

Referenced by Init(), ReadEvent(), and SetThetaRange().

Bool_t PndVolGenerator::fVolVtxIsSet

private

Definition at line 124 of file PndVolGenerator.h.

Referenced by Init(), ReadEvent(), and SetVolXYZ().

Double32_t PndVolGenerator::fX

private

Definition at line 112 of file PndVolGenerator.h.

Referenced by ReadEvent(), and SetXYZ().

Double32_t PndVolGenerator::fX1

private

Definition at line 113 of file PndVolGenerator.h.

Referenced by ReadEvent(), and SetVolXYZ().

Double32_t PndVolGenerator::fX2

private

Definition at line 113 of file PndVolGenerator.h.

Referenced by ReadEvent(), and SetVolXYZ().

Double32_t PndVolGenerator::fY

private

Definition at line 112 of file PndVolGenerator.h.

Referenced by ReadEvent(), and SetXYZ().

Double32_t PndVolGenerator::fY1

private

Definition at line 113 of file PndVolGenerator.h.

Referenced by ReadEvent(), and SetVolXYZ().

Double32_t PndVolGenerator::fY2

private

Definition at line 113 of file PndVolGenerator.h.

Referenced by ReadEvent(), and SetVolXYZ().

Double32_t PndVolGenerator::fYMax

private

Definition at line 109 of file PndVolGenerator.h.

Referenced by ReadEvent(), and SetYRange().

Double32_t PndVolGenerator::fYMin

private

Definition at line 109 of file PndVolGenerator.h.

Referenced by ReadEvent(), and SetYRange().

Bool_t PndVolGenerator::fYRangeIsSet

private

Definition at line 117 of file PndVolGenerator.h.

Referenced by Init(), ReadEvent(), and SetYRange().

Double32_t PndVolGenerator::fZ

private

Definition at line 112 of file PndVolGenerator.h.

Referenced by ReadEvent(), SetVolXYZ(), and SetXYZ().

---

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

• PndVolGenerator.h
• PndVolGenerator.cxx