Bug in weight setting is fixed (Boris Polishchuk)

diff --git a/PWGGA/PHOSTasks/PHOS_PbPb/AliAnalysisTaskPi0FlowMC.cxx b/PWGGA/PHOSTasks/PHOS_PbPb/AliAnalysisTaskPi0FlowMC.cxx
index 8a0119fd64553db107c1a9ac840ebd63d89bf11d..704aba07da63a00ec1267c8bfed67ca0499e6dd0 100644 (file)
--- a/PWGGA/PHOSTasks/PHOS_PbPb/AliAnalysisTaskPi0FlowMC.cxx
+++ b/PWGGA/PHOSTasks/PHOS_PbPb/AliAnalysisTaskPi0FlowMC.cxx
@@ -403,9 +403,7 @@ void AliAnalysisTaskPi0FlowMC::ConsiderPi0s()
       TLorentzVector pv12 = *(ph1->GetMomV2()) + *(ph2->GetMomV2());
       
       const Double_t w2 = ph2->GetWeight();
-      Double_t w = w1*w2;
-      if( FindCommonParent(ph1->GetPrimary(), ph1->GetPrimary()) )
-       w = w1*w2;
+      Double_t w = TMath::Sqrt(w1*w2);
       
       FillHistogram("hPHOSphi",fCentralityV0M,p12.Pt(),p12.Phi(), w) ;
       Double_t dphiA=p12.Phi()-fRPV0A ;
@@ -669,7 +667,7 @@ void AliAnalysisTaskPi0FlowMC::ConsiderPi0sMix()
        TLorentzVector pv12 = *(ph1->GetMomV2()) + *(ph2->GetMomV2());
        
        const Double_t w2 = ph2->GetWeight();
-       Double_t w = w1*w2;
+       Double_t w = TMath::Sqrt(w1*w2);
 
        Double_t dphiA=p12.Phi()-fRPV0A ;
        while(dphiA<0)dphiA+=TMath::Pi() ;