TauDecay Class Reference

#include <TauDataFormat/TauNtuple/src/TauDecay_CMSSW.cc>

Inheritance diagram for TauDecay:

Public Types
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 }

Public Member Functions
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 ()

Static Public Member Functions
static std::string	DecayMode (unsigned int &MODE_ID)

Private Member Functions
void	ClassifyDecayResonance (unsigned int &TauBitMask)

Private Attributes
unsigned int	n_a1
unsigned int	n_a10
unsigned int	n_e
unsigned int	n_eta
unsigned int	n_gamma
unsigned int	n_K
unsigned int	n_K0L
unsigned int	n_K0S
unsigned int	n_Kstar
unsigned int	n_Kstar0
unsigned int	n_mu
unsigned int	n_nu
unsigned int	n_omega
unsigned int	n_pi
unsigned int	n_pi0
unsigned int	n_rho
unsigned int	n_rho0
unsigned int	unknown

Detailed Description

Description: This class reconstructs the MODE modes of the Tauola decays and provides a bit mask of the decay structure for the tau

Description: Bridge class for TauDecay when using HepMC

Implementation: [Notes on implementation]

Definition at line 21 of file TauDecay.h.

Member Enumeration Documentation

enum TauDecay::MODE

Enumerator
MODE_UNKNOWN
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

Definition at line 89 of file TauDecay.h.

            {
    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 TauDecay::TauDecayStructure

Enumerator
other
OneProng
ThreeProng
FiveProng
OnePi0
TwoPi0
ThreePi0
Res_a1_pm
Res_a1_0
Res_rho_pm
Res_rho_0
Res_eta
Res_omega
Res_Kstar_pm
Res_Kstar_0
KS0_to_pipi

Definition at line 131 of file TauDecay.h.

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

Constructor & Destructor Documentation

TauDecay::TauDecay

(

)

Definition at line 7 of file TauDecay.cc.

References Reset().

{ Reset(); }

TauDecay::~TauDecay

(

)

Definition at line 9 of file TauDecay.cc.

{}

Member Function Documentation

void TauDecay::ClassifyDecayMode	(	unsigned int &	MODE_ID,
		unsigned int &	TauBitMask
	)

Definition at line 150 of file TauDecay.cc.

References ClassifyDecayResonance(), gather_cfg::cout, FiveProng, MODE_2K0PIPI0, MODE_3PI, MODE_3PI2PI0, MODE_3PI3PI0, MODE_3PIETA, MODE_3PIOMEGA, MODE_3PIPI0, MODE_5PI, MODE_5PIPI0, MODE_ELECTRON, MODE_ETAPIPI0, MODE_K03PI, MODE_K0BK0PI, MODE_K0PI, MODE_K0PI2PI0, MODE_K2PI0, MODE_K2PIPI0, MODE_K3PI0, MODE_KAON, MODE_KETA, MODE_KK0B, MODE_KK0BPI0, MODE_KOMEGA, MODE_KPI0, MODE_KPIK, MODE_KPIKPI0, MODE_KPIPI, MODE_MUON, MODE_PI2PI0, MODE_PI2PI0ETA, MODE_PI2PI0OMEGA, MODE_PI3PI0, MODE_PI4PI0, MODE_PIK0PI0, MODE_PION, MODE_PIPI0, MODE_PIPI0GAM, MODE_UNKNOWN, n_e, n_eta, n_gamma, n_K, n_K0L, n_K0S, n_mu, n_nu, n_omega, n_pi, n_pi0, n_rho, OnePi0, OneProng, ThreePi0, ThreeProng, and TwoPi0.

Referenced by TauDecay_GenParticle::AnalyzeTau(), and TauDecay_CMSSW::AnalyzeTau().

                                                                                {
  //Reset Bits
  MODE_ID = 0;
  TauBitMask = 0;
  // Classify according to MODE and TauDecayStructure
  if (n_pi + n_K + n_e + n_mu == 1)
    TauBitMask = OneProng;
  if (n_pi + n_K == 3)
    TauBitMask = ThreeProng;
  if (n_pi + n_K == 5)
    TauBitMask = FiveProng;
  if (n_pi0 == 1)
    TauBitMask += OnePi0;
  if (n_pi0 == 2)
    TauBitMask += TwoPi0;
  if (n_pi0 == 3)
    TauBitMask += ThreePi0;
  ClassifyDecayResonance(TauBitMask);

  if (n_e == 0 && n_mu == 0 && n_pi == 1 && n_pi0 == 0 && n_K == 0 && n_K0S + n_K0L == 2 && n_nu == 1) {
    MODE_ID = MODE_K0BK0PI;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 1 && n_pi0 == 0 && n_K == 0 && (n_K0L + n_K0S) == 1 && n_nu == 1) {
    MODE_ID = MODE_K0PI;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 0 && n_pi0 == 1 && n_K == 1 && n_K0L == 0 && n_K0S == 0 && n_nu == 1) {
    MODE_ID = MODE_KPI0;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 1 && n_pi0 == 1 && n_K == 0 && n_K0L + n_K0S == 1 && n_nu == 1) {
    MODE_ID = MODE_PIK0PI0;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 0 && n_pi0 == 0 && n_K == 1 && (n_K0L + n_K0S) == 1 && n_nu == 1) {
    MODE_ID = MODE_KK0B;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 0 && n_pi0 == 1 && n_K == 1 && n_K0L + n_K0S == 1 && n_nu == 1) {
    MODE_ID = MODE_KK0BPI0;
    return;
  }
  if (n_e == 1 && n_mu == 0 && n_pi == 0 && n_pi0 == 0 && n_K == 0 && n_K0L == 0 && n_K0S == 0 && n_nu == 2) {
    MODE_ID = MODE_ELECTRON;
    return;
  }
  if (n_e == 0 && n_mu == 1 && n_pi == 0 && n_pi0 == 0 && n_K == 0 && n_K0L == 0 && n_K0S == 0 && n_nu == 2) {
    MODE_ID = MODE_MUON;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 1 && n_pi0 == 0 && n_K == 0 && n_K0L == 0 && n_K0S == 0 && n_nu == 1) {
    MODE_ID = MODE_PION;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 1 && n_pi0 == 1 && n_K == 0 && n_K0L == 0 && n_K0S == 0 && n_nu == 1) {
    MODE_ID = MODE_PIPI0;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 1 && n_pi0 == 2 && n_K == 0 && n_K0L == 0 && n_K0S == 0 && n_nu == 1) {
    MODE_ID = MODE_PI2PI0;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 3 && n_pi0 == 0 && n_K == 0 && n_K0L == 0 && n_K0S == 0 && n_nu == 1) {
    MODE_ID = MODE_3PI;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 0 && n_pi0 == 0 && n_K == 1 && n_K0L == 0 && n_K0S == 0 && n_nu == 1) {
    MODE_ID = MODE_KAON;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 3 && n_pi0 == 1 && n_K == 0 && n_K0L == 0 && n_K0S == 0 && n_nu == 1) {
    MODE_ID = MODE_3PIPI0;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 1 && n_pi0 == 3 && n_K == 0 && n_K0L == 0 && n_K0S == 0 && n_nu == 1) {
    MODE_ID = MODE_PI3PI0;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 3 && n_pi0 == 2 && n_K == 0 && n_K0L == 0 && n_K0S == 0 && n_nu == 1) {
    MODE_ID = MODE_3PI2PI0;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 5 && n_pi0 == 0 && n_K == 0 && n_K0L == 0 && n_K0S == 0 && n_nu == 1) {
    MODE_ID = MODE_5PI;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 5 && n_pi0 == 1 && n_K == 0 && n_K0L == 0 && n_K0S == 0 && n_nu == 1) {
    MODE_ID = MODE_5PIPI0;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 3 && n_pi0 == 3 && n_K == 0 && n_K0L == 0 && n_K0S == 0 && n_nu == 1) {
    MODE_ID = MODE_3PI3PI0;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 1 && n_pi0 == 0 && n_K == 2 && n_K0L == 0 && n_K0S == 0 && n_nu == 1) {
    MODE_ID = MODE_KPIK;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 0 && n_pi0 == 2 && n_K == 1 && n_K0L == 0 && n_K0S == 0 && n_nu == 1) {
    MODE_ID = MODE_K2PI0;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 2 && n_pi0 == 0 && n_K == 1 && n_K0L == 0 && n_K0S == 0 && n_nu == 1) {
    MODE_ID = MODE_KPIPI;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 1 && n_pi0 == 1 && n_K == 0 && n_K0L == 0 && n_K0S == 0 && n_nu == 1 &&
      n_gamma >= 1 && n_rho == 0) {
    MODE_ID = MODE_PIPI0GAM;
    return;  // Obsolete should not be called
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 1 && n_pi0 == 4 && n_K == 0 && n_K0L == 0 && n_K0S == 0 && n_nu == 1) {
    MODE_ID = MODE_PI4PI0;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 3 && n_pi0 == 0 && n_K == 0 && n_K0L == 0 && n_K0S == 0 && n_nu == 1 &&
      n_eta == 1) {
    MODE_ID = MODE_3PIETA;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 1 && n_pi0 == 2 && n_K == 0 && n_K0L == 0 && n_K0S == 0 && n_nu == 1 &&
      n_eta == 1) {
    MODE_ID = MODE_PI2PI0ETA;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 1 && n_pi0 == 2 && n_K == 0 && n_K0L == 0 && n_K0S == 0 && n_nu == 1 &&
      n_omega == 1) {
    MODE_ID = MODE_PI2PI0OMEGA;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 3 && n_pi0 == 0 && n_K == 0 && n_K0L == 0 && n_K0S == 0 && n_nu == 1 &&
      n_omega == 1) {
    MODE_ID = MODE_3PIOMEGA;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 0 && n_pi0 == 0 && n_K == 1 && n_K0L == 0 && n_K0S == 0 && n_nu == 1 &&
      n_omega == 1) {
    MODE_ID = MODE_KOMEGA;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 0 && n_pi0 == 3 && n_K == 1 && n_K0L == 0 && n_K0S == 0 && n_nu == 1) {
    MODE_ID = MODE_K3PI0;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 2 && n_pi0 == 1 && n_K == 1 && n_K0L == 0 && n_K0S == 0 && n_nu == 1) {
    MODE_ID = MODE_K2PIPI0;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 0 && n_pi0 == 0 && n_K == 1 && n_K0L == 0 && n_K0S == 0 && n_nu == 1 &&
      n_eta == 1) {
    MODE_ID = MODE_KETA;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 1 && n_pi0 == 2 && n_K == 0 && (n_K0L + n_K0S) == 0 && n_nu == 1) {
    MODE_ID = MODE_K0PI2PI0;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 3 && n_pi0 == 0 && n_K == 0 && (n_K0L + n_K0S) == 1 && n_nu == 1) {
    MODE_ID = MODE_K03PI;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 1 && n_pi0 == 1 && n_K == 0 && (n_K0L + n_K0S) == 2 && n_nu == 1) {
    MODE_ID = MODE_2K0PIPI0;
    return;
  }
  if (n_e == 0 && n_mu == 0 && n_pi == 1 && n_pi0 == 1 && n_K == 2 && n_K0L == 0 && n_K0S == 0 && n_nu == 1) {
    MODE_ID = MODE_KPIKPI0;
    return;
  }
  if (n_pi == 1 && n_pi0 == 1 && n_nu == 1 && n_eta == 1) {  // eta modes
    MODE_ID = MODE_ETAPIPI0;
    return;
  }

  std::cout << "Tau Mode not found: n_e " << n_e << " n_mu " << n_mu << " n_pi " << n_pi << " n_pi0 " << n_pi0
            << " n_K " << n_K << "  n_K0L " << n_K0L << "  n_K0S " << n_K0S << " n_nu  " << n_nu << " n_gamma "
            << n_gamma << std::endl;
  MODE_ID = MODE_UNKNOWN;
}

void TauDecay::ClassifyDecayResonance ( unsigned int &  TauBitMask )

private

Definition at line 373 of file TauDecay.cc.

References n_a1, n_a10, n_eta, n_Kstar, n_Kstar0, n_omega, n_rho, n_rho0, Res_a1_0, Res_a1_pm, Res_eta, Res_Kstar_0, Res_Kstar_pm, Res_omega, Res_rho_0, and Res_rho_pm.

Referenced by ClassifyDecayMode().

                                                              {
  // Add Resonance info to TauBitMask
  if (n_a1 > 0)
    TauBitMask += Res_a1_pm;
  if (n_a10 > 0)
    TauBitMask += Res_a1_0;
  if (n_rho > 0)
    TauBitMask += Res_rho_pm;
  if (n_rho0 > 0)
    TauBitMask += Res_rho_0;
  if (n_eta > 0)
    TauBitMask += Res_eta;
  if (n_omega > 0)
    TauBitMask += Res_omega;
  if (n_Kstar > 0)
    TauBitMask += Res_Kstar_pm;
  if (n_Kstar0 > 0)
    TauBitMask += Res_Kstar_0;
}

std::string TauDecay::DecayMode ( unsigned int &  MODE_ID )

static

Definition at line 393 of file TauDecay.cc.

References MODE_2K0PIPI0, MODE_3PI, MODE_3PI2PI0, MODE_3PI3PI0, MODE_3PIETA, MODE_3PIOMEGA, MODE_3PIPI0, MODE_5PI, MODE_5PIPI0, MODE_ELECTRON, MODE_ETAPIPI0, MODE_K03PI, MODE_K0BK0PI, MODE_K0PI, MODE_K0PI2PI0, MODE_K2PI0, MODE_K2PIPI0, MODE_K3PI0, MODE_KAON, MODE_KETA, MODE_KK0B, MODE_KK0BPI0, MODE_KOMEGA, MODE_KPI0, MODE_KPIK, MODE_KPIKPI0, MODE_KPIPI, MODE_MUON, MODE_PI2PI0, MODE_PI2PI0ETA, MODE_PI2PI0OMEGA, MODE_PI3PI0, MODE_PI4PI0, MODE_PIK0PI0, MODE_PION, MODE_PIPI0, and MODE_PIPI0GAM.

Referenced by TauValidation::bookHistograms().

                                                   {
  if (MODE_ID == MODE_ELECTRON)
    return "#tau^{#pm} #rightarrow e^{#pm}#nu#nu";
  else if (MODE_ID == MODE_MUON)
    return "#tau^{#pm} #rightarrow #mu^{#pm}#nu#nu";
  else if (MODE_ID == MODE_PION)
    return "#tau^{#pm} #rightarrow #pi^{#pm}#nu";
  else if (MODE_ID == MODE_PIPI0)
    return "#tau^{#pm} #rightarrow #pi^{#pm}#pi^{0}#nu";
  else if (MODE_ID == MODE_3PI)
    return "#tau^{#pm} #rightarrow #pi^{#pm}#pi^{#pm}#pi^{#mp}#nu";
  else if (MODE_ID == MODE_PI2PI0)
    return "#tau^{#pm} #rightarrow #pi^{#pm}#pi^{0}#pi^{0}#nu";
  else if (MODE_ID == MODE_KAON)
    return "#tau^{#pm} #rightarrow K^{#pm}#nu";
  else if (MODE_ID == MODE_KPI0)
    return "#tau^{#pm} #rightarrow K^{#pm}#pi^{0}#nu";
  else if (MODE_ID == MODE_K0PI)
    return "#tau^{#pm} #rightarrow K^{0}#pi^{#pm}#nu";
  else if (MODE_ID == MODE_3PIPI0)
    return "#tau^{#pm} #rightarrow #pi^{#pm}#pi^{#pm}#pi^{#mp}#pi^{0}#nu";
  else if (MODE_ID == MODE_PI3PI0)
    return "#tau^{#pm} #rightarrow #pi^{#pm}#pi^{0}#pi^{0}#pi^{0}#nu";
  else if (MODE_ID == MODE_3PI2PI0)
    return "#tau^{#pm} #rightarrow #pi^{#pm}#pi^{#pm}#pi^{#mp}#pi^{0}#pi^{0}#nu";
  else if (MODE_ID == MODE_5PI)
    return "#tau^{#pm} #rightarrow #pi^{#pm}#pi^{#pm}#pi^{#pm}#pi^{#mp}#pi^{#mp}#nu";
  else if (MODE_ID == MODE_5PIPI0)
    return "#tau^{#pm} #rightarrow #pi^{#pm}#pi^{#pm}#pi^{#pm}#pi^{#mp}#pi^{#mp}#pi^{0}#nu";
  else if (MODE_ID == MODE_3PI3PI0)
    return "#tau^{#pm} #rightarrow #pi^{#pm}#pi^{#pm}#pi^{#mp}#pi^{0}#pi^{0}#pi^{0}#nu";
  else if (MODE_ID == MODE_KPIK)
    return "#tau^{#pm} #rightarrow K^{#pm}#pi^{#pm}K^{#mp}#nu";
  else if (MODE_ID == MODE_K0BK0PI)
    return "#tau^{#pm} #rightarrow #bar{K}^{0}K^{0}#pi^{#pm}#nu";
  else if (MODE_ID == MODE_KK0BPI0)
    return "#tau^{#pm} #rightarrow #bar{K}^{0}K^{0}#pi^{#pm}#pi^{0}#nu";
  else if (MODE_ID == MODE_K2PI0)
    return "#tau^{#pm} #rightarrow K^{#pm}#pi^{0}#pi^{0}#nu";
  else if (MODE_ID == MODE_KPIPI)
    return "#tau^{#pm} #rightarrow K^{#pm}#pi^{#pm}#pi^{#mp}#nu";
  else if (MODE_ID == MODE_PIK0PI0)
    return "#tau^{#pm} #rightarrow K^{0}#pi^{#pm}#pi^{0}#nu";
  else if (MODE_ID == MODE_ETAPIPI0)
    return "#tau^{#pm} #rightarrow #pi^{#pm}#eta#pi^{0}#nu";
  else if (MODE_ID == MODE_PIPI0GAM)
    return "#tau^{#pm} #rightarrow #pi^{#pm}#pi^{0}#nu#gamma (obsolete)";
  else if (MODE_ID == MODE_KK0B)
    return "#tau^{#pm} #rightarrow K^{#pm}#bar{K}^{0}#nu";
  else if (MODE_ID == MODE_PI4PI0)
    return "#tau^{#pm} #rightarrow #pi^{#pm}#pi^{0}#pi^{0}#pi^{0}#pi^{0}#nu";
  else if (MODE_ID == MODE_3PIETA)
    return "#tau^{#pm} #rightarrow #pi^{#pm}#pi^{#pm}#pi^{#mp}#eta#nu";
  else if (MODE_ID == MODE_PI2PI0ETA)
    return "#tau^{#pm} #rightarrow #pi^{#pm}#pi^{0}#pi^{0}#eta#nu";
  else if (MODE_ID == MODE_PI2PI0OMEGA)
    return "#tau^{#pm} #rightarrow #pi^{#pm}#pi^{0}#pi^{0}#omega#nu";
  else if (MODE_ID == MODE_3PIOMEGA)
    return "#tau^{#pm} #rightarrow #pi^{#pm}#pi^{#pm}#pi^{#mp}#omega#nu";
  else if (MODE_ID == MODE_KOMEGA)
    return "#tau^{#pm} #rightarrow K^{#pm}#omega#nu";
  else if (MODE_ID == MODE_K3PI0)
    return "#tau^{#pm} #rightarrow K#pi^{0}#pi^{0}#pi^{0}#nu";
  else if (MODE_ID == MODE_K2PIPI0)
    return "#tau^{#pm} #rightarrow K^{#pm}#pi^{#pm}#pi^{#mp}#pi^{0}#nu";
  else if (MODE_ID == MODE_KETA)
    return "#tau^{#pm} #rightarrow K^{#pm}#eta#nu";
  else if (MODE_ID == MODE_K0PI2PI0)
    return "#tau^{#pm} #rightarrow K^{0}#pi^{#pm}#pi^{0}#pi^{0}#nu";
  else if (MODE_ID == MODE_K03PI)
    return "#tau^{#pm} #rightarrow K^{0}#pi^{#pm}#pi^{#pm}#pi^{#mp}#nu";
  else if (MODE_ID == MODE_2K0PIPI0)
    return "#tau^{#pm} #rightarrow K^{0}#bar{K}^{0}#pi^{#pm}#pi^{-}#nu";
  else if (MODE_ID == MODE_KPIKPI0)
    return "#tau^{#pm} #rightarrow K^{#pm}#pi^{#pm}K^{#mp}#pi^{0}#nu";
  return "UnKnown";
}

bool TauDecay::hasResonance	(	unsigned int &	TauBitMask,
		int	pdgid
	)

Definition at line 350 of file TauDecay.cc.

References PdtPdgMini::a_10, PdtPdgMini::a_1_plus, funct::abs(), PdtPdgMini::eta, PdtPdgMini::K_S0, PdtPdgMini::K_star0, PdtPdgMini::K_star_plus, KS0_to_pipi, PdtPdgMini::omega, AlCaHLTBitMon_ParallelJobs::p, Res_a1_0, Res_a1_pm, Res_eta, Res_Kstar_0, Res_Kstar_pm, Res_omega, Res_rho_0, Res_rho_pm, PdtPdgMini::rho0, and PdtPdgMini::rho_plus.

                                                               {
  int p = abs(pdgid);
  if (p == PdtPdgMini::a_1_plus && Res_a1_pm & TauBitMask)
    return true;
  if (p == PdtPdgMini::a_10 && Res_a1_0 & TauBitMask)
    return true;
  if (p == PdtPdgMini::rho_plus && Res_rho_pm & TauBitMask)
    return true;
  if (p == PdtPdgMini::rho0 && Res_rho_0 & TauBitMask)
    return true;
  if (p == PdtPdgMini::eta && Res_eta & TauBitMask)
    return true;
  if (p == PdtPdgMini::omega && Res_omega & TauBitMask)
    return true;
  if (p == PdtPdgMini::K_star0 && Res_Kstar_pm & TauBitMask)
    return true;
  if (p == PdtPdgMini::K_star_plus && Res_Kstar_0 & TauBitMask)
    return true;
  if (p == PdtPdgMini::K_S0 && KS0_to_pipi & TauBitMask)
    return true;
  return false;
}

bool TauDecay::isTauFinalStateParticle

(

int

pdgid

)

Definition at line 32 of file TauDecay.cc.

References funct::abs(), PdtPdgMini::e_minus, PdtPdgMini::eta, PdtPdgMini::gamma, PdtPdgMini::K_L0, PdtPdgMini::K_plus, PdtPdgMini::K_S0, PdtPdgMini::mu_minus, PdtPdgMini::nu_e, PdtPdgMini::nu_mu, PdtPdgMini::nu_tau, PdtPdgMini::omega, PdtPdgMini::pi0, and PdtPdgMini::pi_plus.

Referenced by TauValidation::analyze(), TauDecay_GenParticle::Analyze(), and TauDecay_CMSSW::Analyze().

                                                {
  int id = abs(pdgid);
  if (id == PdtPdgMini::e_minus)
    return true;  // e+-
  if (id == PdtPdgMini::nu_e)
    return true;  // nu_e
  if (id == PdtPdgMini::mu_minus)
    return true;  // mu+-
  if (id == PdtPdgMini::nu_mu)
    return true;  // nu_mu
  if (id == PdtPdgMini::nu_tau)
    return true;  // nu_tau
  if (id == PdtPdgMini::gamma)
    return true;  // gamma happends in generator
  if (id == PdtPdgMini::pi0)
    return true;  // pi0
  if (id == PdtPdgMini::pi_plus)
    return true;  // pi+-
  if (id == PdtPdgMini::K_L0)
    return true;  // K0L
  if (id == PdtPdgMini::K_S0)
    return true;  // KS
  if (id == PdtPdgMini::eta)
    return true;
  if (id == PdtPdgMini::omega)
    return true;

  if (id == PdtPdgMini::K_plus)
    return true;  // K+-
  return false;
}

bool TauDecay::isTauParticleCounter

(

int

pdgid

)

Definition at line 64 of file TauDecay.cc.

References funct::abs(), gather_cfg::cout, PdtPdgMini::e_minus, PdtPdgMini::eta, PdtPdgMini::gamma, PdtPdgMini::K0, PdtPdgMini::K_L0, PdtPdgMini::K_plus, PdtPdgMini::K_S0, PdtPdgMini::mu_minus, n_e, n_eta, n_gamma, n_K, n_K0L, n_K0S, n_mu, n_nu, n_omega, n_pi, n_pi0, PdtPdgMini::nu_e, PdtPdgMini::nu_mu, PdtPdgMini::nu_tau, PdtPdgMini::omega, PdtPdgMini::pi0, and PdtPdgMini::pi_plus.

Referenced by TauDecay_GenParticle::Analyze(), and TauDecay_CMSSW::Analyze().

                                             {
  int id = abs(pdgid);
  //count particles
  if (id == PdtPdgMini::pi_plus) {
    n_pi++;
    return true;
  }
  if (id == PdtPdgMini::pi0) {
    n_pi0++;
    return true;
  }
  if (id == PdtPdgMini::K_plus) {
    n_K++;
    return true;
  }
  if (id == PdtPdgMini::K_L0) {
    n_K0L++;
    return true;
  }
  if (id == PdtPdgMini::K_S0) {
    n_K0S++;
    return true;
  }
  if (id == PdtPdgMini::eta) {
    n_eta++;
    return true;
  }
  if (id == PdtPdgMini::omega) {
    n_omega++;
    return true;
  }
  if (id == PdtPdgMini::gamma) {
    n_gamma++;
    return true;
  }
  if (id == PdtPdgMini::nu_tau || id == PdtPdgMini::nu_e || id == PdtPdgMini::nu_mu) {
    n_nu++;
    return true;
  }
  if (id == PdtPdgMini::e_minus) {
    n_e++;
    return true;
  }
  if (id == PdtPdgMini::mu_minus) {
    n_mu++;
    return true;
  }
  if (abs(id) == PdtPdgMini::K0)
    std::cout << "TauDecay::isTauParticleCounter: ERROR unidentified Particle: " << id << std::endl;
  return false;
}

bool TauDecay::isTauResonanceCounter

(

int

pdgid

)

Definition at line 116 of file TauDecay.cc.

References PdtPdgMini::a_10, PdtPdgMini::a_1_plus, funct::abs(), PdtPdgMini::K_star0, PdtPdgMini::K_star_plus, n_a1, n_a10, n_Kstar, n_Kstar0, n_rho, n_rho0, PdtPdgMini::rho0, PdtPdgMini::rho_plus, unknown, and PdtPdgMini::W_plus.

Referenced by TauDecay_GenParticle::Analyze(), and TauDecay_CMSSW::Analyze().

                                              {
  int id = abs(pdgid);
  //count resonances
  if (id == PdtPdgMini::a_1_plus) {
    n_a1++;
    return true;
  }
  if (id == PdtPdgMini::a_10) {
    n_a10++;
    return true;
  }
  if (id == PdtPdgMini::rho_plus) {
    n_rho++;
    return true;
  }
  if (id == PdtPdgMini::rho0) {
    n_rho0++;
    return true;
  }
  if (id == PdtPdgMini::K_star0) {
    n_Kstar0++;
    return true;
  }
  if (id == PdtPdgMini::K_star_plus) {
    n_Kstar++;
    return true;
  }
  if (id == PdtPdgMini::W_plus) {
    return true;
  }
  unknown++;
  return false;
}

unsigned int TauDecay::nPi0

(

unsigned int &

TauBitMask

)

Definition at line 340 of file TauDecay.cc.

References OnePi0, ThreePi0, and TwoPi0.

                                                    {
  if (OnePi0 & TauBitMask)
    return 1;
  if (TwoPi0 & TauBitMask)
    return 2;
  if (ThreePi0 & TauBitMask)
    return 3;
  return 0;
}

unsigned int TauDecay::nProng

(

unsigned int &

TauBitMask

)

Definition at line 331 of file TauDecay.cc.

References FiveProng, OneProng, and ThreeProng.

Referenced by TauValidation::analyze().

                                                      {
  if (OneProng & TauBitMask)
    return 1;
  if (ThreeProng & TauBitMask)
    return 3;
  if (FiveProng & TauBitMask)
    return 5;
  return 7;
}

void TauDecay::Reset

(

)

Definition at line 11 of file TauDecay.cc.

References n_a1, n_a10, n_e, n_eta, n_gamma, n_K, n_K0L, n_K0S, n_Kstar, n_Kstar0, n_mu, n_nu, n_omega, n_pi, n_pi0, n_rho, n_rho0, and unknown.

Referenced by TauDecay_GenParticle::AnalyzeTau(), TauDecay_CMSSW::AnalyzeTau(), and TauDecay().

                     {
  n_pi = 0;
  n_pi0 = 0;
  n_K = 0;
  n_K0L = 0;
  n_K0S = 0;
  n_gamma = 0;
  n_nu = 0;
  n_e = 0;
  n_mu = 0;
  n_a1 = 0;
  n_a10 = 0;
  n_rho = 0;
  n_rho0 = 0;
  n_eta = 0;
  n_omega = 0;
  n_Kstar0 = 0;
  n_Kstar = 0;
  unknown = 0;
}

Member Data Documentation

unsigned int TauDecay::n_a1

private

Definition at line 169 of file TauDecay.h.

Referenced by ClassifyDecayResonance(), isTauResonanceCounter(), and Reset().

unsigned int TauDecay::n_a10

private

Definition at line 169 of file TauDecay.h.

Referenced by ClassifyDecayResonance(), isTauResonanceCounter(), and Reset().

unsigned int TauDecay::n_e

private

Definition at line 168 of file TauDecay.h.

Referenced by ClassifyDecayMode(), isTauParticleCounter(), and Reset().

unsigned int TauDecay::n_eta

private

Definition at line 169 of file TauDecay.h.

Referenced by ClassifyDecayMode(), ClassifyDecayResonance(), isTauParticleCounter(), and Reset().

unsigned int TauDecay::n_gamma

private

Definition at line 168 of file TauDecay.h.

Referenced by ClassifyDecayMode(), isTauParticleCounter(), and Reset().

unsigned int TauDecay::n_K

private

Definition at line 168 of file TauDecay.h.

Referenced by ClassifyDecayMode(), isTauParticleCounter(), and Reset().

unsigned int TauDecay::n_K0L

private

Definition at line 168 of file TauDecay.h.

Referenced by ClassifyDecayMode(), isTauParticleCounter(), and Reset().

unsigned int TauDecay::n_K0S

private

Definition at line 168 of file TauDecay.h.

Referenced by ClassifyDecayMode(), isTauParticleCounter(), and Reset().

unsigned int TauDecay::n_Kstar

private

Definition at line 169 of file TauDecay.h.

Referenced by ClassifyDecayResonance(), isTauResonanceCounter(), and Reset().

unsigned int TauDecay::n_Kstar0

private

Definition at line 169 of file TauDecay.h.

Referenced by ClassifyDecayResonance(), isTauResonanceCounter(), and Reset().

unsigned int TauDecay::n_mu

private

Definition at line 168 of file TauDecay.h.

Referenced by ClassifyDecayMode(), isTauParticleCounter(), and Reset().

unsigned int TauDecay::n_nu

private

Definition at line 168 of file TauDecay.h.

Referenced by ClassifyDecayMode(), isTauParticleCounter(), and Reset().

unsigned int TauDecay::n_omega

private

Definition at line 169 of file TauDecay.h.

Referenced by ClassifyDecayMode(), ClassifyDecayResonance(), isTauParticleCounter(), and Reset().

unsigned int TauDecay::n_pi

private

Definition at line 168 of file TauDecay.h.

Referenced by ClassifyDecayMode(), isTauParticleCounter(), and Reset().

unsigned int TauDecay::n_pi0

private

Definition at line 168 of file TauDecay.h.

Referenced by ClassifyDecayMode(), isTauParticleCounter(), and Reset().

unsigned int TauDecay::n_rho

private

Definition at line 169 of file TauDecay.h.

Referenced by ClassifyDecayMode(), ClassifyDecayResonance(), isTauResonanceCounter(), and Reset().

unsigned int TauDecay::n_rho0

private

Definition at line 169 of file TauDecay.h.

Referenced by ClassifyDecayResonance(), isTauResonanceCounter(), and Reset().

unsigned int TauDecay::unknown

private

Definition at line 169 of file TauDecay.h.

Referenced by isTauResonanceCounter(), and Reset().