// Checking if Dplus is in fiducial acceptance region
//
+ if(fMaxRapidityCand>-998.){
+ if(TMath::Abs(y) > fMaxRapidityCand) return kFALSE;
+ else return kTRUE;
+ }
+
if(pt > 5.) {
// applying cut for pt > 5 GeV
AliDebug(2,Form("pt of D+ = %f (> 5), cutting at |y| < 0.8",pt));
if(pt<fMinPtCand) return 0;
if(pt>fMaxPtCand) return 0;
- if(d->HasBadDaughters()) return 0;
+ if(fUseTrackSelectionWithFilterBits && d->HasBadDaughters()) return 0;
// selection on candidate
if(selectionLevel==AliRDHFCuts::kAll ||
//recalculate vertex w/o daughters
AliAODVertex *origownvtx=0x0;
- if(fRemoveDaughtersFromPrimary) {
+ if(fRemoveDaughtersFromPrimary && !fUseMCVertex) {
if(d->GetOwnPrimaryVtx()) origownvtx=new AliAODVertex(*d->GetOwnPrimaryVtx());
if(!RecalcOwnPrimaryVtx(d,aod)) {
CleanOwnPrimaryVtx(d,aod,origownvtx);
}
}
+ if(fUseMCVertex) {
+ if(d->GetOwnPrimaryVtx()) origownvtx=new AliAODVertex(*d->GetOwnPrimaryVtx());
+ if(!SetMCPrimaryVtx(d,aod)) {
+ CleanOwnPrimaryVtx(d,aod,origownvtx);
+ return 0;
+ }
+ }
+
Int_t ptbin=PtBin(pt);
if (ptbin==-1) {
CleanOwnPrimaryVtx(d,aod,origownvtx);
SetPtBins(nptbins+1,ptbins);
SetCuts(nvars,nptbins,anacutsval);
SetUsePID(kTRUE);
+ fPidHF->SetOldPid(kTRUE);
SetRemoveDaughtersFromPrim(kTRUE);
PrintAll();
delete [] anacutsval;
anacutsval=NULL;
+ delete esdTrackCuts;
+ esdTrackCuts=NULL;
+
return;
}
SetPtBins(nptbins+1,ptbins);
SetCuts(nvars,nptbins,anacutsval);
SetUsePID(kTRUE);
+ fPidHF->SetOldPid(kTRUE);
SetMinCentrality(1E-10);
SetMaxCentrality(20.);
SetUseCentrality(AliRDHFCuts::kCentV0M);
delete [] ptbins;
ptbins=NULL;
+ delete esdTrackCuts;
+ esdTrackCuts=NULL;
+
return;
}
EnableCentralTrigger();
EnableSemiCentralTrigger();
+ // new PID
+ fPidHF->SetOldPid(kFALSE);
+
}
//--------------------------------------------------------------------------