fFlavorSelect = 4;
break;
case kPyBeauty:
+ case kPyBeautyJets:
case kPyBeautyPbPbMNR:
case kPyBeautypPbMNR:
case kPyBeautyppMNR:
case kPyJpsi:
fParentSelect[0] = 443;
break;
+ case kPyMbAtlasTuneMC09:
case kPyMbDefault:
case kPyMb:
case kPyMbWithDirectPhoton:
}
fNpartons = fPythia->GetNumberOfParticles();
} else {
- printf("Loading Event %d\n",AliRunLoader::GetRunLoader()->GetEventNumber());
- fRL->GetEvent(AliRunLoader::GetRunLoader()->GetEventNumber());
+ printf("Loading Event %d\n",AliRunLoader::Instance()->GetEventNumber());
+ fRL->GetEvent(AliRunLoader::Instance()->GetEventNumber());
fPythia->SetNumberOfParticles(0);
fPythia->LoadEvent(fRL->Stack(), 0 , 1);
fPythia->EditEventList(21);
if(mi<0) continue;
mother = (TParticle*)fParticles.At(mi);
mpdg=TMath::Abs(mother->GetPdgCode());
- mpdgUpperFamily=(mpdg>1000 ? mpdg-1000 : mpdg-100); // keep e from c from b
+ mpdgUpperFamily=(mpdg>1000 ? mpdg+1000 : mpdg+100); // keep e from c from b
if ( ParentSelected(mpdg) ||
(fFlavorSelect==5 && ParentSelected(mpdgUpperFamily))) {
if (KinematicSelection(partCheck,1)) {
if (fStack) {
particle = fStack->Particle(nt);
} else {
- particle = gAlice->Stack()->Particle(nt);
+ particle = AliRunLoader::Instance()->Stack()->Particle(nt);
}
// particle->SetFirstDaughter(fPythia->GetK(2, i));
// particle->SetLastDaughter(fPythia->GetK(3, i));