ATLAS Offline Software
PFOClusterMomentPlots.cxx
Go to the documentation of this file.
1 /*
2  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
3 */
4 
5 #include "PFOHistUtils/PFOClusterMomentPlots.h"
6 
7 namespace PFO {
8 
9  PFOClusterMomentPlots::PFOClusterMomentPlots(PlotBase* pParent, std::string sDir, std::string sFEContainerName) : PlotBase(pParent, sDir), m_sFEContainerName(sFEContainerName){
10 
11  m_FE_SECOND_R = nullptr;
12  m_FE_CENTER_LAMBDA = nullptr;
13  m_FE_ISOLATION = nullptr;
14  m_FE_ENG_BAD_CELLS = nullptr;
15  m_FE_N_BAD_CELLS = nullptr;
16  m_FE_BADLARQ_FRAC = nullptr;
17  m_FE_ENG_POS = nullptr;
18  m_FE_AVG_LAR_Q = nullptr;
19  m_FE_AVG_TILE_Q = nullptr;
20  m_FE_EM_PROBABILITY = nullptr;
21  m_FE_SECOND_LAMBDA = nullptr;
22 
23  m_FE_SECOND_R_etaBinA = nullptr;
24  m_FE_CENTER_LAMBDA_etaBinA = nullptr;
25  m_FE_ISOLATION_etaBinA = nullptr;
26  m_FE_ENG_BAD_CELLS_etaBinA = nullptr;
27  m_FE_N_BAD_CELLS_etaBinA = nullptr;
28  m_FE_BADLARQ_FRAC_etaBinA = nullptr;
29  m_FE_ENG_POS_etaBinA = nullptr;
30  m_FE_AVG_LAR_Q_etaBinA = nullptr;
31  m_FE_AVG_TILE_Q_etaBinA = nullptr;
32  m_FE_EM_PROBABILITY_etaBinA = nullptr;
33  m_FE_SECOND_LAMBDA_etaBinA = nullptr;
34 
35  m_FE_SECOND_R_etaBinB = nullptr;
36  m_FE_CENTER_LAMBDA_etaBinB = nullptr;
37  m_FE_ISOLATION_etaBinB = nullptr;
38  m_FE_ENG_BAD_CELLS_etaBinB = nullptr;
39  m_FE_N_BAD_CELLS_etaBinB = nullptr;
40  m_FE_BADLARQ_FRAC_etaBinB = nullptr;
41  m_FE_ENG_POS_etaBinB = nullptr;
42  m_FE_AVG_LAR_Q_etaBinB = nullptr;
43  m_FE_AVG_TILE_Q_etaBinB = nullptr;
44  m_FE_EM_PROBABILITY_etaBinB = nullptr;
45  m_FE_SECOND_LAMBDA_etaBinB = nullptr;
46 
47  m_FE_SECOND_R_etaBinC = nullptr;
48  m_FE_CENTER_LAMBDA_etaBinC = nullptr;
49  m_FE_ISOLATION_etaBinC = nullptr;
50  m_FE_ENG_BAD_CELLS_etaBinC = nullptr;
51  m_FE_N_BAD_CELLS_etaBinC = nullptr;
52  m_FE_BADLARQ_FRAC_etaBinC = nullptr;
53  m_FE_ENG_POS_etaBinC = nullptr;
54  m_FE_AVG_LAR_Q_etaBinC = nullptr;
55  m_FE_AVG_TILE_Q_etaBinC = nullptr;
56  m_FE_EM_PROBABILITY_etaBinC = nullptr;
57  m_FE_SECOND_LAMBDA_etaBinC = nullptr;
58 
59  m_FE_SECOND_R_etaBinD = nullptr;
60  m_FE_CENTER_LAMBDA_etaBinD = nullptr;
61  m_FE_ISOLATION_etaBinD = nullptr;
62  m_FE_ENG_BAD_CELLS_etaBinD = nullptr;
63  m_FE_N_BAD_CELLS_etaBinD = nullptr;
64  m_FE_BADLARQ_FRAC_etaBinD = nullptr;
65  m_FE_ENG_POS_etaBinD = nullptr;
66  m_FE_AVG_LAR_Q_etaBinD = nullptr;
67  m_FE_AVG_TILE_Q_etaBinD = nullptr;
68  m_FE_EM_PROBABILITY_etaBinD = nullptr;
69  m_FE_SECOND_LAMBDA_etaBinD = nullptr;
70  }
71 
72  void PFOClusterMomentPlots::initializePlots(){
73 
74  // FlowElement
75  if(!m_sFEContainerName.empty()){
76  m_FE_SECOND_R = Book1D("_SECOND_R",m_sFEContainerName + "_SECOND_R",60,-1.0,50.0);
77  m_FE_CENTER_LAMBDA = Book1D("_CENTER_LAMBDA",m_sFEContainerName + "_CENTER_LAMBDA",60,-50.0,3000.0);
78  m_FE_ISOLATION = Book1D("_ISOLATION",m_sFEContainerName + "_ISOLATION",60,-1.0,2.0);
79  m_FE_ENG_BAD_CELLS = Book1D("_ENG_BAD_CELLS",m_sFEContainerName + "_ENG_BAD_CELLS",60,-1.0,5);
80  m_FE_N_BAD_CELLS = Book1D("_N_BAD_CELLS",m_sFEContainerName + "_N_BAD_CELLS",30,-1.0,2.0);
81  m_FE_BADLARQ_FRAC = Book1D("_BADLARQ_FRAC",m_sFEContainerName + "_BADLARQ_FRAC",25,-1.0,1.5);
82  m_FE_ENG_POS = Book1D("_ENG_POS",m_sFEContainerName + "_ENG_POS",60,-100.0,10000.0);
83  m_FE_AVG_LAR_Q = Book1D("_AVG_LAR_Q",m_sFEContainerName + "_AVG_LAR_Q",31,-1000.0,30000.0);
84  m_FE_AVG_TILE_Q = Book1D("_AVG_TILE_Q",m_sFEContainerName + "_AVG_TILE_Q",21,-10.0,200.0);
85  m_FE_EM_PROBABILITY = Book1D("_EM_PROBABILITY",m_sFEContainerName + "_EM_PROBABILITY",21,-1.0,1.0);
86  m_FE_SECOND_LAMBDA = Book1D("_SECOND_LAMBDA",m_sFEContainerName + "_SECOND_LAMBDA",60,-1.0,3000.0);
87 
88  m_FE_SECOND_R_etaBinA = Book1D("_SECOND_R_A",m_sFEContainerName + "_SECOND_R (|eta| < 1.5)",60,-1.0,50.0);
89  m_FE_CENTER_LAMBDA_etaBinA = Book1D("_CENTER_LAMBDA_A",m_sFEContainerName + "_CENTER_LAMBDA (|eta| < 1.5)",60,-50.0,3000.0);
90  m_FE_ISOLATION_etaBinA = Book1D("_ISOLATION_A",m_sFEContainerName + "_ISOLATION (|eta| < 1.5)",60,-1.0,2.0);
91  m_FE_ENG_BAD_CELLS_etaBinA = Book1D("_ENG_BAD_CELLS_A",m_sFEContainerName + "_ENG_BAD_CELLS (|eta| < 1.5)",60,-1.0,5);
92  m_FE_N_BAD_CELLS_etaBinA = Book1D("_N_BAD_CELLS_A",m_sFEContainerName + "_N_BAD_CELLS (|eta| < 1.5)",30,-1.0,2.0);
93  m_FE_BADLARQ_FRAC_etaBinA = Book1D("_BADLARQ_FRAC_A",m_sFEContainerName + "_BADLARQ_FRAC (|eta| < 1.5)",25,-1.0,1.5);
94  m_FE_ENG_POS_etaBinA = Book1D("_ENG_POS_A",m_sFEContainerName + "_ENG_POS (|eta| < 1.5)",60,-100.0,10000.0);
95  m_FE_AVG_LAR_Q_etaBinA = Book1D("_AVG_LAR_Q_A",m_sFEContainerName + "_AVG_LAR_Q (|eta| < 1.5)",31,-1000.0,30000.0);
96  m_FE_AVG_TILE_Q_etaBinA = Book1D("_AVG_TILE_Q_A",m_sFEContainerName + "_AVG_TILE_Q (|eta| < 1.5)",21,-10.0,200.0);
97  m_FE_EM_PROBABILITY_etaBinA = Book1D("_EM_PROBABILITY_A",m_sFEContainerName + "_EM_PROBABILITY (|eta| < 1.5)",21,-1.0,1.0);
98  m_FE_SECOND_LAMBDA_etaBinA = Book1D("_SECOND_LAMBDA_A",m_sFEContainerName + "_SECOND_LAMBDA (|eta| < 1.5)",60,-1.0,3000.0);
99 
100  m_FE_SECOND_R_etaBinB = Book1D("_SECOND_R_B",m_sFEContainerName + "_SECOND_R (1.5 <= |eta| < 2.5)",60,-1.0,50.0);
101  m_FE_CENTER_LAMBDA_etaBinB = Book1D("_CENTER_LAMBDA_B",m_sFEContainerName + "_CENTER_LAMBDA (1.5 <= |eta| < 2.5)",60,-50.0,3000.0);
102  m_FE_ISOLATION_etaBinB = Book1D("_ISOLATION_B",m_sFEContainerName + "_ISOLATION (1.5 <= |eta| < 2.5)",60,-1.0,2.0);
103  m_FE_ENG_BAD_CELLS_etaBinB = Book1D("_ENG_BAD_CELLS_B",m_sFEContainerName + "_ENG_BAD_CELLS (1.5 <= |eta| < 2.5)",60,-1.0,5);
104  m_FE_N_BAD_CELLS_etaBinB = Book1D("_N_BAD_CELLS_B",m_sFEContainerName + "_N_BAD_CELLS (1.5 <= |eta| < 2.5)",30,-1.0,2.0);
105  m_FE_BADLARQ_FRAC_etaBinB = Book1D("_BADLARQ_FRAC_B",m_sFEContainerName + "_BADLARQ_FRAC (1.5 <= |eta| < 2.5)",25,-1.0,1.5);
106  m_FE_ENG_POS_etaBinB = Book1D("_ENG_POS_B",m_sFEContainerName + "_ENG_POS (1.5 <= |eta| < 2.5)",60,-100.0,10000.0);
107  m_FE_AVG_LAR_Q_etaBinB = Book1D("_AVG_LAR_Q_B",m_sFEContainerName + "_AVG_LAR_Q (1.5 <= |eta| < 2.5)",31,-1000.0,30000.0);
108  m_FE_AVG_TILE_Q_etaBinB = Book1D("_AVG_TILE_Q_B",m_sFEContainerName + "_AVG_TILE_Q (1.5 <= |eta| < 2.5)",21,-10.0,200.0);
109  m_FE_EM_PROBABILITY_etaBinB = Book1D("_EM_PROBABILITY_B",m_sFEContainerName + "_EM_PROBABILITY (1.5 <= |eta| < 2.5)",21,-1.0,1.0);
110  m_FE_SECOND_LAMBDA_etaBinB = Book1D("_SECOND_LAMBDA_B",m_sFEContainerName + "_SECOND_LAMBDA (1.5 <= |eta| < 2.5)",60,-1.0,3000.0);
111 
112  m_FE_SECOND_R_etaBinC = Book1D("_SECOND_R_C",m_sFEContainerName + "_SECOND_R (2.5 <= |eta| < 3.2)",60,-1.0,50.0);
113  m_FE_CENTER_LAMBDA_etaBinC = Book1D("_CENTER_LAMBDA_C",m_sFEContainerName + "_CENTER_LAMBDA (2.5 <= |eta| < 3.2)",60,-50.0,3000.0);
114  m_FE_ISOLATION_etaBinC = Book1D("_ISOLATION_C",m_sFEContainerName + "_ISOLATION (2.5 <= |eta| < 3.2)",60,-1.0,2.0);
115  m_FE_ENG_BAD_CELLS_etaBinC = Book1D("_ENG_BAD_CELLS_C",m_sFEContainerName + "_ENG_BAD_CELLS (2.5 <= |eta| < 3.2)",60,-1.0,5);
116  m_FE_N_BAD_CELLS_etaBinC = Book1D("_N_BAD_CELLS_C",m_sFEContainerName + "_N_BAD_CELLS (2.5 <= |eta| < 3.2)",30,-1.0,2.0);
117  m_FE_BADLARQ_FRAC_etaBinC = Book1D("_BADLARQ_FRAC_C",m_sFEContainerName + "_BADLARQ_FRAC (2.5 <= |eta| < 3.2)",25,-1.0,1.5);
118  m_FE_ENG_POS_etaBinC = Book1D("_ENG_POS_C",m_sFEContainerName + "_ENG_POS (2.5 <= |eta| < 3.2)",60,-100.0,10000.0);
119  m_FE_AVG_LAR_Q_etaBinC = Book1D("_AVG_LAR_Q_C",m_sFEContainerName + "_AVG_LAR_Q (2.5 <= |eta| < 3.2)",31,-1000.0,30000.0);
120  m_FE_AVG_TILE_Q_etaBinC = Book1D("_AVG_TILE_Q_C",m_sFEContainerName + "_AVG_TILE_Q (2.5 <= |eta| < 3.2)",21,-10.0,200.0);
121  m_FE_EM_PROBABILITY_etaBinC = Book1D("_EM_PROBABILITY_C",m_sFEContainerName + "_EM_PROBABILITY (2.5 <= |eta| < 3.2)",21,-1.0,1.0);
122  m_FE_SECOND_LAMBDA_etaBinC = Book1D("_SECOND_LAMBDA_C",m_sFEContainerName + "_SECOND_LAMBDA (2.5 <= |eta| < 3.2)",60,-1.0,3000.0);
123 
124  m_FE_SECOND_R_etaBinD = Book1D("_SECOND_R_D",m_sFEContainerName + "_SECOND_R (|eta| >= 3.2)",60,-1.0,50.0);
125  m_FE_CENTER_LAMBDA_etaBinD = Book1D("_CENTER_LAMBDA_D",m_sFEContainerName + "_CENTER_LAMBDA (|eta| >= 3.2)",60,-50.0,3000.0);
126  m_FE_ISOLATION_etaBinD = Book1D("_ISOLATION_D",m_sFEContainerName + "_ISOLATION (|eta| >= 3.2)",60,-1.0,2.0);
127  m_FE_ENG_BAD_CELLS_etaBinD = Book1D("_ENG_BAD_CELLS_D",m_sFEContainerName + "_ENG_BAD_CELLS (|eta| >= 3.2)",60,-1.0,5);
128  m_FE_N_BAD_CELLS_etaBinD = Book1D("_N_BAD_CELLS_D",m_sFEContainerName + "_N_BAD_CELLS (|eta| >= 3.2)",30,-1.0,2.0);
129  m_FE_BADLARQ_FRAC_etaBinD = Book1D("_BADLARQ_FRAC_D",m_sFEContainerName + "_BADLARQ_FRAC (|eta| >= 3.2)",25,-1.0,1.5);
130  m_FE_ENG_POS_etaBinD = Book1D("_ENG_POS_D",m_sFEContainerName + "_ENG_POS (|eta| >= 3.2)",60,-100.0,10000.0);
131  m_FE_AVG_LAR_Q_etaBinD = Book1D("_AVG_LAR_Q_D",m_sFEContainerName + "_AVG_LAR_Q (|eta| >= 3.2)",31,-1000.0,30000.0);
132  m_FE_AVG_TILE_Q_etaBinD = Book1D("_AVG_TILE_Q_D",m_sFEContainerName + "_AVG_TILE_Q (|eta| >= 3.2)",21,-10.0,200.0);
133  m_FE_EM_PROBABILITY_etaBinD = Book1D("_EM_PROBABILITY_D",m_sFEContainerName + "_EM_PROBABILITY (|eta| >= 3.2)",21,-1.0,1.0);
134  m_FE_SECOND_LAMBDA_etaBinD = Book1D("_SECOND_LAMBDA_D",m_sFEContainerName + "_SECOND_LAMBDA (|eta| >= 3.2)",60,-1.0,3000.0);
135  }
136  }
137 
138 
139  void PFOClusterMomentPlots::fill(const xAOD::FlowElement& FE, const xAOD::EventInfo& eventInfo){
140  float moment_SECOND_R = -1.0;
141  float moment_CENTER_LAMBDA = -1.0;
142  float moment_ISOLATION = -1.0;
143  float moment_ENG_BAD_CELLS=-1.0;
144  float moment_N_BAD_CELLS = -1.0;
145  float moment_BADLARQ_FRAC = -1.0;
146  float moment_ENG_POS = -1.0;
147  float moment_AVG_LAR_Q = -1.0;
148  float moment_AVG_TILE_Q = -1.0;
149  float moment_EM_PROBABILITY = -1.0;
150  float moment_SECOND_LAMBDA = -1.0;
151 
152  //as opposed to PFO which uses specific functions to grab the cluster moments, the auxdata is used
153  static const SG::AuxElement::ConstAccessor<float>acc_FE_moment_SECOND_R("SECOND_R");
154  static const SG::AuxElement::ConstAccessor<float>acc_FE_moment_CENTER_LAMBDA("CENTER_LAMBDA");
155  static const SG::AuxElement::ConstAccessor<float>acc_FE_moment_ISOLATION("ISOLATION");
156  static const SG::AuxElement::ConstAccessor<float>acc_FE_moment_ENG_BAD_CELLS("ENG_BAD_CELLS");
157  static const SG::AuxElement::ConstAccessor<float>acc_FE_moment_N_BAD_CELLS("N_BAD_CELLS");
158  static const SG::AuxElement::ConstAccessor<float>acc_FE_moment_BADLARQ_FRAC("BADLARQ_FRAC");
159  static const SG::AuxElement::ConstAccessor<float>acc_FE_moment_ENG_POS("ENG_POS");
160  static const SG::AuxElement::ConstAccessor<float>acc_FE_moment_AVG_LAR_Q("AVG_LAR_Q");
161  static const SG::AuxElement::ConstAccessor<float>acc_FE_moment_AVG_TILE_Q("AVG_TILE_Q");
162  static const SG::AuxElement::ConstAccessor<float>acc_FE_moment_EM_PROBABILITY("EM_PROBABILITY");
163  static const SG::AuxElement::ConstAccessor<float>acc_FE_moment_SECOND_LAMBDA("SECOND_LAMBDA");
164 
165  //use accessors to retrieve the auxvars
166  if(acc_FE_moment_SECOND_R.isAvailable(FE))
167  moment_SECOND_R=acc_FE_moment_SECOND_R(FE);
168 
169  if(acc_FE_moment_CENTER_LAMBDA.isAvailable(FE))
170  moment_CENTER_LAMBDA=acc_FE_moment_CENTER_LAMBDA(FE);
171 
172  if(acc_FE_moment_ISOLATION.isAvailable(FE))
173  moment_ISOLATION=acc_FE_moment_ISOLATION(FE);
174 
175  if(acc_FE_moment_ENG_BAD_CELLS.isAvailable(FE))
176  moment_ENG_BAD_CELLS=acc_FE_moment_ENG_BAD_CELLS(FE);
177 
178  if(acc_FE_moment_N_BAD_CELLS.isAvailable(FE))
179  moment_N_BAD_CELLS=acc_FE_moment_N_BAD_CELLS(FE);
180 
181  if(acc_FE_moment_BADLARQ_FRAC.isAvailable(FE))
182  moment_BADLARQ_FRAC=acc_FE_moment_BADLARQ_FRAC(FE);
183 
184  if(acc_FE_moment_ENG_POS.isAvailable(FE))
185  moment_ENG_POS=acc_FE_moment_ENG_POS(FE);
186 
187  if(acc_FE_moment_AVG_LAR_Q.isAvailable(FE))
188  moment_AVG_LAR_Q=acc_FE_moment_AVG_LAR_Q(FE);
189 
190  if(acc_FE_moment_AVG_TILE_Q.isAvailable(FE))
191  moment_AVG_TILE_Q=acc_FE_moment_AVG_TILE_Q(FE);
192 
193  if(acc_FE_moment_EM_PROBABILITY.isAvailable(FE))
194  moment_EM_PROBABILITY=acc_FE_moment_EM_PROBABILITY(FE);
195 
196  if(acc_FE_moment_SECOND_LAMBDA.isAvailable(FE))
197  moment_SECOND_LAMBDA=acc_FE_moment_SECOND_LAMBDA(FE);
198 
199  float FE_eta=FE.eta();
200  m_FE_SECOND_R->Fill(moment_SECOND_R,eventInfo.beamSpotWeight());
201  m_FE_CENTER_LAMBDA->Fill(moment_CENTER_LAMBDA,eventInfo.beamSpotWeight());
202  m_FE_ISOLATION->Fill(moment_ISOLATION,eventInfo.beamSpotWeight());
203  m_FE_ENG_BAD_CELLS->Fill(moment_ENG_BAD_CELLS,eventInfo.beamSpotWeight());
204  m_FE_N_BAD_CELLS->Fill(moment_N_BAD_CELLS,eventInfo.beamSpotWeight());
205  m_FE_BADLARQ_FRAC->Fill(moment_BADLARQ_FRAC,eventInfo.beamSpotWeight());
206  m_FE_ENG_POS->Fill(moment_ENG_POS,eventInfo.beamSpotWeight());
207  m_FE_AVG_LAR_Q->Fill(moment_AVG_LAR_Q,eventInfo.beamSpotWeight());
208  m_FE_AVG_TILE_Q->Fill(moment_AVG_TILE_Q,eventInfo.beamSpotWeight());
209  m_FE_EM_PROBABILITY->Fill(moment_EM_PROBABILITY,eventInfo.beamSpotWeight());
210  m_FE_SECOND_LAMBDA->Fill(moment_SECOND_LAMBDA,eventInfo.beamSpotWeight());
211 
212  if (fabs(FE_eta) < 1.5){
213  m_FE_SECOND_R_etaBinA->Fill(moment_SECOND_R,eventInfo.beamSpotWeight());
214  m_FE_CENTER_LAMBDA_etaBinA->Fill(moment_CENTER_LAMBDA,eventInfo.beamSpotWeight());
215  m_FE_ISOLATION_etaBinA->Fill(moment_ISOLATION,eventInfo.beamSpotWeight());
216  m_FE_ENG_BAD_CELLS_etaBinA->Fill(moment_ENG_BAD_CELLS,eventInfo.beamSpotWeight());
217  m_FE_N_BAD_CELLS_etaBinA->Fill(moment_N_BAD_CELLS,eventInfo.beamSpotWeight());
218  m_FE_BADLARQ_FRAC_etaBinA->Fill(moment_BADLARQ_FRAC,eventInfo.beamSpotWeight());
219  m_FE_ENG_POS_etaBinA->Fill(moment_ENG_POS,eventInfo.beamSpotWeight());
220  m_FE_AVG_LAR_Q_etaBinA->Fill(moment_AVG_LAR_Q,eventInfo.beamSpotWeight());
221  m_FE_AVG_TILE_Q_etaBinA->Fill(moment_AVG_TILE_Q,eventInfo.beamSpotWeight());
222  m_FE_EM_PROBABILITY_etaBinA->Fill(moment_EM_PROBABILITY,eventInfo.beamSpotWeight());
223  m_FE_SECOND_LAMBDA_etaBinA->Fill(moment_SECOND_LAMBDA,eventInfo.beamSpotWeight());
224  }//|eta| < 1.5
225  else if (fabs(FE_eta) < 2.5){
226  m_FE_SECOND_R_etaBinB->Fill(moment_SECOND_R,eventInfo.beamSpotWeight());
227  m_FE_CENTER_LAMBDA_etaBinB->Fill(moment_CENTER_LAMBDA,eventInfo.beamSpotWeight());
228  m_FE_ISOLATION_etaBinB->Fill(moment_ISOLATION,eventInfo.beamSpotWeight());
229  m_FE_ENG_BAD_CELLS_etaBinB->Fill(moment_ENG_BAD_CELLS,eventInfo.beamSpotWeight());
230  m_FE_N_BAD_CELLS_etaBinB->Fill(moment_N_BAD_CELLS,eventInfo.beamSpotWeight());
231  m_FE_BADLARQ_FRAC_etaBinB->Fill(moment_BADLARQ_FRAC,eventInfo.beamSpotWeight());
232  m_FE_ENG_POS_etaBinB->Fill(moment_ENG_POS,eventInfo.beamSpotWeight());
233  m_FE_AVG_LAR_Q_etaBinB->Fill(moment_AVG_LAR_Q,eventInfo.beamSpotWeight());
234  m_FE_AVG_TILE_Q_etaBinB->Fill(moment_AVG_TILE_Q,eventInfo.beamSpotWeight());
235  m_FE_EM_PROBABILITY_etaBinB->Fill(moment_EM_PROBABILITY,eventInfo.beamSpotWeight());
236  m_FE_SECOND_LAMBDA_etaBinB->Fill(moment_SECOND_LAMBDA,eventInfo.beamSpotWeight());
237  }
238  else if (fabs(FE_eta) < 3.2){
239  m_FE_SECOND_R_etaBinC->Fill(moment_SECOND_R,eventInfo.beamSpotWeight());
240  m_FE_CENTER_LAMBDA_etaBinC->Fill(moment_CENTER_LAMBDA,eventInfo.beamSpotWeight());
241  m_FE_ISOLATION_etaBinC->Fill(moment_ISOLATION,eventInfo.beamSpotWeight());
242  m_FE_ENG_BAD_CELLS_etaBinC->Fill(moment_ENG_BAD_CELLS,eventInfo.beamSpotWeight());
243  m_FE_N_BAD_CELLS_etaBinC->Fill(moment_N_BAD_CELLS,eventInfo.beamSpotWeight());
244  m_FE_BADLARQ_FRAC_etaBinC->Fill(moment_BADLARQ_FRAC,eventInfo.beamSpotWeight());
245  m_FE_ENG_POS_etaBinC->Fill(moment_ENG_POS,eventInfo.beamSpotWeight());
246  m_FE_AVG_LAR_Q_etaBinC->Fill(moment_AVG_LAR_Q,eventInfo.beamSpotWeight());
247  m_FE_AVG_TILE_Q_etaBinC->Fill(moment_AVG_TILE_Q,eventInfo.beamSpotWeight());
248  m_FE_EM_PROBABILITY_etaBinC->Fill(moment_EM_PROBABILITY,eventInfo.beamSpotWeight());
249  m_FE_SECOND_LAMBDA_etaBinC->Fill(moment_SECOND_LAMBDA,eventInfo.beamSpotWeight());
250  }
251  else{
252  m_FE_SECOND_R_etaBinD->Fill(moment_SECOND_R,eventInfo.beamSpotWeight());
253  m_FE_CENTER_LAMBDA_etaBinD->Fill(moment_CENTER_LAMBDA,eventInfo.beamSpotWeight());
254  m_FE_ISOLATION_etaBinD->Fill(moment_ISOLATION,eventInfo.beamSpotWeight());
255  m_FE_ENG_BAD_CELLS_etaBinD->Fill(moment_ENG_BAD_CELLS,eventInfo.beamSpotWeight());
256  m_FE_N_BAD_CELLS_etaBinD->Fill(moment_N_BAD_CELLS,eventInfo.beamSpotWeight());
257  m_FE_BADLARQ_FRAC_etaBinD->Fill(moment_BADLARQ_FRAC,eventInfo.beamSpotWeight());
258  m_FE_ENG_POS_etaBinD->Fill(moment_ENG_POS,eventInfo.beamSpotWeight());
259  m_FE_AVG_LAR_Q_etaBinD->Fill(moment_AVG_LAR_Q,eventInfo.beamSpotWeight());
260  m_FE_AVG_TILE_Q_etaBinD->Fill(moment_AVG_TILE_Q,eventInfo.beamSpotWeight());
261  m_FE_EM_PROBABILITY_etaBinD->Fill(moment_EM_PROBABILITY,eventInfo.beamSpotWeight());
262  m_FE_SECOND_LAMBDA_etaBinD->Fill(moment_SECOND_LAMBDA,eventInfo.beamSpotWeight());
263  }
264 
265 
266  }
267 }
PFO::PFOClusterMomentPlots::m_FE_EM_PROBABILITY
TH1 * m_FE_EM_PROBABILITY
Definition: PFOClusterMomentPlots.h:33
PFO::PFOClusterMomentPlots::fill
void fill(const xAOD::FlowElement &FE, const xAOD::EventInfo &eventInfo)
Definition: PFOClusterMomentPlots.cxx:139
PFO::PFOClusterMomentPlots::m_FE_ENG_BAD_CELLS_etaBinD
TH1 * m_FE_ENG_BAD_CELLS_etaBinD
Definition: PFOClusterMomentPlots.h:76
PFO::PFOClusterMomentPlots::PFOClusterMomentPlots
PFOClusterMomentPlots(PlotBase *pParent, std::string sDir, std::string sFEContainerName)
Definition: PFOClusterMomentPlots.cxx:9
PFO::PFOClusterMomentPlots::m_FE_BADLARQ_FRAC_etaBinA
TH1 * m_FE_BADLARQ_FRAC_etaBinA
Definition: PFOClusterMomentPlots.h:42
PFO::PFOClusterMomentPlots::m_FE_EM_PROBABILITY_etaBinC
TH1 * m_FE_EM_PROBABILITY_etaBinC
Definition: PFOClusterMomentPlots.h:70
PFO::PFOClusterMomentPlots::m_FE_ENG_POS
TH1 * m_FE_ENG_POS
Definition: PFOClusterMomentPlots.h:30
PlotBase::Book1D
TH1F * Book1D(const std::string &name, const std::string &labels, int nBins, float start, float end, bool prependDir=true)
Book a TH1D histogram.
Definition: PlotBase.cxx:88
PFO::PFOClusterMomentPlots::m_FE_ENG_BAD_CELLS
TH1 * m_FE_ENG_BAD_CELLS
Definition: PFOClusterMomentPlots.h:27
PFO::PFOClusterMomentPlots::m_FE_AVG_TILE_Q
TH1 * m_FE_AVG_TILE_Q
Definition: PFOClusterMomentPlots.h:32
PFO::PFOClusterMomentPlots::m_FE_CENTER_LAMBDA_etaBinA
TH1 * m_FE_CENTER_LAMBDA_etaBinA
Definition: PFOClusterMomentPlots.h:38
PlotBase
Definition: PlotBase.h:33
PFO::PFOClusterMomentPlots::m_FE_N_BAD_CELLS_etaBinB
TH1 * m_FE_N_BAD_CELLS_etaBinB
Definition: PFOClusterMomentPlots.h:53
PFO::PFOClusterMomentPlots::m_FE_N_BAD_CELLS_etaBinC
TH1 * m_FE_N_BAD_CELLS_etaBinC
Definition: PFOClusterMomentPlots.h:65
PFO::PFOClusterMomentPlots::m_FE_BADLARQ_FRAC_etaBinB
TH1 * m_FE_BADLARQ_FRAC_etaBinB
Definition: PFOClusterMomentPlots.h:54
python.copyTCTOutput.sDir
sDir
Definition: copyTCTOutput.py:60
PFO::PFOClusterMomentPlots::m_FE_AVG_LAR_Q_etaBinC
TH1 * m_FE_AVG_LAR_Q_etaBinC
Definition: PFOClusterMomentPlots.h:68
PFO::PFOClusterMomentPlots::m_FE_N_BAD_CELLS
TH1 * m_FE_N_BAD_CELLS
Definition: PFOClusterMomentPlots.h:28
PFO::PFOClusterMomentPlots::m_FE_ENG_BAD_CELLS_etaBinC
TH1 * m_FE_ENG_BAD_CELLS_etaBinC
Definition: PFOClusterMomentPlots.h:64
PFO::PFOClusterMomentPlots::m_FE_AVG_TILE_Q_etaBinA
TH1 * m_FE_AVG_TILE_Q_etaBinA
Definition: PFOClusterMomentPlots.h:45
PFO::PFOClusterMomentPlots::m_FE_SECOND_LAMBDA_etaBinA
TH1 * m_FE_SECOND_LAMBDA_etaBinA
Definition: PFOClusterMomentPlots.h:47
PFO::PFOClusterMomentPlots::m_FE_SECOND_R_etaBinB
TH1 * m_FE_SECOND_R_etaBinB
Definition: PFOClusterMomentPlots.h:49
PFO::PFOClusterMomentPlots::m_FE_BADLARQ_FRAC_etaBinD
TH1 * m_FE_BADLARQ_FRAC_etaBinD
Definition: PFOClusterMomentPlots.h:78
SG::ConstAccessor
Helper class to provide constant type-safe access to aux data.
Definition: ConstAccessor.h:54
PFO::PFOClusterMomentPlots::m_FE_CENTER_LAMBDA_etaBinB
TH1 * m_FE_CENTER_LAMBDA_etaBinB
Definition: PFOClusterMomentPlots.h:50
PFO::PFOClusterMomentPlots::m_FE_ENG_POS_etaBinB
TH1 * m_FE_ENG_POS_etaBinB
Definition: PFOClusterMomentPlots.h:55
PFO::PFOClusterMomentPlots::m_FE_AVG_TILE_Q_etaBinD
TH1 * m_FE_AVG_TILE_Q_etaBinD
Definition: PFOClusterMomentPlots.h:81
PFO::PFOClusterMomentPlots::m_FE_BADLARQ_FRAC
TH1 * m_FE_BADLARQ_FRAC
Definition: PFOClusterMomentPlots.h:29
PFO::PFOClusterMomentPlots::m_FE_SECOND_R
TH1 * m_FE_SECOND_R
Histograms inclusive in eta.
Definition: PFOClusterMomentPlots.h:24
PFOClusterMomentPlots.h
PFO::PFOClusterMomentPlots::m_FE_ISOLATION_etaBinA
TH1 * m_FE_ISOLATION_etaBinA
Definition: PFOClusterMomentPlots.h:39
PFO::PFOClusterMomentPlots::m_FE_N_BAD_CELLS_etaBinA
TH1 * m_FE_N_BAD_CELLS_etaBinA
Definition: PFOClusterMomentPlots.h:41
PFO::PFOClusterMomentPlots::m_FE_CENTER_LAMBDA_etaBinD
TH1 * m_FE_CENTER_LAMBDA_etaBinD
Definition: PFOClusterMomentPlots.h:74
PFO::PFOClusterMomentPlots::m_FE_N_BAD_CELLS_etaBinD
TH1 * m_FE_N_BAD_CELLS_etaBinD
Definition: PFOClusterMomentPlots.h:77
PFO::PFOClusterMomentPlots::m_FE_CENTER_LAMBDA_etaBinC
TH1 * m_FE_CENTER_LAMBDA_etaBinC
Definition: PFOClusterMomentPlots.h:62
PFO::PFOClusterMomentPlots::m_FE_ISOLATION_etaBinC
TH1 * m_FE_ISOLATION_etaBinC
Definition: PFOClusterMomentPlots.h:63
PFO::PFOClusterMomentPlots::m_FE_ENG_POS_etaBinA
TH1 * m_FE_ENG_POS_etaBinA
Definition: PFOClusterMomentPlots.h:43
PFO::PFOClusterMomentPlots::m_sFEContainerName
std::string m_sFEContainerName
Definition: PFOClusterMomentPlots.h:87
PFO::PFOClusterMomentPlots::m_FE_ENG_POS_etaBinD
TH1 * m_FE_ENG_POS_etaBinD
Definition: PFOClusterMomentPlots.h:79
PFO::PFOClusterMomentPlots::initializePlots
void initializePlots()
Definition: PFOClusterMomentPlots.cxx:72
PFO::PFOClusterMomentPlots::m_FE_EM_PROBABILITY_etaBinB
TH1 * m_FE_EM_PROBABILITY_etaBinB
Definition: PFOClusterMomentPlots.h:58
PFO::PFOClusterMomentPlots::m_FE_CENTER_LAMBDA
TH1 * m_FE_CENTER_LAMBDA
Definition: PFOClusterMomentPlots.h:25
TH1::Fill
int Fill(double)
Definition: rootspy.cxx:285
PFO::PFOClusterMomentPlots::m_FE_SECOND_LAMBDA_etaBinD
TH1 * m_FE_SECOND_LAMBDA_etaBinD
Definition: PFOClusterMomentPlots.h:83
PFO::PFOClusterMomentPlots::m_FE_ISOLATION_etaBinB
TH1 * m_FE_ISOLATION_etaBinB
Definition: PFOClusterMomentPlots.h:51
PFO
Definition: ClusterMomentPlots.h:13
PFO::PFOClusterMomentPlots::m_FE_ENG_POS_etaBinC
TH1 * m_FE_ENG_POS_etaBinC
Definition: PFOClusterMomentPlots.h:67
PFO::PFOClusterMomentPlots::m_FE_SECOND_LAMBDA_etaBinB
TH1 * m_FE_SECOND_LAMBDA_etaBinB
Definition: PFOClusterMomentPlots.h:59
PFO::PFOClusterMomentPlots::m_FE_SECOND_R_etaBinD
TH1 * m_FE_SECOND_R_etaBinD
Definition: PFOClusterMomentPlots.h:73
PFO::PFOClusterMomentPlots::m_FE_AVG_LAR_Q
TH1 * m_FE_AVG_LAR_Q
Definition: PFOClusterMomentPlots.h:31
PFO::PFOClusterMomentPlots::m_FE_AVG_LAR_Q_etaBinA
TH1 * m_FE_AVG_LAR_Q_etaBinA
Definition: PFOClusterMomentPlots.h:44
xAOD::EventInfo_v1
Class describing the basic event information.
Definition: EventInfo_v1.h:43
PFO::PFOClusterMomentPlots::m_FE_EM_PROBABILITY_etaBinD
TH1 * m_FE_EM_PROBABILITY_etaBinD
Definition: PFOClusterMomentPlots.h:82
PFO::PFOClusterMomentPlots::m_FE_BADLARQ_FRAC_etaBinC
TH1 * m_FE_BADLARQ_FRAC_etaBinC
Definition: PFOClusterMomentPlots.h:66
PFO::PFOClusterMomentPlots::m_FE_AVG_LAR_Q_etaBinB
TH1 * m_FE_AVG_LAR_Q_etaBinB
Definition: PFOClusterMomentPlots.h:56
PFO::PFOClusterMomentPlots::m_FE_AVG_TILE_Q_etaBinB
TH1 * m_FE_AVG_TILE_Q_etaBinB
Definition: PFOClusterMomentPlots.h:57
xAOD::FlowElement_v1::eta
virtual double eta() const override
The pseudorapidity ( ) of the particle.
PFO::PFOClusterMomentPlots::m_FE_SECOND_R_etaBinA
TH1 * m_FE_SECOND_R_etaBinA
Same histograms binned in eta regions.
Definition: PFOClusterMomentPlots.h:37
PFO::PFOClusterMomentPlots::m_FE_ENG_BAD_CELLS_etaBinA
TH1 * m_FE_ENG_BAD_CELLS_etaBinA
Definition: PFOClusterMomentPlots.h:40
PFO::PFOClusterMomentPlots::m_FE_EM_PROBABILITY_etaBinA
TH1 * m_FE_EM_PROBABILITY_etaBinA
Definition: PFOClusterMomentPlots.h:46
PFO::PFOClusterMomentPlots::m_FE_SECOND_LAMBDA
TH1 * m_FE_SECOND_LAMBDA
Definition: PFOClusterMomentPlots.h:34
PFO::PFOClusterMomentPlots::m_FE_SECOND_R_etaBinC
TH1 * m_FE_SECOND_R_etaBinC
Definition: PFOClusterMomentPlots.h:61
SG::ConstAccessor::isAvailable
bool isAvailable(const ELT &e) const
Test to see if this variable exists in the store.
PFO::PFOClusterMomentPlots::m_FE_ISOLATION
TH1 * m_FE_ISOLATION
Definition: PFOClusterMomentPlots.h:26
PFO::PFOClusterMomentPlots::m_FE_ISOLATION_etaBinD
TH1 * m_FE_ISOLATION_etaBinD
Definition: PFOClusterMomentPlots.h:75
PFO::PFOClusterMomentPlots::m_FE_AVG_TILE_Q_etaBinC
TH1 * m_FE_AVG_TILE_Q_etaBinC
Definition: PFOClusterMomentPlots.h:69
PFO::PFOClusterMomentPlots::m_FE_ENG_BAD_CELLS_etaBinB
TH1 * m_FE_ENG_BAD_CELLS_etaBinB
Definition: PFOClusterMomentPlots.h:52
PFO::PFOClusterMomentPlots::m_FE_AVG_LAR_Q_etaBinD
TH1 * m_FE_AVG_LAR_Q_etaBinD
Definition: PFOClusterMomentPlots.h:80
xAOD::EventInfo_v1::beamSpotWeight
float beamSpotWeight() const
Weight for beam spot size reweighting.
Definition: EventInfo_v1.cxx:970
PFO::PFOClusterMomentPlots::m_FE_SECOND_LAMBDA_etaBinC
TH1 * m_FE_SECOND_LAMBDA_etaBinC
Definition: PFOClusterMomentPlots.h:71
xAOD::FlowElement_v1
A detector object made of other lower level object(s)
Definition: FlowElement_v1.h:25
Generated on Fri Jun 28 2024 21:24:30 for ATLAS Offline Software by   doxygen 1.8.18