AliAODRecoDecayHF3Prong *dd = (AliAODRecoDecayHF3Prong*)d;
- Int_t iter=-1;
+ Int_t iter=-1;
if(fVarsForOpt[0]){
iter++;
vars[iter]=dd->InvMassLcpKpi();
//
if(!fCutsRD){
- cout<<"Cut matrice not inizialized. Exit..."<<endl;
+ AliError("Cut matrice not inizialized. Exit...\n");
return 0;
}
//PrintAll();
AliAODRecoDecayHF3Prong* d=(AliAODRecoDecayHF3Prong*)obj;
if(!d){
- cout<<"AliAODRecoDecayHF3Prong null"<<endl;
+ AliError("AliAODRecoDecayHF3Prong null \n");
return 0;
}
if(fUseTrackSelectionWithFilterBits && d->HasBadDaughters()) return 0;
+
+ // selection on daughter tracks
+ if(selectionLevel==AliRDHFCuts::kAll ||
+ selectionLevel==AliRDHFCuts::kTracks) {
+ if(!AreDaughtersSelected(d)) return 0;
+ }
+
+
+ // PID selection
if(selectionLevel==AliRDHFCuts::kAll ||
selectionLevel==AliRDHFCuts::kCandidate||
selectionLevel==AliRDHFCuts::kPID) {
- switch (fPIDStrategy) {
- case kNSigma:
- returnvaluePID = IsSelectedPID(d);
-
+ switch (fPIDStrategy) {
+ case kNSigma:
+ returnvaluePID = IsSelectedPID(d);
break;
- case kCombined:
- returnvaluePID = IsSelectedCombinedPID(d);
+ case kCombined:
+ returnvaluePID = IsSelectedCombinedPID(d);
break;
- case kCombinedSoft:
- returnvaluePID = IsSelectedCombinedPIDSoft(d);
+ case kCombinedSoft:
+ returnvaluePID = IsSelectedCombinedPIDSoft(d);
break;
- case kNSigmaStrong:
- returnvaluePID = IsSelectedPIDStrong(d);
- }
- fIsSelectedPID=returnvaluePID;
+ case kNSigmaStrong:
+ returnvaluePID = IsSelectedPIDStrong(d);
+ break;
+ }
+ fIsSelectedPID=returnvaluePID;
}
// if(fUsePID || selectionLevel==AliRDHFCuts::kPID) returnvaluePID = IsSelectedCombinedPID(d); // to test!!
if(returnvaluePID==0) return 0;
+
+
+
// selection on candidate
if(selectionLevel==AliRDHFCuts::kAll ||
selectionLevel==AliRDHFCuts::kCandidate) {
Int_t ptbin=PtBin(pt);
- Double_t mLcpKpi,mLcpiKp;
+ Double_t mLcpKpi=0.,mLcpiKp=0.;
Int_t okLcpKpi=1,okLcpiKp=1;
Double_t mLcPDG = TDatabasePDG::Instance()->GetParticle(4122)->Mass();
}
- // selection on daughter tracks
- if(selectionLevel==AliRDHFCuts::kAll ||
- selectionLevel==AliRDHFCuts::kTracks) {
- if(!AreDaughtersSelected(d)) return 0;
- }
-
Int_t returnvalueTot=CombinePIDCuts(returnvalue,returnvaluePID);
return returnvalueTot;
delete [] ptbins;
ptbins=NULL;
+ delete pidObjK;
+ pidObjK=NULL;
+ delete pidObjpi;
+ pidObjpi=NULL;
delete pidObjp;
pidObjp=NULL;
esdTrackCuts->SetEtaRange(-0.8,0.8);
esdTrackCuts->SetMaxDCAToVertexXY(1.);
esdTrackCuts->SetMaxDCAToVertexZ(1.);
- esdTrackCuts->SetPtRange(0.8,1.e10);
+ esdTrackCuts->SetPtRange(0.49,1.e10);
AddTrackCuts(esdTrackCuts);
const Int_t nptbins=4;
delete [] ptbins;
ptbins=NULL;
+ delete pidObjK;
+ pidObjK=NULL;
+ delete pidObjpi;
+ pidObjpi=NULL;
delete pidObjp;
pidObjp=NULL;