16 m_h_hit_x =
new TH1D(
"h_hit_x",
"hit_x", 100,-150.,150.);
20 m_h_hit_y =
new TH1D(
"h_hit_y",
"hit_y", 100,-150.,150.);
24 m_h_hit_z =
new TH1D(
"h_hit_z",
"hit_z", 100,-20000.,20000.);
28 m_h_xy =
new TH2D(
"h_xy",
"hit_xy", 100,-150.,150.,100,-150,150);
32 m_h_zr =
new TH2D(
"h_zr",
"hit_zr", 100,-20000.,20000.,100,0,250);
36 m_h_hit_post_x =
new TH1D(
"h_hit_post_x",
"hit_post_x", 100,-150.,150.);
40 m_h_hit_post_y =
new TH1D(
"h_hit_post_y",
"hit_post_y", 100,-150,150.);
44 m_h_hit_post_z =
new TH1D(
"h_hit_post_z",
"hit_post_z", 100,-15000,15000.);
48 m_h_hit_edep =
new TH1D(
"h_hit_edep",
"hit_edep", 100,0.,20.);
52 m_h_hit_pdgid =
new TH1D(
"h_hit_pdgid",
"hit_pdgid", 100,0.,7e6);
56 m_h_hit_pretime =
new TH1D(
"h_hit_pretime",
"hit_pretime", 100,0.,100.);
64 m_h_genvolume =
new TH1D(
"h_genvolume",
"genvolume", 20,0.,5.);
68 m_h_wavelength =
new TH1D(
"m_wavelength",
"wavelength", 150,0.,800.);
73 m_tree =
new TTree(
"Lucid",
"Lucid");
90 return StatusCode::SUCCESS;
110 const EventContext& ctx{Gaudi::Hive::currentContext()};
114 double x = i_hit->GetX();
115 double y = i_hit->GetY();
116 double z = i_hit->GetZ();
117 double r = sqrt(
x*
x+
y*
y);
134 m_hit_x->push_back(i_hit->GetX());
135 m_hit_y->push_back(i_hit->GetY());
136 m_hit_z->push_back(i_hit->GetZ());
150 return StatusCode::SUCCESS;