message: Add further track cut for PbPb2018 in MCAODrec option for MCClosure
diff: @@ -3162,6 +3162,11 @@ TObjArray* AliAnalysisTaskBFPsi::GetAcceptedTracks(AliVEvent *event, Double_t gC
  return tracksAccepted;
  }
  + if (fUseRaaGeoCut){
  + fESDtrackCuts = new AliESDtrackCuts();
  + fESDtrackCuts->SetCutGeoNcrNcl(fDeadZoneWidth, fCutGeoNcrNclLength, fCutGeoNcrNclGeom1Pt, fCutGeoNcrNclFractionNcr, fCutGeoNcrNclFractionNcl);
  + }
  +
  for (Int_t iTracks = 0; iTracks < event->GetNumberOfTracks(); iTracks++) {
  AliAODTrack* aodTrack = dynamic_cast<AliAODTrack *>(event->GetTrack(iTracks));
  if (!aodTrack) {
  @@ -3186,6 +3191,29 @@ TObjArray* AliAnalysisTaskBFPsi::GetAcceptedTracks(AliVEvent *event, Double_t gC
  }
  }
  + if(fUseRaaGeoCut){
  + if (!fESDtrackCuts->IsSelected(aodTrack))
  + continue;
  + }
  +
  + if (fUseTOFBCPileUpCut) {
  + if (!aodTrack->GetTOFBunchCrossing()==0)
  + continue;
  + }
  +
  + if (fUseTPCInOutRowsCut) {
  + const TBits& bmap = aodTrack->GetTPCClusterMap();
  + // require at least 20 out of 25 and 3 out of 5 innermost rows
  + int nset25 = 0, nset5 = 0;
  + for (int i=0;i<25; i++) {
  + if (!bmap.TestBitNumber(i)) continue;
  + nset25++;
  + if (i<5) nset5++;
  + }
  + if((nset5<fInRows) || (nset25<fOutRows))
  + continue;
  + }
  +
  vCharge = aodTrack->Charge();
  vEta = aodTrack->Eta();
  @@ -3594,6 +3622,7 @@ TObjArray* AliAnalysisTaskBFPsi::GetAcceptedTracks(AliVEvent *event, Double_t gC
  }
  }//track loop
  + if (fUseRaaGeoCut) delete fESDtrackCuts;
  }//MCAODrec
  //==============================================================================================================
  