Skip to content

Commit

Permalink
Merge pull request #4 from creetz16/master
Browse files Browse the repository at this point in the history 
merge master into kfstudy
  • Loading branch information
@creetz16
creetz16 authored Dec 11, 2023
2 parents aa14943 + cd9ca62 commit 83e0b74
Show file tree
Hide file tree
Showing 97 changed files with 5,660 additions and 1,596 deletions.
4 changes: 2 additions & 2 deletions ALICE3/TableProducer/OTF/onTheFlyTOFPID.cxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -276,10 +276,10 @@ struct OnTheFlyTOFPID {

if (scalarProduct1 > scalarProduct2) {
modulus = std::hypot(point1[0] - trcCircle.xC, point1[1] - trcCircle.yC) * std::hypot(startPoint[0] - trcCircle.xC, startPoint[1] - trcCircle.yC);
cosAngle = (point1[0] - trcCircle.xC) * (startPoint[0] - trcCircle.xC) + (point1[1] - trcCircle.yC) * (startPoint[0] - trcCircle.yC);
cosAngle = (point1[0] - trcCircle.xC) * (startPoint[0] - trcCircle.xC) + (point1[1] - trcCircle.yC) * (startPoint[1] - trcCircle.yC);
} else {
modulus = std::hypot(point2[0] - trcCircle.xC, point2[1] - trcCircle.yC) * std::hypot(startPoint[0] - trcCircle.xC, startPoint[1] - trcCircle.yC);
cosAngle = (point2[0] - trcCircle.xC) * (startPoint[0] - trcCircle.xC) + (point2[1] - trcCircle.yC) * (startPoint[0] - trcCircle.yC);
cosAngle = (point2[0] - trcCircle.xC) * (startPoint[0] - trcCircle.xC) + (point2[1] - trcCircle.yC) * (startPoint[1] - trcCircle.yC);
}
cosAngle /= modulus;
length = trcCircle.rC * TMath::ACos(cosAngle);
Expand Down
2 changes: 1 addition & 1 deletion CODEOWNERS
View file
Original file line number Diff line number Diff line change
Expand Up @@ -29,7 +29,7 @@
/PWGDQ @alibuild @iarsene @dsekihat @feisenhu @lucamicheletti93
/PWGEM @alibuild @mikesas @rbailhac @feisenhu
/PWGEM/Dilepton @alibuild @mikesas @rbailhac @dsekihat @ivorobye @feisenhu
/PWGEM/PhotonMeson @alibuild @mikesas @rbailhac @m-c-danisch @novitzky @mhemmer-cern
/PWGEM/PhotonMeson @alibuild @mikesas @rbailhac @m-c-danisch @novitzky @mhemmer-cern @dsekihat
/PWGHF @alibuild @vkucera @fcolamar @fgrosa @fcatalan92 @mfaggin @mmazzilli @deepathoms @nzardosh @NicoleBastid
/PWGLF @alibuild @lramona @alcaliva @lbariogl @chiarapinto @BongHwi @lbarnby @mbombara @iravasen @njacazio @ChiaraDeMartin95 @skundu692
/PWGMM @alibuild @aalkin
Expand Down
9 changes: 8 additions & 1 deletion Common/DataModel/PIDResponse.h
View file
Original file line number Diff line number Diff line change
Expand Up @@ -46,6 +46,13 @@ void packInTable(const float& valueToBin, T& table)
}
}

// Function to unpack a binned value into a float
template <typename binningType>
float unPackInTable(const typename binningType::binned_t& valueToUnpack)
{
return binningType::bin_width * static_cast<float>(valueToUnpack);
}

// Checkers for TOF PID hypothesis availability (runtime)
template <class T>
using hasTOFEl = decltype(std::declval<T&>().tofNSigmaEl());
Expand Down Expand Up @@ -481,7 +488,7 @@ DECLARE_SOA_COLUMN(TOFNSigmaAl, tofNSigmaAl, float); //! Nsigma separation with
// Macro to convert the stored Nsigmas to floats
#define DEFINE_UNWRAP_NSIGMA_COLUMN(COLUMN, COLUMN_NAME) \
DECLARE_SOA_DYNAMIC_COLUMN(COLUMN, COLUMN_NAME, \
[](binning::binned_t nsigma_binned) -> float { return binning::bin_width * static_cast<float>(nsigma_binned); });
[](binning::binned_t nsigma_binned) -> float { return o2::aod::pidutils::unPackInTable<binning>(nsigma_binned); });

namespace pidtof_tiny
{
Expand Down
2 changes: 1 addition & 1 deletion Common/TableProducer/trackPropagation.cxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -151,7 +151,7 @@ struct TrackPropagation {

for (auto& track : tracks) {
if constexpr (fillCovMat) {
if (fillTracksDCACov) {
if (fillTracksDCA || fillTracksDCACov) {
mDcaInfoCov.set(999, 999, 999, 999, 999);
}
setTrackParCov(track, mTrackParCov);
Expand Down
102 changes: 90 additions & 12 deletions DPG/Tasks/AOTTrack/perfK0sResolution.cxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -46,6 +46,16 @@ struct perfK0sResolution {
const AxisSpec etaAxis{etaBins, "#eta"};
const AxisSpec phiAxis{phiBins, "#phi"};

int nProc = 0;
if (doprocessData) {
LOG(info) << "processData enabled";
nProc++;
}
if (doprocessMC) {
LOG(info) << "processMC enabled";
nProc++;
}

registry.add("h2_masspT", "h2_masspT", {HistType::kTH2F, {mAxis, pTAxis}});
registry.add("h2_masseta", "h2_masseta", {HistType::kTH2F, {mAxis, etaAxis}});
registry.add("h2_massphi", "h2_massphi", {HistType::kTH2F, {mAxis, phiAxis}});
Expand All @@ -62,26 +72,70 @@ struct perfK0sResolution {
Configurable<float> nSigTPC{"nSigTPC", 10., "nSigTPC"};
Configurable<int> trdSelectionPos{"trdSelectionPos", 0, "Flag for the TRD selection on positive daughters: -1 no TRD, 0 no selection, 1 TRD"};
Configurable<int> trdSelectionNeg{"trdSelectionNeg", 0, "Flag for the TRD selection on negative daughters: -1 no TRD, 0 no selection, 1 TRD"};
Configurable<int> tofSelectionPos{"tofSelectionPos", 0, "Flag for the TOF selection on positive daughters: -1 no TOF, 0 no selection, 1 TOF"};
Configurable<int> tofSelectionNeg{"tofSelectionNeg", 0, "Flag for the TOF selection on negative daughters: -1 no TOF, 0 no selection, 1 TOF"};
Configurable<bool> eventSelection{"eventSelection", true, "event selection"};

template <typename T1, typename T2, typename C>
bool acceptV0(const T1& v0, const T2& ntrack, const T2& ptrack, const C& collision)
{
// Apply selections on V0
if (TMath::Abs(v0.yK0Short()) > rapidity)
if (TMath::Abs(v0.yK0Short()) > rapidity) {
return false;
if (v0.v0cosPA() < v0setting_cospa)
}
if (v0.v0cosPA() < v0setting_cospa) {
return false;
if (v0.v0radius() < v0setting_radius)
}
if (v0.v0radius() < v0setting_radius) {
return false;
if (v0.distovertotmom(collision.posX(), collision.posY(), collision.posZ()) * pid_constants::sMasses[PID::K0] > 2.684 * v0lifetime)
}
if (v0.distovertotmom(collision.posX(), collision.posY(), collision.posZ()) * pid_constants::sMasses[PID::K0] > 2.684 * v0lifetime) {
return false;
}

// Apply selections on V0 daughters
if (!ntrack.hasTPC() || !ptrack.hasTPC())
if (!ntrack.hasTPC() || !ptrack.hasTPC()) {
return false;
if (ntrack.tpcNSigmaPi() > nSigTPC || ptrack.tpcNSigmaPi() > nSigTPC)
}
if (ntrack.tpcNSigmaPi() > nSigTPC || ptrack.tpcNSigmaPi() > nSigTPC) {
return false;
}
// TOF selection
switch (tofSelectionPos) {
case -1:
if (ptrack.hasTOF()) {
return false;
}
break;
case 0:
break;
case 1:
if (!ptrack.hasTOF()) {
return false;
}
break;
default:
LOG(fatal) << "Invalid TRD selection for positive daughter";
break;
}
switch (tofSelectionNeg) {
case -1:
if (ntrack.hasTOF()) {
return false;
}
break;
case 0:
break;
case 1:
if (!ntrack.hasTOF()) {
return false;
}
break;
default:
LOG(fatal) << "Invalid TRD selection for negative daughter";
break;
}
// TRD selection
switch (trdSelectionPos) {
case -1:
if (ptrack.hasTRD()) {
Expand Down Expand Up @@ -121,7 +175,7 @@ struct perfK0sResolution {

Filter v0Filter = nabs(aod::v0data::dcapostopv) > v0setting_dcapostopv&& nabs(aod::v0data::dcanegtopv) > v0setting_dcanegtopv&& aod::v0data::dcaV0daughters < v0setting_dcav0dau;

void process(SelectedCollisions::iterator const& collision, soa::Filtered<aod::V0Datas> const& fullV0s, PIDTracks const& tracks)
void processData(SelectedCollisions::iterator const& collision, soa::Filtered<aod::V0Datas> const& fullV0s, PIDTracks const& tracks)
{
if (eventSelection && !collision.sel8())
return;
Expand All @@ -138,10 +192,34 @@ struct perfK0sResolution {
registry.fill(HIST("h2_massphi"), v0.mK0Short(), v0.phi());
}
}
PROCESS_SWITCH(perfK0sResolution, processData, "Process data", true);

void processMC(SelectedCollisions::iterator const& collision, soa::Filtered<aod::V0Datas> const& fullV0s, soa::Join<PIDTracks, aod::McTrackLabels> const& tracks, aod::McParticles const&)
{
if (eventSelection && !collision.sel8())
return;

for (auto& v0 : fullV0s) {

const auto& posTrack = v0.posTrack_as<soa::Join<PIDTracks, aod::McTrackLabels>>();
const auto& negTrack = v0.negTrack_as<soa::Join<PIDTracks, aod::McTrackLabels>>();
if (!acceptV0(v0, negTrack, posTrack, collision))
continue;
if (!posTrack.has_mcParticle()) {
continue;
}
if (!negTrack.has_mcParticle()) {
continue;
}
if (posTrack.mcParticle().pdgCode() != PID::Pion || negTrack.mcParticle().pdgCode() != PID::Pion) {
continue;
}
registry.fill(HIST("h2_masspT"), v0.mK0Short(), v0.pt());
registry.fill(HIST("h2_masseta"), v0.mK0Short(), v0.eta());
registry.fill(HIST("h2_massphi"), v0.mK0Short(), v0.phi());
}
}
PROCESS_SWITCH(perfK0sResolution, processMC, "Process MC", false);
};

WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)
{
return WorkflowSpec{
adaptAnalysisTask<perfK0sResolution>(cfgc)};
}
WorkflowSpec defineDataProcessing(ConfigContext const& cfgc) { return WorkflowSpec{adaptAnalysisTask<perfK0sResolution>(cfgc)}; }
136 changes: 63 additions & 73 deletions DPG/Tasks/AOTTrack/qaEfficiencyV0s.cxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -36,45 +36,53 @@
#include "THashList.h"

using namespace o2::framework;
static constexpr int nSpecies = 2; // One per PDG
// static constexpr const char* particleTitle[nSpecies] = {"K0s", "-K0s"};
static constexpr int PDGs[nSpecies] = {kK0Short, -kK0Short};
int pdgToIndex(int pdg)
{
for (int i = 0; i < nSpecies; i++) {
if (pdg == PDGs[i]) {
return i;
}
}
return -1;
}
static constexpr int kHistoPtNum = 0;
static constexpr int kHistoPtDen = 1;
static constexpr int kHistoTot = 2;

std::shared_ptr<TH1> histograms[nSpecies][kHistoTot] = {{nullptr}};
std::shared_ptr<TH1> histogramsPrm[nSpecies][kHistoTot] = {{nullptr}};
std::shared_ptr<TH1> histogramsPrmRap[nSpecies][kHistoTot] = {{nullptr}};

struct QaEfficiencyV0s {
// Particle information
static constexpr int nSpecies = 1; // One per PDG
static constexpr const char* particleTitle[nSpecies] = {"K0s"};
static constexpr int PDGs[nSpecies] = {kK0Short};
// Particle only selection
Configurable<bool> doK0s{"do-k0s", false, "Flag to run with the PDG code of k0s"};
Configurable<float> rapidityCut{"rapidityCut", 0.5, "Rapidity cut"};
ConfigurableAxis ptBins{"ptBins", {200, 0.f, 5.f}, "Pt binning"};
OutputObj<THashList> listEfficiencyMC{"EfficiencyMC"};
// Histograms

HistogramRegistry registry{"Histos", {}, OutputObjHandlingPolicy::AnalysisObject};

void init(InitContext&)
{
const AxisSpec axisPt{ptBins, "#it{p}_{T} GeV/#it{c}"};
auto h = registry.add<TH1>("Pos/PtNum_310", "Pos/Num_310", kTH1F, {axisPt});
registry.add("Pos/PtDen_310", "Pos/Den_310", kTH1F, {axisPt});
registry.add("Neg/PtNum_310", "Neg/Num_310", kTH1F, {axisPt});
registry.add("Neg/PtDen_310", "Neg/Den_310", kTH1F, {axisPt});
registry.add("Prm/Pos/PtNum_310", "Prm/Pos/Num_310", kTH1F, {axisPt});
registry.add("Prm/Pos/PtDen_310", "Prm/Pos/Den_310", kTH1F, {axisPt});
registry.add("Prm/Neg/PtNum_310", "Prm/Neg/Num_310", kTH1F, {axisPt});
registry.add("Prm/Neg/PtDen_310", "Prm/Neg/Den_310", kTH1F, {axisPt});
registry.add("PrmRap/Pos/PtNum_310", "PrmRap/Pos/Num_310", kTH1F, {axisPt});
registry.add("PrmRap/Pos/PtDen_310", "PrmRap/Pos/Den_310", kTH1F, {axisPt});
registry.add("PrmRap/Neg/PtNum_310", "PrmRap/Neg/Num_310", kTH1F, {axisPt});
registry.add("PrmRap/Neg/PtDen_310", "PrmRap/Neg/Den_310", kTH1F, {axisPt});

TAxis* axis = h->GetXaxis();
listEfficiencyMC.setObject(new THashList);
listEfficiencyMC->Add(new TEfficiency(Form("efficiencyPt_pdg%d", PDGs[0]), Form("efficiencyPt_pdg%d", PDGs[0]), axis->GetNbins(), axis->GetXmin(), axis->GetXmax()));
listEfficiencyMC->Add(new TEfficiency(Form("efficiencyPt_pdg-%d", PDGs[0]), Form("efficiencyPt_pdg-%d", PDGs[0]), axis->GetNbins(), axis->GetXmin(), axis->GetXmax()));
listEfficiencyMC->Add(new TEfficiency(Form("efficiencyPtPrm_pdg%d", PDGs[0]), Form("efficiencyPtPrm_pdg%d", PDGs[0]), axis->GetNbins(), axis->GetXmin(), axis->GetXmax()));
listEfficiencyMC->Add(new TEfficiency(Form("efficiencyPtPrm_pdg-%d", PDGs[0]), Form("efficiencyPtPrm_pdg-%d", PDGs[0]), axis->GetNbins(), axis->GetXmin(), axis->GetXmax()));
listEfficiencyMC->Add(new TEfficiency(Form("efficiencyPtPrmRap_pdg%d", PDGs[0]), Form("efficiencyPtPrmRap_pdg%d", PDGs[0]), axis->GetNbins(), axis->GetXmin(), axis->GetXmax()));
listEfficiencyMC->Add(new TEfficiency(Form("efficiencyPtPrmRap_pdg-%d", PDGs[0]), Form("efficiencyPtPrmRap_pdg-%d", PDGs[0]), axis->GetNbins(), axis->GetXmin(), axis->GetXmax()));

for (int i = 0; i < nSpecies; i++) {
// MC efficiency (PDG code
histograms[i][kHistoPtNum] = registry.add<TH1>(Form("Pt/Num_%i", PDGs[i]), Form("Num %i", PDGs[i]), kTH1F, {axisPt});
histograms[i][kHistoPtDen] = registry.add<TH1>(Form("Pt/Den_%i", PDGs[i]), Form("Den %i", PDGs[i]), kTH1F, {axisPt});

histogramsPrm[i][kHistoPtNum] = registry.add<TH1>(Form("Pt/Prm/Num_%i", PDGs[i]), Form("Pt Prm Num %i", PDGs[i]), kTH1F, {axisPt});
histogramsPrm[i][kHistoPtDen] = registry.add<TH1>(Form("Pt/Prm/Den_%i", PDGs[i]), Form("Pt Prm Den %i", PDGs[i]), kTH1F, {axisPt});

histogramsPrmRap[i][kHistoPtNum] = registry.add<TH1>(Form("Pt/Prm/Rap/PtNum_%i", PDGs[i]), Form("Pt Prm Rap Num %i", PDGs[i]), kTH1F, {axisPt});
histogramsPrmRap[i][kHistoPtDen] = registry.add<TH1>(Form("Pt/Prm/Rap/PtDen_%i", PDGs[i]), Form("Pt Prm Rap Den %i", PDGs[i]), kTH1F, {axisPt});

TAxis* axis = histograms[i][0]->GetXaxis();
listEfficiencyMC->Add(new TEfficiency(Form("efficiencyPt_pdg%d", PDGs[i]), Form("efficiencyPt_pdg%d", PDGs[i]), axis->GetNbins(), axis->GetXmin(), axis->GetXmax()));
listEfficiencyMC->Add(new TEfficiency(Form("efficiencyPtPrm_pdg%d", PDGs[i]), Form("efficiencyPtPrm_pdg%d", PDGs[i]), axis->GetNbins(), axis->GetXmin(), axis->GetXmax()));
listEfficiencyMC->Add(new TEfficiency(Form("efficiencyPtPrmRap_pdg%d", PDGs[i]), Form("efficiencyPtPrmRap_pdg%d", PDGs[i]), axis->GetNbins(), axis->GetXmin(), axis->GetXmax()));
}
}

// MC process
Expand All @@ -83,62 +91,44 @@ struct QaEfficiencyV0s {
void process(o2::aod::McV0Labels const& V0s,
o2::aod::McParticles const& mcParticles)
{
for (auto const& v0 : V0s) {
for (auto const& v0 : V0s) { // Numerator
if (v0.has_mcParticle()) {
auto mcparticle = v0.mcParticle();
if (mcparticle.pdgCode() == PDGs[0]) {
registry.fill(HIST("Pos/PtNum_310"), mcparticle.pt());
if (mcparticle.isPhysicalPrimary()) {
registry.fill(HIST("Prm/Pos/PtNum_310"), mcparticle.pt());
if (TMath::Abs(mcparticle.y()) < rapidityCut) {
registry.fill(HIST("PrmRap/Pos/PtNum_310"), mcparticle.pt());
}
}
} else if (mcparticle.pdgCode() == -PDGs[0]) {
registry.fill(HIST("Neg/PtNum_310"), mcparticle.pt());
if (mcparticle.isPhysicalPrimary()) {
registry.fill(HIST("Prm/Neg/PtNum_310"), mcparticle.pt());
if (TMath::Abs(mcparticle.y()) < rapidityCut) {
registry.fill(HIST("PrmRap/Neg/PtNum_310"), mcparticle.pt());
}
}
const auto index = pdgToIndex(mcparticle.pdgCode());
if (index < 0) {
continue;
}
}
}
for (auto const& mcparticle : mcParticles) {
if (mcparticle.pdgCode() == PDGs[0]) {
registry.fill(HIST("Pos/PtDen_310"), mcparticle.pt());
histograms[index][kHistoPtNum]->Fill(mcparticle.pt());
if (mcparticle.isPhysicalPrimary()) {
registry.fill(HIST("Prm/Pos/PtDen_310"), mcparticle.pt());
histogramsPrm[index][kHistoPtNum]->Fill(mcparticle.pt());
if (TMath::Abs(mcparticle.y()) < rapidityCut) {
registry.fill(HIST("PrmRap/Pos/PtDen_310"), mcparticle.pt());
histogramsPrmRap[index][kHistoPtNum]->Fill(mcparticle.pt());
}
}
} else if (mcparticle.pdgCode() == -PDGs[0]) {
registry.fill(HIST("Neg/PtDen_310"), mcparticle.pt());
if (mcparticle.isPhysicalPrimary()) {
registry.fill(HIST("Prm/Neg/PtDen_310"), mcparticle.pt());
if (TMath::Abs(mcparticle.y()) < rapidityCut) {
registry.fill(HIST("PrmRap/Neg/PtDen_310"), mcparticle.pt());
}
}
}
for (auto const& mcparticle : mcParticles) { // Denominator
const auto index = pdgToIndex(mcparticle.pdgCode());
if (index < 0) {
continue;
}
histograms[index][kHistoPtDen]->Fill(mcparticle.pt());
if (mcparticle.isPhysicalPrimary()) {
histogramsPrm[index][kHistoPtDen]->Fill(mcparticle.pt());
if (TMath::Abs(mcparticle.y()) < rapidityCut) {
histogramsPrmRap[index][kHistoPtDen]->Fill(mcparticle.pt());
}
}
}

static_cast<TEfficiency*>(listEfficiencyMC->At(0))->SetPassedHistogram(*registry.get<TH1>(HIST("Pos/PtNum_310")), "f");
static_cast<TEfficiency*>(listEfficiencyMC->At(0))->SetTotalHistogram(*registry.get<TH1>(HIST("Pos/PtDen_310")), "f");
static_cast<TEfficiency*>(listEfficiencyMC->At(1))->SetPassedHistogram(*registry.get<TH1>(HIST("Neg/PtNum_310")), "f");
static_cast<TEfficiency*>(listEfficiencyMC->At(1))->SetTotalHistogram(*registry.get<TH1>(HIST("Neg/PtDen_310")), "f");

static_cast<TEfficiency*>(listEfficiencyMC->At(2))->SetPassedHistogram(*registry.get<TH1>(HIST("Prm/Pos/PtNum_310")), "f");
static_cast<TEfficiency*>(listEfficiencyMC->At(2))->SetTotalHistogram(*registry.get<TH1>(HIST("Prm/Pos/PtDen_310")), "f");
static_cast<TEfficiency*>(listEfficiencyMC->At(3))->SetPassedHistogram(*registry.get<TH1>(HIST("Prm/Neg/PtNum_310")), "f");
static_cast<TEfficiency*>(listEfficiencyMC->At(3))->SetTotalHistogram(*registry.get<TH1>(HIST("Prm/Neg/PtDen_310")), "f");

static_cast<TEfficiency*>(listEfficiencyMC->At(4))->SetPassedHistogram(*registry.get<TH1>(HIST("PrmRap/Pos/PtNum_310")), "f");
static_cast<TEfficiency*>(listEfficiencyMC->At(4))->SetTotalHistogram(*registry.get<TH1>(HIST("PrmRap/Pos/PtDen_310")), "f");
static_cast<TEfficiency*>(listEfficiencyMC->At(5))->SetPassedHistogram(*registry.get<TH1>(HIST("PrmRap/Neg/PtNum_310")), "f");
static_cast<TEfficiency*>(listEfficiencyMC->At(5))->SetTotalHistogram(*registry.get<TH1>(HIST("PrmRap/Neg/PtDen_310")), "f");
for (int i = 0; i < nSpecies; i++) {
static_cast<TEfficiency*>(listEfficiencyMC->At(i * nSpecies))->SetPassedHistogram(*histograms[i][kHistoPtNum].get(), "f");
static_cast<TEfficiency*>(listEfficiencyMC->At(i * nSpecies))->SetTotalHistogram(*histograms[i][kHistoPtDen].get(), "f");
static_cast<TEfficiency*>(listEfficiencyMC->At(i * nSpecies + 1))->SetPassedHistogram(*histogramsPrm[i][kHistoPtNum].get(), "f");
static_cast<TEfficiency*>(listEfficiencyMC->At(i * nSpecies + 1))->SetTotalHistogram(*histogramsPrm[i][kHistoPtDen].get(), "f");
static_cast<TEfficiency*>(listEfficiencyMC->At(i * nSpecies + 2))->SetPassedHistogram(*histogramsPrmRap[i][kHistoPtNum].get(), "f");
static_cast<TEfficiency*>(listEfficiencyMC->At(i * nSpecies + 2))->SetTotalHistogram(*histogramsPrmRap[i][kHistoPtDen].get(), "f");
}
}
};

Expand Down
Loading

0 comments on commit 83e0b74

Please sign in to comment.