Add a protection to avoid crash in QA (Julian)

[u/mrichter/AliRoot.git] / JETAN / AliCdfJetFinder.h

diff --git a/JETAN/AliCdfJetFinder.h b/JETAN/AliCdfJetFinder.h
index 7d9842b39e08cfbe7c263c92195baca6d167cfc0..9f05e53c36f4f8508205f4abedb9d55e2e84f1c8 100644 (file)
--- a/JETAN/AliCdfJetFinder.h
+++ b/JETAN/AliCdfJetFinder.h
@@ -1,17 +1,28 @@
 #ifndef ALICDFJETFINDER_H
 #define ALICDFJETFINDER_H
 
-/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
- * See cxx source for full Copyright notice                               */
+/*
+ * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved.
+ * See cxx source for full Copyright notice
+ *
+*/
 
-
-//---------------------------------------------------------------------
-//
-//
-//---------------------------------------------------------------------
+//  Definition of constants, structures and functions for jet finder
 
 #include "AliJetFinder.h"
 
+using namespace std ;
+
+
+// structure of jet and particles container
+struct varContainer
+  {
+  Double_t  pt;  // pt of particle/jet
+  Double_t eta;  // eta of particle/jet
+  Double_t phi;  // phi of particle/jet
+  Int_t    njet; // njet is multiplicity of jet or if particles are stored , njet is index number of jet 
+  } ;
+
 class AliCdfJetHeader;
 
 class AliCdfJetFinder : public AliJetFinder
@@ -21,15 +32,48 @@ class AliCdfJetFinder : public AliJetFinder
     AliCdfJetFinder();
     virtual ~AliCdfJetFinder();
 
-    void           CreateOutputObjects(TList *histos);
+    void           CreateOutputObjects(TList * const histos);
     void           FindJets();
+    void           InitData();
+    void           FindCones();
+    void           ComputeConesWeight();
+    void           WriteJets() ;
+    void           AnalizeJets();
+    void           Clean();
+    
     virtual void   FinishRun();
 
+    void           SetAnalyseJets(Bool_t flag = kTRUE) {fAnalyseJets = flag;}
+    
   protected:
     AliCdfJetFinder ( const AliCdfJetFinder& jf );
     AliCdfJetFinder& operator = ( const AliCdfJetFinder& jf );
 
     TList         *fHistos;    // List of histograms
-    ClassDef ( AliCdfJetFinder, 1 )
-  };//
+
+    Bool_t fFromAod ;          // is the data taken from AOD (or ESD filter) 
+    Bool_t fAODwrite ;         // write jets to AOD
+    Bool_t fAODtracksWrite ;   // write jet tracks to AOD
+    Bool_t fAnalyseJets;       // analyse jets
+       
+    TRefArray *fRefArr ;       //! pointer to references array of tracks from AOD
+
+    Int_t         fNJets;     //! counter of number of jets
+    Int_t         fNPart;     //! number of particles in event
+
+    Double_t      fRadius ;   // radius of jet 
+
+    Int_t fMinJetParticles;   // leading jet must contain AT LEAST fMinJetParticles
+    Double_t fJetPtCut;       // leading jet must have AT LEAST fJetPtCut
+
+    varContainer **fVectParticle; //! container for Particles
+    varContainer **fVectJet;      //! container for Jets
+
+    Double_t *fPtArray;  //! momentum array
+    Int_t    *fIdxArray;  //! index array of sorted pts
+
+
+    ClassDef(AliCdfJetFinder,2)
+   };
 #endif
+