#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 ()

Public Member Functions inherited from reco::GenParticle
GenParticle *	clone () const
	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 ()

virtual	~GenParticle ()
	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 ()
	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...

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

virtual Candidate *	daughter (size_type)
	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...

virtual const Candidate *	mother (size_type=0) const
	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...

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

virtual size_t	numberOfMothers () const
	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...

virtual	~CompositeRefCandidateT ()
	destructor More...

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

virtual 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)

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

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

virtual double	energy () const final
	energy More...

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

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

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

virtual void	fillVertexCovariance (CovarianceMatrix &v) const
	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...

virtual bool	hasMasterClone () const

virtual bool	hasMasterClonePtr () const

virtual bool	isCaloMuon () const

virtual bool	isConvertedPhoton () const

virtual bool	isElectron () const

virtual bool	isGlobalMuon () const

virtual bool	isJet () const

virtual bool	isMuon () const

virtual bool	isPhoton () const

virtual bool	isStandAloneMuon () const

virtual bool	isTrackerMuon () const

	LeafCandidate ()

	LeafCandidate (const Candidate &c)

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

	LeafCandidate (LeafCandidate &rh)

	LeafCandidate (LeafCandidate &&)=default

	LeafCandidate (LeafCandidate const &)=default

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

virtual double	mass () const final
	mass More...

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

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

virtual const CandidateBaseRef &	masterClone () const

virtual const CandidatePtr &	masterClonePtr () const

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

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

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

virtual 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...

virtual size_t	numberOfSourceCandidatePtrs () const

LeafCandidate &	operator= (LeafCandidate &&)=default

LeafCandidate &	operator= (LeafCandidate const &)=default

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

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

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

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

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

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

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

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

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

virtual double	rapidity () const final
	rapidity More...

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

virtual void	setLongLived () final

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

virtual void	setMassConstraint () final

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

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

virtual void	setPdgId (int pdgId) final

virtual void	setPz (double pz) final

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

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

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

virtual CandidatePtr	sourceCandidatePtr (size_type i) const

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

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

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

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

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

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

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

virtual double	vertexNdof () const

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

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

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

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

virtual double	y () const final
	rapidity More...

virtual	~LeafCandidate ()
	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 NULL, and theLeadTrack_.

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

GeneratorTau::GeneratorTau ( )

Definition at line 183 of file GeneratorTau.cc.

184 {

185 }

GeneratorTau::~GeneratorTau ( )

inline

Definition at line 42 of file GeneratorTau.h.

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, NULL, patZpeak::numMuons, convertSQLitetoXML_cfg::output, reco::LeafCandidate::pdgId(), and reco::LeafCandidate::status().

Referenced by init().

                                    {
   //return kUndefined if not a tau
   if (theTau == NULL || 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().

 {
   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 )

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(), and getVisibleFourVectors().

                                                             {
   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().

 {
   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_, NULL, convertSQLitetoXML_cfg::output, and theLeadTrack_.

Referenced by init().

                                                    {
   std::vector<const reco::GenParticle*>::const_iterator thePion =
     genChargedPions_.begin();
   double maxPt = 0;
   const reco::GenParticle* output = NULL;
   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().

                                                  {
   return getOpeningAngle(genChargedPions_);
 }

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

Definition at line 287 of file GeneratorTau.cc.

References convertMCVectorToLorentzVectors(), and genChargedPions_.

                                                              {
   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().

                                                {
   return getOpeningAngle(genGammas_);
 }

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

Definition at line 291 of file GeneratorTau.cc.

References convertMCVectorToLorentzVectors(), and genGammas_.

                                                        {
   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().

                                        {
   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.

                                  {
   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().

                                        {
   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.

                                                              {
   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_.

                                                       {
   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(), and getGammaOpeningAngle().

                                                                   {
   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_.

                                            {
   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().

                                                        {
   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().

                                                                    {
   return convertMCVectorToLorentzVectors(visibleDecayProducts_);
 }

float GeneratorTau::getVisNuAngle ( ) const

Definition at line 3 of file GeneratorTau.cc.

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

                                         {
   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, stableDecayProducts_, theDecayMode_, theLeadTrack_, and visibleDecayProducts_.

Referenced by 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();
   }
 }