CMS 3D CMS Logo
List of all members | Public Member Functions | Private Attributes
BiasedTauDecayer Class Reference

#include <BiasedTauDecayer.h>

Inheritance diagram for BiasedTauDecayer:

Public Member Functions

 BiasedTauDecayer (Pythia8::Info *infoPtr, Pythia8::Settings *settingsPtr, Pythia8::ParticleData *particleDataPtr, Pythia8::Rndm *rndmPtr, Pythia8::Couplings *couplingsPtr)
 
bool decay (std::vector< int > &idProd, std::vector< double > &mProd, std::vector< Pythia8::Vec4 > &pProd, int iDec, const Pythia8::Event &event) override
 

Private Attributes

Pythia8::TauDecays decayer
 
bool eDecays_
 
bool filter_
 
std::vector< int > idProdSave
 
std::vector< double > mProdSave
 
bool muDecays_
 
std::vector< Pythia8::Vec4 > pProdSave
 

Detailed Description

Definition at line 9 of file BiasedTauDecayer.h.

Constructor & Destructor Documentation

BiasedTauDecayer::BiasedTauDecayer ( Pythia8::Info *  infoPtr,
Pythia8::Settings *  settingsPtr,
Pythia8::ParticleData *  particleDataPtr,
Pythia8::Rndm *  rndmPtr,
Pythia8::Couplings *  couplingsPtr 
)

Definition at line 11 of file BiasedTauDecayer.cc.

12  {
13  decayer = TauDecays();
14  decayer.init(infoPtr, settingsPtr, particleDataPtr, rndmPtr, couplingsPtr);
15  filter_ = settingsPtr->flag("BiasedTauDecayer:filter");
16  eDecays_ = settingsPtr->flag("BiasedTauDecayer:eDecays");
17  muDecays_ = settingsPtr->flag("BiasedTauDecayer:muDecays");
18 }
BiasedTauDecayer::decayer
Pythia8::TauDecays decayer
Definition: BiasedTauDecayer.h:24

Member Function Documentation

bool BiasedTauDecayer::decay ( std::vector< int > &  idProd,
std::vector< double > &  mProd,
std::vector< Pythia8::Vec4 > &  pProd,
int  iDec,
const Pythia8::Event &  event 
)
override

Definition at line 20 of file BiasedTauDecayer.cc.

References funct::abs(), SelectingProcedure_cff::decay, event(), mps_fire::i, reco::isLepton(), funct::m, and AlCaHLTBitMon_ParallelJobs::p.

21  {
22  if (!filter_)
23  return false;
24  if (idProd[0] != 15 && idProd[0] != -15)
25  return false;
26  int iStart = event[iDec].iTopCopyId();
27  int iMom = event[iStart].mother1();
28  int idMom = event[iMom].idAbs();
29  if (idMom != 23 && idMom != 24 && idMom != 25)
30  return false;
31  int iDau1 = event[iMom].daughter1();
32  int iDau2 = event[iMom].daughter2();
33  int iBot1 = event[iDau1].iBotCopyId();
34  int iBot2 = event[iDau2].iBotCopyId();
35  int iDecSis = (iDec == iBot1) ? iBot2 : iBot1;
36  // Check if sister has been decayed
37  // Since taus decays are correlated, use one decay, store the other
38  bool notDecayed = event[iDecSis].status() > 0 ? true : false;
39  if (notDecayed) {
40  // bias for leptonic decays
41  bool hasLepton = (eDecays_ || muDecays_) ? false : true;
42  Event decay;
43  int i1 = -1;
44  int i2 = -1;
45  while (!hasLepton) {
46  decay = event;
47  decayer.decay(iDec, decay);
48  // check for lepton in first decay
49  i1 = decay[iDec].daughter1();
50  i2 = decay[iDec].daughter2();
51  for (int i = i1; i < i2 + 1; ++i) {
52  if (decay[i].isLepton() && decay[i].isCharged()) {
53  if ((eDecays_ && abs(decay[i].id()) == 11) || (muDecays_ && abs(decay[i].id()) == 13)) {
54  hasLepton = true;
55  break;
56  }
57  }
58  }
59  if (hasLepton)
60  break;
61  // check for lepton in second decay
62  i1 = decay[iDecSis].daughter1();
63  i2 = decay[iDecSis].daughter2();
64  for (int i = i1; i < i2 + 1; ++i) {
65  if (decay[i].isLepton() && decay[i].isCharged()) {
66  if ((eDecays_ && abs(decay[i].id()) == 11) || (muDecays_ && abs(decay[i].id()) == 13)) {
67  hasLepton = true;
68  break;
69  }
70  }
71  }
72  }
73  // Return decay products
74  i1 = decay[iDec].daughter1();
75  i2 = decay[iDec].daughter2();
76  for (int i = i1; i < i2 + 1; ++i) {
77  idProd.push_back(decay[i].id());
78  mProd.push_back(decay[i].m());
79  pProd.push_back(decay[i].p());
80  }
81  // Store correlated decay products
82  i1 = decay[iDecSis].daughter1();
83  i2 = decay[iDecSis].daughter2();
84  idProdSave.clear();
85  mProdSave.clear();
86  pProdSave.clear();
87  for (int i = i1; i < i2 + 1; ++i) {
88  idProdSave.push_back(decay[i].id());
89  mProdSave.push_back(decay[i].m());
90  pProdSave.push_back(decay[i].p());
91  }
92  } else {
93  // Return stored decay products
94  for (size_t i = 0; i < idProdSave.size(); ++i) {
95  idProd.push_back(idProdSave[i]);
96  mProd.push_back(mProdSave[i]);
97  pProd.push_back(pProdSave[i]);
98  }
99  }
100 
101  return true;
102 }
mps_fire.i
i
Definition: mps_fire.py:338
AlCaHLTBitMon_ParallelJobs.p
p
Definition: AlCaHLTBitMon_ParallelJobs.py:153
reco::isLepton
bool isLepton(const Candidate &part)
Definition: pdgIdUtils.h:19
BiasedTauDecayer::pProdSave
std::vector< Pythia8::Vec4 > pProdSave
Definition: BiasedTauDecayer.h:30
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
event
How EventSelector::AcceptEvent() decides whether to accept an event for output otherwise it is excluding the probing of A single or multiple positive and the trigger will pass if any such matching triggers are PASS or EXCEPTION[A criterion thatmatches no triggers at all is detected and causes a throw.] A single negative with an expectation of appropriate bit checking in the decision and the trigger will pass if any such matching triggers are FAIL or EXCEPTION A wildcarded negative criterion that matches more than one trigger in the trigger but the state exists so we define the behavior If all triggers are the negative crieriion will lead to accepting the event(this again matches the behavior of"!*"before the partial wildcard feature was incorporated).The per-event"cost"of each negative criterion with multiple relevant triggers is about the same as!*was in the past
BiasedTauDecayer::idProdSave
std::vector< int > idProdSave
Definition: BiasedTauDecayer.h:28
BiasedTauDecayer::mProdSave
std::vector< double > mProdSave
Definition: BiasedTauDecayer.h:29
BiasedTauDecayer::decay
bool decay(std::vector< int > &idProd, std::vector< double > &mProd, std::vector< Pythia8::Vec4 > &pProd, int iDec, const Pythia8::Event &event) override
Definition: BiasedTauDecayer.cc:20
BiasedTauDecayer::decayer
Pythia8::TauDecays decayer
Definition: BiasedTauDecayer.h:24

Member Data Documentation

Pythia8::TauDecays BiasedTauDecayer::decayer
private

Definition at line 24 of file BiasedTauDecayer.h.

bool BiasedTauDecayer::eDecays_
private

Definition at line 26 of file BiasedTauDecayer.h.

bool BiasedTauDecayer::filter_
private

Definition at line 25 of file BiasedTauDecayer.h.

std::vector<int> BiasedTauDecayer::idProdSave
private

Definition at line 28 of file BiasedTauDecayer.h.

std::vector<double> BiasedTauDecayer::mProdSave
private

Definition at line 29 of file BiasedTauDecayer.h.

bool BiasedTauDecayer::muDecays_
private

Definition at line 27 of file BiasedTauDecayer.h.

std::vector<Pythia8::Vec4> BiasedTauDecayer::pProdSave
private

Definition at line 30 of file BiasedTauDecayer.h.

Generated for CMSSW Reference Manual by doxygen 1.8.11