#include <BPHKinematicFit.h>

Inheritance diagram for BPHKinematicFit:

Classes
struct	FlyingParticle

Public Member Functions
	BPHKinematicFit (const BPHKinematicFit &x)=delete

double	constrMass () const
	retrieve the constraint More...

double	constrSigma () const

virtual const RefCountedKinematicVertex	currentDecayVertex () const

virtual const RefCountedKinematicParticle	currentParticle () const
	get current particle More...

bool	getIndependentFit (const std::string &name, double &mass, double &sigma) const
	retrieve independent fit flag More...

double	getMassSigma (const reco::Candidate *cand) const
	retrieve particle mass sigma More...

virtual bool	isEmpty () const
	get fit status More...

virtual bool	isValidFit () const

virtual const RefCountedKinematicTree &	kinematicTree () const
	perform the kinematic fit and get the result More...

virtual const RefCountedKinematicTree &	kinematicTree (const std::string &name, double mass, double sigma) const

virtual const RefCountedKinematicTree &	kinematicTree (const std::string &name, double mass) const

virtual const RefCountedKinematicTree &	kinematicTree (const std::string &name) const

virtual const RefCountedKinematicTree &	kinematicTree (const std::string &name, KinematicConstraint *kc) const

virtual const RefCountedKinematicTree &	kinematicTree (const std::string &name, MultiTrackKinematicConstraint *kc) const

virtual const std::vector< RefCountedKinematicParticle > &	kinParticles () const
	get kinematic particles More...

virtual std::vector< RefCountedKinematicParticle >	kinParticles (const std::vector< std::string > &names) const

virtual ParticleMass	mass () const

BPHKinematicFit &	operator= (const BPHKinematicFit &x)=delete

virtual const math::XYZTLorentzVector &	p4 () const
	compute total momentum after the fit More...

virtual void	resetKinematicFit () const
	reset the kinematic fit More...

void	setConstraint (double mass, double sigma)
	apply a mass constraint More...

void	setIndependentFit (const std::string &name, bool flag=true, double mass=-1.0, double sigma=-1.0)
	set a decaying daughter as an unique particle fitted independently More...

virtual const RefCountedKinematicVertex	topDecayVertex () const

virtual const RefCountedKinematicParticle	topParticle () const
	get top particle More...

	~BPHKinematicFit () override=default

Public Member Functions inherited from BPHDecayVertex
	BPHDecayVertex (const BPHDecayVertex &x)=delete

const BPHEventSetupWrapper *	getEventSetup () const
	retrieve EventSetup More...

char	getTMode (const reco::Candidate *cand) const
	get Track mode for a daughter More...

const reco::Track *	getTrack (const reco::Candidate *cand) const
	get Track for a daughter More...

const std::string &	getTrackSearchList (const reco::Candidate *cand) const
	retrieve track search list More...

reco::TransientTrack *	getTransientTrack (const reco::Candidate *cand) const
	get TransientTrack for a daughter More...

BPHDecayVertex &	operator= (const BPHDecayVertex &x)=delete

const std::vector< const reco::Track * > &	tracks () const
	get list of Tracks More...

const std::vector< reco::TransientTrack > &	transientTracks () const
	get list of TransientTracks More...

virtual bool	validTracks () const
	check for valid reconstructed vertex More...

virtual bool	validVertex () const

virtual const reco::Vertex &	vertex (VertexFitter< 5 > fitter=nullptr, const reco::BeamSpot bs=nullptr, const GlobalPoint priorPos=nullptr, const GlobalError priorError=nullptr) const
	get reconstructed vertex More...

	~BPHDecayVertex () override

Public Member Functions inherited from BPHDecayMomentum
	BPHDecayMomentum (const BPHDecayMomentum &x)=delete

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

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

const std::map< std::string, const reco::Candidate * > &	daugMap () 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

BPHDecayMomentum &	operator= (const BPHDecayMomentum &x)=delete

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	addK (const std::string &name, const reco::Candidate *daug, double mass=-1.0, double sigma=-1.0)

virtual void	addK (const std::string &name, const reco::Candidate *daug, const std::string &searchList, double mass=-1.0, double sigma=-1.0)

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

	BPHKinematicFit (int daugNum=2, int compNum=2)

	BPHKinematicFit (const BPHKinematicFit *ptr)

void	setNotUpdated () const override

Protected Member Functions inherited from BPHDecayVertex
virtual void	addV (const std::string &name, const reco::Candidate *daug, const std::string &searchList, double mass)

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

	BPHDecayVertex (const BPHEventSetupWrapper *es, int daugNum=2, int compNum=2)

	BPHDecayVertex (const BPHDecayVertex ptr, const BPHEventSetupWrapper es)

Protected Member Functions inherited from BPHDecayMomentum
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)

	BPHDecayMomentum (int daugNum=2, int compNum=2)

	BPHDecayMomentum (const std::map< std::string, Component > &daugMap, int compNum=2)

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

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

virtual void	fill (BPHRecoCandidate *ptr, int level) const =0

Private Member Functions
virtual void	addParticles (std::vector< RefCountedKinematicParticle > &kl, std::map< const reco::Candidate , RefCountedKinematicParticle > &km, std::map< const BPHRecoCandidate , RefCountedKinematicParticle > &cm) const

virtual void	buildParticles () const

virtual void	fitMomentum () const

virtual void	getParticles (const std::string &moth, const std::string &daug, std::vector< RefCountedKinematicParticle > &kl, std::set< RefCountedKinematicParticle > &ks, const BPHKinematicFit *curr) const

virtual void	getParticles (const std::string &moth, const std::vector< std::string > &daug, std::vector< RefCountedKinematicParticle > &kl, std::set< RefCountedKinematicParticle > &ks, const BPHKinematicFit *curr) const

virtual const RefCountedKinematicTree &	kinematicTree (const std::vector< RefCountedKinematicParticle > &kPart, MultiTrackKinematicConstraint *kc) const

virtual unsigned int	numParticles (const BPHKinematicFit *cand=nullptr) const

virtual const BPHKinematicFit *	splitKP (const std::string &name, std::vector< RefCountedKinematicParticle > kComp, std::vector< RefCountedKinematicParticle > kTail=nullptr) const

Static Private Member Functions
static void	insertParticle (RefCountedKinematicParticle &kp, std::vector< RefCountedKinematicParticle > &kl, std::set< RefCountedKinematicParticle > &ks)

Private Attributes
std::vector< RefCountedKinematicParticle >	allParticles

std::map< const BPHRecoCandidate *, FlyingParticle >	cKinP

std::map< const reco::Candidate *, double >	dMSig

std::map< const BPHRecoCandidate *, RefCountedKinematicParticle >	kCDMap

std::map< const reco::Candidate *, RefCountedKinematicParticle >	kinMap

RefCountedKinematicTree	kinTree

double	massConst

double	massSigma

bool	oldFit

bool	oldKPs

bool	oldMom

std::vector< BPHRecoConstCandPtr >	tmpList

math::XYZTLorentzVector	totalMomentum

Detailed Description

Description: Highest-level base class to encapsulate kinematic fit operations

Author: Paolo Ronchese INFN Padova high-level base class to perform a kinematic fit

Definition at line 38 of file BPHKinematicFit.h.

Constructor & Destructor Documentation

◆ BPHKinematicFit() [1/3]

BPHKinematicFit::BPHKinematicFit ( const BPHKinematicFit & x )

delete

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

◆ ~BPHKinematicFit()

BPHKinematicFit::~BPHKinematicFit ( )

overridedefault

Destructor

◆ BPHKinematicFit() [2/3]

BPHKinematicFit::BPHKinematicFit	(	int	daugNum = `2`,
		int	compNum = `2`
	)

protected

Definition at line 39 of file BPHKinematicFit.cc.

     : BPHDecayMomentum(daugNum, compNum),
       BPHDecayVertex(nullptr),
       massConst(-1.0),
       massSigma(-1.0),
       oldKPs(true),
       oldFit(true),
       oldMom(true),
       kinTree(nullptr) {}

◆ BPHKinematicFit() [3/3]

BPHKinematicFit::BPHKinematicFit ( const BPHKinematicFit * ptr )

protected

Definition at line 49 of file BPHKinematicFit.cc.

References HltBtagPostValidation_cff::c, cKinP, BPHDecayMomentum::cMap, BPHDecayMomentum::componentList(), BPHDecayMomentum::daughComp(), BPHDecayMomentum::daughters(), BPHDecayMomentum::dMap, dMSig, mps_fire::i, dqmiolumiharvest::j, visualization-live-secondInstance_cfg::m, dqmiodumpmetadata::n, and BPHDecayMomentum::originalReco().

     : BPHDecayVertex(ptr, nullptr),
       massConst(-1.0),
       massSigma(-1.0),
       oldKPs(true),
       oldFit(true),
       oldMom(true),
       kinTree(nullptr) {
   map<const reco::Candidate*, const reco::Candidate*> iMap;
   const vector<const reco::Candidate*>& daug = daughters();
   const vector<Component>& list = ptr->componentList();
   int i;
   int n = daug.size();
   for (i = 0; i < n; ++i) {
     const reco::Candidate* cand = daug[i];
     iMap[originalReco(cand)] = cand;
   }
   for (i = 0; i < n; ++i) {
     const Component& c = list[i];
     dMSig[iMap[c.cand]] = c.msig;
   }
   const vector<BPHRecoConstCandPtr>& dComp = daughComp();
   int j;
   int m = dComp.size();
   for (j = 0; j < m; ++j) {
     const BPHRecoCandidate* rc = dComp[j].get();
     const map<const reco::Candidate*, double>& dMap = rc->dMSig;
     const map<const BPHRecoCandidate*, FlyingParticle>& cMap = rc->cKinP;
     dMSig.insert(dMap.begin(), dMap.end());
     cKinP.insert(cMap.begin(), cMap.end());
     cKinP[rc];
   }
 }

Member Function Documentation

◆ addK() [1/3]

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

protectedvirtual

Referenced by BPHRecoCandidate::add().

◆ addK() [2/3]

virtual void BPHKinematicFit::addK	(	const std::string &	name,
		const reco::Candidate *	daug,
		const std::string &	searchList,
		double	mass = `-1.0`,
		double	sigma = `-1.0`
	)

protectedvirtual

◆ addK() [3/3]

virtual void BPHKinematicFit::addK	(	const std::string &	name,
		const BPHRecoConstCandPtr &	comp
	)

protectedvirtual

◆ addParticles()

void BPHKinematicFit::addParticles	(	std::vector< RefCountedKinematicParticle > &	kl,
		std::map< const reco::Candidate *, RefCountedKinematicParticle > &	km,
		std::map< const BPHRecoCandidate *, RefCountedKinematicParticle > &	cm
	)		const

privatevirtual

Definition at line 371 of file BPHKinematicFit.cc.

References cKinP, BPHRecoCandidate::clone(), AlCaHLTBitMon_QueryRunRegistry::comp, BPHDecayMomentum::daughComp(), BPHDecayMomentum::daughters(), dMSig, personalPlayback::fp, BPHDecayVertex::getTransientTrack(), isEmpty(), isValidFit(), visualization-live-secondInstance_cfg::m, mass(), dqmiodumpmetadata::n, KinematicParticleFactoryFromTransientTrack::particle(), setConstraint(), tmpList, and topParticle().

Referenced by buildParticles().

                                                                                                         {
   const vector<const reco::Candidate*>& daug = daughters();
   KinematicParticleFactoryFromTransientTrack pFactory;
   int n = daug.size();
   float chi = 0.0;
   float ndf = 0.0;
   while (n--) {
     const reco::Candidate* cand = daug[n];
     ParticleMass mass = cand->mass();
     float sigma = dMSig.find(cand)->second;
     if (sigma < 0)
       sigma = 1.0e-7;
     reco::TransientTrack* tt = getTransientTrack(cand);
     if (tt != nullptr)
       kl.push_back(km[cand] = pFactory.particle(*tt, mass, chi, ndf, sigma));
   }
   const vector<BPHRecoConstCandPtr>& comp = daughComp();
   int m = comp.size();
   while (m--) {
     const BPHRecoCandidate* cptr = comp[m].get();
     const FlyingParticle& fp = cKinP.at(cptr);
     if (fp.flag) {
       BPHRecoCandidate* tptr = cptr->clone();
       double mass = fp.mass;
       double sigma = fp.sigma;
       if (mass > 0.0)
         tptr->setConstraint(mass, sigma);
       tmpList.push_back(BPHRecoConstCandPtr(tptr));
       if (tptr->isEmpty())
         return;
       if (!tptr->isValidFit())
         return;
       kl.push_back(cm[cptr] = tptr->topParticle());
     } else {
       cptr->addParticles(kl, km, cm);
     }
   }
   return;
 }

◆ buildParticles()

void BPHKinematicFit::buildParticles ( ) const

privatevirtual

Definition at line 361 of file BPHKinematicFit.cc.

References addParticles(), allParticles, BPHDecayMomentum::daughFull(), kCDMap, kinMap, oldKPs, and findQualityFiles::size.

Referenced by kinParticles().

                                            {
   kinMap.clear();
   kCDMap.clear();
   allParticles.clear();
   allParticles.reserve(daughFull().size());
   addParticles(allParticles, kinMap, kCDMap);
   oldKPs = false;
   return;
 }

◆ constrMass()

double BPHKinematicFit::constrMass ( ) const

retrieve the constraint

Definition at line 95 of file BPHKinematicFit.cc.

References massConst.

Referenced by BPHRecoCandidate::fill().

95 { return massConst; }

BPHKinematicFit::massConst

double massConst

Definition: BPHKinematicFit.h:125

◆ constrSigma()

double BPHKinematicFit::constrSigma ( ) const

Definition at line 97 of file BPHKinematicFit.cc.

References massSigma.

Referenced by BPHRecoCandidate::fill().

97 { return massSigma; }

BPHKinematicFit::massSigma

double massSigma

Definition: BPHKinematicFit.h:126

◆ currentDecayVertex()

const RefCountedKinematicVertex BPHKinematicFit::currentDecayVertex ( ) const

virtual

Definition at line 272 of file BPHKinematicFit.cc.

References isEmpty(), and kinTree.

                                                                           {
   if (isEmpty())
     return RefCountedKinematicVertex(nullptr);
   return kinTree->currentDecayVertex();
 }

◆ currentParticle()

const RefCountedKinematicParticle BPHKinematicFit::currentParticle ( ) const

virtual

get current particle

Definition at line 266 of file BPHKinematicFit.cc.

References isEmpty(), and kinTree.

                                                                          {
   if (isEmpty())
     return RefCountedKinematicParticle(nullptr);
   return kinTree->currentParticle();
 }

◆ fitMomentum()

void BPHKinematicFit::fitMomentum ( ) const

privatevirtual

Definition at line 540 of file BPHKinematicFit.cc.

References AlCaHLTBitMon_QueryRunRegistry::comp, BPHDecayMomentum::daughComp(), BPHDecayMomentum::daughters(), MillePedeFileConverter_cfg::e, KinematicState::globalMomentum(), isValidFit(), visualization-live-secondInstance_cfg::m, KinematicState::mass(), dqmiodumpmetadata::n, oldMom, mathSSE::sqrt(), topParticle(), totalMomentum, x, PV3DBase< T, PVType, FrameType >::x(), y, PV3DBase< T, PVType, FrameType >::y(), z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by p4().

                                         {
   if (isValidFit()) {
     const KinematicState& ks = topParticle()->currentState();
     GlobalVector tm = ks.globalMomentum();
     double x = tm.x();
     double y = tm.y();
     double z = tm.z();
     double m = ks.mass();
     double e = sqrt((x * x) + (y * y) + (z * z) + (m * m));
     totalMomentum.SetPxPyPzE(x, y, z, e);
   } else {
     edm::LogPrint("FitNotFound") << "BPHKinematicFit::fitMomentum: "
                                  << "simple momentum sum computed";
     math::XYZTLorentzVector tm;
     const vector<const reco::Candidate*>& daug = daughters();
     int n = daug.size();
     while (n--)
       tm += daug[n]->p4();
     const vector<BPHRecoConstCandPtr>& comp = daughComp();
     int m = comp.size();
     while (m--)
       tm += comp[m]->p4();
     totalMomentum = tm;
   }
   oldMom = false;
   return;
 }

◆ getIndependentFit()

bool BPHKinematicFit::getIndependentFit	(	const std::string &	name,
		double &	mass,
		double &	sigma
	)		const

retrieve independent fit flag

Definition at line 316 of file BPHKinematicFit.cc.

References cKinP, AlCaHLTBitMon_QueryRunRegistry::comp, BPHDecayMomentum::getComp(), mass(), and Skims_PA_cff::name.

Referenced by BPHRecoCandidate::fill().

                                                                                              {
   const BPHRecoCandidate* comp = getComp(name).get();
   map<const BPHRecoCandidate*, FlyingParticle>::const_iterator iter = cKinP.find(comp);
   if ((iter != cKinP.end()) && iter->second.flag) {
     mass = iter->second.mass;
     sigma = iter->second.sigma;
     return true;
   }
   return false;
 }

◆ getMassSigma()

double BPHKinematicFit::getMassSigma ( const reco::Candidate * cand ) const

retrieve particle mass sigma

Definition at line 310 of file BPHKinematicFit.cc.

References dMSig.

Referenced by BPHRecoCandidate::fill().

                                                                     {
   map<const reco::Candidate*, double>::const_iterator iter = dMSig.find(cand);
   return (iter != dMSig.end() ? iter->second : -1);
 }

◆ getParticles() [1/2]

virtual void BPHKinematicFit::getParticles	(	const std::string &	moth,
		const std::string &	daug,
		std::vector< RefCountedKinematicParticle > &	kl,
		std::set< RefCountedKinematicParticle > &	ks,
		const BPHKinematicFit *	curr
	)		const

privatevirtual

◆ getParticles() [2/2]

virtual void BPHKinematicFit::getParticles	(	const std::string &	moth,
		const std::vector< std::string > &	daug,
		std::vector< RefCountedKinematicParticle > &	kl,
		std::set< RefCountedKinematicParticle > &	ks,
		const BPHKinematicFit *	curr
	)		const

privatevirtual

◆ insertParticle()

void BPHKinematicFit::insertParticle	(	RefCountedKinematicParticle &	kp,
		std::vector< RefCountedKinematicParticle > &	kl,
		std::set< RefCountedKinematicParticle > &	ks
	)

staticprivate

Definition at line 479 of file BPHKinematicFit.cc.

                                                                            {
   if (ks.find(kp) != ks.end())
     return;
   kl.push_back(kp);
   ks.insert(kp);
   return;
 }

◆ isEmpty()

bool BPHKinematicFit::isEmpty ( void ) const

virtual

get fit status

Definition at line 251 of file BPHKinematicFit.cc.

References kinematicTree(), and kinTree.

Referenced by addParticles(), plotting.Plot::clone(), currentDecayVertex(), currentParticle(), topDecayVertex(), and topParticle().

                                     {
   kinematicTree();
   if (kinTree.get() == nullptr)
     return true;
   return kinTree->isEmpty();
 }

◆ isValidFit()

bool BPHKinematicFit::isValidFit ( ) const

virtual

Definition at line 258 of file BPHKinematicFit.cc.

References topParticle().

Referenced by addParticles(), and fitMomentum().

                                        {
   const RefCountedKinematicParticle kPart = topParticle();
   if (kPart.get() == nullptr)
     return false;
   return kPart->currentState().isValid();
 }

◆ kinematicTree() [1/7]

const RefCountedKinematicTree & BPHKinematicFit::kinematicTree ( ) const

virtual

perform the kinematic fit and get the result

Definition at line 156 of file BPHKinematicFit.cc.

References kinTree, massConst, massSigma, and oldFit.

Referenced by isEmpty().

                                                                     {
   if (oldFit)
     return kinematicTree("", massConst, massSigma);
   return kinTree;
 }

◆ kinematicTree() [2/7]

virtual const RefCountedKinematicTree& BPHKinematicFit::kinematicTree	(	const std::string &	name,
		double	mass,
		double	sigma
	)		const

virtual

◆ kinematicTree() [3/7]

virtual const RefCountedKinematicTree& BPHKinematicFit::kinematicTree	(	const std::string &	name,
		double	mass
	)		const

virtual

◆ kinematicTree() [4/7]

virtual const RefCountedKinematicTree& BPHKinematicFit::kinematicTree ( const std::string & name ) const

virtual

◆ kinematicTree() [5/7]

virtual const RefCountedKinematicTree& BPHKinematicFit::kinematicTree	(	const std::string &	name,
		KinematicConstraint *	kc
	)		const

virtual

◆ kinematicTree() [6/7]

virtual const RefCountedKinematicTree& BPHKinematicFit::kinematicTree	(	const std::string &	name,
		MultiTrackKinematicConstraint *	kc
	)		const

virtual

◆ kinematicTree() [7/7]

const RefCountedKinematicTree & BPHKinematicFit::kinematicTree	(	const std::vector< RefCountedKinematicParticle > &	kPart,
		MultiTrackKinematicConstraint *	kc
	)		const

privatevirtual

Definition at line 527 of file BPHKinematicFit.cc.

References cppFunctionSkipper::exception, KinematicConstrainedVertexFitter::fit(), and kinTree.

                                                                                                        {
   try {
     KinematicConstrainedVertexFitter cvf;
     kinTree = cvf.fit(kPart, kc);
   } catch (std::exception const&) {
     edm::LogPrint("FitFailed") << "BPHKinematicFit::kinematicTree: "
                                << "kin fit reset";
     kinTree = RefCountedKinematicTree(nullptr);
   }
   return kinTree;
 }

◆ kinParticles() [1/2]

const vector< RefCountedKinematicParticle > & BPHKinematicFit::kinParticles ( ) const

virtual

get kinematic particles

Definition at line 121 of file BPHKinematicFit.cc.

References allParticles, buildParticles(), and oldKPs.

Referenced by splitKP().

                                                                                {
   if (oldKPs)
     buildParticles();
   return allParticles;
 }

◆ kinParticles() [2/2]

virtual std::vector<RefCountedKinematicParticle> BPHKinematicFit::kinParticles ( const std::vector< std::string > & names ) const

virtual

◆ mass()

ParticleMass BPHKinematicFit::mass ( ) const

virtual

Definition at line 292 of file BPHKinematicFit.cc.

References KinematicState::isValid(), KinematicState::mass(), and topParticle().

Referenced by Particle.Particle::__str__(), DiObject.DiMuon::__str__(), BPHRecoCandidate::add(), addParticles(), BPHRecoCandidate::build(), BPHDecayToV0DiffMassBuilder::buildCandidate(), BPHRecoCandidate::fill(), getIndependentFit(), setConstraint(), and setIndependentFit().

                                          {
   const RefCountedKinematicParticle kPart = topParticle();
   if (kPart.get() == nullptr)
     return -1.0;
   const KinematicState kStat = kPart->currentState();
   if (kStat.isValid())
     return kStat.mass();
   return -1.0;
 }

◆ numParticles()

unsigned int BPHKinematicFit::numParticles ( const BPHKinematicFit * cand = nullptr ) const

privatevirtual

Definition at line 463 of file BPHKinematicFit.cc.

References AlCaHLTBitMon_QueryRunRegistry::comp, mps_fire::i, and dqmiodumpmetadata::n.

Referenced by splitKP().

                                                                             {
   if (cand == nullptr)
     cand = this;
   unsigned int n = cand->daughters().size();
   const vector<string>& cnames = cand->compNames();
   int i = cnames.size();
   while (i) {
     const BPHRecoCandidate* comp = cand->getComp(cnames[--i]).get();
     if (cand->cKinP.at(comp).flag)
       ++n;
     else
       n += numParticles(comp);
   }
   return n;
 }

◆ operator=()

BPHKinematicFit& BPHKinematicFit::operator= ( const BPHKinematicFit & x )

delete

◆ p4()

const math::XYZTLorentzVector & BPHKinematicFit::p4 ( ) const

virtual

compute total momentum after the fit

Definition at line 303 of file BPHKinematicFit.cc.

References fitMomentum(), oldMom, and totalMomentum.

Referenced by Tau.Tau::dxy_approx(), Tau.Tau::dz(), and Lepton.Lepton::p4WithFSR().

                                                        {
   if (oldMom)
     fitMomentum();
   return totalMomentum;
 }

◆ resetKinematicFit()

void BPHKinematicFit::resetKinematicFit ( ) const

virtual

reset the kinematic fit

Definition at line 245 of file BPHKinematicFit.cc.

References oldFit, oldKPs, and oldMom.

Referenced by setNotUpdated().

                                               {
   oldKPs = oldFit = oldMom = true;
   return;
 }

◆ setConstraint()

void BPHKinematicFit::setConstraint	(	double	mass,
		double	sigma
	)

apply a mass constraint

Operations

Definition at line 87 of file BPHKinematicFit.cc.

References mass(), massConst, massSigma, oldFit, oldKPs, and oldMom.

Referenced by addParticles(), and BPHRecoCandidate::fill().

                                                              {
   oldKPs = oldFit = oldMom = true;
   massConst = mass;
   massSigma = sigma;
   return;
 }

◆ setIndependentFit()

void BPHKinematicFit::setIndependentFit	(	const std::string &	name,
		bool	flag = `true`,
		double	mass = `-1.0`,
		double	sigma = `-1.0`
	)

set a decaying daughter as an unique particle fitted independently

Definition at line 100 of file BPHKinematicFit.cc.

References cKinP, AlCaHLTBitMon_QueryRunRegistry::comp, RemoveAddSevLevel::flag, personalPlayback::fp, BPHDecayMomentum::getComp(), mass(), Skims_PA_cff::name, oldFit, oldKPs, and oldMom.

Referenced by BPHRecoCandidate::fill().

                                                                                                 {
   string::size_type pos = name.find('/');
   if (pos != string::npos) {
     edm::LogPrint("WrongRequest") << "BPHKinematicFit::setIndependentFit: "
                                   << "cascade decay specification not admitted " << name;
     return;
   }
   const BPHRecoCandidate* comp = getComp(name).get();
   if (comp == nullptr) {
     edm::LogPrint("ParticleNotFound") << "BPHKinematicFit::setIndependentFit: " << name << " daughter not found";
     return;
   }
   oldKPs = oldFit = oldMom = true;
   FlyingParticle& fp = cKinP[comp];
   fp.flag = flag;
   fp.mass = mass;
   fp.sigma = sigma;
   return;
 }

◆ setNotUpdated()

void BPHKinematicFit::setNotUpdated ( ) const

overrideprotectedvirtual

Reimplemented from BPHDecayVertex.

Reimplemented in BPHPlusMinusCandidate.

Definition at line 354 of file BPHKinematicFit.cc.

References resetKinematicFit(), and BPHDecayVertex::setNotUpdated().

Referenced by BPHPlusMinusCandidate::setNotUpdated().

                                           {
   BPHDecayVertex::setNotUpdated();
   resetKinematicFit();
   return;
 }

◆ splitKP()

const BPHKinematicFit * BPHKinematicFit::splitKP	(	const std::string &	name,
		std::vector< RefCountedKinematicParticle > *	kComp,
		std::vector< RefCountedKinematicParticle > *	kTail = `nullptr`
	)		const

privatevirtual

Definition at line 489 of file BPHKinematicFit.cc.

References allParticles, cKinP, AlCaHLTBitMon_QueryRunRegistry::comp, BPHDecayMomentum::getComp(), kinParticles(), kinTree, Skims_PA_cff::name, numParticles(), and oldFit.

                                                                                                   {
   kinTree = RefCountedKinematicTree(nullptr);
   oldFit = false;
   kinParticles();
   if (allParticles.size() != numParticles())
     return nullptr;
   kComp->clear();
   if (kTail == nullptr)
     kTail = kComp;
   else
     kTail->clear();
   if (name.empty()) {
     *kComp = allParticles;
     return this;
   }
   const BPHRecoCandidate* comp = getComp(name).get();
   int ns;
   if (comp != nullptr) {
     ns = (cKinP.at(comp).flag ? 1 : numParticles(comp));
   } else {
     edm::LogPrint("ParticleNotFound") << "BPHKinematicFit::splitKP: " << name << " daughter not found";
     *kTail = allParticles;
     return nullptr;
   }
   vector<string> nfull(2);
   nfull[0] = name;
   nfull[1] = "*";
   vector<RefCountedKinematicParticle> kPart = kinParticles(nfull);
   vector<RefCountedKinematicParticle>::const_iterator iter = kPart.begin();
   vector<RefCountedKinematicParticle>::const_iterator imid = iter + ns;
   vector<RefCountedKinematicParticle>::const_iterator iend = kPart.end();
   kComp->insert(kComp->end(), iter, imid);
   kTail->insert(kTail->end(), imid, iend);
   return comp;
 }

◆ topDecayVertex()

const RefCountedKinematicVertex BPHKinematicFit::topDecayVertex ( ) const

virtual

Definition at line 285 of file BPHKinematicFit.cc.

References isEmpty(), and kinTree.

                                                                       {
   if (isEmpty())
     return RefCountedKinematicVertex(nullptr);
   kinTree->movePointerToTheTop();
   return kinTree->currentDecayVertex();
 }

◆ topParticle()

const RefCountedKinematicParticle BPHKinematicFit::topParticle ( ) const

virtual

get top particle

Definition at line 279 of file BPHKinematicFit.cc.

References isEmpty(), and kinTree.

Referenced by addParticles(), fitMomentum(), isValidFit(), and mass().

                                                                      {
   if (isEmpty())
     return RefCountedKinematicParticle(nullptr);
   return kinTree->topParticle();
 }

Member Data Documentation

◆ allParticles

std::vector<RefCountedKinematicParticle> BPHKinematicFit::allParticles

mutableprivate

Definition at line 148 of file BPHKinematicFit.h.

Referenced by buildParticles(), kinParticles(), and splitKP().

◆ cKinP

std::map<const BPHRecoCandidate*, FlyingParticle> BPHKinematicFit::cKinP

private

Definition at line 137 of file BPHKinematicFit.h.

Referenced by addParticles(), BPHKinematicFit(), getIndependentFit(), setIndependentFit(), and splitKP().

◆ dMSig

std::map<const reco::Candidate*, double> BPHKinematicFit::dMSig

private

Definition at line 129 of file BPHKinematicFit.h.

Referenced by addParticles(), BPHKinematicFit(), and getMassSigma().

◆ kCDMap

std::map<const BPHRecoCandidate*, RefCountedKinematicParticle> BPHKinematicFit::kCDMap

mutableprivate

Definition at line 147 of file BPHKinematicFit.h.

Referenced by buildParticles().

◆ kinMap

std::map<const reco::Candidate*, RefCountedKinematicParticle> BPHKinematicFit::kinMap

mutableprivate

Definition at line 146 of file BPHKinematicFit.h.

Referenced by buildParticles().

◆ kinTree

RefCountedKinematicTree BPHKinematicFit::kinTree

mutableprivate

Definition at line 149 of file BPHKinematicFit.h.

Referenced by currentDecayVertex(), currentParticle(), isEmpty(), kinematicTree(), splitKP(), topDecayVertex(), and topParticle().

◆ massConst

double BPHKinematicFit::massConst

private

Definition at line 125 of file BPHKinematicFit.h.

Referenced by constrMass(), kinematicTree(), and setConstraint().

◆ massSigma

double BPHKinematicFit::massSigma

private

Definition at line 126 of file BPHKinematicFit.h.

Referenced by constrSigma(), kinematicTree(), and setConstraint().

◆ oldFit

bool BPHKinematicFit::oldFit

mutableprivate

Definition at line 144 of file BPHKinematicFit.h.

Referenced by kinematicTree(), resetKinematicFit(), setConstraint(), setIndependentFit(), and splitKP().

◆ oldKPs

bool BPHKinematicFit::oldKPs

mutableprivate

Definition at line 143 of file BPHKinematicFit.h.

Referenced by buildParticles(), kinParticles(), resetKinematicFit(), setConstraint(), and setIndependentFit().

◆ oldMom

bool BPHKinematicFit::oldMom

mutableprivate

Definition at line 145 of file BPHKinematicFit.h.

Referenced by fitMomentum(), p4(), resetKinematicFit(), setConstraint(), and setIndependentFit().

◆ tmpList

std::vector<BPHRecoConstCandPtr> BPHKinematicFit::tmpList

mutableprivate

Definition at line 140 of file BPHKinematicFit.h.

Referenced by addParticles().

◆ totalMomentum

math::XYZTLorentzVector BPHKinematicFit::totalMomentum

mutableprivate

Definition at line 150 of file BPHKinematicFit.h.

Referenced by fitMomentum(), and p4().

Classes

Public Member Functions

Protected Member Functions

Private Member Functions

Static Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ BPHKinematicFit() [1/3]

◆ ~BPHKinematicFit()

◆ BPHKinematicFit() [2/3]

◆ BPHKinematicFit() [3/3]

Member Function Documentation

◆ addK() [1/3]

◆ addK() [2/3]

◆ addK() [3/3]

◆ addParticles()

◆ buildParticles()

◆ constrMass()

◆ constrSigma()

◆ currentDecayVertex()

◆ currentParticle()

◆ fitMomentum()

◆ getIndependentFit()

◆ getMassSigma()

◆ getParticles() [1/2]

◆ getParticles() [2/2]

◆ insertParticle()

◆ isEmpty()

◆ isValidFit()

◆ kinematicTree() [1/7]

◆ kinematicTree() [2/7]

◆ kinematicTree() [3/7]

◆ kinematicTree() [4/7]

◆ kinematicTree() [5/7]

◆ kinematicTree() [6/7]

◆ kinematicTree() [7/7]

◆ kinParticles() [1/2]

◆ kinParticles() [2/2]

◆ mass()

◆ numParticles()

◆ operator=()

◆ p4()

◆ resetKinematicFit()

◆ setConstraint()

◆ setIndependentFit()

◆ setNotUpdated()

◆ splitKP()

◆ topDecayVertex()

◆ topParticle()

Member Data Documentation

◆ allParticles

◆ cKinP

◆ dMSig

◆ kCDMap

◆ kinMap

◆ kinTree

◆ massConst

◆ massSigma

◆ oldFit

◆ oldKPs

◆ oldMom

◆ tmpList

◆ totalMomentum