diff --git a/PWGJE/Tasks/bjetTreeCreator.cxx b/PWGJE/Tasks/bjetTreeCreator.cxx index ee059aa9626..884a73d0833 100644 --- a/PWGJE/Tasks/bjetTreeCreator.cxx +++ b/PWGJE/Tasks/bjetTreeCreator.cxx @@ -245,6 +245,7 @@ struct BJetTreeCreator { Configurable> jetRadii{"jetRadii", std::vector{0.4}, "jet resolution parameters"}; + Configurable produceSVTree{"produceSVTree", true, "produce the SV-correlated jet TTreeā€"}; Configurable produceTree{"produceTree", true, "produce the jet TTree"}; Configurable vtxRes{"vtxRes", 0.01, "Vertex position resolution (cluster size) for GNN vertex predictions (cm)"}; @@ -415,7 +416,7 @@ struct BJetTreeCreator { double energySV = candSV.e(); if (svIndices.size() < (svReductionFactor * myJet.template tracks_as().size()) && svIndices.size() < maxConstSV) { - if (produceTree) { + if (produceSVTree) { bjetSVParamsTable(bjetParamsTable.lastIndex() + 1, candSV.pt(), deltaRJetSV, massSV, energySV / myJet.energy(), candSV.impactParameterXY(), candSV.cpa(), candSV.chi2PCA(), candSV.dispersion(), candSV.decayLengthXY(), candSV.errorDecayLengthXY(), candSV.decayLength(), candSV.errorDecayLength()); } svIndices.push_back(bjetSVParamsTable.lastIndex()); @@ -755,6 +756,7 @@ struct BJetTreeCreator { } std::vector indicesTracks; + std::vector indicesSVs; int16_t jetFlavor = analysisJet.origin(); @@ -791,7 +793,6 @@ struct BJetTreeCreator { } if (produceTree) { - std::vector indicesSVs; bjetConstituentsTable(bjetParamsTable.lastIndex() + 1, indicesTracks, indicesSVs); bjetParamsTable(analysisJet.pt(), analysisJet.eta(), analysisJet.phi(), indicesTracks.size(), nVertices, analysisJet.mass(), jetFlavor, analysisJet.r()); }