//
//____________________________________________________________--------------
- return PtScal(*px,3); // 3==> Eta in the PtScal function
+ return PtScal(*px,2); // 2==> Eta in the PtScal function
}
Double_t AliGenGSIlib::YEtaPHOS(const Double_t *py, const Double_t */*dummy*/)
//
//____________________________________________________________--------------
- return PtScal(*px,5); // 5==> Etaprime in the PtScal function
+ return PtScal(*px,4); // 4==> Etaprime in the PtScal function
}
Double_t AliGenGSIlib::YEtaprimePHOS(const Double_t *py, const Double_t */*dummy*/)
//
//____________________________________________________________--------------
- return PtScal(*px,4); // 4==> Omega in the PtScal function
+ return PtScal(*px,3); // 3==> Omega in the PtScal function
}
Double_t AliGenGSIlib::YOmega(const Double_t *py, const Double_t */*dummy*/)
//
//____________________________________________________________--------------
- return PtScal(*px,11); // 11==> Rho in the PtScal function
+ return PtScal(*px,10); // 10==> Rho in the PtScal function
}
Double_t AliGenGSIlib::YRho(const Double_t *py, const Double_t */*dummy*/)
//
//____________________________________________________________--------------
- return PtScal(*px,2); // 2==> Kaon in the PtScal function
+ return PtScal(*px,1); // 1==> Kaon in the PtScal function
}
Double_t AliGenGSIlib::YKaonPHOS(const Double_t *py, const Double_t */*dummy*/)
//
//____________________________________________________________--------------
- return PtScal(*px,6); // 6==> Phi in the PtScal function
+ return PtScal(*px,5); // 5==> Phi in the PtScal function
}
Double_t AliGenGSIlib::YPhiPHOS(const Double_t *py, const Double_t */*dummy*/)
//
//____________________________________________________________--------------
- return PtScal(*px,7); // 7==> Baryon in the PtScal function
+ return PtScal(*px,6); // 6==> Baryon in the PtScal function
}
Double_t AliGenGSIlib::YBaryons(const Double_t *py, const Double_t */*dummy*/)
// The present adds the Rhos
-// MASS 1=>PI, 2=>K, 3=>ETA, 4=>OMEGA, 5=>ETA', 6=>PHI
-// 7=>BARYON-BARYONBAR, 11==>RHO
+// MASS 0=>PI, 1=>K, 2=>ETA, 3=>OMEGA, 4=>ETA', 5=>PHI
+// 6=>BARYON-BARYONBAR, 10==>RHO
const Double_t khm[11] = {0.1396, 0.494, 0.547, 0.782, 0.957, 1.02,
0.938, 0. , 0., 0., 0.769};
const Double_t kfmax[11]={1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1.};
- np--;
Double_t f5=TMath::Power(((sqrt(100.018215)+2.)/(sqrt(100.+khm[np]*khm[np])+2.0)),12.3);
Double_t kfmax2=f5/kfmax[np];
// PIONS