#include <EvtInclusiveDecay.h>

Inheritance diagram for EvtInclusiveDecay:

Collaboration diagram for EvtInclusiveDecay:

Public Member Functions
	EvtInclusiveDecay (const std::string &name, ISvcLocator *pSvcLocator)

virtual	~EvtInclusiveDecay ()

StatusCode	initialize ()

StatusCode	execute ()

StatusCode	finalize ()

std::string	xmlpath (void)

virtual StatusCode	sysInitialize () override
	Override sysInitialize. More...

virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies. More...

ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`. More...

const ServiceHandle< StoreGateSvc > &	evtStore () const
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`. More...

const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard `StoreGateSvc/DetectorStore` Returns (kind of) a pointer to the `StoreGateSvc`. More...

virtual StatusCode	sysStart () override
	Handle START transition. More...

virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles. More...

virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles. More...

Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T > &t)

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc="none")
	Declare a new Gaudi property. More...

void	updateVHKA (Gaudi::Details::PropertyBase &)

MsgStream &	msg () const

MsgStream &	msg (const MSG::Level lvl) const

bool	msgLvl (const MSG::Level lvl) const

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution More...

std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void >	renounce (T &h)

void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed. More...

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
	specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
	specialization for handling Gaudi::Property<SG::VarHandleBase> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
	specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...

Private Attributes
DataObjIDColl	m_extendedExtraObjects

StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default) More...

StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default) More...

std::vector< SG::VarHandleKeyArray * >	m_vhka

bool	m_varHandleArraysDeclared

Features for derived classes to use internally
ServiceHandle< IAthRNGSvc >	m_rndmSvc {this, "RndmSvc", "AthRNGSvc"}

IntegerProperty	m_dsid {this, "Dsid", 999999}

IntegerProperty	m_randomSeed {this, "RandomSeed", 1234567, "Random seed for the built-in random engine"}
	Seed for random number engine. More...

McEventCollection *	m_mcEvtColl {}

EvtInclusiveAtRndmGen *	m_evtAtRndmGen {}

EvtGen *	m_myEvtGen {}

std::string	m_pdtFile

std::string	m_decayFile

std::string	m_userDecayFile

std::string	m_randomStreamName

std::string	m_inputKeyName

std::string	m_outputKeyName

bool	m_readExisting

bool	m_prohibitFinalStateDecay

bool	m_prohibitReDecay

bool	m_prohibitUnDecay

bool	m_prohibitRemoveSelfDecay

std::vector< int >	m_blackList

std::set< int >	m_blackListSet

bool	m_allowAllKnownDecays

bool	m_allowDefaultBDecays

std::vector< int >	m_whiteList

std::set< int >	m_whiteListSet

bool	m_printHepMCBeforeEvtGen

bool	m_printHepMCAfterEvtGen

bool	m_printHepMCHighlighted

bool	m_printHepMCHighLightTopLevelDecays

bool	m_checkDecayTree

bool	m_checkDecayChannels

std::map< int, long >	m_noDecayChannels

int	m_nRepeatedDecays

int	m_maxNRepeatedDecays

bool	m_applyUserSelection

bool	m_userSelRequireOppositeSignedMu

double	m_userSelMu1MinPt

double	m_userSelMu2MinPt

double	m_userSelMu1MaxEta

double	m_userSelMu2MaxEta

double	m_userSelMinDimuMass

double	m_userSelMaxDimuMass

bool	m_isfHerwig

bool	m_setVMtransversePol

void	reseedRandomEngine (const std::string &streamName, const EventContext &ctx)

CLHEP::HepRandomEngine *	getRandomEngine (const std::string &streamName, unsigned long int randomSeedOffset, const EventContext &ctx) const

CLHEP::HepRandomEngine *	getRandomEngineDuringInitialize (const std::string &streamName, unsigned long int randomSeedOffset, unsigned int conditionsRun=1, unsigned int lbn=1) const

StatusCode	traverseDecayTree (HepMC::GenParticlePtr p, bool isToBeRemoved, std::set< HepMC::GenVertexPtr > &visited, std::set< int > &toBeDecayed)

void	removeDecayTree (HepMC::GenEvent *hepMC, HepMC::GenParticlePtr p)

void	decayParticle (HepMC::GenEvent *hepMC, HepMC::GenParticlePtr p)

void	addEvtGenDecayTree (HepMC::GenEvent hepMC, HepMC::GenParticlePtr part, EvtParticle evtPart, EvtVector4R treeStart, double momentumScaleFactor=1.0)

bool	isToBeDecayed (HepMC::ConstGenParticlePtr p, bool doCrossChecks)

bool	isDefaultB (const int pId) const

bool	passesUserSelection (HepMC::GenEvent *hepMC)

double	invMass (HepMC::ConstGenParticlePtr p1, HepMC::ConstGenParticlePtr p2)

void	printHepMC (HepMC::GenEvent hepMC, std::set< int > barcodeList=nullptr)

unsigned int	printTree (HepMC::GenParticlePtr p, std::set< HepMC::GenVertexPtr > &visited, int level, std::set< int > *barcodeList=0)

std::string	pdgName (HepMC::ConstGenParticlePtr p, bool statusHighlighting=false, std::set< int > *barcodeList=nullptr)

Utility event-mangling functions
Todo: Replace with HepMC units when available
ServiceHandle< IPartPropSvc >	m_ppSvc {this, "PartPropSvc", "PartPropSvc"}
	Handle on the particle property service. More...

SG::ReadHandleKey< McEventCollection >	m_mcevents_const { this, "McEventKey", "GEN_EVENT", "StoreGate key of the MC event collection" }
	Const handle to the MC event collection. More...

void	GeVToMeV (HepMC::GenEvent *evt)
	Scale event energies/momenta by x 1000. More...

void	MeVToGeV (HepMC::GenEvent *evt)
	Scale event energies/momenta by x 1/1000. More...

void	cmTomm (HepMC::GenEvent *evt)
	Scale event lengths by x 10. More...

void	mmTocm (HepMC::GenEvent *evt)
	Scale event lengths by x 1/10. More...

Detailed Description

Definition at line 58 of file EvtInclusiveDecay.h.

Member Typedef Documentation

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< Algorithm > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ EvtInclusiveDecay()

EvtInclusiveDecay::EvtInclusiveDecay	(	const std::string &	name,
		ISvcLocator *	pSvcLocator
	)

Definition at line 58 of file EvtInclusiveDecay.cxx.

                                                                                    :
   GenBase( name, pSvcLocator ) {
  
   // Basic EvtGen configuration: decay and particle definition files, random number stream
   declareProperty("pdtFile", m_pdtFile = "inclusive.pdt");
   declareProperty("decayFile", m_decayFile = "2014inclusive.dec");
   declareProperty("userDecayFile", m_userDecayFile = "");
   declareProperty("randomStreamName", m_randomStreamName = "EVTGEN");
   declareProperty("inputKeyName", m_inputKeyName = "GEN_EVENT");
   declareProperty("outputKeyName",m_outputKeyName = "GEN_EVENT_EVTGEN");
   declareProperty("readExisting",m_readExisting=false);
  
   // Selection of particles to be decayed
   declareProperty("prohibitFinalStateDecay", m_prohibitFinalStateDecay=false);
   declareProperty("prohibitReDecay", m_prohibitReDecay=false);
   declareProperty("prohibitUnDecay", m_prohibitUnDecay=true);
   declareProperty("prohibitRemoveSelfDecay", m_prohibitRemoveSelfDecay=false);
   declareProperty("blackList",m_blackList);
   declareProperty("allowAllKnownDecays", m_allowAllKnownDecays=true);
   declareProperty("allowDefaultBDecays", m_allowDefaultBDecays=true);
   declareProperty("whiteList",m_whiteList);
  
   // Level of output
   declareProperty("printHepMCBeforeEvtGen", m_printHepMCBeforeEvtGen=false);
   declareProperty("printHepMCAfterEvtGen", m_printHepMCAfterEvtGen=false);
   declareProperty("printHepMCHighlighted", m_printHepMCHighlighted=true);
   declareProperty("printHepMCHighLightTopLevelDecays", m_printHepMCHighLightTopLevelDecays=true);
  
   // Optional checks
   declareProperty("checkDecayTree", m_checkDecayTree=false);
   declareProperty("checkDecayChannels", m_checkDecayChannels=false);
  
   // Repeated decays
   declareProperty("maxNRepeatedDecays", m_maxNRepeatedDecays=1);
  
   // User selection
   declareProperty("applyUserSelection", m_applyUserSelection=false);
   declareProperty("userSelRequireOppositeSignedMu", m_userSelRequireOppositeSignedMu=true);
   declareProperty("userSelMu1MinPt", m_userSelMu1MinPt=0.);
   declareProperty("userSelMu2MinPt", m_userSelMu2MinPt=0.);
   declareProperty("userSelMu1MaxEta", m_userSelMu1MaxEta=102.5);
   declareProperty("userSelMu2MaxEta", m_userSelMu2MaxEta=102.5);
   declareProperty("userSelMinDimuMass", m_userSelMinDimuMass=0.);
   declareProperty("userSelMaxDimuMass", m_userSelMaxDimuMass=-1.); // set to negative to not apply cut
   declareProperty("isfHerwig", m_isfHerwig=false);
   declareProperty("setVMtransversePol", m_setVMtransversePol=false);
  
   // We have decided to blacklist Tau decays because we are not sure whether the polarization
   // would be properly passed to EvtGen
   m_blackList.push_back(15);
 }

◆ ~EvtInclusiveDecay()

EvtInclusiveDecay::~EvtInclusiveDecay ( )

virtual

Definition at line 112 of file EvtInclusiveDecay.cxx.

                                       {
   delete m_myEvtGen;
   delete m_evtAtRndmGen;
 }

Member Function Documentation

◆ addEvtGenDecayTree()

void EvtInclusiveDecay::addEvtGenDecayTree	(	HepMC::GenEvent *	hepMC,
		HepMC::GenParticlePtr	part,
		EvtParticle *	evtPart,
		EvtVector4R	treeStart,
		double	momentumScaleFactor = `1.0`
	)

private

Definition at line 532 of file EvtInclusiveDecay.cxx.

                                                                                                                     {
   if(evtPart->getNDaug()!=0) {
     // Add decay vertex, starting from production vertex of particle
     double ct=(evtPart->getDaug(0)->get4Pos()).get(0)+treeStart.get(0);
     double x=(evtPart->getDaug(0)->get4Pos()).get(1)+treeStart.get(1);
     double y=(evtPart->getDaug(0)->get4Pos()).get(2)+treeStart.get(2);
     double z=(evtPart->getDaug(0)->get4Pos()).get(3)+treeStart.get(3);
  
     HepMC::GenVertexPtr end_vtx = HepMC::newGenVertexPtr(HepMC::FourVector(x,y,z,ct));
  
     hepMC->add_vertex(end_vtx);
     end_vtx->add_particle_in(part);
  
     // Add decay daughter with their own decay trees
     for(uint it=0; it<evtPart->getNDaug(); it++) {
       double e=(evtPart->getDaug(it)->getP4Lab()).get(0) * momentumScaleFactor;
       double px=(evtPart->getDaug(it)->getP4Lab()).get(1) * momentumScaleFactor;
       double py=(evtPart->getDaug(it)->getP4Lab()).get(2) * momentumScaleFactor;
       double pz=(evtPart->getDaug(it)->getP4Lab()).get(3) * momentumScaleFactor;
       int id=EvtPDL::getStdHep(evtPart->getDaug(it)->getId());
       int status=1;
       if(evtPart->getDaug(it)->getNDaug() != 0) status=2;
       HepMC::GenParticlePtr daughter = HepMC::newGenParticlePtr(HepMC::FourVector(px,py,pz,e),id,status);
       end_vtx->add_particle_out(daughter);
       addEvtGenDecayTree(hepMC, daughter, evtPart->getDaug(it), treeStart, momentumScaleFactor);
     }
   }
 }

◆ ATLAS_NOT_CONST_THREAD_SAFE() [1/2]

HepMC::GenEvent* event GenBase::ATLAS_NOT_CONST_THREAD_SAFE ( )

inlineinherited

Access the current signal event (first in the McEventCollection)

Note: This function will make a new McEventCollection if there is not already a valid one and MakeMcEvent=True.

Definition at line 76 of file GenBase.h.

                                                       {
     if (events()->empty())
       ATH_MSG_ERROR("McEventCollection is empty during first event access");
     return *(events()->begin());
   }

◆ ATLAS_NOT_CONST_THREAD_SAFE() [2/2]

McEventCollection* events GenBase::ATLAS_NOT_CONST_THREAD_SAFE ( )

inherited

Access the current event's McEventCollection.

Note: This function will make a new McEventCollection if there is not already a valid one and MakeMcEvent=True.

◆ cmTomm()

void GenBase::cmTomm ( HepMC::GenEvent * evt )

protectedinherited

Scale event lengths by x 10.

Definition at line 81 of file GenBase.cxx.

                                        {
   for (HepMC::GenEvent::vertex_iterator vtx = evt->vertices_begin(); vtx != evt->vertices_end(); ++vtx) {
     const HepMC::FourVector fv((*vtx)->position().x() * 10,
                                (*vtx)->position().y() * 10,
                                (*vtx)->position().z() * 10,
                                (*vtx)->position().t() * 10);
     (*vtx)->set_position(fv);
   }
 }

◆ decayParticle()

void EvtInclusiveDecay::decayParticle	(	HepMC::GenEvent *	hepMC,
		HepMC::GenParticlePtr	p
	)

private

Definition at line 493 of file EvtInclusiveDecay.cxx.

                                                                                     {
 // TODO the printout below crashes with segfault for HepMC3 - needs expert action
   ATH_MSG_DEBUG("Decaying particle " << pdgName(part) << " " << part);
   if (msgLvl(MSG::VERBOSE)) HepMC::Print::line(std::cout,part);
  
   // Remove existing decay tree, if any, and flag particle as being decayed by EvtGen
   removeDecayTree(hepMC,part);
   part->set_status(HepMC::EVTGENUNDECAYEDSTATUS);
  
   // Create EvtGen version of part and have EvtGen decay it.
   // Since EvtGen uses GeV, convert particles momentum from MeV to GeV.
   int id = part->pdg_id();
   EvtId evtId=EvtPDL::evtIdFromStdHep(id);
   double en =(part->momentum()).e()/1000.;
   double px=(part->momentum()).px()/1000.;
   double py=(part->momentum()).py()/1000.;
   double pz=(part->momentum()).pz()/1000.;
   EvtVector4R evtP(en,px,py,pz);
   EvtParticle* evtPart = EvtParticleFactory::particleFactory(evtId,evtP);
   
   // set transverse polarization to vector mesons (relevant for coherent production of J/Psi etc in UPC)
   if(m_setVMtransversePol && (id==113 || id== 443 || id==100443 || id==553 || id==100553 || id==200553) )evtPart->setVectorSpinDensity();
   
   m_myEvtGen->generateDecay(evtPart);
   if (msgLvl(MSG::VERBOSE)) evtPart->printTree();
   double ct_s = part->production_vertex()->position().t();
   double x_s = part->production_vertex()->position().x();
   double y_s = part->production_vertex()->position().y();
   double z_s = part->production_vertex()->position().z();
  
   EvtVector4R treeStart(ct_s,x_s,y_s,z_s);
   // Add new decay tree to hepMC, converting back from GeV to MeV.
   addEvtGenDecayTree(hepMC, part, evtPart, treeStart, 1000.);
   if(evtPart->getNDaug() !=0) part->set_status(2);
   evtPart->deleteTree();
 }

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty	(	Gaudi::Property< T > &	hndl,
		const SG::VarHandleKeyArrayType &
	)

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKeyArray>

Definition at line 170 of file AthCommonDataStore.h.

   {
     return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                        hndl.value(), 
                                                        hndl.documentation());
  
   }

◆ declareGaudiProperty() [2/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty	(	Gaudi::Property< T > &	hndl,
		const SG::VarHandleKeyType &
	)

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

   {
     return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                        hndl.value(), 
                                                        hndl.documentation());
  
   }

◆ declareGaudiProperty() [3/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty	(	Gaudi::Property< T > &	hndl,
		const SG::VarHandleType &
	)

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleBase>

Definition at line 184 of file AthCommonDataStore.h.

   {
     return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                        hndl.value(), 
                                                        hndl.documentation());
   }

◆ declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty	(	Gaudi::Property< T > &	t,
		const SG::NotHandleType &
	)

inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

   {
     return PBASE::declareProperty(t);
   }

◆ declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty	(	const std::string &	name,
		SG::VarHandleBase &	hndl,
		const std::string &	doc,
		const SG::VarHandleType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

   {
     this->declare(hndl.vhKey());
     hndl.vhKey().setOwner(this);
  
     return PBASE::declareProperty(name,hndl,doc);
   }

◆ declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty	(	const std::string &	name,
		SG::VarHandleKey &	hndl,
		const std::string &	doc,
		const SG::VarHandleKeyType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

   {
     this->declare(hndl);
     hndl.setOwner(this);
  
     return PBASE::declareProperty(name,hndl,doc);
   }

◆ declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty	(	const std::string &	name,
		SG::VarHandleKeyArray &	hndArr,
		const std::string &	doc,
		const SG::VarHandleKeyArrayType &
	)

inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

   {
  
     // std::ostringstream ost;
     // ost << Algorithm::name() << " VHKA declareProp: " << name 
     //     << " size: " << hndArr.keys().size() 
     //     << " mode: " << hndArr.mode() 
     //     << "  vhka size: " << m_vhka.size()
     //     << "\n";
     // debug() << ost.str() << endmsg;
  
     hndArr.setOwner(this);
     m_vhka.push_back(&hndArr);
  
     Gaudi::Details::PropertyBase* p =  PBASE::declareProperty(name, hndArr, doc);
     if (p != 0) {
       p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
     } else {
       ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray " 
                     << name);
     }
  
     return p;
  
   }

◆ declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty	(	const std::string &	name,
		T &	property,
		const std::string &	doc,
		const SG::NotHandleType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

   {
     return PBASE::declareProperty(name, property, doc);
   }

◆ declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty	(	const std::string &	name,
		T &	property,
		const std::string &	doc = `"none"`
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

   {
     typedef typename SG::HandleClassifier<T>::type htype;
     return declareProperty (name, property, doc, htype());
   }

◆ declareProperty() [6/6]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( Gaudi::Property< T > & t )

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                  {
     typedef typename SG::HandleClassifier<T>::type htype;
     return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
   }

◆ detStore()

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< Algorithm > >::detStore ( ) const

inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

95 { return m_detStore; }

◆ event_const()

const HepMC::GenEvent* GenBase::event_const ( ) const

inlineinherited

Access the current signal event (const)

Definition at line 83 of file GenBase.h.

                                            {
     if (events_const()->empty())
       ATH_MSG_ERROR("Const McEventCollection is empty during first event access");
     return *(events_const()->begin());
   }

◆ events_const() [1/2]

const McEventCollection* GenBase::events_const ( ) const

inlineinherited

Access the current event's McEventCollection (const)

Definition at line 96 of file GenBase.h.

                                                 {
     return events_const( getContext() );
   }

◆ events_const() [2/2]

const McEventCollection* GenBase::events_const ( const EventContext & ctx ) const

inlineinherited

Definition at line 99 of file GenBase.h.

                                                                          {
     SG::ReadHandle<McEventCollection> ret = SG::makeHandle(m_mcevents_const, ctx);
     if (!ret.isValid())
       ATH_MSG_ERROR("No McEventCollection found in StoreGate with key " << m_mcevents_const.key());
     return ret.cptr();
   }

◆ evtStore() [1/2]

ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< Algorithm > >::evtStore ( )

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

85 { return m_evtStore; }

◆ evtStore() [2/2]

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< Algorithm > >::evtStore ( ) const

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ execute()

StatusCode EvtInclusiveDecay::execute ( )

virtual

Reimplemented from GenBase.

Definition at line 220 of file EvtInclusiveDecay.cxx.

                                       {
   ATH_MSG_DEBUG("EvtInclusiveDecay executing");
  
   const EventContext& ctx = Gaudi::Hive::currentContext();
   reseedRandomEngine(m_randomStreamName, ctx);
  
   std::string   key = m_inputKeyName;
   // retrieve event from Transient Store (Storegate)
  
   // Load HepMC info
   // FIXME should be using Read/WriteHandles here
   const McEventCollection* oldmcEvtColl{};
   if(m_readExisting) {
     CHECK(evtStore()->retrieve(oldmcEvtColl, key));
     // Fill the new McEventCollection with a copy of the initial HepMC::GenEvent
     m_mcEvtColl = new McEventCollection(*oldmcEvtColl);
   }
   else {CHECK(evtStore()->retrieve(m_mcEvtColl, key));}
  
   if(m_readExisting) {
     if(m_outputKeyName!=key) {
      CHECK(evtStore()->record( m_mcEvtColl,m_outputKeyName));
     }
   }
  
   McEventCollection::iterator mcItr;
   for( mcItr = m_mcEvtColl->begin(); mcItr != m_mcEvtColl->end(); ++mcItr )   {
     HepMC::GenEvent* hepMC = *mcItr;
  
     // Search HepMC record for particles to be decayed by EvtGen
     // NOTE: In order to ensure repeatability, we use a std::set of barcodes to obtain
     //       an ordered list of particles to be decayed by EvtGen.
     std::set<HepMC::GenVertexPtr> visited;
 #ifdef HEPMC3
     std::set<HepMC::GenParticlePtr> toBeDecayed;
     for (auto p: hepMC->particles()) {
       if ( (!p->production_vertex()) ||
            (p->production_vertex()->particles_in().size() == 0) ) {
         StatusCode sc = traverseDecayTree(p,false,visited,toBeDecayed);
         if (sc.isFailure())
           return StatusCode::FAILURE;
       }
     }
     // Print HepMC in tree format if desired (before doing anything)
     if (m_printHepMCBeforeEvtGen) {
       msg(MSG::INFO) << "Printing HepMC record at " << hepMC << " BEFORE running EvtGen:" << endmsg;
       if (m_printHepMCHighLightTopLevelDecays)
         printHepMC(hepMC,&toBeDecayed);
       else
         printHepMC(hepMC);
     }
 #else
     std::set<int> toBeDecayed;
     for (HepMC::GenEvent::particle_iterator itp = hepMC->particles_begin(); itp != hepMC->particles_end(); ++itp) {
       HepMC::GenParticle* p = *itp;
       if ( (!p->production_vertex()) ||
            (p->production_vertex()->particles_in_size() == 0) ) {
         StatusCode sc = traverseDecayTree(p,false,visited,toBeDecayed);
         if (sc.isFailure())
           return StatusCode::FAILURE;
       }
     }
     // Print HepMC in tree format if desired (before doing anything)
     if (m_printHepMCBeforeEvtGen) {
       msg(MSG::INFO) << "Printing HepMC record at " << hepMC << " BEFORE running EvtGen:" << endmsg;
       if (m_printHepMCHighLightTopLevelDecays)
         printHepMC(hepMC,&toBeDecayed);
       else
         printHepMC(hepMC);
     }
 #endif
  
     // Decay selected particles
     bool eventPassesCuts(false);
     int loopCounter(0);
     while( !eventPassesCuts && loopCounter < m_maxNRepeatedDecays ) {
 #ifdef HEPMC3
       for (auto p: toBeDecayed) {
         if (p==0) {
           msg(MSG::ERROR ) << "Overlapping decay tree for particle" << p <<endmsg;
           return StatusCode::FAILURE;
         }
         decayParticle(hepMC,p);
         HepMC::fillBarcodesAttribute(hepMC);
       }
 #else
       for (std::set<int>::iterator itb = toBeDecayed.begin(); itb!=toBeDecayed.end(); ++itb) {
         auto p = hepMC->barcode_to_particle(*itb);
         if (p==0) {
           msg(MSG::ERROR ) << "Overlapping decay tree encountered for barcode " << *itb << endmsg;
           return StatusCode::FAILURE;
         }
         decayParticle(hepMC,p);
         HepMC::fillBarcodesAttribute(hepMC);
       }
 #endif
  
       if(m_applyUserSelection)
         eventPassesCuts = passesUserSelection(hepMC);
       else
         eventPassesCuts = true;
  
       m_nRepeatedDecays++;
       loopCounter++;
     }
  
     // Store the number of decay attempts in event weights std::map, only if repeated decays enabled
 #ifdef HEPMC3
     if(m_maxNRepeatedDecays > 1)
       hepMC->weight("nEvtGenDecayAttempts") = loopCounter;
 #else
     if(m_maxNRepeatedDecays > 1)
       hepMC->weights()["nEvtGenDecayAttempts"] = loopCounter;
 #endif
     // Print HepMC in tree format if desired (after finishing all EvtGen decays)
     if (m_printHepMCAfterEvtGen) {
       msg(MSG::INFO)  << "Printing HepMC record at " << hepMC << " AFTER running EvtGen:" << endmsg;
       if (m_printHepMCHighLightTopLevelDecays)
         printHepMC(hepMC,&toBeDecayed);
       else
         printHepMC(hepMC);
     }
   }
  
   if(m_readExisting && m_outputKeyName==key) {
     CHECK(evtStore()->overwrite(m_mcEvtColl, m_outputKeyName));
   }
  
   return StatusCode::SUCCESS;
 }

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase & ExtraDeps )

protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

◆ extraOutputDeps()

const DataObjIDColl & AthAlgorithm::extraOutputDeps ( ) const

overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 50 of file AthAlgorithm.cxx.

 {
   // If we didn't find any symlinks to add, just return the collection
   // from the base class.  Otherwise, return the extended collection.
   if (!m_extendedExtraObjects.empty()) {
     return m_extendedExtraObjects;
   }
   return Algorithm::extraOutputDeps();
 }

◆ finalize()

StatusCode EvtInclusiveDecay::finalize ( )

Definition at line 352 of file EvtInclusiveDecay.cxx.

                                        {
  
   if (m_checkDecayChannels) {
     ATH_MSG_INFO("The following particles were checked and didn't have any decay channels:");
     if (msgLvl(MSG::INFO)) {
       std::cout << std::endl;
       std::cout << " Particle code    Name from HepPDT        # Occurences" << std::endl;
       std::cout << "------------------------------------------------------"  << std::endl;
       for (std::map<int,long>::iterator p = m_noDecayChannels.begin(); p!=m_noDecayChannels.end(); ++p) {
         int id = p->first;
         int count = p->second;
         std::cout << std::setw(14) << id
                   << std::setw(20) << HepPID::particleName(id)
                   << std::setw(20) << count
                   << std::endl;
       }
       std::cout << std::endl;
     }
   }
   ATH_MSG_INFO("Total number of repeated decays: " << m_nRepeatedDecays);
   ATH_MSG_INFO("EvtInclusiveDecay finalized");
   return StatusCode::SUCCESS;
 }

◆ getRandomEngine()

CLHEP::HepRandomEngine * EvtInclusiveDecay::getRandomEngine	(	const std::string &	streamName,
		unsigned long int	randomSeedOffset,
		const EventContext &	ctx
	)		const

private

Definition at line 194 of file EvtInclusiveDecay.cxx.

 {
   ATHRNG::RNGWrapper* rngWrapper = m_rndmSvc->getEngine(this, streamName);
   rngWrapper->setSeed( streamName, ctx.slot(), randomSeedOffset, ctx.eventID().run_number() );
   return rngWrapper->getEngine(ctx);
 }

◆ getRandomEngineDuringInitialize()

CLHEP::HepRandomEngine * EvtInclusiveDecay::getRandomEngineDuringInitialize	(	const std::string &	streamName,
		unsigned long int	randomSeedOffset,
		unsigned int	conditionsRun = `1`,
		unsigned int	lbn = `1`
	)		const

private

Definition at line 203 of file EvtInclusiveDecay.cxx.

 {
   const size_t slot=0;
   EventContext ctx;
   ctx.setSlot( slot );
   ctx.setEventID (EventIDBase (conditionsRun,
                EventIDBase::UNDEFEVT,  // event
                EventIDBase::UNDEFNUM,  // timestamp
                EventIDBase::UNDEFNUM,  // timestamp ns
                lbn));
   Atlas::setExtendedEventContext(ctx,
                                  Atlas::ExtendedEventContext( evtStore()->hiveProxyDict(),
                                                               conditionsRun) );
   return getRandomEngine(streamName, randomSeedOffset, ctx);
 }

◆ GeVToMeV()

void GenBase::GeVToMeV ( HepMC::GenEvent * evt )

protectedinherited

Scale event energies/momenta by x 1000.

Todo:: Add HepMC units awareness and do it differently when HepMC provides this functionality directly (and reference-based FourVector accessors)

Definition at line 61 of file GenBase.cxx.

                                          {
   for (HepMC::GenEvent::particle_iterator p = evt->particles_begin(); p != evt->particles_end(); ++p) {
     const HepMC::FourVector fv((*p)->momentum().px() * 1000,
                                (*p)->momentum().py() * 1000,
                                (*p)->momentum().pz() * 1000,
                                (*p)->momentum().e()  * 1000);
     (*p)->set_momentum(fv);
     (*p)->set_generated_mass(1000 * (*p)->generated_mass());
   }
 }

◆ initialize()

StatusCode EvtInclusiveDecay::initialize ( )

virtual

Reimplemented from GenBase.

Definition at line 119 of file EvtInclusiveDecay.cxx.

                                          {
  
   ATH_CHECK( GenBase::initialize() );
   // Get the random number service
   CHECK(m_rndmSvc.retrieve());
  
   msg(MSG::INFO) << "EvtInclusiveDecay initialize" << endmsg;
   msg(MSG::INFO) << "Particle properties definition file = " << m_pdtFile << endmsg;
   msg(MSG::INFO) << "Main decay file                     = " << m_decayFile << endmsg;
   msg(MSG::INFO) << "User decay file                     = " << m_userDecayFile << endmsg;
   msg(MSG::INFO) << "Max number of repeated decays       = " << m_maxNRepeatedDecays << endmsg;
   msg(MSG::INFO) << "EvtInclusiveDecay selection parameters:" << endmsg;
   msg(MSG::INFO) << "* prohibitFinalStateDecay = " << m_prohibitFinalStateDecay << endmsg;
   msg(MSG::INFO) << "* prohibitReDecay         = " << m_prohibitReDecay << endmsg;
   msg(MSG::INFO) << "* prohibitUnDecay         = " << m_prohibitUnDecay << endmsg;
   msg(MSG::INFO) << "* prohibitRemoveSelfDecay = " << m_prohibitRemoveSelfDecay << endmsg;
   msg(MSG::INFO) << "* allowAllKnownDecays     = " << m_allowAllKnownDecays << endmsg;
   msg(MSG::INFO) << "* allowDefaultBDecays     = " << m_allowDefaultBDecays << endmsg;
   msg(MSG::INFO) << "User selection parameters:" << endmsg;
   msg(MSG::INFO) << "* applyUserSelection             = " << m_applyUserSelection << endmsg;
   msg(MSG::INFO) << "* userSelRequireOppositeSignedMu = " << m_userSelRequireOppositeSignedMu << endmsg;
   msg(MSG::INFO) << "* userSelMu1MinPt                = " << m_userSelMu1MinPt << endmsg;
   msg(MSG::INFO) << "* userSelMu2MinPt                = " << m_userSelMu2MinPt << endmsg;
   msg(MSG::INFO) << "* userSelMu1MaxEta               = " << m_userSelMu1MaxEta << endmsg;
   msg(MSG::INFO) << "* userSelMu2MaxEta               = " << m_userSelMu2MaxEta << endmsg;
   msg(MSG::INFO) << "* userSelMinDimuMass             = " << m_userSelMinDimuMass << endmsg;
   msg(MSG::INFO) << "* userSelMaxDimuMass             = " << m_userSelMaxDimuMass << endmsg;
  
   // Initialize and print blackList
   m_blackListSet.insert(m_blackList.begin(),m_blackList.end());
   msg(MSG::INFO) << "* blackList;               = ";
   for (std::set<int>::iterator i = m_blackListSet.begin(); i!=m_blackListSet.end(); ++i)
      msg(MSG::INFO) << (*i) << " ";
   msg(MSG::INFO)<< endmsg;
  
   // Initialize and print whiteList
   m_whiteListSet.insert(m_whiteList.begin(),m_whiteList.end());
   msg(MSG::INFO) << "* whiteList               = ";
   for (std::set<int>::iterator i = m_whiteListSet.begin(); i!=m_whiteListSet.end(); ++i)
     msg(MSG::INFO) << (*i) << " ";
   msg(MSG::INFO) << endmsg;
  
   CLHEP::HepRandomEngine* rndmEngine = getRandomEngineDuringInitialize(m_randomStreamName, m_randomSeed, m_dsid);
   // Obtain random number generator for EvtGen
   m_evtAtRndmGen = new EvtInclusiveAtRndmGen(rndmEngine);
  
   // Create an instance of EvtGen and read particle properties and decay files
   EvtExternalGenList genList(true,xmlpath(),"gamma");
   EvtAbsRadCorr* radCorrEngine = genList.getPhotosModel();
   std::list<EvtDecayBase*> extraModels = genList.getListOfModels();
  
    // Create the EvtGen generator object
   // EvtGen myGenerator("decayFile.dec", "evt.pdl", randomEnginePointer,
   //                    radCorrEngine, &extraModels);
  
  
   m_myEvtGen = new EvtGen( m_decayFile.c_str(), m_pdtFile.c_str(), m_evtAtRndmGen, radCorrEngine, &extraModels);
   if(!m_userDecayFile.empty())
     m_myEvtGen->readUDecay(m_userDecayFile.c_str());
  
   m_nRepeatedDecays = 0;
  
   return StatusCode::SUCCESS;
 }

◆ inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< Algorithm > >::inputHandles ( ) const

overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ invMass()

double EvtInclusiveDecay::invMass	(	HepMC::ConstGenParticlePtr	p1,
		HepMC::ConstGenParticlePtr	p2
	)

private

Definition at line 699 of file EvtInclusiveDecay.cxx.

                                                                                           {
   double p1Px = p1->momentum().px();
   double p1Py = p1->momentum().py();
   double p1Pz = p1->momentum().pz();
   double p1E  = p1->momentum().e();
   double p2Px = p2->momentum().px();
   double p2Py = p2->momentum().py();
   double p2Pz = p2->momentum().pz();
   double p2E  = p2->momentum().e();
   double dimuE = p2E + p1E;
   double dimuPx = p2Px + p1Px;
   double dimuPy = p2Py + p1Py;
   double dimuPz = p2Pz + p1Pz;
   double invMass = std::sqrt(dimuE*dimuE - dimuPx*dimuPx - dimuPy*dimuPy - dimuPz*dimuPz);
  
   return invMass;
 }

◆ isDefaultB()

bool EvtInclusiveDecay::isDefaultB ( const int pId ) const

private

Definition at line 647 of file EvtInclusiveDecay.cxx.

                                                       {
   int id = std::abs(pId);
   if ( id == 511   ||
        id == 521   ||
        id == 531   ||
        id == 541   ||
        id == 5122  ||
        id == 5132  ||
        id == 5232  ||
        id == 5112  ||
        id == 5212  ||
        id == 5222 )
     return true;
   else
     return false;
 }

◆ isToBeDecayed()

bool EvtInclusiveDecay::isToBeDecayed	(	HepMC::ConstGenParticlePtr	p,
		bool	doCrossChecks
	)

private

Definition at line 570 of file EvtInclusiveDecay.cxx.

                                                                                     {
   int id = p->pdg_id();
   int nDaughters = 0;
   auto v = p->end_vertex();
   if (v) nDaughters = v->particles_out_size();
  
   // Ignore documentation lines
   if (p->status() == 3) return false;
   // And any particles that aren't stable or decayed
   if(!m_isfHerwig && !MC::isPhysical(p)) return false;
  
   // Particularly for Herwig, try to ignore particles that really should
   // be flagged as documentation lines
   double m2 = p->momentum().m2();
   if (m2 < -1.0E-3) {
     ATH_MSG_DEBUG("Ignoring particle " << pdgName(p) << " with m^2 = " << m2);
     return false;
   }
  
   // Check whether EvtGen has any decay channels defined for this particle
   EvtId evtId = EvtPDL::evtIdFromStdHep(id);
   // std::cout << "EVTID: " << evtId.getId() << " alias " << evtId.getAlias() << std::endl;
   int nModes = 0;
   if (evtId.getId()>=0)
     //    nModes = EvtDecayTable::getNMode(evtId.getAlias());
     nModes = EvtDecayTable::getInstance()->getNMode(evtId.getAlias());
   if (doCrossChecks) {
     ATH_MSG_VERBOSE("Checking particle " << pdgName(p)
         << " (status = " << p->status()
                     <<") -- " << nModes << " decay modes found");
     if (m_checkDecayChannels && nModes==0) {
       std::map<int,long>::iterator pos = m_noDecayChannels.find(id);
       if (pos != m_noDecayChannels.end())
         (pos->second)++;
       else
         m_noDecayChannels[id] = 1;
     }
   }
  
   // Check prohibit* settings
   if (m_prohibitFinalStateDecay && MC::isStable(p)) return false;
   if (m_prohibitReDecay && nDaughters>0) return false;
   if (m_prohibitUnDecay && nModes==0) return false;
   if (m_prohibitRemoveSelfDecay && nDaughters>0) {
     // For now, check only children - this should be sufficient and checking all
     // descendants would be very expensive.
 #ifdef HEPMC3
     for (auto  itd: v->particles_out()) {
       if (std::abs(itd->pdg_id()) == std::abs(id)) return false;
     }
 #else
     for (HepMC::GenVertex::particle_iterator itd = v->particles_begin(HepMC::children);
                                              itd != v->particles_end(HepMC::children);
                                              ++itd) {
       if (std::abs((*itd)->pdg_id()) == std::abs(id)) return false;
     }
 #endif
   }
  
   // Check blackList
   if (m_blackListSet.count(std::abs(id))>0) return false;
  
   // Check allow* settings
   if (m_allowAllKnownDecays && nModes>0) return true;
   if (m_allowDefaultBDecays && isDefaultB(id)) return true;
  
   // Check whiteList
   if (m_whiteListSet.count(std::abs(id))>0) return true;
  
   return false;   // Default is NOT to decay through EvtGen
 }

◆ MeVToGeV()

void GenBase::MeVToGeV ( HepMC::GenEvent * evt )

protectedinherited

Scale event energies/momenta by x 1/1000.

Definition at line 71 of file GenBase.cxx.

                                          {
   for (HepMC::GenEvent::particle_iterator p = evt->particles_begin(); p != evt->particles_end(); ++p) {
     const HepMC::FourVector fv((*p)->momentum().px() / 1000,
                                (*p)->momentum().py() / 1000,
                                (*p)->momentum().pz() / 1000,
                                (*p)->momentum().e()  / 1000);
     (*p)->set_momentum(fv);
     (*p)->set_generated_mass((*p)->generated_mass() / 1000);
   }
 }

◆ mmTocm()

void GenBase::mmTocm ( HepMC::GenEvent * evt )

protectedinherited

Scale event lengths by x 1/10.

Definition at line 90 of file GenBase.cxx.

                                        {
   for (HepMC::GenEvent::vertex_iterator vtx = evt->vertices_begin(); vtx != evt->vertices_end(); ++vtx) {
     const HepMC::FourVector fv((*vtx)->position().x() / 10,
                                (*vtx)->position().y() / 10,
                                (*vtx)->position().z() / 10,
                                (*vtx)->position().t() / 10);
     (*vtx)->set_position(fv);
   }
 }

◆ msg() [1/2]

MsgStream& AthCommonMsg< Algorithm >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                 {
     return this->msgStream();
   }

◆ msg() [2/2]

MsgStream& AthCommonMsg< Algorithm >::msg ( const MSG::Level lvl ) const

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                   {
     return this->msgStream(lvl);
   }

◆ msgLvl()

bool AthCommonMsg< Algorithm >::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                                {
     return this->msgLevel(lvl);
   }

◆ outputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< Algorithm > >::outputHandles ( ) const

overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ particleData()

const HepPDT::ParticleData* GenBase::particleData ( int pid ) const

inlineinherited

Access an element in the particle data table.

Definition at line 126 of file GenBase.h.

                                                         {
     return pdt().particle(HepPDT::ParticleID(std::abs(pid)));
   }

◆ particleTable()

const HepPDT::ParticleDataTable& GenBase::particleTable ( ) const

inlineinherited

Get a particle data table.

Definition at line 118 of file GenBase.h.

                                                        {
     return *(m_ppSvc->PDT());
   }

◆ partPropSvc()

const ServiceHandle<IPartPropSvc> GenBase::partPropSvc ( ) const

inlineinherited

Access the particle property service.

Definition at line 113 of file GenBase.h.

                                                         {
     return m_ppSvc;
   }

◆ passesUserSelection()

bool EvtInclusiveDecay::passesUserSelection ( HepMC::GenEvent * hepMC )

private

Definition at line 669 of file EvtInclusiveDecay.cxx.

                                                                 {
   bool passed(false);
   std::vector<HepMC::GenParticlePtr> *muons = new std::vector<HepMC::GenParticlePtr>;
  
   for ( const auto& p:  *hepMC) {
     if( std::abs(p->pdg_id()) == 13 )
       muons->push_back(p);
   }
  
   for (auto muItr1 = muons->begin(); muItr1 != muons->end(); ++muItr1) {
     for (auto muItr2 = muItr1+1; muItr2 != muons->end(); ++muItr2) {
       if( m_userSelRequireOppositeSignedMu && (*muItr1)->pdg_id() * (*muItr2)->pdg_id() > 0)
         continue;
       if( !( (*muItr1)->momentum().perp() > m_userSelMu1MinPt && std::abs((*muItr1)->momentum().pseudoRapidity()) < m_userSelMu1MaxEta &&
              (*muItr2)->momentum().perp() > m_userSelMu2MinPt && std::abs((*muItr2)->momentum().pseudoRapidity()) < m_userSelMu2MaxEta ) &&
           !( (*muItr2)->momentum().perp() > m_userSelMu1MinPt && std::abs((*muItr2)->momentum().pseudoRapidity()) < m_userSelMu1MaxEta &&
              (*muItr1)->momentum().perp() > m_userSelMu2MinPt && std::abs((*muItr1)->momentum().pseudoRapidity()) < m_userSelMu2MaxEta ) )
         continue;
       double dimuMass = invMass((*muItr1),(*muItr2));
       if( !( dimuMass > m_userSelMinDimuMass && (dimuMass < m_userSelMaxDimuMass || m_userSelMaxDimuMass < 0.) ) )
         continue;
       passed = true;
     }
   }
  
   delete muons;
  
   return passed;
 }

◆ pdgName()

std::string EvtInclusiveDecay::pdgName	(	HepMC::ConstGenParticlePtr	p,
		bool	statusHighlighting = `false`,
		std::set< int > *	barcodeList = `nullptr`
	)

private

Definition at line 872 of file EvtInclusiveDecay.cxx.

                                                                                                                   {
   std::ostringstream buf;
   if (statusHighlighting) {
     if ( ((barcodeList!=0) && (barcodeList->find(HepMC::barcode(p)) != barcodeList->end())) ||
          ((barcodeList==0) && isToBeDecayed(p,false)) )
       buf << "\033[7m";   // reverse
     if (!MC::isStable(p)) {
       if (MC::isDecayed(p))
         buf << "\033[33m";   // yellow
       else
         buf << "\033[31m";   // red
     }
   }
   buf << p->pdg_id();
   buf << "/" << HepPID::particleName(p->pdg_id());
   if (statusHighlighting) {
     buf << "\033[0m";   // revert color attributes
   }
   return buf.str();
 }

◆ pdt()

const HepPDT::ParticleDataTable& GenBase::pdt ( ) const

inlineinherited

Shorter alias to get a particle data table.

Definition at line 123 of file GenBase.h.

123 { return particleTable(); }

◆ printHepMC()

void EvtInclusiveDecay::printHepMC	(	HepMC::GenEvent *	hepMC,
		std::set< int > *	barcodeList = `nullptr`
	)

private

Definition at line 741 of file EvtInclusiveDecay.cxx.

                                                                                  {
   std::set<HepMC::GenVertexPtr> visited;
   unsigned int nParticlesFound = 0;
   unsigned int nTreesFound = 0;
   for (HepMC::GenEvent::particle_iterator itp = hepMC->particles_begin(); itp != hepMC->particles_end(); ++itp) {
     HepMC::GenParticle* p = *itp;
     if ( (!p->production_vertex()) ||
          (p->production_vertex()->particles_in_size() == 0) ) {
       nTreesFound++;
       std::cout << "\n    Found new partial decay tree:\n" << std::endl;
       unsigned int nParticlesVisited = printTree(p,visited,1,barcodeList);
       std::cout << "\n    " << nParticlesVisited << " particles in this subtree" << std::endl;
       nParticlesFound += nParticlesVisited;
     }
   }
   std::cout << "\n    Total of " << nParticlesFound << " particles found in "
             << nTreesFound << " decay subtrees in HepMC event record\n" << std::endl;
 }

◆ printTree()

unsigned int EvtInclusiveDecay::printTree	(	HepMC::GenParticlePtr	p,
		std::set< HepMC::GenVertexPtr > &	visited,
		int	level,
		std::set< int > *	barcodeList = `0`
	)

private

Definition at line 791 of file EvtInclusiveDecay.cxx.

                                                                                                           {
   unsigned int nParticlesVisited = 1;
   for (int i=0; i<level; i++) std::cout << "    ";
   std::cout << pdgName(p,m_printHepMCHighlighted,barcodeList);
   auto v = p->end_vertex();
 #ifdef HEPMC3
   if (v) {
     if (v->particles_in().size() > 1)
       std::cout << " [interaction: " << v->particles_in().size() << " particles, vertex " << v << "]    -->   ";
     else
       std::cout << "   -->   ";
     if (visited.insert(v).second) {
       for (auto itp: v->particles_out()) {
         std::cout << pdgName(itp,m_printHepMCHighlighted,barcodeList) << "   ";
       }
       std::cout << std::endl;
       for (auto itp: v->particles_out()) {
         if (itp->end_vertex())
           nParticlesVisited += printTree(itp, visited, level+1, barcodeList);
         else
           nParticlesVisited++;
       }
     } else
       std:: cout << "see above" << std::endl;
   } else
     std::cout << "   no decay vertex\n" << std::endl;
 #else
   if (v) {
     if (v->particles_in_size() > 1)
       std::cout << " [interaction: " << v->particles_in_size() << " particles, vertex " << v << "]    -->   ";
     else
       std::cout << "   -->   ";
     if (visited.insert(v).second) {
       for (HepMC::GenVertex::particle_iterator itp = v->particles_begin(HepMC::children);
                                                itp != v->particles_end(HepMC::children);
                                                ++itp) {
         std::cout << pdgName(*itp,m_printHepMCHighlighted,barcodeList) << "   ";
       }
       std::cout << std::endl;
       for (HepMC::GenVertex::particle_iterator itp = v->particles_begin(HepMC::children);
                                                itp != v->particles_end(HepMC::children);
                                                ++itp) {
         if ((*itp)->end_vertex())
           nParticlesVisited += printTree(*itp, visited, level+1, barcodeList);
         else
           nParticlesVisited++;
       }
     } else
       std:: cout << "see above" << std::endl;
   } else
     std::cout << "   no decay vertex\n" << std::endl;
 #endif
   return nParticlesVisited;
 }

◆ removeDecayTree()

void EvtInclusiveDecay::removeDecayTree	(	HepMC::GenEvent *	hepMC,
		HepMC::GenParticlePtr	p
	)

private

Definition at line 449 of file EvtInclusiveDecay.cxx.

                                                                                    {
   auto v = p->end_vertex();
   if (v) {
 #ifdef HEPMC3
     //This is recursive in HepMC3. But explicit deletion is allowed as well.
     hepMC->remove_vertex(v);
     p->set_status(1);   // For now, flag particle as undecayed (stable)
     ATH_MSG_DEBUG("Removed existing " << pdgName(p) << " " << p  );
 #else
     std::set<int> vtxBarCodesToDelete;
     vtxBarCodesToDelete.insert(v->barcode());
     for (HepMC::GenVertex::vertex_iterator itv = v->vertices_begin(HepMC::descendants);
                                            itv != v->vertices_end(HepMC::descendants);
                                            ++itv)
       vtxBarCodesToDelete.insert((*itv)->barcode());
     for (std::set<int>::iterator itb = vtxBarCodesToDelete.begin(); itb != vtxBarCodesToDelete.end(); ++itb) {
       auto vdel = hepMC->barcode_to_vertex(*itb);
       hepMC->remove_vertex(vdel);
       delete vdel;
     }
     p->set_status(1);   // For now, flag particle as undecayed (stable)
     ATH_MSG_DEBUG("Removed existing " << pdgName(p) << " (barcode " << p->barcode() << ")"
                   << " decay tree with " << vtxBarCodesToDelete.size() << " vertices");
 #endif
   }
 }

◆ renounce()

std::enable_if_t<std::is_void_v<std::result_of_t<decltype(&T::renounce)(T)> > && !std::is_base_of_v<SG::VarHandleKeyArray, T> && std::is_base_of_v<Gaudi::DataHandle, T>, void> AthCommonDataStore< AthCommonMsg< Algorithm > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

   {
     h.renounce();
     PBASE::renounce (h);
   }

◆ renounceArray()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::renounceArray ( SG::VarHandleKeyArray & handlesArray )

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                           {
     handlesArray.renounce();
   }

◆ reseedRandomEngine()

void EvtInclusiveDecay::reseedRandomEngine	(	const std::string &	streamName,
		const EventContext &	ctx
	)

private

Definition at line 185 of file EvtInclusiveDecay.cxx.

 {
   long seeds[7];
   ATHRNG::calculateSeedsMC21(seeds, streamName,  ctx.eventID().event_number(), m_dsid, m_randomSeed);
   m_evtAtRndmGen->getEngine()->setSeeds(seeds, 0); // NOT THREAD-SAFE
 }

◆ sysInitialize()

StatusCode AthAlgorithm::sysInitialize ( )

overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc

Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc

Reimplemented from AthCommonDataStore< AthCommonMsg< Algorithm > >.

Reimplemented in AthAnalysisAlgorithm, AthFilterAlgorithm, PyAthena::Alg, and AthHistogramAlgorithm.

Definition at line 66 of file AthAlgorithm.cxx.

                                        {
   StatusCode sc=AthCommonDataStore<AthCommonMsg<Algorithm>>::sysInitialize();
  
   if (sc.isFailure()) {
     return sc;
   }
   ServiceHandle<ICondSvc> cs("CondSvc",name());
   for (auto h : outputHandles()) {
     if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
       // do this inside the loop so we don't create the CondSvc until needed
       if ( cs.retrieve().isFailure() ) {
         ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
         return StatusCode::SUCCESS;
       }
       if (cs->regHandle(this,*h).isFailure()) {
         sc = StatusCode::FAILURE;
         ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
                       << " with CondSvc");
       }
     }
   }
   return sc;
 }

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Algorithm > >::sysStart ( )

overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

◆ traverseDecayTree()

StatusCode EvtInclusiveDecay::traverseDecayTree	(	HepMC::GenParticlePtr	p,
		bool	isToBeRemoved,
		std::set< HepMC::GenVertexPtr > &	visited,
		std::set< int > &	toBeDecayed
	)

private

Definition at line 391 of file EvtInclusiveDecay.cxx.

                                                                           {
 #endif
   ATH_MSG_VERBOSE("Inspecting: " << pdgName(p) << " " << p);
   if (!isToBeRemoved) {
     if (isToBeDecayed(p,true)) {
 #ifdef HEPMC3
       toBeDecayed.insert(p);
 #else
       toBeDecayed.insert(HepMC::barcode(p));
 #endif
       isToBeRemoved = true;
       ATH_MSG_VERBOSE("Selected particle for decay: " << pdgName(p) << " " << p );
  
       // In principle we could stop the recursion here. However, to prevent
       // pathological cases in certain decay trees (in particular from Herwig),
       // we continue in order to mark all descendants of this particle
       // as visited. Thus none of these descendants can be flagged for further
       // decay, even if it has several mothers.
     }
   }
   auto v = p->end_vertex();
 #ifdef HEPMC3
   if (v) {
     if (visited.insert(v).second) {
       if ( isToBeRemoved && (v->particles_in().size()>1) && m_checkDecayTree ) {
         ATH_MSG_WARNING("Found particle to be decayed with vertex with >1 incoming mother particles in decay tree");
         ATH_MSG_WARNING( ([&p, &v](){  std::stringstream ss;   HepMC::Print::line(ss,p); HepMC::Print::line(ss,v);  return ss.str();})());
       }
       for (auto itp: v->particles_out()) {
         ATH_CHECK(traverseDecayTree(itp,isToBeRemoved,visited,toBeDecayed) );
       }
     }
   }
 #else
   if (v) {
     if (visited.insert(v).second) {
       if ( isToBeRemoved && (v->particles_in_size()>1) && m_checkDecayTree ) {
         // This is normal for Herwig but should not occur for Pythia
         ATH_MSG_WARNING("Found particle to be decayed with vertex with >1 incoming mother particles in decay tree");
         ATH_MSG_WARNING( ([&p, &v](){  std::stringstream ss;   HepMC::Print::line(ss,p); HepMC::Print::line(ss,v);  return ss.str();})());
       }
       for (auto itp = v->particles_begin(HepMC::children); itp != v->particles_end(HepMC::children); ++itp) {
         ATH_CHECK(traverseDecayTree(*itp,isToBeRemoved,visited,toBeDecayed) );
       }
     }
   }
 #endif
   return StatusCode::SUCCESS;
 }

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::updateVHKA ( Gaudi::Details::PropertyBase & )

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                                {
     // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
     //         << "  size: " << m_vhka.size() << endmsg;
     for (auto &a : m_vhka) {
       std::vector<SG::VarHandleKey*> keys = a->keys();
       for (auto k : keys) {
         k->setOwner(this);
       }
     }
   }

◆ xmlpath()

std::string EvtInclusiveDecay::xmlpath ( void )

Definition at line 920 of file EvtInclusiveDecay.cxx.

                                     {
  
   char *cmtpath = getenv("CMTPATH");
   char *cmtconfig = getenv("CMTCONFIG");
  
   std::string foundpath = "";
  
   if(cmtpath != 0 && cmtconfig != 0){
  
     std::vector<std::string> cmtpaths;
     local_split(cmtpaths, cmtpath,':');
  
     std::string installPath = "/InstallArea/" + std::string(cmtconfig) + "/share/Pythia8/xmldoc";
  
     for(std::vector<std::string>::const_iterator path = cmtpaths.begin();
         path != cmtpaths.end() && foundpath == ""; ++path){
       std::string testPath = *path + installPath;
       std::ifstream testFile(testPath.c_str());
       if(testFile.good()) foundpath = testPath;
       testFile.close();
     }
    }
   else {
 // If the CMT environment is missing, try to find the xmldoc directory
 // using PathResolver:
     foundpath = PathResolverFindCalibDirectory( "Pythia8/xmldoc" );
  
   }
  
  
   return foundpath;
 }

Member Data Documentation

◆ m_allowAllKnownDecays

bool EvtInclusiveDecay::m_allowAllKnownDecays

private

Definition at line 143 of file EvtInclusiveDecay.h.

◆ m_allowDefaultBDecays

bool EvtInclusiveDecay::m_allowDefaultBDecays

private

Definition at line 144 of file EvtInclusiveDecay.h.

◆ m_applyUserSelection

bool EvtInclusiveDecay::m_applyUserSelection

private

Definition at line 161 of file EvtInclusiveDecay.h.

◆ m_blackList

std::vector<int> EvtInclusiveDecay::m_blackList

private

Definition at line 140 of file EvtInclusiveDecay.h.

◆ m_blackListSet

std::set<int> EvtInclusiveDecay::m_blackListSet

private

Definition at line 141 of file EvtInclusiveDecay.h.

◆ m_checkDecayChannels

bool EvtInclusiveDecay::m_checkDecayChannels

private

Definition at line 154 of file EvtInclusiveDecay.h.

◆ m_checkDecayTree

bool EvtInclusiveDecay::m_checkDecayTree

private

Definition at line 153 of file EvtInclusiveDecay.h.

◆ m_decayFile

std::string EvtInclusiveDecay::m_decayFile

private

Definition at line 129 of file EvtInclusiveDecay.h.

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_dsid

IntegerProperty EvtInclusiveDecay::m_dsid {this, "Dsid", 999999}

private

Definition at line 116 of file EvtInclusiveDecay.h.

◆ m_evtAtRndmGen

EvtInclusiveAtRndmGen* EvtInclusiveDecay::m_evtAtRndmGen {}

private

Definition at line 124 of file EvtInclusiveDecay.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthAlgorithm::m_extendedExtraObjects

privateinherited

Definition at line 79 of file AthAlgorithm.h.

◆ m_inputKeyName

std::string EvtInclusiveDecay::m_inputKeyName

private

Definition at line 132 of file EvtInclusiveDecay.h.

◆ m_isfHerwig

bool EvtInclusiveDecay::m_isfHerwig

private

Definition at line 169 of file EvtInclusiveDecay.h.

◆ m_maxNRepeatedDecays

int EvtInclusiveDecay::m_maxNRepeatedDecays

private

Definition at line 159 of file EvtInclusiveDecay.h.

◆ m_mcEventKey

std::string GenBase::m_mcEventKey {}

protectedinherited

StoreGate key for the MC event collection (defaults to GEN_EVENT)

Definition at line 137 of file GenBase.h.

◆ m_mcevents_const

SG::ReadHandleKey<McEventCollection> GenBase::m_mcevents_const { this, "McEventKey", "GEN_EVENT", "StoreGate key of the MC event collection" }

privateinherited

Const handle to the MC event collection.

Definition at line 163 of file GenBase.h.

◆ m_mcEvtColl

McEventCollection* EvtInclusiveDecay::m_mcEvtColl {}

private

Definition at line 121 of file EvtInclusiveDecay.h.

◆ m_mkMcEvent

BooleanProperty GenBase::m_mkMcEvent {this, "MakeMcEvent", false, "Create a new MC event collection if it doesn't exist"}

protectedinherited

Flag to determine if a new MC event collection should be made if it doesn't exist.

Definition at line 139 of file GenBase.h.

◆ m_myEvtGen

EvtGen* EvtInclusiveDecay::m_myEvtGen {}

private

Definition at line 125 of file EvtInclusiveDecay.h.

◆ m_noDecayChannels

std::map<int,long> EvtInclusiveDecay::m_noDecayChannels

private

Definition at line 155 of file EvtInclusiveDecay.h.

◆ m_nRepeatedDecays

int EvtInclusiveDecay::m_nRepeatedDecays

private

Definition at line 157 of file EvtInclusiveDecay.h.

◆ m_outputKeyName

std::string EvtInclusiveDecay::m_outputKeyName

private

Definition at line 133 of file EvtInclusiveDecay.h.

◆ m_pdtFile

std::string EvtInclusiveDecay::m_pdtFile

private

Definition at line 128 of file EvtInclusiveDecay.h.

◆ m_ppSvc

ServiceHandle<IPartPropSvc> GenBase::m_ppSvc {this, "PartPropSvc", "PartPropSvc"}

privateinherited

Handle on the particle property service.

Definition at line 160 of file GenBase.h.

◆ m_printHepMCAfterEvtGen

bool EvtInclusiveDecay::m_printHepMCAfterEvtGen

private

Definition at line 149 of file EvtInclusiveDecay.h.

◆ m_printHepMCBeforeEvtGen

bool EvtInclusiveDecay::m_printHepMCBeforeEvtGen

private

Definition at line 148 of file EvtInclusiveDecay.h.

◆ m_printHepMCHighlighted

bool EvtInclusiveDecay::m_printHepMCHighlighted

private

Definition at line 150 of file EvtInclusiveDecay.h.

◆ m_printHepMCHighLightTopLevelDecays

bool EvtInclusiveDecay::m_printHepMCHighLightTopLevelDecays

private

Definition at line 151 of file EvtInclusiveDecay.h.

◆ m_prohibitFinalStateDecay

bool EvtInclusiveDecay::m_prohibitFinalStateDecay

private

Definition at line 136 of file EvtInclusiveDecay.h.

◆ m_prohibitReDecay

bool EvtInclusiveDecay::m_prohibitReDecay

private

Definition at line 137 of file EvtInclusiveDecay.h.

◆ m_prohibitRemoveSelfDecay

bool EvtInclusiveDecay::m_prohibitRemoveSelfDecay

private

Definition at line 139 of file EvtInclusiveDecay.h.

◆ m_prohibitUnDecay

bool EvtInclusiveDecay::m_prohibitUnDecay

private

Definition at line 138 of file EvtInclusiveDecay.h.

◆ m_randomSeed

IntegerProperty EvtInclusiveDecay::m_randomSeed {this, "RandomSeed", 1234567, "Random seed for the built-in random engine"}

private

Seed for random number engine.

Definition at line 119 of file EvtInclusiveDecay.h.

◆ m_randomStreamName

std::string EvtInclusiveDecay::m_randomStreamName

private

Definition at line 131 of file EvtInclusiveDecay.h.

◆ m_readExisting

bool EvtInclusiveDecay::m_readExisting

private

Definition at line 135 of file EvtInclusiveDecay.h.

◆ m_rndmSvc

ServiceHandle<IAthRNGSvc> EvtInclusiveDecay::m_rndmSvc {this, "RndmSvc", "AthRNGSvc"}

private

Definition at line 113 of file EvtInclusiveDecay.h.

◆ m_setVMtransversePol

bool EvtInclusiveDecay::m_setVMtransversePol

private

Definition at line 170 of file EvtInclusiveDecay.h.

◆ m_userDecayFile

std::string EvtInclusiveDecay::m_userDecayFile

private

Definition at line 130 of file EvtInclusiveDecay.h.

◆ m_userSelMaxDimuMass

double EvtInclusiveDecay::m_userSelMaxDimuMass

private

Definition at line 168 of file EvtInclusiveDecay.h.

◆ m_userSelMinDimuMass

double EvtInclusiveDecay::m_userSelMinDimuMass

private

Definition at line 167 of file EvtInclusiveDecay.h.

◆ m_userSelMu1MaxEta

double EvtInclusiveDecay::m_userSelMu1MaxEta

private

Definition at line 165 of file EvtInclusiveDecay.h.

◆ m_userSelMu1MinPt

double EvtInclusiveDecay::m_userSelMu1MinPt

private

Definition at line 163 of file EvtInclusiveDecay.h.

◆ m_userSelMu2MaxEta

double EvtInclusiveDecay::m_userSelMu2MaxEta

private

Definition at line 166 of file EvtInclusiveDecay.h.

◆ m_userSelMu2MinPt

double EvtInclusiveDecay::m_userSelMu2MinPt

private

Definition at line 164 of file EvtInclusiveDecay.h.

◆ m_userSelRequireOppositeSignedMu

bool EvtInclusiveDecay::m_userSelRequireOppositeSignedMu

private

Definition at line 162 of file EvtInclusiveDecay.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Algorithm > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< Algorithm > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

◆ m_whiteList

std::vector<int> EvtInclusiveDecay::m_whiteList

private

Definition at line 145 of file EvtInclusiveDecay.h.

◆ m_whiteListSet

std::set<int> EvtInclusiveDecay::m_whiteListSet

private

Definition at line 146 of file EvtInclusiveDecay.h.

The documentation for this class was generated from the following files:

Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Features for derived classes to use internally

Utility event-mangling functions

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ EvtInclusiveDecay()

◆ ~EvtInclusiveDecay()

Member Function Documentation

◆ addEvtGenDecayTree()

◆ ATLAS_NOT_CONST_THREAD_SAFE() [1/2]

◆ ATLAS_NOT_CONST_THREAD_SAFE() [2/2]

◆ cmTomm()

◆ decayParticle()

◆ declareGaudiProperty() [1/4]

◆ declareGaudiProperty() [2/4]

◆ declareGaudiProperty() [3/4]

◆ declareGaudiProperty() [4/4]

◆ declareProperty() [1/6]

◆ declareProperty() [2/6]

◆ declareProperty() [3/6]

◆ declareProperty() [4/6]

◆ declareProperty() [5/6]

◆ declareProperty() [6/6]

◆ detStore()

◆ event_const()

◆ events_const() [1/2]

◆ events_const() [2/2]

◆ evtStore() [1/2]

◆ evtStore() [2/2]

◆ execute()

◆ extraDeps_update_handler()

◆ extraOutputDeps()

◆ finalize()

◆ getRandomEngine()

◆ getRandomEngineDuringInitialize()

◆ GeVToMeV()

◆ initialize()

◆ inputHandles()

◆ invMass()

◆ isDefaultB()

◆ isToBeDecayed()

◆ MeVToGeV()

◆ mmTocm()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ outputHandles()

◆ particleData()

◆ particleTable()

◆ partPropSvc()

◆ passesUserSelection()

◆ pdgName()

◆ pdt()

◆ printHepMC()

◆ printTree()

◆ removeDecayTree()

◆ renounce()

◆ renounceArray()

◆ reseedRandomEngine()

◆ sysInitialize()

◆ sysStart()

◆ traverseDecayTree()

◆ updateVHKA()

◆ xmlpath()

Member Data Documentation

◆ m_allowAllKnownDecays

◆ m_allowDefaultBDecays

◆ m_applyUserSelection

◆ m_blackList

◆ m_blackListSet

◆ m_checkDecayChannels

◆ m_checkDecayTree

◆ m_decayFile

◆ m_detStore