Quick Links: ROOT Homepage Class Index Class Hierarchy

Source: header file source file inheritance tree (.pdf)

Sections: class description function members data members class charts

class AliEMCALUnfolding: public TObject

  Base class for the cluster unfolding algorithm
*-- Author: Adam Matyja (SUBATECH)
  Based on unfolding in clusterizerv1 done by Cynthia Hadjidakis
-- Unfolding for eta~0: Cynthia Hadjidakis - still in AliEMCALCLusterizerv1
-- Unfolding extension for whole EMCAL: Adam Matyja (SUBATECH & INP PAN)

  unfolds the clusters having several local maxima.

Function Members (Methods)

public:

	AliEMCALUnfolding()
	AliEMCALUnfolding(AliEMCALGeometry* geometry)
	AliEMCALUnfolding(AliEMCALGeometry* geometry, Float_t ECALocMaxCut, Double_t* SSPars, Double_t* Par5, Double_t* Par6)
virtual	~AliEMCALUnfolding()
void	TObject::AbstractMethod(const char* method) const
virtual void	TObject::AppendPad(Option_t* option = "")
virtual void	TObject::Browse(TBrowser* b)
static TClass*	Class()
virtual const char*	TObject::ClassName() const
virtual void	TObject::Clear(Option_t* = "")
virtual TObject*	TObject::Clone(const char* newname = "") const
virtual Int_t	TObject::Compare(const TObject* obj) const
virtual void	TObject::Copy(TObject& object) const
virtual void	TObject::Delete(Option_t* option = "")MENU
virtual Int_t	TObject::DistancetoPrimitive(Int_t px, Int_t py)
virtual void	TObject::Draw(Option_t* option = "")
virtual void	TObject::DrawClass() constMENU
virtual TObject*	TObject::DrawClone(Option_t* option = "") constMENU
virtual void	TObject::Dump() constMENU
virtual void	TObject::Error(const char* method, const char* msgfmt) const
virtual void	TObject::Execute(const char* method, const char* params, Int_t* error = 0)
virtual void	TObject::Execute(TMethod* method, TObjArray* params, Int_t* error = 0)
virtual void	TObject::ExecuteEvent(Int_t event, Int_t px, Int_t py)
virtual void	TObject::Fatal(const char* method, const char* msgfmt) const
virtual TObject*	TObject::FindObject(const char* name) const
virtual TObject*	TObject::FindObject(const TObject* obj) const
virtual TClonesArray*	GetDigitsArr() const
virtual Option_t*	TObject::GetDrawOption() const
static Long_t	TObject::GetDtorOnly()
virtual Float_t	GetECALocalMaxCut() const
virtual const char*	TObject::GetIconName() const
virtual const char*	TObject::GetName() const
virtual Int_t	GetNumberOfECAClusters() const
virtual char*	TObject::GetObjectInfo(Int_t px, Int_t py) const
static Bool_t	TObject::GetObjectStat()
virtual Option_t*	TObject::GetOption() const
virtual Double_t*	GetPar5() const
virtual Double_t*	GetPar6() const
virtual TObjArray*	GetRecPoints() const
virtual Bool_t	GetRejectBelowThreshold() const
virtual Double_t*	GetShowerShapeParams() const
virtual Float_t	GetThreshold() const
virtual const char*	TObject::GetTitle() const
virtual UInt_t	TObject::GetUniqueID() const
virtual Bool_t	TObject::HandleTimer(TTimer* timer)
virtual ULong_t	TObject::Hash() const
virtual void	TObject::Info(const char* method, const char* msgfmt) const
virtual Bool_t	TObject::InheritsFrom(const char* classname) const
virtual Bool_t	TObject::InheritsFrom(const TClass* cl) const
virtual void	Init()
virtual void	TObject::Inspect() constMENU
void	TObject::InvertBit(UInt_t f)
virtual TClass*	IsA() const
virtual Bool_t	TObject::IsEqual(const TObject* obj) const
virtual Bool_t	TObject::IsFolder() const
Bool_t	TObject::IsOnHeap() const
virtual Bool_t	TObject::IsSortable() const
Bool_t	TObject::IsZombie() const
virtual void	TObject::ls(Option_t* option = "") const
virtual void	MakeUnfolding()
void	TObject::MayNotUse(const char* method) const
virtual Bool_t	TObject::Notify()
void	TObject::Obsolete(const char* method, const char* asOfVers, const char* removedFromVers) const
static void	TObject::operator delete(void* ptr)
static void	TObject::operator delete(void* ptr, void* vp)
static void	TObject::operator delete[](void* ptr)
static void	TObject::operator delete[](void* ptr, void* vp)
void*	TObject::operator new(size_t sz)
void*	TObject::operator new(size_t sz, void* vp)
void*	TObject::operator new[](size_t sz)
void*	TObject::operator new[](size_t sz, void* vp)
virtual void	TObject::Paint(Option_t* option = "")
virtual void	TObject::Pop()
virtual void	TObject::Print(Option_t* option = "") const
virtual Int_t	TObject::Read(const char* name)
virtual void	TObject::RecursiveRemove(TObject* obj)
void	TObject::ResetBit(UInt_t f)
virtual void	TObject::SaveAs(const char* filename = "", Option_t* option = "") constMENU
virtual void	TObject::SavePrimitive(basic_ostream<char,char_traits<char> >& out, Option_t* option = "")
void	TObject::SetBit(UInt_t f)
void	TObject::SetBit(UInt_t f, Bool_t set)
virtual void	SetDigitsArr(TClonesArray* digit)
virtual void	TObject::SetDrawOption(Option_t* option = "")MENU
static void	TObject::SetDtorOnly(void* obj)
virtual void	SetECALocalMaxCut(Float_t cut)
virtual void	SetInput(Int_t numberOfECAClusters, TObjArray* recPoints, TClonesArray* digitsArr)
virtual void	SetNumberOfECAClusters(Int_t n)
static void	TObject::SetObjectStat(Bool_t stat)
virtual void	SetPar5(Double_t* pars)
virtual void	SetPar6(Double_t* pars)
virtual void	SetRecPoints(TObjArray* rec)
virtual void	SetRejectBelowThreshold(Bool_t reject)
virtual void	SetShowerShapeParams(Double_t* pars)
virtual void	SetThreshold(Float_t energy)
virtual void	TObject::SetUniqueID(UInt_t uid)
static Double_t	ShowerShapeV2(Double_t x, Double_t y)
virtual void	ShowMembers(TMemberInspector&)
virtual void	Streamer(TBuffer&)
void	StreamerNVirtual(TBuffer& ClassDef_StreamerNVirtual_b)
virtual void	TObject::SysError(const char* method, const char* msgfmt) const
Bool_t	TObject::TestBit(UInt_t f) const
Int_t	TObject::TestBits(UInt_t f) const
static void	UnfoldingChiSquareV2(Int_t& nPar, Double_t* Grad, Double_t& fret, Double_t* x, Int_t iflag)
virtual Int_t	UnfoldOneCluster(AliEMCALRecPoint* iniTower, Int_t nMax, AliEMCALDigit** maxAt, Float_t* maxAtEnergy, TObjArray* list)
virtual void	TObject::UseCurrentStyle()
virtual void	TObject::Warning(const char* method, const char* msgfmt) const
virtual Int_t	TObject::Write(const char* name = 0, Int_t option = 0, Int_t bufsize = 0)
virtual Int_t	TObject::Write(const char* name = 0, Int_t option = 0, Int_t bufsize = 0) const

protected:

virtual void	TObject::DoError(int level, const char* location, const char* fmt, va_list va) const
void	TObject::MakeZombie()

private:

	AliEMCALUnfolding(const AliEMCALUnfolding&)
static void	EvalPar5(Double_t phi)
static void	EvalPar6(Double_t phi)
static void	EvalParsPhiDependence(Int_t absId, const AliEMCALGeometry* geom)
Bool_t	FindFitV2(AliEMCALRecPoint* emcRP, AliEMCALDigit** MaxAt, const Float_t* maxAtEnergy, Int_t NPar, Float_t* FitParametres) const
AliEMCALUnfolding&	operator=(const AliEMCALUnfolding&)
Bool_t	UnfoldClusterV2(AliEMCALRecPoint* iniEmc, Int_t Nmax, AliEMCALDigit** maxAt, Float_t* maxAtEnergy)
Bool_t	UnfoldClusterV2old(AliEMCALRecPoint* iniEmc, Int_t Nmax, AliEMCALDigit** maxAt, Float_t* maxAtEnergy)

Data Members

public:

enum TObject::EStatusBits {	kCanDelete
	kMustCleanup
	kObjInCanvas
	kIsReferenced
	kHasUUID
	kCannotPick
	kNoContextMenu
	kInvalidObject
};
enum TObject::[unnamed] {	kIsOnHeap
	kNotDeleted
	kZombie
	kBitMask
	kSingleKey
	kOverwrite
	kWriteDelete
};

protected:

TClonesArray*	fDigitsArr	Array with EMCAL digits
Float_t	fECALocMaxCut	minimum energy difference to distinguish local maxima in a cluster
AliEMCALGeometry*	fGeom	! pointer to geometry for utilities
Int_t	fNumberOfECAClusters	number of clusters found in EC section
TObjArray*	fRecPoints	Array with EMCAL clusters
Bool_t	fRejectBelowThreshold	split (false) or reject (true) cell energy below threshold after UF
Float_t	fThreshold	minimum energy for cell to be joined to a cluster

private:

static Double_t	fgPar5[3]	! UF SSPar nr 5 = p0 + phip1 + phi^2 p2
static Double_t	fgPar6[3]	! UF SSPar nr 6 = p0 + phip1 + phi^2 p2
static Double_t	fgSSPars[8]	! Unfolding shower shape parameters

Class Charts

Inheritance Chart:

TObject

←

AliEMCALUnfolding

Function documentation

AliEMCALUnfolding()

 ctor with the indication of the file where header Tree and digits Tree are stored

AliEMCALUnfolding(AliEMCALGeometry* geometry)

 ctor with the indication of the file where header Tree and digits Tree are stored
 use this contructor to avoid usage of Init() which uses runloader
 change needed by HLT - MP

AliEMCALUnfolding(AliEMCALGeometry* geometry, Float_t ECALocMaxCut, Double_t* SSPars, Double_t* Par5, Double_t* Par6)

 ctor with the indication of the file where header Tree and digits Tree are stored
 use this contructor to avoid usage of Init() which uses runloader
 change needed by HLT - MP

void Init()

 Make all memory allocations which can not be done in default constructor.
 Attach the Clusterizer task to the list of EMCAL tasks

~AliEMCALUnfolding()

 dtor

void SetInput(Int_t numberOfECAClusters, TObjArray* recPoints, TClonesArray* digitsArr)

Set input for unfolding purposes

void MakeUnfolding()

 Unfolds clusters using the shape of an ElectroMagnetic shower
 Performs unfolding of all clusters

Int_t UnfoldOneCluster(AliEMCALRecPoint* iniTower, Int_t nMax, AliEMCALDigit** maxAt, Float_t* maxAtEnergy, TObjArray* list)

 Input one cluster
 Output list of clusters
 returns number of clusters
 if fit failed or unfolding is not applicable returns 0 and empty list

Bool_t UnfoldClusterV2(AliEMCALRecPoint* iniEmc, Int_t Nmax, AliEMCALDigit** maxAt, Float_t* maxAtEnergy)

 Extended to whole EMCAL
 Performs the unfolding of a cluster with nMax overlapping showers
 Returns true if success (1->several clusters), otherwise false (fit failed)

Bool_t UnfoldClusterV2old(AliEMCALRecPoint* iniEmc, Int_t Nmax, AliEMCALDigit** maxAt, Float_t* maxAtEnergy)

 Extended to whole EMCAL
 Performs the unfolding of a cluster with nMax overlapping showers

Bool_t FindFitV2(AliEMCALRecPoint* emcRP, AliEMCALDigit** MaxAt, const Float_t* maxAtEnergy, Int_t NPar, Float_t* FitParametres) const

 Calls TMinuit to fit the energy distribution of a cluster with several maxima
 The initial values for fitting procedure are set equal to the
 positions of local maxima.
 Cluster will be fitted as a superposition of nPar/3
 electromagnetic showers

Double_t ShowerShapeV2(Double_t x, Double_t y)

 extended to whole EMCAL
 Shape of the shower
 If you change this function, change also the gradient evaluation in ChiSquare()

void UnfoldingChiSquareV2(Int_t& nPar, Double_t* Grad, Double_t& fret, Double_t* x, Int_t iflag)

 Calculates the Chi square for the cluster unfolding minimization
 Number of parameters, Gradient, Chi squared, parameters, what to do

void SetShowerShapeParams(Double_t* pars)

void SetPar5(Double_t* pars)

void SetPar6(Double_t* pars)

void EvalPar5(Double_t phi)

Evaluate the 5th parameter of the shower shape function
phi in degrees range (-10,10)

fSSPars[5] = 12.31 - phi*0.007381 - phi*phi*0.06936;

void EvalPar6(Double_t phi)

Evaluate the 6th parameter of the shower shape function
phi in degrees range (-10,10)

fSSPars[6] = 0.05452 + phi*0.0001228 + phi*phi*0.001361;

void EvalParsPhiDependence(Int_t absId, const AliEMCALGeometry* geom)

 calculate params p5 and p6 depending on the phi angle in global coordinate
 for the cell with given absId index

AliEMCALUnfolding()

void SetNumberOfECAClusters(Int_t n)

setters and getters

{ fNumberOfECAClusters = n; }

Int_t GetNumberOfECAClusters() const

{ return fNumberOfECAClusters; }

void SetRecPoints(TObjArray* rec)

{ fRecPoints = rec; }

TObjArray * GetRecPoints() const

{ return fRecPoints; }

void SetDigitsArr(TClonesArray* digit)

{ fDigitsArr = digit; }

TClonesArray * GetDigitsArr() const

{ return fDigitsArr; }

void SetECALocalMaxCut(Float_t cut)

{ fECALocMaxCut = cut ; }

Float_t GetECALocalMaxCut() const

{ return fECALocMaxCut; }

void SetThreshold(Float_t energy)

{ fThreshold = energy; }

Float_t GetThreshold() const

{ return fThreshold; }

void SetRejectBelowThreshold(Bool_t reject)

{ fRejectBelowThreshold = reject; }

Bool_t GetRejectBelowThreshold() const

{ return fRejectBelowThreshold; }

Double_t* GetShowerShapeParams() const

{ return fgSSPars ; }

Double_t* GetPar5() const

{ return fgPar5 ; }

Double_t* GetPar6() const

{ return fgPar6 ; }