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 2.0 2004/09/21 20:50:54 kusanagi Exp $ |
* $Id: CalibTrk1Reader.cpp,v 3.0 2005/03/04 15:54:11 kusanagi Exp $ |
4 |
* $Author: kusanagi $ |
* $Author: kusanagi $ |
5 |
* |
* |
6 |
* Implementation of the LogReader class. |
* Implementation of the LogReader class. |
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 2.0 2004/09/21 20:50:54 kusanagi Exp $\n"; |
"$Header: /home/cvsmanager/yoda/techmodel/CalibTrk1Reader.cpp,v 3.0 2005/03/04 15:54:11 kusanagi Exp $\n"; |
71 |
} |
} |
72 |
|
|
73 |
/** |
/** |
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; |
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 |
|
oss.str(""); |
96 |
|
oss << getenv("YODA_DATA") << "/" << time(NULL) << "trc1.dat"; |
97 |
|
FILE *pfile; |
98 |
|
pfile = fopen(oss.str().c_str(), "wb"); |
99 |
|
fwrite(subData, 1, dataLength, pfile); |
100 |
|
fclose(pfile); |
101 |
|
|
102 |
|
//Call to the FORTRAN routin that unpack tracker events |
103 |
|
trkcalibpkt_(&ERROR, (char*)oss.str().c_str()); |
104 |
|
|
105 |
|
calibTrk1->unpackError = ERROR; |
106 |
|
remove(oss.str().c_str()); |
107 |
|
if (ERROR != 0) { |
108 |
|
oss.str(""); |
109 |
|
oss << "Fortran77 function trkcalibpkt error code = " << ERROR; |
110 |
|
logger->warn(oss.str().c_str()); |
111 |
|
} |
112 |
|
|
113 |
|
remove((char*)oss.str().c_str()); |
114 |
|
|
115 |
//took the final CRC to compare it with the previous calculated CRC of the data |
//Store the unpacked data |
116 |
InputFile->read(eventCRC, sizeof(eventCRC)); |
memcpy(calibTrk1->DAQmode, trkcalib_.DAQmode, sizeof(calibTrk1->DAQmode)); |
117 |
readCRC = (((UINT16)(eventCRC[0]<<8))&0xFF00) + (((UINT16)eventCRC[1])&0x00FF); |
memcpy(calibTrk1->DSPnumber, trkcalib_.DSPnumber, sizeof(calibTrk1->DSPnumber)); |
118 |
|
memcpy(calibTrk1->calibnumber, trkcalib_.calibnumber, sizeof(calibTrk1->calibnumber)); |
119 |
if (subCRC == readCRC){ |
memcpy(calibTrk1->DSPnumber, trkcalib_.ncalib_event, sizeof(calibTrk1->ncalib_event)); |
120 |
int ERROR; |
memcpy(calibTrk1->ped_l1, trkcalib_.ped_l1, sizeof(calibTrk1->ped_l1)); |
121 |
//Scrivo un file temporaneo per passarlo alla routine |
memcpy(calibTrk1->ped_l2, trkcalib_.ped_l2, sizeof(calibTrk1->ped_l2)); |
122 |
//Speriamo di cambiare la routine per passargli un buffer..... |
memcpy(calibTrk1->ped_l3, trkcalib_.ped_l3, sizeof(calibTrk1->ped_l3)); |
123 |
DIR *dirp; |
memcpy(calibTrk1->sig_l1, trkcalib_.sig_l1, sizeof(calibTrk1->sig_l1)); |
124 |
std::string pathDir((char*)getenv("YODA_DATA")); |
memcpy(calibTrk1->sig_l2, trkcalib_.sig_l2, sizeof(calibTrk1->sig_l2)); |
125 |
pathDir = pathDir + "/todatemp.dat";; |
memcpy(calibTrk1->sig_l3, trkcalib_.sig_l3, sizeof(calibTrk1->sig_l3)); |
126 |
FILE *pfile; |
memcpy(calibTrk1->nbad_l1, trkcalib_.nbad_l1, sizeof(calibTrk1->nbad_l1)); |
127 |
pfile = fopen((char*)pathDir.c_str(), "wb"); |
memcpy(calibTrk1->nbad_l2, trkcalib_.nbad_l2, sizeof(calibTrk1->nbad_l2)); |
128 |
fwrite(subData, 1, dataLength, pfile); |
memcpy(calibTrk1->nbad_l3, trkcalib_.nbad_l3, sizeof(calibTrk1->nbad_l3)); |
129 |
fclose(pfile); |
memcpy(calibTrk1->cal_flag, trkcalib_.cal_flag, sizeof(calibTrk1->cal_flag)); |
|
|
|
|
//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)); |
|
130 |
// memcpy(calibTrk1->DSPbad_par, trkcalib_.DSPbad_par, sizeof(calibTrk1->DSPbad_par)); |
// memcpy(calibTrk1->DSPbad_par, trkcalib_.DSPbad_par, sizeof(calibTrk1->DSPbad_par)); |
131 |
// memcpy(calibTrk1->DSPped_par, trkcalib_.DSPped_par, sizeof(calibTrk1->DSPped_par)); |
// memcpy(calibTrk1->DSPped_par, trkcalib_.DSPped_par, sizeof(calibTrk1->DSPped_par)); |
132 |
// memcpy(calibTrk1->DSPsig_par, trkcalib_.DSPsig_par, sizeof(calibTrk1->DSPsig_par)); |
// memcpy(calibTrk1->DSPsig_par, trkcalib_.DSPsig_par, sizeof(calibTrk1->DSPsig_par)); |
133 |
//--------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 |
134 |
int tempBad_par[3072][6]; |
int tempBad_par[3072][6]; |
135 |
int tempPed_par[3072][6]; |
int tempPed_par[3072][6]; |
136 |
int tempSig_par[3072][6]; |
int tempSig_par[3072][6]; |
137 |
memcpy(calibTrk1->DSPbad_par,trkcalib_.DSPbad_par, sizeof(calibTrk1->DSPbad_par)); |
memcpy(calibTrk1->DSPbad_par,trkcalib_.DSPbad_par, sizeof(calibTrk1->DSPbad_par)); |
138 |
memcpy(calibTrk1->DSPped_par,trkcalib_.DSPped_par, sizeof(calibTrk1->DSPped_par)); |
memcpy(calibTrk1->DSPped_par,trkcalib_.DSPped_par, sizeof(calibTrk1->DSPped_par)); |
139 |
memcpy(calibTrk1->DSPsig_par,trkcalib_.DSPsig_par, sizeof(calibTrk1->DSPsig_par)); |
memcpy(calibTrk1->DSPsig_par,trkcalib_.DSPsig_par, sizeof(calibTrk1->DSPsig_par)); |
140 |
|
|
141 |
for (int i = 0; i < 6; i++){ |
for (int i = 0; i < 6; i++){ |
142 |
for (int j = 0; j < 3072; j++){ |
for (int j = 0; j < 3072; j++){ |
143 |
calibTrk1->DSPbad_par[i][j] = tempBad_par[j][i]; |
calibTrk1->DSPbad_par[i][j] = tempBad_par[j][i]; |
144 |
calibTrk1->DSPped_par[i][j] = tempPed_par[j][i]; |
calibTrk1->DSPped_par[i][j] = tempPed_par[j][i]; |
145 |
calibTrk1->DSPsig_par[i][j] = tempSig_par[j][i]; |
calibTrk1->DSPsig_par[i][j] = tempSig_par[j][i]; |
|
} |
|
146 |
} |
} |
|
//----------------------------------------------------------------------------------------- |
|
|
oss.flush(); |
|
|
oss << "Fortran77 function trkcalibpkt error code = " << ERROR; |
|
|
logger->warn(oss.str().c_str()); |
|
|
} else { |
|
|
logger->debug(_T("Wrong CRC for CalibTrk1 Packet ")); |
|
147 |
} |
} |
148 |
delete [] subData; |
//----------------------------------------------------------------------------------------- |
149 |
|
//delete [] subData; |
150 |
} |
} |
151 |
|
|
152 |
|
|