From: lmilano <lmilano@cern.ch>
Date: Thu, 4 Sep 2014 09:01:58 +0000 (+0200)
Subject: track cuts updated
X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=commitdiff_plain;h=16e0a56f4eab48e8a9c6093c54ce69c342aefd29

track cuts updated
---

diff --git a/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraAODEventCuts.cxx b/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraAODEventCuts.cxx
index 23cb5b7e4f1..5eb16e1e9d2 100644
--- a/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraAODEventCuts.cxx
+++ b/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraAODEventCuts.cxx
@@ -379,20 +379,15 @@ Double_t AliSpectraAODEventCuts::CalculateQVector(){
   //V0 info    
   Double_t Qxa2 = 0, Qya2 = 0;
   Double_t Qxc2 = 0, Qyc2 = 0;
-  Double_t sumMc = 0, sumMa = 0;
   
   AliAODVZERO* aodV0 = fAOD->GetVZEROData();
   
   for (Int_t iv0 = 0; iv0 < 64; iv0++) {
     
     Float_t multv0 = aodV0->GetMultiplicity(iv0);
+  
     ((TH2F*)fOutput->FindObject("fV0M"))->Fill(iv0,multv0);
     
-    if (iv0 < 32)   
-      sumMc += aodV0->GetMultiplicity(iv0);   
-    else
-      sumMa += aodV0->GetMultiplicity(iv0);
-    
   }
 
   fPsiV0A = fAOD->GetEventplane()->CalculateVZEROEventPlane(fAOD,8,2,Qxa2,Qya2); // V0A
@@ -401,12 +396,12 @@ Double_t AliSpectraAODEventCuts::CalculateQVector(){
   ((TH2F*)fOutput->FindObject("fPsiACor"))->Fill((Float_t)fAOD->GetCentrality()->GetCentralityPercentile("V0M"), fPsiV0A);
   ((TH2F*)fOutput->FindObject("fPsiCCor"))->Fill((Float_t)fAOD->GetCentrality()->GetCentralityPercentile("V0M"), fPsiV0C);
   
-  fqV0A = TMath::Sqrt((Qxa2*Qxa2 + Qya2*Qya2)/sumMa);
-  fqV0C = TMath::Sqrt((Qxc2*Qxc2 + Qyc2*Qyc2)/sumMc);
-  fqV0Ax = Qxa2*TMath::Sqrt(1./sumMa);
-  fqV0Cx = Qxc2*TMath::Sqrt(1./sumMc);
-  fqV0Ay = Qya2*TMath::Sqrt(1./sumMa);
-  fqV0Cy = Qyc2*TMath::Sqrt(1./sumMc);
+  fqV0A = TMath::Sqrt((Qxa2*Qxa2 + Qya2*Qya2));
+  fqV0C = TMath::Sqrt((Qxc2*Qxc2 + Qyc2*Qyc2));
+  fqV0Ax = Qxa2;
+  fqV0Cx = Qxc2;
+  fqV0Ay = Qya2;
+  fqV0Cy = Qyc2;
   
   ((TH2F*)fOutput->FindObject("fQVecACor"))->Fill((Float_t)fAOD->GetCentrality()->GetCentralityPercentile("V0M"), fqV0A);
   ((TH2F*)fOutput->FindObject("fQVecCCor"))->Fill((Float_t)fAOD->GetCentrality()->GetCentralityPercentile("V0M"), fqV0C);