CbmPhsdGenerator.h

Go to the documentation of this file.

// -------------------------------------------------------------------------
// -----                CbmPHSDGenerator header file                  -----
// -----          Created 24/01/14, updated on 17/11/16  by V. Vovchenko           -----
// -------------------------------------------------------------------------
 
 
#ifndef CBMPHSDGENERATOR_H
#define CBMPHSDGENERATOR_H
 
 
#include "FairGenerator.h"
 
#include <fstream>
#include <map>
 
class TVirtualMCStack;
class FairPrimaryGenerator;
 
// This structure holds information about hadron from HSD output file
struct Hadron {
  bool init;  // Are parameters initialised?
  int id;     // HSD particle id
  int charge;
  int ISUB;               // Number of subsequent parallel run
  int IRUN;               // Number of run in subsequent run
  int globalEvent;        // = (ISUB-1)*IRUNS + IRUN
  double px, py, pz, p0;  // 4-momentum
  double b;               // impact parameter of event
};
 
 
class CbmPhsdGenerator : public FairGenerator {
 
public:
  CbmPhsdGenerator();
 
 
  CbmPhsdGenerator(const char* fileNameInput,
                   const char* fileNameBaryons,
                   const char* fileNameMesons);
 
  CbmPhsdGenerator(const char* fileNameInput, const char* fileNameDat);
 
 
  ~CbmPhsdGenerator();
 
 
  Bool_t ReadEvent(FairPrimaryGenerator* primGen);
 
  Bool_t ReadEvent300(FairPrimaryGenerator* primGen);
 
  Bool_t ReadEventDat(FairPrimaryGenerator* primGen);
 
  Bool_t SkipEvents(Int_t count);
 
private:
  Bool_t fReadDat;  
 
  FILE *fBaryonsFile, *fMesonsFile;  
 
  FILE* fDatFile;  
 
  std::map<Int_t, Int_t> fParticleTable;  
 
  const Char_t* fFileName;  
 
  // Necessary to have because of HSD output structure
  Hadron nextBaryon, nextMeson;  // Baryons and mesons read from next event
  int nextEvent;                 // Id of the next event
 
 
  // Data about HSD simulation, is read from HSD input file
  int Ap, Zp, At, Zt;
  int ISUBS;    // Number of subsequent parallel runs in HSD
  int IRUNS;    // Number of runs in one parallel batch in HSD
  double ekin;  // Kinetic energy per projectile nucleon in lab frame
 
 
  void ReadConversionTable();
 
  void ReadCollisionData(const char* fileNameInput);
 
  CbmPhsdGenerator(const CbmPhsdGenerator&);
  CbmPhsdGenerator& operator=(const CbmPhsdGenerator&);
 
  ClassDef(CbmPhsdGenerator, 1);
};
 
#endif