00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016 #include "pdscalibdata/CalibParsCenterCorrV1.h"
00017
00018
00019
00020
00021 #include <string.h>
00022
00023
00024
00025
00026 #include "MsgLogger/MsgLogger.h"
00027
00028
00029
00030
00031
00032 using namespace std;
00033
00034
00035
00036
00037
00038 namespace pdscalibdata {
00039
00040
00041
00042
00043 CalibParsCenterCorrV1::CalibParsCenterCorrV1 ()
00044 {
00045 int NPars = NQuad * NSect;
00046 std::fill_n(&m_center_corr_x[0][0], NPars, double(0));
00047 std::fill_n(&m_center_corr_y[0][0], NPars, double(0));
00048 std::fill_n(&m_center_corr_z[0][0], NPars, double(0));
00049 }
00050
00051
00052
00053 CalibParsCenterCorrV1::CalibParsCenterCorrV1( const std::vector<double> v_parameters )
00054 {
00055 if (v_parameters.size() != NUMBER_OF_PARAMETERS) {
00056 WithMsgLog("CalibParsCenterCorrV1", error, str) {
00057 str << "Expected number of parameters is " << NUMBER_OF_PARAMETERS ;
00058 str << ", read from file " << v_parameters.size() ;
00059 str << ": check the file.\n" ;
00060 }
00061 abort();
00062 }
00063 size_t NPars = NQuad * NSect;
00064 size_t arr_size = sizeof( double ) * v_parameters.size()/3;
00065 memcpy( &m_center_corr_x, &v_parameters[0], arr_size );
00066 memcpy( &m_center_corr_y, &v_parameters[NPars], arr_size );
00067 memcpy( &m_center_corr_z, &v_parameters[NPars*2], arr_size );
00068
00069 }
00070
00071 void CalibParsCenterCorrV1::print()
00072 {
00073 cout << endl << "Center correction:" << endl;
00074 cout << "X:" << endl;
00075 for( int quad=0; quad<NQuad; ++quad ) {
00076 for( int sect=0; sect<NSect; ++sect ) { cout << " " << m_center_corr_x[quad][sect]; }
00077 cout << endl;
00078 }
00079 cout << "Y:" << endl;
00080 for( int quad=0; quad<NQuad; ++quad ) {
00081 for( int sect=0; sect<NSect; ++sect ) { cout << " " << m_center_corr_y[quad][sect]; }
00082 cout << endl;
00083 }
00084 cout << "Z:" << endl;
00085 for( int quad=0; quad<NQuad; ++quad ) {
00086 for( int sect=0; sect<NSect; ++sect ) { cout << " " << m_center_corr_z[quad][sect]; }
00087 cout << endl;
00088 }
00089 }
00090
00091
00092
00093
00094
00095
00096 CalibParsCenterCorrV1::~CalibParsCenterCorrV1 ()
00097 {
00098 }
00099
00100 }