TtDilepLRSignalSelObservables Class Reference

#include <TtDilepLRSignalSelObservables.h>

Public Types
typedef std::pair< unsigned int, bool >	IntBoolPair

Public Member Functions
std::vector< IntBoolPair >	operator() (TtDilepEvtSolution &, const edm::Event &iEvent, bool matchOnly=false)
	TtDilepLRSignalSelObservables (edm::ConsumesCollector &&iC, const edm::EDGetTokenT< std::vector< pat::Jet > > &jetSourceToken)
	~TtDilepLRSignalSelObservables ()

Private Types
typedef std::pair< unsigned int, double >	IntDblPair

Private Member Functions
double	delta (double phi1, double phi2)
void	fillMinMax (double v1, double v2, int obsNbr, std::vector< IntDblPair > &varList, bool match1, bool match2, std::vector< IntBoolPair > &matchList)

Private Attributes
int	count1
int	count2
int	count3
int	count4
int	count5
std::vector< IntBoolPair >	evtselectVarMatch
std::vector< IntDblPair >	evtselectVarVal
edm::EDGetTokenT< TtGenEvent >	genEvtToken_
edm::EDGetTokenT< std::vector < pat::Jet > >	jetSourceToken_

Detailed Description

Definition at line 16 of file TtDilepLRSignalSelObservables.h.

Member Typedef Documentation

typedef std::pair<unsigned int, bool> TtDilepLRSignalSelObservables::IntBoolPair

Definition at line 22 of file TtDilepLRSignalSelObservables.h.

typedef std::pair<unsigned int, double> TtDilepLRSignalSelObservables::IntDblPair

private

Definition at line 26 of file TtDilepLRSignalSelObservables.h.

Constructor & Destructor Documentation

TtDilepLRSignalSelObservables::TtDilepLRSignalSelObservables	(	edm::ConsumesCollector &&	iC,
		const edm::EDGetTokenT< std::vector< pat::Jet > > &	jetSourceToken
	)

Definition at line 11 of file TtDilepLRSignalSelObservables.cc.

References count1, count2, count3, count4, and count5.

    : jetSourceToken_(jetSourceToken), genEvtToken_(iC.consumes<TtGenEvent>(edm::InputTag("genEvt"))) {
  count1 = 0;
  count2 = 0;
  count3 = 0;
  count4 = 0;
  count5 = 0;
  count3 = 0;
}

TtDilepLRSignalSelObservables::~TtDilepLRSignalSelObservables

(

)

Definition at line 23 of file TtDilepLRSignalSelObservables.cc.

{}

Member Function Documentation

double TtDilepLRSignalSelObservables::delta ( double  phi1, double  phi2  )

private

Definition at line 219 of file TtDilepLRSignalSelObservables.cc.

References srCondWrite_cfg::deltaPhi, and M_PI.

Referenced by operator()().

                                                                    {
  double deltaPhi = phi1 - phi2;
  while (deltaPhi > M_PI)
    deltaPhi -= 2 * M_PI;
  while (deltaPhi <= -M_PI)
    deltaPhi += 2 * M_PI;
  return deltaPhi;
}

void TtDilepLRSignalSelObservables::fillMinMax ( double  v1, double  v2, int  obsNbr, std::vector< IntDblPair > &  varList, bool  match1, bool  match2, std::vector< IntBoolPair > &  matchList  )

private

Definition at line 198 of file TtDilepLRSignalSelObservables.cc.

Referenced by operator()().

                                                                                  {
  if (v1 < v2) {
    varList.push_back(IntDblPair(obsNbr, v1));
    varList.push_back(IntDblPair(obsNbr + 1, v2));
    matchList.push_back(IntBoolPair(obsNbr, match1));
    matchList.push_back(IntBoolPair(obsNbr + 1, match2));

  } else {
    varList.push_back(IntDblPair(obsNbr, v2));
    varList.push_back(IntDblPair(obsNbr + 1, v1));
    matchList.push_back(IntBoolPair(obsNbr, match2));
    matchList.push_back(IntBoolPair(obsNbr + 1, match1));
  }
}

std::vector< TtDilepLRSignalSelObservables::IntBoolPair > TtDilepLRSignalSelObservables::operator()	(	TtDilepEvtSolution &	solution,
		const edm::Event &	iEvent,
		bool	matchOnly = false
	)

Definition at line 25 of file TtDilepLRSignalSelObservables.cc.

References funct::abs(), count1, count2, count3, count4, count5, delta(), srCondWrite_cfg::deltaPhi, runTauDisplay::dr, reco::LeafCandidate::et(), evtselectVarMatch, evtselectVarVal, edm::HandleBase::failedToGet(), fillMinMax(), MCTruth::genEvent, genEvtToken_, edm::Event::getByToken(), TtDilepEvtSolution::getCalJetB(), TtDilepEvtSolution::getCalJetBbar(), TtDilepEvtSolution::getGenLepm(), TtDilepEvtSolution::getGenLepp(), TtDilepEvtSolution::getJetB(), TtDilepEvtSolution::getJetBbar(), TtDilepEvtSolution::getLeptNeg(), TtDilepEvtSolution::getLeptPos(), TtDilepEvtSolution::getWmDecay(), TtDilepEvtSolution::getWpDecay(), mps_fire::i, edm::HandleBase::isValid(), fwrapper::jets, jetSourceToken_, M_PI, ResonanceBuilder::mass, reco::Particle::p4(), reco::LeafCandidate::p4(), pat::Jet::partonFlavour(), reco::LeafCandidate::pdgId(), createTree::pp, and TtDilepEvtSolution::setLRSignalEvtObservables().

                                                                          {
  evtselectVarVal.clear();
  evtselectVarMatch.clear();

  // Check whether the objects are matched:
  bool matchB1 = false;
  bool matchB2 = false;
  bool matchB = false;
  bool matchBbar = false;
  bool matchLeptPos = false;
  bool matchLeptNeg = false;

  edm::Handle<TtGenEvent> genEvent;
  iEvent.getByToken(genEvtToken_, genEvent);

  if (!genEvent.failedToGet() && genEvent.isValid()) {
    // std::cout <<std::endl;
    double dr, dr1, dr2;

    if (genEvent->isFullLeptonic()) {
      // Match the leptons, by type and deltaR
      dr = DeltaR<reco::Particle, reco::GenParticle>()(solution.getLeptPos(), *(solution.getGenLepp()));
      matchLeptPos = ((((solution.getWpDecay() == "electron") && (std::abs(solution.getGenLepp()->pdgId()) == 11)) ||
                       ((solution.getWpDecay() == "muon") && (std::abs(solution.getGenLepp()->pdgId()) == 13))) &&
                      (dr < 0.1));

      dr = DeltaR<reco::Particle, reco::GenParticle>()(solution.getLeptNeg(), *(solution.getGenLepm()));
      matchLeptNeg = ((((solution.getWmDecay() == "electron") && (std::abs(solution.getGenLepm()->pdgId()) == 11)) ||
                       ((solution.getWmDecay() == "muon") && (std::abs(solution.getGenLepm()->pdgId()) == 13))) &&
                      (dr < 0.1));
    }

    if (genEvent->isSemiLeptonic()) {
      int id = genEvent->singleLepton()->pdgId();

      if (id > 0) {
        dr = DeltaR<reco::Particle, reco::GenParticle>()(solution.getLeptNeg(), *(genEvent->singleLepton()));
        matchLeptNeg = ((((solution.getWmDecay() == "electron") && (id == 11)) ||
                         ((solution.getWmDecay() == "muon") && (id == 13))) &&
                        (dr < 0.1));
      } else {
        dr = DeltaR<reco::Particle, reco::GenParticle>()(solution.getLeptPos(), *(genEvent->singleLepton()));
        matchLeptPos = ((((solution.getWpDecay() == "electron") && (id == -11)) ||
                         ((solution.getWpDecay() == "muon") && (id == -13))) &&
                        (dr < 0.1));
      }
    }

    if (genEvent->isTtBar() && genEvent->numberOfBQuarks() > 1) {
      if (solution.getJetB().partonFlavour() == 5)
        ++count1;
      if (solution.getJetBbar().partonFlavour() == 5)
        ++count1;

      dr1 = DeltaR<pat::Jet, reco::GenParticle>()(solution.getCalJetB(), *(genEvent->b()));
      dr2 = DeltaR<pat::Jet, reco::GenParticle>()(solution.getCalJetB(), *(genEvent->bBar()));

      matchB1 = ((dr1 < 0.4) || (dr2 < 0.4));
      matchB = ((solution.getJetB().partonFlavour() == 5) && (dr1 < 0.4));
      if (matchB)
        ++count3;
      matchB = ((dr1 < 0.4));
      if (dr1 < 0.5)
        ++count2;
      if (dr1 < 0.4)
        ++count4;
      if (dr1 < 0.3)
        ++count5;

      dr1 = DeltaR<pat::Jet, reco::GenParticle>()(solution.getCalJetBbar(), *(genEvent->b()));
      dr2 = DeltaR<pat::Jet, reco::GenParticle>()(solution.getCalJetBbar(), *(genEvent->bBar()));

      matchBbar = ((solution.getJetBbar().partonFlavour() == 5) && (dr2 < 0.4));
      if (matchBbar)
        ++count3;
      matchBbar = ((dr2 < 0.4));
      matchB2 = ((dr1 < 0.4) || (dr2 < 0.4));
      if (dr2 < 0.5)
        ++count2;
      if (dr2 < 0.4)
        ++count4;
      if (dr2 < 0.3)
        ++count5;
    }
  }

  edm::Handle<std::vector<pat::Jet> > jets;
  iEvent.getByToken(jetSourceToken_, jets);

  //  Lower / Higher of both jet angles

  double v1 = std::abs(solution.getJetB().p4().theta() - M_PI / 2);
  double v2 = std::abs(solution.getJetBbar().p4().theta() - M_PI / 2);
  fillMinMax(v1, v2, 1, evtselectVarVal, matchB1, matchB2, evtselectVarMatch);

  //  Lower / Higher of both jet pT

  double pt1 = solution.getJetB().p4().pt();
  double pt2 = solution.getJetBbar().p4().pt();
  fillMinMax(pt1, pt2, 3, evtselectVarVal, matchB1, matchB2, evtselectVarMatch);

  //  Lower / Higher of both lepton pT

  pt1 = solution.getLeptPos().p4().pt();
  pt2 = solution.getLeptNeg().p4().pt();
  fillMinMax(pt1, pt2, 5, evtselectVarVal, matchLeptPos, matchLeptNeg, evtselectVarMatch);

  // delta theta btw the b-jets

  double deltaPhi = std::abs(delta(solution.getJetB().p4().phi(), solution.getJetBbar().p4().phi()));

  evtselectVarVal.push_back(IntDblPair(7, deltaPhi));
  evtselectVarMatch.push_back(IntBoolPair(7, matchB1 && matchB2));

  // delta phi btw the b-jets

  double deltaTheta = std::abs(delta(solution.getJetBbar().p4().theta(), solution.getJetB().p4().theta()));

  evtselectVarVal.push_back(IntDblPair(8, deltaTheta));
  evtselectVarMatch.push_back(IntBoolPair(8, matchB1 && matchB2));

  //  Lower / Higher of phi difference between the b and associated lepton

  double deltaPhi1 = std::abs(delta(solution.getJetB().p4().phi(), solution.getLeptPos().p4().phi()));
  double deltaPhi2 = std::abs(delta(solution.getJetBbar().p4().phi(), solution.getLeptNeg().p4().phi()));

  fillMinMax(
      deltaPhi1, deltaPhi2, 9, evtselectVarVal, matchB && matchLeptPos, matchBbar && matchLeptNeg, evtselectVarMatch);

  //  Lower / Higher of theta difference between the b and associated lepton

  double deltaTheta1 = std::abs(solution.getJetB().p4().theta() - solution.getLeptPos().p4().theta());
  double deltaTheta2 = std::abs(solution.getJetBbar().p4().theta() - solution.getLeptNeg().p4().theta());
  fillMinMax(deltaTheta1,
             deltaTheta2,
             11,
             evtselectVarVal,
             matchB && matchLeptPos,
             matchBbar && matchLeptNeg,
             evtselectVarMatch);

  // Invariant Mass of lepton pair

  math::XYZTLorentzVector pp = solution.getLeptPos().p4() + solution.getLeptNeg().p4();
  double mass = pp.mass();
  evtselectVarVal.push_back(IntDblPair(13, mass));
  evtselectVarMatch.push_back(IntBoolPair(13, matchLeptNeg && matchLeptPos));

  evtselectVarVal.push_back(IntDblPair(13, mass));
  evtselectVarMatch.push_back(IntBoolPair(13, matchLeptNeg && matchLeptPos));

  std::vector<pat::Jet> jet3;
  for (int i = 0; i < (int)jets->size(); ++i) {
    if (((*jets)[i].et() < solution.getJetB().et()) && ((*jets)[i].et() < solution.getJetBbar().et())) {
      jet3.push_back((*jets)[i]);
    }
  }
  double jet1Ratio = 0., jet2Ratio = 0.;
  if (!jet3.empty()) {
    jet1Ratio = jet3[0].et() / solution.getJetB().et();
    jet2Ratio = jet3[0].et() / solution.getJetBbar().et();
  }
  fillMinMax(jet1Ratio, jet2Ratio, 14, evtselectVarVal, matchB1, matchB2, evtselectVarMatch);

  evtselectVarVal.push_back(IntDblPair(16, jets->size()));
  evtselectVarMatch.push_back(IntBoolPair(16, matchB && matchBbar));

  if (!matchOnly)
    solution.setLRSignalEvtObservables(evtselectVarVal);
  return evtselectVarMatch;
}

Member Data Documentation

int TtDilepLRSignalSelObservables::count1

private

Definition at line 42 of file TtDilepLRSignalSelObservables.h.

Referenced by operator()(), and TtDilepLRSignalSelObservables().

int TtDilepLRSignalSelObservables::count2

private

Definition at line 42 of file TtDilepLRSignalSelObservables.h.

Referenced by operator()(), and TtDilepLRSignalSelObservables().

int TtDilepLRSignalSelObservables::count3

private

Definition at line 42 of file TtDilepLRSignalSelObservables.h.

Referenced by operator()(), and TtDilepLRSignalSelObservables().

int TtDilepLRSignalSelObservables::count4

private

Definition at line 42 of file TtDilepLRSignalSelObservables.h.

Referenced by operator()(), and TtDilepLRSignalSelObservables().

int TtDilepLRSignalSelObservables::count5

private

Definition at line 42 of file TtDilepLRSignalSelObservables.h.

Referenced by operator()(), and TtDilepLRSignalSelObservables().

std::vector<IntBoolPair> TtDilepLRSignalSelObservables::evtselectVarMatch

private

Definition at line 41 of file TtDilepLRSignalSelObservables.h.

Referenced by operator()().

std::vector<IntDblPair> TtDilepLRSignalSelObservables::evtselectVarVal

private

Definition at line 40 of file TtDilepLRSignalSelObservables.h.

Referenced by operator()().

edm::EDGetTokenT<TtGenEvent> TtDilepLRSignalSelObservables::genEvtToken_

private

Definition at line 38 of file TtDilepLRSignalSelObservables.h.

Referenced by operator()().

edm::EDGetTokenT<std::vector<pat::Jet> > TtDilepLRSignalSelObservables::jetSourceToken_

private

Definition at line 37 of file TtDilepLRSignalSelObservables.h.

Referenced by operator()().