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
34 changes: 16 additions & 18 deletions PWGCF/TwoParticleCorrelations/Tasks/neutronProtonCorrZdc.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -73,12 +73,12 @@
const AxisSpec axisZPSignal{2 * cfgNBinsZP, cfgZPmin, 1.5 * cfgZPmax, "ZP (a.u.)"};
const AxisSpec axisAlphaZ{cfgNBinsAlpha, cfgAlphaZmin, cfgAlphaZmax, "#alpha_{spec}"};
const AxisSpec axisZDiffSignal{cfgNBinsZN, cfgDiffZmin, cfgDiffZmax, "#Delta E"};
const AxisSpec AxisMultiplicityF0A{cfgNBinsMultiplicity, 0, 200000, "F0A"};

Check warning on line 76 in PWGCF/TwoParticleCorrelations/Tasks/neutronProtonCorrZdc.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/function-variable]

Use lowerCamelCase for names of functions and variables.
const AxisSpec AxisMultiplicityF0C{cfgNBinsMultiplicity, 0, 100000, "F0C"};

Check warning on line 77 in PWGCF/TwoParticleCorrelations/Tasks/neutronProtonCorrZdc.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/function-variable]

Use lowerCamelCase for names of functions and variables.
const AxisSpec AxisMultiplicityF0M{cfgNBinsMultiplicity, 0, 300000, "F0M"};

Check warning on line 78 in PWGCF/TwoParticleCorrelations/Tasks/neutronProtonCorrZdc.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/function-variable]

Use lowerCamelCase for names of functions and variables.
const AxisSpec AxisMultiplicityFDD{cfgNBinsMultiplicity, 0, 50000, "FDD"};

Check warning on line 79 in PWGCF/TwoParticleCorrelations/Tasks/neutronProtonCorrZdc.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/function-variable]

Use lowerCamelCase for names of functions and variables.
const AxisSpec AxisMultiplicityTPC{cfgNBinsMultiplicity, 0, 100000, "TPC"};

Check warning on line 80 in PWGCF/TwoParticleCorrelations/Tasks/neutronProtonCorrZdc.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/function-variable]

Use lowerCamelCase for names of functions and variables.
const AxisSpec AxisMultiplicityMultNGlobal{cfgNBinsMultiplicity, 0, 3500, "MultsNGlobal"};

Check warning on line 81 in PWGCF/TwoParticleCorrelations/Tasks/neutronProtonCorrZdc.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/function-variable]

Use lowerCamelCase for names of functions and variables.

HistogramConfigSpec defaultZNSectorHist({HistType::kTH2F, {cfgAxisCent, axisZNSectorSignal}});
HistogramConfigSpec defaultZPSectorHist({HistType::kTH2F, {cfgAxisCent, axisZPSectorSignal}});
Expand Down Expand Up @@ -134,24 +134,23 @@
histos.add("MultiplicityHistograms/CentvsFT0M", "CentvsFT0M", kTH2F, {cfgAxisCent, AxisMultiplicityF0M});
histos.add("MultiplicityHistograms/CentvsFV0A", "CentvsFV0A", kTH2F, {cfgAxisCent, AxisMultiplicityF0A});
histos.add("MultiplicityHistograms/CentvsNGlobal", "CentvsNGlobal", kTH2F, {cfgAxisCent, AxisMultiplicityMultNGlobal});

}

template<int mult, typename C>
template <int mult, typename C>
void fillMultHistosRun3(const C& col)
{
static constexpr std::string_view multLabels[] = {"FT0C", "FT0A", "FV0A", "FDDC", "FDDA", "TPC", "NGlobal"};

Check warning on line 142 in PWGCF/TwoParticleCorrelations/Tasks/neutronProtonCorrZdc.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/constexpr-constant]

Use UpperCamelCase for names of constexpr constants. Names of special constants may be prefixed with "k".
std::array<float, 7> multarray = {col.multFT0C(), col.multFT0A(), col.multFV0A(), col.multFDDC(), col.multFDDA(), float(col.multTPC()), float(col.multNTracksGlobal())};

histos.fill(HIST("MultiplicityHistograms/") + HIST(multLabels[mult]), multarray[mult]);
}

template<int cent, typename C>
template <int cent, typename C>
void fillCentHistosRun3(const C& col)
{
static constexpr std::string_view centLabels[] = {"CentvsFT0C", "CentvsFT0CVar1", "CentvsFT0M", "CentvsFV0A", "CentvsNGlobal"};

Check warning on line 151 in PWGCF/TwoParticleCorrelations/Tasks/neutronProtonCorrZdc.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/constexpr-constant]

Use UpperCamelCase for names of constexpr constants. Names of special constants may be prefixed with "k".
std::array<float, 5> centarray = {col.centFT0C(), col.centFT0CVariant1(), col.centFT0M(), col.centFV0A(), col.centNGlobal()};
std::array<float, 5> multarray = {col.multFT0C(), col.multFT0C(), col.multFT0C()+col.multFT0A(), col.multFV0A(), float(col.multNTracksGlobal())};
std::array<float, 5> multarray = {col.multFT0C(), col.multFT0C(), col.multFT0C() + col.multFT0A(), col.multFV0A(), float(col.multNTracksGlobal())};

histos.fill(HIST("MultiplicityHistograms/") + HIST(centLabels[cent]), centarray[cent], multarray[cent]);
}
Expand Down Expand Up @@ -186,10 +185,9 @@

std::array<std::array<float, 4>, 2> znEnergyResponse = {zdc.energySectorZNA(), zdc.energySectorZNC()};
std::array<std::array<float, 4>, 2> zpEnergyResponse = {zdc.energySectorZPA(), zdc.energySectorZPC()};

histos.fill(HIST(SubDir[side]) + HIST(ZNSector[sector]), centr, znEnergyResponse[side][sector]);
histos.fill(HIST(SubDir[side]) + HIST(ZPSector[sector]), centr, zpEnergyResponse[side][sector]);

}

void processRun3(soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::Mults, aod::MultsGlobal, CentralitiesRun3>>::iterator const& collision, BCsRun3 const&, aod::Zdcs const&)
Expand All @@ -213,18 +211,18 @@
histos.fill(HIST("eventCounter"), EventCounter::kZDCSelection);
histos.fill(HIST("CentralityPercentile"), cent);

static_for<0,6>([&](auto i){
static_for<0, 6>([&](auto i) {
fillMultHistosRun3<i>(collision); // Fill multiplicity histograms
});
static_for<0,4>([&](auto i){

static_for<0, 4>([&](auto i) {
fillCentHistosRun3<i>(collision); // Fill centrality histograms
});
});

static_for<0,1>([&](auto i){
static_for<0, 1>([&](auto i) {
fillZDCSideCommonHistos<i>(cent, zdcread); // Fill i-side common histograms
static_for<0,3>([&](auto j){
fillZDCSideSectorHistos<i,j>(cent, zdcread); // Fill i-side sector j histograms
static_for<0, 3>([&](auto j) {
fillZDCSideSectorHistos<i, j>(cent, zdcread); // Fill i-side sector j histograms
});
});

Expand All @@ -251,7 +249,7 @@
histos.fill(HIST("CentvsZPAvsZNC"), cent, sumZPA, sumZNC);
histos.fill(HIST("CentvsZPAvsZPC"), cent, sumZPA, sumZPC);
histos.fill(HIST("CentvsZNCvsZPC"), cent, sumZNC, sumZPC);
histos.fill(HIST("CentvsZNvsZP"), cent, sumZNA+sumZNC, sumZPA+sumZPC);
histos.fill(HIST("CentvsZNvsZP"), cent, sumZNA + sumZNC, sumZPA + sumZPC);
}
}
PROCESS_SWITCH(NeutronProtonCorrZdc, processRun3, "Process analysis for Run 3 data", true);
Expand All @@ -275,10 +273,10 @@
histos.fill(HIST("eventCounter"), EventCounter::kZDCSelection);
histos.fill(HIST("CentralityPercentile"), cent);

static_for<0,1>([&](auto i){
static_for<0, 1>([&](auto i) {
fillZDCSideCommonHistos<i>(cent, zdcread); // Fill i-side common channels
static_for<0,3>([&](auto j){
fillZDCSideSectorHistos<i,j>(cent, zdcread); // Fill i-side sector j
static_for<0, 3>([&](auto j) {
fillZDCSideSectorHistos<i, j>(cent, zdcread); // Fill i-side sector j
});
});

Expand All @@ -305,7 +303,7 @@
histos.fill(HIST("CentvsZPAvsZNC"), cent, sumZPA, sumZNC);
histos.fill(HIST("CentvsZPAvsZPC"), cent, sumZPA, sumZPC);
histos.fill(HIST("CentvsZNCvsZPC"), cent, sumZNC, sumZPC);
histos.fill(HIST("CentvsZNvsZP"), cent, sumZNA+sumZNC, sumZPA+sumZPC);
histos.fill(HIST("CentvsZNvsZP"), cent, sumZNA + sumZNC, sumZPA + sumZPC);
}
}
PROCESS_SWITCH(NeutronProtonCorrZdc, processRun2, "Process analysis for Run 2 converted data", false);
Expand Down
Loading