![Logo](../../ATLAS-Logo-Square-Blue-RGB.png) |
ATLAS Offline Software
|
Go to the documentation of this file.
7 #include "GaudiKernel/MsgStream.h"
14 #include "GaudiKernel/SmartDataPtr.h"
15 #include "GaudiKernel/NTuple.h"
16 #include "GaudiKernel/INTupleSvc.h"
29 m_h_PixEC_Align_Disk{},
30 m_h_PixEC_Align_first{},
32 m_AlignResults_nModules{},
58 ATH_MSG_DEBUG(
"got ID helpers from detector store (relying on GeoModel to put them)" );
64 return StatusCode::SUCCESS;
71 NTupleFilePtr file1(
ntupleSvc(),
"/NTUPLES/FILE1");
73 NTuplePtr
nt(
ntupleSvc(),
"/NTUPLES/FILE1/InitialAlignment");
75 nt =
ntupleSvc()->book(
"/NTUPLES/FILE1/InitialAlignment", CLID_ColumnWiseTuple,
"InitialAlignment");
99 std::map<Identifier, SurveyConstraintModule*, std::less<Identifier> > ModuleMap;
101 std::map<Identifier, SurveyConstraintModule*, std::less<Identifier> >
::iterator it;
102 for (
it = ModuleMap.begin();
it != ModuleMap.end(); ++
it) {
136 sc =
ntupleSvc()->writeRecord(
"NTUPLES/FILE1/InitialAlignment");
137 if (
sc.isFailure()) {
154 double deltachisq = 0;
159 bool NewDisk =
true, NewSector =
true;
160 int previous_disk = -1, previous_sector = -1;
187 ATH_MSG_DEBUG(
"alignment parameters = (" << dparams[0] <<
"," << dparams[1] <<
"," << dparams[2] <<
"," << dparams[3] <<
"," << dparams[4] <<
"," << dparams[5] <<
")");
189 ATH_MSG_DEBUG(
"DOCA_Vector = (" << DOCA_Vector[0] <<
"," << DOCA_Vector[1] <<
"," << DOCA_Vector[2] <<
"," << DOCA_Vector[3] <<
"," << DOCA_Vector[4] <<
"," << DOCA_Vector[5] <<
")");
190 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) <<
")");
191 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) <<
")");
192 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) <<
")");
207 ATH_MSG_DEBUG(
"Pixel Barrel ModuleID = " << Pixel_ModuleID);
208 ATH_MSG_DEBUG(
"alignment parameters = (" << dparams[0] <<
"," << dparams[1] <<
"," << dparams[2] <<
"," << dparams[3] <<
"," << dparams[4] <<
"," << dparams[5] <<
")");
210 ATH_MSG_DEBUG(
"DOCA_Vector = (" << DOCA_Vector[0] <<
"," << DOCA_Vector[1] <<
"," << DOCA_Vector[2] <<
"," << DOCA_Vector[3] <<
"," << DOCA_Vector[4] <<
"," << DOCA_Vector[5] <<
")");
211 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) <<
")");
212 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) <<
")");
213 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) <<
")");
229 ATH_MSG_DEBUG(
"alignment parameters = (" << dparams[0] <<
"," << dparams[1] <<
"," << dparams[2] <<
"," << dparams[3] <<
"," << dparams[4] <<
"," << dparams[5] <<
")");
231 ATH_MSG_DEBUG(
"DOCA_Vector = (" << DOCA_Vector[0] <<
"," << DOCA_Vector[1] <<
"," << DOCA_Vector[2] <<
"," << DOCA_Vector[3] <<
"," << DOCA_Vector[4] <<
"," << DOCA_Vector[5] <<
")");
232 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) <<
")");
233 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) <<
")");
234 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) <<
")");
250 ATH_MSG_DEBUG(
"alignment parameters = (" << dparams[0] <<
"," << dparams[1] <<
"," << dparams[2] <<
"," << dparams[3] <<
"," << dparams[4] <<
"," << dparams[5] <<
")");
252 ATH_MSG_DEBUG(
"DOCA_Vector = (" << DOCA_Vector[0] <<
"," << DOCA_Vector[1] <<
"," << DOCA_Vector[2] <<
"," << DOCA_Vector[3] <<
"," << DOCA_Vector[4] <<
"," << DOCA_Vector[5] <<
")");
253 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) <<
")");
254 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) <<
")");
255 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) <<
")");
260 return StatusCode::SUCCESS;
275 return StatusCode::SUCCESS;
290 const static std::string
prefix =
"/stat/MisAlign/";
315 return StatusCode::SUCCESS;
def retrieve(aClass, aKey=None)
JetConstituentVector::iterator iterator
This is an Identifier helper class for the Pixel subdetector. This class is a factory for creating co...
NTuple::Item< long > m_AlignResults_Identifier_LayerDisc
SurveyConstraintTestAlg(const std::string &name, ISvcLocator *pSvcLocator)
NTuple::Item< long > m_AlignResults_Identifier_Eta
This is an Identifier helper class for the SCT subdetector. This class is a factory for creating comp...
std::vector< Identifier >::const_iterator const_id_iterator
Eigen::Matrix< double, Eigen::Dynamic, 1 > VectorX
Dynamic Vector - dynamic allocation.
NTuple::Item< double > m_AlignResults_y
Eigen::Matrix< double, Eigen::Dynamic, Eigen::Dynamic > MatrixX
Dynamic Matrix - dynamic allocation.
const_id_iterator wafer_end(void) const
Scalar phi() const
phi method
int barrel_ec(const Identifier &id) const
Values of different levels (failure returns 0)
Amg::VectorX DOCA_Vector() const
Scalar eta() const
pseudorapidity method
NTuple::Item< double > m_AlignResults_z
const_id_iterator wafer_begin(void) const
Iterators over full set of ids. Wafer iterator is sorted.
const_id_iterator wafer_begin(void) const
Iterators over full set of ids. Wafer iterator is sorted.
NTuple::Item< double > m_AlignResults_gamma
ISurveyConstraint * m_SurvConstr
TH1 * m_h_PixEC_Align_Disk[6]
std::vector< Identifier >::const_iterator const_id_iterator
IToolSvc * m_toolsvc
Pointer to tool service.
int barrel_ec(const Identifier &id) const
Values of different levels (failure returns 0)
int phi_module(const Identifier &id) const
NTuple::Item< double > m_AlignResults_x
const ServiceHandle< StoreGateSvc > & detStore() const
The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
void CreateMisAlignNtuple()
T * book(const std::string &n, const std::string &t, unsigned nbins)
::StatusCode StatusCode
StatusCode definition for legacy code.
NTuple::Item< double > m_AlignResults_beta
NTuple::Item< long > m_AlignResults_Identifier_Phi
int layer_disk(const Identifier &id) const
virtual int SectorNumber(int)=0
int eta_module(const Identifier &id) const
int layer_disk(const Identifier &id) const
TH1 * m_h_PixEC_Align_first[6]
NTuple::Item< long > m_AlignResults_Identifier_SCT
NTuple::Item< long > m_AlignResults_Identifier_BarrelEC
NTuple::Item< long > m_AlignResults_Identifier_ID
int eta_module(const Identifier &id) const
int m_AlignResults_nModules
int side(const Identifier &id) const
const_id_iterator wafer_end(void) const
def TH1F(name, title, nxbins, bins_par2, bins_par3=None, path='', **kwargs)
int phi_module(const Identifier &id) const
NTuple::Item< double > m_AlignResults_alpha