#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

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

void	setCollisionId (int s)

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

virtual const_iterator	begin () const
	first daughter const_iterator More...

virtual iterator	begin ()
	first daughter iterator 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_iterator	end () const
	last daughter const_iterator More...

virtual iterator	end ()
	last daughter iterator 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
template<typename S >
daughter_iterator< S >::type	beginFilter (const S &s) const

virtual Vector	boostToCM () const

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

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

template<typename S >
daughter_iterator< S >::type	endFilter (const S &s) const

virtual double	energy () const
	energy More...

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

virtual double	eta () const
	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 ()
	default constructor More...

	LeafCandidate (const Candidate &c)

	LeafCandidate (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...

	LeafCandidate (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...

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

virtual double	mass () const
	mass More...

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

virtual double	massSqr () const
	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
	spatial momentum vector More...

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

virtual double	mtSqr () const
	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

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

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

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

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

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

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

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

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

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

virtual double	rapidity () const
	rapidity More...

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

virtual void	setLongLived ()
	set long lived flag More...

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

virtual void	setMassConstraint ()
	set mass constraint flag More...

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

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

virtual void	setPdgId (int pdgId)

virtual void	setPz (double pz)

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

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

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

virtual CandidatePtr	sourceCandidatePtr (size_type i) const

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

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

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

virtual const Point &	vertex () const
	vertex position 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
	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
	rapidity More...

virtual	~LeafCandidate ()
	destructor More...

Public Member Functions inherited from reco::Candidate
template<typename S >
daughter_iterator< S >::type	beginFilter (const S &s) const

	Candidate ()
	default constructor More...

template<typename S >
daughter_iterator< S >::type	endFilter (const S &s) const

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 Attributes inherited from reco::LeafCandidate
static const unsigned int	longLivedTag = 65536
	long lived flag More...

static const unsigned int	massConstraintTag = 131072
	do mass constraint flag More...

Protected Member Functions inherited from reco::LeafCandidate
void	cacheCartesian () const
	set internal cache More...

void	cachePolar () const
	set internal cache More...

void	clearCache () const
	clear internal cache More...

Protected Attributes inherited from reco::LeafCandidate
bool	cacheCartesianFixed_

bool	cachePolarFixed_
	has cache been set? More...

float	eta_

float	mass_

LorentzVector	p4Cartesian_
	internal cache for p4 More...

PolarLorentzVector	p4Polar_
	internal cache for p4 More...

int	pdgId_
	PDG identifier. More...

float	phi_

float	pt_
	four-momentum Lorentz vector More...

Charge	qx3_
	electric charge More...

int	status_
	status word More...

Point	vertex_
	vertex position More...

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