Skip to content
Merged
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
15 changes: 6 additions & 9 deletions PWGHF/HFC/Tasks/taskCorrelationDplusHadrons.cxx
Original file line number Diff line number Diff line change
@@ -1,3 +1,3 @@
// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
Expand All @@ -10,7 +10,7 @@
// or submit itself to any jurisdiction.

/// \file taskCorrelationDplusHadrons.cxx
/// \author Shyam Kumar <shyam.kumar@cern.ch>

Check warning on line 13 in PWGHF/HFC/Tasks/taskCorrelationDplusHadrons.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[doc/file]

Documentation for \brief is missing, incorrect or misplaced.
#include <memory> // std::shared_ptr
#include <string>
#include <vector>
Expand Down Expand Up @@ -48,15 +48,15 @@

const int npTBinsCorrelations = 8;
const double pTBinsCorrelations[npTBinsCorrelations + 1] = {0., 2., 4., 6., 8., 12., 16., 24., 99.};
auto pTBinsCorrelations_v = std::vector<double>{pTBinsCorrelations, pTBinsCorrelations + npTBinsCorrelations + 1};

Check warning on line 51 in PWGHF/HFC/Tasks/taskCorrelationDplusHadrons.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/function-variable]

Use lowerCamelCase for names of functions and variables.
const double signalRegionInnerDefault[npTBinsCorrelations] = {1.8490, 1.8490, 1.8490, 1.8490, 1.8490, 1.8490, 1.8490, 1.8490};
const double signalRegionOuterDefault[npTBinsCorrelations] = {1.8890, 1.8890, 1.8890, 1.8890, 1.8890, 1.8890, 1.8890, 1.8890};
const double sidebandLeftOuterDefault[npTBinsCorrelations] = {1.7690, 1.7690, 1.7690, 1.7690, 1.7690, 1.7690, 1.7690, 1.7690};
const double sidebandLeftInnerDefault[npTBinsCorrelations] = {1.8250, 1.8250, 1.8250, 1.8250, 1.8250, 1.8250, 1.8250, 1.8250};
const double sidebandRightInnerDefault[npTBinsCorrelations] = {1.9130, 1.9130, 1.9130, 1.9130, 1.9130, 1.9130, 1.9130, 1.9130};
const double sidebandRightOuterDefault[npTBinsCorrelations] = {1.9690, 1.9690, 1.9690, 1.9690, 1.9690, 1.9690, 1.9690, 1.9690};
auto signalRegionInner_v = std::vector<double>{signalRegionInnerDefault, signalRegionInnerDefault + npTBinsCorrelations};

Check warning on line 58 in PWGHF/HFC/Tasks/taskCorrelationDplusHadrons.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/function-variable]

Use lowerCamelCase for names of functions and variables.
auto signalRegionOuter_v = std::vector<double>{signalRegionOuterDefault, signalRegionOuterDefault + npTBinsCorrelations};

Check warning on line 59 in PWGHF/HFC/Tasks/taskCorrelationDplusHadrons.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/function-variable]

Use lowerCamelCase for names of functions and variables.
auto sidebandLeftInner_v = std::vector<double>{sidebandLeftInnerDefault, sidebandLeftInnerDefault + npTBinsCorrelations};
auto sidebandLeftOuter_v = std::vector<double>{sidebandLeftOuterDefault, sidebandLeftOuterDefault + npTBinsCorrelations};
auto sidebandRightInner_v = std::vector<double>{sidebandRightInnerDefault, sidebandRightInnerDefault + npTBinsCorrelations};
Expand Down Expand Up @@ -276,7 +276,7 @@
}
auto effD = isPromptAnalysis ? mEfficiencyPrompt : mEfficiencyFD;
idxBdtScore = isPromptAnalysis ? 1 : 2;

if (activateQA) {
const int regionLimits = 6;
std::string labels[regionLimits] = {"SigReg Left", "SigReg Right", "Left SB Low", "Left SB Up", "Right SB Low", "Right SB Up"};
Expand Down Expand Up @@ -305,7 +305,7 @@
float bdtScoreFd = candidate.mlScoreFd();
float bdtScoreBkg = candidate.mlScoreBkg();
int effBinD = o2::analysis::findBin(binsPtEfficiencyD, ptD);
float bdtScorePromptorFd = isPromptAnalysis ? bdtScorePrompt: bdtScoreFd;
float bdtScorePromptorFd = isPromptAnalysis ? bdtScorePrompt : bdtScoreFd;

// reject entries outside pT ranges of interest
if (ptD < binsPtEfficiencyD->front() || ptD > binsPtEfficiencyD->back()) {
Expand Down Expand Up @@ -344,8 +344,7 @@
double massD = pairEntry.mD();
int effBinD = o2::analysis::findBin(binsPtEfficiencyD, ptD);
int pTBinD = o2::analysis::findBin(binsPtCorrelations, ptD);
float bdtScorePromptorFd = isPromptAnalysis ? bdtScorePrompt: bdtScoreFd;

float bdtScorePromptorFd = isPromptAnalysis ? bdtScorePrompt : bdtScoreFd;

// reject entries outside pT ranges of interest
if (ptD < binsPtEfficiencyD->front() || ptD > binsPtEfficiencyD->back()) {
Expand Down Expand Up @@ -409,8 +408,7 @@
float bdtScoreBkg = candidate.mlScoreBkg();
int effBinD = o2::analysis::findBin(binsPtEfficiencyD, ptD);
bool isDplusPrompt = candidate.isPrompt();
float bdtScorePromptorFd = isPromptAnalysis ? bdtScorePrompt: bdtScoreFd;

float bdtScorePromptorFd = isPromptAnalysis ? bdtScorePrompt : bdtScoreFd;

// reject entries outside pT ranges of interest
if (ptD < binsPtEfficiencyD->front() || ptD > binsPtEfficiencyD->back())
Expand Down Expand Up @@ -453,7 +451,7 @@
float ptHadron = pairEntry.ptHadron();
float massD = pairEntry.mD();
float bdtScorePrompt = pairEntry.mlScorePrompt();
float bdtScoreFd = pairEntry.mlScoreFd();
float bdtScoreFd = pairEntry.mlScoreFd();
float bdtScoreBkg = pairEntry.mlScoreBkg();
bool isPhysicalPrimary = pairEntry.isPhysicalPrimary();
float trackDcaXY = pairEntry.trackDcaXY();
Expand All @@ -464,8 +462,7 @@
int poolBin = pairEntry.poolBin();
int effBinD = o2::analysis::findBin(binsPtEfficiencyD, ptD);
int pTBinD = o2::analysis::findBin(binsPtCorrelations, ptD);
float bdtScorePromptorFd = isPromptAnalysis ? bdtScorePrompt: bdtScoreFd;

float bdtScorePromptorFd = isPromptAnalysis ? bdtScorePrompt : bdtScoreFd;

// reject entries outside pT ranges of interest
if (ptD < binsPtEfficiencyD->front() || ptD > binsPtEfficiencyD->back())
Expand Down Expand Up @@ -512,7 +509,7 @@
registry.fill(HIST("hCorrel2DVsPtPhysicalPrimaryMcRec"), deltaPhi, deltaEta, ptD, ptHadron, statusDplusPrompt, poolBin, efficiencyWeight);
if (statusDplusPrompt == 1 && statusPromptHadron == 1) {
registry.fill(HIST("hCorrel2DVsPtSignalRegionPromptDplusPromptHadronMcRec"), deltaPhi, deltaEta, ptD, ptHadron, poolBin, efficiencyWeight);
} else if (statusDplusPrompt == 0 && statusPromptHadron == 2) {

Check warning on line 512 in PWGHF/HFC/Tasks/taskCorrelationDplusHadrons.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number]

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
registry.fill(HIST("hCorrel2DVsPtSignalRegionNonPromptDplusNonPromptHadronMcRec"), deltaPhi, deltaEta, ptD, ptHadron, poolBin, efficiencyWeight);
}
}
Expand Down Expand Up @@ -564,7 +561,7 @@
}
} else {
registry.fill(HIST("hCorrel2DVsPtMcGenNonPrompt"), deltaPhi, deltaEta, ptD, ptHadron, poolBin);
if (statusPromptHadron == 2) {

Check warning on line 564 in PWGHF/HFC/Tasks/taskCorrelationDplusHadrons.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number]

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
registry.fill(HIST("hCorrel2DVsPtMcGenNonPromptDNonPromptHadron"), deltaPhi, deltaEta, ptD, ptHadron, poolBin);
}
}
Expand Down
Loading