Pythia8Generator.h

Go to the documentation of this file.

/********************************************************************************
 *    Copyright (C) 2014 GSI Helmholtzzentrum fuer Schwerionenforschung GmbH    *
 *                                                                              *
 *              This software is distributed under the terms of the             *
 *         GNU Lesser General Public Licence version 3 (LGPL) version 3,        *
 *                  copied verbatim in the file "LICENSE"                       *
 ********************************************************************************/
// -------------------------------------------------------------------------
// -----                  M. Al-Turany   June 2014                     -----
// -------------------------------------------------------------------------
 
#ifndef PNDP8GENERAtOR_H
#define PNDP8GENERAtOR_H 1
 
#include <FairGenerator.h>
 
#include <TROOT.h>
#include <TRandom1.h>
#include <TRandom3.h>
 
#include "Pythia.h"
 
class FairPrimaryGenerator;
using namespace Pythia8;
 
class PyTr1Rng : public RndmEngine {
public:
   PyTr1Rng() { rng = new TRandom1(gRandom->GetSeed()); };
   virtual ~PyTr1Rng(){};
 
   Double_t flat() { return rng->Rndm(); };
 
private:
   TRandom1 *rng; 
};
 
class PyTr3Rng : public RndmEngine {
public:
   PyTr3Rng() { rng = new TRandom3(gRandom->GetSeed()); };
   virtual ~PyTr3Rng(){};
 
   Double_t flat() { return rng->Rndm(); };
 
private:
   TRandom3 *rng; 
};
 
class Pythia8Generator : public FairGenerator {
public:
   Pythia8Generator();
 
   virtual ~Pythia8Generator();
 
   Bool_t ReadEvent(FairPrimaryGenerator *);
   void SetParameters(char *);
   void Print(); 
 
   virtual Bool_t Init(); 
 
   void SetMom(Double_t mom) { fMom = mom; };
   void SetId(Double_t id) { fId = id; };
   void SetHNLId(Int_t id) { fHNL = id; };
   void UseRandom1()
   {
      fUseRandom1 = kTRUE;
      fUseRandom3 = kFALSE;
   };
   void UseRandom3()
   {
      fUseRandom1 = kFALSE;
      fUseRandom3 = kTRUE;
   };
   void GetPythiaInstance(int);
 
private:
   Pythia fPythia;            
   RndmEngine *fRandomEngine; 
 
protected:
   Double_t fMom;      // proton momentum
   Int_t fHNL;         // HNL ID
   Int_t fId;          // target type
   Bool_t fUseRandom1; // flag to use TRandom1
   Bool_t fUseRandom3; // flag to use TRandom3 (default)
 
   ClassDef(Pythia8Generator, 1);
};
 
#endif /* !PNDP8GENERAtOR_H */