pythiadecayer.cpp

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file pythiadecayer.cpp

#include "pythiadecayer.h"
#include "reportingUtils.h"
#include "starlightconfig.h"
using namespace Pythia8;



pythiaDecayer::pythiaDecayer() :
    _pythia(PYTHIA8_SETTINGS_DIR)
{}
pythiaDecayer::~pythiaDecayer()
{}
void pythiaDecayer::init()
{
  _pythia.readString("ProcessLevel:all = off"); 
//  _pythia.readString("Standalone:allowResDec = on");//Option removed from Pythia8 JB05192015
  _pythia.readString("Next:numberShowEvent = 0");
  _pythia.init();                               
  _pythia.event.reset();
}

void pythiaDecayer::addParticle(const starlightParticle &p)
{
  
  Event &pyEvent = _pythia.event;
  int status = 23; // Outgoing particle from the hardest sub-process
  int col = 0;
  int acol = 0;
  int code = p.getPdgCode();

  pyEvent.append(code, status, col, acol, p.GetPx(), p.GetPy(), p.GetPz(), p.GetE(), p.M());
  
}

eXEvent pythiaDecayer::execute()
{
  eXEvent slEvent;
  
  Event &pyEvent = _pythia.event;
  _pythia.forceTimeShower(1, 2, 100000.0);
  if(!_pythia.next())
  {
    printWarn << "Pythia::next() failed" << std::endl;
    return eXEvent();
  }
  
  for(int i = 0; i < pyEvent.size(); ++i)
  {
    
      Particle p = pyEvent[i];
      starlightParticle slPart(p.px(), p.py(), p.pz(), p.e(), p.m(), p.idAbs()*(p.charge()<0?-1:1), p.charge(),
             p.xProd(), p.yProd(), p.zProd(), p.tProd(),
             p.mother1(), p.mother2(), p.daughter1(), p.daughter2(), p.status());
      slEvent.addParticle(slPart);
  }
  pyEvent.clear();
  pyEvent.reset(); 
  return slEvent;
  
}