small fix
authorsnelling <snelling@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 30 Jan 2009 17:49:24 +0000 (17:49 +0000)
committersnelling <snelling@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 30 Jan 2009 17:49:24 +0000 (17:49 +0000)
PWG2/FLOW/macros/runFlowAnalysis.C

index 7e7accd..d5bdf5a 100644 (file)
@@ -171,6 +171,13 @@ int runFlowAnalysis(Int_t aRuns = 100, const char*
     AliFittingQDistribution* fqd = new AliFittingQDistribution();
     fqd->CreateOutputObjects();
   }
+
+  //SP = Scalar Product 
+  if(SP) {
+    AliFlowAnalysisWithScalarProduct* sp = new AliFlowAnalysisWithScalarProduct();
+    sp->Init();
+  }
+
   //LYZ1 = Lee-Yang Zeroes first run
   if(LYZ1) {
     AliFlowAnalysisWithLeeYangZeros* lyz1 = new AliFlowAnalysisWithLeeYangZeros();
@@ -223,12 +230,6 @@ int runFlowAnalysis(Int_t aRuns = 100, const char*
       }
     }
   }
-  
-  //SP = Scalar Product 
-  if(SP) {
-    AliFlowAnalysisWithScalarProduct* sp = new AliFlowAnalysisWithScalarProduct();
-    sp->Init();
-  }
   //------------------------------------------------------------------------
   
   //standard code
@@ -326,7 +327,8 @@ int runFlowAnalysis(Int_t aRuns = 100, const char*
              AliFlowEventSimple* fEvent = fEventMaker->FillTracks(kTree, cutsInt, cutsDiff);
              
              //pass the flow event to flow methods for analysis 
-             if (MCEP) //mcep->Make(fEvent,fEP);  //fix fEP
+             Double_t fEP = 0.; // temporary number need true value
+             if(MCEP) mcep->Make(fEvent,fEP);  //fix fEP
              if(QC) qc->Make(fEvent);
              if(GFC) gfc->Make(fEvent);
              if(FQD) fqd->Make(fEvent);
@@ -348,7 +350,7 @@ int runFlowAnalysis(Int_t aRuns = 100, const char*
 
   //--------------------------------------------------------------
   //calculating and storing the final results of flow analysis
-  if (MCEP) {mcep->Finish(); mcep->WriteHistograms("outputMCEPanalysis.root");}
+  if(MCEP) {mcep->Finish(); mcep->WriteHistograms("outputMCEPanalysis.root");}
   if(QC) {qc->Finish(); qc->WriteHistograms("outputQCanalysis.root");}
   if(GFC) {gfc->Finish(); gfc->WriteHistograms("outputGFCanalysis.root");}
   if(FQD) {fqd->Finish(); fqd->WriteHistograms("outputFQDanalysis.root");}