* provided "as is" without express or implied warranty. *
**************************************************************************/
+/* $Id$ */
+
/////////////////////////////////////////////////////////////
//
// AliAnalysisTaskSE for the comparison of heavy flavor
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;
}
// load D*+ candidates
inputArrayDstar = (TClonesArray*)aodFromExt->GetList()->FindObject("Dstar");
}
- } else {
+ } else if(aod) {
// load HF vertices
inputArrayVertices = (TClonesArray*)aod->GetList()->FindObject("VerticesHF");
// load D0->Kpi candidates
}
- if(!inputArrayVertices) {
+ if(!inputArrayVertices || !aod) {
printf("AliAnalysisTaskSECompareHF::UserExec: Vertices branch not found!\n");
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;