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

Annotation of /yoda/techmodel/McmdReader.cpp

Parent Directory Parent Directory | Revision Log Revision Log

Revision 5.1 - (hide annotations) (download)
Sat Feb 4 12:37:45 2006 UTC (18 years, 9 months ago) by kusanagi
Branch: MAIN
Changes since 5.0: +3 -3 lines 
Several new features in this release:
a) all the packets are conform to the Mass Memory Format specifications (http://people.roma2.infn.it/~cpu/Mass_Memory_Format.html)
b) unpacking either using the old files structure OR the new one file unpacking.
c) parametrized root files compression factor
d) deleting of the following packet: TofTest, TrkTest, TrkEvent.
e) the Tracker routines now work without the use of temp files.

The point a) allow Yoda to unpack in the root file all the packets generated by the CPU. According to the MassMemoryFormat; that is three possible data are available:

1) almost explicit structure of the packet (like for Log, Tracker, Mcmd, etc....);
2) dummy data collection structure (InitHeader, InitTrailer, CalibHeader, CalibTrailer);
3) just the data of the packet (almost all Alarm and Init procedures). The class regarding this packets have only one parameters, a TArrayC class, which contain the data-block included in the packet (tat is the data below the packet Header).

The point b) has been implemented as a consequence of an agreement about a more compact structure of the unpacked data. Up to now the structure of each unpacked data consisted of a folder, named after the packet type, and three files: xxx.Header.root, xxx.NamePacket.root, xxx.Registry.root.
Starting from this release YODA, by default will unpack the data in a unique root file. The structure of this file will consist of:
- several TTree(s) named after the packet type;
- into each TTree are foreseen three TBranche(s):
    - 'Header'  (the old xxx.Header.root file)
    - 'NameOfThePacket' (the old xxx.Event.root file or the xxx.Event.DETECTOR.root)
    - 'Registry' (the old xxx.Registry.root file)

Anyway is still possible, but deprecated, to unpack using the old structure, passing to the "yoda" command the optional parameter "-multifile"

The point c) has been implemented because is well know that writing time in a TTree is as much fast as much lower is the compression factor for the root file; anyway for a PAMELA dat file, a compression equal to 0 will generate a root file which will be more than two times the original size. To modify the compression parameter just add the optional parameter "-c [0-9]" to the yoda command line.
1 kusanagi 1.1 /** @file
2 kusanagi 1.2 * $Source: /home/cvsmanager/yoda/techmodel/McmdReader.cpp,v $
3 kusanagi 5.1 * $Id: McmdReader.cpp,v 5.0 2005/08/29 09:46:13 Maurizio Nagni Exp $
4     * $Author: Maurizio Nagni $
5 kusanagi 1.1 *
6     * Implementation of the McmdReader class.
7     */
8    
9     #define BYTE unsigned char
10     #include <string>
11 kusanagi 1.5 #include <log4cxx/logger.h>
12 kusanagi 1.1 extern "C" {
13     #include <sys/time.h>
14     #include "CRC.h"
15     }
16    
17     #include <fstream>
18     #include "stdio.h"
19     #include "ReaderAlgorithms.h"
20    
21     #include "event/mcmd/McmdRecord.h"
22    
23     using namespace pamela;
24     using namespace pamela::techmodel;
25    
26 kusanagi 1.5 static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.McmdReader"));
27     static std::stringstream oss;
28 kusanagi 1.1 /**
29     * Constructor.
30     */
31     McmdReader::McmdReader(void):
32     TechmodelAlgorithm(PacketType::Mcmd, "Mcmd") {
33 kusanagi 1.5 logger->debug(_T("Constructor"));
34 kusanagi 1.1 Mcmd = new McmdEvent();
35     }
36    
37     /**
38     * Get a string with the version info of the algorithm.
39     */
40     std::string McmdReader::GetVersionInfo(void) const {
41     return
42 kusanagi 5.1 "$Header: /home/cvsmanager/yoda/techmodel/McmdReader.cpp,v 5.0 2005/08/29 09:46:13 Maurizio Nagni Exp $\n";
43 kusanagi 1.1 }
44    
45     /**
46     * Initialize the algorithm with a special run. This will initialize the
47     * event reader routines for all packet types.
48     */
49     void McmdReader::Init(PamelaRun *run) {
50 kusanagi 1.5 logger->debug(_T("Initialize"));
51 kusanagi 1.1 SetInputStream(run);
52     run->WriteSubPacket(this, &Mcmd, Mcmd->Class());
53     }
54    
55     /**
56     * Unpack the Mcmd event from an input file.
57     * Each subpacket is prafaceded by 4 OBT bytes.
58     * The elementar structure is a kind of
59     * --------CPU - OBT---------------------------------
60     * OBT - 4 Bytes
61     * --------Start Sub-Packet---------------------------------
62     * SeqID - 2 Bytes
63     * T - 1 bit | spare - 7 bits | Mcmd-ID - 1 byte
64     * Spare - 4 bits | Mcmd Lenght 12 bits
65     * Time Tag - 4 Bytes
66     * DATA - McmdLength * 2 Bytes
67     * don't care - 1 Byte
68     * END-ID - 1 Byte is the sum of subPacket Bytes module 256
69     * --------End Sub-Packet---------------------------------
70     * --------CPU CRC-16 on OBT + subPacket---------------------------------
71     * subCRC - 2 Bytes
72     * this structure is repeated one or more times.
73     * The last 2 Bytes are CRC16 computed by CPU on all previous (OBT + subpacket data + subCRC) repetitions
74    
75     * @Event Number
76     * @length is the size in bytes of the event (or packet)
77     */
78 kusanagi 2.1 void McmdReader::RunEvent(int EventNumber, long int length) throw (WrongCRCException){
79 kusanagi 1.1
80     int i = 0;
81     char OBT[4];
82     char subHeader[10];
83     char *subData;
84     char subTrailer[2];
85     char subPckCRC[2];
86     char eventCRC[2];
87     long int start;
88     unsigned short headerCRC;
89     UINT16 subCRC; //CRC of the subpacket (updated as subPckt is read)
90     UINT16 readCRC; //CRC read from the subpacket
91     UINT16 partialCRC; //partial CRC updated as mcmd packet is read (to compare with CRC read in the total Mcmd header)
92    
93     long int lastPosition, dataLength;
94     McmdRecord *rec;
95    
96     //the 2 bytes subtracted belong to the final event CRC bytes
97     lastPosition = (long int)InputFile->tellg() + length - 2;
98     Mcmd->Records->Clear();
99     TClonesArray &recs = *(Mcmd->Records);
100     partialCRC = 0;
101     while(InputFile->tellg() < lastPosition) {
102    
103     //Just to test crc
104     /* char tempbuff[length];
105     InputFile->read(tempbuff, sizeof(tempbuff));
106     subCRC = 0;
107     for (int jj = 0; jj < length ; jj++){
108     subCRC = CM_Compute_CRC16(subCRC, (BYTE*)(tempbuff+jj), 1);
109     }*/
110    
111     //Read the OBT preceeding the subpacket then calculate a partial CRC for it
112     //and update the partialCRC
113     start = InputFile->tellg();
114     InputFile->read(OBT, sizeof(OBT));
115     subCRC = CM_Compute_CRC16(0, (BYTE*)&OBT, 4);
116     partialCRC = CM_Compute_CRC16(partialCRC, (BYTE*)&OBT, 4);
117    
118     //Read the header for the subPacket and read mcmdLength
119     //12 is the total lenght of subHeader + subTrailer
120     InputFile->read(subHeader, sizeof(subHeader));
121     subCRC = CM_Compute_CRC16(subCRC, (BYTE*)&subHeader, 10);
122     partialCRC = CM_Compute_CRC16(partialCRC, (BYTE*)&subHeader, 10);
123     dataLength = (((0x0fff)&((UINT16)subHeader[4]))<<8|((UINT16)subHeader[5])*2) - 12;
124    
125     if (dataLength < 0) break; //it should throw an exception ***TBD***
126     //read subpacket data according to data length then calculate partial CRC for data
127     //and update the partialCRC
128    
129     // If dataLength > 0 reads the data and compute again a partial CRC
130     // else set subData to 0
131     //if(dataLength > 0){
132     subData = new char[dataLength];
133     InputFile->read(subData, sizeof(unsigned char)*dataLength);
134     for (int jj = 0; jj < dataLength ; jj++){
135     subCRC = CM_Compute_CRC16(subCRC, (BYTE*)(subData+jj), 1);
136     partialCRC = CM_Compute_CRC16(partialCRC, (BYTE*)(subData+jj), 1);
137     }
138     //subCRC = CM_Compute_CRC16(subCRC, (BYTE*)&subData, dataLength);
139    
140     //} else {
141     // subData = "0";
142     //}
143    
144     //Read the CRC inside of MCMD
145     //To check sum of all sub packet data, module 256 ****TBD****
146     InputFile->read(subTrailer, sizeof(subTrailer));
147     subCRC = CM_Compute_CRC16(subCRC, (BYTE*)&subTrailer, 2);
148     partialCRC = CM_Compute_CRC16(partialCRC, (BYTE*)&subTrailer, 2);
149    
150     //Read the CRC of subPacket
151     InputFile->read(subPckCRC, sizeof(subPckCRC));
152     partialCRC = CM_Compute_CRC16(partialCRC, (BYTE*)&subPckCRC, 2);
153     readCRC = (((BYTE)subPckCRC[0])<<8) + ((BYTE)subPckCRC[1]);
154    
155     //finally check if the RecordCRC is correct
156     //and finally update the partialCRC
157     //TO DO - if one CRC is wrong also the total one will be corrupted
158     if (subCRC == readCRC){
159     rec = new(recs[i++]) McmdRecord(); //aggiungo un nuovo McmdRecord all'evento
160     rec->MCMD_RECORD_OBT = (((UINT32)OBT[0]<<24)&0xFF000000) + (((UINT32)OBT[1]<<16)&0x00FF0000) + (((UINT32)OBT[2]<<8)&0x0000FF00) + (((UINT32)OBT[3])&0x000000FF);
161     rec->SeqID = (((BYTE)subHeader[0]<<8)&0xFF00) + (((BYTE)subHeader[1])&0x00FF);
162     rec->Tbit = ((BYTE)((subHeader[2]&0x80))>>7);
163     rec->ID1 = (BYTE)subHeader[3];
164     rec->McmdLength = (0x0fff)&(((UINT16)(subHeader[4]<<8)) + ((UINT16)subHeader[5]));
165     rec->TimeTag = (((UINT32)OBT[6]<<24)&0xFF000000) + (((UINT32)OBT[7]<<16)&0x00FF0000) + (((UINT32)OBT[8]<<8)&0x0000FF00) + (((UINT32)OBT[9])&0x000000FF);
166     rec->endID = (BYTE)subTrailer[1];
167 kusanagi 1.4
168     rec->McmdData = new TArrayC(dataLength, subData);
169 kusanagi 1.3 delete [] subData;
170 kusanagi 1.1 } else {
171 kusanagi 2.1 oss.str("");
172 kusanagi 1.5 oss << "Wrong CRC on Subpacket internal to TMTC Packet starting at position"
173     << start;
174     logger->warn(oss.str().c_str());
175 kusanagi 1.1 }
176     }
177     //in the end compare the calculated partial CRC with the MCMD packet CRC
178     InputFile->read(eventCRC, sizeof(eventCRC));
179     readCRC = (((UINT16)(eventCRC[0]<<8))&0xFF00) + (((UINT16)eventCRC[1])&0x00FF);
180    
181 kusanagi 2.1 if(partialCRC != readCRC) throw WrongCRCException(" Wrong Global CRC for MCMD Packet ");
182     /*if(!(partialCRC == readCRC)) {
183     throw WrongCRCException();
184     }*/
185 kusanagi 1.1 }
186    
187    
  ViewVC Help
Powered by ViewVC 1.1.23