#include <TauDecay_CMSSW.h>

Inheritance diagram for TauDecay_CMSSW:

Public Member Functions
bool	AnalyzeTau (HepMC::GenParticle *Tau, unsigned int &JAK_ID, unsigned int &TauBitMask, bool dores=true, bool dopi0=true)

std::vector< unsigned int >	Get_MotherIdx ()

std::vector< HepMC::GenParticle * >	Get_TauDecayProducts ()

	TauDecay_CMSSW ()

	~TauDecay_CMSSW ()

Public Member Functions inherited from TauDecay
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	AddPi0Info (HepMC::GenParticle *Particle, unsigned int midx)

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

Private Attributes
unsigned int	JAK_ID

std::vector< unsigned int >	MotherIdx

unsigned int	TauBitMask

std::vector< HepMC::GenParticle * >	TauDecayProducts

Additional Inherited Members
Public Types inherited from TauDecay
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 }

Detailed Description

Definition at line 30 of file TauDecay_CMSSW.h.

Constructor & Destructor Documentation

TauDecay_CMSSW::TauDecay_CMSSW ( )

Definition at line 7 of file TauDecay_CMSSW.cc.

                               :
   TauDecay()
 {
 
 }

TauDecay_CMSSW::~TauDecay_CMSSW ( )

Definition at line 13 of file TauDecay_CMSSW.cc.

                                {
 
 } 

Member Function Documentation

void TauDecay_CMSSW::AddPi0Info	(	HepMC::GenParticle *	Particle,
		unsigned int	midx
	)

private

Definition at line 61 of file TauDecay_CMSSW.cc.

References MotherIdx, and TauDecayProducts.

                                                                            {
   if(Particle->status()==1){
     TauDecayProducts.push_back(Particle);
     MotherIdx.push_back(midx);
     return;
   }
   HepMC::GenVertex::particle_iterator des;
   for(des = Particle->end_vertex()->particles_begin(HepMC::children);
       des!= Particle->end_vertex()->particles_end(HepMC::children);++des ) {
     AddPi0Info((*des),midx);
   }
 }

void TauDecay_CMSSW::Analyze	(	HepMC::GenParticle *	Particle,
		unsigned int	midx,
		bool	dores,
		bool	dopi0
	)

private

Definition at line 39 of file TauDecay_CMSSW.cc.

References funct::abs(), gather_cfg::cout, TauDecay::isTauFinalStateParticle(), TauDecay::isTauParticleCounter(), TauDecay::isTauResonanceCounter(), MotherIdx, and TauDecayProducts.

Referenced by AnalyzeTau().

                                                                                                 {
   unsigned int pdgid=abs(Particle->pdg_id());
   isTauResonanceCounter(pdgid);
   if(isTauFinalStateParticle(pdgid)){
     if(!isTauParticleCounter(pdgid)) std::cout << "TauDecay_CMSSW::Analyze WARNING: Unknow Final State Particle in Tau Decay... " << pdgid << std::endl;
     TauDecayProducts.push_back(Particle);
     MotherIdx.push_back(midx);
     return;
   }
   HepMC::GenVertex::particle_iterator des;
   if(Particle->end_vertex()){
     for(des = Particle->end_vertex()->particles_begin(HepMC::children);
     des!= Particle->end_vertex()->particles_end(HepMC::children) && Particle->end_vertex()-> particles_out_size()>0;++des ) {
       Analyze((*des),midx,dores,dopi0);
     }
   }
   else {
     std::cout << "Unstable particle that is undecayed in Tau decay tree. PDG ID: " << pdgid << std::endl;
   }
 }

bool TauDecay_CMSSW::AnalyzeTau	(	HepMC::GenParticle *	Tau,
		unsigned int &	JAK_ID,
		unsigned int &	TauBitMask,
		bool	dores = `true`,
		bool	dopi0 = `true`
	)

Definition at line 17 of file TauDecay_CMSSW.cc.

References funct::abs(), Analyze(), TauDecay::ClassifyDecayMode(), MotherIdx, TauDecay::Reset(), PdtPdgMini::tau_minus, and TauDecayProducts.

                                                                                                                          {
   Reset();
   MotherIdx.clear();
   TauDecayProducts.clear();
   if(abs(Tau->pdg_id())==PdtPdgMini::tau_minus){ // check that it is a tau
     unsigned int Tauidx=TauDecayProducts.size();
     HepMC::GenVertex::particle_iterator des;
     if( Tau->end_vertex()){
       for(des = Tau->end_vertex()->particles_begin(HepMC::children);
       des!= Tau->end_vertex()->particles_end(HepMC::children);++des ) {
     Analyze((*des),Tauidx,dores,dopi0);
       }
       ClassifyDecayMode(JAK_ID,TauBitMask);
       return true;
     }
   }
   return false;
 }