Filling the QA histograms also for the MC case

diff --git a/PWG2/EBYE/AliAnalysisTaskBF.cxx b/PWG2/EBYE/AliAnalysisTaskBF.cxx
index 2e43652ee9b7a3ece4fc3c3a219e217c7d054958..9b86aac6a1ab74fdd936b3035ed765562927cfaf 100755 (executable)
--- a/PWG2/EBYE/AliAnalysisTaskBF.cxx
+++ b/PWG2/EBYE/AliAnalysisTaskBF.cxx
@@ -754,6 +754,8 @@ void AliAnalysisTaskBF::UserExec(Option_t *) {
        gImpactParameter = headerH->ImpactParameter();\r
        fCentrality = gImpactParameter;\r
       }\r
+      fCentrality = gImpactParameter;\r
+\r
       // take only events inside centrality class (DIDN'T CHANGE THIS UP TO NOW)\r
       if((fImpactParameterMin > gImpactParameter) || (fImpactParameterMax < gImpactParameter))\r
        return;\r
@@ -844,6 +846,10 @@ void AliAnalysisTaskBF::UserExec(Option_t *) {
            track->PxPyPz(v_p);\r
            //Printf("phi (before): %lf",v_phi);\r
 \r
+           fHistPt->Fill(v_pt);\r
+           fHistEta->Fill(v_eta);\r
+           fHistPhi->Fill(v_phi);\r
+\r
            //Flow after burner\r
            if(fUseFlowAfterBurner) {\r
              Double_t precisionPhi = 0.001;\r
@@ -852,6 +858,7 @@ void AliAnalysisTaskBF::UserExec(Option_t *) {
              Double_t phi0 = v_phi;\r
              Double_t gV2 = fDifferentialV2->Eval(v_pt);\r
              fHistPhiBefore->Fill(v_phi);\r
+             gReactionPlane = 0.0;\r
 \r
              for (Int_t j = 0; j < maxNumberOfIterations; j++) {\r
                Double_t phiprev = v_phi;\r

diff --git a/PWG2/EBYE/AliBalance.cxx b/PWG2/EBYE/AliBalance.cxx
index 1e6b1679b5293c125070d0e0472e6dfca8710258..b9400068572e9d7db297f6f541ffd77a18b628b3 100644 (file)
--- a/PWG2/EBYE/AliBalance.cxx
+++ b/PWG2/EBYE/AliBalance.cxx
@@ -683,9 +683,9 @@ void AliBalance::PrintResults(Int_t iAnalysisType, TH1D *gHistBalance) {
   cout<<"=================================================="<<endl;
   for(Int_t i = 1; i <= fNumberOfBins[iAnalysisType]; i++) { 
     x[i-1] = fP2Start[iAnalysisType] + fP2Step[iAnalysisType]*i + fP2Step[iAnalysisType]/2;
-    cout<<"B: "<<gHistBalance->GetBinContent(i)<<"\t Error: "<<gHistBalance->GetBinError(i)<<"\t bin: "<<gHistBalance->GetBinCenter(i)<<endl;
+    //cout<<"B: "<<gHistBalance->GetBinContent(i)<<"\t Error: "<<gHistBalance->GetBinError(i)<<"\t bin: "<<gHistBalance->GetBinCenter(i)<<endl;
   } 
-  cout<<"=================================================="<<endl;
+  //cout<<"=================================================="<<endl;
   for(Int_t i = 2; i <= fNumberOfBins[iAnalysisType]; i++) {
     gSumXi += gHistBalance->GetBinCenter(i);
     gSumBi += gHistBalance->GetBinContent(i);
@@ -705,7 +705,7 @@ void AliBalance::PrintResults(Int_t iAnalysisType, TH1D *gHistBalance) {
   
   Double_t delta = gSumBiXi / gSumBi;
   Double_t deltaError = (gSumBiXi / gSumBi) * TMath::Sqrt(TMath::Power((TMath::Sqrt(gSumXi2DeltaBi2)/gSumBiXi),2) + TMath::Power((gSumDeltaBi2/gSumBi),2) );
-  
+  cout<<"Analysis type: "<<gBFAnalysisType[iAnalysisType].Data()<<endl;
   cout<<"Width: "<<delta<<"\t Error: "<<deltaError<<endl;
   cout<<"New error: "<<deltaErrorNew<<endl;
   cout<<"Integral: "<<integral<<"\t Error: "<<integralError<<endl;