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;
207 int simId = (*i_hit).MMId();
209 int sim_side = hitHelper->
GetSide(simId);
211 static const std::string s_m1s1(
"M1S1");
212 static const std::string s_m2s1(
"M2S1");
213 static const std::string s_m1s2(
"M1S2");
214 static const std::string s_m2s2(
"M2S2");
215 static const std::string s_m1l1(
"M1L1");
216 static const std::string s_m2l1(
"M2L1");
217 static const std::string s_m1l2(
"M1L2");
218 static const std::string s_m2l2(
"M2L2");
225 if ((sim_stationName==s_m1s1 || sim_stationName==s_m2s1) && sim_side==1){
226 double r_S1_A = std::sqrt(p.x()*p.x()+p.y()*p.y());
233 if ((sim_stationName==s_m1s1 || sim_stationName==s_m2s1) && sim_side==-1){
234 double r_S1_C = std::sqrt(p.x()*p.x()+p.y()*p.y());
244 if ((sim_stationName==s_m1s2 || sim_stationName==s_m2s2) && sim_side==1){
245 double r_S2_A = std::sqrt(p.x()*p.x()+p.y()*p.y());
252 if ((sim_stationName==s_m1s2 || sim_stationName==s_m2s2) && sim_side==-1){
253 double r_S2_C = std::sqrt(p.x()*p.x()+p.y()*p.y());
263 if ((sim_stationName==s_m1s2 || sim_stationName==s_m2s2 || sim_stationName==s_m1s1 || sim_stationName==s_m2s1) && sim_side==1){
264 double r_S_A = std::sqrt(p.x()*p.x()+p.y()*p.y());
270 if ((sim_stationName==s_m1s2 || sim_stationName==s_m2s2 || sim_stationName==s_m1s1 || sim_stationName==s_m2s1) && sim_side==-1){
271 double r_S_C = std::sqrt(p.x()*p.x()+p.y()*p.y());
281 double r_A = std::sqrt(p.x()*p.x()+p.y()*p.y());
288 double r_C = std::sqrt(p.x()*p.x()+p.y()*p.y());
299 if ((sim_stationName==s_m1l1 || sim_stationName==s_m2l1) && sim_side==1){
300 double r_L1_A = std::sqrt(p.x()*p.x()+p.y()*p.y());
307 if ((sim_stationName==s_m1l1 || sim_stationName==s_m2l1) && sim_side==-1){
308 double r_L1_C = std::sqrt(p.x()*p.x()+p.y()*p.y());
320 if ((sim_stationName==s_m1l2 || sim_stationName==s_m2l2) && sim_side==1){
321 double r_L2_A = std::sqrt(p.x()*p.x()+p.y()*p.y());
328 if ((sim_stationName==s_m1l2 || sim_stationName==s_m2l2) && sim_side==-1){
329 double r_L2_C = std::sqrt(p.x()*p.x()+p.y()*p.y());
339 if ((sim_stationName==s_m1l2 || sim_stationName==s_m2l2 || sim_stationName==s_m1l1 || sim_stationName==s_m2l1) && sim_side==1){
340 double r_L_A = std::sqrt(p.x()*p.x()+p.y()*p.y());
346 if ((sim_stationName==s_m1l2 || sim_stationName==s_m2l2 || sim_stationName==s_m1l1 || sim_stationName==s_m2l1) && sim_side==-1){
347 double r_L_C = std::sqrt(p.x()*p.x()+p.y()*p.y());
354 return StatusCode::SUCCESS;