Skip to content

Commit 8b000bb

Browse files
torkjellsdatterIda Torkjellsdatter Storehaug
andauthored
[PWGDQ] [PWGDQ-120] debugging of undefined histogram in AnalysisDileptonTrack (dqEfficiency_withAssoc.cxx) (#8470)
Co-authored-by: Ida Torkjellsdatter Storehaug <idats@eduroam-193-157-237-37.wlan.uio.no>
1 parent 2dfa3d8 commit 8b000bb

File tree

1 file changed

+10
-6
lines changed

1 file changed

+10
-6
lines changed

PWGDQ/Tasks/dqEfficiency_withAssoc.cxx

Lines changed: 10 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -1954,6 +1954,12 @@ struct AnalysisDileptonTrack {
19541954
if (objArraySingleCuts->FindObject(fConfigTrackCut.value.data()) == nullptr) {
19551955
LOG(fatal) << " Track cut chosen for the correlation task was not computed in the single-track task! Check it out!";
19561956
}
1957+
for (int icut = 0; icut < objArraySingleCuts->GetEntries(); ++icut) {
1958+
TString tempStr = objArraySingleCuts->At(icut)->GetName();
1959+
if (tempStr.CompareTo(fConfigTrackCut.value.data()) == 0) {
1960+
fTrackCutBit = icut; // the bit corresponding to the track to be combined with dileptons
1961+
}
1962+
}
19571963
// get the cuts employed for same-event pairing
19581964
string tempCutsPair;
19591965
if (isBarrel) {
@@ -1965,23 +1971,21 @@ struct AnalysisDileptonTrack {
19651971
if (!tempCutsSingleStr.IsNull() && !tempCutsPairStr.IsNull()) {
19661972
std::unique_ptr<TObjArray> objArray(tempCutsPairStr.Tokenize(","));
19671973
fNCuts = objArray->GetEntries();
1968-
for (int icut = 0; icut < objArraySingleCuts->GetEntries(); ++icut) {
1969-
TString tempStr = objArraySingleCuts->At(icut)->GetName();
1974+
for (int icut = 0; icut < objArray->GetEntries(); ++icut) {
1975+
TString tempStr = objArray->At(icut)->GetName();
19701976
if (objArray->FindObject(tempStr.Data()) != nullptr) {
19711977
fHistNamesDileptonTrack[icut] = Form("DileptonTrack_%s_%s", tempStr.Data(), fConfigTrackCut.value.data());
19721978
fHistNamesDileptons[icut] = Form("DileptonsSelected_%s", tempStr.Data());
19731979
DefineHistograms(fHistMan, fHistNamesDileptonTrack[icut], fConfigHistogramSubgroups.value.data()); // define dilepton-track histograms
19741980
DefineHistograms(fHistMan, fHistNamesDileptons[icut], "barrel,vertexing"); // define dilepton histograms
19751981
std::vector<TString> mcHistNames;
1982+
int isig = 0;
19761983
for (auto& sig : fRecMCSignals) {
19771984
mcHistNames.push_back(Form("DileptonTrackMCMatched_%s_%s_%s", tempStr.Data(), fConfigTrackCut.value.data(), sig.GetName()));
1978-
DefineHistograms(fHistMan, mcHistNames[mcHistNames.size() - 1], "");
1985+
DefineHistograms(fHistMan, mcHistNames[mcHistNames.size() - 1], fConfigHistogramSubgroups.value.data());
19791986
}
19801987
fHistNamesDileptonTrackMCmatched[icut] = mcHistNames;
19811988
}
1982-
if (tempStr.CompareTo(fConfigTrackCut.value.data()) == 0) {
1983-
fTrackCutBit = icut; // the bit corresponding to the track to be combined with dileptons
1984-
}
19851989
}
19861990
}
19871991
// Add histogram classes for each specified MCsignal at the generator level

0 commit comments

Comments
 (0)