35 m_h_S1_xy_A =
new TH2D(
"h_mm_S1_xy_A",
"S1_xy_A", 1200,-6000., 6000., 1200, -6000., 6000.);
39 m_h_S1_rz_A =
new TH2D(
"h_mm_S1_rz_A",
"S1_rz_A", 1000,7000., 8000., 5000, 0., 5000.);
43 m_h_S1_r_A =
new TH1D(
"h_mm_S1_r_A",
"S1_r_A", 1000,0., 14000.);
47 m_h_S1_xy_C =
new TH2D(
"h_mm_S1_xy_C",
"S1_xy_C", 1200,-6000., 6000., 1200, -6000., 6000.);
51 m_h_S1_rz_C =
new TH2D(
"h_mm_S1_rz_C",
"S1_rz_C", 1000,-8000., -7000., 5000, 0., 5000.);
55 m_h_S1_r_C =
new TH1D(
"h_mm_S1_r_C",
"S1_r_C", 1000,0., 14000.);
61 m_h_S2_xy_A =
new TH2D(
"h_mm_S2_xy_A",
"S2_xy_A", 1200,-6000., 6000., 1200, -6000., 6000.);
65 m_h_S2_rz_A =
new TH2D(
"h_mm_S2_rz_A",
"S2_rz_A", 1000,7000., 8000., 5000, 0., 5000.);
69 m_h_S2_r_A =
new TH1D(
"h_mm_S2_r_A",
"S2_r_A", 1000,0., 14000.);
73 m_h_S2_xy_C =
new TH2D(
"h_mm_S2_xy_C",
"S2_xy_C", 1200,-6000., 6000., 1200, -6000., 6000.);
77 m_h_S2_rz_C =
new TH2D(
"h_mm_S2_rz_C",
"S2_rz_C", 1000,-8000., -7000., 5000, 0., 5000.);
81 m_h_S2_r_C =
new TH1D(
"h_mm_S2_r_C",
"S2_r_C", 1000,0., 14000.);
87 m_h_S_xy_A =
new TH2D(
"h_mm_S_xy_A",
"S_xy_A", 1200,-6000., 6000., 1200, -6000., 6000.);
91 m_h_S_xy_C =
new TH2D(
"h_mm_S_xy_C",
"S_xy_C", 1200,-6000., 6000., 1200, -6000., 6000.);
95 m_h_S_rz_A =
new TH2D(
"h_mm_S_rz_A",
"S_rz_A", 1000,7000., 8000., 5000, 0., 5000.);
99 m_h_S_rz_C =
new TH2D(
"h_mm_S_rz_C",
"S_rz_C", 1000,-8000., -7000., 5000, 0., 5000.);
105 m_h_xy_A =
new TH2D(
"h_mm_xy_A",
"xy_A", 1200,-6000., 6000., 1200, -6000., 6000.);
109 m_h_xy_C =
new TH2D(
"h_mm_xy_C",
"xy_C", 1200,-6000., 6000., 1200, -6000., 6000.);
113 m_h_rz_A =
new TH2D(
"h_mm_rz_A",
"rz_A", 1000,7000., 8000., 5000, 0., 5000.);
117 m_h_rz_C =
new TH2D(
"h_mm_rz_C",
"rz_C", 1000,-8000., -7000., 5000, 0., 5000.);
123 m_h_L1_xy_A =
new TH2D(
"h_mm_L1_xy_A",
"L1_xy_A", 1200,-6000., 6000., 1200, -6000., 6000.);
127 m_h_L1_rz_A =
new TH2D(
"h_mm_L1_rz_A",
"L1_rz_A", 1000,7000., 8000., 5000, 0., 5000.);
131 m_h_L1_r_A =
new TH1D(
"h_mm_L1_r_A",
"L1_r_A", 1000,0., 14000.);
135 m_h_L1_xy_C =
new TH2D(
"h_mm_L1_xy_C",
"L1_xy_C", 1200,-6000., 6000., 1200, -6000., 6000.);
139 m_h_L1_rz_C =
new TH2D(
"h_mm_L1_rz_C",
"L1_rz_C", 1000,-8000., -7000., 5000, 0., 5000.);
143 m_h_L1_r_C =
new TH1D(
"h_mm_L1_r_C",
"L1_r_C", 1000,0., 14000.);
149 m_h_L2_xy_A =
new TH2D(
"h_mm_L2_xy_A",
"L2_xy_A", 1200,-6000., 6000., 1200, -6000., 6000.);
153 m_h_L2_rz_A =
new TH2D(
"h_mm_L2_rz_A",
"L2_rz_A", 1000,7000., 8000., 5000, 0., 5000.);
157 m_h_L2_r_A =
new TH1D(
"h_mm_L2_r_A",
"L2_r_A", 1000,0., 14000.);
161 m_h_L2_xy_C =
new TH2D(
"h_mm_L2_xy_C",
"L2_xy_C", 1200,-6000., 6000., 1200, -6000., 6000.);
165 m_h_L2_rz_C =
new TH2D(
"h_mm_L2_rz_C",
"L2_rz_C", 1000,-8000., -7000., 5000, 0., 5000.);
169 m_h_L2_r_C =
new TH1D(
"h_mm_L2_r_C",
"L2_r_C", 1000,0., 14000.);
175 m_h_L_xy_A =
new TH2D(
"h_mm_L_xy_A",
"L_xy_A", 1200,-6000., 6000., 1200, -6000., 6000.);
179 m_h_L_rz_A =
new TH2D(
"h_mm_L_rz_A",
"L_rz_A", 1000,7000., 8000., 5000, 0., 5000.);
183 m_h_L_xy_C =
new TH2D(
"h_mm_L_xy_C",
"L_xy_C", 1200,-6000., 6000., 1200, -6000., 6000.);
187 m_h_L_rz_C =
new TH2D(
"h_mm_L_rz_C",
"L_rz_C", 1000,-8000., -7000., 5000, 0., 5000.);
193 return StatusCode::SUCCESS;
200 const EventContext& ctx{Gaudi::Hive::currentContext()};
208 int simId = (*i_hit).MMId();
210 int sim_side = hitHelper->
GetSide(simId);
212 static const std::string s_m1s1(
"M1S1");
213 static const std::string s_m2s1(
"M2S1");
214 static const std::string s_m1s2(
"M1S2");
215 static const std::string s_m2s2(
"M2S2");
216 static const std::string s_m1l1(
"M1L1");
217 static const std::string s_m2l1(
"M2L1");
218 static const std::string s_m1l2(
"M1L2");
219 static const std::string s_m2l2(
"M2L2");
226 if ((sim_stationName==s_m1s1 || sim_stationName==s_m2s1) && sim_side==1){
227 double r_S1_A = std::sqrt(p.x()*p.x()+p.y()*p.y());
234 if ((sim_stationName==s_m1s1 || sim_stationName==s_m2s1) && sim_side==-1){
235 double r_S1_C = std::sqrt(p.x()*p.x()+p.y()*p.y());
245 if ((sim_stationName==s_m1s2 || sim_stationName==s_m2s2) && sim_side==1){
246 double r_S2_A = std::sqrt(p.x()*p.x()+p.y()*p.y());
253 if ((sim_stationName==s_m1s2 || sim_stationName==s_m2s2) && sim_side==-1){
254 double r_S2_C = std::sqrt(p.x()*p.x()+p.y()*p.y());
264 if ((sim_stationName==s_m1s2 || sim_stationName==s_m2s2 || sim_stationName==s_m1s1 || sim_stationName==s_m2s1) && sim_side==1){
265 double r_S_A = std::sqrt(p.x()*p.x()+p.y()*p.y());
271 if ((sim_stationName==s_m1s2 || sim_stationName==s_m2s2 || sim_stationName==s_m1s1 || sim_stationName==s_m2s1) && sim_side==-1){
272 double r_S_C = std::sqrt(p.x()*p.x()+p.y()*p.y());
282 double r_A = std::sqrt(p.x()*p.x()+p.y()*p.y());
289 double r_C = std::sqrt(p.x()*p.x()+p.y()*p.y());
300 if ((sim_stationName==s_m1l1 || sim_stationName==s_m2l1) && sim_side==1){
301 double r_L1_A = std::sqrt(p.x()*p.x()+p.y()*p.y());
308 if ((sim_stationName==s_m1l1 || sim_stationName==s_m2l1) && sim_side==-1){
309 double r_L1_C = std::sqrt(p.x()*p.x()+p.y()*p.y());
321 if ((sim_stationName==s_m1l2 || sim_stationName==s_m2l2) && sim_side==1){
322 double r_L2_A = std::sqrt(p.x()*p.x()+p.y()*p.y());
329 if ((sim_stationName==s_m1l2 || sim_stationName==s_m2l2) && sim_side==-1){
330 double r_L2_C = std::sqrt(p.x()*p.x()+p.y()*p.y());
340 if ((sim_stationName==s_m1l2 || sim_stationName==s_m2l2 || sim_stationName==s_m1l1 || sim_stationName==s_m2l1) && sim_side==1){
341 double r_L_A = std::sqrt(p.x()*p.x()+p.y()*p.y());
347 if ((sim_stationName==s_m1l2 || sim_stationName==s_m2l2 || sim_stationName==s_m1l1 || sim_stationName==s_m2l1) && sim_side==-1){
348 double r_L_C = std::sqrt(p.x()*p.x()+p.y()*p.y());
355 return StatusCode::SUCCESS;