Class with analysis utils for simulations. More...

#include <AliMCAnalysisUtils.h>

Inheritance diagram for AliMCAnalysisUtils:

Public Types
enum	mcTypes { kMCPhoton, kMCPrompt, kMCFragmentation, kMCISR, kMCPi0Decay, kMCEtaDecay, kMCOtherDecay, kMCDecayPairLost, kMCDecayPairInCalo, kMCDecayDalitz, kMCConversion, kMCElectron, kMCEFromCFromB, kMCEFromC, kMCEFromB, kMCZDecay, kMCWDecay, kMCMuon, kMCPion, kMCPi0, kMCKaon, kMCEta, kMCProton, kMCAntiProton, kMCNeutron, kMCAntiNeutron, kMCUnknown, kMCBadLabel }

enum	generator { kPythia = 0, kHerwig = 1, kHijing = 2, kBoxLike = 3 }

Public Member Functions
	AliMCAnalysisUtils ()
	Constructor. More...

virtual	~AliMCAnalysisUtils ()
	Destructor. More...

Int_t	CheckCommonAncestor (Int_t index1, Int_t index2, const AliMCEvent *mcevent, Int_t &ancPDG, Int_t &ancStatus, TLorentzVector &momentum, TVector3 &prodVertex)

Int_t	CheckOrigin (Int_t label, AliMCEvent *mcevent, TString selectHeaderName)

Int_t	CheckOrigin (const Int_t labels, Int_t nlabels, AliMCEvent mcevent, TString selectHeaderName, const TObjArray *arrayCluster=0x0)

void	CheckOverlapped2GammaDecay (const Int_t labels, Int_t nlabels, Int_t mesonIndex, const AliMCEvent mcevent, Int_t &tag)

void	CheckLostDecayPair (const TObjArray arrayCluster, Int_t iMom, Int_t iParent, const AliMCEvent mcevent, Int_t &tag)
	Check on AODs if the current decay photon has the second photon companion lost. More...

TLorentzVector	GetMother (Int_t label, const AliMCEvent *mcevent, Bool_t &ok)

TLorentzVector	GetMother (Int_t label, const AliMCEvent *mcevent, Int_t &pdg, Int_t &status, Bool_t &ok)

TLorentzVector	GetMother (Int_t label, const AliMCEvent *mcevent, Int_t &pdg, Int_t &status, Bool_t &ok, Int_t &momLabel)

TLorentzVector	GetGrandMother (Int_t label, const AliMCEvent *mcevent, Int_t &pdg, Int_t &status, Bool_t &ok, Int_t &grandMomLabel, Int_t &greatMomLabel)

TLorentzVector	GetMotherWithPDG (Int_t label, Int_t pdg, const AliMCEvent *mcevent, Bool_t &ok, Int_t &momLabel)

TLorentzVector	GetFirstMotherWithPDG (Int_t label, Int_t pdg, const AliMCEvent *mcevent, Bool_t &ok, Int_t &momLabel, Int_t &gparentlabel)

void	GetMCDecayAsymmetryAngleForPDG (Int_t label, Int_t pdg, const AliMCEvent *mcevent, Float_t &asy, Float_t &angle, Bool_t &ok)
	In case of an eta or pi0 decay into 2 photons, get the asymmetry in the energy of the photons. More...

Int_t	GetNDaughters (Int_t label, const AliMCEvent *mcevent, Bool_t &ok)

TLorentzVector	GetDaughter (Int_t daughter, Int_t label, const AliMCEvent *mcevent, Int_t &pdg, Int_t &status, Bool_t &ok, Int_t &daugLabel, TVector3 &prodVertex)

Int_t	GetNOverlaps (const Int_t label, UInt_t nlabels, Int_t mctag, Int_t mesonLabel, AliMCEvent mcevent, Int_t overpdg, Int_t overlabel)

void	SetTagBit (Int_t &tag, UInt_t set) const

Bool_t	CheckTagBit (Int_t tag, UInt_t test) const

TList *	GetJets (AliMCEvent *mcevent, Bool_t check)

AliGenPythiaEventHeader *	CheckAndGetPythiaEventHeader (AliMCEvent *mcevent, TString selecHeaderName)

TString	GetPythiaHeaderName () const

TString	GetPythiaProcessName () const

Int_t	GetPythiaProcess () const

Int_t	GetPythiaFirstParticle () const

Int_t	GetPythiaVersion () const

void	SetDebug (Int_t deb)

Int_t	GetDebug () const

void	SetMCGenerator (Int_t mcgen)
	Set the generator type. More...

void	SetMCGenerator (TString mcgen)
	Set the generator type. More...

Int_t	GetMCGenerator () const

TString	GetMCGeneratorString () const

void	Print (const Option_t *opt) const
	Print some relevant parameters set for the analysis. More...

void	PrintAncestry (AliMCEvent *mcevent, Int_t label, Int_t nGenerMax=1000) const

void	PrintMCTag (Int_t tag) const
	Print the assigned origins to this particle. More...

Static Public Member Functions
static AliGenPythiaEventHeader *	GetPythiaEventHeader (AliMCEvent *mcevent, TString selecHeaderName, TString &genName, TString &processName, Int_t &process, Int_t &firstParticle, Int_t &pythiaVersion)

Private Member Functions
AliMCAnalysisUtils &	operator= (const AliMCAnalysisUtils &mcu)
	Copy constructor not implemented. More...

	AliMCAnalysisUtils (const AliMCAnalysisUtils &mcu)
	Assignment operator not implemented. More...

Private Attributes
Int_t	fCurrentEvent
	Current Event number - GetJets() More...

Int_t	fDebug
	Debug level. More...

TList *	fJetsList
	List of jets - GetJets() More...

Int_t	fMCGenerator
	MC generator used to generate data in simulation. More...

TString	fMCGeneratorString
	MC generator used to generate data in simulation. More...

TLorentzVector	fDaughMom
	! particle momentum More...

TLorentzVector	fDaughMom2
	! particle momentum More...

TLorentzVector	fMotherMom
	! particle momentum More...

TLorentzVector	fGMotherMom
	! particle momentum More...

AliGenPythiaEventHeader *	fPyGenHead
	! pythia event header of current event More...

TString	fPyGenName
	Pythia header assigned name. More...

TString	fPyProcessName
	Pythia process name, Gamma-Jet or Jet-Jet. More...

Int_t	fPyProcess
	Pythia process code. More...

Int_t	fPyFirstParticle
	First Pythia generated particle in array. More...

Int_t	fPyVersion
	Pythia guessed version. More...

Int_t	fMinPartonicParent
	Minimum label of partonic parent of direct photon. More...

Int_t	fMaxPartonicParent
	Minimum label of partonic parent of direct photon. More...

Detailed Description

Class with analysis utils for simulations.

Class containing utility methods for analysis of simulations, ESD or AOD format. Contains:

method to check the origin of a given track/cluster, in several levels of ancestry
methods to check if a cluster contains several particles, particularly decay photons pairs form same neutral meson.
method to obtain the generated jets

More information can be found in this twiki.

Author: Gustavo Conesa Balbastre Gusta.nosp@m.vo.C.nosp@m.onesa.nosp@m..Bal.nosp@m.bastr.nosp@m.e@ce.nosp@m.rn.ch, LPSC-IN2P3-CNRS

Definition at line 35 of file AliMCAnalysisUtils.h.

Member Enumeration Documentation

enum AliMCAnalysisUtils::generator

Enumerator
kPythia
kHerwig
kHijing
kBoxLike

Definition at line 133 of file AliMCAnalysisUtils.h.

enum AliMCAnalysisUtils::mcTypes

"Mostly" photon parent types on line 1,

then "mostly" electron parent types on line 3, (e.g. electrons can come from pi0 decay) then charged particles on line 4, followed by other and unknown on line 5

Enumerator
kMCPhoton
kMCPrompt
kMCFragmentation
kMCISR
kMCPi0Decay
kMCEtaDecay
kMCOtherDecay
kMCDecayPairLost
kMCDecayPairInCalo
kMCDecayDalitz
kMCConversion
kMCElectron
kMCEFromCFromB
kMCEFromC
kMCEFromB
kMCZDecay
kMCWDecay
kMCMuon
kMCPion
kMCPi0
kMCKaon
kMCEta
kMCProton
kMCAntiProton
kMCNeutron
kMCAntiNeutron
kMCUnknown
kMCBadLabel

Definition at line 53 of file AliMCAnalysisUtils.h.

Constructor & Destructor Documentation

AliMCAnalysisUtils::AliMCAnalysisUtils ( )

Constructor.

Definition at line 40 of file AliMCAnalysisUtils.cxx.

AliMCAnalysisUtils::~AliMCAnalysisUtils ( )

virtual

Destructor.

Definition at line 62 of file AliMCAnalysisUtils.cxx.

AliMCAnalysisUtils::AliMCAnalysisUtils ( const AliMCAnalysisUtils & mcu )

private

Assignment operator not implemented.

Member Function Documentation

AliGenPythiaEventHeader * AliMCAnalysisUtils::CheckAndGetPythiaEventHeader	(	AliMCEvent *	mcevent,
		TString	selectHeaderName
	)

Recover the event header if it is from Pythia. Check it only once per event since it is called many times.

Get and set internally the generated process pythia version and first generated particle in case of cocktail events

Parameters

mcevent	Access to AliVMCEvent
selectHeaderName	String that must have the header name

Definition at line 1513 of file AliMCAnalysisUtils.cxx.

Referenced by CheckOrigin(), CheckTagBit(), and GetNOverlaps().

Int_t AliMCAnalysisUtils::CheckCommonAncestor	(	Int_t	index1,
		Int_t	index2,
		const AliMCEvent *	mcevent,
		Int_t &	ancPDG,
		Int_t &	ancStatus,
		TLorentzVector &	momentum,
		TVector3 &	prodVertex
	)

Check the first common ancestor of 2 clusters, given the most likely labels of the primaries generating such clusters.

Definition at line 75 of file AliMCAnalysisUtils.cxx.

Referenced by AliAnaInsideClusterInvariantMass::CheckLocalMaximaMCOrigin(), AliAnaElectron::FillShowerShapeHistograms(), GetNOverlaps(), AliAnaPhotonConvInCalo::MakeAnalysisFillHistograms(), and AliAnaPi0::MakeAnalysisFillHistograms().

void AliMCAnalysisUtils::CheckLostDecayPair	(	const TObjArray *	arrayCluster,
		Int_t	iMom,
		Int_t	iParent,
		const AliMCEvent *	mcevent,
		Int_t &	tag
	)

Check on AODs if the current decay photon has the second photon companion lost.

Definition at line 789 of file AliMCAnalysisUtils.cxx.

Referenced by CheckOrigin().

Int_t AliMCAnalysisUtils::CheckOrigin	(	Int_t	label,
		AliMCEvent *	mcevent,
		TString	selectHeaderName
	)

Returns: tag with primary particle at the origin of the cluster/track. Here we have only one input MC label not multiple.

Parameters

labels	list of MC labels of cluster
mcevent	pointer to MCEvent()
selectHeaderName	String that must have the header name, for pythia event header/process selection

Definition at line 192 of file AliMCAnalysisUtils.cxx.

Referenced by AliAnaClusterShapeCorrelStudies::ClusterLoopHistograms(), AliAnaCalorimeterQA::ClusterMCHistograms(), AliAnaPhoton::FillAcceptanceHistograms(), AliAnaParticleIsolation::FillAcceptanceHistograms(), AliAnaInsideClusterInvariantMass::GetMCIndex(), AliAnaPhoton::MakeAnalysisFillAOD(), AliAnaElectron::MakeAnalysisFillAOD(), and AliAnaPi0EbE::MakeShowerShapeIdentification().

Int_t AliMCAnalysisUtils::CheckOrigin	(	const Int_t *	labels,
		Int_t	nlabels,
		AliMCEvent *	mcevent,
		TString	selectHeaderName,
		const TObjArray *	arrayCluster = `0x0`
	)

Returns: tag with primary particle(S) at the origin of the cluster/track.

Generally speaking, label is the MC label of a reconstructed entity (track, cluster, etc) for which we want to know something about its heritage, but one can also use it directly with stack particles not connected to reconstructed entities.

Array of clusters needed in case we want to check if the cluster originated from a pi0/eta meson has the companion decay photon in the list of clusters.

Parameters

labels	list of MC labels of cluster
nlabels	total number of labels attached to cluster
mcevent	pointer to MCEvent()
selectHeaderName	String that must have the header name, for pythia event header/process selection
arrayCluster	list of calorimeter clusters, needed to check lost meson decays

Definition at line 216 of file AliMCAnalysisUtils.cxx.

void AliMCAnalysisUtils::CheckOverlapped2GammaDecay	(	const Int_t *	labels,
		Int_t	nlabels,
		Int_t	mesonIndex,
		const AliMCEvent *	mcevent,
		Int_t &	tag
	)

Check if cluster is formed from the contribution of 2 decay photons from pi0 or eta. Input are AOD AliVParticles.

Definition at line 616 of file AliMCAnalysisUtils.cxx.

Referenced by CheckOrigin().

Bool_t AliMCAnalysisUtils::CheckTagBit	(	Int_t	tag,
		UInt_t	test
	)		const

inline

Definition at line 103 of file AliMCAnalysisUtils.h.

Referenced by AliAnaPhoton::ActivityNearCluster(), CheckOrigin(), CheckOverlapped2GammaDecay(), AliAnaClusterShapeCorrelStudies::ClusterLoopHistograms(), AliAnaParticleIsolation::FillAcceptanceHistograms(), AliAnaPi0::FillMCVersusRecDataHistograms(), AliAnaPhoton::FillTrackMatchingResidualHistograms(), AliAnaParticleIsolation::FillTrackMatchingShowerShapeControlHistograms(), AliAnaInsideClusterInvariantMass::GetMCIndex(), GetNOverlaps(), AliAnaElectron::MakeAnalysisFillAOD(), AliAnaParticleHadronCorrelation::MakeAnalysisFillHistograms(), AliAnaPi0::MakeAnalysisFillHistograms(), AliAnaPi0EbE::MakeInvMassInCalorimeter(), AliAnaPi0EbE::MakeInvMassInCalorimeterAndCTS(), PrintMCTag(), and AliAnaClusterShapeCorrelStudies::WeightHistograms().

TLorentzVector AliMCAnalysisUtils::GetDaughter	(	Int_t	idaugh,
		Int_t	label,
		const AliMCEvent *	mcevent,
		Int_t &	pdg,
		Int_t &	status,
		Bool_t &	ok,
		Int_t &	daughlabel,
		TVector3 &	prodVertex
	)

Returns: the kinematics of the particle that generated the signal, its pdg and its status and its label mother.

Definition at line 1071 of file AliMCAnalysisUtils.cxx.

Referenced by AliAnaInsideClusterInvariantMass::CheckLocalMaximaMCOrigin(), AliAnaParticleIsolation::FillAcceptanceHistograms(), AliAnaPhoton::MakeAnalysisFillHistograms(), and AliAnaParticleIsolation::StudyMCConversionRadius().

Int_t AliMCAnalysisUtils::GetDebug ( ) const

inline

Definition at line 131 of file AliMCAnalysisUtils.h.

Referenced by AliAnaCaloTrackCorrBaseClass::InitDebug().

TLorentzVector AliMCAnalysisUtils::GetFirstMotherWithPDG	(	Int_t	label,
		Int_t	pdg,
		const AliMCEvent *	mcevent,
		Bool_t &	ok,
		Int_t &	momlabel,
		Int_t &	gparentlabel
	)

Returns: the kinematics of the particle ancestor for the first time with a given pdg Prompt/fragmentation photons can have multiple times a photon as its parent

Definition at line 1234 of file AliMCAnalysisUtils.cxx.

Referenced by CheckOrigin().

TLorentzVector AliMCAnalysisUtils::GetGrandMother	(	Int_t	label,
		const AliMCEvent *	mcevent,
		Int_t &	pdg,
		Int_t &	status,
		Bool_t &	ok,
		Int_t &	grandMomLabel,
		Int_t &	greatMomLabel
	)

Returns: the kinematics of the particle that generated the signal.

Definition at line 1288 of file AliMCAnalysisUtils.cxx.

Referenced by AliAnaInsideClusterInvariantMass::CheckLocalMaximaMCOrigin(), and GetNOverlaps().

TList * AliMCAnalysisUtils::GetJets	(	AliMCEvent *	mcevent,
		Bool_t	check
	)

Returns: list of jets (TParticles) and index of most likely parton that originated it.

Parameters

mcevent	pointer to AliMCEvent
check	get the list from the pythia header or if already done return just return it.

Definition at line 912 of file AliMCAnalysisUtils.cxx.

Referenced by CheckTagBit().

void AliMCAnalysisUtils::GetMCDecayAsymmetryAngleForPDG	(	Int_t	label,
		Int_t	pdg,
		const AliMCEvent *	mcevent,
		Float_t &	asy,
		Float_t &	angle,
		Bool_t &	ok
	)

In case of an eta or pi0 decay into 2 photons, get the asymmetry in the energy of the photons.

Definition at line 1334 of file AliMCAnalysisUtils.cxx.

Referenced by AliAnaInsideClusterInvariantMass::GetMCPrimaryKine().

Int_t AliMCAnalysisUtils::GetMCGenerator ( ) const

inline

Definition at line 136 of file AliMCAnalysisUtils.h.

TString AliMCAnalysisUtils::GetMCGeneratorString ( ) const

inline

Definition at line 137 of file AliMCAnalysisUtils.h.

TLorentzVector AliMCAnalysisUtils::GetMother	(	Int_t	label,
		const AliMCEvent *	mcevent,
		Bool_t &	ok
	)

Returns: the kinematics of the particle that generated the signal.

Definition at line 1116 of file AliMCAnalysisUtils.cxx.

Referenced by AliAnaInsideClusterInvariantMass::CheckLocalMaximaMCOrigin(), AliAnaPhoton::CocktailGeneratorsClusterOverlaps(), AliAnaInsideClusterInvariantMass::GetMCPrimaryKine(), GetMother(), GetNOverlaps(), AliAnaPi0EbE::MakeAnalysisFillHistograms(), AliAnaPhoton::MakeAnalysisFillHistograms(), AliAnaPi0EbE::MakeShowerShapeIdentification(), and AliAnaParticleIsolation::StudyMCConversionRadius().

TLorentzVector AliMCAnalysisUtils::GetMother	(	Int_t	label,
		const AliMCEvent *	mcevent,
		Int_t &	pdg,
		Int_t &	status,
		Bool_t &	ok
	)

Definition at line 1126 of file AliMCAnalysisUtils.cxx.

TLorentzVector AliMCAnalysisUtils::GetMother	(	Int_t	label,
		const AliMCEvent *	mcevent,
		Int_t &	pdg,
		Int_t &	status,
		Bool_t &	ok,
		Int_t &	momlabel
	)

Returns: the kinematics of the particle that generated the signal, its pdg and its status and its label mother.

Definition at line 1137 of file AliMCAnalysisUtils.cxx.

TLorentzVector AliMCAnalysisUtils::GetMotherWithPDG	(	Int_t	label,
		Int_t	pdg,
		const AliMCEvent *	mcevent,
		Bool_t &	ok,
		Int_t &	momlabel
	)

Returns: the kinematics of the particle that generated the signal.

Definition at line 1171 of file AliMCAnalysisUtils.cxx.

Referenced by AliAnaParticleIsolation::FillAcceptanceHistograms(), AliAnaInsideClusterInvariantMass::GetMCPrimaryKine(), AliAnaPi0EbE::MakeAnalysisFillHistograms(), and AliAnaPi0EbE::MakeShowerShapeIdentification().

Int_t AliMCAnalysisUtils::GetNDaughters	(	Int_t	label,
		const AliMCEvent *	mcevent,
		Bool_t &	ok
	)

Returns: the the number of daughters of a given MC particle.

Definition at line 1394 of file AliMCAnalysisUtils.cxx.

Referenced by AliAnaInsideClusterInvariantMass::CheckLocalMaximaMCOrigin(), and AliAnaParticleIsolation::FillAcceptanceHistograms().

Int_t AliMCAnalysisUtils::GetNOverlaps	(	const Int_t *	label,
		UInt_t	nlabels,
		Int_t	mctag,
		Int_t	mesonLabel,
		AliMCEvent *	mcevent,
		Int_t *	overpdg,
		Int_t *	overlabel
	)

Compare the primary depositing more energy with the rest, if no photon/electron (conversion) or neutral meson as comon ancestor, consider it as other particle contributing. Give as input the meson label in case it was a pi0 or eta merged cluster.

Definition at line 1424 of file AliMCAnalysisUtils.cxx.

Referenced by AliAnaClusterShapeCorrelStudies::ClusterLoopHistograms(), AliAnaPhoton::FillShowerShapeHistograms(), AliAnaParticleIsolation::FillTrackMatchingShowerShapeControlHistograms(), AliAnaCaloTrackCorrBaseClass::GetCocktailGeneratorBackgroundTag(), AliAnaInsideClusterInvariantMass::GetMCPrimaryKine(), and AliAnaPi0EbE::MakeShowerShapeIdentification().

AliGenPythiaEventHeader * AliMCAnalysisUtils::GetPythiaEventHeader	(	AliMCEvent *	mcevent,
		TString	selectHeaderName,
		TString &	genName,
		TString &	processName,
		Int_t &	process,
		Int_t &	firstParticle,
		Int_t &	pythiaVersion
	)

static

Recover the event header if it is from Pythia. Get the generated process pythia version and first generated particle in case of cocktail events

Parameters

mcevent	Access to AliVMCEvent
selectHeaderName	String that must have the header name
genName	String with name of generator
processName	String with type opf event Jet-Jet or Gamma-Jet
process	Integer process number id
firstParticle	Integer of first generated particle in case of cocktail simulations
pythiaVersion	6 or 8

Definition at line 1564 of file AliMCAnalysisUtils.cxx.

Referenced by CheckAndGetPythiaEventHeader(), CheckTagBit(), AliCaloTrackMCReader::FillInputEvent(), and AliCaloTrackReader::FillInputEvent().

Int_t AliMCAnalysisUtils::GetPythiaFirstParticle ( ) const

inline

Definition at line 127 of file AliMCAnalysisUtils.h.

TString AliMCAnalysisUtils::GetPythiaHeaderName ( ) const

inline

Definition at line 124 of file AliMCAnalysisUtils.h.

Int_t AliMCAnalysisUtils::GetPythiaProcess ( ) const

inline

Definition at line 126 of file AliMCAnalysisUtils.h.

TString AliMCAnalysisUtils::GetPythiaProcessName ( ) const

inline

Definition at line 125 of file AliMCAnalysisUtils.h.

Int_t AliMCAnalysisUtils::GetPythiaVersion ( ) const

inline

Definition at line 128 of file AliMCAnalysisUtils.h.

AliMCAnalysisUtils& AliMCAnalysisUtils::operator= ( const AliMCAnalysisUtils & mcu )

private

Copy constructor not implemented.

void AliMCAnalysisUtils::Print ( const Option_t * opt ) const

Print some relevant parameters set for the analysis.

Definition at line 1698 of file AliMCAnalysisUtils.cxx.

Referenced by GetMCGeneratorString().

void AliMCAnalysisUtils::PrintAncestry	(	AliMCEvent *	mcevent,
		Int_t	label,
		Int_t	nGenerMax = `1000`
	)		const

Print info of generated particles, for different generations If no generation specified, all ancestry is printed

Parameters

mcevent	access to AliVMCEvent
label	index of generated particle under investigation
nGenerMax	limit to the number of generations back to the particle under investigation

Definition at line 1719 of file AliMCAnalysisUtils.cxx.

Referenced by GetMCGeneratorString().

void AliMCAnalysisUtils::PrintMCTag ( Int_t tag ) const

Print the assigned origins to this particle.

Definition at line 1766 of file AliMCAnalysisUtils.cxx.

Referenced by GetMCGeneratorString().

void AliMCAnalysisUtils::SetDebug ( Int_t deb )

inline

Definition at line 130 of file AliMCAnalysisUtils.h.

void AliMCAnalysisUtils::SetMCGenerator ( Int_t mcgen )

Set the generator type.

Definition at line 1805 of file AliMCAnalysisUtils.cxx.