BPHDecayToTkpTknSymChargeBuilder Class Reference

#include <BPHDecayToTkpTknSymChargeBuilder.h>

Inheritance diagram for BPHDecayToTkpTknSymChargeBuilder:

Classes
class	Particle

Public Member Functions
	BPHDecayToTkpTknSymChargeBuilder (const edm::EventSetup &es, const std::string &daug1Name, double daug1Mass, double daug1Sigma, const std::string &daug2Name, double daug2Mass, double daug2Sigma, const BPHRecoBuilder::BPHGenericCollection *posCollection, const BPHRecoBuilder::BPHGenericCollection *negCollection, double expectedMass)
	BPHDecayToTkpTknSymChargeBuilder (const BPHDecayToTkpTknSymChargeBuilder &x)=delete
std::vector < BPHPlusMinusConstCandPtr >	build ()
	build candidates More...
double	getDzMax () const
double	getTrk1EtaMax () const
double	getTrk1PtMin () const
	get current cuts More...
double	getTrk2EtaMax () const
double	getTrk2PtMin () const
BPHDecayToTkpTknSymChargeBuilder &	operator= (const BPHDecayToTkpTknSymChargeBuilder &x)=delete
void	setDzMax (double dz)
void	setTrk1EtaMax (double eta)
void	setTrk1PtMin (double pt)
	set cuts More...
void	setTrk2EtaMax (double eta)
void	setTrk2PtMin (double pt)
	~BPHDecayToTkpTknSymChargeBuilder () override
Public Member Functions inherited from BPHDecayGenericBuilder
	BPHDecayGenericBuilder (const edm::EventSetup &es, BPHMassFitSelect *mfs=nullptr)
	BPHDecayGenericBuilder (const BPHDecayGenericBuilder &x)=delete
double	getMassFitMax () const
double	getMassFitMin () const
double	getMassMax () const
double	getMassMin () const
	get current cuts More...
double	getMinPDiff ()
double	getProbMin () const
BPHDecayGenericBuilder &	operator= (const BPHDecayGenericBuilder &x)=delete
void	setMassFitMax (double m)
void	setMassFitMin (double m)
void	setMassFitRange (double mMin, double mMax)
void	setMassMax (double m)
void	setMassMin (double m)
	set cuts More...
void	setMassRange (double mMin, double mMax)
void	setMinPDiff (double mpd)
	track min p difference More...
void	setProbMin (double p)
virtual	~BPHDecayGenericBuilder ()

Private Member Functions
void	addParticle (const BPHRecoBuilder::BPHGenericCollection *collection, int charge, std::vector< Particle * > &list)

Private Attributes
double	d1Mass
std::string	d1Name
double	d1Sigma
double	d2Mass
std::string	d2Name
double	d2Sigma
double	dzMax
double	eMass
double	eta1Max
double	eta2Max
const BPHRecoBuilder::BPHGenericCollection *	nCollection
const BPHRecoBuilder::BPHGenericCollection *	pCollection
double	pt1Min
double	pt2Min
std::vector < BPHPlusMinusConstCandPtr >	recList

Additional Inherited Members
Protected Attributes inherited from BPHDecayGenericBuilder
BPHChi2Select *	chi2Sel
const edm::EventSetup *	evSetup
BPHMassSelect *	massSel
BPHMassFitSelect *	mFitSel
double	minPDiff
bool	updated

Detailed Description

Description: Class to build a decay to an oppositely charged particle pair with different masses, choosing the best mass assignment on the reconstructed mass basis

Author

Paolo Ronchese INFN Padova

Definition at line 43 of file BPHDecayToTkpTknSymChargeBuilder.h.

Constructor & Destructor Documentation

BPHDecayToTkpTknSymChargeBuilder::BPHDecayToTkpTknSymChargeBuilder	(	const edm::EventSetup &	es,
		const std::string &	daug1Name,
		double	daug1Mass,
		double	daug1Sigma,
		const std::string &	daug2Name,
		double	daug2Mass,
		double	daug2Sigma,
		const BPHRecoBuilder::BPHGenericCollection *	posCollection,
		const BPHRecoBuilder::BPHGenericCollection *	negCollection,
		double	expectedMass
	)

Constructor

Definition at line 40 of file BPHDecayToTkpTknSymChargeBuilder.cc.

    : BPHDecayGenericBuilder(es),
      d1Name(daug1Name),
      d1Mass(daug1Mass),
      d1Sigma(daug1Sigma),
      d2Name(daug2Name),
      d2Mass(daug2Mass),
      d2Sigma(daug2Sigma),
      eMass(expectedMass),
      pCollection(posCollection),
      nCollection(negCollection),
      pt1Min(-1.0),
      pt2Min(-1.0),
      eta1Max(10.0),
      eta2Max(10.0),
      dzMax(1.0) {}

BPHDecayToTkpTknSymChargeBuilder::BPHDecayToTkpTknSymChargeBuilder ( const BPHDecayToTkpTknSymChargeBuilder &  x )

delete

BPHDecayToTkpTknSymChargeBuilder::~BPHDecayToTkpTknSymChargeBuilder ( )

override

Destructor

Definition at line 70 of file BPHDecayToTkpTknSymChargeBuilder.cc.

{}

Member Function Documentation

void BPHDecayToTkpTknSymChargeBuilder::addParticle ( const BPHRecoBuilder::BPHGenericCollection *  collection, int  charge, std::vector< Particle * > &  list  )

private

Definition at line 206 of file BPHDecayToTkpTknSymChargeBuilder.cc.

References reco::Candidate::charge(), d1Mass, d2Mass, eta1Max, eta2Max, BPHRecoBuilder::BPHGenericCollection::get(), BPHTrackReference::getTrack(), mps_fire::i, dqmiodumpmetadata::n, fireworks::p2, reco::Candidate::p4(), pt1Min, pt2Min, submitPVResolutionJobs::q, BPHRecoBuilder::BPHGenericCollection::size(), and mathSSE::sqrt().

Referenced by build().

                                                                            {
  int i;
  int n = collection->size();
  list.reserve(n);
  for (i = 0; i < n; ++i) {
    const reco::Candidate& cand = collection->get(i);
    int q = cand.charge();
    if ((charge > 0) && (q <= 0))
      continue;
    if ((charge < 0) && (q >= 0))
      continue;
    const reco::Candidate::LorentzVector p4 = cand.p4();
    const reco::Track* tk = BPHTrackReference::getTrack(cand, "cfhp");
    if (tk == nullptr)
      continue;
    double px = p4.px();
    double py = p4.py();
    double pz = p4.pz();
    double p2 = (px * px) + (py * py) + (pz * pz);
    double e1 = ((p4.pt() >= pt1Min) && (p4.eta() <= eta1Max) ? sqrt(p2 + (d1Mass * d1Mass)) : -1.0);
    double e2 = ((p4.pt() >= pt2Min) && (p4.eta() <= eta2Max) ? sqrt(p2 + (d2Mass * d2Mass)) : -1.0);
    if ((e1 > 0.0) || (e2 > 0.0))
      list.push_back(new Particle(&cand, tk, px, py, pz, e1, e2));
  }
  return;
}

vector< BPHPlusMinusConstCandPtr > BPHDecayToTkpTknSymChargeBuilder::build

(

)

build candidates

Operations

Definition at line 75 of file BPHDecayToTkpTknSymChargeBuilder.cc.

References BPHChi2Select::accept(), addParticle(), BPHDecayToTkpTknSymChargeBuilder::Particle::cand, BPHDecayGenericBuilder::chi2Sel, BPHGenericPtr< T >::create(), d1Mass, d1Name, d1Sigma, d2Mass, d2Name, d2Sigma, reco::TrackBase::dz(), dzMax, BPHDecayToTkpTknSymChargeBuilder::Particle::e1, BPHDecayToTkpTknSymChargeBuilder::Particle::e2, eMass, BPHDecayGenericBuilder::evSetup, BPHDecayGenericBuilder::getMassMax(), BPHDecayGenericBuilder::getMassMin(), ResonanceBuilder::mass, nCollection, nt, fireworks::p1, fireworks::p2, pCollection, DiDispStaMuonMonitor_cfi::pt, BPHDecayToTkpTknSymChargeBuilder::Particle::px, BPHDecayToTkpTknSymChargeBuilder::Particle::py, BPHDecayToTkpTknSymChargeBuilder::Particle::pz, recList, findQualityFiles::rr, BPHDecayToTkpTknSymChargeBuilder::Particle::track, and BPHDecayGenericBuilder::updated.

Referenced by BPHDecayToV0DiffMassBuilder::buildFromBPHGenericCollection(), and BPHWriteSpecificDecay::fill().

                                                                         {
  if (updated)
    return recList;

  recList.clear();

  // extract basic informations from input collections

  vector<Particle*> pList;
  vector<Particle*> nList;

  addParticle(pCollection, +1, pList);
  addParticle(nCollection, -1, nList);
  int iPos;
  int iNeg;
  int nPos = pList.size();
  int nNeg = nList.size();
  double massMin = getMassMin();
  double massMax = getMassMax();
  double mSqMin = massMin * massMin * 0.9;
  double mSqMax = massMax * massMax * 1.1;
  if (mSqMin < 0.0)
    mSqMin = 0.0;

  for (iPos = 0; iPos < nPos; ++iPos) {
    Particle* pc = pList[iPos];
    const reco::Track* pt = pc->track;
    double px = pc->px;
    double py = pc->py;
    double pz = pc->pz;
    double p1 = pc->e1;
    double p2 = pc->e2;
    for (iNeg = 0; iNeg < nNeg; ++iNeg) {
      Particle* nc = nList[iNeg];
      const reco::Track* nt = nc->track;
      if (fabs(nt->dz() - pt->dz()) > dzMax)
        continue;
      double nx = nc->px;
      double ny = nc->py;
      double nz = nc->pz;
      double n1 = nc->e1;
      double n2 = nc->e2;
      const float tx = px + nx;
      const float ty = py + ny;
      const float tz = pz + nz;
      const float ta = ((p1 > 0.0) && (n2 > 0.0) ? p1 + n2 : -1.0);
      const float tb = ((p2 > 0.0) && (n1 > 0.0) ? p2 + n1 : -1.0);
      float ma = (ta > 0 ? (ta * ta) - ((tx * tx) + (ty * ty) + (tz * tz)) : -1.0);
      float mb = (tb > 0 ? (tb * tb) - ((tx * tx) + (ty * ty) + (tz * tz)) : -1.0);
      //      if ( ma > 0.0 ) ma = sqrt( ma );
      //      else            ma = -1.0;
      //      if ( mb > 0.0 ) mb = sqrt( mb );
      //      else            mb = -1.0;
      if (((ma < mSqMin) || (ma > mSqMax)) && ((mb < mSqMin) || (mb > mSqMax)))
        continue;
      BPHPlusMinusCandidatePtr rc(nullptr);
      float rcMass = -1.0;
      if (ma > 0.0) {
        rc = BPHPlusMinusCandidateWrap::create(evSetup);
        rc->add(d1Name, pc->cand, d1Mass, d1Sigma);
        rc->add(d2Name, nc->cand, d2Mass, d2Sigma);
        rcMass = rc->composite().mass();
      }
      BPHPlusMinusCandidatePtr rb(nullptr);
      float rbMass = -1.0;
      if (mb > 0.0) {
        rb = BPHPlusMinusCandidateWrap::create(evSetup);
        rb->add(d1Name, nc->cand, d1Mass, d1Sigma);
        rb->add(d2Name, pc->cand, d2Mass, d2Sigma);
        rbMass = rb->composite().mass();
      }
      BPHPlusMinusCandidatePtr* rp(nullptr);
      double mass = -1.0;
      if ((rc.get() != nullptr) && ((rb.get() == nullptr) || (fabs(rcMass - eMass) < fabs(rbMass - eMass)))) {
        mass = rcMass;
        rp = &rc;
      } else {
        mass = rbMass;
        rp = &rb;
      }
      BPHPlusMinusCandidate* rr = rp->get();
      if (mass < massMin)
        continue;
      if (mass > massMax)
        continue;
      if (!chi2Sel->accept(*rr))
        continue;
      recList.push_back(*rp);
    }
  }

  for (iPos = 0; iPos < nPos; ++iPos)
    delete pList[iPos];
  for (iNeg = 0; iNeg < nNeg; ++iNeg)
    delete nList[iNeg];

  updated = true;
  return recList;
}

double BPHDecayToTkpTknSymChargeBuilder::getDzMax ( ) const

inline

Definition at line 83 of file BPHDecayToTkpTknSymChargeBuilder.h.

References dzMax.

{ return dzMax; }

double BPHDecayToTkpTknSymChargeBuilder::getTrk1EtaMax ( ) const

inline

Definition at line 81 of file BPHDecayToTkpTknSymChargeBuilder.h.

References eta1Max.

Referenced by BPHKx0ToKPiBuilder::getPiEtaMax().

{ return eta1Max; }

double BPHDecayToTkpTknSymChargeBuilder::getTrk1PtMin ( ) const

inline

get current cuts

Definition at line 79 of file BPHDecayToTkpTknSymChargeBuilder.h.

References pt1Min.

Referenced by BPHKx0ToKPiBuilder::getPiPtMin().

{ return pt1Min; }

double BPHDecayToTkpTknSymChargeBuilder::getTrk2EtaMax ( ) const

inline

Definition at line 82 of file BPHDecayToTkpTknSymChargeBuilder.h.

References eta2Max.

Referenced by BPHKx0ToKPiBuilder::getKEtaMax().

{ return eta2Max; }

double BPHDecayToTkpTknSymChargeBuilder::getTrk2PtMin ( ) const

inline

Definition at line 80 of file BPHDecayToTkpTknSymChargeBuilder.h.

References pt2Min.

Referenced by BPHKx0ToKPiBuilder::getKPtMin().

{ return pt2Min; }

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

delete

void BPHDecayToTkpTknSymChargeBuilder::setDzMax

(

double

dz

)

Definition at line 200 of file BPHDecayToTkpTknSymChargeBuilder.cc.

References PVValHelper::dz, dzMax, and BPHDecayGenericBuilder::updated.

                                                         {
  updated = false;
  dzMax = dz;
  return;
}

void BPHDecayToTkpTknSymChargeBuilder::setTrk1EtaMax

(

double

eta

)

Definition at line 188 of file BPHDecayToTkpTknSymChargeBuilder.cc.

References PVValHelper::eta, eta1Max, and BPHDecayGenericBuilder::updated.

Referenced by BPHKx0ToKPiBuilder::BPHKx0ToKPiBuilder(), BPHDecayToV0DiffMassBuilder::buildFromBPHGenericCollection(), BPHKx0ToKPiBuilder::setEtaMax(), and BPHKx0ToKPiBuilder::setPiEtaMax().

                                                               {
  updated = false;
  eta1Max = eta;
  return;
}

void BPHDecayToTkpTknSymChargeBuilder::setTrk1PtMin

(

double

pt

)

set cuts

Definition at line 176 of file BPHDecayToTkpTknSymChargeBuilder.cc.

References DiDispStaMuonMonitor_cfi::pt, pt1Min, and BPHDecayGenericBuilder::updated.

Referenced by BPHKx0ToKPiBuilder::BPHKx0ToKPiBuilder(), BPHDecayToV0DiffMassBuilder::buildFromBPHGenericCollection(), and BPHKx0ToKPiBuilder::setPtMin().

                                                             {
  updated = false;
  pt1Min = pt;
  return;
}

void BPHDecayToTkpTknSymChargeBuilder::setTrk2EtaMax

(

double

eta

)

Definition at line 194 of file BPHDecayToTkpTknSymChargeBuilder.cc.

References PVValHelper::eta, eta2Max, and BPHDecayGenericBuilder::updated.

Referenced by BPHKx0ToKPiBuilder::BPHKx0ToKPiBuilder(), BPHDecayToV0DiffMassBuilder::buildFromBPHGenericCollection(), BPHKx0ToKPiBuilder::setEtaMax(), and BPHKx0ToKPiBuilder::setKEtaMax().

                                                               {
  updated = false;
  eta2Max = eta;
  return;
}

void BPHDecayToTkpTknSymChargeBuilder::setTrk2PtMin

(

double

pt

)

Definition at line 182 of file BPHDecayToTkpTknSymChargeBuilder.cc.

References DiDispStaMuonMonitor_cfi::pt, pt2Min, and BPHDecayGenericBuilder::updated.

Referenced by BPHKx0ToKPiBuilder::BPHKx0ToKPiBuilder(), BPHDecayToV0DiffMassBuilder::buildFromBPHGenericCollection(), BPHKx0ToKPiBuilder::setKPtMin(), and BPHKx0ToKPiBuilder::setPtMin().

                                                             {
  updated = false;
  pt2Min = pt;
  return;
}

Member Data Documentation

double BPHDecayToTkpTknSymChargeBuilder::d1Mass

private

Definition at line 87 of file BPHDecayToTkpTknSymChargeBuilder.h.

Referenced by addParticle(), and build().

std::string BPHDecayToTkpTknSymChargeBuilder::d1Name

private

Definition at line 86 of file BPHDecayToTkpTknSymChargeBuilder.h.

Referenced by build().

double BPHDecayToTkpTknSymChargeBuilder::d1Sigma

private

Definition at line 88 of file BPHDecayToTkpTknSymChargeBuilder.h.

Referenced by build().

double BPHDecayToTkpTknSymChargeBuilder::d2Mass

private

Definition at line 90 of file BPHDecayToTkpTknSymChargeBuilder.h.

Referenced by addParticle(), and build().

std::string BPHDecayToTkpTknSymChargeBuilder::d2Name

private

Definition at line 89 of file BPHDecayToTkpTknSymChargeBuilder.h.

Referenced by build().

double BPHDecayToTkpTknSymChargeBuilder::d2Sigma

private

Definition at line 91 of file BPHDecayToTkpTknSymChargeBuilder.h.

Referenced by build().

double BPHDecayToTkpTknSymChargeBuilder::dzMax

private

Definition at line 101 of file BPHDecayToTkpTknSymChargeBuilder.h.

Referenced by build(), getDzMax(), and setDzMax().

double BPHDecayToTkpTknSymChargeBuilder::eMass

private

Definition at line 92 of file BPHDecayToTkpTknSymChargeBuilder.h.

Referenced by build().

double BPHDecayToTkpTknSymChargeBuilder::eta1Max

private

Definition at line 99 of file BPHDecayToTkpTknSymChargeBuilder.h.

Referenced by addParticle(), getTrk1EtaMax(), and setTrk1EtaMax().

double BPHDecayToTkpTknSymChargeBuilder::eta2Max

private

Definition at line 100 of file BPHDecayToTkpTknSymChargeBuilder.h.

Referenced by addParticle(), getTrk2EtaMax(), and setTrk2EtaMax().

const BPHRecoBuilder::BPHGenericCollection* BPHDecayToTkpTknSymChargeBuilder::nCollection

private

Definition at line 95 of file BPHDecayToTkpTknSymChargeBuilder.h.

Referenced by build().

const BPHRecoBuilder::BPHGenericCollection* BPHDecayToTkpTknSymChargeBuilder::pCollection

private

Definition at line 94 of file BPHDecayToTkpTknSymChargeBuilder.h.

Referenced by build().

double BPHDecayToTkpTknSymChargeBuilder::pt1Min

private

Definition at line 97 of file BPHDecayToTkpTknSymChargeBuilder.h.

Referenced by addParticle(), getTrk1PtMin(), and setTrk1PtMin().

double BPHDecayToTkpTknSymChargeBuilder::pt2Min

private

Definition at line 98 of file BPHDecayToTkpTknSymChargeBuilder.h.

Referenced by addParticle(), getTrk2PtMin(), and setTrk2PtMin().

std::vector<BPHPlusMinusConstCandPtr> BPHDecayToTkpTknSymChargeBuilder::recList

private

Definition at line 103 of file BPHDecayToTkpTknSymChargeBuilder.h.

Referenced by build().