BiasedTauDecayer Class Reference

#include <BiasedTauDecayer.h>

Inheritance diagram for BiasedTauDecayer:

Public Member Functions
	BiasedTauDecayer (Pythia8::Info *infoPtr, Pythia8::Settings *settingsPtr)
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::BiasedTauDecayer	(	Pythia8::Info *	infoPtr,
		Pythia8::Settings *	settingsPtr
	)

Definition at line 11 of file BiasedTauDecayer.cc.

                                                                       {
  decayer = TauDecays();
  decayer.initInfoPtr(*infoPtr);
  decayer.init();
  filter_ = settingsPtr->flag("BiasedTauDecayer:filter");
  eDecays_ = settingsPtr->flag("BiasedTauDecayer:eDecays");
  muDecays_ = settingsPtr->flag("BiasedTauDecayer:muDecays");
}

Member Function Documentation

decay()

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.

                                                                                                            {
  if (!filter_)
    return false;
  if (idProd[0] != 15 && idProd[0] != -15)
    return false;
  int iStart = event[iDec].iTopCopyId();
  int iMom = event[iStart].mother1();
  int idMom = event[iMom].idAbs();
  if (idMom != 23 && idMom != 24 && idMom != 25)
    return false;
  int iDau1 = event[iMom].daughter1();
  int iDau2 = event[iMom].daughter2();
  int iBot1 = event[iDau1].iBotCopyId();
  int iBot2 = event[iDau2].iBotCopyId();
  int iDecSis = (iDec == iBot1) ? iBot2 : iBot1;
  // Check if sister has been decayed
  // Since taus decays are correlated, use one decay, store the other
  bool notDecayed = event[iDecSis].status() > 0 ? true : false;
  if (notDecayed) {
    // bias for leptonic decays
    bool hasLepton = (eDecays_ || muDecays_) ? false : true;
    Event decay;
    int i1 = -1;
    int i2 = -1;
    while (!hasLepton) {
      decay = event;
      decayer.decay(iDec, decay);
      // check for lepton in first decay
      i1 = decay[iDec].daughter1();
      i2 = decay[iDec].daughter2();
      for (int i = i1; i < i2 + 1; ++i) {
        if (decay[i].isLepton() && decay[i].isCharged()) {
          if ((eDecays_ && abs(decay[i].id()) == 11) || (muDecays_ && abs(decay[i].id()) == 13)) {
            hasLepton = true;
            break;
          }
        }
      }
      if (hasLepton)
        break;
      // check for lepton in second decay
      i1 = decay[iDecSis].daughter1();
      i2 = decay[iDecSis].daughter2();
      for (int i = i1; i < i2 + 1; ++i) {
        if (decay[i].isLepton() && decay[i].isCharged()) {
          if ((eDecays_ && abs(decay[i].id()) == 11) || (muDecays_ && abs(decay[i].id()) == 13)) {
            hasLepton = true;
            break;
          }
        }
      }
    }
    // Return decay products
    i1 = decay[iDec].daughter1();
    i2 = decay[iDec].daughter2();
    for (int i = i1; i < i2 + 1; ++i) {
      idProd.push_back(decay[i].id());
      mProd.push_back(decay[i].m());
      pProd.push_back(decay[i].p());
    }
    // Store correlated decay products
    i1 = decay[iDecSis].daughter1();
    i2 = decay[iDecSis].daughter2();
    idProdSave.clear();
    mProdSave.clear();
    pProdSave.clear();
    for (int i = i1; i < i2 + 1; ++i) {
      idProdSave.push_back(decay[i].id());
      mProdSave.push_back(decay[i].m());
      pProdSave.push_back(decay[i].p());
    }
  } else {
    // Return stored decay products
    for (size_t i = 0; i < idProdSave.size(); ++i) {
      idProd.push_back(idProdSave[i]);
      mProd.push_back(mProdSave[i]);
      pProd.push_back(pProdSave[i]);
    }
  }
 
  return true;
}

References funct::abs(), PA_ZEESkim_cff::decay, edmPickEvents::event, mps_fire::i, testProducerWithPsetDescEmpty_cfi::i1, testProducerWithPsetDescEmpty_cfi::i2, reco::isLepton(), visualization-live-secondInstance_cfg::m, AlCaHLTBitMon_ParallelJobs::p, and funct::true.

Member Data Documentation

decayer

Pythia8::TauDecays BiasedTauDecayer::decayer

private

Definition at line 20 of file BiasedTauDecayer.h.

eDecays_

bool BiasedTauDecayer::eDecays_

private

Definition at line 22 of file BiasedTauDecayer.h.

filter_

bool BiasedTauDecayer::filter_

private

Definition at line 21 of file BiasedTauDecayer.h.

idProdSave

std::vector<int> BiasedTauDecayer::idProdSave

private

Definition at line 24 of file BiasedTauDecayer.h.

mProdSave

std::vector<double> BiasedTauDecayer::mProdSave

private

Definition at line 25 of file BiasedTauDecayer.h.

muDecays_

bool BiasedTauDecayer::muDecays_

private

Definition at line 23 of file BiasedTauDecayer.h.

pProdSave

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

private

Definition at line 26 of file BiasedTauDecayer.h.