#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) 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) const

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

virtual const RefCountedKinematicTree &	kinematicTree (const std::string &name, double mass, double sigma) 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

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

const edm::EventSetup *	getEventSetup () const
	retrieve EventSetup 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

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

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 BPHRecoConstCandPtr &comp)
	add a previously reconstructed particle giving it a name More...

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 reco::Candidate *daug, double mass=-1.0, double sigma=-1.0)

	BPHKinematicFit ()

	BPHKinematicFit (const BPHKinematicFit *ptr)

void	setNotUpdated () const override

Protected Member Functions inherited from BPHDecayVertex
virtual void	addV (const std::string &name, const BPHRecoConstCandPtr &comp)
	add a previously reconstructed particle giving it a name More...

virtual void	addV (const std::string &name, const reco::Candidate *daug, const std::string &searchList, double mass)

	BPHDecayVertex (const BPHDecayVertex ptr, const edm::EventSetup es)

	BPHDecayVertex (const edm::EventSetup *es)

Protected Member Functions inherited from BPHDecayMomentum
virtual void	addP (const std::string &name, const BPHRecoConstCandPtr &comp)
	add a previously reconstructed particle giving it a name More...

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

	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	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

virtual void	getParticles (const std::string &moth, const std::vector< std::string > &daug, std::vector< RefCountedKinematicParticle > &kl, std::set< RefCountedKinematicParticle > &ks) 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 ( )

override

Destructor

Definition at line 85 of file BPHKinematicFit.cc.

85 {}

◆ BPHKinematicFit() [2/3]

BPHKinematicFit::BPHKinematicFit ( )

protected

Definition at line 39 of file BPHKinematicFit.cc.

     : 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 48 of file BPHKinematicFit.cc.

     : 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];
   }
 }

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().

Member Function Documentation

◆ addK() [1/3]

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

protectedvirtual

add a previously reconstructed particle giving it a name

◆ 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

add a simple particle and specify a criterion to search for the associated track

◆ addK() [3/3]

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

protectedvirtual

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

Referenced by BPHRecoCandidate::add().

◆ 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.

                                                                                                         {
   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) && (sigma > 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;
 }

References addParticles(), 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, topParticle(), and groupFilesInBlocks::tt.

Referenced by addParticles(), and buildParticles().

◆ buildParticles()

void BPHKinematicFit::buildParticles ( ) const

privatevirtual

Definition at line 361 of file BPHKinematicFit.cc.

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

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

Referenced by kinParticles().

◆ constrMass()

double BPHKinematicFit::constrMass ( ) const

retrieve the constraint

Definition at line 99 of file BPHKinematicFit.cc.

99 { return massConst; }

References massConst.

Referenced by BPHRecoCandidate::fill().

◆ constrSigma()

double BPHKinematicFit::constrSigma ( ) const

Definition at line 101 of file BPHKinematicFit.cc.

101 { return massSigma; }

References massSigma.

Referenced by BPHRecoCandidate::fill().

◆ currentDecayVertex()

const RefCountedKinematicVertex BPHKinematicFit::currentDecayVertex ( ) const

virtual

Definition at line 277 of file BPHKinematicFit.cc.

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

References isEmpty(), and kinTree.

◆ currentParticle()

const RefCountedKinematicParticle BPHKinematicFit::currentParticle ( ) const

virtual

get current particle

Definition at line 271 of file BPHKinematicFit.cc.

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

References isEmpty(), and kinTree.

◆ fitMomentum()

void BPHKinematicFit::fitMomentum ( ) const

privatevirtual

Definition at line 538 of file BPHKinematicFit.cc.

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

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().

◆ getIndependentFit()

bool BPHKinematicFit::getIndependentFit ( const std::string & name ) const

retrieve independent fit flag

Definition at line 321 of file BPHKinematicFit.cc.

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

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

Referenced by BPHRecoCandidate::fill().

◆ getMassSigma()

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

retrieve particle mass sigma

Definition at line 315 of file BPHKinematicFit.cc.

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

References dMSig.

Referenced by BPHRecoCandidate::fill().

◆ getParticles() [1/2]

virtual void BPHKinematicFit::getParticles	(	const std::string &	moth,
		const std::string &	daug,
		std::vector< RefCountedKinematicParticle > &	kl,
		std::set< RefCountedKinematicParticle > &	ks
	)		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

privatevirtual

◆ insertParticle()

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

staticprivate

Definition at line 477 of file BPHKinematicFit.cc.

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

References kp.

◆ isEmpty()

bool BPHKinematicFit::isEmpty ( void ) const

virtual

get fit status

Definition at line 256 of file BPHKinematicFit.cc.

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

References kinematicTree(), and kinTree.

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

◆ isValidFit()

bool BPHKinematicFit::isValidFit ( ) const

virtual

Definition at line 263 of file BPHKinematicFit.cc.

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

References topParticle().

Referenced by addParticles(), and fitMomentum().

◆ kinematicTree() [1/7]

const RefCountedKinematicTree & BPHKinematicFit::kinematicTree ( ) const

virtual

perform the kinematic fit and get the result

Definition at line 160 of file BPHKinematicFit.cc.

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

References kinTree, massConst, massSigma, and oldFit.

Referenced by isEmpty().

◆ kinematicTree() [2/7]

virtual const RefCountedKinematicTree& BPHKinematicFit::kinematicTree ( const std::string & name ) 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,
		double	mass,
		double	sigma
	)		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 525 of file BPHKinematicFit.cc.

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

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

◆ kinParticles() [1/2]

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

virtual

get kinematic particles

Definition at line 125 of file BPHKinematicFit.cc.

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

References allParticles, buildParticles(), and oldKPs.

Referenced by splitKP().

◆ 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 297 of file BPHKinematicFit.cc.

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

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

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

◆ numParticles()

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

privatevirtual

Definition at line 461 of file BPHKinematicFit.cc.

                                                                             {
   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 (cKinP.at(comp).flag)
       ++n;
     else
       n += numParticles(comp);
   }
   return n;
 }

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

Referenced by splitKP().

◆ 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 308 of file BPHKinematicFit.cc.

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

References fitMomentum(), oldMom, and totalMomentum.

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

◆ resetKinematicFit()

void BPHKinematicFit::resetKinematicFit ( ) const

virtual

reset the kinematic fit

Definition at line 250 of file BPHKinematicFit.cc.

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

References oldFit, oldKPs, and oldMom.

Referenced by BPHDecayToResFlyingBuilder::build(), and setNotUpdated().

◆ setConstraint()

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

apply a mass constraint

Operations

Definition at line 91 of file BPHKinematicFit.cc.

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

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

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

◆ 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 104 of file BPHKinematicFit.cc.

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

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

Referenced by BPHDecayToResFlyingBuilder::build(), and BPHRecoCandidate::fill().

◆ setNotUpdated()

void BPHKinematicFit::setNotUpdated ( ) const

overrideprotectedvirtual

Reimplemented from BPHDecayVertex.

Reimplemented in BPHPlusMinusCandidate.

Definition at line 354 of file BPHKinematicFit.cc.

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

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

Referenced by BPHPlusMinusCandidate::setNotUpdated().

◆ splitKP()

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

privatevirtual

Definition at line 487 of file BPHKinematicFit.cc.

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

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

◆ topDecayVertex()

const RefCountedKinematicVertex BPHKinematicFit::topDecayVertex ( ) const

virtual

Definition at line 290 of file BPHKinematicFit.cc.

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

References isEmpty(), and kinTree.

Referenced by BPHDecayToResFlyingBuilder::build().

◆ topParticle()

const RefCountedKinematicParticle BPHKinematicFit::topParticle ( ) const

virtual

get top particle

Definition at line 284 of file BPHKinematicFit.cc.

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

References isEmpty(), and kinTree.

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

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(), numParticles(), 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