InitialState Class Reference

#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 fDatabase, 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 fDatabase.

                          {
    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(), alignCSCRings::e, fDatabase, GetDecayTime(), and edm::shift.

Referenced by gen::Hydjet2Hadronizer::generatePartonsAndHadronize().

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

virtual double InitialState::GetWeakDecayLimit ( )

pure virtual

Implemented in gen::Hydjet2Hadronizer.

virtual bool InitialState::RunDecays ( )

pure virtual

Implemented in gen::Hydjet2Hadronizer.

Member Data Documentation

DatabasePDG* InitialState::fDatabase

protected

Definition at line 16 of file InitialState.h.

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