#include <InitialState.h>

Inheritance diagram for InitialState:

Public Member Functions
virtual void	Evolve (List_t &secondaries, ParticleAllocator &allocator, double weakDecayLimit)

virtual double	GetWeakDecayLimit ()=0

	InitialState ()

virtual bool	RunDecays ()=0

virtual	~InitialState ()

Protected Attributes
DatabasePDG *	fDatabase

Detailed Description

Definition at line 14 of file InitialState.h.

Constructor & Destructor Documentation

InitialState::InitialState ( )

inline

Definition at line 18 of file InitialState.h.

References DatabasePDG::LoadData(), DatabasePDG::SetMassRange(), and DatabasePDG::SetWidthRange().

                  {
     fDatabase = new DatabasePDG();
     fDatabase->LoadData();
     fDatabase->SetMassRange(0.0, 200.);
     fDatabase->SetWidthRange(0., 10.);
   };

virtual InitialState::~InitialState ( )

inlinevirtual

Definition at line 24 of file InitialState.h.

References Evolve(), fDatabase, GetWeakDecayLimit(), and RunDecays().

                           {
     delete fDatabase;
   };

Member Function Documentation

void InitialState::Evolve	(	List_t &	secondaries,
		ParticleAllocator &	allocator,
		double	weakDecayLimit
	)

virtual

Definition at line 13 of file InitialState.cc.

References Decay(), MillePedeFileConverter_cfg::e, fDatabase, GetDecayTime(), and edm::shift.

Referenced by gen::Hydjet2Hadronizer::generatePartonsAndHadronize(), and ~InitialState().

                                                                                                   {
   // Particle indexes are set for primaries already from InitialStateHydjet::Initialize()
    
   // particle list iterators
   LPIT_t it;
   LPIT_t e;
 
   // Decay loop
   // Note that the decay products are always added at the end of list so the unstable products are
   // decayed when the iterator reaches them (used for cascade decays)
   for(it = secondaries.begin(), e = secondaries.end(); it != e; ++it) {
 
     // if the decay procedure was applied already skip ... (e.g. particles from pythia history information)
     if(it->GetDecayed()) {
       continue;
     }
 
     // generate the decay time; if particle is stable or set to stable decayTime=0
     double decayTime = GetDecayTime(*it, weakDecayLimit);
     it->SetLastInterTime(it->T() + decayTime);
     TVector3 shift(it->Mom().BoostVector());
     shift *= decayTime; 
     it->SetDecayed();
 
     // if decayTime>0 then apply the decay procedure (only 2 or 3 body decays)
     if(decayTime > 0.) {
 
       it->Pos(it->Pos() += TLorentzVector(shift, 0.));
       it->T(it->T() + decayTime);
       Decay(secondaries, *it, allocator, fDatabase);
     }
     // if particle is stable just continue
   }
 }