fZDCCutRefDelta(-2.1),
fZDCCutSigmaSum(3.25),
fZDCCutSigmaDelta(2.25),
+ fZDCCutRefSumCorr(-65.5),
+ fZDCCutRefDeltaCorr(-2.1),
+ fZDCCutSigmaSumCorr(6.0),
+ fZDCCutSigmaDeltaCorr(1.2),
fDoFMD(kTRUE),
fFMDLowCut(0.2),
fFMDHitCut(0.5),
continue;
if (track->GetTPCNcls() < 30) continue;
- if (track->GetTPCchi2()/track->GetTPCNcls() > 0.3) continue;
if (TMath::Abs(track->Eta()) > 0.005) continue;
if (track->Pt() < 4) continue;
if (track->GetKinkIndex(0) > 0) continue;
UInt_t status = track->GetStatus();
- if ((status&AliESDtrack::kITSrefit)==1) continue; // explicitly ask for tracks without ITS refit
- if ((status&AliESDtrack::kTPCrefit)==0) continue;
+ if ((status&AliESDtrack::kITSrefit)==AliESDtrack::kITSrefit) continue; // explicitly ask for tracks without ITS refit
+ if ((status&AliESDtrack::kTPCrefit)!=AliESDtrack::kTPCrefit) continue;
if (track->GetTPCsignal() > 10) continue; // explicitly ask for tracks without dE/dx
- if (TMath::Abs(track->GetZ()) < 50) continue;
-
trackCounter++;
}
if (trackCounter > 15)
zdcAccept = (znaFired | zncFired);
}
else {
-
- // Cuts for the new corrected TDC values
- const Float_t refSumCorr = -65.5;
- const Float_t sigmaSumCorr = 6.0;
- const Float_t sigmaDeltaCorr = 1.2;
-
for(Int_t i = 0; i < 4; ++i) {
if (esdZDC->GetZDCTDCData(10,i) != 0) {
Float_t tdcC = 0.025*(esdZDC->GetZDCTDCData(10,i)-esdZDC->GetZDCTDCData(14,i));
fHistTimeCorrZDC->Fill(tdcCcorr-tdcAcorr,tdcCcorr+tdcAcorr);
}
if (esdZDC->TestBit(AliESDZDC::kCorrectedTDCFilled)) {
- if (((tdcCcorr-tdcAcorr-fZDCCutRefDelta)*(tdcCcorr-tdcAcorr-fZDCCutRefDelta)/(sigmaDeltaCorr*sigmaDeltaCorr) +
- (tdcCcorr+tdcAcorr-refSumCorr)*(tdcCcorr+tdcAcorr-refSumCorr)/(sigmaSumCorr*sigmaSumCorr))< 1.0)
+ if (((tdcCcorr-tdcAcorr-fZDCCutRefDeltaCorr)*(tdcCcorr-tdcAcorr-fZDCCutRefDeltaCorr)/(fZDCCutSigmaDeltaCorr*fZDCCutSigmaDeltaCorr) +
+ (tdcCcorr+tdcAcorr-fZDCCutRefSumCorr)*(tdcCcorr+tdcAcorr-fZDCCutRefSumCorr)/(fZDCCutSigmaSumCorr*fZDCCutSigmaSumCorr))< 1.0)
zdcAccept = kTRUE;
}
else {