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

Inheritance diagram for TauDecay:

Public Types
enum	JAK { JAK_UNKNOWN = 0, JAK_ELECTRON = 1, JAK_MUON = 2, JAK_PION = 3, JAK_RHO_PIPI0 = 4, JAK_A1_3PI = 5, JAK_KAON = 6, JAK_KSTAR = 7, JAK_3PIPI0 = 8, JAK_PI3PI0 = 9, JAK_3PI2PI0 = 10, JAK_5PI = 11, JAK_5PIPI0 = 12, JAK_3PI3PI0 = 13, JAK_KPIK = 14, JAK_K0BK0PI = 15, JAK_KK0BPI0 = 16, JAK_K2PI0 = 17, JAK_KPIPI = 18, JAK_PIK0PI0 = 19, JAK_ETAPIPI0 = 20, JAK_PIPI0GAM = 21, JAK_KK0B = 22, NJAKID = 23 }
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 &JAK_ID, unsigned int &TauBitMask)
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	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 JAK modes of the Tauola decays and provides a bit mask of the decay structure for the tau

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 23 of file TauDecay.h.

Member Enumeration Documentation

enum TauDecay::JAK

Enumerator:

JAK_UNKNOWN
JAK_ELECTRON
JAK_MUON
JAK_PION
JAK_RHO_PIPI0
JAK_A1_3PI
JAK_KAON
JAK_KSTAR
JAK_3PIPI0
JAK_PI3PI0
JAK_3PI2PI0
JAK_5PI
JAK_5PIPI0
JAK_3PI3PI0
JAK_KPIK
JAK_K0BK0PI
JAK_KK0BPI0
JAK_K2PI0
JAK_KPIPI
JAK_PIK0PI0
JAK_ETAPIPI0
JAK_PIPI0GAM
JAK_KK0B
NJAKID

Definition at line 51 of file TauDecay.h.

          {JAK_UNKNOWN=0,
           JAK_ELECTRON=1,
           JAK_MUON=2,
           JAK_PION=3,
           JAK_RHO_PIPI0=4,
           JAK_A1_3PI=5,
           JAK_KAON=6,
           JAK_KSTAR=7,
           JAK_3PIPI0=8,
           JAK_PI3PI0=9,
           JAK_3PI2PI0=10,
           JAK_5PI=11,
           JAK_5PIPI0=12,
           JAK_3PI3PI0=13,
           JAK_KPIK=14,
           JAK_K0BK0PI=15,
           JAK_KK0BPI0=16,
           JAK_K2PI0=17,
           JAK_KPIPI=18,
           JAK_PIK0PI0=19,
           JAK_ETAPIPI0=20,
           JAK_PIPI0GAM=21,
           JAK_KK0B=22,
           NJAKID=23
  };

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

References Reset().

                  {
  Reset();
}

TauDecay::~TauDecay ( )

Definition at line 13 of file TauDecay.cc.

                   {

}

Member Function Documentation

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

Definition at line 88 of file TauDecay.cc.

References ClassifyDecayResonance(), gather_cfg::cout, FiveProng, JAK_3PI2PI0, JAK_3PI3PI0, JAK_3PIPI0, JAK_5PI, JAK_5PIPI0, JAK_A1_3PI, JAK_ELECTRON, JAK_ETAPIPI0, JAK_K0BK0PI, JAK_K2PI0, JAK_KAON, JAK_KK0B, JAK_KK0BPI0, JAK_KPIK, JAK_KPIPI, JAK_KSTAR, JAK_MUON, JAK_PI3PI0, JAK_PIK0PI0, JAK_PION, JAK_RHO_PIPI0, JAK_UNKNOWN, KS0_to_pipi, n_e, n_eta, n_gamma, n_K, n_K0L, n_K0S, n_mu, n_nu, n_pi, n_pi0, OnePi0, OneProng, ThreePi0, ThreeProng, and TwoPi0.

Referenced by TauDecay_CMSSW::AnalyzeTau().

                                                                             {
  //Reset Bits
  JAK_ID=0;
  TauBitMask=0;
  // Classify according to JAK and TauDecayStructure
  //
  // Exlusive modes remove first
  //
  if(n_pi>=1 && n_pi0>=1 && n_nu==1 && n_eta==1){ // eta modes
    JAK_ID=JAK_ETAPIPI0;
    TauBitMask=OneProng;
    if(n_pi0==1)TauBitMask+=OnePi0;
    if(n_pi0==2)TauBitMask+=TwoPi0;
    if(n_pi0==3)TauBitMask+=ThreePi0;
    ClassifyDecayResonance(TauBitMask);
    return;
  }
  //JAK_K0BK0PI
  if ((n_K0S+n_K0L==2 && n_pi==1 && n_pi0==0 && n_K==0 && n_nu==1) || (n_pi>=1 && n_pi0==0 && n_K0L+n_K0S==2 && n_K==0 && n_nu==1) ){
    JAK_ID=JAK_K0BK0PI;
    TauBitMask=OneProng;
    if(n_pi>1)TauBitMask+=KS0_to_pipi;
    ClassifyDecayResonance(TauBitMask);
    return;
  }
  //if(n_Kstar==1){ //JAK_KSTAR K0SPi && KPi0 
    if(n_e==0 && n_mu==0 && n_pi>=1 && n_pi0==0 && n_K==0 && n_K0L+n_K0S==1 && n_nu==1){
      JAK_ID=JAK_KSTAR;
      TauBitMask=OneProng;
      if(n_pi==3)TauBitMask+=KS0_to_pipi;
      ClassifyDecayResonance(TauBitMask);
      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){
      JAK_ID=JAK_KSTAR;
      TauBitMask=OneProng;
      ClassifyDecayResonance(TauBitMask);
      return;
    }
    //}
  //JAK_PIK0PI0
  if(n_e==0 && n_mu==0 && n_pi>=1 && n_pi0==1 && n_K==0 && n_K0L+n_K0S==1  && n_nu==1){
    JAK_ID=JAK_PIK0PI0;
    TauBitMask=OneProng;
    TauBitMask+=OnePi0;
    if(n_pi==3)TauBitMask+=KS0_to_pipi;
    ClassifyDecayResonance(TauBitMask);
    return;
  }
  //JAK_KK0B
  if(n_e==0 && n_mu==0 && n_pi0==0 && n_K==1 && n_K0L+n_K0S==1 && n_nu==1){
    JAK_ID=JAK_KK0B;
    TauBitMask=OneProng;
    if(n_pi==2)TauBitMask+=KS0_to_pipi;
    return;
  }
  //JAK_ID=JAK_KK0BPI0
  if(n_e==0 && n_mu==0 && n_pi0==1 && n_K==1 && n_K0L+n_K0S==1  && n_nu==1){
    JAK_ID=JAK_KK0BPI0;
    TauBitMask=OneProng;
    TauBitMask+=OnePi0;
    if(n_pi==2)TauBitMask+=KS0_to_pipi;
    ClassifyDecayResonance(TauBitMask);
    return;
  }
  
  
  //Safty handelling for exlusive modes
  if (n_K0L!=0){
    std::cout << "Unknown mode with KL0: 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;
    return;
  }
  if (n_K0S!=0){
    std::cout << "Unknown mode with KS0: 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;
    return;
  }
  if (n_eta!=0){
    std::cout << "Unknown mode with eta: 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;
    return;
  }



  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);
  //
  // Standard modes
  //
  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){
    JAK_ID=JAK_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){
    JAK_ID=JAK_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){
    JAK_ID=JAK_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){// && n_rho==1){ removing intermediate resoance to be compatible with pythia8
    JAK_ID=JAK_RHO_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){
    JAK_ID=JAK_A1_3PI;
    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){
    JAK_ID=JAK_A1_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){
    JAK_ID=JAK_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){
    JAK_ID=JAK_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){
    JAK_ID=JAK_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){
    JAK_ID=JAK_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){
    JAK_ID=JAK_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){
    JAK_ID=JAK_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){
    JAK_ID=JAK_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){
    JAK_ID=JAK_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){
    JAK_ID=JAK_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){
    JAK_ID=JAK_KPIPI;
    return;
  }
  // removing JAKID 21 to allow for compatibility with Pythia8 
  /*  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){
    JAK_ID=JAK_PIPI0GAM;
    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;
  JAK_ID=JAK_UNKNOWN;
}

void TauDecay::ClassifyDecayResonance ( unsigned int & TauBitMask ) [private]

Definition at line 259 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;
}

bool TauDecay::isTauFinalStateParticle ( int pdgid )

Definition at line 38 of file TauDecay.cc.

References abs, PdtPdgMini::e_minus, PdtPdgMini::gamma, PdtPdgMini::K_L0, PdtPdgMini::K_plus, PdtPdgMini::K_S0, PdtPdgMini::mu_minus, PdtPdgMini::nu_e, PdtPdgMini::nu_mu, PdtPdgMini::nu_tau, PdtPdgMini::pi0, and PdtPdgMini::pi_plus.

Referenced by 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::K_plus)    return true;  // K+-
  return false;
}

bool TauDecay::isTauParticleCounter ( int pdgid )

Definition at line 54 of file TauDecay.cc.

References abs, PdtPdgMini::e_minus, PdtPdgMini::gamma, PdtPdgMini::K_L0, PdtPdgMini::K_plus, PdtPdgMini::K_S0, PdtPdgMini::mu_minus, n_e, n_gamma, n_K, n_K0L, n_K0S, n_mu, n_nu, n_pi, n_pi0, PdtPdgMini::nu_e, PdtPdgMini::nu_mu, PdtPdgMini::nu_tau, PdtPdgMini::pi0, and PdtPdgMini::pi_plus.

Referenced by 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::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;}
  return false;
}

bool TauDecay::isTauResonanceCounter ( int pdgid )

Definition at line 71 of file TauDecay.cc.

References PdtPdgMini::a_10, PdtPdgMini::a_1_plus, abs, PdtPdgMini::eta, PdtPdgMini::K_star0, PdtPdgMini::K_star_plus, n_a1, n_a10, n_eta, n_Kstar, n_Kstar0, n_omega, n_rho, n_rho0, PdtPdgMini::omega, PdtPdgMini::rho0, PdtPdgMini::rho_plus, unknown, and PdtPdgMini::W_plus.

Referenced by 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::eta)        { n_eta++;     return true;}
  if(id==PdtPdgMini::omega)      { n_omega++;   return true;}
  //if(id==PdtPdgMini::K_S0)       { n_K0S++;     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 ) [inline]

Definition at line 110 of file TauDecay.h.

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 ) [inline]

Definition at line 104 of file TauDecay.h.

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