Skip to content

Commit 8ec6fe3

Browse files
committed
fixed Sumw2 issue + other minor bugs
1 parent c264435 commit 8ec6fe3

File tree

1 file changed

+4
-4
lines changed

1 file changed

+4
-4
lines changed

PWGJE/Tasks/fullJetSpectra.cxx

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -135,7 +135,7 @@ struct FullJetSpectra {
135135
// Instantiate the Zorro processor for skimmed data and define an output object
136136
Zorro zorro;
137137
OutputObj<ZorroSummary> zorroSummary{"zorroSummary"};
138-
const bool doSumw2 = doMBGapTrigger;
138+
bool doSumw2 = false;
139139

140140
// Multiplicity Utilities
141141
// struct CentClass {
@@ -334,7 +334,7 @@ struct FullJetSpectra {
334334
triggerMaskBits = jetderiveddatautilities::initialiseTriggerMaskBits(triggerMasks);
335335
particleSelection = static_cast<std::string>(particleSelections);
336336
jetRadiiValues = (std::vector<double>)jetRadii;
337-
337+
doSumw2 = doMBGapTrigger;
338338
/* if (doMcClosure) {
339339
// randGen.SetSeed(mcSplitSeed);
340340
// randGen.SetSeed(static_cast<UInt_t>(std::time(nullptr)));
@@ -508,7 +508,7 @@ struct FullJetSpectra {
508508
registry.add("h2_full_jet_energyscaleChargedVsFullPart", "Jet Energy Scale (charged part, vs. full jet pt); p_{T,part} (GeV/c); (p_{T,det} - p_{T,part})/p_{T,part}", {HistType::kTH2F, {{400, 0., 400.}, {200, -1., 1.}}}, doSumw2);
509509
registry.add("h2_full_jet_energyscaleNeutralVsFullPart", "Jet Energy Scale (neutral part, vs. full jet pt); p_{T,part} (GeV/c); (p_{T,det} - p_{T,part})/p_{T,part}", {HistType::kTH2F, {{400, 0., 400.}, {200, -1., 1.}}}, doSumw2);
510510
registry.add("h2_full_fakemcdjets", "Fake MCD Jets; p_{T,det} (GeV/c); NCounts", {HistType::kTH2F, {{350, 0., 350.}, {100, 0., 100.}}}, doSumw2);
511-
registry.add("h2FullfakeMcpJets", "Fake MCP Jets; p_{T,part} (GeV/c); NCounts", {HistType::kTH2F, {{350, 0., 350.}, {100, 0., 100.}}}, doSumw2);
511+
registry.add("h2_full_fakemcpjets", "Fake MCP Jets; p_{T,part} (GeV/c); NCounts", {HistType::kTH2F, {{350, 0., 350.}, {100, 0., 100.}}}, doSumw2);
512512
registry.add("h2_full_matchedmcpjet_pt", "Matched MCP jet in EMC Fiducial Acceptance #it{p}_{T,part};#it{p}_{T,part} (GeV/c); Ncounts", {HistType::kTH2F, {{350, 0., 350.}, {10000, 0., 10000.}}}, doSumw2);
513513

514514
// Response Matrix
@@ -1740,7 +1740,7 @@ struct FullJetSpectra {
17401740
// apply emcal fiducial cuts to the matched particle level jets - if the matched mcp jet lies outside of the EMCAL fiducial, flag it as a fake jet
17411741
if (mcpjet.eta() > jetEtaMax || mcpjet.eta() < jetEtaMin || mcpjet.phi() > jetPhiMax || mcpjet.phi() < jetPhiMin) {
17421742
fakeMcpJet++;
1743-
registry.fill(HIST("h2FullfakeMcpJets"), mcpjet.pt(), fakeMcpJet, eventWeight);
1743+
registry.fill(HIST("h2_full_fakemcpjets"), mcpjet.pt(), fakeMcpJet, eventWeight);
17441744
continue;
17451745
} else {
17461746
NPartJetFid++;

0 commit comments

Comments
 (0)