Streamlining. Jet finder should be used with AliAnalysisTaskJets.
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 23 Jul 2009 10:27:29 +0000 (10:27 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 23 Jul 2009 10:27:29 +0000 (10:27 +0000)
JETAN/AliCdfJetFinder.cxx
JETAN/AliCdfJetFinder.h

index 4a88372..9007581 100644 (file)
@@ -64,6 +64,7 @@ AliCdfJetFinder::AliCdfJetFinder():
     fFromAod(0),
     fAODwrite(0),
     fAODtracksWrite(0),
+    fAnalyseJets(0),
     fRefArr (NULL),
     fNJets(-9999),
     fNPart(-9999),
@@ -366,7 +367,7 @@ ComputeConesWeight();
 
 if (fAODwrite) { cout << "Writing AOD" << endl ; WriteJets(); }
 
-AnalizeJets();
+if (fAnalyseJets) AnalizeJets();
 
 Clean();
 
@@ -378,7 +379,7 @@ void AliCdfJetFinder::InitData()
 // initialisation of variables and data members
 
   TClonesArray * vectArray = fReader->GetMomentumArray() ;
-    if ( vectArray == 0 ) { cout << "Could not get the momentum array" << endl; return; }
+  if ( vectArray == 0 ) { cout << "Could not get the momentum array" << endl; return; }
 
   fNPart = vectArray->GetEntries()  ; // n particles in this event;
 
@@ -583,17 +584,17 @@ for(  Int_t jetnr = 0 ; jetnr < fNJets ; jetnr++ )
   en = TMath::Sqrt ( px * px + py * py + pz * pz );
 
   AliAODJet jet (px, py, pz, en);
-  AddJet(jet);
+
 
   if (fDebug) jet.Print("");
 
   if (fFromAod && fAODtracksWrite)
-    {
-    for (  Int_t jetTrack = 0; jetTrack < fNPart; jetTrack++ )
+  {
+      for (  Int_t jetTrack = 0; jetTrack < fNPart; jetTrack++ )
       { if ( fVectParticle[jetTrack]->njet == jetnr ) { jet.AddTrack(fRefArr->At(jetTrack)) ; } }
-    }
-    // tracks REFs written in AOD
-
+  }
+  // tracks REFs written in AOD
+  AddJet(jet);
   }
 //jets vector parsed and written to AOD
 }
@@ -604,6 +605,8 @@ void AliCdfJetFinder::AnalizeJets()
 {
 // analyzing of jets and filling of histograms
 
+    const Double_t kPI = TMath::Pi();
+    
   //persistent pointer to histo20
   TH1F *hR = (TH1F*)fHistos->FindObject("histo20");
 
index 09896dd..add15c7 100644 (file)
 
 using namespace std ;
 
-Double_t const  kPI        = TMath::Pi();
-// Double_t const  kTWOPI     = 2.*kPI;
 
 // 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 
+  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;
@@ -42,9 +40,11 @@ class AliCdfJetFinder : public AliJetFinder
     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 );
@@ -55,6 +55,8 @@ class AliCdfJetFinder : public AliJetFinder
     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
@@ -69,9 +71,10 @@ class AliCdfJetFinder : public AliJetFinder
     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