TauDecay_GenParticle Class Reference

#include <TauDecay_GenParticle.h>

Inheritance diagram for TauDecay_GenParticle:

Public Member Functions
bool	AddRadInfo (const reco::GenParticle *Particle, unsigned int midx)
bool	AnalyzeTau (const reco::GenParticle *Tau, unsigned int &MODE_ID, unsigned int &TauBitMask, bool dores, bool dopi0)
std::vector< unsigned int >	Get_MotherIdx ()
std::vector< const reco::GenParticle * >	Get_TauDecayProducts ()
	TauDecay_GenParticle ()
	~TauDecay_GenParticle ()
Public Member Functions inherited from TauDecay
void	ClassifyDecayMode (unsigned int &MODE_ID, unsigned int &TauBitMask)
bool	hasResonance (unsigned int &TauBitMask, int pdgid)
bool	isTauFinalStateParticle (int pdgid)
bool	isTauParticleCounter (int pdgid)
bool	isTauResonanceCounter (int pdgid)
unsigned int	nPi0 (unsigned int &TauBitMask)
unsigned int	nProng (unsigned int &TauBitMask)
void	Reset ()
	TauDecay ()
	~TauDecay ()

Private Member Functions
void	AddPi0Info (const reco::GenParticle *Particle, unsigned int midx)
void	Analyze (const reco::GenParticle *Particle, unsigned int midx, bool dores, bool dopi0)

Private Attributes
unsigned int	MODE_ID
std::vector< unsigned int >	MotherIdx
unsigned int	TauBitMask
std::vector< const reco::GenParticle * >	TauDecayProducts

Additional Inherited Members
Public Types inherited from TauDecay
enum	MODE {   MODE_UNKNOWN = 0, MODE_ELECTRON, MODE_MUON, MODE_PION,   MODE_PIPI0, MODE_3PI, MODE_PI2PI0, MODE_KAON,   MODE_K0PI, MODE_KPI0, MODE_3PIPI0, MODE_PI3PI0,   MODE_3PI2PI0, MODE_5PI, MODE_5PIPI0, MODE_3PI3PI0,   MODE_KPIK, MODE_K0BK0PI, MODE_KK0BPI0, MODE_K2PI0,   MODE_KPIPI, MODE_PIK0PI0, MODE_ETAPIPI0, MODE_PIPI0GAM,   MODE_KK0B, MODE_PI4PI0, MODE_3PIETA, MODE_PI2PI0ETA,   MODE_PI2PI0OMEGA, MODE_3PIOMEGA, MODE_KOMEGA, MODE_K3PI0,   MODE_K2PIPI0, MODE_KETA, MODE_K0PI2PI0, MODE_K03PI,   MODE_2K0PIPI0, MODE_KPIKPI0, NMODEID }
enum	TauDecayStructure {   other = 0, OneProng = 1, ThreeProng = 2, FiveProng = 4,   OnePi0 = 8, TwoPi0 = 32, ThreePi0 = 64, Res_a1_pm = 128,   Res_a1_0 = 256, Res_rho_pm = 512, Res_rho_0 = 1024, Res_eta = 2048,   Res_omega = 4096, Res_Kstar_pm = 8192, Res_Kstar_0 = 16384, KS0_to_pipi = 32768 }
Static Public Member Functions inherited from TauDecay
static std::string	DecayMode (unsigned int &MODE_ID)

Detailed Description

Definition at line 23 of file TauDecay_GenParticle.h.

Constructor & Destructor Documentation

TauDecay_GenParticle()

TauDecay_GenParticle::TauDecay_GenParticle

(

)

Definition at line 7 of file TauDecay_GenParticle.cc.

: TauDecay() {}

~TauDecay_GenParticle()

TauDecay_GenParticle::~TauDecay_GenParticle

(

)

Definition at line 9 of file TauDecay_GenParticle.cc.

{}

Member Function Documentation

AddPi0Info()

void TauDecay_GenParticle::AddPi0Info ( const reco::GenParticle *  Particle, unsigned int  midx  )

private

Definition at line 61 of file TauDecay_GenParticle.cc.

References mps_fire::i, MotherIdx, and TauDecayProducts.

Referenced by Analyze().

                                                                                        {
  if (Particle->status() == 1) {
    TauDecayProducts.push_back(Particle);
    MotherIdx.push_back(midx);
    return;
  }
  for (unsigned int i = 0; i < Particle->numberOfDaughters(); i++) {
    const reco::Candidate *dau = Particle->daughter(i);
    AddPi0Info(static_cast<const reco::GenParticle *>(dau), midx);
  }
}

AddRadInfo()

bool TauDecay_GenParticle::AddRadInfo	(	const reco::GenParticle *	Particle,
		unsigned int	midx
	)

Definition at line 73 of file TauDecay_GenParticle.cc.

References PdtPdgMini::gamma, mps_fire::i, MotherIdx, dqmiodumpmetadata::n, reco::Candidate::pdgId(), displacedMuons_cfi::photon, and TauDecayProducts.

Referenced by Analyze().

                                                                                        {
  bool same(false), photon(false);
  int n = 0;
  for (unsigned int i = 0; i < Particle->numberOfDaughters(); i++) {
    const reco::Candidate *dau = Particle->daughter(i);
    if (Particle->pdgId() == dau->pdgId()) {
      same = true;
    }
    if (dau->pdgId() == PdtPdgMini::gamma)
      photon = true;
    n++;
  }
  if (same && photon && n == 2) {
    for (unsigned int i = 0; i < Particle->numberOfDaughters(); i++) {
      const reco::Candidate *dau = Particle->daughter(i);
      if (Particle->pdgId() == dau->pdgId()) {
        if (!AddRadInfo(static_cast<const reco::GenParticle *>(dau), midx)) {
          TauDecayProducts.push_back(static_cast<const reco::GenParticle *>(dau));
          MotherIdx.push_back(midx);
        }
      }
      if (dau->pdgId() == PdtPdgMini::gamma) {
        TauDecayProducts.push_back(static_cast<const reco::GenParticle *>(dau));
        MotherIdx.push_back(midx);
      }
    }
    return true;
  }
  return false;
}

Analyze()

void TauDecay_GenParticle::Analyze ( const reco::GenParticle *  Particle, unsigned int  midx, bool  dores, bool  dopi0  )

private

Definition at line 30 of file TauDecay_GenParticle.cc.

References funct::abs(), AddPi0Info(), AddRadInfo(), gather_cfg::cout, mps_fire::i, TauDecay::isTauFinalStateParticle(), TauDecay::isTauParticleCounter(), TauDecay::isTauResonanceCounter(), MotherIdx, EgammaValidation_cff::pdgid, PdtPdgMini::pi0, and TauDecayProducts.

Referenced by AnalyzeTau().

                                                                                                             {
  unsigned int pdgid = abs(Particle->pdgId());
  if (isTauFinalStateParticle(pdgid)) {
    if (!isTauParticleCounter(pdgid))
      std::cout << "TauDecay_GenParticle::Analyze WARNING: Unknow Final State Particle in Tau Decay... " << std::endl;

    if (!AddRadInfo(static_cast<const reco::GenParticle *>(Particle), midx)) {
      TauDecayProducts.push_back(Particle);
      MotherIdx.push_back(midx);
    }
    // store information on pi0 decay products even though a pi0 is a finsal state particle (for 3PiPi0 studies)
    if (pdgid == PdtPdgMini::pi0 && dopi0) {
      midx = MotherIdx.size() - 1;
      for (unsigned int i = 0; i < Particle->numberOfDaughters(); i++) {
        const reco::Candidate *dau = Particle->daughter(i);
        AddPi0Info(static_cast<const reco::GenParticle *>(dau), midx);
      }
    }
    return;
  }
  if (Particle->status() == 1 || (isTauResonanceCounter(pdgid) && dores)) {
    TauDecayProducts.push_back(Particle);
    MotherIdx.push_back(midx);
    midx = MotherIdx.size() - 1;
  }
  for (unsigned int i = 0; i < Particle->numberOfDaughters(); i++) {
    const reco::Candidate *dau = Particle->daughter(i);
    Analyze(static_cast<const reco::GenParticle *>(dau), midx, dores, dopi0);
  }
}

AnalyzeTau()

bool TauDecay_GenParticle::AnalyzeTau	(	const reco::GenParticle *	Tau,
		unsigned int &	MODE_ID,
		unsigned int &	TauBitMask,
		bool	dores,
		bool	dopi0
	)

Definition at line 11 of file TauDecay_GenParticle.cc.

References funct::abs(), Analyze(), TauDecay::ClassifyDecayMode(), mps_fire::i, MODE_ID, MotherIdx, TauDecay::Reset(), PdtPdgMini::tau_minus, TauBitMask, and TauDecayProducts.

Referenced by TauValidation::analyze(), and TauValidation::spinEffectsZH().

                                                                                                         {
  Reset();
  MotherIdx.clear();
  TauDecayProducts.clear();
  if (abs(Tau->pdgId()) == PdtPdgMini::tau_minus) {  // check that it is a tau
    unsigned int Tauidx = TauDecayProducts.size();
    TauDecayProducts.push_back(Tau);
    MotherIdx.push_back(Tauidx);
    for (unsigned int i = 0; i < Tau->numberOfDaughters(); i++) {
      const reco::Candidate *dau = Tau->daughter(i);
      Analyze(static_cast<const reco::GenParticle *>(dau), Tauidx, dores, dopi0);
    }
    ClassifyDecayMode(MODE_ID, TauBitMask);
    return true;
  }
  return false;
}

Get_MotherIdx()

std::vector<unsigned int> TauDecay_GenParticle::Get_MotherIdx ( )

inline

Definition at line 32 of file TauDecay_GenParticle.h.

References MotherIdx.

{ return MotherIdx; }

Get_TauDecayProducts()

std::vector<const reco::GenParticle *> TauDecay_GenParticle::Get_TauDecayProducts ( )

inline

Definition at line 31 of file TauDecay_GenParticle.h.

References TauDecayProducts.

Referenced by TauValidation::analyze(), and TauValidation::spinEffectsZH().

{ return TauDecayProducts; }

Member Data Documentation

MODE_ID

unsigned int TauDecay_GenParticle::MODE_ID

private

Definition at line 42 of file TauDecay_GenParticle.h.

Referenced by AnalyzeTau().

MotherIdx

std::vector<unsigned int> TauDecay_GenParticle::MotherIdx

private

Definition at line 41 of file TauDecay_GenParticle.h.

Referenced by AddPi0Info(), AddRadInfo(), Analyze(), AnalyzeTau(), and Get_MotherIdx().

TauBitMask

unsigned int TauDecay_GenParticle::TauBitMask

private

Definition at line 42 of file TauDecay_GenParticle.h.

Referenced by AnalyzeTau().

TauDecayProducts

std::vector<const reco::GenParticle *> TauDecay_GenParticle::TauDecayProducts

private

Definition at line 40 of file TauDecay_GenParticle.h.

Referenced by AddPi0Info(), AddRadInfo(), Analyze(), AnalyzeTau(), and Get_TauDecayProducts().