]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - JETAN/AliDAJetFinder.h
Updated macro to plot SDD injector performance vs. time
[u/mrichter/AliRoot.git] / JETAN / AliDAJetFinder.h
index 82a79e4619ea3fd1673fc77a9664366450cdbe83..4e5a0c65ada13d5cabebc7beaf61c1d646f9599a 100644 (file)
@@ -19,18 +19,16 @@ class AliDAJetFinder : public AliJetFinder
 public:
     AliDAJetFinder();
     virtual  ~AliDAJetFinder();
-    
+
     void FindJets      ();
-    void SetJetHeader  (AliDAJetHeader *h) {fHeader = h;}
-    
  private:
-    void InitDetAnn    (Double_t &dEtSum);
-    void Annealing     (Int_t nk);
-    void NumCl         (Int_t &nc,Int_t &nk);
-    void ReduceClusters(Int_t **iSame,Int_t nc,Int_t &ncout,Int_t **cont,Int_t *nSameOut);
-    void DoubleClusters(Int_t nc,Int_t &nk);
-    void EndDetAnn     (Int_t &nk,Int_t *xx,Double_t etx);
-    void StoreJets     (Int_t nk,Int_t *xx);
+    void InitDetAnn    (Double_t &dEtSum, Double_t **xData, TVectorD *px, TVectorD *py, TMatrixD *pyx, TMatrixD *y);
+    void Annealing     (Int_t nk, Double_t **xData, TVectorD *vPx, TVectorD *vPy, TMatrixD *mPyx, TMatrixD *mY);
+    void NumCl         (Int_t &nc, Int_t &nk, TVectorD *vPy, TMatrixD *mPyx, TMatrixD *mY);
+    void ReduceClusters(Int_t **iSame, Int_t nc, Int_t &ncout, Int_t **cont, Int_t *nSameOut) const;
+    void DoubleClusters(Int_t nc, Int_t &nk,  TVectorD *vPy,  TMatrixD *mY) const;
+    void EndDetAnn     (Int_t &nk, Double_t **xData, Int_t *xx, Double_t etx, TVectorD *px, TVectorD *py, TMatrixD *pyx, TMatrixD *y);
+    void StoreJets     (Int_t nk, Double_t **xData, Int_t *xx, TMatrixD *mY);
 
 protected:
     AliDAJetFinder(const AliDAJetFinder &jf);
@@ -43,16 +41,9 @@ protected:
     Int_t      fNloopMax;                              // maximum number of loops at a fixed beta
     Double_t   fBeta;                                  // increasing multiplier of entropy
     Int_t      fNclustMax;                             // maximum number of clusters to find
-    TMatrixD  *fPyx;                                   // conditional probability matrix
-    TMatrixD  *fY;                                      // clusters matrix
-    TVectorD  *fPx;                                     // input data weights
-    TVectorD  *fPy;                                     // clusters' probability
-    Double_t  *fXEta;                                  // input data eta coordinate
-    Double_t  *fXPhi;                                  // input data phi coordinate
     Int_t      fNin;                                   // number of input data
-     
-    AliDAJetHeader *fHeader;                            // the header  
-    
-    ClassDef(AliDAJetFinder,1)
-};// 
+    Int_t      fNeff;                                  // total input data, including fakes
+
+    ClassDef(AliDAJetFinder,3)
+};
 #endif