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

AliRoot » MUON » AliMUONClusterFinderMLEM

class AliMUONClusterFinderMLEM: public AliMUONVClusterFinder

 \class AliMUONClusterFinderMLEM

 Clusterizer class based on the Expectation-Maximization algorithm

 Pre-clustering is handled by AliMUONPreClusterFinder
 From a precluster a pixel array is built, and from this pixel array
 a list of clusters is output (using AliMUONClusterSplitterMLEM).

 \author Laurent Aphecetche (for the "new" C++ structure) and
 Alexander Zinchenko, JINR Dubna, for the hardcore of it ;-)

Function Members (Methods)

public:

	AliMUONClusterFinderMLEM(Bool_t plot, AliMUONVClusterFinder* clusterFinder)
virtual	~AliMUONClusterFinderMLEM()
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
static Int_t	GetCoupledFlag()
virtual Option_t*	TObject::GetDrawOption() const
static Long_t	TObject::GetDtorOnly()
virtual const char*	TObject::GetIconName() const
static Int_t	GetModifiedFlag()
static Int_t	GetMustKeepFlag()
virtual const char*	TObject::GetName() const
virtual char*	TObject::GetObjectInfo(Int_t px, Int_t py) const
static Bool_t	TObject::GetObjectStat()
virtual Option_t*	TObject::GetOption() const
static Int_t	GetOverFlag()
virtual const char*	TObject::GetTitle() const
virtual UInt_t	TObject::GetUniqueID() const
static Int_t	GetUseForFitFlag()
static Int_t	GetZeroFlag()
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	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
void	TObject::MayNotUse(const char* method) const
virtual Bool_t	NeedSegmentation() const
virtual AliMUONCluster*	NextCluster()
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	Paint(Option_t* opt = "")
virtual void	TObject::Pop()
virtual Bool_t	Prepare(Int_t detElemId, TObjArray pads, const AliMpArea& area, const AliMpVSegmentation segmentations)
virtual void	Print(Option_t* opt = "") 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	SetChargeHints(Double_t lowestPadCharge, Double_t lowestClusterCharge)
virtual void	TObject::SetDrawOption(Option_t* option = "")MENU
static void	TObject::SetDtorOnly(void* obj)
static void	TObject::SetObjectStat(Bool_t stat)
virtual void	TObject::SetUniqueID(UInt_t uid)
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
virtual void	TObject::UseCurrentStyle()
virtual Bool_t	AliMUONVClusterFinder::UsePad(const AliMUONPad& pad)
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:

	AliMUONClusterFinderMLEM(const AliMUONClusterFinderMLEM& rhs)
void	AddBinSimple(TH2D* mlem, Int_t ic, Int_t jc)
void	AddVirtualPad(AliMUONCluster& cluster)
void	BuildPixArray(AliMUONCluster& cluster)
void	BuildPixArrayOneCathode(AliMUONCluster& cluster)
void	CheckOverlaps()
AliMUONCluster*	CheckPrecluster(const AliMUONCluster& cluster)
AliMUONCluster*	CheckPreclusterTwoCathodes(AliMUONCluster* cluster)
void	ComputeCoefficients(AliMUONCluster& cluster, Double_t* coef, Double_t* probi)
void	FindCluster(AliMUONCluster& cluster, const Int_t* localMax, Int_t iMax)
void	FindCOG(Double_t* xyc)
Int_t	FindLocalMaxima(TObjArray* pixArray, Int_t* localMax, Double_t* maxVal)
Int_t	FindNearest(const AliMUONPad* pixPtr0)
void	FlagLocalMax(TH2D* hist, Int_t i, Int_t j, Int_t* isLocalMax)
Bool_t	MainLoop(AliMUONCluster& cluster, Int_t iSimple)
void	MaskPeaks(Int_t mask)
void	Mlem(AliMUONCluster& cluster, const Double_t* coef, Double_t* probi, Int_t nIter)
AliMUONClusterFinderMLEM&	operator=(const AliMUONClusterFinderMLEM& rhs)
Bool_t	Overlap(const AliMUONPad& pad, const AliMUONPad& pixel)
void	PadOverHist(Int_t idir, Int_t ix0, Int_t iy0, AliMUONPad* pad, TH2D* hist1, TH2D* hist2)
void	PadsInXandY(AliMUONCluster& cluster, Int_t& nInX, Int_t& nInY) const
AliMUONPad*	Pixel(Int_t i) const
void	Plot(const char* outputfile)
void	RemovePixel(Int_t i)
void	Simple(AliMUONCluster& cluster)
Bool_t	WorkOnPreCluster()

Data Members

public:

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

private:

TObjArray	fClusterList	!< clusters corresponding to the current pre-cluster
Int_t	fClusterNumber	!< current cluster number
Int_t	fDebug	!< debug level
Int_t	fDetElemId	!< current DE being processed
Int_t	fEventNumber	!< current event being processed
TH2D*	fHistAnode	!< histogram for local maxima search
TH2D*	fHistMlem	!< histogram for MLEM procedure
Double_t	fLowestClusterCharge	!< see AliMUONRecoParam
Double_t	fLowestPadCharge	!< see AliMUONRecoParam
Double_t	fLowestPixelCharge	!< see AliMUONRecoParam
Int_t	fNAddVirtualPads	!< number of clusters for which we added virtual pads
Int_t	fNClusters	!< total number of clusters
TObjArray*	fPixArray	!< collection of pixels
Bool_t	fPlot	!< whether we should plot thing (for debug only, quite slow!)
AliMUONCluster*	fPreCluster	!< current pre-cluster
AliMUONVClusterFinder*	fPreClusterFinder	!< the pre-clustering worker
AliMUONClusterSplitterMLEM*	fSplitter	!< helper class to go from pixel arrays to clusters
static const Int_t	fgkCoupled	/< coupled pad
static const TVector2	fgkDecreaseSize	/< idem
static const Double_t	fgkDistancePrecision	/< used to check overlaps and so on
static const TVector2	fgkIncreaseSize	/< idem
static const Int_t	fgkModified	/< modified pad charge
static const Int_t	fgkMustKeep	/< do not kill (for pixels)
static const Int_t	fgkOver	/< processing is over
static const Int_t	fgkUseForFit	/< should be used for fit
static const Int_t	fgkZero	/< pad "basic" state
const AliMpVSegmentation*	fkSegmentation[2]	!< new segmentation

Class Charts

Inheritance Chart:

TObject

←

AliMUONVClusterFinder

←

AliMUONClusterFinderMLEM

Function documentation

AliMUONClusterFinderMLEM(Bool_t plot, AliMUONVClusterFinder* clusterFinder)

 Constructor

~AliMUONClusterFinderMLEM()

 Destructor

Prepare(Int_t detElemId, TObjArray** pads, const AliMpArea& area, const AliMpVSegmentation** segmentations)

 Prepare for clustering
  AliCodeTimerAuto("",0)

NextCluster()

 Return next cluster
  AliCodeTimerAuto("",0)

WorkOnPreCluster()

 Starting from a precluster, builds a pixel array, and then
 extract clusters from this array

Overlap(const AliMUONPad& pad, const AliMUONPad& pixel)

 Check if the pad and the pixel overlaps

CheckPrecluster(const AliMUONCluster& cluster)

 Check precluster in order to attempt to simplify it (mostly for
 two-cathode preclusters)

CheckPreclusterTwoCathodes(AliMUONCluster* cluster)

 Check two-cathode cluster

CheckOverlaps()

 For debug only : check if some pixels overlap...

void BuildPixArray(AliMUONCluster& cluster)

 Build pixel array for MLEM method

void BuildPixArrayOneCathode(AliMUONCluster& cluster)

 Build the pixel array

void PadOverHist(Int_t idir, Int_t ix0, Int_t iy0, AliMUONPad* pad, TH2D* hist1, TH2D* hist2)

 "Span" pad over histogram in the direction idir

Plot(const char* outputfile)

 Make a plot and save it as png

ComputeCoefficients(AliMUONCluster& cluster, Double_t* coef, Double_t* probi)

 Compute coefficients needed for MLEM algorithm

Bool_t MainLoop(AliMUONCluster& cluster, Int_t iSimple)

 Repeat MLEM algorithm until pixel size becomes sufficiently small

void MaskPeaks(Int_t mask)

 Mask/unmask pixels corresponding to local maxima (add/subtract 10000 to their charge
 - to avoid loosing low charge pixels after sorting)

void Mlem(AliMUONCluster& cluster, const Double_t* coef, Double_t* probi, Int_t nIter)

 Use MLEM to find pixel charges

void FindCOG(Double_t* xyc)

 Calculate position of the center-of-gravity around the maximum pixel

Int_t FindNearest(const AliMUONPad* pixPtr0)

 Find the pixel nearest to the given one
 (algorithm may be not very efficient)

Paint(Option_t* opt = "")

 Paint cluster and pixels

Int_t FindLocalMaxima(TObjArray* pixArray, Int_t* localMax, Double_t* maxVal)

 Find local maxima in pixel space for large preclusters in order to
 try to split them into smaller pieces (to speed up the MLEM procedure)
 or to find additional fitting seeds if clusters were not completely resolved

void FlagLocalMax(TH2D* hist, Int_t i, Int_t j, Int_t* isLocalMax)

 Flag pixels (whether or not local maxima)

void FindCluster(AliMUONCluster& cluster, const Int_t* localMax, Int_t iMax)

 Find pixel cluster around local maximum \a iMax and pick up pads
 overlapping with it

AddBinSimple(TH2D* mlem, Int_t ic, Int_t jc)

 Add adjacent bins (+-1 in X and Y) to the cluster

operator=(const AliMUONClusterFinderMLEM& rhs)

 Protected assignement operator

void AddVirtualPad(AliMUONCluster& cluster)

 Add virtual pad (with small charge) to improve fit for clusters
 with number of pads == 2 per direction

void PadsInXandY(AliMUONCluster& cluster, Int_t& nInX, Int_t& nInY) const

 Find number of pads in X and Y-directions (excluding virtual ones and
 overflows)

void RemovePixel(Int_t i)

 Remove pixel at index i

void Simple(AliMUONCluster& cluster)

 Process simple cluster (small number of pads) without EM-procedure

Pixel(Int_t i) const

 Returns pixel at index i

Print(Option_t* opt = "") const

 printout

SetChargeHints(Double_t lowestPadCharge, Double_t lowestClusterCharge)

 Set some thresholds we use later on in the algorithm

AliMUONClusterFinderMLEM(Bool_t plot, AliMUONVClusterFinder* clusterFinder)

Bool_t NeedSegmentation() const

{ return kTRUE; }

Int_t GetZeroFlag()

 Status flags for pads
 Return pad "basic" state flag

{ return fgkZero; }

Int_t GetMustKeepFlag()

 Return do not kill flag

{ return fgkMustKeep; }

Int_t GetUseForFitFlag()

 Return should be used for fit flag

{ return fgkUseForFit; }

Int_t GetOverFlag()

 Return processing is over flag

{ return fgkOver; }

Int_t GetModifiedFlag()

 Return modified pad charge flag

{ return fgkModified; }

Int_t GetCoupledFlag()

 Return coupled pad flag

{ return fgkCoupled; }