const Int_t kNLeadingPhiBins = 40;
Double_t leadingPhiBins[kNLeadingPhiBins+1];
for (Int_t i=0; i<=kNLeadingPhiBins; i++)
- leadingPhiBins[i] = -1.5 * TMath::Pi() + 1.0 / 40 * i * TMath::TwoPi();
+ leadingPhiBins[i] = -0.5 * TMath::Pi() + 1.0 / 40 * i * TMath::TwoPi();
// multiplicity
const Int_t kNMultiplicityBins = 15;
vars[2] = leading->Pt();
vars[3] = multiplicity;
vars[4] = leading->Phi() - particle->Phi();
- if (vars[4]>TMath::Pi()/2.)vars[4]-=TMath::TwoPi();
- if (vars[4]< - 3.*TMath::Pi()/2.)vars[4]+=TMath::TwoPi();
+ if (vars[4] > 1.5 * TMath::Pi())
+ vars[4] -= TMath::TwoPi();
+ if (vars[4] < -0.5 * TMath::Pi())
+ vars[4] += TMath::TwoPi();
fNumberDensitypT->GetTrackHist(region)->Fill(vars, step);
fSumpT->GetTrackHist(region)->Fill(vars, step, particle->Pt());