#include <TauDecay_GenParticle.h>

Inheritance diagram for TauDecay_GenParticle:

Public Member Functions
bool	AddRadInfo (const reco::GenParticle *Particle, unsigned int midx)

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

std::vector< unsigned int >	Get_MotherIdx ()

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

	TauDecay_GenParticle ()

	~TauDecay_GenParticle ()

Public Member Functions inherited from TauDecay
void	ClassifyDecayMode (unsigned int &MODE_ID, unsigned int &TauBitMask)

bool	hasResonance (unsigned int &TauBitMask, int pdgid)

bool	isTauFinalStateParticle (int pdgid)

bool	isTauParticleCounter (int pdgid)

bool	isTauResonanceCounter (int pdgid)

unsigned int	nPi0 (unsigned int &TauBitMask)

unsigned int	nProng (unsigned int &TauBitMask)

void	Reset ()

	TauDecay ()

	~TauDecay ()

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

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

Private Attributes
unsigned int	MODE_ID

std::vector< unsigned int >	MotherIdx

unsigned int	TauBitMask

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

Additional Inherited Members
Public Types inherited from TauDecay
enum	MODE { MODE_UNKNOWN =0, MODE_ELECTRON, MODE_MUON, MODE_PION, MODE_PIPI0, MODE_3PI, MODE_PI2PI0, MODE_KAON, MODE_K0PI, MODE_KPI0, MODE_3PIPI0, MODE_PI3PI0, MODE_3PI2PI0, MODE_5PI, MODE_5PIPI0, MODE_3PI3PI0, MODE_KPIK, MODE_K0BK0PI, MODE_KK0BPI0, MODE_K2PI0, MODE_KPIPI, MODE_PIK0PI0, MODE_ETAPIPI0, MODE_PIPI0GAM, MODE_KK0B, MODE_PI4PI0, MODE_3PIETA, MODE_PI2PI0ETA, MODE_PI2PI0OMEGA, MODE_3PIOMEGA, MODE_KOMEGA, MODE_K3PI0, MODE_K2PIPI0, MODE_KETA, MODE_K0PI2PI0, MODE_K03PI, MODE_2K0PIPI0, MODE_KPIKPI0, NMODEID }

enum	TauDecayStructure { other =0, OneProng =1, ThreeProng =2, FiveProng =4, OnePi0 =8, TwoPi0 =32, ThreePi0 =64, Res_a1_pm =128, Res_a1_0 =256, Res_rho_pm =512, Res_rho_0 =1024, Res_eta =2048, Res_omega =4096, Res_Kstar_pm =8192, Res_Kstar_0 =16384, KS0_to_pipi =32768 }

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

Detailed Description

Definition at line 23 of file TauDecay_GenParticle.h.

Constructor & Destructor Documentation

TauDecay_GenParticle::TauDecay_GenParticle ( )

Definition at line 7 of file TauDecay_GenParticle.cc.

TauDecay_GenParticle::~TauDecay_GenParticle ( )

Definition at line 13 of file TauDecay_GenParticle.cc.

                                            {
 
 }

Member Function Documentation

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

private

Definition at line 68 of file TauDecay_GenParticle.cc.

References reco::CompositeRefCandidateT< D >::daughter(), i, MotherIdx, reco::CompositeRefCandidateT< D >::numberOfDaughters(), reco::LeafCandidate::status(), and TauDecayProducts.

Referenced by Analyze().

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

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

Definition at line 80 of file TauDecay_GenParticle.cc.

References reco::CompositeRefCandidateT< D >::daughter(), PdtPdgMini::gamma, i, MotherIdx, n, reco::CompositeRefCandidateT< D >::numberOfDaughters(), reco::Candidate::pdgId(), reco::LeafCandidate::pdgId(), and TauDecayProducts.

Referenced by Analyze().

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

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

private

Definition at line 38 of file TauDecay_GenParticle.cc.

References funct::abs(), AddPi0Info(), AddRadInfo(), gather_cfg::cout, reco::CompositeRefCandidateT< D >::daughter(), i, TauDecay::isTauFinalStateParticle(), TauDecay::isTauParticleCounter(), TauDecay::isTauResonanceCounter(), MotherIdx, reco::CompositeRefCandidateT< D >::numberOfDaughters(), reco::LeafCandidate::pdgId(), PdtPdgMini::pi0, reco::LeafCandidate::status(), and TauDecayProducts.

Referenced by AnalyzeTau().

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

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

Definition at line 17 of file TauDecay_GenParticle.cc.

References funct::abs(), Analyze(), TauDecay::ClassifyDecayMode(), reco::CompositeRefCandidateT< D >::daughter(), i, MotherIdx, reco::CompositeRefCandidateT< D >::numberOfDaughters(), reco::LeafCandidate::pdgId(), TauDecay::Reset(), PdtPdgMini::tau_minus, and TauDecayProducts.

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

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