00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016 #include "pdscalibdata/CalibParsCenterV1.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 CalibParsCenterV1::CalibParsCenterV1 ()
00044 {
00045 double arr_x[] = {198., 198., 310., 98., 627., 628., 711., 498.,
00046 198., 198., 310., 98., 627., 628., 711., 498.,
00047 198., 198., 310., 98., 627., 628., 711., 498.,
00048 198., 198., 310., 98., 627., 628., 711., 498.};
00049
00050 double arr_y[] = {307., 95., 625., 625., 515., 727., 198., 199.,
00051 307., 95., 625., 625., 515., 727., 198., 199.,
00052 307., 95., 625., 625., 515., 727., 198., 199.,
00053 307., 95., 625., 625., 515., 727., 198., 199.};
00054
00055 int NPars = NQuad * NSect;
00056 size_t arr_size = sizeof( double ) * NPars;
00057 memcpy ( &m_center_x[0][0], arr_x, arr_size);
00058 memcpy ( &m_center_y[0][0], arr_y, arr_size);
00059 std::fill_n( &m_center_z[0][0], NPars, double(0));
00060 }
00061
00062
00063
00064 CalibParsCenterV1::CalibParsCenterV1( const std::vector<double> v_parameters )
00065 {
00066 if (v_parameters.size() != NUMBER_OF_PARAMETERS) {
00067 WithMsgLog("CalibParsCenterV1", error, str) {
00068 str << "Expected number of parameters is " << NUMBER_OF_PARAMETERS ;
00069 str << ", read from file " << v_parameters.size() ;
00070 str << ": check the file.\n" ;
00071 }
00072 abort();
00073 }
00074 size_t NPars = NQuad * NSect;
00075 size_t arr_size = sizeof( double ) * v_parameters.size()/3;
00076 memcpy( &m_center_x, &v_parameters[0], arr_size );
00077 memcpy( &m_center_y, &v_parameters[NPars], arr_size );
00078 memcpy( &m_center_z, &v_parameters[NPars*2], arr_size );
00079
00080 }
00081
00082 void CalibParsCenterV1::print()
00083 {
00084 cout << endl << "Center:" << endl;
00085 cout << "X:" << endl;
00086 for( int quad=0; quad<NQuad; ++quad ) {
00087 for( int sect=0; sect<NSect; ++sect ) { cout << " " << m_center_x[quad][sect]; }
00088 cout << endl;
00089 }
00090 cout << "Y:" << endl;
00091 for( int quad=0; quad<NQuad; ++quad ) {
00092 for( int sect=0; sect<NSect; ++sect ) { cout << " " << m_center_y[quad][sect]; }
00093 cout << endl;
00094 }
00095 cout << "Z:" << endl;
00096 for( int quad=0; quad<NQuad; ++quad ) {
00097 for( int sect=0; sect<NSect; ++sect ) { cout << " " << m_center_z[quad][sect]; }
00098 cout << endl;
00099 }
00100 }
00101
00102
00103
00104
00105
00106 CalibParsCenterV1::~CalibParsCenterV1 ()
00107 {
00108 }
00109
00110 }