// Set b-quark mass
SetPMAS(5,1,4.75);
break;
+ case kPyBeautyJets:
case kPyBeautyppMNRwmi:
// Tuning of Pythia parameters aimed to get a resonable agreement
// between with the NLO calculation by Mangano, Nason, Ridolfi for the
//
// Isotropic decay ????
Double_t cost = 2. * gRandom->Rndm() - 1.;
- Double_t sint = TMath::Sqrt(1. - cost * cost);
+ Double_t sint = TMath::Sqrt((1.-cost)*(1.+cost));
Double_t phis = 2. * TMath::Pi() * gRandom->Rndm();
Double_t pz1 = pst * cost;