| 1 |
/** @file |
| 2 |
* $Source: /home/cvsmanager/yoda/techmodel/CalibTrk2Reader.cpp,v $ |
| 3 |
* $Id: CalibTrk2Reader.cpp,v 2.0 2004/09/21 20:50:54 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.0 2004/09/21 20:50:54 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 length) { |
| 86 |
std::stringstream oss; |
| 87 |
char *subData; |
| 88 |
char eventCRC[2]; |
| 89 |
UINT16 subCRC; //CRC of the data |
| 90 |
UINT16 readCRC; //CRC read from the end of the subpacket |
| 91 |
long int dataLength; |
| 92 |
|
| 93 |
//the 2 bytes subtracted belong to the final event CRC bytes |
| 94 |
dataLength = length - (long int)2; |
| 95 |
|
| 96 |
subData = new char[dataLength]; |
| 97 |
InputFile->read(subData, sizeof(unsigned char)*dataLength); |
| 98 |
subCRC = CM_Compute_CRC16(0, (BYTE*)subData, dataLength); |
| 99 |
|
| 100 |
//took the final CRC to compare it with the previous calculated CRC of the data |
| 101 |
InputFile->read(eventCRC, sizeof(eventCRC)); |
| 102 |
readCRC = (((UINT16)(eventCRC[0]<<8))&0xFF00) + (((UINT16)eventCRC[1])&0x00FF); |
| 103 |
|
| 104 |
if (subCRC == readCRC){ |
| 105 |
int ERROR; |
| 106 |
|
| 107 |
//Scrivo un file temporaneo per passarlo alla routine |
| 108 |
//Speriamo di cambiare la routine per passargli un buffer..... |
| 109 |
DIR *dirp; |
| 110 |
std::string pathDir((char*)getenv("YODA_DATA")); |
| 111 |
pathDir = pathDir + "/todatemp.dat"; |
| 112 |
FILE *pfile; |
| 113 |
pfile = fopen((char*)pathDir.c_str(), "wb"); |
| 114 |
fwrite(subData, 1, dataLength, pfile); |
| 115 |
fclose(pfile); |
| 116 |
|
| 117 |
//Call to the FORTRAN routin that unpack tracker events |
| 118 |
trkcalibpkt_(&ERROR, (char*)pathDir.c_str()); |
| 119 |
|
| 120 |
//delete the temporary file |
| 121 |
remove((char*)pathDir.c_str()); |
| 122 |
|
| 123 |
//Store the unpacked data |
| 124 |
memcpy(calibTrk2->DAQmode, trkcalib_.DAQmode, sizeof(calibTrk2->DAQmode)); |
| 125 |
memcpy(calibTrk2->DSPnumber, trkcalib_.DSPnumber, sizeof(calibTrk2->DSPnumber)); |
| 126 |
memcpy(calibTrk2->calibnumber, trkcalib_.calibnumber, sizeof(calibTrk2->calibnumber)); |
| 127 |
memcpy(calibTrk2->DSPnumber, trkcalib_.ncalib_event, sizeof(calibTrk2->ncalib_event)); |
| 128 |
memcpy(calibTrk2->ped_l1, trkcalib_.ped_l1, sizeof(calibTrk2->ped_l1)); |
| 129 |
memcpy(calibTrk2->ped_l2, trkcalib_.ped_l2, sizeof(calibTrk2->ped_l2)); |
| 130 |
memcpy(calibTrk2->ped_l3, trkcalib_.ped_l3, sizeof(calibTrk2->ped_l3)); |
| 131 |
memcpy(calibTrk2->sig_l1, trkcalib_.sig_l1, sizeof(calibTrk2->sig_l1)); |
| 132 |
memcpy(calibTrk2->sig_l2, trkcalib_.sig_l2, sizeof(calibTrk2->sig_l2)); |
| 133 |
memcpy(calibTrk2->sig_l3, trkcalib_.sig_l3, sizeof(calibTrk2->sig_l3)); |
| 134 |
memcpy(calibTrk2->nbad_l1, trkcalib_.nbad_l1, sizeof(calibTrk2->nbad_l1)); |
| 135 |
memcpy(calibTrk2->nbad_l2, trkcalib_.nbad_l2, sizeof(calibTrk2->nbad_l2)); |
| 136 |
memcpy(calibTrk2->nbad_l3, trkcalib_.nbad_l3, sizeof(calibTrk2->nbad_l3)); |
| 137 |
memcpy(calibTrk2->cal_flag, trkcalib_.cal_flag, sizeof(calibTrk2->cal_flag)); |
| 138 |
// memcpy(calibTrk2->DSPbad_par, trkcalib_.DSPbad_par, sizeof(calibTrk2->DSPbad_par)); |
| 139 |
// memcpy(calibTrk2->DSPped_par, trkcalib_.DSPped_par, sizeof(calibTrk2->DSPped_par)); |
| 140 |
// memcpy(calibTrk2->DSPsig_par, trkcalib_.DSPsig_par, sizeof(calibTrk2->DSPsig_par)); |
| 141 |
//--------have to invert array because of FORTRAN <-> C different management of the indexes |
| 142 |
int tempBad_par[3072][6]; |
| 143 |
int tempPed_par[3072][6]; |
| 144 |
int tempSig_par[3072][6]; |
| 145 |
memcpy(tempBad_par,trkcalib_.DSPbad_par, sizeof(tempBad_par)); |
| 146 |
memcpy(tempPed_par,trkcalib_.DSPped_par, sizeof(tempPed_par)); |
| 147 |
memcpy(tempSig_par,trkcalib_.DSPsig_par, sizeof(tempSig_par)); |
| 148 |
|
| 149 |
for (int i = 0; i < 6; i++){ |
| 150 |
for (int j = 0; j < 3072; j++){ |
| 151 |
calibTrk2->DSPbad_par[i][j] = tempBad_par[j][i]; |
| 152 |
calibTrk2->DSPped_par[i][j] = tempPed_par[j][i]; |
| 153 |
calibTrk2->DSPsig_par[i][j] = tempSig_par[j][i]; |
| 154 |
} |
| 155 |
} |
| 156 |
//----------------------------------------------------------------------------------------- |
| 157 |
oss.flush(); |
| 158 |
oss << "Fortran77 function trkcalibpkt error code = " << ERROR; |
| 159 |
logger->warn(oss.str().c_str()); |
| 160 |
} else { |
| 161 |
logger->debug(_T("Wrong CRC for CalibTrk2 Packet ")); |
| 162 |
} |
| 163 |
delete [] subData; |
| 164 |
} |
| 165 |
|
Parent Directory
| 
Revision Log