/[PAMELA software]/yoda/techmodel/CalibCalPedReader.cpp
ViewVC logotype

Diff of /yoda/techmodel/CalibCalPedReader.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.2 by kusanagi, Thu Jul 8 12:31:42 2004 UTC revision 1.4 by kusanagi, Thu Aug 19 15:24:46 2004 UTC
# Line 1  Line 1 
 /** @file  
  * $Source: /home/cvsmanager/yoda/techmodel/CalibCalPedReader.cpp,v $  
  * $Id: CalibCalPedReader.cpp,v 1.1 2004/07/06 14:07:27 kusanagi Exp $  
  * $Author: kusanagi $  
  *  
  * Implementation of the LogReader class.  
 * ToBeDone:  
 * Control the CRC for the entire data Packet not just for single records  
  */  
1    
 //#define UINT unsigned int  
2  #define BYTE  unsigned char  #define BYTE  unsigned char
3  #include <string>  #include <string>
4  #include <log4cpp/Category.hh>  #include <log4cpp/Category.hh>
5  #include <fstream>  #include <fstream>
6  #include "stdio.h"  #include "stdio.h"
7  extern "C" {  extern "C" {
8      #include "CRC.h"  
9      //Struct per il passaggio di dati da e verso la chiamata fortran      //Struct per il passaggio di dati da e verso la chiamata fortran
10      extern struct {      extern struct {
11          int IEV2;          int IEV2;
12          int calped[96][11][4];          int calped[4][11][96];
13          int calgood[96][11][4];          int calgood[4][11][96];
14          int calthr[96][11][4];          int calthr[4][11][6];
15          int calrms[96][11][4];          int calrms[4][11][96];
16          int calbase[96][11][4];          int calbase[4][11][6];
17          int calvar[96][11][4];          int calvar[4][11][96];
18          int calpuls[96][11][4];          int calpuls[4][11][96];
19      } calib_;      } calib_;
20            
21      //external declaration of the Fortran function      //external declaration of the Fortran function
22  //    void calpedestal_(short[], int, int*);      void calpedestal_(short[], long int*, int*);
23  }  }
24    
25  #include "ReaderAlgorithms.h"  #include "ReaderAlgorithms.h"
26    
27  #include "event/CalibCalPedEvent.h"  #include "event/CalibCalPedEvent.h"
28    
29  using namespace pamela;  using namespace pamela;
# Line 48  CalibCalPedReader::CalibCalPedReader(voi Line 39  CalibCalPedReader::CalibCalPedReader(voi
39    cat <<  log4cpp::Priority::DEBUG    cat <<  log4cpp::Priority::DEBUG
40        <<  "Constructor "        <<  "Constructor "
41        <<  "\n " << log4cpp::CategoryStream::ENDLINE;        <<  "\n " << log4cpp::CategoryStream::ENDLINE;
42    CalibCalPed = new CalibCalPedEvent();    calibCalPed = new CalibCalPedEvent();
43  }  }
44    
45  /**  /**
# Line 56  CalibCalPedReader::CalibCalPedReader(voi Line 47  CalibCalPedReader::CalibCalPedReader(voi
47   */   */
48  std::string CalibCalPedReader::GetVersionInfo(void) const {  std::string CalibCalPedReader::GetVersionInfo(void) const {
49    return    return
50      "$Header: /home/cvsmanager/yoda/techmodel/CalibCalPedReader.cpp,v 1.1 2004/07/06 14:07:27 kusanagi Exp $\n";      "$Header: /home/cvsmanager/yoda/techmodel/CalibCalPedReader.cpp,v 1.3 2004/07/17 20:03:38 kusanagi Exp $\n";
51  }  }
52    
53  /**  /**
# Line 65  std::string CalibCalPedReader::GetVersio Line 56  std::string CalibCalPedReader::GetVersio
56   */   */
57  void CalibCalPedReader::Init(PamelaRun *run) {  void CalibCalPedReader::Init(PamelaRun *run) {
58    SetInputStream(run);    SetInputStream(run);
59    run->WriteSubPacket(this, &CalibCalPed, CalibCalPed->Class());    run->WriteSubPacket(this, &calibCalPed, calibCalPed->Class());
60  }  }
61    
62  /**  /**
# Line 73  void CalibCalPedReader::Init(PamelaRun * Line 64  void CalibCalPedReader::Init(PamelaRun *
64   */   */
65  void CalibCalPedReader::RunEvent(int EventNumber, long int length) {  void CalibCalPedReader::RunEvent(int EventNumber, long int length) {
66    
67    char     *subData;      char        packetData[length-2];
68    long int  dataLength;      char        CRCevent[2];
69    int ERROR;      UINT16      calculatedCRC    = 0;   //calculated CRC
70          UINT16      readCRC          = 0;   //read CRC
71      //the 2 bytes subtracted belong to the final event CRC bytes      long int    dataLength;
72      dataLength = length - (long int)2;      int         ERROR;
73    
74      subData = new char[dataLength];      dataLength = length - 2;
75      InputFile->read(subData, sizeof(unsigned char)*dataLength);      InputFile->read(packetData, sizeof(packetData));
76        InputFile->read(CRCevent, sizeof(CRCevent));
77      //Skip the last two crc bytes already checked in UnpackPscu  
78      //This part have to be refactored!!!! too bad......      calculatedCRC = CM_Compute_CRC16(0, (BYTE*)packetData, dataLength);
79      InputFile->seekg((long int)2, std::ios::cur);      readCRC = ((UINT16)(CRCevent[0]<<8)&0xFF00) + ((UINT16)(CRCevent[1])&0x00FF);
80                      
81      //Chiamata alla funzione fortran per la lettura dei piedistalli      if (calculatedCRC == readCRC) {
82    //  calpedestal_((short*)subData, dataLength, &ERROR);              calpedestal_((short*)packetData, &dataLength, &ERROR);
83          cat <<  log4cpp::Priority::ERROR    //Store the unpacked data
84              <<  "Fortran77 function calpedestal error code = " << ERROR      calibCalPed->IEV2 = calib_.IEV2;
85              <<  "\n " << log4cpp::CategoryStream::ENDLINE;  //--------have to invert array because of FORTRAN <-> C different management of the indexes
86      free(subData);      int tempCalped[96][11][4];
87    }      int tempCalgood[96][11][4];
88        int tempCalthr[6][11][4];
89        int tempCalrms[96][11][4];
90        int tempCalbase[6][11][4];
91        int tempCalvar[6][11][4];
92        int tempCalpuls[96][11][4];
93    
94        memcpy(tempCalped,  calib_.calped,  sizeof(tempCalped));
95        memcpy(tempCalgood, calib_.calgood, sizeof(tempCalgood));
96        memcpy(tempCalthr,  calib_.calthr,  sizeof(tempCalthr));
97        memcpy(tempCalrms,  calib_.calrms,  sizeof(tempCalrms));
98        memcpy(tempCalbase, calib_.calbase, sizeof(tempCalbase));
99        memcpy(tempCalvar,  calib_.calvar,  sizeof(tempCalvar));
100        memcpy(tempCalpuls, calib_.calpuls, sizeof(tempCalpuls));
101    
102        for (int i = 0; i < 4; i++){
103            for (int j = 0; j <11; j++){
104                for (int z = 0; z < 96; z++){
105                    calibCalPed->calped[i][j][z]    = tempCalped[z][j][i];
106                    calibCalPed->calgood[i][j][z]   = tempCalgood[z][j][i];
107                    calibCalPed->calrms[i][j][z]    = tempCalrms[z][j][i];
108                    calibCalPed->calpuls[i][j][z]    = tempCalpuls[z][j][i];
109                }
110            }
111        }
112    
113        for (int i = 0; i < 4; i++){
114            for (int j = 0; j <11; j++){
115                for (int z = 0; z < 6; z++){
116                    calibCalPed->calthr[i][j][z]    = tempCalthr[z][j][i];
117                    calibCalPed->calbase[i][j][z]   = tempCalbase[z][j][i];
118                    calibCalPed->calvar[i][j][z]    = tempCalvar[z][j][i];
119                }
120            }
121        }
122    //-----------------------------------------------------------------------------------------
123    
124        cat <<  log4cpp::Priority::ERROR
125            <<  "Fortran77 function calpedestal error code = " << ERROR
126            <<  "\n " << log4cpp::CategoryStream::ENDLINE;
127        } else {
128                cat <<  log4cpp::Priority::ERROR
129                <<  "The test of calculated CRC with one wrote on file FAILED!!"
130                <<  "\n " << log4cpp::CategoryStream::ENDLINE;    
131        }
132        free(packetData);
133    }
134    
135    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.23