void Three_Particle_Correlation::
Evaluate(const Vec4D & mom1,const Vec4D & mom2,const Vec4D & mom3,double weight, double ncount)
{
double eta1 = mom1.Eta();
double eta2 = mom2.Eta();
double eta3 = mom3.Eta();
if (ATOOLS::dabs(eta1)>2. || ATOOLS::dabs(eta2)>2. || ATOOLS::dabs(eta3)>2.) return;
double phi1 = mom1.Phi();
double phi2 = mom2.Phi();
double phi3 = mom3.Phi();
double delta12 = ATOOLS::dabs(phi1-phi2); delta12 = ATOOLS::Min(delta12,2.*M_PI-delta12);
double delta23 = ATOOLS::dabs(phi2-phi3); delta23 = ATOOLS::Min(delta23,2.*M_PI-delta23);
double delta31 = ATOOLS::dabs(phi3-phi1); delta31 = ATOOLS::Min(delta31,2.*M_PI-delta31);
if (ATOOLS::dabs(eta1-eta2)>2.) p_histo->Insert(delta12,weight,ncount);
if (ATOOLS::dabs(eta2-eta3)>2.) p_histo->Insert(delta23,weight,ncount);
if (ATOOLS::dabs(eta3-eta1)>2.) p_histo->Insert(delta31,weight,ncount);
if (mom1.PPerp()>10. || mom2.PPerp()>10. || mom3.PPerp()>10.) {
if (ATOOLS::dabs(eta1-eta2)>2.) p_histo10->Insert(delta12,weight,ncount);
if (ATOOLS::dabs(eta2-eta3)>2.) p_histo10->Insert(delta23,weight,ncount);
if (ATOOLS::dabs(eta3-eta1)>2.) p_histo10->Insert(delta31,weight,ncount);
}
if (mom1.PPerp()>20. || mom2.PPerp()>20. || mom3.PPerp()>20.) {
if (ATOOLS::dabs(eta1-eta2)>2.) p_histo20->Insert(delta12,weight,ncount);
if (ATOOLS::dabs(eta2-eta3)>2.) p_histo20->Insert(delta23,weight,ncount);
if (ATOOLS::dabs(eta3-eta1)>2.) p_histo20->Insert(delta31,weight,ncount);
}
if (mom1.PPerp()>50. || mom2.PPerp()>50. || mom3.PPerp()>50.) {
if (ATOOLS::dabs(eta1-eta2)>2.) p_histo50->Insert(delta12,weight,ncount);
if (ATOOLS::dabs(eta2-eta3)>2.) p_histo50->Insert(delta23,weight,ncount);
if (ATOOLS::dabs(eta3-eta1)>2.) p_histo50->Insert(delta31,weight,ncount);
}
if (mom1.PPerp()>100. || mom2.PPerp()>100. || mom3.PPerp()>100.) {
if (ATOOLS::dabs(eta1-eta2)>2.) p_histo100->Insert(delta12,weight,ncount);
if (ATOOLS::dabs(eta2-eta3)>2.) p_histo100->Insert(delta23,weight,ncount);
if (ATOOLS::dabs(eta3-eta1)>2.) p_histo100->Insert(delta31,weight,ncount);
}
}
void EtaTracksAsym::Evaluate(const Vec4D & mom,double weight,double ncount)
{
double eta = 0.;
eta=mom.Eta();
if (eta<0.) {
p_histo->Insert(-eta,-weight,ncount);
}
else {
p_histo->Insert(eta,weight,ncount);
}
}
static bool OrderV(const Vec4D &a,const Vec4D &b)
{ return dabs(a.Eta())>dabs(b.Eta()); }
double EV_SC_Eta::Calc(const Vec4D &mom)
{
return mom.Eta();
}
void Three_Particle_DEta::Evaluate(const Vec4D & mom1,const Vec4D & mom2,const Vec4D & mom3,double weight, double ncount)
{
Vec4D mother = mom1+mom2;
double deta = abs((mother.Eta()-mom3.Eta()));
p_histo->Insert(deta,weight,ncount);
}