[PWGCF] new group of QA histos + minor fixes (AliceO2Group#8998)

Author: abilandz

PWGCF/MultiparticleCorrelations/Core/MuPa-Configurables.h

@@ -49,6 +49,9 @@ struct : ConfigurableGroup {
   Configurable<vector<string>> cfBookQAEventHistograms2D{"cfBookQAEventHistograms2D", {"1-Multiplicity_vs_ReferenceMultiplicity", "1-Multiplicity_vs_NContributors", "1-Multiplicity_vs_Centrality", "1-Multiplicity_vs_Vertex_z", "1-Multiplicity_vs_Occupancy", "1-ReferenceMultiplicity_vs_NContributors", "1-ReferenceMultiplicity_vs_Centrality", "1-ReferenceMultiplicity_vs_Vertex_z", "1-ReferenceMultiplicity_vs_Occupancy", "1-NContributors_vs_Centrality", "1-NContributors_vs_Vertex_z", "1-NContributors_vs_Occupancy", "1-Centrality_vs_Vertex_z", "1-Centrality_vs_Occupancy", "0-Centrality_vs_ImpactParameter", "1-Vertex_z_vs_Occupancy", "0-MultNTracksPV_vs_MultNTracksGlobal", "0-CentFT0C_vs_CentNTPV", "0-CentFT0M_vs_CentNTPV", "0-CentRun2V0M_vs_CentRun2SPDTracklets", "1-TrackOccupancyInTimeRange_vs_FT0COccupancyInTimeRange", "1-CurrentRunDuration_vs_InteractionRate"}, "book (1) or do not book (0) this QA 2D event histogram"};
   Configurable<bool> cfFillQAParticleHistograms2D{"cfFillQAParticleHistograms2D", false, "if false, all QA 2D particle histograms are not filled. if true, only the ones for which fBookQAParticleHistograms2D[...] is true, are filled"};
   Configurable<vector<string>> cfBookQAParticleHistograms2D{"cfBookQAParticleHistograms2D", {"1-Pt_vs_dcaXY"}, "book (1) or do not book (0) this QA 2D particle histogram"};
+  Configurable<bool> cfFillQAParticleEventHistograms2D{"cfFillQAParticleEventHistograms2D", false, "if false, all QA 2D particle event histograms are not filled. if true, only the ones for which fBookQAParticleEventHistograms2D[...] is true, are filled"};
+  Configurable<vector<string>> cfBookQAParticleEventHistograms2D{"cfBookQAParticleEventHistograms2D", {"1-CurrentRunDuration_vs_itsNCls", "1-CurrentRunDuration_vs_itsNClsNegEtaEbyE", "1-CurrentRunDuration_vs_itsNClsPosEtaEbyE", "1-CurrentRunDuration_vs_Eta0804EbyE", "1-CurrentRunDuration_vs_Eta0400EbyE", "1-CurrentRunDuration_vs_Eta0004EbyE", "1-CurrentRunDuration_vs_Eta0408EbyE", "1-CurrentRunDuration_vs_Pt0005EbyE", "1-CurrentRunDuration_vs_Pt0510EbyE", "1-CurrentRunDuration_vs_Pt1050EbyE"}, "book (1) or do not book (0) this QA 2D particle event histogram"};
+
 } cf_qa;
 
 // *) Event histograms:

PWGCF/MultiparticleCorrelations/Core/MuPa-DataMembers.h

@@ -84,27 +84,43 @@ struct EventByEventQuantities {
 //    Remark 1: I keep new histograms in this group, until I need them permanently in the analysis. Then, they are moved to EventHistograms or ParticleHistograms (yes, even if they are 2D).
 //    Remark 2: All 2D histograms book as TH2F, due to "stmem error" in terminate (see .cxx for further details)
 struct QualityAssurance {
-  TList* fQAList = NULL;                                                                  //!<! base list to hold all QA output object
-  TProfile* fQAHistogramsPro = NULL;                                                      //!<! keeps flags relevant for the QA histograms
-  Bool_t fCheckUnderflowAndOverflow = kFALSE;                                             // check and bail out if in event and particle histograms there are entries which went to underflow or overflow bins
-  Int_t fRebin = 1;                                                                       // number of bins of selected heavy 2D histograms are devided with this number, there is a configurable cfRebin
-  TH2F* fQAEventHistograms2D[eQAEventHistograms2D_N][2][2] = {{{NULL}}};                  //! [ type - see enum eQAEventHistograms2D ][reco,sim][before, after particle cuts]
-  Bool_t fFillQAEventHistograms2D = kTRUE;                                                // if kFALSE, all 2D event histograms are not filled. if kTRUE, the ones for which fBookQAEventHistograms2D[...] is kTRUE, are filled
-  Bool_t fBookQAEventHistograms2D[eQAEventHistograms2D_N] = {kTRUE};                      // book or not this 2D histogram, see configurable cfBookQAEventHistograms2D
-  Float_t fEventHistogramsBins2D[eQAEventHistograms2D_N][2][3] = {{{0.}}};                // [type - see enum][x,y][nBins,min,max]
-  TString fEventHistogramsName2D[eQAEventHistograms2D_N] = {""};                          // name of fQAEventHistograms2D, determined programatically from other 1D names, to ease bookkeeping
-  TH2F* fQAParticleHistograms2D[eQAParticleHistograms2D_N][2][2] = {{{NULL}}};            //! [ type - see enum eQAParticleHistograms2D ][reco,sim][before, after particle cuts]
-  Bool_t fFillQAParticleHistograms2D = kTRUE;                                             // if kFALSE, all 2D particle histograms are not filled. if kTRUE, the ones for which fBookQAParticleHistograms2D[...] is kTRUE, are filled
-  Bool_t fBookQAParticleHistograms2D[eQAParticleHistograms2D_N] = {kTRUE};                // book or not this 2D histogram, see configurable cfBookQAParticleHistograms2D
-  Float_t fParticleHistogramsBins2D[eQAParticleHistograms2D_N][2][3] = {{{0.}}};          // [type - see enum][x,y][nBins,min,max]
-  TString fParticleHistogramsName2D[eQAParticleHistograms2D_N] = {""};                    // name of fQAParticleHistograms2D, determined programatically from other 1D names, to ease bookkeeping
+  TList* fQAList = NULL;                      //!<! base list to hold all QA output object
+  TProfile* fQAHistogramsPro = NULL;          //!<! keeps flags relevant for the QA histograms
+  Bool_t fCheckUnderflowAndOverflow = kFALSE; // check and bail out if in event and particle histograms there are entries which went to underflow or overflow bins
+  Int_t fRebin = 1;                           // number of bins of selected heavy 2D histograms are devided with this number, there is a configurable cfRebin
+  // ...
+
+  TList* fQAEventList = NULL;                                              //!<! base list to hold all QA event output object
+  TH2F* fQAEventHistograms2D[eQAEventHistograms2D_N][2][2] = {{{NULL}}};   //! [ type - see enum eQAEventHistograms2D ][reco,sim][before, after particle cuts]
+  Bool_t fFillQAEventHistograms2D = kTRUE;                                 // if kFALSE, all 2D event histograms are not filled. if kTRUE, the ones for which fBookQAEventHistograms2D[...] is kTRUE, are filled
+  Bool_t fBookQAEventHistograms2D[eQAEventHistograms2D_N] = {kTRUE};       // book or not this 2D histogram, see configurable cfBookQAEventHistograms2D
+  Float_t fEventHistogramsBins2D[eQAEventHistograms2D_N][2][3] = {{{0.}}}; // [type - see enum][x,y][nBins,min,max]
+  TString fEventHistogramsName2D[eQAEventHistograms2D_N] = {""};           // name of fQAEventHistograms2D, determined programatically from other 1D names, to ease bookkeeping
+
+  TList* fQAParticleList = NULL;                                                 //!<! base list to hold all QA particle output object
+  TH2F* fQAParticleHistograms2D[eQAParticleHistograms2D_N][2][2] = {{{NULL}}};   //! [ type - see enum eQAParticleHistograms2D ][reco,sim][before, after particle cuts]
+  Bool_t fFillQAParticleHistograms2D = kTRUE;                                    // if kFALSE, all 2D histograms in this category are not filled. If kTRUE, the ones for which fBookQAParticleHistograms2D[...] is kTRUE, are filled
+  Bool_t fBookQAParticleHistograms2D[eQAParticleHistograms2D_N] = {kTRUE};       // book or not this 2D histogram, see configurable cfBookQAParticleHistograms2D
+  Float_t fParticleHistogramsBins2D[eQAParticleHistograms2D_N][2][3] = {{{0.}}}; // [type - see enum][x,y][nBins,min,max]
+  TString fParticleHistogramsName2D[eQAParticleHistograms2D_N] = {""};           // name of fQAParticleHistograms2D, determined programatically from other 1D names, to ease bookkeeping
+
+  TList* fQAParticleEventList = NULL;                                                        //!<! base list to hold all QA particle event output object
+  TH2F* fQAParticleEventHistograms2D[eQAParticleEventHistograms2D_N][2][2] = {{{NULL}}};     //! [ type - see enum eQAParticleEventHistograms2D ][reco,sim][before, after cuts]
+  bool fFillQAParticleEventHistograms2D = true;                                              // if false, all 2D histograms in this category are not filled. If true, the ones for which fBookQAParticleEventHistograms2D[...] is true, are filled
+  Bool_t fBookQAParticleEventHistograms2D[eQAParticleEventHistograms2D_N] = {kTRUE};         // book or not this 2D histogram, see configurable cfBookQAParticleEventHistograms2D
+  Float_t fQAParticleEventHistogramsBins2D[eQAParticleEventHistograms2D_N][2][3] = {{{0.}}}; // [type - see enum][x,y][nBins,min,max]
+  TString fQAParticleEventHistogramsName2D[eQAParticleEventHistograms2D_N] = {""};           // name of fQAParticleEventHistograms2D, determined programatically from other 1D names, to ease bookkeeping
+  TProfile* fQAParticleEventProEbyE[2][2] = {{NULL}};                                        // helper profile to calculate <some-particle-property> event-by-event
+                                                                                             // [reco, sim][before, after]. Type dimension is bin.
+
   Float_t fReferenceMultiplicity[eReferenceMultiplicityEstimators_N] = {0.};              // used mostly in QA correlation plots
   TString fReferenceMultiplicityEstimatorName[eReferenceMultiplicityEstimators_N] = {""}; // TBI 20241123 add comment
   Float_t fCentrality[eCentralityEstimators_N] = {0.};                                    // used mostly in QA correlation plots
   TString fCentralityEstimatorName[eCentralityEstimators_N] = {""};                       // TBI 20241123 add comment
   Float_t fOccupancy[eOccupancyEstimators_N] = {0.};                                      // used mostly in QA correlation plots
   TString fOccupancyEstimatorName[eOccupancyEstimators_N] = {""};                         // TBI 20241123 add comment
-} qa;                                                                                     // "qa" is a common label for objects in this struct
+
+} qa; // "qa" is a common label for objects in this struct
 
 // *) Event histograms:
 struct EventHistograms {

PWGCF/MultiparticleCorrelations/Core/MuPa-Enums.h

@@ -278,6 +278,38 @@ enum eQAParticleHistograms2D {
   eQAParticleHistograms2D_N
 };
 
+enum eQAParticleEventHistograms2D {
+  // In this category I do correlation <some-particle-property> vs. some-event-property.
+  // The < ... > goes over all particles in that event.
+  // All < ... > over particles are calculated with helper TProfile
+  // For instance: <nITScls> vs. current run duration
+  eCurrentRunDuration_vs_itsNClsEbyE,
+  eCurrentRunDuration_vs_itsNClsNegEtaEbyE,
+  eCurrentRunDuration_vs_itsNClsPosEtaEbyE,
+  eCurrentRunDuration_vs_Eta0804EbyE,
+  eCurrentRunDuration_vs_Eta0400EbyE,
+  eCurrentRunDuration_vs_Eta0004EbyE,
+  eCurrentRunDuration_vs_Eta0408EbyE,
+  eCurrentRunDuration_vs_Pt0005EbyE,
+  eCurrentRunDuration_vs_Pt0510EbyE,
+  eCurrentRunDuration_vs_Pt1050EbyE,
+  eQAParticleEventHistograms2D_N
+};
+
+enum eQAParticleEventProEbyE {
+  eitsNClsEbyE = 1,   // Labels average <itsNCls> in a given event (therefore "EbyE" is appended). Yes, from one, because it runs over bin content and entries in TProfile for most of the time.
+  eitsNClsNegEtaEbyE, // <itsNCls> in a given event for eta < 0
+  eitsNClsPosEtaEbyE, // <itsNCls> in a given event for eta > 0
+  eEta0804EbyE,       // <eta> in a given event for -0.8 < eta < -0.4
+  eEta0400EbyE,       // <eta> in a given event for -0.4 < eta <  0.0
+  eEta0004EbyE,       // <eta> in a given event for  0.0 < eta <  0.4
+  eEta0408EbyE,       // <eta> in a given event for  0.4 < eta <  0.8
+  ePt0005EbyE,        // <pt> in a given event for  0.0 < pt < 0.5
+  ePt0510EbyE,        // <pt> in a given event for  0.5 < pt < 1.0
+  ePt1050EbyE,        // <pt> in a given event for  1.0 < pt < 5.0
+  eQAParticleEventProEbyE_N
+};
+
 enum eReferenceMultiplicityEstimators {
   // Run 3:
   eMultTPC = 0,