#include <BPHDecayMomentum.h>

Inheritance diagram for BPHDecayMomentum:

Classes
struct	Component

Public Member Functions
virtual const std::vector< std::string > &	compNames () const

virtual const pat::CompositeCandidate &	composite () const
	get a composite by the simple sum of simple particles More...

virtual const std::vector< BPHRecoConstCandPtr > &	daughComp () const

virtual const std::vector< const reco::Candidate * > &	daughFull () const

virtual const std::vector< const reco::Candidate * > &	daughters () const

virtual const std::vector< std::string > &	daugNames () const

virtual BPHRecoConstCandPtr	getComp (const std::string &name) const

virtual const reco::Candidate *	getDaug (const std::string &name) const

virtual const reco::Candidate *	originalReco (const reco::Candidate *daug) const
	get the original particle from the clone More...

virtual	~BPHDecayMomentum ()

Protected Member Functions
virtual void	addP (const std::string &name, const reco::Candidate *daug, double mass=-1.0)

virtual void	addP (const std::string &name, const BPHRecoConstCandPtr &comp)
	add a previously reconstructed particle giving it a name More...

	BPHDecayMomentum ()

	BPHDecayMomentum (const std::map< std::string, Component > &daugMap)

	BPHDecayMomentum (const std::map< std::string, Component > &daugMap, const std::map< std::string, BPHRecoConstCandPtr > compMap)

const std::vector< Component > &	componentList () const

virtual void	setNotUpdated () const

Private Member Functions
void	clonesList (const std::map< std::string, Component > &daugMap)

virtual void	computeMomentum () const

void	dCompList ()

virtual void	fillDaug (std::vector< const reco::Candidate * > &ad) const

virtual void	sumMomentum (const std::vector< const reco::Candidate * > dl) const

Private Attributes
std::vector< BPHRecoConstCandPtr >	cList

std::map< const reco::Candidate , const reco::Candidate >	clonesMap

std::map< std::string, BPHRecoConstCandPtr >	cMap

pat::CompositeCandidate	compCand

std::vector< Component >	compList

std::vector< const reco::Candidate * >	dFull

std::vector< const reco::Candidate * >	dList

std::map< std::string, const reco::Candidate * >	dMap

std::vector< std::string >	nComp

std::vector< std::string >	nList

bool	oldMom

Friends
class	BPHRecoBuilder

Detailed Description

Description: Lowest-level base class to contain decay products and compute total momentum

Author: Paolo Ronchese INFN Padova

Definition at line 36 of file BPHDecayMomentum.h.

Constructor & Destructor Documentation

BPHDecayMomentum::~BPHDecayMomentum ( )

virtual

Constructors are protected this object can exist only as part of a derived classDestructor

Definition at line 64 of file BPHDecayMomentum.cc.

References dList, and gen::n.

                                     {
   // delete all clones
   int n = dList.size();
   while ( n-- ) delete dList[n];
 }

BPHDecayMomentum::BPHDecayMomentum ( )

protected

Definition at line 35 of file BPHDecayMomentum.cc.

References clonesList(), cMap, dCompList(), dList, and oldMom.

                                   :
  oldMom( true ) {
   dList.reserve( 2 );
 }

BPHDecayMomentum::BPHDecayMomentum ( const std::map< std::string, Component > & daugMap )

protected

BPHDecayMomentum::BPHDecayMomentum	(	const std::map< std::string, Component > &	daugMap,
		const std::map< std::string, BPHRecoConstCandPtr >	compMap
	)

protected

Member Function Documentation

virtual void BPHDecayMomentum::addP	(	const std::string &	name,
		const reco::Candidate *	daug,
		double	mass = `-1.0`
	)

protectedvirtual

add a simple particle giving it a name particles are cloned, eventually specifying a different mass

Referenced by componentList().

virtual void BPHDecayMomentum::addP	(	const std::string &	name,
		const BPHRecoConstCandPtr &	comp
	)

protectedvirtual

add a previously reconstructed particle giving it a name

void BPHDecayMomentum::clonesList ( const std::map< std::string, Component > & daugMap )

private

Definition at line 190 of file BPHDecayMomentum.cc.

References reco::Candidate::begin(), BPHDecayMomentum::Component::cand, reco::Candidate::clone(), clonesMap, AlCaHLTBitMon_QueryRunRegistry::comp, compList, dList, dMap, mps_splice::entry, mps_fire::i, ResonanceBuilder::mass, BPHDecayMomentum::Component::mass, gen::n, dataset::name, nList, and reco::Candidate::setMass().

Referenced by BPHDecayMomentum().

                                                                         {
   int n = daugMap.size();
   dList.resize( n );
   nList.resize( n );
   // reset and fill a list
   // to be used in the creation of other bases of BPHRecoCandidate
   compList.clear();
   compList.reserve( n );
   // loop over daughters
   int i = 0;
   double mass;
   reco::Candidate* dnew;
   map<string,Component>::const_iterator iter = daugMap.begin();
   map<string,Component>::const_iterator iend = daugMap.end();
   while ( iter != iend ) {
     const pair<string,Component>& entry = *iter++;
     const Component& comp = entry.second;
     const reco::Candidate* cand = comp.cand;
     // store component for usage
     // in the creation of other bases of BPHRecoCandidate
     compList.push_back( comp );
     // clone particle and store it with its name
     dList[i] = dnew = cand->clone();
     const string& name = nList[i++] = entry.first;
     dMap[name] = dnew;
     clonesMap[dnew] = cand;
     // set daughter mass if requested
     mass = comp.mass;
     if ( mass > 0 ) dnew->setMass( mass );
   }
   return;
 }

const vector< string > & BPHDecayMomentum::compNames ( ) const

virtual

get the list of names of previously reconstructed particles e.g. in JPsi -> mu+mu- returns an empty list in B+ -> JPsi K+ returns the name used for the JPsi in Bs -> JPsi Phi returns the two names used for JPsi and Phi

Definition at line 84 of file BPHDecayMomentum.cc.

References nComp.

                                                         {
   return nComp;
 }

const vector< BPHDecayMomentum::Component > & BPHDecayMomentum::componentList ( ) const

protected

Definition at line 147 of file BPHDecayMomentum.cc.

References addP(), cList, reco::Candidate::clone(), clonesMap, cMap, AlCaHLTBitMon_QueryRunRegistry::comp, compList, dList, dMap, ResonanceBuilder::mass, dataset::name, nComp, nList, reco::Candidate::setMass(), and setNotUpdated().

Referenced by BPHDecayVertex::BPHDecayVertex(), and BPHKinematicFit::BPHKinematicFit().

                                                    {
   // return an object filled in the constructor
   // to be used in the creation of other bases of BPHRecoCandidate
   return compList;
 }

const pat::CompositeCandidate & BPHDecayMomentum::composite ( ) const

virtual

get a composite by the simple sum of simple particles

Operations

Reimplemented in BPHPlusMinusCandidate.

Definition at line 73 of file BPHDecayMomentum.cc.

References compCand, computeMomentum(), and oldMom.

Referenced by BPHPlusMinusCandidate::composite(), and BPHMassSelect::~BPHMassSelect().

                                                                {
   if ( oldMom ) computeMomentum();
   return compCand;
 }

void BPHDecayMomentum::computeMomentum ( ) const

privatevirtual

Definition at line 263 of file BPHDecayMomentum.cc.

References reco::CompositeCandidate::clearDaughters(), compCand, dFull, fillDaug(), oldMom, AddFourMomenta::set(), and sumMomentum().

Referenced by composite(), and daughFull().

                                              {
   // reset full list of daughters
   dFull.clear();
   fillDaug( dFull );
   // reset and fill pat::CompositeCandidate
   compCand.clearDaughters();
   sumMomentum( dFull );
   // compute the total momentum
   AddFourMomenta addP4;
   addP4.set( compCand );
   oldMom = false;
   return;
 }

const vector< BPHRecoConstCandPtr > & BPHDecayMomentum::daughComp ( ) const

virtual

get the list of previously reconstructed particles e.g. in JPsi -> mu+mu- returns an empty list in B+ -> JPsi K+ returns the JPsi in Bs -> JPsi Phi returns the JPsi and Phi

Definition at line 111 of file BPHDecayMomentum.cc.

References cList.

Referenced by BPHDecayVertex::BPHDecayVertex(), BPHKinematicFit::BPHKinematicFit(), and BPHKinematicFit::fitMomentum().

                                                                      {
   // return the list of previously reconstructed particles
   return cList;
 }

const vector< const reco::Candidate * > & BPHDecayMomentum::daughFull ( ) const

virtual

get the full list of simple particles produced in the decay, directly or through cascade decays e.g. in B+ -> JPsi K+ ; JPsi -> mu+mu- returns the mu+, mu-, K+ in Bs -> JPsi Phi; JPsi -> mu+mu-; Phi -> K+K- returns the mu+, mu-, K+, K- (clones are actually returned)

Definition at line 94 of file BPHDecayMomentum.cc.

References computeMomentum(), dFull, and oldMom.

Referenced by BPHKinematicFit::buildParticles(), BPHKinematicFit::kinematicTree(), BPHKinematicFit::kinParticles(), and BPHDecayVertex::tTracks().

                                                                       {
   // compute total momentum to update the full list of decay products
   if ( oldMom ) computeMomentum();
   return dFull;
 }

const vector< const reco::Candidate * > & BPHDecayMomentum::daughters ( ) const

virtual

get the list of simple particles directly produced in the decay e.g. in JPsi -> mu+mu- returns the two muons in B+ -> JPsi K+ returns the K+ in Bs -> JPsi Phi returns an empty list (clones are actually returned)

Definition at line 89 of file BPHDecayMomentum.cc.

References dList.

Referenced by BPHDecayVertex::BPHDecayVertex(), BPHKinematicFit::BPHKinematicFit(), BPHPlusMinusVertex::cAppInRPhi(), BPHWriteSpecificDecay::fill(), BPHKinematicFit::fitMomentum(), BPHKinematicFit::p4(), BPHPlusMinusCandidate::phiDiff(), and BPHPlusMinusCandidate::~BPHPlusMinusCandidate().

                                                                       {
   return dList;
 }

const vector< string > & BPHDecayMomentum::daugNames ( ) const

virtual

get the list of names of simple particles directly produced in the decay e.g. in JPsi -> mu+mu- returns the two names used for muons in B+ -> JPsi K+ returns the name used for the K+ in Bs -> JPsi Phi returns an empty list

Definition at line 79 of file BPHDecayMomentum.cc.

References nList.

Referenced by BPHKinematicFit::kinematicTree().

                                                         {
   return nList;
 }

void BPHDecayMomentum::dCompList ( )

private

Definition at line 224 of file BPHDecayMomentum.cc.

References cList, clonesMap, cMap, AlCaHLTBitMon_QueryRunRegistry::comp, mps_splice::entry, mps_fire::i, gen::n, and nComp.

Referenced by BPHDecayMomentum().

                                  {
   // fill lists of previously reconstructed particles and their names
   // and retrieve cascade decay products
   int n = cMap.size();
   cList.resize( n );
   nComp.resize( n );
   int i = 0;
   map<string,BPHRecoConstCandPtr>::const_iterator iter = cMap.begin();
   map<string,BPHRecoConstCandPtr>::const_iterator iend = cMap.end();
   while ( iter != iend ) {
     const pair<string,BPHRecoConstCandPtr>& entry = *iter++;
     nComp[i] = entry.first;
     BPHRecoConstCandPtr comp = entry.second;
     cList[i++] = comp;
     clonesMap.insert( comp->clonesMap.begin(), comp->clonesMap.end() );
   }
   return;
 }

void BPHDecayMomentum::fillDaug ( std::vector< const reco::Candidate * > & ad ) const

privatevirtual

Definition at line 253 of file BPHDecayMomentum.cc.

References cList, dList, and gen::n.

Referenced by computeMomentum().

                                                                         {
   // recursively fill the list of simple particles, produced
   // directly or in cascade decays
   ad.insert( ad.end(), dList.begin(), dList.end() );
   int n = cList.size();
   while ( n-- ) cList[n]->fillDaug( ad );
   return;
 }

BPHRecoConstCandPtr BPHDecayMomentum::getComp ( const std::string & name ) const

virtual

get a previously reconstructed particle from the name return null pointer if not found

Definition at line 132 of file BPHDecayMomentum.cc.

References cMap, AlCaHLTBitMon_QueryRunRegistry::comp, getComp(), genParticles_cff::map, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by getComp(), getDaug(), and BPHKinematicFit::kinematicTree().

                                                                         {
   // return a previously reconstructed particle from the name
   // return null pointer if not found
   string::size_type pos = name.find( "/" );
   if ( pos != string::npos ) {
     const BPHRecoCandidate* comp = getComp( name.substr( 0, pos ) ).get();
     return ( comp == 0 ? 0 : comp->getComp( name.substr( pos + 1 ) ) );
   }
   map<const string,
       BPHRecoConstCandPtr>::const_iterator iter = cMap.find( name );
   return ( iter != cMap.end() ? iter->second : 0 );
 }

const reco::Candidate * BPHDecayMomentum::getDaug ( const std::string & name ) const

virtual

get a simple particle from the name return null pointer if not found

Definition at line 117 of file BPHDecayMomentum.cc.

References AlCaHLTBitMon_QueryRunRegistry::comp, dMap, getComp(), getDaug(), genParticles_cff::map, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by BPHKx0ToKPiBuilder::build(), BPHWriteSpecificDecay::fill(), getDaug(), BPHOniaToMuMuBuilder::getOriginalCandidate(), BPHKinematicFit::kinParticles(), and BPHMassSymSelect::~BPHMassSymSelect().

                                                   {
   // return a simple particle from the name
   // return null pointer if not found
   string::size_type pos = name.find( "/" );
   if ( pos != string::npos ) {
     const BPHRecoCandidate* comp = getComp( name.substr( 0, pos ) ).get();
     return ( comp == 0 ? 0 : comp->getDaug( name.substr( pos + 1 ) ) );
   }
   map<const string,
       const reco::Candidate*>::const_iterator iter = dMap.find( name );
   return ( iter != dMap.end() ? iter->second : 0 );
 }

const reco::Candidate * BPHDecayMomentum::originalReco ( const reco::Candidate * daug ) const

virtual

get the original particle from the clone

Definition at line 101 of file BPHDecayMomentum.cc.

References clonesMap, and genParticles_cff::map.

Referenced by BPHDecayVertex::BPHDecayVertex(), BPHKinematicFit::BPHKinematicFit(), BPHKx0ToKPiBuilder::build(), BPHWriteSpecificDecay::fill(), BPHOniaToMuMuBuilder::getOriginalCandidate(), and BPHMassSymSelect::~BPHMassSymSelect().

                                                          {
   // return the original particle for a given clone
   // return null pointer if not found
   map<const reco::Candidate*,
       const reco::Candidate*>::const_iterator iter = clonesMap.find( daug );
   return ( iter != clonesMap.end() ? iter->second : 0 );
 }

void BPHDecayMomentum::setNotUpdated ( ) const

protectedvirtual

Reimplemented in BPHKinematicFit, BPHDecayVertex, BPHPlusMinusCandidate, and BPHPlusMinusVertex.

Definition at line 184 of file BPHDecayMomentum.cc.

References oldMom.

Referenced by componentList(), and BPHDecayVertex::setNotUpdated().

                                            {
   oldMom = true;
   return;
 }

void BPHDecayMomentum::sumMomentum ( const std::vector< const reco::Candidate * > dl ) const

privatevirtual

Definition at line 244 of file BPHDecayMomentum.cc.

References reco::CompositeCandidate::addDaughter(), compCand, and gen::n.

Referenced by computeMomentum().

                                                                      {
   // add the particles to pat::CompositeCandidate
   int n = dl.size();
   while ( n-- ) compCand.addDaughter( *dl[n] );  
   return;
 }