BPHX3872ToJPsiPiPiBuilder Class Reference

#include <BPHX3872ToJPsiPiPiBuilder.h>

Public Member Functions
	BPHX3872ToJPsiPiPiBuilder (const BPHX3872ToJPsiPiPiBuilder &x)=delete
	BPHX3872ToJPsiPiPiBuilder (const edm::EventSetup &es, const std::vector< BPHPlusMinusConstCandPtr > &jpsiCollection, const BPHRecoBuilder::BPHGenericCollection *posCollection, const BPHRecoBuilder::BPHGenericCollection *negCollection)
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() [1/2]

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.

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

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

BPHX3872ToJPsiPiPiBuilder() [2/2]

BPHX3872ToJPsiPiPiBuilder::BPHX3872ToJPsiPiPiBuilder ( const BPHX3872ToJPsiPiPiBuilder &  x )

delete

~BPHX3872ToJPsiPiPiBuilder()

BPHX3872ToJPsiPiPiBuilder::~BPHX3872ToJPsiPiPiBuilder ( )

virtual

Destructor

Definition at line 61 of file BPHX3872ToJPsiPiPiBuilder.cc.

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

References chi2Sel, jpsiSel, massSel, and mFitSel.

Member Function Documentation

build()

vector< BPHRecoConstCandPtr > BPHX3872ToJPsiPiPiBuilder::build

(

)

build X3872 candidates

Operations

Definition at line 71 of file BPHX3872ToJPsiPiPiBuilder.cc.

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

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

Referenced by BPHWriteSpecificDecay::fill().

getConstr()

bool BPHX3872ToJPsiPiPiBuilder::getConstr

(

)

const

Definition at line 306 of file BPHX3872ToJPsiPiPiBuilder.cc.

{ return massConstr; }

References massConstr.

getJPsiMassMax()

double BPHX3872ToJPsiPiPiBuilder::getJPsiMassMax

(

)

const

Definition at line 290 of file BPHX3872ToJPsiPiPiBuilder.cc.

{ return jpsiSel->getMassMax(); }

References BPHMassCuts::getMassMax(), and jpsiSel.

getJPsiMassMin()

double BPHX3872ToJPsiPiPiBuilder::getJPsiMassMin

(

)

const

get current cuts

Definition at line 288 of file BPHX3872ToJPsiPiPiBuilder.cc.

{ return jpsiSel->getMassMin(); }

References BPHMassCuts::getMassMin(), and jpsiSel.

getMassFitMax()

double BPHX3872ToJPsiPiPiBuilder::getMassFitMax

(

)

const

Definition at line 304 of file BPHX3872ToJPsiPiPiBuilder.cc.

{ return mFitSel->getMassMax(); }

References BPHMassCuts::getMassMax(), and mFitSel.

getMassFitMin()

double BPHX3872ToJPsiPiPiBuilder::getMassFitMin

(

)

const

Definition at line 302 of file BPHX3872ToJPsiPiPiBuilder.cc.

{ return mFitSel->getMassMin(); }

References BPHMassCuts::getMassMin(), and mFitSel.

getMassMax()

double BPHX3872ToJPsiPiPiBuilder::getMassMax

(

)

const

Definition at line 298 of file BPHX3872ToJPsiPiPiBuilder.cc.

{ return massSel->getMassMax(); }

References BPHMassCuts::getMassMax(), and massSel.

Referenced by build().

getMassMin()

double BPHX3872ToJPsiPiPiBuilder::getMassMin

(

)

const

Definition at line 296 of file BPHX3872ToJPsiPiPiBuilder.cc.

{ return massSel->getMassMin(); }

References BPHMassCuts::getMassMin(), and massSel.

getPiEtaMax()

double BPHX3872ToJPsiPiPiBuilder::getPiEtaMax

(

)

const

Definition at line 294 of file BPHX3872ToJPsiPiPiBuilder.cc.

{ return etaMax; }

References etaMax.

getPiPtMin()

double BPHX3872ToJPsiPiPiBuilder::getPiPtMin

(

)

const

get current cuts

Definition at line 292 of file BPHX3872ToJPsiPiPiBuilder.cc.

{ return ptMin; }

References ptMin.

getProbMin()

double BPHX3872ToJPsiPiPiBuilder::getProbMin

(

)

const

Definition at line 300 of file BPHX3872ToJPsiPiPiBuilder.cc.

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

References chi2Sel, and BPHChi2Select::getProbMin().

operator=()

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

delete

setConstr()

void BPHX3872ToJPsiPiPiBuilder::setConstr

(

bool

flag

)

Definition at line 281 of file BPHX3872ToJPsiPiPiBuilder.cc.

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

References RemoveAddSevLevel::flag, massConstr, and updated.

Referenced by BPHWriteSpecificDecay::fill().

setJPsiMassMax()

void BPHX3872ToJPsiPiPiBuilder::setJPsiMassMax

(

double

m

)

Definition at line 232 of file BPHX3872ToJPsiPiPiBuilder.cc.

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

References jpsiSel, visualization-live-secondInstance_cfg::m, BPHMassCuts::setMassMax(), and updated.

Referenced by BPHWriteSpecificDecay::fill().

setJPsiMassMin()

void BPHX3872ToJPsiPiPiBuilder::setJPsiMassMin

(

double

m

)

set cuts

Definition at line 226 of file BPHX3872ToJPsiPiPiBuilder.cc.

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

References jpsiSel, visualization-live-secondInstance_cfg::m, BPHMassCuts::setMassMin(), and updated.

Referenced by BPHWriteSpecificDecay::fill().

setMassFitMax()

void BPHX3872ToJPsiPiPiBuilder::setMassFitMax

(

double

m

)

Definition at line 275 of file BPHX3872ToJPsiPiPiBuilder.cc.

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

References visualization-live-secondInstance_cfg::m, mFitSel, BPHMassCuts::setMassMax(), and updated.

Referenced by BPHWriteSpecificDecay::fill().

setMassFitMin()

void BPHX3872ToJPsiPiPiBuilder::setMassFitMin

(

double

m

)

Definition at line 269 of file BPHX3872ToJPsiPiPiBuilder.cc.

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

References visualization-live-secondInstance_cfg::m, mFitSel, BPHMassCuts::setMassMin(), and updated.

Referenced by BPHWriteSpecificDecay::fill().

setMassMax()

void BPHX3872ToJPsiPiPiBuilder::setMassMax

(

double

m

)

Definition at line 256 of file BPHX3872ToJPsiPiPiBuilder.cc.

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

References visualization-live-secondInstance_cfg::m, massSel, BPHMassCuts::setMassMax(), and updated.

Referenced by BPHWriteSpecificDecay::fill().

setMassMin()

void BPHX3872ToJPsiPiPiBuilder::setMassMin

(

double

m

)

Definition at line 250 of file BPHX3872ToJPsiPiPiBuilder.cc.

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

References visualization-live-secondInstance_cfg::m, massSel, BPHMassCuts::setMassMin(), and updated.

Referenced by BPHWriteSpecificDecay::fill().

setPiEtaMax()

void BPHX3872ToJPsiPiPiBuilder::setPiEtaMax

(

double

eta

)

Definition at line 244 of file BPHX3872ToJPsiPiPiBuilder.cc.

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

References PVValHelper::eta, etaMax, and updated.

Referenced by BPHWriteSpecificDecay::fill().

setPiPtMin()

void BPHX3872ToJPsiPiPiBuilder::setPiPtMin

(

double

pt

)

set cuts

Definition at line 238 of file BPHX3872ToJPsiPiPiBuilder.cc.

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

References DiDispStaMuonMonitor_cfi::pt, ptMin, and updated.

Referenced by BPHWriteSpecificDecay::fill().

setProbMin()

void BPHX3872ToJPsiPiPiBuilder::setProbMin

(

double

p

)

Definition at line 262 of file BPHX3872ToJPsiPiPiBuilder.cc.

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

References chi2Sel, AlCaHLTBitMon_ParallelJobs::p, and updated.

Referenced by BPHWriteSpecificDecay::fill().

Member Data Documentation

chi2Sel

BPHChi2Select* BPHX3872ToJPsiPiPiBuilder::chi2Sel

private

Definition at line 103 of file BPHX3872ToJPsiPiPiBuilder.h.

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

etaMax

double BPHX3872ToJPsiPiPiBuilder::etaMax

private

Definition at line 100 of file BPHX3872ToJPsiPiPiBuilder.h.

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

evSetup

const edm::EventSetup* BPHX3872ToJPsiPiPiBuilder::evSetup

private

Definition at line 93 of file BPHX3872ToJPsiPiPiBuilder.h.

Referenced by build().

jCollection

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

private

Definition at line 94 of file BPHX3872ToJPsiPiPiBuilder.h.

Referenced by build().

jPsiName

std::string BPHX3872ToJPsiPiPiBuilder::jPsiName

private

Definition at line 89 of file BPHX3872ToJPsiPiPiBuilder.h.

Referenced by BPHX3872ToJPsiPiPiBuilder(), and build().

jpsiSel

BPHMassSelect* BPHX3872ToJPsiPiPiBuilder::jpsiSel

private

Definition at line 98 of file BPHX3872ToJPsiPiPiBuilder.h.

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

massConstr

bool BPHX3872ToJPsiPiPiBuilder::massConstr

private

Definition at line 106 of file BPHX3872ToJPsiPiPiBuilder.h.

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

massSel

BPHMassSelect* BPHX3872ToJPsiPiPiBuilder::massSel

private

Definition at line 102 of file BPHX3872ToJPsiPiPiBuilder.h.

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

mFitSel

BPHMassFitSelect* BPHX3872ToJPsiPiPiBuilder::mFitSel

private

Definition at line 104 of file BPHX3872ToJPsiPiPiBuilder.h.

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

minPDiff

float BPHX3872ToJPsiPiPiBuilder::minPDiff

private

Definition at line 107 of file BPHX3872ToJPsiPiPiBuilder.h.

Referenced by BPHX3872ToJPsiPiPiBuilder().

nCollection

const BPHRecoBuilder::BPHGenericCollection* BPHX3872ToJPsiPiPiBuilder::nCollection

private

Definition at line 96 of file BPHX3872ToJPsiPiPiBuilder.h.

Referenced by build().

pCollection

const BPHRecoBuilder::BPHGenericCollection* BPHX3872ToJPsiPiPiBuilder::pCollection

private

Definition at line 95 of file BPHX3872ToJPsiPiPiBuilder.h.

Referenced by build().

pionNegName

std::string BPHX3872ToJPsiPiPiBuilder::pionNegName

private

Definition at line 91 of file BPHX3872ToJPsiPiPiBuilder.h.

Referenced by build().

pionPosName

std::string BPHX3872ToJPsiPiPiBuilder::pionPosName

private

Definition at line 90 of file BPHX3872ToJPsiPiPiBuilder.h.

Referenced by build().

ptMin

double BPHX3872ToJPsiPiPiBuilder::ptMin

private

Definition at line 99 of file BPHX3872ToJPsiPiPiBuilder.h.

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

updated

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

x3872List

std::vector<BPHRecoConstCandPtr> BPHX3872ToJPsiPiPiBuilder::x3872List

private

Definition at line 110 of file BPHX3872ToJPsiPiPiBuilder.h.

Referenced by build().