Skip to content

Commit

Permalink
[PWGDQ] use float type for statistic histograms (AliceO2Group#6802)
Browse files Browse the repository at this point in the history 
* [PWGDQ] use float type for statistic histograms

* [PWGDQ] use double type for statistic histograms

* [PWGDQ] use double type for statistic histogram in filterPbPb
  • Loading branch information
@feisenhu
feisenhu authored Jul 10, 2024
1 parent fca5c25 commit 2bd0e14
Show file tree
Hide file tree
Showing 5 changed files with 59 additions and 61 deletions.
70 changes: 35 additions & 35 deletions PWGDQ/TableProducer/tableMaker.cxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -416,10 +416,10 @@ struct TableMaker {
// fill stats information, before selections
for (int i = 0; i < kNaliases; i++) {
if (triggerAliases & (uint32_t(1) << i)) {
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(2.0, static_cast<float>(i));
(reinterpret_cast<TH2F*>(fStatsList->At(0)))->Fill(2.0, static_cast<float>(i));
}
}
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(2.0, static_cast<float>(kNaliases));
(reinterpret_cast<TH2F*>(fStatsList->At(0)))->Fill(2.0, static_cast<float>(kNaliases));

if (fConfigRunZorro) {
for (int i = 0; i < kNTriggersDQ; ++i) {
Expand All @@ -437,10 +437,10 @@ struct TableMaker {
// fill stats information, after selections
for (int i = 0; i < kNaliases; i++) {
if (triggerAliases & (uint32_t(1) << i)) {
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(3.0, static_cast<float>(i));
(reinterpret_cast<TH2F*>(fStatsList->At(0)))->Fill(3.0, static_cast<float>(i));
}
}
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(3.0, static_cast<float>(kNaliases));
(reinterpret_cast<TH2F*>(fStatsList->At(0)))->Fill(3.0, static_cast<float>(kNaliases));

fHistMan->FillHistClass("Event_AfterCuts", VarManager::fgValues);

Expand Down Expand Up @@ -519,7 +519,7 @@ struct TableMaker {
fHistMan->FillHistClass(Form("Ambiguous_TrackBarrel_%s", (*cut).GetName()), VarManager::fgValues);
}
}
(reinterpret_cast<TH1I*>(fStatsList->At(1)))->Fill(static_cast<float>(i));
(reinterpret_cast<TH1D*>(fStatsList->At(1)))->Fill(static_cast<float>(i));
}
}
if (!trackTempFilterMap) {
Expand All @@ -537,7 +537,7 @@ struct TableMaker {
trackFilteringTag = uint64_t(track.pidbit());
for (int iv0 = 0; iv0 < 5; iv0++) {
if (track.pidbit() & (uint8_t(1) << iv0)) {
(reinterpret_cast<TH1I*>(fStatsList->At(1)))->Fill(fTrackCuts.size() + static_cast<float>(iv0));
(reinterpret_cast<TH1D*>(fStatsList->At(1)))->Fill(fTrackCuts.size() + static_cast<float>(iv0));
}
}
if (fConfigIsOnlyforMaps) {
Expand Down Expand Up @@ -719,7 +719,7 @@ struct TableMaker {
fHistMan->FillHistClass(Form("Ambiguous_Muons_%s", (*cut).GetName()), VarManager::fgValues);
}
}
(reinterpret_cast<TH1I*>(fStatsList->At(2)))->Fill(static_cast<float>(i));
(reinterpret_cast<TH1D*>(fStatsList->At(2)))->Fill(static_cast<float>(i));
}
}
if (!trackTempFilterMap) {
Expand Down Expand Up @@ -864,10 +864,10 @@ struct TableMaker {
// fill stats information, before selections
for (int i = 0; i < kNaliases; i++) {
if (triggerAliases & (uint32_t(1) << i)) {
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(2.0, static_cast<float>(i));
(reinterpret_cast<TH2F*>(fStatsList->At(0)))->Fill(2.0, static_cast<float>(i));
}
}
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(2.0, static_cast<float>(kNaliases));
(reinterpret_cast<TH2F*>(fStatsList->At(0)))->Fill(2.0, static_cast<float>(kNaliases));

if (fConfigRunZorro) {
for (int i = 0; i < kNTriggersDQ; ++i) {
Expand All @@ -885,10 +885,10 @@ struct TableMaker {
// fill stats information, after selections
for (int i = 0; i < kNaliases; i++) {
if (triggerAliases & (uint32_t(1) << i)) {
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(3.0, static_cast<float>(i));
(reinterpret_cast<TH2F*>(fStatsList->At(0)))->Fill(3.0, static_cast<float>(i));
}
}
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(3.0, static_cast<float>(kNaliases));
(reinterpret_cast<TH2F*>(fStatsList->At(0)))->Fill(3.0, static_cast<float>(kNaliases));

fHistMan->FillHistClass("Event_AfterCuts", VarManager::fgValues);

Expand Down Expand Up @@ -953,7 +953,7 @@ struct TableMaker {
fHistMan->FillHistClass(Form("Ambiguous_TrackBarrel_%s", (*cut).GetName()), VarManager::fgValues);
}
}
(reinterpret_cast<TH1I*>(fStatsList->At(1)))->Fill(static_cast<float>(i));
(reinterpret_cast<TH1D*>(fStatsList->At(1)))->Fill(static_cast<float>(i));
}
}
if (!trackTempFilterMap) {
Expand All @@ -971,7 +971,7 @@ struct TableMaker {
trackFilteringTag |= (uint64_t(track.pidbit()) << 2);
for (int iv0 = 0; iv0 < 5; iv0++) {
if (track.pidbit() & (uint8_t(1) << iv0)) {
(reinterpret_cast<TH1I*>(fStatsList->At(1)))->Fill(fTrackCuts.size() + static_cast<float>(iv0));
(reinterpret_cast<TH1D*>(fStatsList->At(1)))->Fill(fTrackCuts.size() + static_cast<float>(iv0));
}
}
if (fConfigIsOnlyforMaps) {
Expand Down Expand Up @@ -1101,7 +1101,7 @@ struct TableMaker {
fHistMan->FillHistClass(Form("Ambiguous_Muons_%s", (*cut).GetName()), VarManager::fgValues);
}
}
(reinterpret_cast<TH1I*>(fStatsList->At(2)))->Fill(static_cast<float>(i));
(reinterpret_cast<TH1D*>(fStatsList->At(2)))->Fill(static_cast<float>(i));
}
}
if (!trackTempFilterMap) {
Expand Down Expand Up @@ -1216,7 +1216,7 @@ struct TableMaker {
fStatsList.setObject(new TList());
fStatsList->SetOwner(kTRUE);
std::vector<TString> eventLabels{"BCs", "Collisions before filtering", "Before cuts", "After cuts"};
TH2I* histEvents = new TH2I("EventStats", "Event statistics", eventLabels.size(), -0.5, eventLabels.size() - 0.5, kNaliases + 1, -0.5, +kNaliases + 0.5);
TH2F* histEvents = new TH2F("EventStats", "Event statistics", eventLabels.size(), -0.5, eventLabels.size() - 0.5, kNaliases + 1, -0.5, +kNaliases + 0.5);
int ib = 1;
for (auto label = eventLabels.begin(); label != eventLabels.end(); label++, ib++) {
histEvents->GetXaxis()->SetBinLabel(ib, (*label).Data());
Expand All @@ -1228,7 +1228,7 @@ struct TableMaker {
fStatsList->Add(histEvents);

// Track statistics: one bin for each track selection and 5 bins for V0 tags (gamma, K0s, Lambda, anti-Lambda, Omega)
TH1I* histTracks = new TH1I("TrackStats", "Track statistics", fTrackCuts.size() + 5.0, -0.5, fTrackCuts.size() - 0.5 + 5.0);
TH1D* histTracks = new TH1D("TrackStats", "Track statistics", fTrackCuts.size() + 5.0, -0.5, fTrackCuts.size() - 0.5 + 5.0);
ib = 1;
for (auto cut = fTrackCuts.begin(); cut != fTrackCuts.end(); cut++, ib++) {
histTracks->GetXaxis()->SetBinLabel(ib, (*cut).GetName());
Expand All @@ -1238,7 +1238,7 @@ struct TableMaker {
histTracks->GetXaxis()->SetBinLabel(fTrackCuts.size() + 1 + ib, v0TagNames[ib]);
}
fStatsList->Add(histTracks);
TH1I* histMuons = new TH1I("MuonStats", "Muon statistics", fMuonCuts.size(), -0.5, fMuonCuts.size() - 0.5);
TH1D* histMuons = new TH1D("MuonStats", "Muon statistics", fMuonCuts.size(), -0.5, fMuonCuts.size() - 0.5);
ib = 1;
for (auto cut = fMuonCuts.begin(); cut != fMuonCuts.end(); cut++, ib++) {
histMuons->GetXaxis()->SetBinLabel(ib, (*cut).GetName());
Expand Down Expand Up @@ -1286,10 +1286,10 @@ struct TableMaker {
{
for (int i = 0; i < kNaliases; i++) {
if (collision.alias_bit(i) > 0) {
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(i));
(reinterpret_cast<TH2F*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(i));
}
}
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(kNaliases));
(reinterpret_cast<TH2F*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(kNaliases));
if (collision.eventFilter()) {
fullSkimming<gkEventFillMap, gkTrackFillMapWithCov, gkMuonFillMapWithCov>(collision, bcs, tracksBarrel, tracksMuon, nullptr, nullptr);
}
Expand All @@ -1301,10 +1301,10 @@ struct TableMaker {
{
for (int i = 0; i < kNaliases; i++) {
if (collision.alias_bit(i) > 0) {
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(i));
(reinterpret_cast<TH2F*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(i));
}
}
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(kNaliases));
(reinterpret_cast<TH2F*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(kNaliases));
if (collision.eventFilter()) {
fullSkimming<gkEventFillMapWithMult, gkTrackFillMapWithCov, gkMuonFillMapWithCov>(collision, bcs, tracksBarrel, tracksMuon, nullptr, nullptr);
}
Expand Down Expand Up @@ -1367,10 +1367,10 @@ struct TableMaker {
{
for (int i = 0; i < kNaliases; i++) {
if (collision.alias_bit(i) > 0) {
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(i));
(reinterpret_cast<TH2F*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(i));
}
}
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(kNaliases));
(reinterpret_cast<TH2F*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(kNaliases));
if (collision.eventFilter()) {
fullSkimming<gkEventFillMap, gkTrackFillMap, 0u>(collision, bcs, tracksBarrel, nullptr, nullptr, nullptr);
}
Expand All @@ -1389,10 +1389,10 @@ struct TableMaker {
{
for (int i = 0; i < kNaliases; i++) {
if (collision.alias_bit(i) > 0) {
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(i));
(reinterpret_cast<TH2F*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(i));
}
}
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(kNaliases));
(reinterpret_cast<TH2F*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(kNaliases));
if (collision.eventFilter()) {
fullSkimming<gkEventFillMapWithMultsAndEventFilter, gkTrackFillMap, 0u>(collision, bcs, tracksBarrel, nullptr, nullptr, nullptr);
}
Expand All @@ -1404,10 +1404,10 @@ struct TableMaker {
{
for (int i = 0; i < kNaliases; i++) {
if (collision.alias_bit(i) > 0) {
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(i));
(reinterpret_cast<TH2F*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(i));
}
}
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(kNaliases));
(reinterpret_cast<TH2F*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(kNaliases));
if (collision.eventFilter()) {
fullSkimming<gkEventFillMap, gkTrackFillMapWithCov, 0u>(collision, bcs, tracksBarrel, nullptr, nullptr, nullptr);
}
Expand Down Expand Up @@ -1475,10 +1475,10 @@ struct TableMaker {
{
for (int i = 0; i < kNaliases; i++) {
if (collision.alias_bit(i) > 0) {
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(i));
(reinterpret_cast<TH2F*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(i));
}
}
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(kNaliases));
(reinterpret_cast<TH2F*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(kNaliases));
if (collision.eventFilter()) {
fullSkimming<gkEventFillMap, 0u, gkMuonFillMapWithCov>(collision, bcs, nullptr, tracksMuon, nullptr, nullptr);
}
Expand Down Expand Up @@ -1583,10 +1583,10 @@ struct TableMaker {
{
for (int i = 0; i < kNaliases; i++) {
if (bc.alias_bit(i) > 0) {
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(0.0, static_cast<float>(i));
(reinterpret_cast<TH2F*>(fStatsList->At(0)))->Fill(0.0, static_cast<float>(i));
}
}
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(0.0, static_cast<float>(kNaliases));
(reinterpret_cast<TH2F*>(fStatsList->At(0)))->Fill(0.0, static_cast<float>(kNaliases));
}

// List of process functions removed because they are not using cov matrix
Expand Down Expand Up @@ -1625,10 +1625,10 @@ struct TableMaker {
{
for (int i = 0; i < kNaliases; i++) {
if (collision.alias_bit(i) > 0) {
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(i));
(reinterpret_cast<TH2F*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(i));
}
}
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(kNaliases));
(reinterpret_cast<TH2F*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(kNaliases));
if (collision.eventFilter()) {
fullSkimming<gkEventFillMap, 0u, gkMuonFillMap>(collision, bcs, nullptr, tracksMuon, nullptr, nullptr);
}
Expand All @@ -1640,10 +1640,10 @@ struct TableMaker {
{
for (int i = 0; i < kNaliases; i++) {
if (collision.alias_bit(i) > 0) {
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(i));
(reinterpret_cast<TH2F*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(i));
}
}
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(kNaliases));
(reinterpret_cast<TH2F*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(kNaliases));
if (collision.eventFilter()) {
fullSkimming<gkEventFillMapWithMultsAndEventFilter, 0u, gkMuonFillMap>(collision, bcs, nullptr, tracksMuon, nullptr, nullptr);
}
Expand Down
28 changes: 14 additions & 14 deletions PWGDQ/TableProducer/tableMaker_withAssoc.cxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -412,7 +412,7 @@ struct TableMaker {
fStatsList.setObject(new TList());
fStatsList->SetOwner(kTRUE);
std::vector<TString> eventLabels{"BCs", "Collisions before filtering", "Before cuts", "After cuts"};
TH2I* histEvents = new TH2I("EventStats", "Event statistics", eventLabels.size(), -0.5, eventLabels.size() - 0.5, o2::aod::evsel::kNsel + 1, -0.5, static_cast<double>(o2::aod::evsel::kNsel) + 0.5);
TH2D* histEvents = new TH2D("EventStats", "Event statistics", eventLabels.size(), -0.5, eventLabels.size() - 0.5, o2::aod::evsel::kNsel + 1, -0.5, static_cast<double>(o2::aod::evsel::kNsel) + 0.5);
int ib = 1;
for (auto label = eventLabels.begin(); label != eventLabels.end(); label++, ib++) {
histEvents->GetXaxis()->SetBinLabel(ib, (*label).Data());
Expand All @@ -424,7 +424,7 @@ struct TableMaker {
fStatsList->Add(histEvents);

// Track statistics: one bin for each track selection and 5 bins for V0 tags (gamma, K0s, Lambda, anti-Lambda, Omega)
TH1I* histTracks = new TH1I("TrackStats", "Track statistics", fTrackCuts.size() + 5.0, -0.5, fTrackCuts.size() - 0.5 + 5.0);
TH1D* histTracks = new TH1D("TrackStats", "Track statistics", fTrackCuts.size() + 5.0, -0.5, fTrackCuts.size() - 0.5 + 5.0);
ib = 1;
for (auto cut = fTrackCuts.begin(); cut != fTrackCuts.end(); cut++, ib++) {
histTracks->GetXaxis()->SetBinLabel(ib, (*cut).GetName());
Expand All @@ -434,7 +434,7 @@ struct TableMaker {
histTracks->GetXaxis()->SetBinLabel(fTrackCuts.size() + 1 + ib, v0TagNames[ib]);
}
fStatsList->Add(histTracks);
TH1I* histMuons = new TH1I("MuonStats", "Muon statistics", fMuonCuts.size(), -0.5, fMuonCuts.size() - 0.5);
TH1D* histMuons = new TH1D("MuonStats", "Muon statistics", fMuonCuts.size(), -0.5, fMuonCuts.size() - 0.5);
ib = 1;
for (auto cut = fMuonCuts.begin(); cut != fMuonCuts.end(); cut++, ib++) {
histMuons->GetXaxis()->SetBinLabel(ib, (*cut).GetName());
Expand Down Expand Up @@ -467,10 +467,10 @@ struct TableMaker {

for (int i = 0; i < o2::aod::evsel::kNsel; i++) {
if (collision.selection_bit(i)) {
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(i));
(reinterpret_cast<TH2D*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(i));
}
}
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(o2::aod::evsel::kNsel));
(reinterpret_cast<TH2D*>(fStatsList->At(0)))->Fill(1.0, static_cast<float>(o2::aod::evsel::kNsel));

// apply the event filter computed by filter-PP
if constexpr ((TEventFillMap & VarManager::ObjTypes::EventFilter) > 0) {
Expand Down Expand Up @@ -509,10 +509,10 @@ struct TableMaker {
// fill stats information, before selections
for (int i = 0; i < o2::aod::evsel::kNsel; i++) {
if (collision.selection_bit(i)) {
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(2.0, static_cast<float>(i));
(reinterpret_cast<TH2D*>(fStatsList->At(0)))->Fill(2.0, static_cast<float>(i));
}
}
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(2.0, static_cast<float>(o2::aod::evsel::kNsel));
(reinterpret_cast<TH2D*>(fStatsList->At(0)))->Fill(2.0, static_cast<float>(o2::aod::evsel::kNsel));

if (fConfigRunZorro) {
for (int i = 0; i < kNTriggersDQ; ++i) {
Expand All @@ -530,10 +530,10 @@ struct TableMaker {
// fill stats information, after selections
for (int i = 0; i < o2::aod::evsel::kNsel; i++) {
if (collision.selection_bit(i)) {
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(3.0, static_cast<float>(i));
(reinterpret_cast<TH2D*>(fStatsList->At(0)))->Fill(3.0, static_cast<float>(i));
}
}
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(3.0, static_cast<float>(o2::aod::evsel::kNsel));
(reinterpret_cast<TH2D*>(fStatsList->At(0)))->Fill(3.0, static_cast<float>(o2::aod::evsel::kNsel));

fHistMan->FillHistClass("Event_AfterCuts", VarManager::fgValues);

Expand Down Expand Up @@ -610,7 +610,7 @@ struct TableMaker {
if (fConfigQA) {
fHistMan->FillHistClass(Form("TrackBarrel_%s", (*cut).GetName()), VarManager::fgValues);
}
(reinterpret_cast<TH1I*>(fStatsList->At(1)))->Fill(static_cast<float>(i));
(reinterpret_cast<TH1D*>(fStatsList->At(1)))->Fill(static_cast<float>(i));
}
}
if (!trackTempFilterMap) {
Expand All @@ -622,7 +622,7 @@ struct TableMaker {
trackFilteringTag |= uint64_t(track.pidbit());
for (int iv0 = 0; iv0 < 5; iv0++) {
if (track.pidbit() & (uint8_t(1) << iv0)) {
(reinterpret_cast<TH1I*>(fStatsList->At(1)))->Fill(fTrackCuts.size() + static_cast<float>(iv0));
(reinterpret_cast<TH1D*>(fStatsList->At(1)))->Fill(fTrackCuts.size() + static_cast<float>(iv0));
}
}
if (fConfigIsOnlyforMaps) {
Expand Down Expand Up @@ -769,7 +769,7 @@ struct TableMaker {
if (fConfigQA) {
fHistMan->FillHistClass(Form("Muons_%s", (*cut).GetName()), VarManager::fgValues);
}
(reinterpret_cast<TH1I*>(fStatsList->At(2)))->Fill(static_cast<float>(i));
(reinterpret_cast<TH1D*>(fStatsList->At(2)))->Fill(static_cast<float>(i));
}
}

Expand Down Expand Up @@ -999,10 +999,10 @@ struct TableMaker {
{
for (int i = 0; i < o2::aod::evsel::kNsel; i++) {
if (bc.selection_bit(i) > 0) {
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(0.0, static_cast<float>(i));
(reinterpret_cast<TH2D*>(fStatsList->At(0)))->Fill(0.0, static_cast<float>(i));
}
}
(reinterpret_cast<TH2I*>(fStatsList->At(0)))->Fill(0.0, static_cast<float>(o2::aod::evsel::kNsel));
(reinterpret_cast<TH2D*>(fStatsList->At(0)))->Fill(0.0, static_cast<float>(o2::aod::evsel::kNsel));
}

PROCESS_SWITCH(TableMaker, processPPWithFilter, "Build full DQ skimmed data model typically for pp/p-Pb and UPC Pb-Pb, w/ event filtering", false);
Expand Down
7 changes: 3 additions & 4 deletions PWGDQ/Tasks/filterPP.cxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -15,8 +15,7 @@
#include <vector>
#include <memory>
#include <cstring>
#include <TH1F.h>
#include <TH2I.h>
#include <TH1.h>
#include <THashList.h>
#include <TString.h>
#include "Framework/AnalysisTask.h"
Expand Down Expand Up @@ -378,7 +377,7 @@ struct DQFilterPPTask {
Produces<aod::DQEventFilter> eventFilter;
Produces<aod::DqFilters> dqtable;
OutputObj<THashList> fOutputList{"output"};
OutputObj<TH1I> fStats{"Statistics"};
OutputObj<TH1D> fStats{"Statistics"};
HistogramManager* fHistMan;

Configurable<std::string> fConfigBarrelSelections{"cfgBarrelSels", "jpsiPID1:pairMassLow:1", "<track-cut>:[<pair-cut>]:<n>,[<track-cut>:[<pair-cut>]:<n>],..."};
Expand Down Expand Up @@ -450,7 +449,7 @@ struct DQFilterPPTask {
VarManager::SetUseVars(AnalysisCut::fgUsedVars);

// setup the Stats histogram
fStats.setObject(new TH1I("Statistics", "Stats for DQ triggers", fNBarrelCuts + fNMuonCuts + 2, -2.5, -0.5 + fNBarrelCuts + fNMuonCuts));
fStats.setObject(new TH1D("Statistics", "Stats for DQ triggers", fNBarrelCuts + fNMuonCuts + 2, -2.5, -0.5 + fNBarrelCuts + fNMuonCuts));
fStats->GetXaxis()->SetBinLabel(1, "Events inspected");
fStats->GetXaxis()->SetBinLabel(2, "Events selected");
if (fNBarrelCuts) {
Expand Down
Loading

0 comments on commit 2bd0e14

Please sign in to comment.