#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 () noexcept(false)

Protected Attributes
DatabasePDG *	fDatabase

Detailed Description

Definition at line 14 of file InitialState.h.

Constructor & Destructor Documentation

◆ InitialState()

InitialState::InitialState ( )

inline

Definition at line 19 of file InitialState.h.

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

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

◆ ~InitialState()

virtual InitialState::~InitialState ( )

inlinevirtualnoexcept

Definition at line 25 of file InitialState.h.

25 { delete fDatabase; };

References fDatabase.

Member Function Documentation

◆ Evolve()

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

virtual

Definition at line 13 of file InitialState.cc.

                                                                                                   {
   // 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
   }
 }