#include <GeneratorTau.h>

Inheritance diagram for GeneratorTau:

Public Types
enum	tauDecayModeEnum { kElectron, kMuon, kOneProng0pi0, kOneProng1pi0, kOneProng2pi0, kThreeProng0pi0, kThreeProng1pi0, kOther, kUndefined }

Public Types inherited from reco::CompositeRefCandidateT< GenParticleRefVector >
typedef GenParticleRefVector	daughters
	collection of references to daughters More...

typedef GenParticleRefVector	mothers
	collection of references to daughters More...

Public Types inherited from reco::LeafCandidate
typedef int	Charge
	electric charge type More...

typedef CandidateCollection	daughters
	collection of daughter candidates More...

typedef unsigned int	index

typedef math::XYZTLorentzVector	LorentzVector
	Lorentz vector. More...

typedef math::XYZPoint	Point
	point in the space More...

typedef math::PtEtaPhiMLorentzVector	PolarLorentzVector
	Lorentz vector. More...

typedef math::XYZVector	Vector
	point in the space More...

Public Types inherited from reco::Candidate
enum	{ dimension = 3 }

enum	{ size = dimension * (dimension + 1)/2 }
	matix size More...

typedef int	Charge
	electric charge type More...

typedef candidate::const_iterator	const_iterator

typedef math::Error< dimension >::type	CovarianceMatrix
	covariance error matrix (3x3) More...

typedef unsigned int	index
	index type More...

typedef candidate::iterator	iterator

typedef math::XYZTLorentzVector	LorentzVector
	Lorentz vector. More...

typedef math::XYZPoint	Point
	point in the space More...

typedef math::PtEtaPhiMLorentzVector	PolarLorentzVector
	Lorentz vector. More...

typedef size_t	size_type

typedef math::XYZVector	Vector
	point in the space More...

Public Member Functions
tauDecayModeEnum	computeDecayMode (const reco::GenParticle *)

void	computeStableDecayProducts (const reco::GenParticle , std::vector< const reco::GenParticle > &)

LorentzVector	convertHepMCFourVec (const reco::GenParticle *theParticle)

std::vector< LorentzVector >	convertMCVectorToLorentzVectors (const std::vector< const reco::GenParticle * > &theList) const

void	decayToPDGClassification (const reco::GenParticle , std::vector< const reco::GenParticle > &)
	Return list of stable & "semi-stable" tau decay products (e.g. decay the rhos) More...

const reco::GenParticle *	findLeadTrack ()

	GeneratorTau (const reco::GenParticle &input)

	GeneratorTau ()

float	getChargedOpeningAngle () const

std::vector< LorentzVector >	getChargedPions () const

tauDecayModeEnum	getDecayType () const

float	getGammaOpeningAngle () const

std::vector< LorentzVector >	getGammas () const

std::vector< const reco::Candidate * >	getGenChargedPions () const

std::vector< const reco::Candidate * >	getGenGammas () const

std::vector< const reco::Candidate * >	getGenNeutralPions () const

std::vector< const reco::Candidate * >	getGenNu () const

const reco::Candidate *	getLeadTrack () const

float	getOpeningAngle (const std::vector< const reco::GenParticle * > &aCollection) const

std::vector< const reco::Candidate * >	getStableDecayProducts () const

LorentzVector	getVisibleFourVector () const

std::vector< LorentzVector >	getVisibleFourVectors () const

float	getVisNuAngle () const

void	init ()

bool	isFinalStateTau () const

	~GeneratorTau () override

Public Member Functions inherited from reco::GenParticle
GenParticle *	clone () const override
	return a clone More...

int	collisionId () const

bool	fromHardProcessBeforeFSR () const

bool	fromHardProcessDecayed () const

bool	fromHardProcessFinalState () const

	GenParticle ()
	default constructor More...

	GenParticle (const LeafCandidate &c)
	default constructor More...

	GenParticle (Charge q, const LorentzVector &p4, const Point &vtx, int pdgId, int status, bool integerCharge)
	constrocturo from values More...

	GenParticle (Charge q, const PolarLorentzVector &p4, const Point &vtx, int pdgId, int status, bool integerCharge)
	constrocturo from values More...

bool	isDirectHardProcessTauDecayProductFinalState () const

bool	isDirectPromptTauDecayProductFinalState () const

bool	isHardProcess () const

bool	isLastCopy () const

bool	isLastCopyBeforeFSR () const

bool	isMostlyLikePythia6Status3 ()

bool	isPromptDecayed () const

bool	isPromptFinalState () const

void	setCollisionId (int s)

const GenStatusFlags &	statusFlags () const

GenStatusFlags &	statusFlags ()

	~GenParticle () override
	destructor More...

Public Member Functions inherited from reco::CompositeRefCandidateT< GenParticleRefVector >
void	addDaughter (const typename daughters::value_type &)
	add a daughter via a reference More...

void	addMother (const typename mothers::value_type &)
	add a daughter via a reference More...

void	clearDaughters ()
	clear daughter references More...

void	clearMothers ()
	clear mother references More...

CompositeRefCandidateT< GenParticleRefVector > *	clone () const override
	returns a clone of the candidate More...

	CompositeRefCandidateT ()
	default constructor More...

	CompositeRefCandidateT (Charge q, const LorentzVector &p4, const Point &vtx=Point(0, 0, 0), int pdgId=0, int status=0, bool integerCharge=true)
	constructor from values More...

	CompositeRefCandidateT (Charge q, const PolarLorentzVector &p4, const Point &vtx=Point(0, 0, 0), int pdgId=0, int status=0, bool integerCharge=true)
	constructor from values More...

	CompositeRefCandidateT (const LeafCandidate &c)
	constructor from a particle More...

const Candidate *	daughter (size_type) const override
	return daughter at a given position, i = 0, ... numberOfDaughters() - 1 (read only mode) More...

Candidate *	daughter (size_type) override
	return daughter at a given position, i = 0, ... numberOfDaughters() - 1 More...

daughters::value_type	daughterRef (size_type i) const
	reference to daughter at given position More...

const daughters &	daughterRefVector () const
	references to daughtes More...

const Candidate *	mother (size_type=0) const override
	return mother at a given position, i = 0, ... numberOfMothers() - 1 (read only mode) More...

daughters::value_type	motherRef (size_type i=0) const
	reference to mother at given position More...

const mothers &	motherRefVector () const
	references to mothers More...

size_t	numberOfDaughters () const override
	number of daughters More...

size_t	numberOfMothers () const override
	number of mothers More...

void	resetDaughters (const edm::ProductID &id)
	set daughters product ID More...

void	resetMothers (const edm::ProductID &id)
	set mother product ID More...

	~CompositeRefCandidateT () override
	destructor More...

Public Member Functions inherited from reco::LeafCandidate
Vector	boostToCM () const final

int	charge () const final
	electric charge More...

void	construct (int qx3, float pt, float eta, float phi, float mass, const Point &vtx, int pdgId, int status)

Candidate *	daughter (const std::string &s) override
	return daughter with a specified role name More...

const Candidate *	daughter (const std::string &s) const override
	return daughter with a specified role name More...

double	energy () const final
	energy More...

double	et () const final
	transverse energy More...

double	et2 () const final
	transverse energy squared (use this for cut!) More...

double	eta () const final
	momentum pseudorapidity More...

void	fillVertexCovariance (CovarianceMatrix &v) const override
	fill SMatrix More...

template<typename T >
T	get () const
	get a component More...

template<typename T , typename Tag >
T	get () const
	get a component More...

template<typename T >
T	get (size_type i) const
	get a component More...

template<typename T , typename Tag >
T	get (size_type i) const
	get a component More...

bool	hasMasterClone () const override

bool	hasMasterClonePtr () const override

bool	isCaloMuon () const override

bool	isConvertedPhoton () const override

bool	isElectron () const override

bool	isGlobalMuon () const override

bool	isJet () const override

bool	isMuon () const override

bool	isPhoton () const override

bool	isStandAloneMuon () const override

bool	isTrackerMuon () const override

	LeafCandidate ()

	LeafCandidate (const Candidate &c)

template<typename... Args>
	LeafCandidate (Args &&...args)

	LeafCandidate (LeafCandidate &rh)

	LeafCandidate (LeafCandidate &&)=default

	LeafCandidate (LeafCandidate const &)=default

bool	longLived () const final
	is long lived? More...

double	mass () const final
	mass More...

bool	massConstraint () const final
	do mass constraint? More...

double	massSqr () const final
	mass squared More...

const CandidateBaseRef &	masterClone () const override

const CandidatePtr &	masterClonePtr () const override

template<typename Ref >
Ref	masterRef () const
	cast master clone reference to a concrete type More...

Vector	momentum () const final
	spatial momentum vector More...

double	mt () const final
	transverse mass More...

double	mtSqr () const final
	transverse mass squared More...

template<typename T >
size_type	numberOf () const
	number of components More...

template<typename T , typename Tag >
size_type	numberOf () const
	number of components More...

size_t	numberOfSourceCandidatePtrs () const override

LeafCandidate &	operator= (LeafCandidate &&)=default

LeafCandidate &	operator= (LeafCandidate const &)=default

double	p () const final
	magnitude of momentum vector More...

const LorentzVector &	p4 () const final
	four-momentum Lorentz vector More...

int	pdgId () const final
	PDG identifier. More...

double	phi () const final
	momentum azimuthal angle More...

const PolarLorentzVector &	polarP4 () const final
	four-momentum Lorentz vector More...

double	pt () const final
	transverse momentum More...

double	px () const final
	x coordinate of momentum vector More...

double	py () const final
	y coordinate of momentum vector More...

double	pz () const final
	z coordinate of momentum vector More...

double	rapidity () const final
	rapidity More...

void	setCharge (Charge q) final
	set electric charge More...

void	setLongLived () final

void	setMass (double m) final
	set particle mass More...

void	setMassConstraint () final

void	setP4 (const LorentzVector &p4) final
	set 4-momentum More...

void	setP4 (const PolarLorentzVector &p4) final
	set 4-momentum More...

void	setPdgId (int pdgId) final

void	setPz (double pz) final

void	setStatus (int status) final
	set status word More...

void	setThreeCharge (Charge qx3) final
	set electric charge More...

void	setVertex (const Point &vertex) override
	set vertex More...

CandidatePtr	sourceCandidatePtr (size_type i) const override

int	status () const final
	status word More...

double	theta () const final
	momentum polar angle More...

int	threeCharge () const final
	electric charge More...

const Point &	vertex () const override
	vertex position (overwritten by PF...) More...

double	vertexChi2 () const override
	chi-squares More...

double	vertexCovariance (int i, int j) const override
	(i, j)-th element of error matrix, i, j = 0, ... 2 More...

CovarianceMatrix	vertexCovariance () const final
	return SMatrix More...

double	vertexNdof () const override

double	vertexNormalizedChi2 () const override
	chi-squared divided by n.d.o.f. More...

double	vx () const override
	x coordinate of vertex position More...

double	vy () const override
	y coordinate of vertex position More...

double	vz () const override
	z coordinate of vertex position More...

double	y () const final
	rapidity More...

	~LeafCandidate () override
	destructor More...

Public Member Functions inherited from reco::Candidate
const_iterator	begin () const
	first daughter const_iterator More...

iterator	begin ()
	first daughter iterator More...

virtual const Track *	bestTrack () const

	Candidate ()
	default constructor More...

virtual float	dxyError () const
	uncertainty on dxy More...

virtual float	dzError () const
	uncertainty on dz More...

const_iterator	end () const
	last daughter const_iterator More...

iterator	end ()
	last daughter iterator More...

template<typename T >
T	get () const
	get a component More...

template<typename T , typename Tag >
T	get () const
	get a component More...

template<typename T >
T	get (size_type i) const
	get a component More...

template<typename T , typename Tag >
T	get (size_type i) const
	get a component More...

template<typename Ref >
Ref	masterRef () const
	cast master clone reference to a concrete type More...

template<typename T >
size_type	numberOf () const
	number of components More...

template<typename T , typename Tag >
size_type	numberOf () const
	number of components More...

virtual void	setSourceCandidatePtr (const CandidatePtr &ptr)
	Set the ptr to the source Candidate. More...

virtual	~Candidate ()
	destructor More...

Private Attributes
int	aFinalStateTau_

Angle< LorentzVector >	angleFinder

DeltaR< LorentzVector >	deltaRFinder

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

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

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

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

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

tauDecayModeEnum	theDecayMode_

const reco::GenParticle *	theLeadTrack_

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

Additional Inherited Members
Static Public Member Functions inherited from reco::CompositeRefCandidateT< GenParticleRefVector >
static short	Class_Version ()

Detailed Description

Definition at line 28 of file GeneratorTau.h.

Member Enumeration Documentation

enum GeneratorTau::tauDecayModeEnum

Enumerator
kElectron
kMuon
kOneProng0pi0
kOneProng1pi0
kOneProng2pi0
kThreeProng0pi0
kThreeProng1pi0
kOther
kUndefined

Definition at line 31 of file GeneratorTau.h.

                             {kElectron, kMuon,
          kOneProng0pi0, kOneProng1pi0, kOneProng2pi0,
          kThreeProng0pi0, kThreeProng1pi0,
          kOther, kUndefined};

Constructor & Destructor Documentation

GeneratorTau::GeneratorTau ( const reco::GenParticle & input )

inline

Definition at line 37 of file GeneratorTau.h.

References GeneratorTau(), and theLeadTrack_.

                                                 :GenParticle(input){
         theLeadTrack_ = nullptr;
       };

GeneratorTau::GeneratorTau ( )

Definition at line 183 of file GeneratorTau.cc.

Referenced by GeneratorTau().

184 {

185 }

GeneratorTau::~GeneratorTau ( )

inlineoverride

Definition at line 42 of file GeneratorTau.h.

References getGenChargedPions(), getGenGammas(), getGenNeutralPions(), getGenNu(), and getStableDecayProducts().

42 {};

Member Function Documentation

GeneratorTau::tauDecayModeEnum GeneratorTau::computeDecayMode ( const reco::GenParticle * theTau )

find the decay products, in terms of the PDG table (eg 1 pi0, 1 pi+, etc)

Definition at line 48 of file GeneratorTau.cc.

References funct::abs(), decayToPDGClassification(), kElectron, kMuon, kOneProng0pi0, kOneProng1pi0, kOneProng2pi0, kOther, kThreeProng0pi0, kThreeProng1pi0, kUndefined, convertSQLitetoXML_cfg::output, HiggsValidation_cfi::pdg_id, reco::LeafCandidate::pdgId(), and reco::LeafCandidate::status().

Referenced by init(), and isFinalStateTau().

                                    {
   //return kUndefined if not a tau
   if (theTau == nullptr || std::abs(theTau->pdgId()) != 15
       || theTau->status() != 2)
     return kUndefined;
 
   tauDecayModeEnum output;
 
   //counters to determine decay type (adapted from Ricardo's code)
   int numElectrons      = 0;
   int numMuons          = 0;
   int numChargedPions   = 0;
   int numNeutralPions   = 0;
   int numNeutrinos      = 0;
   int numOtherParticles = 0;
 
   std::vector<const reco::GenParticle* > pdgDecayProductTypes;
 
   GeneratorTau::decayToPDGClassification(theTau, pdgDecayProductTypes);
 
   for (std::vector<const reco::GenParticle* >::const_iterator decayProduct =
       pdgDecayProductTypes.begin();
       decayProduct != pdgDecayProductTypes.end(); ++decayProduct) {
     int pdg_id = std::abs( (*decayProduct)->pdgId() );
     //edm::LogInfo("GeneratorTau") << "Has decay product w/ PDG ID: " << pdg_id;
     if (pdg_id == 11) numElectrons++;
     else if (pdg_id == 13) numMuons++;
     else if (pdg_id == 211) numChargedPions++;
     else if (pdg_id == 111) numNeutralPions++;
     else if (pdg_id == 12 ||
         pdg_id == 14 ||
         pdg_id == 16)  numNeutrinos++;
     else if (pdg_id != 22)
       numOtherParticles++;
   }
   output = kOther;
 
   //determine tauDecayMode
   if ( numOtherParticles == 0 ){
     if ( numElectrons == 1 ){
       //--- tau decays into electrons
       output = kElectron;
     } else if ( numMuons == 1 ){
       //--- tau decays into muons
       output = kMuon;
     } else {
       //--- hadronic tau decays
       switch ( numChargedPions ){
         case 1 :
           switch ( numNeutralPions ){
             case 0 :
               output = kOneProng0pi0;
               break;
             case 1 :
               output = kOneProng1pi0;
               break;
             case 2 :
               output = kOneProng2pi0;
               break;
           }
           break;
         case 3 :
           switch ( numNeutralPions ){
             case 0 :
               output = kThreeProng0pi0;
               break;
             case 1 :
               output = kThreeProng1pi0;
               break;
           }
           break;
       }
     }
   }
   return output;
 }

void GeneratorTau::computeStableDecayProducts	(	const reco::GenParticle *	theParticle,
		std::vector< const reco::GenParticle * > &	container
	)

Definition at line 161 of file GeneratorTau.cc.

References reco::Candidate::Candidate(), reco::CompositeRefCandidateT< D >::daughter(), reco::CompositeRefCandidateT< GenParticleRefVector >::daughter(), reco::CompositeRefCandidateT< D >::numberOfDaughters(), and reco::LeafCandidate::status().

Referenced by init(), and isFinalStateTau().

 {
   if (theParticle)
   {
     if (theParticle->status() == 1) //status = 1 indicates final state particle
     {
       //edm::LogInfo("GeneratorTau") << "computeStableDecayProducts: Found a final state daughter with status: " << theParticle->status() << " Num stable decay products so far: " << container.size();
       container.push_back(theParticle);
     }
     else
     {
       unsigned int nDaughters = theParticle->numberOfDaughters();
       for (size_t dIter = 0; dIter < nDaughters; ++dIter)
       {
         const Candidate * daughter = theParticle->daughter(dIter);
         //edm::LogInfo("Debug") << "Recursing on daughter with PDG: " << daughter->pdgId();
         GeneratorTau::computeStableDecayProducts(static_cast<const reco::GenParticle*>(daughter), container);
       }
     }
   }
 }

LorentzVector GeneratorTau::convertHepMCFourVec ( const reco::GenParticle * theParticle )

Referenced by isFinalStateTau().

std::vector< LorentzVector > GeneratorTau::convertMCVectorToLorentzVectors ( const std::vector< const reco::GenParticle * > & theList ) const

Definition at line 231 of file GeneratorTau.cc.

References convertSQLitetoXML_cfg::output.

Referenced by getChargedPions(), getGammas(), getVisibleFourVectors(), and isFinalStateTau().

                                                             {
   std::vector<LorentzVector> output;
   std::vector<const reco::GenParticle*>::const_iterator theParticle;
   for (theParticle = theList.begin();
       theParticle != theList.end(); ++theParticle) {
     output.push_back( (*theParticle)->p4() );
   }
   return output;
 }

void GeneratorTau::decayToPDGClassification	(	const reco::GenParticle *	theParticle,
		std::vector< const reco::GenParticle * > &	container
	)

Return list of stable & "semi-stable" tau decay products (e.g. decay the rhos)

Definition at line 129 of file GeneratorTau.cc.

References funct::abs(), reco::Candidate::Candidate(), reco::CompositeRefCandidateT< D >::daughter(), reco::CompositeRefCandidateT< GenParticleRefVector >::daughter(), genNeutralPions_, reco::CompositeRefCandidateT< D >::numberOfDaughters(), reco::LeafCandidate::pdgId(), and reco::LeafCandidate::status().

Referenced by computeDecayMode(), and isFinalStateTau().

 {
   if (theParticle)
   {
     //edm::LogInfo("Debug") << "It's non-null";
     int pdgId = std::abs(theParticle->pdgId());
     //edm::LogInfo("Debug") << "PDGID = " << pdgId << " Status = " << theStatus;
     if (theParticle->status() == 1 || pdgId == 211 || pdgId == 111 || pdgId == 11 || pdgId == 13)
     {
       //edm::LogInfo("Debug") << "Adding to container...";
       container.push_back(theParticle);
       //add neutral pions and this step....
       if (pdgId == 111)
         genNeutralPions_.push_back(theParticle);
     }
     else
     {
       unsigned int nDaughters = theParticle->numberOfDaughters();
       for (size_t dIter = 0; dIter < nDaughters; ++dIter)
       {
         const Candidate * daughter = theParticle->daughter(dIter);
         //edm::LogInfo("Debug") << "Recursing on daughter with PDG: " << daughter->pdgId();
         GeneratorTau::decayToPDGClassification(static_cast<const reco::GenParticle*>(daughter), container);
       }
 
     }
   }
 }

const reco::GenParticle * GeneratorTau::findLeadTrack ( )

Definition at line 13 of file GeneratorTau.cc.

References genChargedPions_, MuonErrorMatrixAnalyzer_cfi::maxPt, convertSQLitetoXML_cfg::output, and theLeadTrack_.

Referenced by init(), and isFinalStateTau().

                                                    {
   std::vector<const reco::GenParticle*>::const_iterator thePion =
     genChargedPions_.begin();
   double maxPt = 0;
   const reco::GenParticle* output = nullptr;
   for (; thePion != genChargedPions_.end(); ++thePion) {
     if ((*thePion)->pt() > maxPt) {
       maxPt = (*thePion)->pt();
       output = (*thePion);
     }
   }
   theLeadTrack_ = output;
   return output;
 }

float GeneratorTau::getChargedOpeningAngle ( ) const

Definition at line 40 of file GeneratorTau.cc.

References genChargedPions_, and getOpeningAngle().

Referenced by isFinalStateTau().

                                                  {
   return getOpeningAngle(genChargedPions_);
 }

std::vector< LorentzVector > GeneratorTau::getChargedPions ( ) const

Definition at line 287 of file GeneratorTau.cc.

References convertMCVectorToLorentzVectors(), and genChargedPions_.

Referenced by isFinalStateTau().

                                                              {
   return convertMCVectorToLorentzVectors(genChargedPions_);
 }

tauDecayModeEnum GeneratorTau::getDecayType ( ) const

inline

Definition at line 50 of file GeneratorTau.h.

References theDecayMode_.

50 {return theDecayMode_;};

GeneratorTau::theDecayMode_

tauDecayModeEnum theDecayMode_

Definition: GeneratorTau.h:88

float GeneratorTau::getGammaOpeningAngle ( ) const

Definition at line 44 of file GeneratorTau.cc.

References genGammas_, and getOpeningAngle().

Referenced by isFinalStateTau().

                                                {
   return getOpeningAngle(genGammas_);
 }

std::vector< LorentzVector > GeneratorTau::getGammas ( ) const

Definition at line 291 of file GeneratorTau.cc.

References convertMCVectorToLorentzVectors(), and genGammas_.

Referenced by isFinalStateTau().

                                                        {
   return convertMCVectorToLorentzVectors(genGammas_);
 }

std::vector< const reco::Candidate * > GeneratorTau::getGenChargedPions ( ) const

Definition at line 243 of file GeneratorTau.cc.

References genChargedPions_, and convertSQLitetoXML_cfg::output.

Referenced by TruthTauDecayModeProducer::produce(), and ~GeneratorTau().

                                        {
   std::vector<const reco::Candidate*> output;
   std::vector<const GenParticle*>::const_iterator iter;
   for (iter = genChargedPions_.begin(); iter != genChargedPions_.end(); ++iter)
     output.push_back(static_cast<const reco::Candidate*>(*iter));
   return output;
 }

std::vector< const reco::Candidate * > GeneratorTau::getGenGammas ( ) const

Definition at line 261 of file GeneratorTau.cc.

References genGammas_, and convertSQLitetoXML_cfg::output.

Referenced by ~GeneratorTau().

                                  {
   std::vector<const reco::Candidate*> output;
   std::vector<const GenParticle*>::const_iterator iter;
   for (iter = genGammas_.begin(); iter != genGammas_.end(); ++iter)
     output.push_back(static_cast<const reco::Candidate*>(*iter));
   return output;
 }

std::vector< const reco::Candidate * > GeneratorTau::getGenNeutralPions ( ) const

Definition at line 252 of file GeneratorTau.cc.

References genNeutralPions_, and convertSQLitetoXML_cfg::output.

Referenced by TruthTauDecayModeProducer::produce(), and ~GeneratorTau().

                                        {
   std::vector<const reco::Candidate*> output;
   std::vector<const GenParticle*>::const_iterator iter;
   for (iter = genNeutralPions_.begin(); iter != genNeutralPions_.end(); ++iter)
     output.push_back(static_cast<const reco::Candidate*>(*iter));
   return output;
 }

std::vector< const reco::Candidate * > GeneratorTau::getGenNu ( ) const

Definition at line 279 of file GeneratorTau.cc.

References genNus_, and convertSQLitetoXML_cfg::output.

Referenced by ~GeneratorTau().

                                                              {
   std::vector<const reco::Candidate*> output;
   std::vector<const GenParticle*>::const_iterator iter;
   for (iter = genNus_.begin(); iter != genNus_.end(); ++iter)
     output.push_back(static_cast<const reco::Candidate*>(*iter));
   return output;
 }

const reco::Candidate * GeneratorTau::getLeadTrack ( ) const

Definition at line 9 of file GeneratorTau.cc.

References reco::Candidate::Candidate(), and theLeadTrack_.

Referenced by isFinalStateTau().

                                                       {
   return static_cast<const Candidate*>(theLeadTrack_);
 }

float GeneratorTau::getOpeningAngle ( const std::vector< const reco::GenParticle * > & aCollection ) const

Definition at line 28 of file GeneratorTau.cc.

References angleFinder, convertSQLitetoXML_cfg::output, reco::LeafCandidate::p4(), and theLeadTrack_.

Referenced by getChargedOpeningAngle(), getGammaOpeningAngle(), and isFinalStateTau().

                                                                   {
   double output = 0;
   std::vector<const reco::GenParticle*>::const_iterator theObject =
     theCollection.begin();
   for (; theObject != theCollection.end(); ++theObject) {
     if (output < angleFinder(theLeadTrack_->p4(), (*theObject)->p4()))
       output = angleFinder(theLeadTrack_->p4(), (*theObject)->p4());
   }
   return output;
 }

std::vector< const reco::Candidate * > GeneratorTau::getStableDecayProducts ( ) const

Definition at line 270 of file GeneratorTau.cc.

References convertSQLitetoXML_cfg::output, and stableDecayProducts_.

Referenced by ~GeneratorTau().

                                            {
   std::vector<const reco::Candidate*> output;
   std::vector<const GenParticle*>::const_iterator iter;
   for (iter = stableDecayProducts_.begin();
       iter != stableDecayProducts_.end(); ++iter)
     output.push_back(static_cast<const reco::Candidate*>(*iter));
   return output;
 }

LorentzVector GeneratorTau::getVisibleFourVector ( ) const

Definition at line 299 of file GeneratorTau.cc.

References getVisibleFourVectors(), and convertSQLitetoXML_cfg::output.

Referenced by getVisNuAngle(), and isFinalStateTau().

                                                        {
   LorentzVector output;
   std::vector<LorentzVector> tempForSum = getVisibleFourVectors();
   for (std::vector<LorentzVector>::iterator iter = tempForSum.begin();
       iter != tempForSum.end(); ++iter)
     output += (*iter);
   return output;
 }

std::vector< LorentzVector > GeneratorTau::getVisibleFourVectors ( ) const

Definition at line 295 of file GeneratorTau.cc.

References convertMCVectorToLorentzVectors(), and visibleDecayProducts_.

Referenced by getVisibleFourVector(), and isFinalStateTau().

                                                                    {
   return convertMCVectorToLorentzVectors(visibleDecayProducts_);
 }

float GeneratorTau::getVisNuAngle ( ) const

Definition at line 3 of file GeneratorTau.cc.

References angleFinder, getVisibleFourVector(), and reco::LeafCandidate::p4().

Referenced by isFinalStateTau().

                                         {
   LorentzVector suckVector = getVisibleFourVector();
   LorentzVector suckNuVector = this->p4() - suckVector;
   return angleFinder(suckVector, suckNuVector);
 }

void GeneratorTau::init ( void )

Definition at line 189 of file GeneratorTau.cc.

References funct::abs(), aFinalStateTau_, computeDecayMode(), computeStableDecayProducts(), findLeadTrack(), genChargedPions_, genGammas_, genNus_, kUndefined, HiggsValidation_cfi::pdg_id, stableDecayProducts_, theDecayMode_, theLeadTrack_, and visibleDecayProducts_.

Referenced by isFinalStateTau(), and TruthTauDecayModeProducer::produce().

 {
   //make sure this tau really decays
   theDecayMode_        = kUndefined;
   aFinalStateTau_      = false;
 
   //get Decaymode
   //edm::LogInfo("GeneratorTau") << "Computing decay mode..";
   theDecayMode_ = computeDecayMode(this);
 
   //make sure it is a real tau decay
   if (theDecayMode_ != kUndefined) {
     aFinalStateTau_ = true;
     //edm::LogInfo("GeneratorTau") << "Found decay type: " << theDecayMode_ << ", computing stable decay products.";
     //get the stable decay products
     computeStableDecayProducts(this, stableDecayProducts_);
     //from the stable products, fill the lists
     //edm::LogInfo("GeneratorTau") << "Found " << stableDecayProducts_.size() << " stable decay products, filtering.";
     for (std::vector<const reco::GenParticle*>::const_iterator iter = stableDecayProducts_.begin();
         iter != stableDecayProducts_.end();
         ++iter)
     {
       //fill vectors
       int pdg_id = std::abs( (*iter)->pdgId() );
       if (pdg_id == 16 || pdg_id == 12 || pdg_id == 14)
         genNus_.push_back( (*iter) );
       else  {
         visibleDecayProducts_.push_back( (*iter) );
         if (pdg_id == 211 || (*iter)->charge() != 0)
           genChargedPions_.push_back( (*iter) );
         else if (pdg_id == 22)
           genGammas_.push_back( (*iter) );
       }
     }
     // find the lead charged object
     theLeadTrack_ = findLeadTrack();
   }
 }