#include <PndHypBupGenerator.h>

Inheritance diagram for PndHypBupGenerator:

Public Member Functions
	PndHypBupGenerator ()

	PndHypBupGenerator (const char *fileName)

virtual	~PndHypBupGenerator ()

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

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

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

void	GammaEmissPar (Bool_t On=kFALSE)

virtual Bool_t	ReadEvent (FairPrimaryGenerator *primGen)

	PndHypBupGenerator ()

	PndHypBupGenerator (const char *fileName)

virtual	~PndHypBupGenerator ()

virtual Bool_t	ReadEvent (FairPrimaryGenerator *primGen)

Private Member Functions
void	CloseInput ()

	ClassDef (PndHypBupGenerator, 2)

void	CloseInput ()
	Particle array from PLUTO. More...

	ClassDef (PndHypBupGenerator, 1)

Private Attributes
Int_t	iEvent
	PDG database. More...

const Char_t *	fFileName
	Event number. More...

TFile *	fInputFile
	Input file name. More...

TTree *	fInputTree
	Pointer to input file. More...

TClonesArray *	fParticles
	Pointer to input tree. More...

Double32_t	fPMin
	Particle array. More...

Double32_t	fPMax

Double32_t	fPhiMin
	Transverse momentum range [GeV]. More...

Double32_t	fPhiMax

Double32_t	fThetaMin
	Azimuth angle range [degree]. More...

Double32_t	fThetaMax

Bool_t	fThetaRangeIsSet
	Polar angle range in lab system [degree]. More...

Bool_t	fPRangeIsSet

Bool_t	fGamOn

Detailed Description

Definition at line 36 of file detectors/hyp/hypTools/PndHypBupGenerator.h.

Constructor & Destructor Documentation

PndHypBupGenerator::PndHypBupGenerator ( )

Default constructor without arguments should not be used.

Definition at line 25 of file detectors/hyp/hypTools/PndHypBupGenerator.cxx.

References fInputFile, fInputTree, and iEvent.

                                       :fPMin(0),fPMax(0),
 fThetaMin(0),fThetaMax(0),fPhiMin(0),fPhiMax(0),
 fThetaRangeIsSet(0),fPRangeIsSet(0),fGamOn(0) {
  iEvent     = 0;
  fInputFile = NULL;
  fInputTree = NULL;
 }

PndHypBupGenerator::PndHypBupGenerator ( const char * fileName )

Standard constructor.

Parameters

fileName The input file name

Definition at line 37 of file detectors/hyp/hypTools/PndHypBupGenerator.cxx.

References fFileName, fInputFile, fInputTree, fParticles, and iEvent.

                                                           :fPMin(0),fPMax(0),
 fThetaMin(0),fThetaMax(0),fPhiMin(0),fPhiMax(0),
 fThetaRangeIsSet(0),fPRangeIsSet(0),fGamOn(0){
 
  
   fFileName  = fileName;
    iEvent     = 0;
    fInputFile = new TFile(fFileName);
   fInputTree = (TTree*) fInputFile->Get("data");
   fParticles = new TClonesArray("THParticle",100);
   fInputTree->SetBranchAddress("Particles", &fParticles);
 
  
    
 }

PndHypBupGenerator::~PndHypBupGenerator ( )

virtual

Destructor.

Definition at line 57 of file detectors/hyp/hypTools/PndHypBupGenerator.cxx.

References CloseInput().

                                         {
   CloseInput();
 }

PndHypBupGenerator::PndHypBupGenerator ( )

Default constructor without arguments should not be used.

PndHypBupGenerator::PndHypBupGenerator ( const char * fileName )

Standard constructor.

Parameters

fileName The input file name

virtual PndHypBupGenerator::~PndHypBupGenerator ( )

virtual

Destructor.

Member Function Documentation

PndHypBupGenerator::ClassDef	(	PndHypBupGenerator	,
		1
	)

private

Private method RegisterIons. Goes through the input file and registers any ion needed. STL map from ion name to CbmIon

PndHypBupGenerator::ClassDef	(	PndHypBupGenerator	,
		2
	)

private

Private method RegisterIons. Goes through the input file and registers any ion needed. STL map from ion name to CbmIon

void PndHypBupGenerator::CloseInput ( )

private

Particle array from PLUTO.

Private method CloseInput. Just for convenience. Closes the input file properly. Called from destructor and from ReadEvent.

void PndHypBupGenerator::CloseInput ( )

private

Private method CloseInput. Just for convenience. Closes the input file properly. Called from destructor and from ReadEvent.

Definition at line 168 of file detectors/hyp/hypTools/PndHypBupGenerator.cxx.

References fFileName, and fInputFile.

Referenced by ReadEvent(), and ~PndHypBupGenerator().

                                     {
  
     if ( fInputFile) {
        cout << "-I- CbmHypBupGenerator: Closing input file " 
             << fFileName << endl;
        fInputFile->Close();
     
        delete fInputFile;
     }
     
     fInputFile = NULL;
   }

void PndHypBupGenerator::GammaEmissPar ( Bool_t On = kFALSE )

inline

Definition at line 65 of file detectors/hyp/hypTools/PndHypBupGenerator.h.

References fGamOn.

Referenced by run(), run_test(), and sim_Gamma().

65 {fGamOn=On;}

PndHypBupGenerator::fGamOn

Bool_t fGamOn

Definition: detectors/hyp/hypTools/PndHypBupGenerator.h:89

virtual Bool_t PndHypBupGenerator::ReadEvent ( FairPrimaryGenerator * primGen )

virtual

Reads on event from the input file and pushes the tracks onto the stack. Abstract method in base class.

Parameters

primGen pointer to the CbmPrimaryGenerator

Bool_t PndHypBupGenerator::ReadEvent ( FairPrimaryGenerator * primGen )

virtual

Reads on event from the input file and pushes the tracks onto the stack. Abstract method in base class.

Parameters

primGen pointer to the CbmPrimaryGenerator

Definition at line 65 of file detectors/hyp/hypTools/PndHypBupGenerator.cxx.

References acos(), CloseInput(), CAMath::Cos(), cos(), Double_t, fGamOn, fInputFile, fInputTree, fParticles, fPMax, fPMin, fThetaMax, fThetaMin, g, THParticle::GetPdgCode(), iEvent, p, phi, pz, CAMath::Sin(), theta, THParticle::Vx(), THParticle::Vy(), and THParticle::Vz().

                                                                   {
 
   // Check for input file
   if ( ! fInputFile ) {
     cout << "-E- PndHypBupGenerator: Input file not open!" << endl;
     return kFALSE;
   }
  // Check for number of events in input file
   if ( iEvent > fInputTree->GetEntries() ) {
     cout << "-E PndDpmGenerator: No more events in input file!" << endl;
     CloseInput();
     return kFALSE;
   }
 
   // Defining kinematics variables, angle is expressed in degrees.
    Double_t phi, p=0.0,theta=0.0;
 
   TFile  *g=gFile;
   fInputFile->cd();
   fInputTree->GetEntry(iEvent++);
   g->cd();
 
    // Get number of particle in TClonesrray
   Int_t nParts = fParticles->GetEntriesFast();
 
   // Loop over particles in TClonesArray
   
   for (Int_t iPart=0; iPart < nParts; iPart++) {
     THParticle* part = (THParticle*) fParticles->At(iPart);
     Int_t pdgType = part->GetPdgCode();
     
     // Check if particle type is known to database
     if ( ! pdgType ) {
       cout << "-W PndDpmGenerator: Unknown type " << part->GetPdgCode()
            << ", skipping particle." << endl;
       continue;
     }
 
     Double_t px = 0.;
     Double_t py = 0.;
     Double_t pz = 0.;
 
     Double_t vx ;
     Double_t vy ;
     Double_t vz ;
 
     if(fGamOn){
  
       phi = gRandom->Uniform(0,360)* TMath::DegToRad();
       p =  gRandom->Uniform(fPMin,fPMax);
       theta = acos(gRandom->Uniform(cos(fThetaMin* TMath::DegToRad()),cos(fThetaMax* TMath::DegToRad())));
 
       pz = p*TMath::Cos(theta);
       py = p*TMath::Sin(theta)*TMath::Sin(phi);
       px = p*TMath::Sin(theta)*TMath::Cos(phi);
 
       vx = part->Vx(); 
       vy = part->Vy();
       vz = part->Vz();
       
       
       cout<<" Secondary vertex for stopping Xi - gammas "<<vx<<" "
           <<vy<<" "<<vz<<endl;
       
       // Give track to PrimaryGenerator
       
       
       primGen->AddTrack(22, px, py, pz, vx, vy,vz);
       
       
       
     }else
       {
         px= part->Px();
         py= part->Py();
         pz= part->Pz();
         
         
         vx = part->Vx(); 
         vy = part->Vy();
         vz = part->Vz();
         
         
         cout<<" Secondary vertex for stopping Xi - "<<pdgType<<" "<<vx<<" "
             <<vy<<" "<<vz<<endl;
         
         // Give track to PrimaryGenerator
         
         
         primGen->AddTrack(pdgType, px, py, pz, vx, vy,vz);
       }
     
     
   }
   
   
   return kTRUE;
 }