18 #include "GaudiKernel/ITHistSvc.h"
19 #include "GaudiKernel/ServiceHandle.h"
37 m_myTree =
new TTree(
"data",
"data");
39 ATH_CHECK( m_jFEXOutputCollectionSGKey.initialize() );
41 m_myTree->Branch (
"smallRJet_eta", &m_smallRJet_eta);
42 m_myTree->Branch (
"smallRJet_phi", &m_smallRJet_phi);
43 m_myTree->Branch (
"smallRJet_ET", &m_smallRJet_ET);
44 m_myTree->Branch (
"smallRJet_clusterET", &m_smallRJet_clusterET);
45 m_myTree->Branch (
"smallRJet_nTOBs", &m_smallRJet_nTOBs,
"nTOBs");
46 m_myTree->Branch (
"smallRJet_isCentralTowerSeed", &m_smallRJet_isCentralTowerSeed);
48 m_myTree->Branch (
"smallRJetTOB_eta", &m_smallRJetTOB_eta);
49 m_myTree->Branch (
"smallRJetTOB_phi", &m_smallRJetTOB_phi);
50 m_myTree->Branch (
"smallRJetTOB_ET", &m_smallRJetTOB_ET);
51 m_myTree->Branch (
"smallRJetTOB_sat", &m_smallRJetTOB_sat);
52 m_myTree->Branch (
"smallRJetTOB_word", &m_smallRJetTOB_word);
53 m_myTree->Branch (
"smallRJetTOB_jfexID", &m_smallRJetTOB_jfexID);
54 m_myTree->Branch (
"smallRJetTOB_fpgaID", &m_smallRJetTOB_fpgaID);
57 m_myTree->Branch (
"largeRJet_eta", &m_largeRJet_eta);
58 m_myTree->Branch (
"largeRJet_phi", &m_largeRJet_phi);
59 m_myTree->Branch (
"largeRJet_ET", &m_largeRJet_ET);
60 m_myTree->Branch (
"largeRJet_nTOBs", &m_largeRJet_nTOBs);
62 m_myTree->Branch (
"largeRJetTOB_ET", &m_largeRJetTOB_ET);
63 m_myTree->Branch (
"largeRJetTOB_eta", &m_largeRJetTOB_eta);
64 m_myTree->Branch (
"largeRJetTOB_phi", &m_largeRJetTOB_phi);
65 m_myTree->Branch (
"largeRJetTOB_sat", &m_largeRJetTOB_sat);
67 m_myTree->Branch (
"largeRJetTOB_word", &m_largeRJetTOB_word);
68 m_myTree->Branch (
"largeRJetTOB_jfexID", &m_largeRJetTOB_jfexID);
69 m_myTree->Branch (
"largeRJetTOB_fpgaID", &m_largeRJetTOB_fpgaID);
71 m_myTree->Branch (
"tau_TT_ID", &m_tau_TT_ID);
72 m_myTree->Branch (
"tau_isLocalMax", &m_tau_isLocalMax);
73 m_myTree->Branch (
"tau_ET", &m_tau_ET);
74 m_myTree->Branch (
"tau_clusterET", &m_tau_clusterET);
75 m_myTree->Branch (
"tau_eta", &m_tau_eta);
76 m_myTree->Branch (
"tau_phi", &m_tau_phi);
77 m_myTree->Branch (
"tau_realeta", &m_tau_realeta);
78 m_myTree->Branch (
"tau_ISO", &m_tau_ISO);
79 m_myTree->Branch (
"tau_jFEXid", &m_tau_jFEXid);
80 m_myTree->Branch (
"tau_FPGAid", &m_tau_FPGAid);
82 m_myTree->Branch (
"tau_TOB_word", &m_tau_TOB_word);
83 m_myTree->Branch (
"tau_TOB_ET", &m_tau_TOB_ET);
84 m_myTree->Branch (
"tau_TOB_eta", &m_tau_TOB_eta);
85 m_myTree->Branch (
"tau_TOB_phi", &m_tau_TOB_phi);
86 m_myTree->Branch (
"tau_TOB_ISO", &m_tau_TOB_ISO);
87 m_myTree->Branch (
"tau_TOB_Sat", &m_tau_TOB_Sat);
90 m_myTree->Branch (
"pileup_FPGAid", &m_pileup_FPGAid);
91 m_myTree->Branch (
"pileup_jFEXid", &m_pileup_jFEXid);
92 m_myTree->Branch (
"pileup_rho_EM", &m_pileup_rho_EM);
93 m_myTree->Branch (
"pileup_rho_HAD1", &m_pileup_rho_HAD1);
94 m_myTree->Branch (
"pileup_rho_HAD2", &m_pileup_rho_HAD2);
95 m_myTree->Branch (
"pileup_rho_HAD3", &m_pileup_rho_HAD3);
96 m_myTree->Branch (
"pileup_rho_FCAL", &m_pileup_rho_FCAL);
98 m_myTree->Branch (
"pileup_map_ID", &m_pileup_map_ID);
99 m_myTree->Branch (
"pileup_map_Et_values_HAD_jet", &m_pileup_map_Et_values_HAD_jet);
100 m_myTree->Branch (
"pileup_map_Et_values_EM_jet", &m_pileup_map_Et_values_EM_jet);
101 m_myTree->Branch (
"pileup_map_Et_values_Total_jet", &m_pileup_map_Et_values_Total_jet);
102 m_myTree->Branch (
"pileup_map_Et_values_HAD_met", &m_pileup_map_Et_values_HAD_met);
103 m_myTree->Branch (
"pileup_map_Et_values_EM_met", &m_pileup_map_Et_values_EM_met);
104 m_myTree->Branch (
"pileup_map_Et_values_Total_met", &m_pileup_map_Et_values_Total_met);
106 return StatusCode::SUCCESS;
111 if(!jFEXOutputCollectionobj.
isValid()) {
112 ATH_MSG_ERROR(
"Could not retrieve jFEXOutputCollection " << m_jFEXOutputCollectionSGKey.key());
113 return StatusCode::FAILURE;
116 return StatusCode::SUCCESS;
118 CHECK(loadsmallRJetAlgoVariables(jFEXOutputCollectionobj));
119 CHECK(loadlargeRJetAlgoVariables(jFEXOutputCollectionobj));
120 CHECK(loadtauAlgoVariables(jFEXOutputCollectionobj));
121 CHECK(loadPileupVariables(jFEXOutputCollectionobj));
122 CHECK(loadPileupEt(jFEXOutputCollectionobj));
125 return StatusCode::SUCCESS;
130 return StatusCode::SUCCESS;
135 m_smallRJet_eta.clear();
136 m_smallRJet_phi.clear();
137 m_smallRJet_isCentralTowerSeed.clear();
138 m_smallRJet_ET.clear();
139 m_smallRJet_clusterET.clear();
140 m_smallRJetTOB_eta.clear();
141 m_smallRJetTOB_phi.clear();
142 m_smallRJetTOB_ET.clear();
143 m_smallRJetTOB_sat.clear();
144 m_smallRJetTOB_word.clear();
145 m_smallRJetTOB_jfexID.clear();
146 m_smallRJetTOB_fpgaID.clear();
148 for (
int i = 0;
i < jFEXOutputCollectionobj->
SRJetSize();
i++)
150 m_smallRJet_isCentralTowerSeed.push_back( jFEXOutputCollectionobj->
get_smallRJet(
i,
"smallRJet_isCentralTowerSeed"));
151 m_smallRJet_phi.push_back( jFEXOutputCollectionobj->
get_smallRJet(
i,
"smallRJet_phi"));
152 m_smallRJet_eta.push_back( jFEXOutputCollectionobj->
get_smallRJet(
i,
"smallRJet_eta"));
153 m_smallRJet_ET.push_back( jFEXOutputCollectionobj->
get_smallRJet(
i,
"smallRJet_ET"));
154 m_smallRJet_clusterET.push_back( jFEXOutputCollectionobj->
get_smallRJet(
i,
"smallRJet_clusterET"));
155 m_smallRJetTOB_eta.push_back( jFEXOutputCollectionobj->
get_smallRJet(
i,
"smallRJetTOB_eta"));
156 m_smallRJetTOB_phi.push_back( jFEXOutputCollectionobj->
get_smallRJet(
i,
"smallRJetTOB_phi"));
157 m_smallRJetTOB_ET.push_back( jFEXOutputCollectionobj->
get_smallRJet(
i,
"smallRJetTOB_ET"));
158 m_smallRJetTOB_sat.push_back( jFEXOutputCollectionobj->
get_smallRJet(
i,
"smallRJetTOB_sat"));
159 m_smallRJetTOB_word.push_back( jFEXOutputCollectionobj->
get_smallRJet(
i,
"smallRJetTOB_word"));
160 m_smallRJetTOB_jfexID.push_back( jFEXOutputCollectionobj->
get_smallRJet(
i,
"smallRJetTOB_jfexID"));
161 m_smallRJetTOB_fpgaID.push_back( jFEXOutputCollectionobj->
get_smallRJet(
i,
"smallRJetTOB_fpgaID"));
163 return StatusCode::SUCCESS;
168 m_largeRJet_ET.clear();
169 m_largeRJet_nTOBs.clear();
170 m_largeRJet_eta.clear();
171 m_largeRJet_phi.clear();
172 m_largeRJetTOB_ET.clear();
173 m_largeRJetTOB_eta.clear();
174 m_largeRJetTOB_phi.clear();
175 m_largeRJetTOB_sat.clear();
176 m_largeRJetTOB_word.clear();
177 m_largeRJetTOB_fpgaID.clear();
178 m_largeRJetTOB_jfexID.clear();
180 for (
int i = 0;
i < jFEXOutputCollectionobj->
LRJetSize();
i++)
182 m_largeRJet_ET.push_back( jFEXOutputCollectionobj->
get_largeRJet(
i,
"largeRJet_ET"));
183 m_largeRJet_nTOBs.push_back( jFEXOutputCollectionobj->
get_largeRJet(
i,
"largeRJet_nTOBs"));
184 m_largeRJet_eta.push_back( jFEXOutputCollectionobj->
get_largeRJet(
i,
"largeRJet_eta"));
185 m_largeRJet_phi.push_back( jFEXOutputCollectionobj->
get_largeRJet(
i,
"largeRJet_phi"));
186 m_largeRJetTOB_ET.push_back( jFEXOutputCollectionobj->
get_largeRJet(
i,
"largeRJetTOB_ET"));
187 m_largeRJetTOB_eta.push_back( jFEXOutputCollectionobj->
get_largeRJet(
i,
"largeRJetTOB_eta"));
188 m_largeRJetTOB_phi.push_back( jFEXOutputCollectionobj->
get_largeRJet(
i,
"largeRJetTOB_phi"));
189 m_largeRJetTOB_sat.push_back( jFEXOutputCollectionobj->
get_largeRJet(
i,
"largeRJetTOB_sat"));
190 m_largeRJetTOB_word.push_back( jFEXOutputCollectionobj->
get_largeRJet(
i,
"largeRJetTOB_word"));
191 m_largeRJetTOB_fpgaID.push_back(jFEXOutputCollectionobj->
get_largeRJet(
i,
"largeRJetTOB_fpgaID"));
192 m_largeRJetTOB_jfexID.push_back(jFEXOutputCollectionobj->
get_largeRJet(
i,
"largeRJetTOB_jfexID"));
195 return StatusCode::SUCCESS;
201 m_tau_jFEXid.clear();
202 m_tau_FPGAid.clear();
203 m_tau_isLocalMax.clear();
205 m_tau_clusterET.clear();
208 m_tau_realeta.clear();
211 m_tau_TOB_word.clear();
212 m_tau_TOB_ET.clear();
213 m_tau_TOB_eta.clear();
214 m_tau_TOB_phi.clear();
215 m_tau_TOB_ISO.clear();
216 m_tau_TOB_Sat.clear();
218 for (
int i = 0;
i < jFEXOutputCollectionobj->
TauSize();
i++)
220 m_tau_isLocalMax.push_back( jFEXOutputCollectionobj->
get_tau(
i,
"tau_isLocalMax"));
221 m_tau_TT_ID.push_back( jFEXOutputCollectionobj->
get_tau(
i,
"tau_TT_ID"));
222 m_tau_jFEXid.push_back( jFEXOutputCollectionobj->
get_tau(
i,
"tau_jFEXid"));
223 m_tau_FPGAid.push_back( jFEXOutputCollectionobj->
get_tau(
i,
"tau_FPGAid"));
224 m_tau_ET.push_back( jFEXOutputCollectionobj->
get_tau(
i,
"tau_ET"));
225 m_tau_clusterET.push_back( jFEXOutputCollectionobj->
get_tau(
i,
"tau_clusterET"));
226 m_tau_eta.push_back( jFEXOutputCollectionobj->
get_tau(
i,
"tau_eta"));
227 m_tau_phi.push_back( jFEXOutputCollectionobj->
get_tau(
i,
"tau_phi"));
228 m_tau_realeta.push_back( jFEXOutputCollectionobj->
get_tau(
i,
"tau_realeta"));
229 m_tau_ISO.push_back( jFEXOutputCollectionobj->
get_tau(
i,
"tau_ISO"));
231 m_tau_TOB_word.push_back( jFEXOutputCollectionobj->
get_tau(
i,
"tau_TOB_word"));
232 m_tau_TOB_ET.push_back( jFEXOutputCollectionobj->
get_tau(
i,
"tau_TOB_ET"));
233 m_tau_TOB_eta.push_back( jFEXOutputCollectionobj->
get_tau(
i,
"tau_TOB_eta"));
234 m_tau_TOB_phi.push_back( jFEXOutputCollectionobj->
get_tau(
i,
"tau_TOB_phi"));
235 m_tau_TOB_ISO.push_back( jFEXOutputCollectionobj->
get_tau(
i,
"tau_TOB_ISO"));
236 m_tau_TOB_Sat.push_back( jFEXOutputCollectionobj->
get_tau(
i,
"tau_TOB_Sat"));
239 return StatusCode::SUCCESS;
244 m_pileup_FPGAid.clear();
245 m_pileup_jFEXid.clear();
246 m_pileup_rho_EM.clear();
247 m_pileup_rho_HAD1.clear();
248 m_pileup_rho_HAD2.clear();
249 m_pileup_rho_HAD3.clear();
250 m_pileup_rho_FCAL.clear();
252 for (
int i = 0;
i < jFEXOutputCollectionobj->
PileupSize();
i++)
255 m_pileup_FPGAid.push_back( jFEXOutputCollectionobj->
get_pileup(
i,
"pileup_FPGAid"));
256 m_pileup_jFEXid.push_back( jFEXOutputCollectionobj->
get_pileup(
i,
"pileup_jFEXid"));
257 m_pileup_rho_EM.push_back( jFEXOutputCollectionobj->
get_pileup(
i,
"pileup_rho_EM"));
258 m_pileup_rho_HAD1.push_back(jFEXOutputCollectionobj->
get_pileup(
i,
"pileup_rho_HAD1"));
259 m_pileup_rho_HAD2.push_back(jFEXOutputCollectionobj->
get_pileup(
i,
"pileup_rho_HAD2"));
260 m_pileup_rho_HAD3.push_back(jFEXOutputCollectionobj->
get_pileup(
i,
"pileup_rho_HAD3"));
261 m_pileup_rho_FCAL.push_back(jFEXOutputCollectionobj->
get_pileup(
i,
"pileup_rho_FCAL"));
263 return StatusCode::SUCCESS;
268 m_pileup_map_ID.clear();
269 m_pileup_map_Et_values_HAD_jet.clear();
270 m_pileup_map_Et_values_EM_jet.clear();
271 m_pileup_map_Et_values_Total_jet.clear();
272 m_pileup_map_Et_values_HAD_met.clear();
273 m_pileup_map_Et_values_EM_met.clear();
274 m_pileup_map_Et_values_Total_met.clear();
278 m_pileup_map_ID.push_back( jFEXOutputCollectionobj->
get_pileup_map(
i,
"pileup_map_ID"));
279 m_pileup_map_Et_values_HAD_jet.push_back( jFEXOutputCollectionobj->
get_pileup_map(
i,
"pileup_map_Et_values_HAD_jet"));
280 m_pileup_map_Et_values_EM_jet.push_back( jFEXOutputCollectionobj->
get_pileup_map(
i,
"pileup_map_Et_values_EM_jet"));
281 m_pileup_map_Et_values_Total_jet.push_back(jFEXOutputCollectionobj->
get_pileup_map(
i,
"pileup_map_Et_values_Total_jet"));
282 m_pileup_map_Et_values_HAD_met.push_back( jFEXOutputCollectionobj->
get_pileup_map(
i,
"pileup_map_Et_values_HAD_met"));
283 m_pileup_map_Et_values_EM_met.push_back( jFEXOutputCollectionobj->
get_pileup_map(
i,
"pileup_map_Et_values_EM_met"));
284 m_pileup_map_Et_values_Total_met.push_back(jFEXOutputCollectionobj->
get_pileup_map(
i,
"pileup_map_Et_values_Total_met"));
289 return StatusCode::SUCCESS;