Page Menu
Home
HEPForge
Search
Configure Global Search
Log In
Files
F10881559
No One
Temporary
Actions
View File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Flag For Later
Size
2 KB
Subscribers
None
View Options
Index: trunk/src/pythiadecayer.cpp
===================================================================
--- trunk/src/pythiadecayer.cpp (revision 140)
+++ trunk/src/pythiadecayer.cpp (revision 141)
@@ -1,73 +1,76 @@
#include "pythiadecayer.h"
#include "reportingUtils.h"
+#include "starlightconfig.h"
using namespace Pythia8;
+
+
pythiaDecayer::pythiaDecayer() :
- _pythia()
+ _pythia(PYTHIA8_SETTINGS_DIR)
{}
/*
pythiaDecayer::pythiaDecayer(const pythiaDecayer &obj) :
_pythia(obj._pythia)
{}
*/
pythiaDecayer::~pythiaDecayer()
{}
/*
pythiaDecayer& pythiaDecayer::operator=(const pythiaDecayer &other)
{
if (this != &other)
{
_pythia = other._pythia;
}
return *this;
}
*/
void pythiaDecayer::init()
{
_pythia.readString("ProcessLevel:all = off");
_pythia.readString("Standalone:allowResDec = on");
_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());
}
upcEvent pythiaDecayer::execute()
{
upcEvent slEvent;
Event &pyEvent = _pythia.event;
_pythia.forceTimeShower(1, 2, 100000.0);
// pyEvent.list();
if(!_pythia.next())
{
printWarn << "Pythia::next() failed" << std::endl;
return upcEvent();
}
for(int i = 0; i < pyEvent.size(); ++i)
{
Particle p = pyEvent[i];
starlightParticle slPart(p.px(), p.py(), p.pz(), p.e(), p.mass(), 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;
}
\ No newline at end of file
Index: trunk/starlightconfig.h.in
===================================================================
--- trunk/starlightconfig.h.in (revision 140)
+++ trunk/starlightconfig.h.in (revision 141)
@@ -1,6 +1,8 @@
// the configured options and settings for Tutorial
#define Starlight_VERSION_MAJOR @Starlight_VERSION_MAJOR@
#define Starlight_VERSION_MINOR @Starlight_VERSION_MINOR@
+#define PYTHIA8_SETTINGS_DIR "@PYTHIA8_SETTINGS_DIR@"
+
#cmakedefine ENABLE_PYTHIA
#cmakedefine ENABLE_PYTHIA6
#cmakedefine ENABLE_DPMJET
File Metadata
Details
Attached
Mime Type
text/x-diff
Expires
Sat, May 3, 6:32 AM (1 d, 17 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4983061
Default Alt Text
(2 KB)
Attached To
rSTARLIGHTSVN starlightsvn
Event Timeline
Log In to Comment