126 std::string
names[162] = {
"eta",
"pt",
"phi",
"pos_x",
"pos_y",
"pos_z",
127 "emb0_cell",
"emb1_cell",
"emb2_cell",
"emb3_cell",
"emec0_cell",
"emec1_cell",
"emec2_cell",
"emec3_cell",
128 "hec0_cell",
"hec1_cell",
"hec2_cell",
"hec3_cell",
"fc1_cell",
"fc2_cell",
"fc3_cell",
129 "emb0_hits",
"emb1_hits",
"emb2_hits",
"emb3_hits",
"emec0_hits",
"emec1_hits",
"emec2_hits",
"emec3_hits",
130 "hec0_hits",
"hec1_hits",
"hec2_hits",
"hec3_hits",
"fc1_hits",
"fc2_hits",
"fc3_hits",
131 "emb0_sumE",
"emb1_sumE",
"emb2_sumE",
"emb3_sumE",
"emec0_sumE",
"emec1_sumE",
"emec2_sumE",
"emec3_sumE",
132 "hec0_sumE",
"hec1_sumE",
"hec2_sumE",
"hec3_sumE",
"fc1_sumE",
"fc2_sumE",
"fc3_sumE",
133 "emb0_dPhi",
"emb1_dPhi",
"emb2_dPhi",
"emb3_dPhi",
"emec0_dPhi",
"emec1_dPhi",
"emec2_dPhi",
"emec3_dPhi",
134 "hec0_dPhi",
"hec1_dPhi",
"hec2_dPhi",
"hec3_dPhi",
"fc1_dX",
"fc2_dX",
"fc3_dX",
135 "emb0_sPhi",
"emb1_sPhi",
"emb2_sPhi",
"emb3_sPhi",
"emec0_sPhi",
"emec1_sPhi",
"emec2_sPhi",
"emec3_sPhi",
136 "hec0_sPhi",
"hec1_sPhi",
"hec2_sPhi",
"hec3_sPhi",
"fc1_sX",
"fc2_sX",
"fc3_sX",
137 "emb0_dEta",
"emb1_dEta",
"emb2_dEta",
"emb3_dEta",
"emec0_dEta",
"emec1_dEta",
"emec2_dEta",
"emec3_dEta",
138 "hec0_dEta",
"hec1_dEta",
"hec2_dEta",
"hec3_dEta",
"fc1_dY",
"fc2_dY",
"fc3_dY",
139 "emb0_sEta",
"emb1_sEta",
"emb2_sEta",
"emb3_sEta",
"emec0_sEta",
"emec1_sEta",
"emec2_sEta",
"emec3_sEta",
140 "hec0_sEta",
"hec1_sEta",
"hec2_sEta",
"hec3_sEta",
"fc1_sY",
"fc2_sY",
"fc3_sY",
141 "emb0_time",
"emb1_time",
"emb2_time",
"emb3_time",
"emec0_time",
"emec1_time",
"emec2_time",
"emec3_time",
142 "hec0_time",
"hec1_time",
"hec2_time",
"hec3_time",
"fc1_time",
"fc2_time",
"fc3_time",
143 "emb0_widthX",
"emb1_widthX",
"emb2_widthX",
"emb3_widthX",
"emec0_widthX",
"emec1_widthX",
"emec2_widthX",
"emec3_widthX",
144 "hec0_widthX",
"hec1_widthX",
"hec2_widthX",
"hec3_widthX",
"fc1_widthX",
"fc2_widthX",
"fc3_widthX",
145 "emb0_widthY",
"emb1_widthY",
"emb2_widthY",
"emb3_widthY",
"emec0_widthY",
"emec1_widthY",
"emec2_widthY",
"emec3_widthY",
146 "hec0_widthY",
"hec1_widthY",
"hec2_widthY",
"hec3_widthY",
"fc1_widthY",
"fc2_widthY",
"fc3_widthY",
147 "cpuTime",
"Energy",
"PDG_ID",
"RunNo",
"EventNo",
"E_Dep" };
149 double lim[162][2] = { {0,5}, {0,100}, {-4,4}, {-1600,1600}, {-1600,1600}, {-4000,4000},
150 {0,0.25}, {0,3}, {0,6}, {0,0.1}, {0,0.25}, {0,2}, {0,3}, {0,0.1}, {0,1}, {0,10}, {0,10}, {0,10}, {0,0.1}, {0,0.1}, {0,0.1},
151 {0,100}, {0,600}, {0,600}, {0,50}, {0,50}, {0,400}, {0,500}, {0,200}, {0,100}, {0,1000}, {0,1000}, {0,100}, {0,150}, {0,50}, {0,10},
152 {0,0.4}, {0,5}, {0,10}, {0,0.2}, {0,0.1}, {0,3}, {0,5}, {0,0.2}, {0,1}, {0,10}, {0,10}, {0,0.1}, {0,1}, {0,0.1}, {0,0.1},
153 {-500,500}, {-100,100}, {-15,15}, {-200,200}, {-3000,3000}, {-60,60}, {-25,25}, {-200,200}, {-50,50}, {-50,50}, {-50,50}, {-50,50}, {-60,60}, {-500,500}, {-200,200},
154 {0,1000}, {0,500}, {0,200}, {0,500}, {0,2000}, {0,250}, {0,300}, {0,500}, {0,200}, {0,200}, {0,200}, {0,200}, {0,100}, {0,100}, {0,50},
155 {-150,150}, {-15,15}, {-20,20}, {-200,200}, {-0,2500}, {-50,20}, {-15,15}, {-150,150}, {-50,50}, {-50,50}, {-50,50}, {-50,50}, {-60,60}, {-500,500}, {-200,200},
156 {0,500}, {0,100}, {0,100}, {0,400}, {0,1000}, {0,150}, {0,100}, {0,400}, {0,200}, {0,200}, {0,200}, {0,200}, {0,60}, {0,100}, {0,50},
157 {0,750}, {0,25}, {0,20}, {0,1000}, {0,10000}, {0,40}, {0,30}, {0,1000}, {0,1000}, {0,100}, {0,100}, {0,200}, {0,10}, {0,500}, {0,100},
158 {-150,150}, {-15,15}, {-20,20}, {-200,200}, {-0,2500}, {-50,20}, {-15,15}, {-150,150}, {-50,50}, {-50,50}, {-50,50}, {-50,50}, {-60,60}, {-500,500}, {-200,200},
159 {-150,150}, {-15,15}, {-20,20}, {-200,200}, {-0,2500}, {-50,20}, {-15,15}, {-150,150}, {-50,50}, {-50,50}, {-50,50}, {-50,50}, {-60,60}, {-500,500}, {-200,200},
160 {0,50}, {0,100}, {-25,25}, {0,10}, {0,1000}, {0,10}};
179 for (
int i=0;
i<162;
i++){
184 ATH_MSG_WARNING(
"Failed to register historam " <<
names[
i] <<
". Not sure what will happen now..." );
190 if (!
file)
throw std::runtime_error (
"Ntuple MGR not open");
191 NTuple::Directory *
col=
ntupleSvc()->createDirectory(
"/NTUPLES/FILE/COL");
192 NTuplePtr
nt(
ntupleSvc(),
"/NTUPLES/FILE/COL/SingleTrackValidation");
193 if (!
nt)
nt=
ntupleSvc()->book(
col, 1, CLID_ColumnWiseTuple,
"SingleTrackValidation");
195 if (
nt->addItem(
"Eta",
m_c->
eta ).isFailure() ||
196 nt->addItem(
"Pt",
m_c->
pt ).isFailure() ||
197 nt->addItem(
"BarrelX",
m_c->
x ).isFailure() ||
198 nt->addItem(
"BarrelY",
m_c->
y ).isFailure() ||
199 nt->addItem(
"BarrelZ",
m_c->
z ).isFailure() ||
200 nt->addItem(
"Phi",
m_c->
phi ).isFailure() ){
201 ATH_MSG_WARNING(
"Registration of some of the ntuple branches failed. No idea what will happen next..." );
210 for (
int i=0;
i<15;
i++){
211 if (
i<12) sprintf(
title,
"S%i_C00",
i);
212 else sprintf(
title,
"FC%i_C00",
i-11);
214 if (
i<12) sprintf(
title,
"S%i_SumE",
i);
215 else sprintf(
title,
"FC%i_SumE",
i-11);
217 if (
i<12) sprintf(
title,
"S%i_Hits",
i);
218 else sprintf(
title,
"FC%i_Hits",
i-11);
220 if (
i<12) sprintf(
title,
"S%i_DeltaPhi",
i);
221 else sprintf(
title,
"FC%i_DeltaX",
i-11);
223 if (
i<12) sprintf(
title,
"S%i_SigmaPhi",
i);
224 else sprintf(
title,
"FC%i_SigmaX",
i-11);
226 if (
i<12) sprintf(
title,
"S%i_DeltaEta",
i);
227 else sprintf(
title,
"FC%i_DeltaY",
i-11);
229 if (
i<12) sprintf(
title,
"S%i_SigmaEta",
i);
230 else sprintf(
title,
"FC%i_SigmaY",
i-11);
232 if (
i<12) sprintf(
title,
"S%i_Time",
i);
233 else sprintf(
title,
"FC%i_Time",
i-11);
235 if (
i<12) sprintf(
title,
"S%i_WidthX",
i);
236 else sprintf(
title,
"FC%i_WidthX",
i-11);
238 if (
i<12) sprintf(
title,
"S%i_WidthY",
i);
239 else sprintf(
title,
"FC%i_WidthY",
i-11);
244 nt->addItem(
"TrackEnergy" ,
m_c->
Energy ).isFailure() ||
245 nt->addItem(
"ParticleID" ,
m_c->
PDG ).isFailure() ||
246 nt->addItem(
"Run#" ,
m_c->
RunNo ).isFailure() ||
249 ATH_MSG_WARNING(
"Registration of some of the ntuple branches failed. No idea what will happen next..." );
259 return StatusCode::SUCCESS;