00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016 #include "pdscalibdata/CsPadCenterGlobalV1.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 CsPadCenterGlobalV1::CsPadCenterGlobalV1 ()
00044 {
00045 double arr_x[] = { 473.38, 685.26, 155.01, 154.08, 266.81, 53.95, 583.04, 582.15,
00046 989.30, 987.12, 1096.93, 884.11, 1413.16, 1414.94, 1500.83, 1288.02,
00047 1142.59, 930.23, 1459.44, 1460.67, 1347.57, 1559.93, 1032.27, 1033.44,
00048 626.78, 627.42, 516.03, 729.15, 198.28, 198.01, 115.31, 327.66};
00049
00050 double arr_y[] = {1028.07, 1026.28, 1139.46, 926.91, 1456.78, 1457.35, 1539.71, 1327.89,
00051 1180.51, 967.36, 1497.74, 1498.54, 1385.08, 1598.19, 1069.65, 1069.93,
00052 664.89, 666.83, 553.60, 765.91, 237.53, 236.06, 152.17, 365.47,
00053 510.38, 722.95, 193.33, 193.41, 308.04, 95.25, 625.28, 624.14};
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 CsPadCenterGlobalV1::CsPadCenterGlobalV1( const std::vector<double> v_parameters )
00065 {
00066 if (v_parameters.size() != NUMBER_OF_PARAMETERS) {
00067 WithMsgLog("CsPadCenterGlobalV1", 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 CsPadCenterGlobalV1::print()
00083 {
00084 cout << endl << "Center global:" << 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 CsPadCenterGlobalV1::~CsPadCenterGlobalV1 ()
00107 {
00108 }
00109
00110 }