BPHX3872ToJPsiPiPiBuilder Class Reference

#include <BPHX3872ToJPsiPiPiBuilder.h>

Public Member Functions
	BPHX3872ToJPsiPiPiBuilder (const edm::EventSetup &es, const std::vector< BPHPlusMinusConstCandPtr > &jpsiCollection, const BPHRecoBuilder::BPHGenericCollection *posCollection, const BPHRecoBuilder::BPHGenericCollection *negCollection)
	BPHX3872ToJPsiPiPiBuilder (const BPHX3872ToJPsiPiPiBuilder &x)=delete
std::vector< BPHRecoConstCandPtr >	build ()
	build X3872 candidates More...
bool	getConstr () const
double	getJPsiMassMax () const
double	getJPsiMassMin () const
	get current cuts More...
double	getMassFitMax () const
double	getMassFitMin () const
double	getMassMax () const
double	getMassMin () const
double	getPiEtaMax () const
double	getPiPtMin () const
	get current cuts More...
double	getProbMin () const
BPHX3872ToJPsiPiPiBuilder &	operator= (const BPHX3872ToJPsiPiPiBuilder &x)=delete
void	setConstr (bool flag)
void	setJPsiMassMax (double m)
void	setJPsiMassMin (double m)
	set cuts More...
void	setMassFitMax (double m)
void	setMassFitMin (double m)
void	setMassMax (double m)
void	setMassMin (double m)
void	setPiEtaMax (double eta)
void	setPiPtMin (double pt)
	set cuts More...
void	setProbMin (double p)
virtual	~BPHX3872ToJPsiPiPiBuilder ()

Private Attributes
BPHChi2Select *	chi2Sel
double	etaMax
const edm::EventSetup *	evSetup
const std::vector < BPHPlusMinusConstCandPtr > *	jCollection
std::string	jPsiName
BPHMassSelect *	jpsiSel
bool	massConstr
BPHMassSelect *	massSel
BPHMassFitSelect *	mFitSel
float	minPDiff
const BPHRecoBuilder::BPHGenericCollection *	nCollection
const BPHRecoBuilder::BPHGenericCollection *	pCollection
std::string	pionNegName
std::string	pionPosName
double	ptMin
bool	updated
std::vector< BPHRecoConstCandPtr >	x3872List

Detailed Description

Description: Class to build X3872 to JPsi pi+ pi- candidates

Author

Paolo Ronchese INFN Padova

Definition at line 42 of file BPHX3872ToJPsiPiPiBuilder.h.

Constructor & Destructor Documentation

BPHX3872ToJPsiPiPiBuilder::BPHX3872ToJPsiPiPiBuilder	(	const edm::EventSetup &	es,
		const std::vector< BPHPlusMinusConstCandPtr > &	jpsiCollection,
		const BPHRecoBuilder::BPHGenericCollection *	posCollection,
		const BPHRecoBuilder::BPHGenericCollection *	negCollection
	)

Constructor

Definition at line 36 of file BPHX3872ToJPsiPiPiBuilder.cc.

References chi2Sel, etaMax, BPHParticleMasses::jPsiMass, BPHParticleMasses::jPsiMWidth, jPsiName, jpsiSel, massConstr, massSel, mFitSel, minPDiff, ptMin, and updated.

    : jPsiName("JPsi"),
      pionPosName("PionPos"),
      pionNegName("PionNeg"),
      evSetup(&es),
      jCollection(&jpsiCollection),
      pCollection(posCollection),
      nCollection(negCollection) {
  jpsiSel = new BPHMassSelect(2.80, 3.40);
  ptMin = 1.0;
  etaMax = 10.0;
  massSel = new BPHMassSelect(3.00, 4.50);
  chi2Sel = new BPHChi2Select(0.02);
  mFitSel = new BPHMassFitSelect(jPsiName, BPHParticleMasses::jPsiMass, BPHParticleMasses::jPsiMWidth, 3.50, 4.20);
  massConstr = true;
  minPDiff = 1.0e-4;
  updated = false;
}

BPHX3872ToJPsiPiPiBuilder::BPHX3872ToJPsiPiPiBuilder ( const BPHX3872ToJPsiPiPiBuilder &  x )

delete

BPHX3872ToJPsiPiPiBuilder::~BPHX3872ToJPsiPiPiBuilder ( )

virtual

Destructor

Definition at line 61 of file BPHX3872ToJPsiPiPiBuilder.cc.

References chi2Sel, jpsiSel, massSel, and mFitSel.

                                                      {
  delete jpsiSel;
  delete massSel;
  delete chi2Sel;
  delete mFitSel;
}

Member Function Documentation

vector< BPHRecoConstCandPtr > BPHX3872ToJPsiPiPiBuilder::build

(

)

build X3872 candidates

Operations

Definition at line 71 of file BPHX3872ToJPsiPiPiBuilder.cc.

References BPHMassSelect::accept(), BPHChi2Select::accept(), BPHMassFitSelect::accept(), BPHRecoCandidate::add(), c, reco::Candidate::charge(), chi2Sel, reco::TrackBase::dz(), etaMax, evSetup, BPHRecoBuilder::BPHGenericCollection::get(), getMassMax(), BPHTrackReference::getTrack(), jCollection, BPHParticleMasses::jPsiMass, jPsiName, jpsiSel, ResonanceBuilder::mass, massSel, mFitSel, nCollection, compareTotals::ne, nt, AlCaHLTBitMon_ParallelJobs::p, fireworks::p2, reco::Candidate::p4(), pCollection, BPHParticleMasses::pionMass, BPHParticleMasses::pionMSigma, pionNegName, pionPosName, createTree::pp, DiDispStaMuonMonitor_cfi::pt, ptMin, BPHRecoBuilder::BPHGenericCollection::size(), mathSSE::sqrt(), HLT_FULL_cff::track, updated, x, x3872List, y, and z.

Referenced by BPHWriteSpecificDecay::fill().

                                                             {
  if (updated)
    return x3872List;

  x3872List.clear();

  // extract basic informations from input collections

  class Particle {
  public:
    Particle(const reco::Candidate* c, const reco::Track* tk, double x, double y, double z, double p)
        : cand(c), track(tk), px(x), py(y), pz(z), ePion(p) {}
    const reco::Candidate* cand;
    const reco::Track* track;
    double px;
    double py;
    double pz;
    double ePion;
  };

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

  int nPos = pCollection->size();
  int nNeg = nCollection->size();

  pList.reserve(nPos);
  nList.reserve(nNeg);

  // filter input collections

  int iPos;
  int iNeg;

  for (iPos = 0; iPos < nPos; ++iPos) {
    const reco::Candidate& cand = pCollection->get(iPos);
    if (cand.charge() != +1)
      continue;
    const reco::Candidate::LorentzVector p4 = cand.p4();
    if (p4.pt() < ptMin)
      continue;
    if (p4.eta() > etaMax)
      continue;
    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);
    pList.push_back(
        Particle(&cand, tk, px, py, pz, sqrt(p2 + (BPHParticleMasses::pionMass * BPHParticleMasses::pionMass))));
  }

  for (iNeg = 0; iNeg < nNeg; ++iNeg) {
    const reco::Candidate& cand = nCollection->get(iNeg);
    if (cand.charge() != -1)
      continue;
    const reco::Candidate::LorentzVector p4 = cand.p4();
    if (p4.pt() < ptMin)
      continue;
    if (p4.eta() > etaMax)
      continue;
    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);
    nList.push_back(
        Particle(&cand, tk, px, py, pz, sqrt(p2 + (BPHParticleMasses::pionMass * BPHParticleMasses::pionMass))));
  }

  // filter basic candidates

  nPos = pList.size();
  nNeg = nList.size();

  double mMax = getMassMax() - (0.8 * BPHParticleMasses::jPsiMass);

  struct PionPair {
    const reco::Candidate* posPion;
    const reco::Candidate* negPion;
  };
  vector<const PionPair*> pionPairs;
  pionPairs.reserve(nPos * nNeg);

  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 pe = pc.ePion;
    for (iNeg = 0; iNeg < nNeg; ++iNeg) {
      Particle& nc = nList[iNeg];
      const reco::Track* nt = nc.track;
      if (fabs(nt->dz() - pt->dz()) > 1.0)
        continue;
      double nx = nc.px;
      double ny = nc.py;
      double nz = nc.pz;
      double ne = nc.ePion;
      const float tx = px + nx;
      const float ty = py + ny;
      const float tz = pz + nz;
      const float te = pe + ne;
      float mass = (te * te) - ((tx * tx) + (ty * ty) + (tz * tz));
      if (mass > mMax)
        continue;
      PionPair* pp = new PionPair;
      pp->posPion = pc.cand;
      pp->negPion = nc.cand;
      pionPairs.push_back(pp);
    }
  }

  vector<BPHPlusMinusConstCandPtr> jPsi;
  int nJPsi = jCollection->size();
  int iJPsi;
  jPsi.reserve(nJPsi);
  for (iJPsi = 0; iJPsi < nJPsi; ++iJPsi) {
    const BPHPlusMinusConstCandPtr& jpCand = jCollection->at(iJPsi);
    if (jpsiSel->accept(*jpCand))
      jPsi.push_back(jpCand);
  }
  nJPsi = jPsi.size();

  int nPair = pionPairs.size();
  int iPair;
  for (iPair = 0; iPair < nPair; ++iPair) {
    const PionPair* pp = pionPairs[iPair];
    for (iJPsi = 0; iJPsi < nJPsi; ++iJPsi) {
      BPHRecoCandidate* x3872 = new BPHRecoCandidate(evSetup);
      BPHRecoCandidatePtr xPtr(x3872);
      x3872->add(jPsiName, jPsi[iJPsi]);
      x3872->add(pionPosName, pp->posPion, BPHParticleMasses::pionMass, BPHParticleMasses::pionMSigma);
      x3872->add(pionNegName, pp->negPion, BPHParticleMasses::pionMass, BPHParticleMasses::pionMSigma);
      if (!massSel->accept(*x3872))
        continue;
      if ((chi2Sel != nullptr) && !chi2Sel->accept(*x3872))
        continue;
      if (!mFitSel->accept(*x3872))
        continue;
      x3872List.push_back(xPtr);
    }
    delete pp;
  }

  updated = true;
  return x3872List;
}

bool BPHX3872ToJPsiPiPiBuilder::getConstr

(

)

const

Definition at line 306 of file BPHX3872ToJPsiPiPiBuilder.cc.

References massConstr.

{ return massConstr; }

double BPHX3872ToJPsiPiPiBuilder::getJPsiMassMax

(

)

const

Definition at line 290 of file BPHX3872ToJPsiPiPiBuilder.cc.

References BPHMassCuts::getMassMax(), and jpsiSel.

{ return jpsiSel->getMassMax(); }

double BPHX3872ToJPsiPiPiBuilder::getJPsiMassMin

(

)

const

get current cuts

Definition at line 288 of file BPHX3872ToJPsiPiPiBuilder.cc.

References BPHMassCuts::getMassMin(), and jpsiSel.

{ return jpsiSel->getMassMin(); }

double BPHX3872ToJPsiPiPiBuilder::getMassFitMax

(

)

const

Definition at line 304 of file BPHX3872ToJPsiPiPiBuilder.cc.

References BPHMassCuts::getMassMax(), and mFitSel.

{ return mFitSel->getMassMax(); }

double BPHX3872ToJPsiPiPiBuilder::getMassFitMin

(

)

const

Definition at line 302 of file BPHX3872ToJPsiPiPiBuilder.cc.

References BPHMassCuts::getMassMin(), and mFitSel.

{ return mFitSel->getMassMin(); }

double BPHX3872ToJPsiPiPiBuilder::getMassMax

(

)

const

Definition at line 298 of file BPHX3872ToJPsiPiPiBuilder.cc.

References BPHMassCuts::getMassMax(), and massSel.

Referenced by build().

{ return massSel->getMassMax(); }

double BPHX3872ToJPsiPiPiBuilder::getMassMin

(

)

const

Definition at line 296 of file BPHX3872ToJPsiPiPiBuilder.cc.

References BPHMassCuts::getMassMin(), and massSel.

{ return massSel->getMassMin(); }

double BPHX3872ToJPsiPiPiBuilder::getPiEtaMax

(

)

const

Definition at line 294 of file BPHX3872ToJPsiPiPiBuilder.cc.

References etaMax.

{ return etaMax; }

double BPHX3872ToJPsiPiPiBuilder::getPiPtMin

(

)

const

get current cuts

Definition at line 292 of file BPHX3872ToJPsiPiPiBuilder.cc.

References ptMin.

{ return ptMin; }

double BPHX3872ToJPsiPiPiBuilder::getProbMin

(

)

const

Definition at line 300 of file BPHX3872ToJPsiPiPiBuilder.cc.

References chi2Sel, and BPHChi2Select::getProbMin().

{ return (chi2Sel == nullptr ? -1.0 : chi2Sel->getProbMin()); }

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

delete

void BPHX3872ToJPsiPiPiBuilder::setConstr

(

bool

flag

)

Definition at line 281 of file BPHX3872ToJPsiPiPiBuilder.cc.

References massConstr, and updated.

Referenced by BPHWriteSpecificDecay::fill().

                                                   {
  updated = false;
  massConstr = flag;
  return;
}

void BPHX3872ToJPsiPiPiBuilder::setJPsiMassMax

(

double

m

)

Definition at line 232 of file BPHX3872ToJPsiPiPiBuilder.cc.

References jpsiSel, BPHMassCuts::setMassMax(), and updated.

Referenced by BPHWriteSpecificDecay::fill().

                                                       {
  updated = false;
  jpsiSel->setMassMax(m);
  return;
}

void BPHX3872ToJPsiPiPiBuilder::setJPsiMassMin

(

double

m

)

set cuts

Definition at line 226 of file BPHX3872ToJPsiPiPiBuilder.cc.

References jpsiSel, BPHMassCuts::setMassMin(), and updated.

Referenced by BPHWriteSpecificDecay::fill().

                                                       {
  updated = false;
  jpsiSel->setMassMin(m);
  return;
}

void BPHX3872ToJPsiPiPiBuilder::setMassFitMax

(

double

m

)

Definition at line 275 of file BPHX3872ToJPsiPiPiBuilder.cc.

References mFitSel, BPHMassCuts::setMassMax(), and updated.

Referenced by BPHWriteSpecificDecay::fill().

                                                      {
  updated = false;
  mFitSel->setMassMax(m);
  return;
}

void BPHX3872ToJPsiPiPiBuilder::setMassFitMin

(

double

m

)

Definition at line 269 of file BPHX3872ToJPsiPiPiBuilder.cc.

References mFitSel, BPHMassCuts::setMassMin(), and updated.

Referenced by BPHWriteSpecificDecay::fill().

                                                      {
  updated = false;
  mFitSel->setMassMin(m);
  return;
}

void BPHX3872ToJPsiPiPiBuilder::setMassMax

(

double

m

)

Definition at line 256 of file BPHX3872ToJPsiPiPiBuilder.cc.

References massSel, BPHMassCuts::setMassMax(), and updated.

Referenced by BPHWriteSpecificDecay::fill().

                                                   {
  updated = false;
  massSel->setMassMax(m);
  return;
}

void BPHX3872ToJPsiPiPiBuilder::setMassMin

(

double

m

)

Definition at line 250 of file BPHX3872ToJPsiPiPiBuilder.cc.

References massSel, BPHMassCuts::setMassMin(), and updated.

Referenced by BPHWriteSpecificDecay::fill().

                                                   {
  updated = false;
  massSel->setMassMin(m);
  return;
}

void BPHX3872ToJPsiPiPiBuilder::setPiEtaMax

(

double

eta

)

Definition at line 244 of file BPHX3872ToJPsiPiPiBuilder.cc.

References PVValHelper::eta, etaMax, and updated.

Referenced by BPHWriteSpecificDecay::fill().

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

void BPHX3872ToJPsiPiPiBuilder::setPiPtMin

(

double

pt

)

set cuts

Definition at line 238 of file BPHX3872ToJPsiPiPiBuilder.cc.

References DiDispStaMuonMonitor_cfi::pt, ptMin, and updated.

Referenced by BPHWriteSpecificDecay::fill().

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

void BPHX3872ToJPsiPiPiBuilder::setProbMin

(

double

p

)

Definition at line 262 of file BPHX3872ToJPsiPiPiBuilder.cc.

References chi2Sel, and updated.

Referenced by BPHWriteSpecificDecay::fill().

                                                   {
  updated = false;
  delete chi2Sel;
  chi2Sel = (p < 0.0 ? nullptr : new BPHChi2Select(p));
  return;
}

Member Data Documentation

BPHChi2Select* BPHX3872ToJPsiPiPiBuilder::chi2Sel

private

Definition at line 103 of file BPHX3872ToJPsiPiPiBuilder.h.

Referenced by BPHX3872ToJPsiPiPiBuilder(), build(), getProbMin(), setProbMin(), and ~BPHX3872ToJPsiPiPiBuilder().

double BPHX3872ToJPsiPiPiBuilder::etaMax

private

Definition at line 100 of file BPHX3872ToJPsiPiPiBuilder.h.

Referenced by BPHX3872ToJPsiPiPiBuilder(), build(), getPiEtaMax(), and setPiEtaMax().

const edm::EventSetup* BPHX3872ToJPsiPiPiBuilder::evSetup

private

Definition at line 93 of file BPHX3872ToJPsiPiPiBuilder.h.

Referenced by build().

const std::vector<BPHPlusMinusConstCandPtr>* BPHX3872ToJPsiPiPiBuilder::jCollection

private

Definition at line 94 of file BPHX3872ToJPsiPiPiBuilder.h.

Referenced by build().

std::string BPHX3872ToJPsiPiPiBuilder::jPsiName

private

Definition at line 89 of file BPHX3872ToJPsiPiPiBuilder.h.

Referenced by BPHX3872ToJPsiPiPiBuilder(), and build().

BPHMassSelect* BPHX3872ToJPsiPiPiBuilder::jpsiSel

private

Definition at line 98 of file BPHX3872ToJPsiPiPiBuilder.h.

Referenced by BPHX3872ToJPsiPiPiBuilder(), build(), getJPsiMassMax(), getJPsiMassMin(), setJPsiMassMax(), setJPsiMassMin(), and ~BPHX3872ToJPsiPiPiBuilder().

bool BPHX3872ToJPsiPiPiBuilder::massConstr

private

Definition at line 106 of file BPHX3872ToJPsiPiPiBuilder.h.

Referenced by BPHX3872ToJPsiPiPiBuilder(), getConstr(), and setConstr().

BPHMassSelect* BPHX3872ToJPsiPiPiBuilder::massSel

private

Definition at line 102 of file BPHX3872ToJPsiPiPiBuilder.h.

Referenced by BPHX3872ToJPsiPiPiBuilder(), build(), getMassMax(), getMassMin(), setMassMax(), setMassMin(), and ~BPHX3872ToJPsiPiPiBuilder().

BPHMassFitSelect* BPHX3872ToJPsiPiPiBuilder::mFitSel

private

Definition at line 104 of file BPHX3872ToJPsiPiPiBuilder.h.

Referenced by BPHX3872ToJPsiPiPiBuilder(), build(), getMassFitMax(), getMassFitMin(), setMassFitMax(), setMassFitMin(), and ~BPHX3872ToJPsiPiPiBuilder().

float BPHX3872ToJPsiPiPiBuilder::minPDiff

private

Definition at line 107 of file BPHX3872ToJPsiPiPiBuilder.h.

Referenced by BPHX3872ToJPsiPiPiBuilder().

const BPHRecoBuilder::BPHGenericCollection* BPHX3872ToJPsiPiPiBuilder::nCollection

private

Definition at line 96 of file BPHX3872ToJPsiPiPiBuilder.h.

Referenced by build().

const BPHRecoBuilder::BPHGenericCollection* BPHX3872ToJPsiPiPiBuilder::pCollection

private

Definition at line 95 of file BPHX3872ToJPsiPiPiBuilder.h.

Referenced by build().

std::string BPHX3872ToJPsiPiPiBuilder::pionNegName

private

Definition at line 91 of file BPHX3872ToJPsiPiPiBuilder.h.

Referenced by build().

std::string BPHX3872ToJPsiPiPiBuilder::pionPosName

private

Definition at line 90 of file BPHX3872ToJPsiPiPiBuilder.h.

Referenced by build().

double BPHX3872ToJPsiPiPiBuilder::ptMin

private

Definition at line 99 of file BPHX3872ToJPsiPiPiBuilder.h.

Referenced by BPHX3872ToJPsiPiPiBuilder(), build(), getPiPtMin(), and setPiPtMin().

bool BPHX3872ToJPsiPiPiBuilder::updated

private

Definition at line 108 of file BPHX3872ToJPsiPiPiBuilder.h.

Referenced by BPHX3872ToJPsiPiPiBuilder(), build(), setConstr(), setJPsiMassMax(), setJPsiMassMin(), setMassFitMax(), setMassFitMin(), setMassMax(), setMassMin(), setPiEtaMax(), setPiPtMin(), and setProbMin().

std::vector<BPHRecoConstCandPtr> BPHX3872ToJPsiPiPiBuilder::x3872List

private

Definition at line 110 of file BPHX3872ToJPsiPiPiBuilder.h.

Referenced by build().