break;
case kHardMuons:
case kBSemiMuonic:
+ case kDSemiMuonic:
case kSemiMuonic:
case kDiMuon:
case kJpsiDiMuon:
Double_t gamma = 1./TMath::Sqrt((1.-beta)*(1.+beta));
Double_t gb = gamma * beta;
- printf("\n Boosting particles to lab frame %f %f %f", fDyBoost, beta, gamma);
+ // printf("\n Boosting particles to lab frame %f %f %f", fDyBoost, beta, gamma);
Int_t i;
Int_t np = fParticles.GetEntriesFast();
// Applies a boost in the y-direction in order to take into account the
// beam crossing angle
- Double_t thetaPr0, phiPr0, pyPr2, pzPr2;
+ Double_t thetaPr0, pyPr2, pzPr2;
TVector3 beta;
thetaPr0 = fXingAngleY / 2.;
- phiPr0 = 0;
// Momentum of the CMS system
pyPr2 = TMath::Sqrt(fEnergyCMS * fEnergyCMS/ 4 - 0.938 * 0.938) * TMath::Sin(thetaPr0);