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);
}
}
double PHASIC::ComputePhi
(Vec4D pijt,Vec4D pkt,Vec4D pi,const int mode)
{
LN_Pair ln(GetLN((mode&1)?-pijt:pijt,(mode&2)?-pkt:pkt,0));
Vec4D ktt(0.0,cross(Vec3D(ln.m_l),Vec3D(ln.m_n)));
Poincare cms(ln.SaneCMS());
cms.Boost(ln.m_l);
cms.Boost(pi);
Poincare zax(ln.m_l,Vec4D::ZVEC);
if (mode!=3 && ktt.PSpat2()>rpa->gen.SqrtAccu()) {
zax.Rotate(ktt);
}
else {
msg_Debugging()<<"Set fixed n_perp\n";
if (IsZero(ln.m_l.PPerp2())) ln.m_l[1]=ln.m_l[2]=0.0;
zax=Poincare(ln.m_l,Vec4D::ZVEC);
ktt=Vec4D(0.0,1.0,1.0,0.0);
}
zax.Rotate(pi);
Poincare xax(ktt,Vec4D::XVEC);
xax.Rotate(pi);
return pi.Phi();
}
void One_Particle_Phi::EvaluateNLOcontrib(const Vec4D & mom,double weight, double ncount)
{
p_histo->Insert(mom.Phi(),weight,ncount);
}