1 |
/** @file |
/** @file |
2 |
* $Source: /home/cvsmanager/yoda/techmodel/CalibTrk1Reader.cpp,v $ |
* $Source: /home/cvsmanager/yoda/techmodel/CalibTrk1Reader.cpp,v $ |
3 |
* $Id: CalibTrk1Reader.cpp,v 1.3 2004/07/17 20:03:38 kusanagi Exp $ |
* $Id: CalibTrk1Reader.cpp,v 2.1 2004/10/17 12:28:38 kusanagi Exp $ |
4 |
* $Author: kusanagi $ |
* $Author: kusanagi $ |
5 |
* |
* |
6 |
* Implementation of the LogReader class. |
* Implementation of the LogReader class. |
11 |
#define UINT unsigned int |
#define UINT unsigned int |
12 |
#define BYTE unsigned char |
#define BYTE unsigned char |
13 |
#include <string> |
#include <string> |
14 |
#include <log4cpp/Category.hh> |
#include <log4cxx/logger.h> |
15 |
extern "C" { |
extern "C" { |
16 |
#include "CRC.h" |
#include "CRC.h" |
17 |
//The int* return the error code |
//The int* return the error code |
38 |
float DSPsig_par[6][3072]; |
float DSPsig_par[6][3072]; |
39 |
float DSPbad_par[6][3072]; |
float DSPbad_par[6][3072]; |
40 |
} trkcalib_; |
} trkcalib_; |
|
|
|
41 |
#include <dirent.h> |
#include <dirent.h> |
42 |
} |
} |
43 |
|
|
51 |
using namespace pamela; |
using namespace pamela; |
52 |
using namespace pamela::techmodel; |
using namespace pamela::techmodel; |
53 |
|
|
54 |
static log4cpp::Category& cat = log4cpp::Category::getInstance("pamela.techmodel.CalibTrk1Reader"); |
static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.CalibTrk1Reader")); |
55 |
|
|
56 |
/** |
/** |
57 |
* Constructor. |
* Constructor. |
58 |
*/ |
*/ |
59 |
CalibTrk1Reader::CalibTrk1Reader(void): |
CalibTrk1Reader::CalibTrk1Reader(void): |
60 |
TechmodelAlgorithm(PacketType::CalibTrk1, "TechmodelCalibTrk1Reader") { |
TechmodelAlgorithm(PacketType::CalibTrk1, "TechmodelCalibTrk1Reader") { |
61 |
cat << log4cpp::Priority::DEBUG |
logger->debug(_T("Constructor")); |
|
<< "Constructor " |
|
|
<< "\n " << log4cpp::CategoryStream::ENDLINE; |
|
62 |
calibTrk1 = new CalibTrk1Event(); |
calibTrk1 = new CalibTrk1Event(); |
63 |
} |
} |
64 |
|
|
67 |
*/ |
*/ |
68 |
std::string CalibTrk1Reader::GetVersionInfo(void) const { |
std::string CalibTrk1Reader::GetVersionInfo(void) const { |
69 |
return |
return |
70 |
"$Header: /home/cvsmanager/yoda/techmodel/CalibTrk1Reader.cpp,v 1.3 2004/07/17 20:03:38 kusanagi Exp $\n"; |
"$Header: /home/cvsmanager/yoda/techmodel/CalibTrk1Reader.cpp,v 2.1 2004/10/17 12:28:38 kusanagi Exp $\n"; |
71 |
} |
} |
72 |
|
|
73 |
/** |
/** |
75 |
* event reader routines for all packet types. |
* event reader routines for all packet types. |
76 |
*/ |
*/ |
77 |
void CalibTrk1Reader::Init(PamelaRun *run) { |
void CalibTrk1Reader::Init(PamelaRun *run) { |
78 |
|
logger->debug(_T("Initialize")); |
79 |
SetInputStream(run); |
SetInputStream(run); |
80 |
run->WriteSubPacket(this, &calibTrk1, calibTrk1->Class()); |
run->WriteSubPacket(this, &calibTrk1, calibTrk1->Class()); |
81 |
} |
} |
84 |
* Unpack the CalibTrk event from an input file. |
* Unpack the CalibTrk event from an input file. |
85 |
*/ |
*/ |
86 |
//void CalibTrkReader::RunEvent(int EventNumber, long int length, char *subData) { |
//void CalibTrkReader::RunEvent(int EventNumber, long int length, char *subData) { |
87 |
void CalibTrk1Reader::RunEvent(int EventNumber, long int length) { |
void CalibTrk1Reader::RunEvent(int EventNumber, long int dataLength) throw (Exception){ |
88 |
|
std::stringstream oss; |
89 |
char *subData; |
char subData[dataLength]; |
90 |
char eventCRC[2]; |
int ERROR; |
|
UINT16 subCRC; //CRC of the data |
|
|
UINT16 readCRC; //CRC read from the end of the subpacket |
|
|
long int dataLength; |
|
|
|
|
|
//the 2 bytes subtracted belong to the final event CRC bytes |
|
|
dataLength = length - (long int)2; |
|
|
|
|
|
subData = new char[dataLength]; |
|
91 |
InputFile->read(subData, sizeof(unsigned char)*dataLength); |
InputFile->read(subData, sizeof(unsigned char)*dataLength); |
92 |
subCRC = CM_Compute_CRC16(0, (BYTE*)subData, dataLength); |
|
93 |
|
//Scrivo un file temporaneo per passarlo alla routine |
94 |
|
//Speriamo di cambiare la routine per passargli un buffer..... |
95 |
|
DIR *dirp; |
96 |
|
std::string pathDir((char*)getenv("YODA_DATA")); |
97 |
|
pathDir = pathDir + "/todatemp.dat";; |
98 |
|
FILE *pfile; |
99 |
|
pfile = fopen((char*)pathDir.c_str(), "wb"); |
100 |
|
fwrite(subData, 1, dataLength, pfile); |
101 |
|
fclose(pfile); |
102 |
|
|
103 |
|
//Call to the FORTRAN routin that unpack tracker events |
104 |
|
trkcalibpkt_(&ERROR, (char*)pathDir.c_str()); |
105 |
|
|
106 |
|
if (ERROR != 0) { |
107 |
|
oss.str(""); |
108 |
|
oss << "Fortran77 function trkcalibpkt error code = " << ERROR; |
109 |
|
logger->warn(oss.str().c_str()); |
110 |
|
} |
111 |
|
|
112 |
//took the final CRC to compare it with the previous calculated CRC of the data |
remove((char*)pathDir.c_str()); |
|
InputFile->read(eventCRC, sizeof(eventCRC)); |
|
|
readCRC = (((UINT16)(eventCRC[0]<<8))&0xFF00) + (((UINT16)eventCRC[1])&0x00FF); |
|
|
|
|
|
if (subCRC == readCRC){ |
|
|
int ERROR; |
|
|
//Scrivo un file temporaneo per passarlo alla routine |
|
|
//Speriamo di cambiare la routine per passargli un buffer..... |
|
|
DIR *dirp; |
|
|
std::string pathDir((char*)getenv("YODA_DATA")); |
|
|
pathDir = pathDir + "/todatemp.dat";; |
|
|
FILE *pfile; |
|
|
pfile = fopen((char*)pathDir.c_str(), "wb"); |
|
|
fwrite(subData, 1, dataLength, pfile); |
|
|
fclose(pfile); |
|
|
|
|
|
//Call to the FORTRAN routin that unpack tracker events |
|
|
trkcalibpkt_(&ERROR, (char*)pathDir.c_str()); |
|
|
|
|
|
remove((char*)pathDir.c_str()); |
|
|
|
|
|
//Store the unpacked data |
|
|
memcpy(calibTrk1->DAQmode, trkcalib_.DAQmode, sizeof(calibTrk1->DAQmode)); |
|
|
memcpy(calibTrk1->DSPnumber, trkcalib_.DSPnumber, sizeof(calibTrk1->DSPnumber)); |
|
|
memcpy(calibTrk1->calibnumber, trkcalib_.calibnumber, sizeof(calibTrk1->calibnumber)); |
|
|
memcpy(calibTrk1->DSPnumber, trkcalib_.ncalib_event, sizeof(calibTrk1->ncalib_event)); |
|
|
memcpy(calibTrk1->ped_l1, trkcalib_.ped_l1, sizeof(calibTrk1->ped_l1)); |
|
|
memcpy(calibTrk1->ped_l2, trkcalib_.ped_l2, sizeof(calibTrk1->ped_l2)); |
|
|
memcpy(calibTrk1->ped_l3, trkcalib_.ped_l3, sizeof(calibTrk1->ped_l3)); |
|
|
memcpy(calibTrk1->sig_l1, trkcalib_.sig_l1, sizeof(calibTrk1->sig_l1)); |
|
|
memcpy(calibTrk1->sig_l2, trkcalib_.sig_l2, sizeof(calibTrk1->sig_l2)); |
|
|
memcpy(calibTrk1->sig_l3, trkcalib_.sig_l3, sizeof(calibTrk1->sig_l3)); |
|
|
memcpy(calibTrk1->nbad_l1, trkcalib_.nbad_l1, sizeof(calibTrk1->nbad_l1)); |
|
|
memcpy(calibTrk1->nbad_l2, trkcalib_.nbad_l2, sizeof(calibTrk1->nbad_l2)); |
|
|
memcpy(calibTrk1->nbad_l3, trkcalib_.nbad_l3, sizeof(calibTrk1->nbad_l3)); |
|
|
memcpy(calibTrk1->cal_flag, trkcalib_.cal_flag, sizeof(calibTrk1->cal_flag)); |
|
113 |
|
|
114 |
|
//Store the unpacked data |
115 |
|
memcpy(calibTrk1->DAQmode, trkcalib_.DAQmode, sizeof(calibTrk1->DAQmode)); |
116 |
|
memcpy(calibTrk1->DSPnumber, trkcalib_.DSPnumber, sizeof(calibTrk1->DSPnumber)); |
117 |
|
memcpy(calibTrk1->calibnumber, trkcalib_.calibnumber, sizeof(calibTrk1->calibnumber)); |
118 |
|
memcpy(calibTrk1->DSPnumber, trkcalib_.ncalib_event, sizeof(calibTrk1->ncalib_event)); |
119 |
|
memcpy(calibTrk1->ped_l1, trkcalib_.ped_l1, sizeof(calibTrk1->ped_l1)); |
120 |
|
memcpy(calibTrk1->ped_l2, trkcalib_.ped_l2, sizeof(calibTrk1->ped_l2)); |
121 |
|
memcpy(calibTrk1->ped_l3, trkcalib_.ped_l3, sizeof(calibTrk1->ped_l3)); |
122 |
|
memcpy(calibTrk1->sig_l1, trkcalib_.sig_l1, sizeof(calibTrk1->sig_l1)); |
123 |
|
memcpy(calibTrk1->sig_l2, trkcalib_.sig_l2, sizeof(calibTrk1->sig_l2)); |
124 |
|
memcpy(calibTrk1->sig_l3, trkcalib_.sig_l3, sizeof(calibTrk1->sig_l3)); |
125 |
|
memcpy(calibTrk1->nbad_l1, trkcalib_.nbad_l1, sizeof(calibTrk1->nbad_l1)); |
126 |
|
memcpy(calibTrk1->nbad_l2, trkcalib_.nbad_l2, sizeof(calibTrk1->nbad_l2)); |
127 |
|
memcpy(calibTrk1->nbad_l3, trkcalib_.nbad_l3, sizeof(calibTrk1->nbad_l3)); |
128 |
|
memcpy(calibTrk1->cal_flag, trkcalib_.cal_flag, sizeof(calibTrk1->cal_flag)); |
129 |
|
// memcpy(calibTrk1->DSPbad_par, trkcalib_.DSPbad_par, sizeof(calibTrk1->DSPbad_par)); |
130 |
|
// memcpy(calibTrk1->DSPped_par, trkcalib_.DSPped_par, sizeof(calibTrk1->DSPped_par)); |
131 |
|
// memcpy(calibTrk1->DSPsig_par, trkcalib_.DSPsig_par, sizeof(calibTrk1->DSPsig_par)); |
132 |
//--------have to invert array because of FORTRAN <-> C different management of the indexes |
//--------have to invert array because of FORTRAN <-> C different management of the indexes |
133 |
int tempBad_par[3072][6]; |
int tempBad_par[3072][6]; |
134 |
int tempPed_par[3072][6]; |
int tempPed_par[3072][6]; |
135 |
int tempSig_par[3072][6]; |
int tempSig_par[3072][6]; |
136 |
memcpy(calibTrk1->DSPbad_par,trkcalib_.DSPbad_par, sizeof(calibTrk1->DSPbad_par)); |
memcpy(calibTrk1->DSPbad_par,trkcalib_.DSPbad_par, sizeof(calibTrk1->DSPbad_par)); |
137 |
memcpy(calibTrk1->DSPped_par,trkcalib_.DSPped_par, sizeof(calibTrk1->DSPped_par)); |
memcpy(calibTrk1->DSPped_par,trkcalib_.DSPped_par, sizeof(calibTrk1->DSPped_par)); |
138 |
memcpy(calibTrk1->DSPsig_par,trkcalib_.DSPsig_par, sizeof(calibTrk1->DSPsig_par)); |
memcpy(calibTrk1->DSPsig_par,trkcalib_.DSPsig_par, sizeof(calibTrk1->DSPsig_par)); |
139 |
|
|
140 |
for (int i = 0; i < 6; i++){ |
for (int i = 0; i < 6; i++){ |
141 |
for (int j = 0; j < 3072; j++){ |
for (int j = 0; j < 3072; j++){ |
142 |
calibTrk1->DSPbad_par[i][j] = tempBad_par[j][i]; |
calibTrk1->DSPbad_par[i][j] = tempBad_par[j][i]; |
143 |
calibTrk1->DSPped_par[i][j] = tempPed_par[j][i]; |
calibTrk1->DSPped_par[i][j] = tempPed_par[j][i]; |
144 |
calibTrk1->DSPsig_par[i][j] = tempSig_par[j][i]; |
calibTrk1->DSPsig_par[i][j] = tempSig_par[j][i]; |
|
} |
|
145 |
} |
} |
|
//----------------------------------------------------------------------------------------- |
|
|
|
|
|
|
|
|
cat << log4cpp::Priority::ERROR |
|
|
<< "Fortran77 function trkcalibpkt error code = " << ERROR |
|
|
<< "\n " << log4cpp::CategoryStream::ENDLINE; |
|
|
free(subData); |
|
|
} else { |
|
|
cat << log4cpp::Priority::ERROR |
|
|
<< "Wrong CRC for CalibTrk1 Packet " |
|
|
<< "\n " << log4cpp::CategoryStream::ENDLINE; |
|
146 |
} |
} |
147 |
|
//----------------------------------------------------------------------------------------- |
148 |
|
delete [] subData; |
149 |
} |
} |
150 |
|
|
151 |
|
|