message: Better QA plot for HM events
diff: @@ -6,6 +6,7 @@ using std::array;
  #include <memory>
  using std::string;
  using std::vector;
  +#include <numeric>
  #include <TClonesArray.h>
  #include <TH1D.h>
  @@ -447,12 +448,25 @@ void AliEventCuts::AddQAplotsToList(TList *qaList, bool addCorrelationPlots) {
  qaList->Add(fNormalisationHist);
  string titles[2] = {"before event cuts","after event cuts"};
  +
  + int nBins{100};
  + std::array<double, 120> centBins;
  + centBins.fill(0);
  + if (fOverrideAutoTriggerMask && fTriggerMask & AliVEvent::kHighMultV0) {
  + nBins = 119;
  + for (int i = 0; i < 20; ++i)
  + centBins[i+1] = centBins[i] + 0.01;
  + std::iota(centBins.begin() + 20, centBins.end(), 1);
  + } else {
  + std::iota(centBins.begin(), centBins.end(), 0);
  + }
  +
  for (int iS = 0; iS < 2; ++iS) {
  fVtz[iS] = new TH1D(Form("Vtz_%s",fkLabels[iS].data()),Form("Vertex z %s; #it{v_{z}} (cm); Events",titles[iS].data()),400,-20.,20.);
  fDeltaTrackSPDvtz[iS] = new TH1D(Form("DeltaVtz_%s",fkLabels[iS].data()),Form("Vertex tracks - Vertex SPD %s; #Delta#it{v_{z}} (cm); Events",titles[iS].data()),400,-2.,2.);
  - fCentrality[iS] = new TH1D(Form("Centrality_%s",fkLabels[iS].data()),Form("Centrality percentile %s; Centrality (%%); Events",titles[iS].data()),100,0.,100.);
  - fEstimCorrelation[iS] = new TH2D(Form("EstimCorrelation_%s",fkLabels[iS].data()),Form("Correlation estimators %s;%s;%s",titles[iS].data(),fCentEstimators[1].data(),fCentEstimators[0].data()),100,0.,100.,100,0.,100.);
  - fMultCentCorrelation[iS] = new TH2D(Form("MultCentCorrelation_%s",fkLabels[iS].data()),Form("Correlation multiplicity-centrality %s;Percentile of %s; Number of tracklets",titles[iS].data(),fCentEstimators[0].data()),100,0.,100.,2000,0.,10000.);
  + fCentrality[iS] = new TH1D(Form("Centrality_%s",fkLabels[iS].data()),Form("Centrality percentile %s; Centrality (%%); Events",titles[iS].data()), nBins, centBins.data());
  + fEstimCorrelation[iS] = new TH2D(Form("EstimCorrelation_%s",fkLabels[iS].data()),Form("Correlation estimators %s;%s;%s",titles[iS].data(),fCentEstimators[1].data(),fCentEstimators[0].data()), nBins, centBins.data(), nBins, centBins.data());
  + fMultCentCorrelation[iS] = new TH2D(Form("MultCentCorrelation_%s",fkLabels[iS].data()),Form("Correlation multiplicity-centrality %s;Percentile of %s; Number of tracklets",titles[iS].data(),fCentEstimators[0].data()), nBins, centBins.data(),2000,0.,10000.);
  qaList->Add(fVtz[iS]);
  qaList->Add(fDeltaTrackSPDvtz[iS]);
  