Page MenuHomeHEPForge
Files F9501329

DecayHandler.cc
No OneTemporary
Actions

DecayHandler.cc
View Options

// -*- C++ -*-
//
// DecayHandler.cc is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 1999-2007 Leif Lonnblad
//
// ThePEG is licenced under version 2 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
//
// This is the implementation of the non-inlined, non-templated member
// functions of the DecayHandler class.
//
#include "DecayHandler.h"
#include "Hint.h"
#include "ThePEG/PDT/ParticleData.h"
#include "ThePEG/PDT/DecayMode.h"
#include "ThePEG/PDT/Decayer.h"
#include "ThePEG/EventRecord/Particle.h"
#include "ThePEG/EventRecord/Step.h"
#include "ThePEG/EventRecord/Collision.h"
#include "ThePEG/Interface/Parameter.h"
#include "ThePEG/Interface/ClassDocumentation.h"
#include "ThePEG/Persistency/PersistentOStream.h"
#include "ThePEG/Persistency/PersistentIStream.h"
#ifdef ThePEG_TEMPLATES_IN_CC_FILE
// #include "DecayHandler.tcc"
#endif
using namespace ThePEG;
DecayHandler::~DecayHandler() {}
void DecayHandler::
handle(EventHandler &, const tPVector & tagged,
const Hint &) ThePEG_THROW_SPEC((Veto, Stop, Exception)) {
// First go through to see which of the tagged particles should be
// decayed: Exit if none are found.
tPVector parents;
for ( int i = 0, N = tagged.size(); i < N; ++i )
if ( tagged[i] && !tagged[i]->data().stable() )
parents.push_back(tagged[i]);
if ( parents.empty() ) return;
// Create a new step, decay all particles and add their children in
// the new step.
for ( int i = 0, N = parents.size(); i < N; ++i )
performDecay(newStep()->find(parents[i]), *newStep());
}
void DecayHandler::
performDecay(tPPtr parent, Step & s) const throw(Veto, Exception) {
ParticleVector children = Decayer::DecayParticle(parent, s, maxLoop());
for ( int i = 0, N = children.size(); i < N; ++i )
if ( !children[i]->data().stable() ) performDecay(children[i], s);
}
void DecayHandler::persistentOutput(PersistentOStream & os) const {
os << theMaxLoop;
}
void DecayHandler::persistentInput(PersistentIStream & is, int) {
is >> theMaxLoop;
}
ClassDescription<DecayHandler> DecayHandler::initDecayHandler;
void DecayHandler::Init() {
static ClassDocumentation<DecayHandler> documentation
("This is the main class handling the decay of unstable particles. Note "
"that the actual decays will be performed by objects of the "
"ThePEG::Decayer class.");
static Parameter<DecayHandler,long> interfaceMaxLoop
("MaxLoop",
"The maximum number of attempts per event when selecting a decay channel.",
&DecayHandler::theMaxLoop, 100000, 100, 100000000, true, false, true);
}

File Metadata

Mime Type
text/x-c
Expires
Sun, Feb 23, 2:19 PM (2 h, 23 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4486525
Default Alt Text
DecayHandler.cc (2 KB)

Event Timeline