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

Contents of /yoda/techmodel/McmdReader.cpp

Parent Directory Parent Directory | Revision Log Revision Log


Revision 6.0 - (show annotations) (download)
Tue Feb 7 17:11:09 2006 UTC (18 years, 10 months ago) by kusanagi
Branch: MAIN
CVS Tags: YODA6_2/01, YODA6_2/00, YODA6_1/00, YODA6_0/00
Changes since 5.1: +3 -3 lines
Several new features in this revision:
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 /** @file
2 * $Source: /home/cvsmanager/yoda/techmodel/McmdReader.cpp,v $
3 * $Id: McmdReader.cpp,v 5.1 2006/02/04 12:37:45 kusanagi Exp $
4 * $Author: kusanagi $
5 *
6 * Implementation of the McmdReader class.
7 */
8
9 #define BYTE unsigned char
10 #include <string>
11 #include <log4cxx/logger.h>
12 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 static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.McmdReader"));
27 static std::stringstream oss;
28 /**
29 * Constructor.
30 */
31 McmdReader::McmdReader(void):
32 TechmodelAlgorithm(PacketType::Mcmd, "Mcmd") {
33 logger->debug(_T("Constructor"));
34 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 "$Header: /home/cvsmanager/yoda/techmodel/McmdReader.cpp,v 5.1 2006/02/04 12:37:45 kusanagi Exp $\n";
43 }
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 logger->debug(_T("Initialize"));
51 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 void McmdReader::RunEvent(int EventNumber, long int length) throw (WrongCRCException){
79
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
168 rec->McmdData = new TArrayC(dataLength, subData);
169 delete [] subData;
170 } else {
171 oss.str("");
172 oss << "Wrong CRC on Subpacket internal to TMTC Packet starting at position"
173 << start;
174 logger->warn(oss.str().c_str());
175 }
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 if(partialCRC != readCRC) throw WrongCRCException(" Wrong Global CRC for MCMD Packet ");
182 /*if(!(partialCRC == readCRC)) {
183 throw WrongCRCException();
184 }*/
185 }
186
187

  ViewVC Help
Powered by ViewVC 1.1.23