PamOffLineSW/techmodel/CalibCalPedReader.cpp


extern "C" {

    //Struct per il passaggio di dati da e verso la chiamata fortran
    extern struct {
        int   iev;
        int   cstwerr[4];
        float cperror[4];
        float calped[4][11][96];
        float calgood[4][11][96];
        float calthr[4][11][6];
        float calrms[4][11][96];
        float calbase[4][11][6];
        float calvar[4][11][6];
        //float calpuls[4][11][96];
    } calib_; 
    //external declaration of the Fortran function 
    void calpedestal_(char*, long int*, int*);
}

#include "ReaderAlgorithms.h"
using namespace pamela::techmodel;


/**
 * Constructor. 
 */
CalibCalPedReader::CalibCalPedReader(void): 
  TechmodelAlgorithm(PacketType::Log, "TechmodelCalibCalPedReader") { 
  calibCalPed = new CalibCalPedEvent();
}


/**
 * Initialize the algorithm with a special run. This will initialize the
 * event reader routines for all packet types.
 */
void CalibCalPedReader::Init(PamelaRun *run) {
  run->WriteSubPacket(this, &calibCalPed, calibCalPed->Class());
}

/**
 * Unpack the CalibCalPed event from an input file.
 */
void CalibCalPedReader::PKT_RunEvent(char* packetData, long int dataLength) throw (CalibException){
    stringstream oss;
        string msg;
    int         ERROR;  
    calpedestal_((char*)packetData, &dataLength, &ERROR);
    calibCalPed->unpackError = ERROR;
    if (ERROR != 0) {
        char *errmsg;
        switch (ERROR){
            case 1: errmsg = "CALORIMETER NOT FOUND";
        }
        oss.str("");
        oss << "CalibCalPed: Fortran77 function calpedestal error code = " << ERROR
            <<  " " << errmsg;
        msg=oss.str();
        PamOffLineSW::mainLogUtil->logWarning(msg);     
     } //else {
      //Store the unpacked data
        calibCalPed->iev = calib_.iev;
        memcpy(calibCalPed->cstwerr, calib_.cstwerr, sizeof(calibCalPed->cstwerr));
        memcpy(calibCalPed->cperror, calib_.cperror, sizeof(calibCalPed->cperror));
    //--------have to invert array because of FORTRAN <-> C different management of the indexes
        float tempCalped[96][11][4];
        float tempCalgood[96][11][4];
        float tempCalthr[6][11][4];
        float tempCalrms[96][11][4];
        float tempCalbase[6][11][4];
        float tempCalvar[6][11][4];
        //float tempCalpuls[96][11][4];

        memcpy(tempCalped,  calib_.calped,  sizeof(tempCalped));
        memcpy(tempCalgood, calib_.calgood, sizeof(tempCalgood));
        memcpy(tempCalthr,  calib_.calthr,  sizeof(tempCalthr));
        memcpy(tempCalrms,  calib_.calrms,  sizeof(tempCalrms));
        memcpy(tempCalbase, calib_.calbase, sizeof(tempCalbase));
        memcpy(tempCalvar,  calib_.calvar,  sizeof(tempCalvar));
        //memcpy(tempCalpuls, calib_.calpuls, sizeof(tempCalpuls));

        for (int i = 0; i < 4; i++){
            for (int j = 0; j <11; j++){
                for (int z = 0; z < 96; z++){
                    calibCalPed->calped[i][j][z]    = tempCalped[z][j][i];
                    calibCalPed->calgood[i][j][z]   = tempCalgood[z][j][i];
                    calibCalPed->calrms[i][j][z]    = tempCalrms[z][j][i];
                    //calibCalPed->calpuls[i][j][z]   = tempCalpuls[z][j][i];
                }
            }
        }

        for (int i = 0; i < 4; i++){
            for (int j = 0; j <11; j++){
                for (int z = 0; z < 6; z++){
                    calibCalPed->calthr[i][j][z]    = tempCalthr[z][j][i];
                    calibCalPed->calbase[i][j][z]   = tempCalbase[z][j][i];
                    calibCalPed->calvar[i][j][z]    = tempCalvar[z][j][i];
                }
            }
        }
    //-----------------------------------------------------------------------------------------
    //}
        
        //TODO: check here ....
        if((calibCalPed->cperror[0])||(calibCalPed->cperror[1])||(calibCalPed->cperror[2])||(calibCalPed->cperror[3]))
        {
                oss.str("");
                oss <<  "CalibCalPed Problem: cperror[0] = " <<calibCalPed->cperror[0]<<" cperror[1] = " <<calibCalPed->cperror[1]<<" cperror[2] = " <<calibCalPed->cperror[2]<<" cperror[3] = " <<calibCalPed->cperror[3];
                msg=oss.str();
                PamOffLineSW::mainLogUtil->logError(msg);
                throw CalibException("Wrong Global CRC for CalibCalPed Packet ");
        }

}

/**
 * Get a string with the version info of the algorithm.
 */
std::string CalibCalPedReader::GetVersionInfo(void) const {
  return 
    "$Header: /repository/PamOffLineSW/techmodel/CalibCalPedReader.cpp,v 1.7 2008-03-21 16:18:46 messineo Exp $\n";
}
1
2	extern "C" {
3
4	//Struct per il passaggio di dati da e verso la chiamata fortran
5	extern struct {
6	int iev;
7	int cstwerr[4];
8	float cperror[4];
9	float calped[4][11][96];
10	float calgood[4][11][96];
11	float calthr[4][11][6];
12	float calrms[4][11][96];
13	float calbase[4][11][6];
14	float calvar[4][11][6];
15	//float calpuls[4][11][96];
16	} calib_;
17	//external declaration of the Fortran function
18	void calpedestal_(char, long int, int*);
19	}
20
21	#include "ReaderAlgorithms.h"
22	using namespace pamela::techmodel;
23
24
25	/**
26	* Constructor.
27	*/
28	CalibCalPedReader::CalibCalPedReader(void):
29	TechmodelAlgorithm(PacketType::Log, "TechmodelCalibCalPedReader") {
30	calibCalPed = new CalibCalPedEvent();
31	}
32
33
34
35	/**
36	* Initialize the algorithm with a special run. This will initialize the
37	* event reader routines for all packet types.
38	*/
39	void CalibCalPedReader::Init(PamelaRun *run) {
40	run->WriteSubPacket(this, &calibCalPed, calibCalPed->Class());
41	}
42
43	/**
44	* Unpack the CalibCalPed event from an input file.
45	*/
46	void CalibCalPedReader::PKT_RunEvent(char* packetData, long int dataLength) throw (CalibException){
47	stringstream oss;
48	string msg;
49	int ERROR;
50	calpedestal_((char*)packetData, &dataLength, &ERROR);
51	calibCalPed->unpackError = ERROR;
52	if (ERROR != 0) {
53	char *errmsg;
54	switch (ERROR){
55	case 1: errmsg = "CALORIMETER NOT FOUND";
56	}
57	oss.str("");
58	oss << "CalibCalPed: Fortran77 function calpedestal error code = " << ERROR
59	<< " " << errmsg;
60	msg=oss.str();
61	PamOffLineSW::mainLogUtil->logWarning(msg);
62	} //else {
63	//Store the unpacked data
64	calibCalPed->iev = calib_.iev;
65	memcpy(calibCalPed->cstwerr, calib_.cstwerr, sizeof(calibCalPed->cstwerr));
66	memcpy(calibCalPed->cperror, calib_.cperror, sizeof(calibCalPed->cperror));
67	//--------have to invert array because of FORTRAN <-> C different management of the indexes
68	float tempCalped[96][11][4];
69	float tempCalgood[96][11][4];
70	float tempCalthr[6][11][4];
71	float tempCalrms[96][11][4];
72	float tempCalbase[6][11][4];
73	float tempCalvar[6][11][4];
74	//float tempCalpuls[96][11][4];
75
76	memcpy(tempCalped, calib_.calped, sizeof(tempCalped));
77	memcpy(tempCalgood, calib_.calgood, sizeof(tempCalgood));
78	memcpy(tempCalthr, calib_.calthr, sizeof(tempCalthr));
79	memcpy(tempCalrms, calib_.calrms, sizeof(tempCalrms));
80	memcpy(tempCalbase, calib_.calbase, sizeof(tempCalbase));
81	memcpy(tempCalvar, calib_.calvar, sizeof(tempCalvar));
82	//memcpy(tempCalpuls, calib_.calpuls, sizeof(tempCalpuls));
83
84	for (int i = 0; i < 4; i++){
85	for (int j = 0; j <11; j++){
86	for (int z = 0; z < 96; z++){
87	calibCalPed->calped[i][j][z] = tempCalped[z][j][i];
88	calibCalPed->calgood[i][j][z] = tempCalgood[z][j][i];
89	calibCalPed->calrms[i][j][z] = tempCalrms[z][j][i];
90	//calibCalPed->calpuls[i][j][z] = tempCalpuls[z][j][i];
91	}
92	}
93	}
94
95	for (int i = 0; i < 4; i++){
96	for (int j = 0; j <11; j++){
97	for (int z = 0; z < 6; z++){
98	calibCalPed->calthr[i][j][z] = tempCalthr[z][j][i];
99	calibCalPed->calbase[i][j][z] = tempCalbase[z][j][i];
100	calibCalPed->calvar[i][j][z] = tempCalvar[z][j][i];
101	}
102	}
103	}
104	//-----------------------------------------------------------------------------------------
105	//}
106
107	//TODO: check here ....
108	if((calibCalPed->cperror[0])\|\|(calibCalPed->cperror[1])\|\|(calibCalPed->cperror[2])\|\|(calibCalPed->cperror[3]))
109	{
110	oss.str("");
111	oss << "CalibCalPed Problem: cperror[0] = " <<calibCalPed->cperror[0]<<" cperror[1] = " <<calibCalPed->cperror[1]<<" cperror[2] = " <<calibCalPed->cperror[2]<<" cperror[3] = " <<calibCalPed->cperror[3];
112	msg=oss.str();
113	PamOffLineSW::mainLogUtil->logError(msg);
114	throw CalibException("Wrong Global CRC for CalibCalPed Packet ");
115	}
116
117	}
118
119	/**
120	* Get a string with the version info of the algorithm.
121	*/
122	std::string CalibCalPedReader::GetVersionInfo(void) const {
123	return
124	"$Header: /repository/PamOffLineSW/techmodel/CalibCalPedReader.cpp,v 1.7 2008-03-21 16:18:46 messineo Exp $\n";
125	}