[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$ */

////////////////////////////////////////////////
//  MUON Cluster Finder Class                 //
////////////////////////////////////////////////
#include <TObject.h>
#include "AliMUONHitMapA1.h"
#include "AliMUONClusterInput.h"
#include "TF1.h"

class AliMUONResponse;
class AliMUONSegmentation;
class TClonesArray;
class AliMUONRawCluster;
class AliMUONDigit;


class AliMUONClusterFinderVS : public TObject 
{
 public:
    AliMUONClusterFinderVS();
    AliMUONClusterFinderVS(const AliMUONClusterFinderVS& clusterFinder);
    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(const AliMUONRawCluster cluster);
//  Set tracks for debugging    
    virtual void SetTracks(Int_t t1, Int_t t2) {fTrack[0]=t1; fTrack[1]=t2;}
    virtual Bool_t TestTrack(Int_t t);
//  Assignment operator
    AliMUONClusterFinderVS & operator = (const AliMUONClusterFinderVS& rhs);

 protected:
    AliMUONClusterInput*    fInput;              // AliMUONClusterInput instance
    AliMUONHitMapA1*        fHitMap[2];          // Hit Map cathode 1
// Configuration    
    Int_t                   fDeclusterFlag;      // flag for declusterin
    Int_t                   fClusterSize;        // cluster size 
    Int_t                   fNperMax;            // Maximum number of pads per peak
// Current decluster result    
    Int_t                   fMul[2];             // current multiplicity
    Int_t                   fNPeaks;             // number of local maxima
    Int_t                   fNRawClusters;       // Number of Raw Clusters
// 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.
    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
// 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,1) //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$ */
	7
	8	////////////////////////////////////////////////
	9	// MUON Cluster Finder Class //
	10	////////////////////////////////////////////////
30aaba74	11	#include <TObject.h>
	12	#include "AliMUONHitMapA1.h"
	13	#include "AliMUONClusterInput.h"
a9e2aefa	14	#include "TF1.h"
9825400f	15
30aaba74	16	class AliMUONResponse;
	17	class AliMUONSegmentation;
	18	class TClonesArray;
	19	class AliMUONRawCluster;
	20	class AliMUONDigit;
	21
	22
	23	class AliMUONClusterFinderVS : public TObject
a9e2aefa	24	{
a9e2aefa	25	public:
a9e2aefa	26	AliMUONClusterFinderVS();
	27	AliMUONClusterFinderVS(const AliMUONClusterFinderVS& clusterFinder);
	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
a9e2aefa	55	virtual void AddRawCluster(const AliMUONRawCluster cluster);
30aaba74	56	// Set tracks for debugging
	57	virtual void SetTracks(Int_t t1, Int_t t2) {fTrack[0]=t1; fTrack[1]=t2;}
	58	virtual Bool_t TestTrack(Int_t t);
	59	// Assignment operator
a9e2aefa	60	AliMUONClusterFinderVS & operator = (const AliMUONClusterFinderVS& rhs);
30aaba74	61
	62	protected:
	63	AliMUONClusterInput* fInput; // AliMUONClusterInput instance
	64	AliMUONHitMapA1* fHitMap[2]; // Hit Map cathode 1
	65	// Configuration
	66	Int_t fDeclusterFlag; // flag for declusterin
	67	Int_t fClusterSize; // cluster size
	68	Int_t fNperMax; // Maximum number of pads per peak
	69	// Current decluster result
	70	Int_t fMul[2]; // current multiplicity
	71	Int_t fNPeaks; // number of local maxima
	72	Int_t fNRawClusters; // Number of Raw Clusters
	73	// Local data store
a9e2aefa	74	AliMUONDigit* fDig[100][2]; // current list of digits
	75	Int_t fIx[100][2]; // current list of x-pad-coord.
	76	Int_t fIy[100][2]; // current list of y-pad-coord.
	77	Float_t fX[100][2]; // current list of x-coord.
	78	Float_t fY[100][2]; // current list of y-coord.
	79	Int_t fIndLocal[100][2]; // indices of local maxima
	80	Int_t fNLocal[2]; // Number of local maxima
	81	Int_t fQ[100][2]; // current list of charges
a9e2aefa	82	// Current Fit
	83	Double_t fXFit[2]; // x-coordinate
	84	Double_t fYFit[2]; // y-coordinate
	85	Double_t fQrFit[2]; // charge ratio
	86	Float_t fChi2[2]; // chi2 of fit
	87	Float_t fXInit[2]; // start values
	88	Float_t fYInit[2]; // start values
	89	Float_t fQrInit[2]; // start values
	90	Int_t fFitStat; // status of fit
	91
	92	// Selected track for debugging
	93	Int_t fTrack[2]; // Only digits with main contributions from these tracks are
	94	// considered
	95	// Return pointer to raw clusters
	96	ClassDef(AliMUONClusterFinderVS,1) //Class for clustering and reconstruction of space points
	97	};
	98	#endif
	99
	100
	101
	102
	103
	104
	105
	106
	107
	108
30aaba74	109
	110
	111
	112
	113