+
+ //Only fill the MC containers if MC information is available
+ if(eventHandler) {
+ Int_t label = TMath::Abs(track->GetLabel());
+ TParticle *particle = stack->Particle(label) ;
+ if(!particle) continue;
+
+ containerInputTPConlyMC[0] = particle->Pt();
+ containerInputTPConlyMC[1] = particle->Phi();
+ containerInputTPConlyMC[2] = particle->Eta();
+
+ //Container with primaries
+ if(stack->IsPhysicalPrimary(label)) {
+ if(particle->GetPDG()->Charge()>0.) {
+ fCFManagerPos->GetParticleContainer()->Fill(containerInputMC,kStepReconstructedTPCOnlyMC);
+ }
+ if(particle->GetPDG()->Charge()<0.) {
+ fCFManagerNeg->GetParticleContainer()->Fill(containerInputMC,kStepReconstructedTPCOnlyMC);
+ }
+ }
+ }
+