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 11 of file TauDecay.cc.

                   {

}

Member Function Documentation

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

Definition at line 89 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, mathSSE::return(), 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 246 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 258 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 231 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 36 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 TauDecay_GenParticle::Analyze(), TauDecay_CMSSW::Analyze(), and TauValidation::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 55 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 75 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 224 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 218 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

(

void

)

Definition at line 15 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 168 of file TauDecay.h.

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

unsigned int TauDecay::n_a10

private

Definition at line 168 of file TauDecay.h.

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

unsigned int TauDecay::n_e

private

Definition at line 167 of file TauDecay.h.

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

unsigned int TauDecay::n_eta

private

Definition at line 168 of file TauDecay.h.

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

unsigned int TauDecay::n_gamma

private

Definition at line 167 of file TauDecay.h.

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

unsigned int TauDecay::n_K

private

Definition at line 167 of file TauDecay.h.

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

unsigned int TauDecay::n_K0L

private

Definition at line 167 of file TauDecay.h.

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

unsigned int TauDecay::n_K0S

private

Definition at line 167 of file TauDecay.h.

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

unsigned int TauDecay::n_Kstar

private

Definition at line 168 of file TauDecay.h.

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

unsigned int TauDecay::n_Kstar0

private

Definition at line 168 of file TauDecay.h.

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

unsigned int TauDecay::n_mu

private

Definition at line 167 of file TauDecay.h.

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

unsigned int TauDecay::n_nu

private

Definition at line 167 of file TauDecay.h.

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

unsigned int TauDecay::n_omega

private

Definition at line 168 of file TauDecay.h.

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

unsigned int TauDecay::n_pi

private

Definition at line 167 of file TauDecay.h.

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

unsigned int TauDecay::n_pi0

private

Definition at line 167 of file TauDecay.h.

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

unsigned int TauDecay::n_rho

private

Definition at line 168 of file TauDecay.h.

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

unsigned int TauDecay::n_rho0

private

Definition at line 168 of file TauDecay.h.

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

unsigned int TauDecay::unknown

private

Definition at line 168 of file TauDecay.h.

Referenced by isTauResonanceCounter(), and Reset().