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 3.0 2005/03/04 15:54:11 kusanagi Exp $ |
* $Id: CalibTrk1Reader.cpp,v 6.1 2006/05/30 19:10:02 kusanagi Exp $ |
4 |
* $Author: kusanagi $ |
* $Author: kusanagi $ |
5 |
* |
* |
6 |
* Implementation of the LogReader class. |
* Implementation of the LogReader class. |
8 |
* Control the CRC for the entire data Packet not just for single records |
* Control the CRC for the entire data Packet not just for single records |
9 |
*/ |
*/ |
10 |
|
|
11 |
#define UINT unsigned int |
#include "ReaderAlgorithms.h" |
12 |
#define BYTE unsigned char |
|
|
#include <string> |
|
|
#include <log4cxx/logger.h> |
|
13 |
extern "C" { |
extern "C" { |
14 |
#include "CRC.h" |
#include "CRC.h" |
15 |
//The int* return the error code |
extern void trkcalibpkt_(int*, unsigned char[], long int*, int*); //(*) |
|
//The char* is the path to the temporary file |
|
|
extern void trkcalibpkt_(int*, char*); |
|
16 |
|
|
17 |
//Struct per il passaggio di dati da e verso la chiamata fortran |
//Struct per il passaggio di dati da e verso la chiamata fortran |
18 |
extern struct { |
extern struct { |
19 |
int DAQmode[6]; |
int good0; |
20 |
int DSPnumber[6]; |
int DAQmode[6]; |
21 |
int calibnumber[6]; |
int DSPnumber[6]; |
22 |
int ncalib_event[6]; |
int calibnumber[6]; |
23 |
int ped_l1[6]; |
int ncalib_event[6]; |
24 |
int ped_l2[6]; |
int ped_l1[6]; |
25 |
int ped_l3[6]; |
int ped_l2[6]; |
26 |
int sig_l1[6]; |
int ped_l3[6]; |
27 |
int sig_l2[6]; |
int sig_l1[6]; |
28 |
int sig_l3[6]; |
int sig_l2[6]; |
29 |
int nbad_l1[6]; |
int sig_l3[6]; |
30 |
int nbad_l2[6]; |
int nbad_l1[6]; |
31 |
int nbad_l3[6]; |
int nbad_l2[6]; |
32 |
int cal_flag[6]; |
int nbad_l3[6]; |
33 |
int DSPped_par[6][3072]; |
int cal_flag[6]; |
34 |
float DSPsig_par[6][3072]; |
float DSPped_par[3072][6]; |
35 |
float DSPbad_par[6][3072]; |
float DSPsig_par[3072][6]; |
36 |
|
int DSPbad_par[3072][6]; |
37 |
|
int crc_cal[3][6]; |
38 |
|
int crc_hcal[6]; |
39 |
} trkcalib_; |
} trkcalib_; |
40 |
#include <dirent.h> |
#include <dirent.h> |
41 |
} |
} |
42 |
|
|
|
#include <fstream> |
|
|
#include "stdio.h" |
|
|
#include "ReaderAlgorithms.h" |
|
|
#include "event/PamelaRun.h" |
|
|
|
|
|
#include "event/CalibTrk1Event.h" |
|
|
|
|
|
using namespace pamela; |
|
43 |
using namespace pamela::techmodel; |
using namespace pamela::techmodel; |
44 |
|
|
45 |
static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.CalibTrk1Reader")); |
static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.CalibTrk1Reader")); |
58 |
*/ |
*/ |
59 |
std::string CalibTrk1Reader::GetVersionInfo(void) const { |
std::string CalibTrk1Reader::GetVersionInfo(void) const { |
60 |
return |
return |
61 |
"$Header: /home/cvsmanager/yoda/techmodel/CalibTrk1Reader.cpp,v 3.0 2005/03/04 15:54:11 kusanagi Exp $\n"; |
"$Header: /home/cvsmanager/yoda/techmodel/CalibTrk1Reader.cpp,v 6.1 2006/05/30 19:10:02 kusanagi Exp $\n"; |
62 |
} |
} |
63 |
|
|
64 |
/** |
/** |
80 |
char subData[dataLength]; |
char subData[dataLength]; |
81 |
int ERROR; |
int ERROR; |
82 |
InputFile->read(subData, sizeof(unsigned char)*dataLength); |
InputFile->read(subData, sizeof(unsigned char)*dataLength); |
|
|
|
|
//Scrivo un file temporaneo per passarlo alla routine |
|
|
//Speriamo di cambiare la routine per passargli un buffer..... |
|
|
oss.str(""); |
|
|
oss << getenv("YODA_DATA") << "/" << time(NULL) << "trc1.dat"; |
|
|
FILE *pfile; |
|
|
pfile = fopen(oss.str().c_str(), "wb"); |
|
|
fwrite(subData, 1, dataLength, pfile); |
|
|
fclose(pfile); |
|
83 |
|
|
84 |
|
char *data = new char[dataLength]; |
85 |
|
memcpy(data, subData, dataLength); |
86 |
|
int curpos=1; |
87 |
|
|
88 |
//Call to the FORTRAN routin that unpack tracker events |
//Call to the FORTRAN routin that unpack tracker events |
89 |
trkcalibpkt_(&ERROR, (char*)oss.str().c_str()); |
trkcalibpkt_( &ERROR,(unsigned char*)data, &dataLength, &curpos); |
90 |
|
|
91 |
calibTrk1->unpackError = ERROR; |
calibTrk1->unpackError = ERROR; |
|
remove(oss.str().c_str()); |
|
92 |
if (ERROR != 0) { |
if (ERROR != 0) { |
93 |
oss.str(""); |
oss.str(""); |
94 |
oss << "Fortran77 function trkcalibpkt error code = " << ERROR; |
oss << "Fortran77 function trkcalibpkt error code = " << ERROR; |
95 |
logger->warn(oss.str().c_str()); |
logger->warn(oss.str().c_str()); |
96 |
} |
} |
97 |
|
|
|
remove((char*)oss.str().c_str()); |
|
|
|
|
98 |
//Store the unpacked data |
//Store the unpacked data |
99 |
|
calibTrk1->good0 = trkcalib_.good0; |
100 |
memcpy(calibTrk1->DAQmode, trkcalib_.DAQmode, sizeof(calibTrk1->DAQmode)); |
memcpy(calibTrk1->DAQmode, trkcalib_.DAQmode, sizeof(calibTrk1->DAQmode)); |
101 |
memcpy(calibTrk1->DSPnumber, trkcalib_.DSPnumber, sizeof(calibTrk1->DSPnumber)); |
memcpy(calibTrk1->DSPnumber, trkcalib_.DSPnumber, sizeof(calibTrk1->DSPnumber)); |
102 |
memcpy(calibTrk1->calibnumber, trkcalib_.calibnumber, sizeof(calibTrk1->calibnumber)); |
memcpy(calibTrk1->calibnumber, trkcalib_.calibnumber, sizeof(calibTrk1->calibnumber)); |
103 |
memcpy(calibTrk1->DSPnumber, trkcalib_.ncalib_event, sizeof(calibTrk1->ncalib_event)); |
memcpy(calibTrk1->ncalib_event, trkcalib_.ncalib_event, sizeof(calibTrk1->ncalib_event)); |
104 |
memcpy(calibTrk1->ped_l1, trkcalib_.ped_l1, sizeof(calibTrk1->ped_l1)); |
memcpy(calibTrk1->ped_l1, trkcalib_.ped_l1, sizeof(calibTrk1->ped_l1)); |
105 |
memcpy(calibTrk1->ped_l2, trkcalib_.ped_l2, sizeof(calibTrk1->ped_l2)); |
memcpy(calibTrk1->ped_l2, trkcalib_.ped_l2, sizeof(calibTrk1->ped_l2)); |
106 |
memcpy(calibTrk1->ped_l3, trkcalib_.ped_l3, sizeof(calibTrk1->ped_l3)); |
memcpy(calibTrk1->ped_l3, trkcalib_.ped_l3, sizeof(calibTrk1->ped_l3)); |
111 |
memcpy(calibTrk1->nbad_l2, trkcalib_.nbad_l2, sizeof(calibTrk1->nbad_l2)); |
memcpy(calibTrk1->nbad_l2, trkcalib_.nbad_l2, sizeof(calibTrk1->nbad_l2)); |
112 |
memcpy(calibTrk1->nbad_l3, trkcalib_.nbad_l3, sizeof(calibTrk1->nbad_l3)); |
memcpy(calibTrk1->nbad_l3, trkcalib_.nbad_l3, sizeof(calibTrk1->nbad_l3)); |
113 |
memcpy(calibTrk1->cal_flag, trkcalib_.cal_flag, sizeof(calibTrk1->cal_flag)); |
memcpy(calibTrk1->cal_flag, trkcalib_.cal_flag, sizeof(calibTrk1->cal_flag)); |
114 |
// memcpy(calibTrk1->DSPbad_par, trkcalib_.DSPbad_par, sizeof(calibTrk1->DSPbad_par)); |
|
115 |
// memcpy(calibTrk1->DSPped_par, trkcalib_.DSPped_par, sizeof(calibTrk1->DSPped_par)); |
int tempBad_par[3072][6]; |
116 |
// memcpy(calibTrk1->DSPsig_par, trkcalib_.DSPsig_par, sizeof(calibTrk1->DSPsig_par)); |
float tempPed_par[3072][6]; |
117 |
//--------have to invert array because of FORTRAN <-> C different management of the indexes |
float tempSig_par[3072][6]; |
118 |
int tempBad_par[3072][6]; |
int tempCrc_cal[3][6]; |
119 |
int tempPed_par[3072][6]; |
memcpy(tempBad_par,trkcalib_.DSPbad_par, sizeof(tempBad_par)); |
120 |
int tempSig_par[3072][6]; |
memcpy(tempPed_par,trkcalib_.DSPped_par, sizeof(tempPed_par)); |
121 |
memcpy(calibTrk1->DSPbad_par,trkcalib_.DSPbad_par, sizeof(calibTrk1->DSPbad_par)); |
memcpy(tempSig_par,trkcalib_.DSPsig_par, sizeof(tempSig_par)); |
122 |
memcpy(calibTrk1->DSPped_par,trkcalib_.DSPped_par, sizeof(calibTrk1->DSPped_par)); |
memcpy(tempCrc_cal,trkcalib_.crc_cal, sizeof(tempCrc_cal)); |
|
memcpy(calibTrk1->DSPsig_par,trkcalib_.DSPsig_par, sizeof(calibTrk1->DSPsig_par)); |
|
123 |
|
|
124 |
for (int i = 0; i < 6; i++){ |
for (int i = 0; i < 6; i++){ |
125 |
for (int j = 0; j < 3072; j++){ |
for (int j = 0; j < 3072; j++){ |
127 |
calibTrk1->DSPped_par[i][j] = tempPed_par[j][i]; |
calibTrk1->DSPped_par[i][j] = tempPed_par[j][i]; |
128 |
calibTrk1->DSPsig_par[i][j] = tempSig_par[j][i]; |
calibTrk1->DSPsig_par[i][j] = tempSig_par[j][i]; |
129 |
} |
} |
130 |
|
for (int k = 0; k < 3; k++){ |
131 |
|
calibTrk1->crc_cal[i][k] = tempCrc_cal[k][i]; |
132 |
|
} |
133 |
} |
} |
134 |
//----------------------------------------------------------------------------------------- |
//----------------------------------------------------------------------------------------- |
135 |
//delete [] subData; |
|
136 |
} |
} |
137 |
|
|
138 |
|
|