TSSAplotter/src/TSSAplotter.h

0001 #include <utility>
0002 #include <string>
0003
0004 #include <TH1.h>
0005 #include <TGraphErrors.h>
0006 #include <TGraphAsymmErrors.h>
0007 #include <TCanvas.h>
0008 #include <TPad.h>
0009 #include <TString.h>
0010
0011 #include "binnedhistset/BinnedHistSet.h"
0012
0013 class TSSAplotter {
0014     public:
0015     TSSAplotter();
0016     ~TSSAplotter();
0017
0018     void GetNMHists();
0019     void GetTSSAHists();
0020     void PlotNMHists();
0021     void PlotTSSAs();
0022
0023     double RelLumiAsym(double Nup, double Ndown, double relLumi);
0024     double RelLumiError(double Nup, double Ndown, double relLumi);
0025     double SqrtAsym(double NLup, double NLdown, double NRup, double NRdown);
0026     double SqrtError(double NLup, double NLdown, double NRup, double NRdown,
0027         double NLupErr, double NLdownErr, double NRupErr, double NRdownErr);
0028     TGraphErrors* RelLumiGraph(TH1* phi_up, TH1* phi_down, double relLumi);
0029     TGraphErrors* SqrtGraph(TH1* phi_up, TH1* phi_down);
0030     std::pair<float, float> RawAsym(std::string type, TGraphErrors* gr);
0031     void PlotOneAsym(std::string type, std::string title, TPad* pad, TGraphErrors* gr);
0032     void CompareOneAsym(std::string title, std::string outfilename, TGraphErrors* rl, TGraphErrors* sqrt);
0033     void PlotAsymsBinned(std::string type, std::string which, BinnedHistSet* bhs_up, BinnedHistSet* bhs_down, std::string outfilename);
0034     void PlotAsymsBinned(std::string which, BinnedHistSet* bhs_up, BinnedHistSet* bhs_down, std::string outfilename);
0035     void PlotAsymsBinnedBlueYellow(std::string type, std::string which, BinnedHistSet* bhs_up_blue, BinnedHistSet* bhs_down_blue, BinnedHistSet* bhs_up_yellow, BinnedHistSet* bhs_down_yellow, std::string outfilename);
0036     std::pair<double, double> FitGraph(std::string type, TGraphErrors* gr);
0037     TGraphErrors* AmplitudeBinnedGraph(std::string type, std::string title, std::string xlabel, BinnedHistSet* bhs_up, BinnedHistSet* bhs_down);
0038     void PlotAmplitudes(std::string type, std::string title, std::string xlabel, BinnedHistSet* bhs_up, BinnedHistSet* bhs_down, std::string outfilename);
0039     void CompareAmplitudes(std::string title, std::string xlabel, BinnedHistSet* bhs_up, BinnedHistSet* bhs_down, std::string outfilename);
0040     void CompareAmplitudesEtaBinned(std::string type, std::string title, std::string xlabel, BinnedHistSet* bhs_up_lowEta, BinnedHistSet* bhs_down_lowEta, BinnedHistSet* bhs_up_highEta, BinnedHistSet* bhs_down_highEta, std::string outfilename);
0041     void CompareAmplitudesBlueYellow(std::string type, std::string title, std::string xlabel, BinnedHistSet* bhs_up_blue, BinnedHistSet* bhs_down_blue, BinnedHistSet* bhs_up_yellow, BinnedHistSet* bhs_down_yellow, std::string outfilename);
0042     std::pair<double, double> GetCorrectedAsymErr(double A_sig, double A_sig_err, double A_bkgr, double A_bkgr_err, double r, double r_err, double pol);
0043     TGraphErrors* CorrectedAmplitudesGraph(std::string type, std::string title, std::string xlabel, BinnedHistSet* bhs_up_sig, BinnedHistSet* bhs_down_sig, BinnedHistSet* bhs_up_bkgr, BinnedHistSet* bhs_down_bkgr);
0044     void PlotCorrectedAmplitudes(std::string type, std::string title, std::string xlabel, BinnedHistSet* bhs_up_sig, BinnedHistSet* bhs_down_sig, BinnedHistSet* bhs_up_bkgr, BinnedHistSet* bhs_down_bkgr, std::string outfilename);
0045     TGraphErrors* AveragedGraphs(TGraphErrors* gr_blue, TGraphErrors* gr_yellow);
0046     void PlotAveragedAmplitudes(std::string type, std::string title, std::string xlabel, std::string outfilename, BinnedHistSet* bhs_up_sig_blue, BinnedHistSet* bhs_down_sig_blue, BinnedHistSet* bhs_up_bkgr_blue, BinnedHistSet* bhs_down_bkgr_blue, BinnedHistSet* bhs_up_sig_yellow, BinnedHistSet* bhs_down_sig_yellow, BinnedHistSet* bhs_up_bkgr_yellow, BinnedHistSet* bhs_down_bkgr_yellow);
0047     void PlotGraph(TGraphErrors* gr, std::string title, std::string xlabel, std::string ylabel, std::string outfilename, bool inset = false, bool ylines = false);
0048     void CompareGraphs(TGraphErrors* gr1, TGraphErrors* gr2, std::string title, std::string xlabel, std::string ylabel, std::string leglabel1, std::string leglabel2, std::string outfilename, bool inset = false);
0049     void CompareGraphs(TGraphErrors* gr1, TGraphAsymmErrors* gr2, std::string title, std::string xlabel, std::string ylabel, std::string leglabel1, std::string leglabel2, std::string outfilename, bool inset = false);
0050     TGraphErrors* PolCorrectedGraph(TGraphErrors* gr, double pol);
0051     void PlotBackgroundFractions(std::string outfilename);
0052     TGraphErrors* tTestGraph(std::string type, TGraphErrors* gr1, TGraphErrors* gr2);
0053     TGraphErrors* tTestZeroGraph(TGraphErrors* gr1);
0054     void PlottTest(TGraphErrors* gr, std::string title, std::string xlabel, std::string outfilename, bool limit_yrange = true);
0055     void PlotPi0(std::string outfilename);
0056     void PlotEta(std::string outfilename);
0057     int main(std::string TSSApdf = "TSSAplots.pdf", std::string NMpdf = "");
0058
0059     bool do_NMhists = false;
0060     std::string NMhistfname = "";
0061     std::string TSSAhistfname = "";
0062     std::string NMpdfname = "NMplots_Aug26.pdf";
0063     std::string TSSApdfname = "TSSAplots_Aug26.pdf";
0064
0065     private:
0066     const float PI = 3.141592;
0067     double relLumiBlue = 0.99822;
0068     double relLumiYellow = 0.99767;
0069     double polBlue = 0.48644;
0070     double polYellow = 0.45455;
0071
0072     // NM hists
0073     TH1* h_nEvents = nullptr; // from NM
0074     TH1* h_vtxz = nullptr;
0075     TH2* h_towerEta_Phi_500MeV = nullptr;
0076     TH2* h_towerEta_Phi_800MeV = nullptr;
0077     // cluster distributions
0078     TH1* h_nAllClusters = nullptr;
0079     TH1* h_nGoodClusters = nullptr;
0080     TH1* h_clusterE = nullptr;
0081     TH1* h_clusterEta = nullptr;
0082     TH2* h_clusterVtxz_Eta = nullptr;
0083     TH1* h_clusterPhi = nullptr;
0084     TH2* h_clusterEta_Phi = nullptr;
0085     TH1* h_clusterpT = nullptr;
0086     TH1* h_clusterxF = nullptr;
0087     TH2* h_clusterpT_xF = nullptr;
0088     TH1* h_clusterChi2 = nullptr;
0089     TH1* h_clusterChi2zoomed = nullptr;
0090     TH1* h_mesonClusterChi2 = nullptr;
0091     TH1* h_goodClusterE = nullptr;
0092     TH1* h_goodClusterEta = nullptr;
0093     TH2* h_goodClusterVtxz_Eta = nullptr;
0094     TH1* h_goodClusterPhi = nullptr;
0095     TH2* h_goodClusterEta_Phi = nullptr;
0096     TH1* h_goodClusterpT = nullptr;
0097     TH1* h_goodClusterxF = nullptr;
0098     TH2* h_goodClusterpT_xF = nullptr;
0099     // diphoton distributions
0100     TH1* h_nDiphotons = nullptr;
0101     TH1* h_nRecoPi0s = nullptr;
0102     TH1* h_nRecoEtas = nullptr;
0103     TH1* h_diphotonMass = nullptr;
0104     TH1* h_diphotonE = nullptr;
0105     TH1* h_diphotonEta = nullptr;
0106     TH2* h_diphotonVtxz_Eta = nullptr;
0107     TH1* h_diphotonPhi = nullptr;
0108     TH2* h_diphotonEta_Phi = nullptr;
0109     TH1* h_diphotonpT = nullptr;
0110     TH1* h_diphotonxF = nullptr;
0111     TH2* h_diphotonpT_xF = nullptr;
0112     TH1* h_diphotonAsym = nullptr;
0113     TH1* h_diphotonDeltaR = nullptr;
0114     TH2* h_diphotonAsym_DeltaR = nullptr;
0115     BinnedHistSet* bhs_diphotonMass_pT = nullptr;
0116     BinnedHistSet* bhs_diphotonMass_xF = nullptr;
0117
0118     // TSSA hists
0119     TH1* h_nClustersMinbiasTrig = nullptr;
0120     TH1* h_nClustersPhotonTrig = nullptr;
0121     TH1* h_Nevents = nullptr;
0122     TH1* h_Npi0s = nullptr;
0123     TH1* h_Netas = nullptr;
0124
0125     BinnedHistSet* bhs_diphotonMass_pT_pi0binning = nullptr;
0126     BinnedHistSet* bhs_diphotonMass_pT_etabinning = nullptr;
0127     BinnedHistSet* bhs_pi0_pT_pT = nullptr;
0128     BinnedHistSet* bhs_eta_pT_pT = nullptr;
0129     BinnedHistSet* bhs_diphotonxF_xF = nullptr;
0130     BinnedHistSet* bhs_diphotoneta_eta = nullptr;
0131     BinnedHistSet* bhs_diphotonvtxz_vtxz = nullptr;
0132     TH2* h_DiphotonMassAsym = nullptr;
0133     TH2* h_DiphotonMassdeltaR = nullptr;
0134     TH2* h_DiphotondeltaRAsym = nullptr;
0135     TH2* h_DiphotondeltaRAsym_pi0 = nullptr;
0136     TH2* h_DiphotondeltaRAsym_pi0_smalldR = nullptr;
0137     TH2* h_DiphotondeltaRAsym_eta = nullptr;
0138     TH2* h_DiphotondeltaRAsym_etabkgr = nullptr;
0139     TH2* h_DiphotonMassAsym_highpT = nullptr;
0140     TH2* h_DiphotonMassdeltaR_highpT = nullptr;
0141     TH2* h_DiphotondeltaRAsym_highpT = nullptr;
0142     TH2* h_DiphotondeltaRAsym_highpT_smalldR = nullptr;
0143     TH2* h_DiphotondeltaRAsym_eta_highpT = nullptr;
0144     TH2* h_DiphotondeltaRAsym_etabkgr_highpT = nullptr;
0145     // Armenteros-Podolanski plots
0146     TH2* h_armen_all = nullptr;
0147     TH2* h_armen_pi0 = nullptr;
0148     TH2* h_armen_pi0bkgr = nullptr;
0149     TH2* h_armen_eta = nullptr;
0150     TH2* h_armen_etabkgr = nullptr;
0151     TH2* h_armen_eta_highpT = nullptr;
0152     TH2* h_armen_etabkgr_highpT = nullptr;
0153     TH1* h_armen_p_L = nullptr;
0154     TH1* h_armen_p_T = nullptr;
0155     TH2* h_armen_alpha_alphaE = nullptr;
0156     TH2* h_armen_alpha_deltaR = nullptr;
0157     TH2* h_armen_alpha_deltaR_pi0 = nullptr;
0158     TH2* h_armen_alpha_deltaR_eta = nullptr;
0159     TH2* h_armen_p_T_deltaR = nullptr;
0160     TH2* h_armen_pT_p_L = nullptr;
0161     TH2* h_armen_pT_p_T = nullptr;
0162
0163     // Phi hists
0164     /* std::vector<double> pTbins = {2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 10.0, 20.0}; */
0165     std::vector<double> pTbins_pi0 = {1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 10.0};
0166     std::vector<double> pTbins_eta = {2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 10.0, 20.0};
0167     std::vector<double> xFbins = {-0.15, -0.10, -0.06, -0.03, 0.0, 0.03, 0.06, 0.10, 0.15};
0168     std::vector<double> etabins = {-3.0, -1.75, -1.15, -0.35, 0.0, 0.35, 1.15, 1.75, 3.0};
0169     std::vector<double> vtxzbins = {200.0, -100.0, -50.0, -30.0, 0.0, 30.0, 50.0, 100.0, 200.0};
0170     std::vector<double> pTmeans_pi0 = {1.63, 2.483, 3.436, 4.41, 5.395, 6.391, 7.39, 8.588};
0171     std::vector<double> pTmeans_eta = {2.432, 3.428, 4.422, 5.414, 6.409, 7.415, 8.717, 11.54};
0172     std::vector<double> xFmeans = {-0.112, -0.0687, -0.0394, -0.0153, 0.0141, 0.0400, 0.0687, 0.112};
0173     std::vector<double> etameans = {-1.915, -1.368, -0.7822, -0.1776, 0.1731, 0.7643, 1.361, 1.923};
0174     std::vector<double> vtxzmeans = {-120.6, -68.79, -39.3, -13.03, 13.44, 39.43, 68.59, 121.4};
0175     std::vector<double> pi0_bkgr_fractions = {0.6831, 0.47138, 0.22116, 0.17628, 0.17482, 0.19737, 0.24736, 0.7249};
0176     std::vector<double> pi0_bkgr_fraction_errs = {0.00018, 0.00014, 0.00017, 0.00026, 0.00047, 0.00091, 0.00188, 0.00184};
0177     std::vector<double> eta_bkgr_fractions = {0.96982, 0.89879, 0.83979, 0.79816, 0.76808, 0.73977, 0.74278, 0.76591};
0178     std::vector<double> eta_bkgr_fraction_errs = {0.00005, 0.00015, 0.00031, 0.00058, 0.00109, 0.00202, 0.00275, 0.00500};
0179     /* std::vector<double> xFbins = {-0.15, -0.10, -0.06, -0.03, 0.0, 0.03, 0.06, 0.10, 0.15}; */
0180     /* std::vector<double> etabins = {-2.0, -1.15, -0.35, 0.0, 0.35, 1.15, 2.0}; */
0181     /* std::vector<double> vtxzbins = {-100.0, -50.0, -30.0, 0.0, 30.0, 50.0, 100.0}; */
0182     int nHistBins_phi = 12;
0183
0184     BinnedHistSet* bhs_pi0_blue_up_phi_pT = nullptr;
0185     BinnedHistSet* bhs_pi0_blue_down_phi_pT = nullptr;
0186     BinnedHistSet* bhs_pi0_yellow_up_phi_pT = nullptr;
0187     BinnedHistSet* bhs_pi0_yellow_down_phi_pT = nullptr;
0188     BinnedHistSet* bhs_pi0_blue_up_forward_phi_pT = nullptr;
0189     BinnedHistSet* bhs_pi0_blue_down_forward_phi_pT = nullptr;
0190     BinnedHistSet* bhs_pi0_yellow_up_forward_phi_pT = nullptr;
0191     BinnedHistSet* bhs_pi0_yellow_down_forward_phi_pT = nullptr;
0192     BinnedHistSet* bhs_pi0_blue_up_backward_phi_pT = nullptr;
0193     BinnedHistSet* bhs_pi0_blue_down_backward_phi_pT = nullptr;
0194     BinnedHistSet* bhs_pi0_yellow_up_backward_phi_pT = nullptr;
0195     BinnedHistSet* bhs_pi0_yellow_down_backward_phi_pT = nullptr;
0196     BinnedHistSet* bhs_pi0_lowEta_blue_up_forward_phi_pT = nullptr;
0197     BinnedHistSet* bhs_pi0_lowEta_blue_down_forward_phi_pT = nullptr;
0198     BinnedHistSet* bhs_pi0_lowEta_yellow_up_forward_phi_pT = nullptr;
0199     BinnedHistSet* bhs_pi0_lowEta_yellow_down_forward_phi_pT = nullptr;
0200     BinnedHistSet* bhs_pi0_highEta_blue_up_forward_phi_pT = nullptr;
0201     BinnedHistSet* bhs_pi0_highEta_blue_down_forward_phi_pT = nullptr;
0202     BinnedHistSet* bhs_pi0_highEta_yellow_up_forward_phi_pT = nullptr;
0203     BinnedHistSet* bhs_pi0_highEta_yellow_down_forward_phi_pT = nullptr;
0204     BinnedHistSet* bhs_pi0_blue_up_phi_xF = nullptr;
0205     BinnedHistSet* bhs_pi0_blue_down_phi_xF = nullptr;
0206     BinnedHistSet* bhs_pi0_yellow_up_phi_xF = nullptr;
0207     BinnedHistSet* bhs_pi0_yellow_down_phi_xF = nullptr;
0208     BinnedHistSet* bhs_pi0_blue_up_phi_eta = nullptr;
0209     BinnedHistSet* bhs_pi0_blue_down_phi_eta = nullptr;
0210     BinnedHistSet* bhs_pi0_yellow_up_phi_eta = nullptr;
0211     BinnedHistSet* bhs_pi0_yellow_down_phi_eta = nullptr;
0212     BinnedHistSet* bhs_pi0_blue_up_phi_vtxz = nullptr;
0213     BinnedHistSet* bhs_pi0_blue_down_phi_vtxz = nullptr;
0214     BinnedHistSet* bhs_pi0_yellow_up_phi_vtxz = nullptr;
0215     BinnedHistSet* bhs_pi0_yellow_down_phi_vtxz = nullptr;
0216
0217     BinnedHistSet* bhs_eta_blue_up_phi_pT = nullptr;
0218     BinnedHistSet* bhs_eta_blue_down_phi_pT = nullptr;
0219     BinnedHistSet* bhs_eta_yellow_up_phi_pT = nullptr;
0220     BinnedHistSet* bhs_eta_yellow_down_phi_pT = nullptr;
0221     BinnedHistSet* bhs_eta_blue_up_forward_phi_pT = nullptr;
0222     BinnedHistSet* bhs_eta_blue_down_forward_phi_pT = nullptr;
0223     BinnedHistSet* bhs_eta_yellow_up_forward_phi_pT = nullptr;
0224     BinnedHistSet* bhs_eta_yellow_down_forward_phi_pT = nullptr;
0225     BinnedHistSet* bhs_eta_lowEta_blue_up_forward_phi_pT = nullptr;
0226     BinnedHistSet* bhs_eta_lowEta_blue_down_forward_phi_pT = nullptr;
0227     BinnedHistSet* bhs_eta_lowEta_yellow_up_forward_phi_pT = nullptr;
0228     BinnedHistSet* bhs_eta_lowEta_yellow_down_forward_phi_pT = nullptr;
0229     BinnedHistSet* bhs_eta_highEta_blue_up_forward_phi_pT = nullptr;
0230     BinnedHistSet* bhs_eta_highEta_blue_down_forward_phi_pT = nullptr;
0231     BinnedHistSet* bhs_eta_highEta_yellow_up_forward_phi_pT = nullptr;
0232     BinnedHistSet* bhs_eta_highEta_yellow_down_forward_phi_pT = nullptr;
0233     BinnedHistSet* bhs_eta_blue_up_backward_phi_pT = nullptr;
0234     BinnedHistSet* bhs_eta_blue_down_backward_phi_pT = nullptr;
0235     BinnedHistSet* bhs_eta_yellow_up_backward_phi_pT = nullptr;
0236     BinnedHistSet* bhs_eta_yellow_down_backward_phi_pT = nullptr;
0237     BinnedHistSet* bhs_eta_blue_up_phi_xF = nullptr;
0238     BinnedHistSet* bhs_eta_blue_down_phi_xF = nullptr;
0239     BinnedHistSet* bhs_eta_yellow_up_phi_xF = nullptr;
0240     BinnedHistSet* bhs_eta_yellow_down_phi_xF = nullptr;
0241     BinnedHistSet* bhs_eta_blue_up_phi_eta = nullptr;
0242     BinnedHistSet* bhs_eta_blue_down_phi_eta = nullptr;
0243     BinnedHistSet* bhs_eta_yellow_up_phi_eta = nullptr;
0244     BinnedHistSet* bhs_eta_yellow_down_phi_eta = nullptr;
0245     BinnedHistSet* bhs_eta_blue_up_phi_vtxz = nullptr;
0246     BinnedHistSet* bhs_eta_blue_down_phi_vtxz = nullptr;
0247     BinnedHistSet* bhs_eta_yellow_up_phi_vtxz = nullptr;
0248     BinnedHistSet* bhs_eta_yellow_down_phi_vtxz = nullptr;
0249
0250     BinnedHistSet* bhs_pi0bkgr_blue_up_phi_pT = nullptr;
0251     BinnedHistSet* bhs_pi0bkgr_blue_down_phi_pT = nullptr;
0252     BinnedHistSet* bhs_pi0bkgr_yellow_up_phi_pT = nullptr;
0253     BinnedHistSet* bhs_pi0bkgr_yellow_down_phi_pT = nullptr;
0254     BinnedHistSet* bhs_pi0bkgr_blue_up_forward_phi_pT = nullptr;
0255     BinnedHistSet* bhs_pi0bkgr_blue_down_forward_phi_pT = nullptr;
0256     BinnedHistSet* bhs_pi0bkgr_yellow_up_forward_phi_pT = nullptr;
0257     BinnedHistSet* bhs_pi0bkgr_yellow_down_forward_phi_pT = nullptr;
0258     BinnedHistSet* bhs_pi0bkgr_blue_up_backward_phi_pT = nullptr;
0259     BinnedHistSet* bhs_pi0bkgr_blue_down_backward_phi_pT = nullptr;
0260     BinnedHistSet* bhs_pi0bkgr_yellow_up_backward_phi_pT = nullptr;
0261     BinnedHistSet* bhs_pi0bkgr_yellow_down_backward_phi_pT = nullptr;
0262     BinnedHistSet* bhs_pi0bkgr_lowEta_blue_up_forward_phi_pT = nullptr;
0263     BinnedHistSet* bhs_pi0bkgr_lowEta_blue_down_forward_phi_pT = nullptr;
0264     BinnedHistSet* bhs_pi0bkgr_lowEta_yellow_up_forward_phi_pT = nullptr;
0265     BinnedHistSet* bhs_pi0bkgr_lowEta_yellow_down_forward_phi_pT = nullptr;
0266     BinnedHistSet* bhs_pi0bkgr_highEta_blue_up_forward_phi_pT = nullptr;
0267     BinnedHistSet* bhs_pi0bkgr_highEta_blue_down_forward_phi_pT = nullptr;
0268     BinnedHistSet* bhs_pi0bkgr_highEta_yellow_up_forward_phi_pT = nullptr;
0269     BinnedHistSet* bhs_pi0bkgr_highEta_yellow_down_forward_phi_pT = nullptr;
0270     BinnedHistSet* bhs_pi0bkgr_blue_up_phi_xF = nullptr;
0271     BinnedHistSet* bhs_pi0bkgr_blue_down_phi_xF = nullptr;
0272     BinnedHistSet* bhs_pi0bkgr_yellow_up_phi_xF = nullptr;
0273     BinnedHistSet* bhs_pi0bkgr_yellow_down_phi_xF = nullptr;
0274     BinnedHistSet* bhs_pi0bkgr_blue_up_phi_eta = nullptr;
0275     BinnedHistSet* bhs_pi0bkgr_blue_down_phi_eta = nullptr;
0276     BinnedHistSet* bhs_pi0bkgr_yellow_up_phi_eta = nullptr;
0277     BinnedHistSet* bhs_pi0bkgr_yellow_down_phi_eta = nullptr;
0278     BinnedHistSet* bhs_pi0bkgr_blue_up_phi_vtxz = nullptr;
0279     BinnedHistSet* bhs_pi0bkgr_blue_down_phi_vtxz = nullptr;
0280     BinnedHistSet* bhs_pi0bkgr_yellow_up_phi_vtxz = nullptr;
0281     BinnedHistSet* bhs_pi0bkgr_yellow_down_phi_vtxz = nullptr;
0282
0283     BinnedHistSet* bhs_etabkgr_blue_up_phi_pT = nullptr;
0284     BinnedHistSet* bhs_etabkgr_blue_down_phi_pT = nullptr;
0285     BinnedHistSet* bhs_etabkgr_yellow_up_phi_pT = nullptr;
0286     BinnedHistSet* bhs_etabkgr_yellow_down_phi_pT = nullptr;
0287     BinnedHistSet* bhs_etabkgr_blue_up_forward_phi_pT = nullptr;
0288     BinnedHistSet* bhs_etabkgr_blue_down_forward_phi_pT = nullptr;
0289     BinnedHistSet* bhs_etabkgr_yellow_up_forward_phi_pT = nullptr;
0290     BinnedHistSet* bhs_etabkgr_yellow_down_forward_phi_pT = nullptr;
0291     BinnedHistSet* bhs_etabkgr_blue_up_backward_phi_pT = nullptr;
0292     BinnedHistSet* bhs_etabkgr_blue_down_backward_phi_pT = nullptr;
0293     BinnedHistSet* bhs_etabkgr_yellow_up_backward_phi_pT = nullptr;
0294     BinnedHistSet* bhs_etabkgr_yellow_down_backward_phi_pT = nullptr;
0295     BinnedHistSet* bhs_etabkgr_blue_up_phi_xF = nullptr;
0296     BinnedHistSet* bhs_etabkgr_blue_down_phi_xF = nullptr;
0297     BinnedHistSet* bhs_etabkgr_yellow_up_phi_xF = nullptr;
0298     BinnedHistSet* bhs_etabkgr_yellow_down_phi_xF = nullptr;
0299     BinnedHistSet* bhs_etabkgr_blue_up_phi_eta = nullptr;
0300     BinnedHistSet* bhs_etabkgr_blue_down_phi_eta = nullptr;
0301     BinnedHistSet* bhs_etabkgr_yellow_up_phi_eta = nullptr;
0302     BinnedHistSet* bhs_etabkgr_yellow_down_phi_eta = nullptr;
0303     BinnedHistSet* bhs_etabkgr_blue_up_phi_vtxz = nullptr;
0304     BinnedHistSet* bhs_etabkgr_blue_down_phi_vtxz = nullptr;
0305     BinnedHistSet* bhs_etabkgr_yellow_up_phi_vtxz = nullptr;
0306     BinnedHistSet* bhs_etabkgr_yellow_down_phi_vtxz = nullptr;
0307 };