PndPythia8Direct.cxx

Go to the documentation of this file.

#include <math.h>
#include "TROOT.h"
#include "Pythia.h"
#include "FairPrimaryGenerator.h"
//#include "FairGenerator.h"

#include "PndPythia8Direct.h"

using namespace Pythia8;

// -----   Default constructor   -------------------------------------------
PndPythia8Direct::PndPythia8Direct()
{
  fPar_it = 0;
  fMom = 15.;
  fUseRandom1 = kFALSE;
  fUseRandom3 = kTRUE;
}
// -------------------------------------------------------------------------

// -----   Default constructor   -------------------------------------------
Bool_t PndPythia8Direct::Init()
{
  if (fUseRandom1) fRandomEngine = new PyTr1Rng();
  if (fUseRandom3) fRandomEngine = new PyTr3Rng();

  fPythia.setRndmEnginePtr(fRandomEngine);

  // Set Parameters
  for(Int_t jj=0; jj<fPar_it; jj++) {
    fPythia.readString(fParams[jj]);
    cout<<"fPythia.readString(\""<<fParams[jj]<<"\")"<<endl;
  }
  cout<<"Beam Momentum "<<fMom<<endl;

  // fPythia.init(-2212, 2212, 0., 0., fMom, 0., 0., 0.);
  // Settings according to the obsolete init function
  fPythia.settings.mode("Beams:idA",  -2212);
  fPythia.settings.mode("Beams:idB",   2212);
  fPythia.settings.mode("Beams:frameType",  3);
  fPythia.settings.parm("Beams:pxA",    0.);
  fPythia.settings.parm("Beams:pyA",    0.);
  fPythia.settings.parm("Beams:pzA",    fMom);
  fPythia.settings.parm("Beams:pxB",    0.);
  fPythia.settings.parm("Beams:pyB",    0.);
  fPythia.settings.parm("Beams:pzB",    0.);
  fPythia.init();

  return kTRUE;
}
// -------------------------------------------------------------------------


// -----   Destructor   ----------------------------------------------------
PndPythia8Direct::~PndPythia8Direct()
{
}
// -------------------------------------------------------------------------

// -----   Passing the event   ---------------------------------------------
Bool_t PndPythia8Direct::ReadEvent(FairPrimaryGenerator* cpg)
{
  Int_t npart = 0;
  while(npart == 0)
  {
    fPythia.next();
    for(int i=0; i<fPythia.event.size(); i++)
    {
      if(fPythia.event[i].isFinal())
      {
        npart++;
      };
    };
  };

  for(Int_t ii=0; ii<fPythia.event.size(); ii++)
    if(fPythia.event[ii].isFinal())
    {
      cpg->AddTrack((Int_t)fPythia.event[ii].id(),
                    (Double_t)fPythia.event[ii].px(),(Double_t)fPythia.event[ii].py(),(Double_t)fPythia.event[ii].pz(),
                    0.,0.,0.);
    };

  return kTRUE;
}
// -------------------------------------------------------------------------
void PndPythia8Direct::SetParameters(char* par)
{
  fParams[fPar_it] = par;
  fPar_it++;
}

// -------------------------------------------------------------------------
void PndPythia8Direct::Print() {
  fPythia.settings.listAll();
}
// -------------------------------------------------------------------------

ClassImp(PndPythia8Direct)