69 const std::string resultPath{
"NTUPLES/FILE1/InitialAlignment"};
70 NTupleFilePtr file1(
ntupleSvc(),
"/NTUPLES/FILE1");
74 nt =
ntupleSvc()->book(resultPath, CLID_ColumnWiseTuple,
"InitialAlignment");
98 std::map<Identifier, SurveyConstraintModule*, std::less<Identifier> > ModuleMap;
100 std::map<Identifier, SurveyConstraintModule*, std::less<Identifier> >
::iterator it;
101 for (it = ModuleMap.begin(); it != ModuleMap.end(); ++it) {
102 const Identifier& ModuleID = (it->second)->moduleID();
106 barrel_ec =
m_pixid->barrel_ec(ModuleID);
107 layer_disk =
m_pixid->layer_disk(ModuleID);
111 if (
m_sctid->side(ModuleID) != 0)
continue;
113 barrel_ec =
m_sctid->barrel_ec(ModuleID);
114 layer_disk =
m_sctid->layer_disk(ModuleID);
136 if (
sc.isFailure()) {
153 double deltachisq = 0;
158 bool NewDisk =
true, NewSector =
true;
159 int previous_disk = -1, previous_sector = -1;
164 if(std::abs(
m_pixid->barrel_ec(Pixel_ModuleID)) == 2){
172 if(previous_disk ==
m_pixid->layer_disk(Pixel_ModuleID) &&
176 if(previous_disk ==
m_pixid->layer_disk(Pixel_ModuleID))
179 if (NewDisk)
for (
unsigned int i=0;i!=6;++i)
m_h_PixEC_Align_Disk[i] -> Fill(multipl*dparams[i]);
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);
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) <<
")");
200 if(
m_pixid->barrel_ec(Pixel_ModuleID) == 0){
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) <<
")");
221 if(std::abs(
m_sctid->barrel_ec(SCT_ModuleID)) == 2){
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) <<
")");
242 if(
m_sctid->barrel_ec(SCT_ModuleID) == 0){
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) <<
")");
259 return StatusCode::SUCCESS;