//____________________________________________________________
AliGenParam::AliGenParam(Int_t npart, Int_t param,
- Double_t (*PtPara) (Double_t*, Double_t*),
- Double_t (*YPara ) (Double_t* ,Double_t*),
+ Double_t (*PtPara) (const Double_t*, const Double_t*),
+ Double_t (*YPara ) (const Double_t* ,const Double_t*),
Int_t (*IpPara) (TRandom *))
:AliGenMC(npart),
"Division by 0: Please check you rapidity range !");
}
- pl=xmt*ty/sqrt(1.-ty*ty);
+ pl=xmt*ty/sqrt((1.-ty)*(1.+ty));
theta=TMath::ATan2(pt,pl);
// Cut on theta
if(theta<fThetaMin || theta>fThetaMax) continue;
//
// children
- if (ChildSelected(TMath::Abs(kf)) || fForceDecay == kAll && trackIt[i])
+ if ((ChildSelected(TMath::Abs(kf)) || fForceDecay == kAll) && trackIt[i])
{
if (fCutOnChild) {
pc[0]=iparticle->Px();