Float_t xDecay=0.;
Float_t yDecay=0.;
Float_t zDecay=0.;
- Bool_t IsInjected = -1;
+ // Bool_t IsInjected = -1;
Bool_t IsLc = 0;
Bool_t IsLcfromLb = 0;
Bool_t IsLcfromc = 0;
Bool_t IsLcFromq = 0;
if(fReadMC){
- AliAODMCHeader *mcHeader2 = (AliAODMCHeader*)aod->GetList()->FindObject(AliAODMCHeader::StdBranchName());
+ // AliAODMCHeader *mcHeader2 = (AliAODMCHeader*)aod->GetList()->FindObject(AliAODMCHeader::StdBranchName());
AliVertexingHFUtils *util = new AliVertexingHFUtils();
- IsInjected = util->IsCandidateInjected(part,mcHeader2,arrayMC)?1:0;
+ // IsInjected = util->IsCandidateInjected(part,mcHeader2,arrayMC)?1:0;
labDp = MatchToMCLambdac(part,arrayMC);
if(labDp>0){
Int_t pdgDgLctopKpi[3]={2212,321,211};
Int_t lab=-9999;
- Bool_t IsLcfromLb=0;
+ // Bool_t IsLcfromLb=0;
Bool_t IsLcfromc=0;
if(fReadMC){
lab=part->MatchToMC(4122,arrMC,3,pdgDgLctopKpi); //return MC particle label if the array corresponds to a Lc, -1 if not (cf. AliAODRecoDecay.cxx)
AliAODMCParticle *partDp = (AliAODMCParticle*)arrMC->At(lab);
AliVertexingHFUtils *util = new AliVertexingHFUtils();
Int_t pdgMom=util->CheckOrigin(arrMC,partDp,kFALSE);
- if(pdgMom == 5) IsLcfromLb =1;
+ // if(pdgMom == 5) IsLcfromLb =1;
if(pdgMom == 4) IsLcfromc =1;
delete util;
}
Double_t theoryRatioExtremeA=1., theoryRatioExtremeB=1.;
Double_t correctionConservativeA=1., correctionConservativeB=1.;
Double_t theoryRatioConservativeA=1., theoryRatioConservativeB=1.;
- Double_t correctionUnc=0.;
+ // Double_t correctionUnc=0.;
Double_t correctionExtremeAUnc=0., correctionExtremeBUnc=0.;
Double_t correctionConservativeAUnc=0., correctionConservativeBUnc=0.;
theoryRatioExtremeA=1.; theoryRatioExtremeB=1.;
correctionConservativeA=1.; correctionConservativeB=1.;
theoryRatioConservativeA=1.; theoryRatioConservativeB=1.;
- correctionUnc=0.;
+ // correctionUnc=0.;
correctionExtremeAUnc=0.; correctionExtremeBUnc=0.;
correctionConservativeAUnc=0.; correctionConservativeBUnc=0.;
correctionConservativeAUncStatEffc=0.; correctionConservativeBUncStatEffc=0.;
(fhDirectEffpt->GetBinError(ibin)/fhDirectEffpt->GetBinContent(ibin))*(fhDirectEffpt->GetBinError(ibin)/fhDirectEffpt->GetBinContent(ibin))
);
- correctionUnc = correction*correction * theoryRatio * effRatio * relEffUnc;
+ // correctionUnc = correction*correction * theoryRatio * effRatio * relEffUnc;
correctionExtremeAUnc = correctionExtremeA*correctionExtremeA * theoryRatioExtremeA * effRatio * relEffUnc;
correctionExtremeBUnc = correctionExtremeB*correctionExtremeB * theoryRatioExtremeB * effRatio * relEffUnc;