--- yoda/techmodel/CalibTrk1Reader.cpp 2004/07/06 12:31:55 1.1 +++ yoda/techmodel/CalibTrk1Reader.cpp 2006/05/30 19:10:02 6.1 @@ -1,6 +1,6 @@ /** @file - * $Source: /home/cvsmanager/yoda/techmodel/CalibTrkReader.cpp,v $ - * $Id: CalibTrkReader.cpp,v 1.1.1.1 2004/07/06 12:20:23 kusanagi Exp $ + * $Source: /home/cvsmanager/yoda/techmodel/CalibTrk1Reader.cpp,v $ + * $Id: CalibTrk1Reader.cpp,v 6.0 2006/02/07 17:11:09 kusanagi Exp $ * $Author: kusanagi $ * * Implementation of the LogReader class. @@ -8,136 +8,131 @@ * Control the CRC for the entire data Packet not just for single records */ -#define UINT unsigned int -#define BYTE unsigned char -#include -#include +#include "ReaderAlgorithms.h" + extern "C" { - //The int* return the error code - //The char* is the path to the temporary file - extern void trkcalibpkt_(int*, char*); + #include "CRC.h" + extern void trkcalibpkt_(int*, unsigned char[], long int*, int*); //(*) //Struct per il passaggio di dati da e verso la chiamata fortran extern struct { - int DAQmode[6]; - int DSPnumber[6]; - int calibnumber[6]; - int ncalib_event[6]; - int ped_l1[6]; - int ped_l2[6]; - int ped_l3[6]; - int sig_l1[6]; - int sig_l2[6]; - int sig_l3[6]; - int nbad_l1[6]; - int nbad_l2[6]; - int nbad_l3[6]; - int cal_flag[6]; - int checksum[6]; - int DSPbad_par[6][3072]; - float DSPped_par[6][3072]; - float DSPsig_par[6][3072]; + int good0; + int DAQmode[6]; + int DSPnumber[6]; + int calibnumber[6]; + int ncalib_event[6]; + int ped_l1[6]; + int ped_l2[6]; + int ped_l3[6]; + int sig_l1[6]; + int sig_l2[6]; + int sig_l3[6]; + int nbad_l1[6]; + int nbad_l2[6]; + int nbad_l3[6]; + int cal_flag[6]; + float DSPped_par[3072][6]; + float DSPsig_par[3072][6]; + int DSPbad_par[3072][6]; + int crc_cal[3][6]; + int crc_hcal[6]; } trkcalib_; - #include } -#include -#include "stdio.h" -#include "ReaderAlgorithms.h" -#include "event/PamelaRun.h" - -#include "event/CalibTrkEvent.h" - -using namespace pamela; using namespace pamela::techmodel; -static log4cpp::Category& cat = log4cpp::Category::getInstance("pamela.techmodel.CalibTrkReader"); +static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.CalibTrk1Reader")); /** * Constructor. */ -CalibTrkReader::CalibTrkReader(void): - TechmodelAlgorithm(PacketType::CalibTrk, "TechmodelCalibTrkReader") { - cat << log4cpp::Priority::DEBUG - << "Constructor " - << "\n " << log4cpp::CategoryStream::ENDLINE; - calibTrk = new CalibTrkEvent(); +CalibTrk1Reader::CalibTrk1Reader(void): + TechmodelAlgorithm(PacketType::CalibTrk1, "TechmodelCalibTrk1Reader") { + logger->debug(_T("Constructor")); + calibTrk1 = new CalibTrk1Event(); } /** * Get a string with the version info of the algorithm. */ -std::string CalibTrkReader::GetVersionInfo(void) const { +std::string CalibTrk1Reader::GetVersionInfo(void) const { return - "$Header: /home/cvsmanager/yoda/techmodel/CalibTrkReader.cpp,v 1.1.1.1 2004/07/06 12:20:23 kusanagi Exp $\n"; + "$Header: /home/cvsmanager/yoda/techmodel/CalibTrk1Reader.cpp,v 6.0 2006/02/07 17:11:09 kusanagi Exp $\n"; } /** * Initialize the algorithm with a special run. This will initialize the * event reader routines for all packet types. */ -void CalibTrkReader::Init(PamelaRun *run) { +void CalibTrk1Reader::Init(PamelaRun *run) { + logger->debug(_T("Initialize")); SetInputStream(run); - run->WriteSubPacket(this, &calibTrk, calibTrk->Class()); + run->WriteSubPacket(this, &calibTrk1, calibTrk1->Class()); } /** * Unpack the CalibTrk event from an input file. */ //void CalibTrkReader::RunEvent(int EventNumber, long int length, char *subData) { -void CalibTrkReader::RunEvent(int EventNumber, long int length) { +void CalibTrk1Reader::RunEvent(int EventNumber, long int dataLength) throw (Exception){ + std::stringstream oss; + char subData[dataLength]; int ERROR; - - //the 2 bytes subtracted belong to the final event CRC bytes - long int dataLength = length;// - (long int)2; - - char *subData = new char[dataLength]; InputFile->read(subData, sizeof(unsigned char)*dataLength); - //Skip the last two crc bytes already checked in UnpackPscu - //This part have to be refactored!!!! too bad...... - InputFile->seekg((long int)2, std::ios::cur); - - //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 + "/"; - pathDir = pathDir + PamelaRun::GetRunName(EventNumber) + "/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()); + char *data = new char[dataLength]; + memcpy(data, subData, dataLength); + int curpos=1; + //Call to the FORTRAN routin that unpack tracker events + trkcalibpkt_( &ERROR,(unsigned char*)data, &dataLength, &curpos); + + calibTrk1->unpackError = ERROR; + if (ERROR != 0) { + oss.str(""); + oss << "Fortran77 function trkcalibpkt error code = " << ERROR; + logger->warn(oss.str().c_str()); + } + //Store the unpacked data - memcpy(calibTrk->DAQmode, trkcalib_.DAQmode, sizeof(calibTrk->DAQmode)); - memcpy(calibTrk->DSPnumber, trkcalib_.DSPnumber, sizeof(calibTrk->DSPnumber)); - memcpy(calibTrk->calibnumber, trkcalib_.calibnumber, sizeof(calibTrk->calibnumber)); - memcpy(calibTrk->DSPnumber, trkcalib_.ncalib_event, sizeof(calibTrk->ncalib_event)); - memcpy(calibTrk->ped_l1, trkcalib_.ped_l1, sizeof(calibTrk->ped_l1)); - memcpy(calibTrk->ped_l2, trkcalib_.ped_l2, sizeof(calibTrk->ped_l2)); - memcpy(calibTrk->ped_l3, trkcalib_.ped_l3, sizeof(calibTrk->ped_l3)); - memcpy(calibTrk->sig_l1, trkcalib_.sig_l1, sizeof(calibTrk->sig_l1)); - memcpy(calibTrk->sig_l2, trkcalib_.sig_l2, sizeof(calibTrk->sig_l2)); - memcpy(calibTrk->sig_l3, trkcalib_.sig_l3, sizeof(calibTrk->sig_l3)); - memcpy(calibTrk->nbad_l1, trkcalib_.nbad_l1, sizeof(calibTrk->nbad_l1)); - memcpy(calibTrk->nbad_l2, trkcalib_.nbad_l2, sizeof(calibTrk->nbad_l2)); - memcpy(calibTrk->nbad_l3, trkcalib_.nbad_l3, sizeof(calibTrk->nbad_l3)); - memcpy(calibTrk->cal_flag, trkcalib_.cal_flag, sizeof(calibTrk->cal_flag)); - memcpy(calibTrk->checksum, trkcalib_.checksum, sizeof(calibTrk->checksum)); - memcpy(calibTrk->DSPbad_par,trkcalib_.DSPbad_par, sizeof(calibTrk->DSPbad_par)); - memcpy(calibTrk->DSPped_par,trkcalib_.DSPped_par, sizeof(calibTrk->DSPped_par)); - memcpy(calibTrk->DSPsig_par,trkcalib_.DSPsig_par, sizeof(calibTrk->DSPsig_par)); - - cat << log4cpp::Priority::ERROR - << "Fortran77 function trkcalibpkt error code = " << ERROR - << "\n " << log4cpp::CategoryStream::ENDLINE; - free(subData); + calibTrk1->good0 = trkcalib_.good0; + 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->ncalib_event, 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)); + + int tempBad_par[3072][6]; + float tempPed_par[3072][6]; + float tempSig_par[3072][6]; + int tempCrc_cal[3][6]; + memcpy(tempBad_par,trkcalib_.DSPbad_par, sizeof(tempBad_par)); + memcpy(tempPed_par,trkcalib_.DSPped_par, sizeof(tempPed_par)); + memcpy(tempSig_par,trkcalib_.DSPsig_par, sizeof(tempSig_par)); + memcpy(tempCrc_cal,trkcalib_.crc_cal, sizeof(tempCrc_cal)); + + for (int i = 0; i < 6; i++){ + for (int j = 0; j < 3072; j++){ + calibTrk1->DSPbad_par[i][j] = tempBad_par[j][i]; + calibTrk1->DSPped_par[i][j] = tempPed_par[j][i]; + calibTrk1->DSPsig_par[i][j] = tempSig_par[j][i]; + } + for (int k = 0; k < 3; k++){ + calibTrk1->crc_cal[i][k] = tempCrc_cal[k][i]; + } + } +//----------------------------------------------------------------------------------------- + } +