AtTPCXSReader.h

Go to the documentation of this file.

// -------------------------------------------------------------------------
// -----                    AtTPCXSReader header file                  -----
// -----               Created 03/07/18 by H. Alvarez Pol              -----
// -------------------------------------------------------------------------
 
#ifndef AtTPCXSREADER_H
#define AtTPCXSREADER_H
 
#include <FairGenerator.h>
 
#include <Rtypes.h>
#include <TString.h>
 
#include <vector>
 
class FairPrimaryGenerator;
class FairIon;
class FairParticle;
class TBuffer;
class TClass;
class TH2F;
class TMemberInspector;
 
class AtTPCXSReader : public FairGenerator {
 
public:
   AtTPCXSReader();
 
   AtTPCXSReader(const char *name, std::vector<Int_t> *z, std::vector<Int_t> *a, std::vector<Int_t> *q, Int_t mult,
                 std::vector<Double_t> *px, std::vector<Double_t> *py, std::vector<Double_t> *pz,
                 std::vector<Double_t> *mass);
 
   AtTPCXSReader &operator=(const AtTPCXSReader &) { return *this; }
 
   virtual ~AtTPCXSReader() = default;
 
   virtual Bool_t ReadEvent(FairPrimaryGenerator *primGen);
 
   void SetXSFileName(TString name = "xs_22Mgp_fusionEvaporation.txt") { fXSFileName = name; }
 
private:
   TString fXSFileName;
 
   static Int_t fgNIon;                 
   Int_t fMult;                         // Multiplicity per event
   std::vector<Double_t> fPx, fPy, fPz; // Momentum components [GeV] per nucleon
   std::vector<Double_t> Masses;        // Masses of the N products
   std::vector<Double_t> fExEnergy;     // Excitation energies of the products
   Double_t fVx, fVy, fVz;              // Vertex coordinates [cm]
   std::vector<FairIon *> fIon;         // Pointer to the FairIon to be generated
   std::vector<TString> fPType;
   std::vector<FairParticle *> fParticle;
   std::vector<Int_t> fQ;  // Electric charge [e]
   Double_t fBeamEnergy{}; // Residual beam energy for phase calculation
   Double_t fPxBeam{};
   Double_t fPyBeam{};
   Double_t fPzBeam{};
   std::vector<Double_t> fWm; // Total mass
 
   TH2F *fh_pdf{};
 
   ClassDef(AtTPCXSReader, 1)
};
 
#endif