first step for error estimation QC{2}
authorsnelling <snelling@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 30 Jun 2009 15:32:24 +0000 (15:32 +0000)
committersnelling <snelling@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 30 Jun 2009 15:32:24 +0000 (15:32 +0000)
PWG2/FLOW/AliFlowCommon/AliFlowAnalysisWithQCumulants.cxx
PWG2/FLOW/macros/compareFlowResults.C

index 3f59f863d945b2e43495c7ebc13d4069c5e863eb..0e99b372ffa5061bfa0ffe0307a40da6c2e4503e 100644 (file)
@@ -5190,6 +5190,14 @@ void AliFlowAnalysisWithQCumulants::CalculateFinalResultsForNoNameIntegratedFlow
  //sixthOrderQCumulant = sixthOrderQCumulant - sixCorrection;
  //eightOrderQCumulant = eightOrderQCumulant - eightCorrection;
  
  //sixthOrderQCumulant = sixthOrderQCumulant - sixCorrection;
  //eightOrderQCumulant = eightOrderQCumulant - eightCorrection;
  
+ // errors: (to be improved i.e. reimplemented)
+ Double_t secondOrderError = 0;
+ if(two > 0. && !(useWeights))
+ { 
+  Double_t nEvtsNoName = (fCommonHists2nd->GetHistMultRP())->GetEntries(); 
+  if(nEvtsNoName>0) secondOrderError = (fQCorrelations->GetBinError(1)/(2.*pow(two,0.5)))/(pow(nEvtsNoName,0.5));
+ }
  if(useWeights) sixthOrderQCumulant = 0.; // to be removed (once 6th order with weights is calculated)
  if(useWeights) eightOrderQCumulant = 0.; // to be removed (once 8th order with weights is calculated)
  
  if(useWeights) sixthOrderQCumulant = 0.; // to be removed (once 6th order with weights is calculated)
  if(useWeights) eightOrderQCumulant = 0.; // to be removed (once 8th order with weights is calculated)
  
@@ -5203,7 +5211,7 @@ void AliFlowAnalysisWithQCumulants::CalculateFinalResultsForNoNameIntegratedFlow
   if(!(useWeights)) 
   { 
    fIntFlowResultsQC->SetBinContent(1,dVn2);
   if(!(useWeights)) 
   { 
    fIntFlowResultsQC->SetBinContent(1,dVn2);
-   fIntFlowResultsQC->SetBinError(1,0.); // to be improved
+   fIntFlowResultsQC->SetBinError(1,secondOrderError); // to be improved
   }
   if(useWeights)
   { 
   }
   if(useWeights)
   { 
@@ -5212,7 +5220,7 @@ void AliFlowAnalysisWithQCumulants::CalculateFinalResultsForNoNameIntegratedFlow
   }
   
   // fill common histogram:
   }
   
   // fill common histogram:
-  fCommonHistsResults2nd->FillIntegratedFlow(dVn2, 0.); // to be improved 
+  fCommonHistsResults2nd->FillIntegratedFlow(dVn2,secondOrderError); // to be improved 
   
  } 
  if(fourthOrderQCumulant<0.)
   
  } 
  if(fourthOrderQCumulant<0.)
index 72a9c4ec4d6f0cdfd2cf1b98181d2b49bab607e7..d1dc8e86181c9ab8c23ae8ecc89d020ef37f446e 100644 (file)
@@ -613,7 +613,7 @@ void compareFlowResults(TString type="ESD",Int_t mode=mLocal)
       qcCommonHistRes2 = dynamic_cast<AliFlowCommonHistResults*> (pListQC->FindObject("AliFlowCommonHistResults2ndOrderQC"));
       if(qcCommonHistRes2)  {
        flowValue[binQC2-1] = (qcCommonHistRes2->GetHistIntFlow())->GetBinContent(1);
       qcCommonHistRes2 = dynamic_cast<AliFlowCommonHistResults*> (pListQC->FindObject("AliFlowCommonHistResults2ndOrderQC"));
       if(qcCommonHistRes2)  {
        flowValue[binQC2-1] = (qcCommonHistRes2->GetHistIntFlow())->GetBinContent(1);
-       //flowError[binQC2-1] = (qcCommonHistRes2->GetHistIntFlow())->GetBinError(1);
+       flowError[binQC2-1] = (qcCommonHistRes2->GetHistIntFlow())->GetBinError(1);
        flowValueRP[binQC2RP-1] = (qcCommonHistRes2->GetHistIntFlowRP())->GetBinContent(1);
        //flowErrorRP[binQC2RP-1] = (qcCommonHistRes2->GetHistIntFlowRP())->GetBinError(1);
        flowValuePOI[binQC2POI-1] = (qcCommonHistRes2->GetHistIntFlowPOI())->GetBinContent(1);
        flowValueRP[binQC2RP-1] = (qcCommonHistRes2->GetHistIntFlowRP())->GetBinContent(1);
        //flowErrorRP[binQC2RP-1] = (qcCommonHistRes2->GetHistIntFlowRP())->GetBinError(1);
        flowValuePOI[binQC2POI-1] = (qcCommonHistRes2->GetHistIntFlowPOI())->GetBinContent(1);