fOutput(0),
fNtupleCmp(0),
fHistMass(0),
-fHistNEvents(0),
-fVHF(0)
+fHistNEvents(0)
{
// Default constructor
fOutput(0),
fNtupleCmp(0),
fHistMass(0),
-fHistNEvents(0),
-fVHF(0)
+fHistNEvents(0)
{
// Standard constructor
delete fOutput;
fOutput = 0;
}
- if (fVHF) {
- delete fVHF;
- fVHF = 0;
- }
}
//________________________________________________________________________
if(fDebug > 1) printf("AnalysisTaskSECompareHF::Init() \n");
- gROOT->LoadMacro("ConfigVertexingHF.C");
-
- fVHF = (AliAnalysisVertexingHF*)gROOT->ProcessLine("ConfigVertexingHF()");
- fVHF->PrintStatus();
-
return;
}
d3->SetOwnPrimaryVtx(vtx1); // needed to compute all variables
unsetvtx=kTRUE;
}
- //if(d3->SelectDplus(fVHF->GetDplusCuts())) {
AliAODMCParticle *dMC = (AliAODMCParticle*)mcArray->At(lab);
pdg = dMC->GetPdgCode();
invmass = d3->InvMassDplus();
(Float_t)d3->CosPointingAngle(),(Float_t)(d3->Getd0Prong(0)*d3->Getd0Prong(1)*d3->Getd0Prong(2))};
fNtupleCmp->Fill(tmp);
PostData(2,fNtupleCmp);
- //}
if(unsetvtx) d3->UnsetOwnPrimaryVtx();
}
break;
unsetvtx=kTRUE;
}
okD0=0; okD0bar=0;
- //if(d4->SelectD0(fVHF->GetD0to4ProngsCuts(),okD0,okD0bar)) {
AliAODMCParticle *dMC = (AliAODMCParticle*)mcArray->At(lab);
pdg = dMC->GetPdgCode();
//invmass = (pdg==421 ? d->InvMassD0() : d->InvMassD0bar());
(Float_t)d4->CosPointingAngle(),(Float_t)(d4->Getd0Prong(0)*d4->Getd0Prong(1)*d4->Getd0Prong(2)*d4->Getd0Prong(3))};
fNtupleCmp->Fill(tmp);
PostData(2,fNtupleCmp);
- //}
if(unsetvtx) d4->UnsetOwnPrimaryVtx();
}
break;
TNtuple *fNtupleCmp; //! output ntuple
TH1F *fHistMass; //! output histogram
TH1F *fHistNEvents; //! output histogram
- AliAnalysisVertexingHF *fVHF; // Vertexer heavy flavour (used to pass the cuts)
- ClassDef(AliAnalysisTaskSECompareHF,4); // AliAnalysisTaskSE for the MC association of heavy-flavour decay candidates
+ ClassDef(AliAnalysisTaskSECompareHF,5); // AliAnalysisTaskSE for the MC association of heavy-flavour decay candidates
};
#endif