BPHKinematicFit Class Reference

#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, 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
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 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)
	add a previously reconstructed particle giving it a name More...
	BPHKinematicFit ()
	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)
	add a previously reconstructed particle giving it a name More...
	BPHDecayVertex (const edm::EventSetup *es)
	BPHDecayVertex (const BPHDecayVertex *ptr, const edm::EventSetup *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)
	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	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::BPHKinematicFit ( const BPHKinematicFit &  x )

delete

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

BPHKinematicFit::~BPHKinematicFit ( )

override

Destructor

Definition at line 85 of file BPHKinematicFit.cc.

{}

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::BPHKinematicFit ( const BPHKinematicFit *  ptr )

protected

Definition at line 48 of file BPHKinematicFit.cc.

References c, BPHDecayMomentum::Component::cand, cKinP, BPHDecayMomentum::cMap, BPHDecayMomentum::componentList(), BPHDecayMomentum::daughComp(), BPHDecayMomentum::daughters(), BPHDecayMomentum::dMap, dMSig, mps_fire::i, dqmiolumiharvest::j, visualization-live-secondInstance_cfg::m, BPHDecayMomentum::Component::msig, 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

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

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

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

protectedvirtual

add a previously reconstructed particle giving it a name

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(), BPHDecayMomentum::daughComp(), BPHDecayMomentum::daughters(), dMSig, BPHKinematicFit::FlyingParticle::flag, personalPlayback::fp, BPHDecayVertex::getTransientTrack(), isEmpty(), isValidFit(), visualization-live-secondInstance_cfg::m, reco::Candidate::mass(), BPHKinematicFit::FlyingParticle::mass, dqmiodumpmetadata::n, KinematicParticleFactoryFromTransientTrack::particle(), setConstraint(), BPHKinematicFit::FlyingParticle::sigma, tmpList, topParticle(), and groupFilesInBlocks::tt.

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

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

double BPHKinematicFit::constrMass

(

)

const

retrieve the constraint

Definition at line 99 of file BPHKinematicFit.cc.

References massConst.

Referenced by BPHRecoCandidate::fill().

{ return massConst; }

double BPHKinematicFit::constrSigma

(

)

const

Definition at line 101 of file BPHKinematicFit.cc.

References massSigma.

Referenced by BPHRecoCandidate::fill().

{ return massSigma; }

const RefCountedKinematicVertex BPHKinematicFit::currentDecayVertex ( ) const

virtual

Definition at line 277 of file BPHKinematicFit.cc.

References isEmpty(), and kinTree.

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

const RefCountedKinematicParticle BPHKinematicFit::currentParticle ( ) const

virtual

get current particle

Definition at line 271 of file BPHKinematicFit.cc.

References isEmpty(), and kinTree.

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

void BPHKinematicFit::fitMomentum ( ) const

privatevirtual

Definition at line 538 of file BPHKinematicFit.cc.

References BPHDecayMomentum::daughComp(), BPHDecayMomentum::daughters(), alignCSCRings::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;
}

bool BPHKinematicFit::getIndependentFit

(

const std::string &

name

)

const

retrieve independent fit flag

Definition at line 321 of file BPHKinematicFit.cc.

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

Referenced by BPHRecoCandidate::fill().

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

double BPHKinematicFit::getMassSigma

(

const reco::Candidate *

cand

)

const

retrieve particle mass sigma

Definition at line 315 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);
}

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

privatevirtual

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

privatevirtual

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

bool BPHKinematicFit::isEmpty ( void  ) const

virtual

get fit status

Definition at line 256 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();
}

bool BPHKinematicFit::isValidFit ( ) const

virtual

Definition at line 263 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();
}

const RefCountedKinematicTree & BPHKinematicFit::kinematicTree ( ) const

virtual

perform the kinematic fit and get the result

Definition at line 160 of file BPHKinematicFit.cc.

References kinTree, massConst, massSigma, and oldFit.

Referenced by isEmpty().

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

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

virtual

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

virtual

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

virtual

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

virtual

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

virtual

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

privatevirtual

Definition at line 525 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;
}

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

virtual

get kinematic particles

Definition at line 125 of file BPHKinematicFit.cc.

References allParticles, buildParticles(), and oldKPs.

Referenced by splitKP().

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

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

virtual

ParticleMass BPHKinematicFit::mass ( ) const

virtual

Definition at line 297 of file BPHKinematicFit.cc.

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

Referenced by Particle.Particle::__str__(), DiObject.DiMuon::__str__(), BPHRecoCandidate::add(), BPHPlusMinusCandidate::build(), BPHRecoCandidate::build(), BPHDecayToV0DiffMassBuilder::buildCandidate(), 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;
}

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

privatevirtual

Definition at line 461 of file BPHKinematicFit.cc.

References cKinP, BPHDecayMomentum::compNames(), BPHDecayMomentum::daughters(), BPHDecayMomentum::getComp(), 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 (cKinP.at(comp).flag)
      ++n;
    else
      n += numParticles(comp);
  }
  return n;
}

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

delete

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

virtual

compute total momentum after the fit

Definition at line 308 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;
}

void BPHKinematicFit::resetKinematicFit ( ) const

virtual

reset the kinematic fit

Definition at line 250 of file BPHKinematicFit.cc.

References oldFit, oldKPs, and oldMom.

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

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

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

apply a mass constraint

Operations

Definition at line 91 of file BPHKinematicFit.cc.

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

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

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

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.

References cKinP, AlCaHLTBitMon_QueryRunRegistry::comp, BPHKinematicFit::FlyingParticle::flag, personalPlayback::fp, BPHDecayMomentum::getComp(), mass(), BPHKinematicFit::FlyingParticle::mass, mergeVDriftHistosByStation::name, oldFit, oldKPs, oldMom, and BPHKinematicFit::FlyingParticle::sigma.

Referenced by BPHDecayToResFlyingBuilder::build(), and 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;
}

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

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.

References allParticles, cKinP, AlCaHLTBitMon_QueryRunRegistry::comp, BPHDecayMomentum::getComp(), kinParticles(), kinTree, mergeVDriftHistosByStation::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;
}

const RefCountedKinematicVertex BPHKinematicFit::topDecayVertex ( ) const

virtual

Definition at line 290 of file BPHKinematicFit.cc.

References isEmpty(), and kinTree.

Referenced by BPHDecayToResFlyingBuilder::build().

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

const RefCountedKinematicParticle BPHKinematicFit::topParticle ( ) const

virtual

get top particle

Definition at line 284 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

std::vector<RefCountedKinematicParticle> BPHKinematicFit::allParticles

mutableprivate

Definition at line 148 of file BPHKinematicFit.h.

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

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

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

private

Definition at line 129 of file BPHKinematicFit.h.

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

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

mutableprivate

Definition at line 147 of file BPHKinematicFit.h.

Referenced by buildParticles().

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

mutableprivate

Definition at line 146 of file BPHKinematicFit.h.

Referenced by buildParticles().

RefCountedKinematicTree BPHKinematicFit::kinTree

mutableprivate

Definition at line 149 of file BPHKinematicFit.h.

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

double BPHKinematicFit::massConst

private

Definition at line 125 of file BPHKinematicFit.h.

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

double BPHKinematicFit::massSigma

private

Definition at line 126 of file BPHKinematicFit.h.

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

bool BPHKinematicFit::oldFit

mutableprivate

Definition at line 144 of file BPHKinematicFit.h.

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

bool BPHKinematicFit::oldKPs

mutableprivate

Definition at line 143 of file BPHKinematicFit.h.

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

bool BPHKinematicFit::oldMom

mutableprivate

Definition at line 145 of file BPHKinematicFit.h.

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

std::vector<BPHRecoConstCandPtr> BPHKinematicFit::tmpList

mutableprivate

Definition at line 140 of file BPHKinematicFit.h.

Referenced by addParticles().

math::XYZTLorentzVector BPHKinematicFit::totalMomentum

mutableprivate

Definition at line 150 of file BPHKinematicFit.h.

Referenced by fitMomentum(), and p4().