for(Int_t i=0; i<9; i++) fD0toKpiCuts[i]=source.fD0toKpiCuts[i];
for(Int_t i=0; i<9; i++) fBtoJPSICuts[i]=source.fBtoJPSICuts[i];
for(Int_t i=0; i<12; i++) fDplusCuts[i]=source.fDplusCuts[i];
- for(Int_t i=0; i<13; i++) fDsCuts[i]=source.fDsCuts[i];
+ for(Int_t i=0; i<14; i++) fDsCuts[i]=source.fDsCuts[i];
for(Int_t i=0; i<12; i++) fLcCuts[i]=source.fLcCuts[i];
for(Int_t i=0; i<5; i++) fDstarCuts[i]=source.fDstarCuts[i];
for(Int_t i=0; i<9; i++) fD0to4ProngsCuts[i]=source.fD0to4ProngsCuts[i];
for(Int_t i=0; i<9; i++) fD0toKpiCuts[i]=source.fD0toKpiCuts[i];
for(Int_t i=0; i<9; i++) fBtoJPSICuts[i]=source.fBtoJPSICuts[i];
for(Int_t i=0; i<12; i++) fDplusCuts[i]=source.fDplusCuts[i];
- for(Int_t i=0; i<13; i++) fDsCuts[i]=source.fDsCuts[i];
+ for(Int_t i=0; i<14; i++) fDsCuts[i]=source.fDsCuts[i];
for(Int_t i=0; i<12; i++) fLcCuts[i]=source.fLcCuts[i];
for(Int_t i=0; i<5; i++) fDstarCuts[i]=source.fDstarCuts[i];
for(Int_t i=0; i<9; i++) fD0to4ProngsCuts[i]=source.fD0to4ProngsCuts[i];
if(f3Prong) {
ok3Prong = kFALSE;
Int_t ok1,ok2;
+ Int_t dum1,dum2;
if(the3Prong->SelectDplus(fDplusCuts)) ok3Prong = kTRUE;
- if(the3Prong->SelectDs(fDsCuts,ok1,ok2)) ok3Prong = kTRUE;
+ if(the3Prong->SelectDs(fDsCuts,ok1,ok2,dum1,dum2)) ok3Prong = kTRUE;
if(the3Prong->SelectLc(fLcCuts,ok1,ok2)) ok3Prong = kTRUE;
}
//if(fDebug) printf("ok3Prong: %d\n",(Int_t)ok3Prong);
printf(" cosThetaPoint > %f\n",fDsCuts[9]);
printf(" Sum d0^2 [cm^2] > %f\n",fDsCuts[10]);
printf(" dca cut [cm] < %f\n",fDsCuts[11]);
- printf(" Inv. Mass phi/K0* [GeV] < %f\n",fDsCuts[12]);
+ printf(" Inv. Mass phi [GeV] < %f\n",fDsCuts[12]);
+ printf(" Inv. Mass K0* [GeV] < %f\n",fDsCuts[13]);
printf(" Lc->pKpi cuts:\n");
printf(" |M-MLc| [GeV] < %f\n",fLcCuts[0]);
printf(" pTP [GeV/c] > %f\n",fLcCuts[1]);
AliVTrack *track;
track = (AliVTrack*)event->GetTrack(i);
+ // skip pure ITS SA tracks
+ //if(track->GetStatus()&AliESDtrack::kITSpureSA) continue;
+
// TEMPORARY: check that the cov matrix is there
Double_t covtest[21];
if(!track->GetCovarianceXYZPxPyPz(covtest)) continue;
}
//-----------------------------------------------------------------------------
void AliAnalysisVertexingHF::SetDsCuts(Double_t cut0,Double_t cut1,
- Double_t cut2,Double_t cut3,Double_t cut4,
- Double_t cut5,Double_t cut6,
- Double_t cut7,Double_t cut8,
- Double_t cut9,Double_t cut10,
- Double_t cut11,Double_t cut12)
+ Double_t cut2,Double_t cut3,
+ Double_t cut4,Double_t cut5,
+ Double_t cut6,Double_t cut7,
+ Double_t cut8,Double_t cut9,
+ Double_t cut10,Double_t cut11,
+ Double_t cut12,Double_t cut13)
{
// Set the cuts for Ds->KKpi selection
fDsCuts[0] = cut0;
fDsCuts[10] = cut10;
fDsCuts[11] = cut11;
fDsCuts[12] = cut12;
+ fDsCuts[13] = cut13;
return;
}
{
// Set the cuts for Ds->KKpi selection
- for(Int_t i=0; i<13; i++) fDsCuts[i] = cuts[i];
+ for(Int_t i=0; i<14; i++) fDsCuts[i] = cuts[i];
return;
}
Double_t cut11=10000000000.);
void SetDplusCuts(const Double_t cuts[12]);
void SetDsCuts(Double_t cut0=1000.,Double_t cut1=0.,
- Double_t cut2=0.,Double_t cut3=0.,Double_t cut4=0.,
- Double_t cut5=0.,Double_t cut6=10000000000.,
- Double_t cut7=0.,Double_t cut8=0.,
- Double_t cut9=-1.1,Double_t cut10=0.,
- Double_t cut11=10000000000., Double_t cut12=1000.);
- void SetDsCuts(const Double_t cuts[13]);
+ Double_t cut2=0.,Double_t cut3=0.,Double_t cut4=0.,
+ Double_t cut5=0.,Double_t cut6=10000000000.,
+ Double_t cut7=0.,Double_t cut8=0.,
+ Double_t cut9=-1.1,Double_t cut10=0.,
+ Double_t cut11=10000000000., Double_t cut12=1000.,
+ Double_t cut13=1000.);
+ void SetDsCuts(const Double_t cuts[14]);
void SetLcCuts(Double_t cut0=1000.,Double_t cut1=0.,
Double_t cut2=0.,Double_t cut3=0.,Double_t cut4=0.,
Double_t cut5=0.,Double_t cut6=10000000000.,
// 9 = cosThetaPoint
// 10 = Sum d0^2 (cm^2)
// 11 = dca cut (cm)
- Double_t fDsCuts[13]; // cuts on Ds candidates
+ Double_t fDsCuts[14]; // cuts on Ds candidates
// (to be passed to AliAODRecoDecayHF3Prong::SelectDs())
// 0 = inv. mass half width [GeV]
// 1 = pTK [GeV/c]
// 9 = cosThetaPoint
// 10 = Sum d0^2 (cm^2)
// 11 = dca cut (cm)
- // 12 = inv. mass cut arounf phi and K0* [GeV]
+ // 12 = inv. mass cut around phi
+ // 13 = inv. mass cut around phi K0* [GeV]
Double_t fLcCuts[12]; // cuts on Lambdac candidates
// (to be passed to AliAODRecoDecayHF3Prong::SelectLc())
// 0 = inv. mass half width [GeV]
void SetPrimaryVertex(AliESDVertex *v1) { fV1 = v1; }
Bool_t SingleTrkCuts(AliESDtrack *trk,Bool_t &okDisplaced,Bool_t &okSoftPi) const;
//
- ClassDef(AliAnalysisVertexingHF,12); // Reconstruction of HF decay candidates
+ ClassDef(AliAnalysisVertexingHF,13); // Reconstruction of HF decay candidates
};