1 |
/** @file |
2 |
* $Source: /home/cvsmanager/yoda/techmodel/CalibTrk2Reader.cpp,v $ |
3 |
* $Id: CalibTrk2Reader.cpp,v 2.4 2004/12/22 11:38:35 kusanagi Exp $ |
4 |
* $Author: kusanagi $ |
5 |
* |
6 |
* Implementation of the LogReader class. |
7 |
* ToBeDone: |
8 |
* Control the CRC for the entire data Packet not just for single records |
9 |
*/ |
10 |
|
11 |
#define UINT unsigned int |
12 |
#define BYTE unsigned char |
13 |
#include <string> |
14 |
#include <log4cxx/logger.h> |
15 |
extern "C" { |
16 |
#include "CRC.h" |
17 |
//Passo il path verso la il file temporaneo |
18 |
extern void trkcalibpkt_(int*, char*); |
19 |
|
20 |
//Struct per il passaggio di dati da e verso la chiamata fortran |
21 |
extern struct { |
22 |
int DAQmode[6]; |
23 |
int DSPnumber[6]; |
24 |
int calibnumber[6]; |
25 |
int ncalib_event[6]; |
26 |
int ped_l1[6]; |
27 |
int ped_l2[6]; |
28 |
int ped_l3[6]; |
29 |
int sig_l1[6]; |
30 |
int sig_l2[6]; |
31 |
int sig_l3[6]; |
32 |
int nbad_l1[6]; |
33 |
int nbad_l2[6]; |
34 |
int nbad_l3[6]; |
35 |
int cal_flag[6]; |
36 |
int checksum[6]; |
37 |
int DSPbad_par[6][3072]; |
38 |
float DSPped_par[6][3072]; |
39 |
float DSPsig_par[6][3072]; |
40 |
} trkcalib_; |
41 |
#include <dirent.h> |
42 |
} |
43 |
|
44 |
#include <fstream> |
45 |
#include "stdio.h" |
46 |
#include "ReaderAlgorithms.h" |
47 |
|
48 |
#include "event/CalibTrk2Event.h" |
49 |
|
50 |
using namespace pamela; |
51 |
using namespace pamela::techmodel; |
52 |
|
53 |
static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.CalibTrk2Reader")); |
54 |
|
55 |
/** |
56 |
* Constructor. |
57 |
*/ |
58 |
CalibTrk2Reader::CalibTrk2Reader(void): |
59 |
TechmodelAlgorithm(PacketType::CalibTrk2, "TechmodelCalibTrk2Reader") { |
60 |
logger->debug(_T("Constructor")); |
61 |
calibTrk2 = new CalibTrk2Event(); |
62 |
} |
63 |
|
64 |
/** |
65 |
* Get a string with the version info of the algorithm. |
66 |
*/ |
67 |
std::string CalibTrk2Reader::GetVersionInfo(void) const { |
68 |
return |
69 |
"$Header: /home/cvsmanager/yoda/techmodel/CalibTrk2Reader.cpp,v 2.4 2004/12/22 11:38:35 kusanagi Exp $\n"; |
70 |
} |
71 |
|
72 |
/** |
73 |
* Initialize the algorithm with a special run. This will initialize the |
74 |
* event reader routines for all packet types. |
75 |
*/ |
76 |
void CalibTrk2Reader::Init(PamelaRun *run) { |
77 |
logger->debug(_T("Initialize")); |
78 |
SetInputStream(run); |
79 |
run->WriteSubPacket(this, &calibTrk2, calibTrk2->Class()); |
80 |
} |
81 |
|
82 |
/** |
83 |
* Unpack the CalibTrk2 event from an input file. |
84 |
*/ |
85 |
void CalibTrk2Reader::RunEvent(int EventNumber, long int dataLength) throw (Exception){ |
86 |
std::stringstream oss; |
87 |
char subData[dataLength]; |
88 |
int ERROR; |
89 |
|
90 |
InputFile->read(subData, sizeof(unsigned char)*dataLength); |
91 |
|
92 |
//Scrivo un file temporaneo per passarlo alla routine |
93 |
//Speriamo di cambiare la routine per passargli un buffer..... |
94 |
oss.str(""); |
95 |
oss << getenv("YODA_DATA") << "/" << time(NULL) << "trc2.dat"; |
96 |
FILE *pfile; |
97 |
pfile = fopen(oss.str().c_str(), "wb"); |
98 |
fwrite(subData, 1, dataLength, pfile); |
99 |
fclose(pfile); |
100 |
|
101 |
//Call to the FORTRAN routin that unpack tracker events |
102 |
trkcalibpkt_(&ERROR, (char*)oss.str().c_str()); |
103 |
remove(oss.str().c_str()); |
104 |
if (ERROR != 0) { |
105 |
oss.str(""); |
106 |
oss << "Fortran77 function trkcalibpkt error code = " << ERROR; |
107 |
logger->warn(oss.str().c_str()); |
108 |
} |
109 |
//delete the temporary file |
110 |
remove((char*)oss.str().c_str()); |
111 |
|
112 |
//Store the unpacked data |
113 |
memcpy(calibTrk2->DAQmode, trkcalib_.DAQmode, sizeof(calibTrk2->DAQmode)); |
114 |
memcpy(calibTrk2->DSPnumber, trkcalib_.DSPnumber, sizeof(calibTrk2->DSPnumber)); |
115 |
memcpy(calibTrk2->calibnumber, trkcalib_.calibnumber, sizeof(calibTrk2->calibnumber)); |
116 |
memcpy(calibTrk2->DSPnumber, trkcalib_.ncalib_event, sizeof(calibTrk2->ncalib_event)); |
117 |
memcpy(calibTrk2->ped_l1, trkcalib_.ped_l1, sizeof(calibTrk2->ped_l1)); |
118 |
memcpy(calibTrk2->ped_l2, trkcalib_.ped_l2, sizeof(calibTrk2->ped_l2)); |
119 |
memcpy(calibTrk2->ped_l3, trkcalib_.ped_l3, sizeof(calibTrk2->ped_l3)); |
120 |
memcpy(calibTrk2->sig_l1, trkcalib_.sig_l1, sizeof(calibTrk2->sig_l1)); |
121 |
memcpy(calibTrk2->sig_l2, trkcalib_.sig_l2, sizeof(calibTrk2->sig_l2)); |
122 |
memcpy(calibTrk2->sig_l3, trkcalib_.sig_l3, sizeof(calibTrk2->sig_l3)); |
123 |
memcpy(calibTrk2->nbad_l1, trkcalib_.nbad_l1, sizeof(calibTrk2->nbad_l1)); |
124 |
memcpy(calibTrk2->nbad_l2, trkcalib_.nbad_l2, sizeof(calibTrk2->nbad_l2)); |
125 |
memcpy(calibTrk2->nbad_l3, trkcalib_.nbad_l3, sizeof(calibTrk2->nbad_l3)); |
126 |
memcpy(calibTrk2->cal_flag, trkcalib_.cal_flag, sizeof(calibTrk2->cal_flag)); |
127 |
// memcpy(calibTrk2->DSPbad_par, trkcalib_.DSPbad_par, sizeof(calibTrk2->DSPbad_par)); |
128 |
// memcpy(calibTrk2->DSPped_par, trkcalib_.DSPped_par, sizeof(calibTrk2->DSPped_par)); |
129 |
// memcpy(calibTrk2->DSPsig_par, trkcalib_.DSPsig_par, sizeof(calibTrk2->DSPsig_par)); |
130 |
//--------have to invert array because of FORTRAN <-> C different management of the indexes |
131 |
int tempBad_par[3072][6]; |
132 |
int tempPed_par[3072][6]; |
133 |
int tempSig_par[3072][6]; |
134 |
memcpy(tempBad_par,trkcalib_.DSPbad_par, sizeof(tempBad_par)); |
135 |
memcpy(tempPed_par,trkcalib_.DSPped_par, sizeof(tempPed_par)); |
136 |
memcpy(tempSig_par,trkcalib_.DSPsig_par, sizeof(tempSig_par)); |
137 |
|
138 |
for (int i = 0; i < 6; i++){ |
139 |
for (int j = 0; j < 3072; j++){ |
140 |
calibTrk2->DSPbad_par[i][j] = tempBad_par[j][i]; |
141 |
calibTrk2->DSPped_par[i][j] = tempPed_par[j][i]; |
142 |
calibTrk2->DSPsig_par[i][j] = tempSig_par[j][i]; |
143 |
} |
144 |
} |
145 |
//----------------------------------------------------------------------------------------- |
146 |
oss.str(""); |
147 |
oss << "Fortran77 function trkcalibpkt error code = " << ERROR; |
148 |
logger->warn(oss.str().c_str()); |
149 |
|
150 |
//delete [] subData; |
151 |
} |
152 |
|