AliAODTrack *aodtrack = fAOD->GetTrack(iTrack);
if(!aodtrack)
continue;
- if( !aodtrack->TestFilterBit(fFilterMask) )
+ if((!aodtrack->TestFilterBit(fFilterMask)) ||
+ ((!fCFManagerPos->CheckParticleCuts(kStepReconstructed,aodtrack)) && (aodtrack->Charge()>0.)) ||
+ ((!fCFManagerNeg->CheckParticleCuts(kStepReconstructed,aodtrack)) && (aodtrack->Charge()<0)) )
continue;
else {
track = static_cast<AliVTrack*>(aodtrack);
//Container with primaries
if(particle->IsPhysicalPrimary()) {
- if(particle->Charge()>0. && fCFManagerPos->CheckParticleCuts(kStepReconstructedMC,particle)) fCFManagerPos->GetParticleContainer()->Fill(containerInputRecMC,kStepReconstructedMC);
- if(particle->Charge()<0. && fCFManagerNeg->CheckParticleCuts(kStepReconstructedMC,particle)) fCFManagerNeg->GetParticleContainer()->Fill(containerInputRecMC,kStepReconstructedMC);
+ if(particle->Charge()>0.) fCFManagerPos->GetParticleContainer()->Fill(containerInputRecMC,kStepReconstructedMC);
+ if(particle->Charge()<0.) fCFManagerNeg->GetParticleContainer()->Fill(containerInputRecMC,kStepReconstructedMC);
//Fill pT resolution plots for primaries
//fPtRelUncertainty1PtPrim->Fill(containerInputRec[0],containerInputRec[0]*TMath::Sqrt(track->GetSigma1Pt2())); //This has not been implemented in AOD analysis, since they are also produced by the AddTaskPWG4HighPtTrackQA.C macro
}
// mcAccCuts->SetMinNHitITS(mintrackrefsITS);
//mcAccCuts->SetMinNHitTPC(mintrackrefsTPC);
- TObjArray* recList = new TObjArray(0);
+ TObjArray* recMCList = new TObjArray(0);
TObjArray* secList = new TObjArray(0) ;
printf("CREATE MC KINE CUTS\n");
- TObjArray* recMCList = new TObjArray(0);
+ TObjArray* recList = new TObjArray(0);
TObjArray* mcList = new TObjArray(0) ;
mcList->AddLast(mcKineCuts);
- recMCList->AddLast(mcKineCuts);
+ recList->AddLast(mcKineCuts);
//mcList->AddLast(mcAccCuts);
//CREATE THE INTERFACE TO CORRECTION FRAMEWORK USED IN THE TASK