#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

◆ MODE

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

◆ TauDecayStructure

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::TauDecay ( )

Definition at line 7 of file TauDecay.cc.

References Reset().

7 { Reset(); }

TauDecay::Reset

void Reset()

Definition: TauDecay.cc:11

◆ ~TauDecay()

TauDecay::~TauDecay ( )

Definition at line 9 of file TauDecay.cc.

9 {}

Member Function Documentation

◆ ClassifyDecayMode()

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

◆ ClassifyDecayResonance()

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

◆ DecayMode()

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

◆ hasResonance()

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, EgammaValidation_cff::pdgid, 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;
 }

◆ isTauFinalStateParticle()

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, EgammaValidation_cff::pdgid, 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;
 }

◆ isTauParticleCounter()

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, EgammaValidation_cff::pdgid, 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;
 }

◆ isTauResonanceCounter()

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, EgammaValidation_cff::pdgid, 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;
 }

◆ nPi0()

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

◆ nProng()

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

◆ Reset()

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

◆ n_a1

unsigned int TauDecay::n_a1

private

Definition at line 169 of file TauDecay.h.

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

◆ n_a10

unsigned int TauDecay::n_a10

private

Definition at line 169 of file TauDecay.h.

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

◆ n_e

unsigned int TauDecay::n_e

private

Definition at line 168 of file TauDecay.h.

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

◆ n_eta

unsigned int TauDecay::n_eta

private

Definition at line 169 of file TauDecay.h.

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

◆ n_gamma

unsigned int TauDecay::n_gamma

private

Definition at line 168 of file TauDecay.h.

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

◆ n_K

unsigned int TauDecay::n_K

private

Definition at line 168 of file TauDecay.h.

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

◆ n_K0L

unsigned int TauDecay::n_K0L

private

Definition at line 168 of file TauDecay.h.

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

◆ n_K0S

unsigned int TauDecay::n_K0S

private

Definition at line 168 of file TauDecay.h.

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

◆ n_Kstar

unsigned int TauDecay::n_Kstar

private

Definition at line 169 of file TauDecay.h.

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

◆ n_Kstar0

unsigned int TauDecay::n_Kstar0

private

Definition at line 169 of file TauDecay.h.

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

◆ n_mu

unsigned int TauDecay::n_mu

private

Definition at line 168 of file TauDecay.h.

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

◆ n_nu

unsigned int TauDecay::n_nu

private

Definition at line 168 of file TauDecay.h.

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

◆ n_omega

unsigned int TauDecay::n_omega

private

Definition at line 169 of file TauDecay.h.

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

◆ n_pi

unsigned int TauDecay::n_pi

private

Definition at line 168 of file TauDecay.h.

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

◆ n_pi0

unsigned int TauDecay::n_pi0

private

Definition at line 168 of file TauDecay.h.

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

◆ n_rho

unsigned int TauDecay::n_rho

private

Definition at line 169 of file TauDecay.h.

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

◆ n_rho0

unsigned int TauDecay::n_rho0

private

Definition at line 169 of file TauDecay.h.

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

◆ unknown

unsigned int TauDecay::unknown

private

Definition at line 169 of file TauDecay.h.

Referenced by isTauResonanceCounter(), and Reset().

Public Types

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Enumeration Documentation

◆ MODE

◆ TauDecayStructure

Constructor & Destructor Documentation

◆ TauDecay()

◆ ~TauDecay()

Member Function Documentation

◆ ClassifyDecayMode()

◆ ClassifyDecayResonance()

◆ DecayMode()

◆ hasResonance()

◆ isTauFinalStateParticle()

◆ isTauParticleCounter()

◆ isTauResonanceCounter()

◆ nPi0()

◆ nProng()

◆ Reset()

Member Data Documentation

◆ 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

◆ unknown