1 #ifndef ALIMUONCLUSTERSPLITTERMLEM_H
2 #define ALIMUONCLUSTERSPLITTERMLEM_H
4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
10 /// \class AliMUONClusterSplitterMLEM
11 /// \brief Splitter class for the MLEM algorithm
13 // Author Alexander Zinchenko, JINR Dubna; Laurent Aphecetche, SUBATECH
20 #include "TMatrixDfwd.h"
26 class AliMUONMathieson;
28 class AliMUONClusterSplitterMLEM : public TObject
31 AliMUONClusterSplitterMLEM(Int_t detElemId, TObjArray* pixArray);
33 virtual ~AliMUONClusterSplitterMLEM();
35 void AddBin(TH2 *mlem,
36 Int_t ic, Int_t jc, Int_t mode,
37 Bool_t *used, TObjArray *pix);
39 void AddCluster(Int_t ic, Int_t nclust,
41 Bool_t *used, Int_t *clustNumb, Int_t &nCoupled);
43 TObject* BinToPix(TH2 *mlem, Int_t jc, Int_t ic);
45 Float_t ChargeIntegration(Double_t x, Double_t y, const AliMUONPad& pad);
47 void Fcn1(const AliMUONCluster& cluster,
48 Int_t & npar, Double_t * gin,
49 Double_t &f, Double_t *par, Int_t iflag);
51 Int_t Fit(const AliMUONCluster& cluster,
52 Int_t iSimple, Int_t nfit,
53 Int_t *clustFit, TObjArray **clusters,
54 Double_t *parOk, TObjArray& clusterList, TH2 *mlem);
56 void Merge(const AliMUONCluster& cluster,
57 Int_t nForFit, Int_t nCoupled,
58 Int_t *clustNumb, Int_t *clustFit,
60 TMatrixD& aijcluclu, TMatrixD& aijclupad);
62 Double_t MinGroupCoupl(Int_t nCoupled, Int_t *clustNumb,
63 TMatrixD& aijcluclu, Int_t *minGroup);
65 Int_t SelectPad(const AliMUONCluster& cluster,
66 Int_t nCoupled, Int_t nForFit,
67 Int_t *clustNumb, Int_t *clustFit,
70 void Split(const AliMUONCluster& cluster,
72 Double_t* coef, TObjArray& clusterList);
75 void UpdatePads(const AliMUONCluster& cluster, Int_t nfit, Double_t *par);
77 void SetDebug (Int_t debug) { fDebug = debug; }
80 /// will not be implemented
81 AliMUONClusterSplitterMLEM(const AliMUONClusterSplitterMLEM&);
82 /// will not be implemented
83 AliMUONClusterSplitterMLEM& operator=(const AliMUONClusterSplitterMLEM&);
84 Double_t Param2Coef(Int_t icand, Double_t coef, Double_t *par);
88 static const Double_t fgkCouplMin; ///< threshold on coupling
90 TObjArray* fPixArray; //!< \todo add comment
91 AliMUONMathieson* fMathieson; //!< Mathieson
92 Int_t fDetElemId; //!< detection element we are working on
93 Int_t fNpar; //!< number of fit parameters
94 Double_t fQtot; //!< total charge
95 Int_t fnCoupled; //!< number of coupled pixels ?
96 Int_t fDebug; //!< debug level
98 ClassDef(AliMUONClusterSplitterMLEM,1) //