#include <PndPythia6Generator.h>

Inheritance diagram for PndPythia6Generator:

Public Member Functions
	PndPythia6Generator ()

	PndPythia6Generator (const char *fileName)

virtual	~PndPythia6Generator ()

virtual Bool_t	ReadEvent (FairPrimaryGenerator *primGen)

void	SetVerbose (Int_t verb)

Private Member Functions
void	CloseInput ()
	Verbose Level. More...

	ClassDef (PndPythia6Generator, 1)

Private Attributes
const Char_t *	fFileName

FILE *	fInputFile
	Input file Name. More...

Int_t	fVerbose
	File. More...

Detailed Description

Definition at line 60 of file PndPythia6Generator.h.

Constructor & Destructor Documentation

PndPythia6Generator::PndPythia6Generator ( )

Default constructor without arguments should not be used.

Definition at line 17 of file PndPythia6Generator.cxx.

17 {}

PndPythia6Generator::PndPythia6Generator ( const char * fileName )

Standard constructor.

Parameters

fileName The input file name

Definition at line 23 of file PndPythia6Generator.cxx.

References fFileName, fInputFile, and fVerbose.

                                                              {
   fFileName  = fileName;
   fVerbose = 0;
   cout << "-I PndPythia6Generator: Opening input file " << fileName << endl;
   if ((fInputFile = fopen(fFileName,"r"))==NULL)
     //  fInputFile = new ifstream(fFileName);
     //  if ( ! fInputFile->is_open() ) 
     Fatal("PndPythia6Generator","Cannot open input file.");
   
   // fPDG=TDatabasePDG::Instance();
 }

PndPythia6Generator::~PndPythia6Generator ( )

virtual

Destructor.

Definition at line 39 of file PndPythia6Generator.cxx.

References CloseInput().

                                           {
   CloseInput();
 }

Member Function Documentation

PndPythia6Generator::ClassDef	(	PndPythia6Generator	,
		1
	)

private

PDG data base

void PndPythia6Generator::CloseInput ( )

private

Verbose Level.

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

Definition at line 111 of file PndPythia6Generator.cxx.

References fFileName, and fInputFile.

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

                                      {
   if ( fInputFile ) {
     //if ( fInputFile->is_open() ) {
     {
       cout << "-I PndPythia6Generator: Closing input file " 
            << fFileName << endl;
       //  fInputFile->close();
       
       fclose(fInputFile);
     }
     delete fInputFile;
     fInputFile = NULL;
   }
 }

Bool_t PndPythia6Generator::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 CbmrimaryGenerator

Definition at line 47 of file PndPythia6Generator.cxx.

References CloseInput(), fInputFile, fT, and fVerbose.

                                                                    {
 
   // Check for input file
   if (!fInputFile) {
     // if ( ! fInputFile->is_open() ) {
     cout << "-E PndPythia6Generator: Input file not open!" << endl;
     return kFALSE;
   }
   
   // Define event variable to be read from file
    Int_t ntracks = 0, eventID = 0, ncols = 0;
  
   // Define track variables to be read from file
   Int_t nLev = 0, pdgID = 0, nM1 = -1, nM2 = -1, nDF = -1, nDL = -1;
   Float_t fPx = 0., fPy = 0., fPz = 0., fM = 0., fE = 0.;
   Float_t fVx = 0., fVy = 0., fVz = 0., fT = 0.;
   
   // Read event header line from input file
 
   //Int_t max_nr = 0; //[R.K. 01/2017] unused variable?
   
   //Text_t buffer[200]; //[R.K. 01/2017] unused variable?
   ncols = fscanf(fInputFile,"%d\t%d", &eventID, &ntracks);
   
   if (ncols && ntracks>0) {
     
     if (fVerbose>0) cout << "Event number: " << eventID << "\tNtracks: " << ntracks << endl;
     
     for (Int_t ll=0; ll<ntracks; ll++) 
       {
         ncols = fscanf(fInputFile,"%d %d %d %d %d %d %f %f %f %f %f %f %f %f %f", &nLev, &pdgID, &nM1, &nM2, &nDF, &nDL, &fPx, &fPy, &fPz, &fE, &fM, &fVx, &fVy, &fVz, &fT);
         if (fVerbose>0) cout << nLev << "\t" << pdgID << "\t" << nM1 << "\t" << nM2 << "\t" << nDF << "\t" << nDL <<
           "\t" << fPx << "\t" << fPy << "\t" << fPz << "\t" << fE << "\t" << fM << "\t" << fVx << "\t" << fVy << "\t" << fVz << "\t" << fT <<  endl;
         if (nLev==1)  
           primGen->AddTrack(pdgID, fPx, fPy, fPz, fVx, fVy, fVz);
       }
   }
   else {
     cout << "-I PndPythia6Generator: End of input file reached " << endl;
     CloseInput();
     return kFALSE;
   }
   
   
   // If end of input file is reached : close it and abort run
   if ( feof(fInputFile) ) {
     cout << "-I PndPythia6Generator: End of input file reached " << endl;
     CloseInput();
     return kFALSE;
   }
   
   /*
     cout << "-I PndPythia6Generator: Event " << eventID << ",  vertex = ("
     << vx << "," << vy << "," << vz << ") cm,  multiplicity "
     << ntracks << endl;
   */
   
   return kTRUE;
 }