#include <PndGiBuuGenerator.h>

Inheritance diagram for PndGiBuuGenerator:

Public Member Functions
	PndGiBuuGenerator (const Char_t fileName, const Char_t pidDataFileName)

virtual	~PndGiBuuGenerator ()

virtual Bool_t	ReadEvent (FairPrimaryGenerator *primGen)

void	RegisterDecayer (Int_t pdgCode, PndGiBuuDecayerAbs *decayer)

Protected Member Functions
Bool_t	Init ()

Private Member Functions
	PndGiBuuGenerator ()

void	FillPidMap ()

Int_t	GetPdgParticleId (Int_t GiBuuId, Int_t charge)

void	WriteoutDecayParticle (Int_t oldEventNr, Int_t oldRunId, FairPrimaryGenerator *primGen)

void	CloseInput ()
	PDG code of particle. More...

	ClassDef (PndGiBuuGenerator, 1)

Private Attributes
std::map< std::pair< Int_t, Int_t >, Int_t >	fMapGiBuuToPdg

Int_t	fEvent

Int_t	fRunId
	Event number. More...

TString	fFileName
	Run number in GiBuu. More...

TString	fPidDataFileName
	Input file name. More...

std::ifstream *	fInputAsciiFile
	Input Pid data file name. More...

std::map< Int_t, PndGiBuuDecayerAbs * >	fDecayerMap

Double_t	fPx

Double_t	fPy

Double_t	fPz

Int_t	fPdg
	Momentum of particle. More...

Detailed Description

Definition at line 30 of file PndGiBuuGenerator.h.

Constructor & Destructor Documentation

PndGiBuuGenerator::PndGiBuuGenerator	(	const Char_t *	fileName,
		const Char_t *	pidDataFileName
	)

Standard constructor

Parameters

fileName	The input root file name
pidDataFileName	Translation between GiBuu particle Id and Pdg particle Id file

Definition at line 11 of file PndGiBuuGenerator.cxx.

                                                                                           :
                 fEvent(0), fRunId(0), fFileName(fileName), fPidDataFileName(pidDataFileName), fInputAsciiFile(0), fPx(0), fPy(0), fPz(0), fPdg(0)
 {
         std::cout << "Constructor Called!" << std::endl;
 }

PndGiBuuGenerator::~PndGiBuuGenerator ( )

virtual

Destructor

Definition at line 20 of file PndGiBuuGenerator.cxx.

References CloseInput(), and fDecayerMap.

                                       {
   CloseInput();
   for (std::map<Int_t, PndGiBuuDecayerAbs*>::const_iterator iter = fDecayerMap.begin(); iter != fDecayerMap.end(); iter++){
           delete(iter->second);
   }
 }

PndGiBuuGenerator::PndGiBuuGenerator ( )

inlineprivate

Default constructor (should not be used)

Definition at line 60 of file PndGiBuuGenerator.h.

60 {};

Member Function Documentation

PndGiBuuGenerator::ClassDef	(	PndGiBuuGenerator	,
		1
	)

private

void PndGiBuuGenerator::CloseInput ( )

private

PDG code of particle.

Definition at line 137 of file PndGiBuuGenerator.cxx.

References fInputAsciiFile.

Referenced by ~PndGiBuuGenerator().

                                    {
 
         fInputAsciiFile->close();
         //delete fInputAsciiFile;
         fInputAsciiFile = NULL;
 
 }

void PndGiBuuGenerator::FillPidMap ( )

private

Definition at line 103 of file PndGiBuuGenerator.cxx.

References fMapGiBuuToPdg, fPidDataFileName, and PidFile.

Referenced by Init().

 {
         std::ifstream* PidFile = new std::ifstream(fPidDataFileName.Data(), std::ifstream::in);
         if (0 == PidFile){
                 Fatal("PndGiBuuGenerator::FillPidMap", "No pid file!");
         }
         Int_t giPid = 0;
         Int_t giCharge = 0;
         Int_t pdgPid = 0;
 
 //      std::cout << "Reading In PidFile: " << fPidDataFileName.Data() <<  " " << PidFile << " " << PidFile->eof() << std::endl;
 
         while (!(PidFile->eof())){
                 *PidFile >> giPid >> giCharge >> pdgPid;
                 std::pair<Int_t, Int_t> GiBuuPid(giPid, giCharge);
                 fMapGiBuuToPdg[GiBuuPid] = pdgPid;
 //              std::cout << giPid << " " << giCharge << " = " << pdgPid << std::endl;
         }
         PidFile->close();
         PidFile = 0;
 
 }

Int_t PndGiBuuGenerator::GetPdgParticleId	(	Int_t	GiBuuId,
		Int_t	charge
	)

private

Definition at line 126 of file PndGiBuuGenerator.cxx.

References fMapGiBuuToPdg.

Referenced by ReadEvent().

 {
         std::pair<Int_t, Int_t> GiBuuPid(GiBuuId, charge);
         if (fMapGiBuuToPdg.count(GiBuuPid) > 0){
                 return fMapGiBuuToPdg[GiBuuPid];
         } else {
                 return 0;
         }
 }

Bool_t PndGiBuuGenerator::Init ( )

protected

Definition at line 28 of file PndGiBuuGenerator.cxx.

References fFileName, FillPidMap(), and fInputAsciiFile.

 {
 
         fInputAsciiFile = new std::ifstream(fFileName.Data(), std::ifstream::in);
         if (0 == fInputAsciiFile) {
                 Fatal("PndGiBuuGenerator", "Cannot open ASCII input file.");
         }
 
         FillPidMap();
 
         return kTRUE;
 }

Bool_t PndGiBuuGenerator::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 FairPrimaryGenerator

Definition at line 46 of file PndGiBuuGenerator.cxx.

References Double_t, energy, fEvent, fInputAsciiFile, fPdg, fPx, fPy, fPz, fRunId, GetPdgParticleId(), and WriteoutDecayParticle().

 {
 
         Int_t GiBuuPid = 0, GiBuuCharge = 0;
         //Int_t runId = 0; //[R.K. 01/2017] unused variable?
         Int_t motherId1 = 0;
         Int_t motherId2 = 0;
         Int_t dummyInt = 0;
         Double_t mass = 0.;
         Double_t energy = 0.;
         Int_t oldEventNr = fEvent;
         Int_t oldRunId = fRunId;
 
         Double_t dummyDouble = 0.;
 
         std::cout << "ReadEvent: " << fEvent << std::endl;
 
         if (fEvent > 0){
                 WriteoutDecayParticle(oldEventNr, oldRunId, primGen);
         } else {
                 oldEventNr = 0;
                 oldRunId = 0;
                 fEvent = 0;
                 fRunId = 0;
         }
 
         while (fEvent == oldEventNr && !(fInputAsciiFile->eof())){
                 *fInputAsciiFile >> GiBuuPid >> GiBuuCharge >> motherId1 >> motherId2 >> dummyInt >> mass >> dummyInt >> energy >> fPx >> fPy >> fPz >> dummyDouble >> fEvent >> fRunId >> dummyDouble;
                 if (0 == oldEventNr){
                         oldEventNr = fEvent;
                         oldRunId = fRunId;
                 }
                 fPdg = GetPdgParticleId(GiBuuPid, GiBuuCharge);
                 WriteoutDecayParticle(oldEventNr, oldRunId, primGen);
         }
 
         return kTRUE;
 }

void PndGiBuuGenerator::RegisterDecayer	(	Int_t	pdgCode,
		PndGiBuuDecayerAbs *	decayer
	)

inline

Definition at line 53 of file PndGiBuuGenerator.h.

References fDecayerMap.

53 { fDecayerMap[pdgCode] = decayer; }

PndGiBuuGenerator::fDecayerMap

std::map< Int_t, PndGiBuuDecayerAbs * > fDecayerMap

Definition: PndGiBuuGenerator.h:74

void PndGiBuuGenerator::WriteoutDecayParticle	(	Int_t	oldEventNr,
		Int_t	oldRunId,
		FairPrimaryGenerator *	primGen
	)

private

Definition at line 85 of file PndGiBuuGenerator.cxx.

References fDecayerMap, fEvent, fPdg, fPx, fPy, fPz, fRunId, i, X, Y, and Z.

Referenced by ReadEvent().

 {
         if (fEvent == oldEventNr && fRunId == oldRunId) {
                 if (fDecayerMap.count(fPdg) > 0) {
                         std::vector<PndGiBuuTrack> tracks = fDecayerMap[fPdg]->DecayTrack(PndGiBuuTrack(fPdg, fPx, fPy, fPz, 0., 0., 0.));
                         for (size_t i = 0; i < tracks.size(); i++) {
                                 std::cout << "Writing decayed particles: " << tracks[i].GetPdgId() << std::endl;
                                 primGen->AddTrack(tracks[i].GetPdgId(),
                                                 tracks[i].GetMomentum().X(), tracks[i].GetMomentum().Y(), tracks[i].GetMomentum().Z(),
                                                 tracks[i].GetVertex().X(), tracks[i].GetVertex().Y(), tracks[i].GetVertex().Z());
                         }
                 } else {
                         std::cout << "Writing particle direct: " << fPdg << std::endl;
                         primGen->AddTrack(fPdg, fPx, fPy, fPz, 0., 0., 0.);
                 }
         }
 }