[u/mrichter/AliRoot.git] / MUON / AliMUONClusterFinderVS.h

#ifndef ALIMUONCLUSTERFINDERVS_H
#define ALIMUONCLUSTERFINDERVS_H
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

/* $Id$ */
// Revision of includes 07/05/2004

/// \ingroup rec
/// \class AliMUONClusterFinderVS
/// \brief Class for clustering and reconstruction of space points

#include <TObject.h>

class TClonesArray;

class AliMUONClusterInput;
class AliMUONDigitMapA1;
class AliMUONGeometrySegmentation;
class AliMUONRawCluster;
class AliMUONDigit;


class AliMUONClusterFinderVS : public TObject 
{
 public:
    AliMUONClusterFinderVS();
    virtual ~AliMUONClusterFinderVS();
// Decluster ?
    virtual void SetDeclusterFlag(Int_t flag=1) {fDeclusterFlag =flag;}
// Set max. cluster size ; bigger clusters will deconvoluted
    virtual void SetClusterSize(Int_t clsize=5) {fClusterSize = clsize;}
// Set max. number of pads per local cluster
    virtual void SetNperMax(Int_t npermax=5) {fNperMax = npermax;}
// Search for raw clusters
    virtual void  FindRawClusters();
// Find cluster    
    virtual void  FindCluster(Int_t i, Int_t j, Int_t cath, AliMUONRawCluster &c);
// Decluster
    virtual void Decluster(AliMUONRawCluster *cluster);
//  Perform split by local maxima  
    virtual void   SplitByLocalMaxima(AliMUONRawCluster *cluster);
    virtual void   FindLocalMaxima(AliMUONRawCluster *cluster);
    virtual void   Split(AliMUONRawCluster * cluster);
//  Perform Double Mathieson Fit
    Bool_t  DoubleMathiesonFit(AliMUONRawCluster *c, Int_t cath);
    Float_t CombiDoubleMathiesonFit(AliMUONRawCluster *c);
    Float_t SingleMathiesonFit(AliMUONRawCluster *c, Int_t cath);
    Float_t CombiSingleMathiesonFit(AliMUONRawCluster *c);    
//  Build up full cluster information    
    virtual void   FillCluster(AliMUONRawCluster *cluster, Int_t flag, Int_t cath);
    virtual void   FillCluster(AliMUONRawCluster *cluster, Int_t cath);
    virtual void   FillCluster(AliMUONRawCluster *cluster) {FillCluster(cluster,1,0);}
// Add a new raw cluster    
    virtual void AddRawCluster(AliMUONRawCluster& cluster);
//  Set tracks for debugging    
    virtual void SetTracks(Int_t t1, Int_t t2) {fTrack[0]=t1; fTrack[1]=t2;}
    void SetGhostChi2Cut(Float_t cut) {fGhostChi2Cut = cut;}
// get raw cluster pointer 
    TClonesArray*  GetRawClusters() {return fRawClusters;}
// reset raw clusters
    void ResetRawClusters();

 protected:
    AliMUONClusterFinderVS(const AliMUONClusterFinderVS& clusterFinder);
//  Assignment operator
    AliMUONClusterFinderVS & operator = (const AliMUONClusterFinderVS& rhs);

    AliMUONClusterInput*          fInput;        //!< AliMUONClusterInput instance
    AliMUONDigitMapA1*            fDigitMap[2];  ///< Hit Maps for cathode 1 and 2
    AliMUONGeometrySegmentation*  fSeg2[2];      ///< New Segmentations for cathode 1 and 2

    
// Configuration    
    Int_t                   fDeclusterFlag;      ///< flag for declusterin
    Int_t                   fClusterSize;        ///< cluster size 
    Int_t                   fNperMax;            ///< Maximum number of pads per peak
    Float_t                 fGhostChi2Cut;       ///< \brief Cut in charge matching chi2
	                                         /// (2 degrees of freedom)
                                                 /// Used by ghost removal
    // Current decluster result    
    Int_t                   fMul[2];             ///< current multiplicity
    Int_t                   fNPeaks;             ///< number of local maxima
    Int_t                   fNRawClusters;       ///< Number of Raw Clusters
    TClonesArray*           fRawClusters;        ///< array of cluster per ch.

    // Local data store    
    AliMUONDigit*           fDig[100][2];        ///< current list of digits 
    Int_t                   fIx[100][2];         ///< current list of x-pad-coord.
    Int_t                   fIy[100][2];         ///< current list of y-pad-coord.
    Float_t                 fX[100][2];          ///< current list of x-coord.
    Float_t                 fY[100][2];          ///< current list of y-coord.
    Float_t                 fZ[100][2];          ///< current list of z-coord.
    Int_t                   fIndLocal[100][2];   ///< indices of local maxima
    Int_t                   fNLocal[2];          ///< Number of local maxima
    Int_t                   fQ[100][2];          ///< current list of charges
    Float_t                 fZPlane;             ///< currenz z-plane position
    Int_t                   fSector;             ///< current sector
    
    // Current Fit
    Double_t                 fXFit[2];         ///< x-coordinate
    Double_t                 fYFit[2];         ///< y-coordinate
    Double_t                 fQrFit[2];        ///< charge ratio
    Float_t                  fChi2[2];         ///< chi2 of fit
    Float_t                  fXInit[2];        ///< start values
    Float_t                  fYInit[2];        ///< start values
    Float_t                  fQrInit[2];       ///< start values
    Int_t                    fFitStat;         ///< status of fit
    
    // Selected track for debugging
    Int_t                    fTrack[2];        ///< Only digits with main contributions from these tracks are
    // considered 
    
    //  Return pointer to raw clusters    
    ClassDef(AliMUONClusterFinderVS,2) //Class for clustering and reconstruction of space points
      };
#endif
Commit	Line	Data
a9e2aefa	1	#ifndef ALIMUONCLUSTERFINDERVS_H
	2	#define ALIMUONCLUSTERFINDERVS_H
	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
	5
	6	/* $Id$ */
30178c30	7	// Revision of includes 07/05/2004
a9e2aefa	8
692de412	9	/// \ingroup rec
	10	/// \class AliMUONClusterFinderVS
	11	/// \brief Class for clustering and reconstruction of space points
c1a185bf	12
30aaba74	13	#include <TObject.h>
9825400f	14
30178c30	15	class TClonesArray;
30178c30	16
c1a185bf	17	class AliMUONClusterInput;
5289cf2f	18	class AliMUONDigitMapA1;
fed772f3	19	class AliMUONGeometrySegmentation;
30aaba74	20	class AliMUONRawCluster;
	21	class AliMUONDigit;
	22
	23
	24	class AliMUONClusterFinderVS : public TObject
a9e2aefa	25	{
a9e2aefa	26	public:
a9e2aefa	27	AliMUONClusterFinderVS();
4da78c65	28	virtual ~AliMUONClusterFinderVS();
30aaba74	29	// Decluster ?
	30	virtual void SetDeclusterFlag(Int_t flag=1) {fDeclusterFlag =flag;}
	31	// Set max. cluster size ; bigger clusters will deconvoluted
	32	virtual void SetClusterSize(Int_t clsize=5) {fClusterSize = clsize;}
	33	// Set max. number of pads per local cluster
	34	virtual void SetNperMax(Int_t npermax=5) {fNperMax = npermax;}
a9e2aefa	35	// Search for raw clusters
30aaba74	36	virtual void FindRawClusters();
a9e2aefa	37	// Find cluster
	38	virtual void FindCluster(Int_t i, Int_t j, Int_t cath, AliMUONRawCluster &c);
	39	// Decluster
	40	virtual void Decluster(AliMUONRawCluster *cluster);
	41	// Perform split by local maxima
	42	virtual void SplitByLocalMaxima(AliMUONRawCluster *cluster);
	43	virtual void FindLocalMaxima(AliMUONRawCluster *cluster);
	44	virtual void Split(AliMUONRawCluster * cluster);
a9e2aefa	45	// Perform Double Mathieson Fit
30aaba74	46	Bool_t DoubleMathiesonFit(AliMUONRawCluster *c, Int_t cath);
a9e2aefa	47	Float_t CombiDoubleMathiesonFit(AliMUONRawCluster *c);
	48	Float_t SingleMathiesonFit(AliMUONRawCluster *c, Int_t cath);
	49	Float_t CombiSingleMathiesonFit(AliMUONRawCluster *c);
	50	// Build up full cluster information
	51	virtual void FillCluster(AliMUONRawCluster *cluster, Int_t flag, Int_t cath);
	52	virtual void FillCluster(AliMUONRawCluster *cluster, Int_t cath);
30aaba74	53	virtual void FillCluster(AliMUONRawCluster *cluster) {FillCluster(cluster,1,0);}
30aaba74	54	// Add a new raw cluster
b137f8b9	55	virtual void AddRawCluster(AliMUONRawCluster& cluster);
30aaba74	56	// Set tracks for debugging
30aaba74	57	virtual void SetTracks(Int_t t1, Int_t t2) {fTrack[0]=t1; fTrack[1]=t2;}
07cfabcf	58	void SetGhostChi2Cut(Float_t cut) {fGhostChi2Cut = cut;}
4da78c65	59	// get raw cluster pointer
	60	TClonesArray* GetRawClusters() {return fRawClusters;}
	61	// reset raw clusters
	62	void ResetRawClusters();
30aaba74	63
30aaba74	64	protected:
30178c30	65	AliMUONClusterFinderVS(const AliMUONClusterFinderVS& clusterFinder);
	66	// Assignment operator
	67	AliMUONClusterFinderVS & operator = (const AliMUONClusterFinderVS& rhs);
	68
829425a5	69	AliMUONClusterInput* fInput; //!< AliMUONClusterInput instance
	70	AliMUONDigitMapA1* fDigitMap[2]; ///< Hit Maps for cathode 1 and 2
	71	AliMUONGeometrySegmentation* fSeg2[2]; ///< New Segmentations for cathode 1 and 2
fed772f3	72
f0d86bc4	73
30aaba74	74	// Configuration
829425a5	75	Int_t fDeclusterFlag; ///< flag for declusterin
	76	Int_t fClusterSize; ///< cluster size
	77	Int_t fNperMax; ///< Maximum number of pads per peak
	78	Float_t fGhostChi2Cut; ///< \brief Cut in charge matching chi2
	79	/// (2 degrees of freedom)
	80	/// Used by ghost removal
b137f8b9	81	// Current decluster result
829425a5	82	Int_t fMul[2]; ///< current multiplicity
	83	Int_t fNPeaks; ///< number of local maxima
	84	Int_t fNRawClusters; ///< Number of Raw Clusters
	85	TClonesArray* fRawClusters; ///< array of cluster per ch.
4da78c65	86
b137f8b9	87	// Local data store
829425a5	88	AliMUONDigit* fDig[100][2]; ///< current list of digits
	89	Int_t fIx[100][2]; ///< current list of x-pad-coord.
	90	Int_t fIy[100][2]; ///< current list of y-pad-coord.
	91	Float_t fX[100][2]; ///< current list of x-coord.
	92	Float_t fY[100][2]; ///< current list of y-coord.
	93	Float_t fZ[100][2]; ///< current list of z-coord.
	94	Int_t fIndLocal[100][2]; ///< indices of local maxima
	95	Int_t fNLocal[2]; ///< Number of local maxima
	96	Int_t fQ[100][2]; ///< current list of charges
	97	Float_t fZPlane; ///< currenz z-plane position
	98	Int_t fSector; ///< current sector
f0d86bc4	99
b137f8b9	100	// Current Fit
829425a5	101	Double_t fXFit[2]; ///< x-coordinate
	102	Double_t fYFit[2]; ///< y-coordinate
	103	Double_t fQrFit[2]; ///< charge ratio
	104	Float_t fChi2[2]; ///< chi2 of fit
	105	Float_t fXInit[2]; ///< start values
	106	Float_t fYInit[2]; ///< start values
	107	Float_t fQrInit[2]; ///< start values
	108	Int_t fFitStat; ///< status of fit
a9e2aefa	109
b137f8b9	110	// Selected track for debugging
829425a5	111	Int_t fTrack[2]; ///< Only digits with main contributions from these tracks are
b137f8b9	112	// considered
b137f8b9	113
b137f8b9	114	// Return pointer to raw clusters
4da78c65	115	ClassDef(AliMUONClusterFinderVS,2) //Class for clustering and reconstruction of space points
b137f8b9	116	};
a9e2aefa	117	#endif
	118
	119
	120
	121
	122
	123
	124
	125
	126
	127
30aaba74	128
	129
	130
	131
	132