1 #ifndef ALIDAJETFINDER_H
2 #define ALIDAJETFINDER_H
4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
7 //---------------------------------------------------------------------
8 //Jet finder based on Deterministic Annealing
9 //Author: Davide Perrino (davide.perrino@ba.infn.it)
10 //---------------------------------------------------------------------
12 #include <AliJetFinder.h>
17 class AliDAJetFinder : public AliJetFinder
21 virtual ~AliDAJetFinder();
25 void InitDetAnn (Double_t &dEtSum,Double_t **xData,TVectorD *px,TVectorD *py,TMatrixD *pyx,TMatrixD *y);
26 void Annealing (Int_t nk,Double_t **xData,TVectorD *vPx,TVectorD *vPy,TMatrixD *mPyx,TMatrixD *mY);
27 void NumCl (Int_t &nc,Int_t &nk,TVectorD *vPy,TMatrixD *mPyx,TMatrixD *mY);
28 void ReduceClusters(Int_t **iSame,Int_t nc,Int_t &ncout,Int_t **cont,Int_t *nSameOut);
29 void DoubleClusters(Int_t nc,Int_t &nk,TVectorD *vPy,TMatrixD *mY);
30 void EndDetAnn (Int_t &nk,Double_t **xData,Int_t *xx,Double_t etx,TVectorD *px,TVectorD *py,TMatrixD *pyx,TMatrixD *y);
31 void StoreJets (Int_t nk,Double_t **xData,Int_t *xx,TMatrixD *mY);
34 AliDAJetFinder(const AliDAJetFinder &jf);
35 AliDAJetFinder& operator=(const AliDAJetFinder &jf);
36 Double_t fAlpha; // beta increment
37 Double_t fDelta; // perturbation proportional to Delta
38 Double_t fAvDist; // minimum distance to distinguish two clusters
39 Double_t fEps; // convergence criterium below max number of loops
40 Double_t fEpsMax; // convergence criterium above max number of loops
41 Int_t fNloopMax; // maximum number of loops at a fixed beta
42 Double_t fBeta; // increasing multiplier of entropy
43 Int_t fNclustMax; // maximum number of clusters to find
44 Int_t fNin; // number of input data
46 ClassDef(AliDAJetFinder,2)