145 {
146
147
149
150 float multipl=1.E3;
151
153 double deltachisq = 0;
156
157
158 bool NewDisk = true, NewSector = true;
159 int previous_disk = -1, previous_sector = -1;
163 const Identifier Pixel_ModuleID = *wafer_it;
164 if(std::abs(
m_pixid->barrel_ec(Pixel_ModuleID)) == 2){
166 dparams,
167 deltachisq,
168 DOCA_Vector,
169 DOCA_Matrix));
170
171
172 if(previous_disk ==
m_pixid->layer_disk(Pixel_ModuleID) &&
174 NewSector=false;
175 else NewSector=true;
176 if(previous_disk ==
m_pixid->layer_disk(Pixel_ModuleID))
177 NewDisk=false;
178 else NewDisk=true;
181 else for (
unsigned int i=0;
i!=6;++
i)
m_h_PixEC_Align[i] -> Fill(multipl*dparams[i]);
182 previous_disk =
m_pixid->layer_disk(Pixel_ModuleID);
184
186 ATH_MSG_DEBUG(
"alignment parameters = (" << dparams[0] <<
"," << dparams[1] <<
"," << dparams[2] <<
"," << dparams[3] <<
"," << dparams[4] <<
"," << dparams[5] <<
")");
188 ATH_MSG_DEBUG(
"DOCA_Vector = (" << DOCA_Vector[0] <<
"," << DOCA_Vector[1] <<
"," << DOCA_Vector[2] <<
"," << DOCA_Vector[3] <<
"," << DOCA_Vector[4] <<
"," << DOCA_Vector[5] <<
")");
189 ATH_MSG_DEBUG(
"DOCA_Matrix = (" <<DOCA_Matrix(0,0) <<
"," << DOCA_Matrix(1,0) <<
"," << DOCA_Matrix(2,0) <<
"," << DOCA_Matrix(3,0) <<
"," << DOCA_Matrix(4,0) <<
"," << DOCA_Matrix(5,0) <<
")");
190 ATH_MSG_DEBUG(
"DOCA_Matrix = (" <<DOCA_Matrix(0,1) <<
"," << DOCA_Matrix(1,1) <<
"," << DOCA_Matrix(2,1) <<
"," << DOCA_Matrix(3,1) <<
"," << DOCA_Matrix(4,1) <<
"," << DOCA_Matrix(5,1) <<
")");
191 ATH_MSG_DEBUG(
"DOCA_Matrix = (" <<DOCA_Matrix(0,0) <<
"," << DOCA_Matrix(1,1) <<
"," << DOCA_Matrix(2,2) <<
"," << DOCA_Matrix(3,2) <<
"," << DOCA_Matrix(4,4) <<
"," << DOCA_Matrix(5,5) <<
")");
192
193 }
194 }
195 }
196
197
199 const Identifier Pixel_ModuleID = *wafer_it;
200 if(
m_pixid->barrel_ec(Pixel_ModuleID) == 0){
202 dparams,
203 deltachisq,
204 DOCA_Vector,
205 DOCA_Matrix));
206 ATH_MSG_DEBUG(
"Pixel Barrel ModuleID = " << Pixel_ModuleID);
207 ATH_MSG_DEBUG(
"alignment parameters = (" << dparams[0] <<
"," << dparams[1] <<
"," << dparams[2] <<
"," << dparams[3] <<
"," << dparams[4] <<
"," << dparams[5] <<
")");
209 ATH_MSG_DEBUG(
"DOCA_Vector = (" << DOCA_Vector[0] <<
"," << DOCA_Vector[1] <<
"," << DOCA_Vector[2] <<
"," << DOCA_Vector[3] <<
"," << DOCA_Vector[4] <<
"," << DOCA_Vector[5] <<
")");
210 ATH_MSG_DEBUG(
"DOCA_Matrix = (" <<DOCA_Matrix(0,0) <<
"," << DOCA_Matrix(1,0) <<
"," << DOCA_Matrix(2,0) <<
"," << DOCA_Matrix(3,0) <<
"," << DOCA_Matrix(4,0) <<
"," << DOCA_Matrix(5,0) <<
")");
211 ATH_MSG_DEBUG(
"DOCA_Matrix = (" <<DOCA_Matrix(0,1) <<
"," << DOCA_Matrix(1,1) <<
"," << DOCA_Matrix(2,1) <<
"," << DOCA_Matrix(3,1) <<
"," << DOCA_Matrix(4,1) <<
"," << DOCA_Matrix(5,1) <<
")");
212 ATH_MSG_DEBUG(
"DOCA_Matrix = (" <<DOCA_Matrix(0,0) <<
"," << DOCA_Matrix(1,1) <<
"," << DOCA_Matrix(2,2) <<
"," << DOCA_Matrix(3,2) <<
"," << DOCA_Matrix(4,4) <<
"," << DOCA_Matrix(5,5) <<
")");
213 break;
214 }
215 }
216
217
218
220 const Identifier SCT_ModuleID = *wafer_it;
221 if(std::abs(
m_sctid->barrel_ec(SCT_ModuleID)) == 2){
223 dparams,
224 deltachisq,
225 DOCA_Vector,
226 DOCA_Matrix));
228 ATH_MSG_DEBUG(
"alignment parameters = (" << dparams[0] <<
"," << dparams[1] <<
"," << dparams[2] <<
"," << dparams[3] <<
"," << dparams[4] <<
"," << dparams[5] <<
")");
230 ATH_MSG_DEBUG(
"DOCA_Vector = (" << DOCA_Vector[0] <<
"," << DOCA_Vector[1] <<
"," << DOCA_Vector[2] <<
"," << DOCA_Vector[3] <<
"," << DOCA_Vector[4] <<
"," << DOCA_Vector[5] <<
")");
231 ATH_MSG_DEBUG(
"DOCA_Matrix = (" <<DOCA_Matrix(0,0) <<
"," << DOCA_Matrix(1,0) <<
"," << DOCA_Matrix(2,0) <<
"," << DOCA_Matrix(3,0) <<
"," << DOCA_Matrix(4,0) <<
"," << DOCA_Matrix(5,0) <<
")");
232 ATH_MSG_DEBUG(
"DOCA_Matrix = (" <<DOCA_Matrix(0,1) <<
"," << DOCA_Matrix(1,1) <<
"," << DOCA_Matrix(2,1) <<
"," << DOCA_Matrix(3,1) <<
"," << DOCA_Matrix(4,1) <<
"," << DOCA_Matrix(5,1) <<
")");
233 ATH_MSG_DEBUG(
"DOCA_Matrix = (" <<DOCA_Matrix(0,0) <<
"," << DOCA_Matrix(1,1) <<
"," << DOCA_Matrix(2,2) <<
"," << DOCA_Matrix(3,2) <<
"," << DOCA_Matrix(4,4) <<
"," << DOCA_Matrix(5,5) <<
")");
234 break;
235 }
236 }
237
238
239
241 const Identifier SCT_ModuleID = *wafer_it;
242 if(
m_sctid->barrel_ec(SCT_ModuleID) == 0){
244 dparams,
245 deltachisq,
246 DOCA_Vector,
247 DOCA_Matrix));
249 ATH_MSG_DEBUG(
"alignment parameters = (" << dparams[0] <<
"," << dparams[1] <<
"," << dparams[2] <<
"," << dparams[3] <<
"," << dparams[4] <<
"," << dparams[5] <<
")");
251 ATH_MSG_DEBUG(
"DOCA_Vector = (" << DOCA_Vector[0] <<
"," << DOCA_Vector[1] <<
"," << DOCA_Vector[2] <<
"," << DOCA_Vector[3] <<
"," << DOCA_Vector[4] <<
"," << DOCA_Vector[5] <<
")");
252 ATH_MSG_DEBUG(
"DOCA_Matrix = (" <<DOCA_Matrix(0,0) <<
"," << DOCA_Matrix(1,0) <<
"," << DOCA_Matrix(2,0) <<
"," << DOCA_Matrix(3,0) <<
"," << DOCA_Matrix(4,0) <<
"," << DOCA_Matrix(5,0) <<
")");
253 ATH_MSG_DEBUG(
"DOCA_Matrix = (" <<DOCA_Matrix(0,1) <<
"," << DOCA_Matrix(1,1) <<
"," << DOCA_Matrix(2,1) <<
"," << DOCA_Matrix(3,1) <<
"," << DOCA_Matrix(4,1) <<
"," << DOCA_Matrix(5,1) <<
")");
254 ATH_MSG_DEBUG(
"DOCA_Matrix = (" <<DOCA_Matrix(0,0) <<
"," << DOCA_Matrix(1,1) <<
"," << DOCA_Matrix(2,2) <<
"," << DOCA_Matrix(3,2) <<
"," << DOCA_Matrix(4,4) <<
"," << DOCA_Matrix(5,5) <<
")");
255 break;
256 }
257 }
258
259 return StatusCode::SUCCESS;
260 }
std::vector< Identifier >::const_iterator const_id_iterator
std::vector< Identifier >::const_iterator const_id_iterator
virtual StatusCode finalize() override
Eigen::Matrix< double, Eigen::Dynamic, Eigen::Dynamic > MatrixX
Dynamic Matrix - dynamic allocation.