protect against invoking MC Handler for AOD

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

diff --git a/JETAN/AliCdfJetFinder.h b/JETAN/AliCdfJetFinder.h
index c3dbdab9191063e92db7b0796bcbde6eec5c1942..add15c724d07401362b1ad2dcd17bf022b62e361 100644 (file)
--- a/JETAN/AliCdfJetFinder.h
+++ b/JETAN/AliCdfJetFinder.h
@@ -1,29 +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
+ *
+*/
 
-#include "AliJetFinder.h"
-#include <math.h>
+//  Definition of constants, structures and functions for jet finder
 
-//
-//  Definition of constants, structures and functions
-//
+#include "AliJetFinder.h"
 
 using namespace std ;
 
-const Double_t pi = TMath::Pi();
 
-struct varContainer // container for Particle and Jets
-  { // variables of container struct
-  Double_t  pt; Double_t eta; Double_t phi;
-  Int_t njet; // if jets are stored in varContainer njet is multiplicity of jet
-              // if particles are stored , njet is index number of jet
+// 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
@@ -33,63 +32,49 @@ 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           WriteJets() ;
     void           AnalizeJets();
     void           Clean();
-
+    
     virtual void   FinishRun();
 
-    Double_t DeltaPhiNorm (Double_t dphi)
-      {
-      if ( dphi < - pi ) { dphi = -dphi - 2.0 * pi ; }
-      if ( dphi >   pi ) { dphi = -dphi + 2.0 * pi ; }
-      return dphi;
-      }
-
-    inline long double Distance (long double x, long double y)
-      { return TMath::Sqrt ( pow(x,2) + pow(y,2) ); }
-
-
-
+    void           SetAnalyseJets(Bool_t flag = kTRUE) {fAnalyseJets = flag;}
+    
   protected:
     AliCdfJetFinder ( const AliCdfJetFinder& jf );
     AliCdfJetFinder& operator = ( const AliCdfJetFinder& jf );
 
     TList         *fHistos;    // List of histograms
+    Bool_t        fDebug;      //  enables debugging
 
-    Bool_t        fDebug;   //  enables debugging
-
-    Bool_t fFromAod;
-    Bool_t fAODwrite;   // write jets to AOD
-    Bool_t fAODtracksWrite;  // write jet tracks to AOD
-    TRefArray *fRefArr ; // pointer to references array of tracks from AOD
-
+    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;
+    Double_t      fRadius ;   // radius of jet 
 
-    Int_t fMinJetParticles;
-    Double_t fJetPtCut;
+    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
-
-
-
-
+    Int_t    *fIdxArray;  // index array of sorted pts
 
 
-    ClassDef ( AliCdfJetFinder, 1 )
-  };//
+    ClassDef(AliCdfJetFinder,1)
+   };
 #endif