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

Annotation of /yoda/techmodel/CalibTrk1Reader.cpp

Parent Directory Parent Directory | Revision Log Revision Log

Revision 5.1 - (hide annotations) (download)
Sat Feb 4 12:37:44 2006 UTC (18 years, 10 months ago) by kusanagi
Branch: MAIN
Changes since 5.0: +44 -48 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/CalibTrk1Reader.cpp,v $
3 kusanagi 5.1 * $Id: CalibTrk1Reader.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 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 kusanagi 1.6 #include <log4cxx/logger.h>
15 kusanagi 5.1 #include <fstream>
16     #include "stdio.h"
17     #include "ReaderAlgorithms.h"
18     #include "event/PamelaRun.h"
19    
20     #include "event/CalibTrk1Event.h"
21    
22 kusanagi 1.1 extern "C" {
23 kusanagi 1.3 #include "CRC.h"
24 kusanagi 5.1 extern void trkcalibpkt_(int*, unsigned char[], long int*, int*); //(*)
25 kusanagi 1.1
26     //Struct per il passaggio di dati da e verso la chiamata fortran
27     extern struct {
28 kusanagi 5.1 int good0;
29     int DAQmode[6];
30     int DSPnumber[6];
31     int calibnumber[6];
32     int ncalib_event[6];
33     int ped_l1[6];
34     int ped_l2[6];
35     int ped_l3[6];
36     int sig_l1[6];
37     int sig_l2[6];
38     int sig_l3[6];
39     int nbad_l1[6];
40     int nbad_l2[6];
41     int nbad_l3[6];
42     int cal_flag[6];
43     float DSPped_par[3072][6];
44     float DSPsig_par[3072][6];
45     int DSPbad_par[3072][6];
46     int crc_cal[3][6];
47     int crc_hcal[6];
48 kusanagi 1.1 } trkcalib_;
49     #include <dirent.h>
50     }
51    
52     using namespace pamela;
53     using namespace pamela::techmodel;
54    
55 kusanagi 1.6 static log4cxx::LoggerPtr logger = log4cxx::Logger::getLogger(_T("pamela.techmodel.CalibTrk1Reader"));
56 kusanagi 1.1
57     /**
58     * Constructor.
59     */
60 kusanagi 1.2 CalibTrk1Reader::CalibTrk1Reader(void):
61     TechmodelAlgorithm(PacketType::CalibTrk1, "TechmodelCalibTrk1Reader") {
62 kusanagi 1.6 logger->debug(_T("Constructor"));
63 kusanagi 1.2 calibTrk1 = new CalibTrk1Event();
64 kusanagi 1.1 }
65    
66     /**
67     * Get a string with the version info of the algorithm.
68     */
69 kusanagi 1.2 std::string CalibTrk1Reader::GetVersionInfo(void) const {
70 kusanagi 1.1 return
71 kusanagi 5.1 "$Header: /home/cvsmanager/yoda/techmodel/CalibTrk1Reader.cpp,v 5.0 2005/08/29 09:46:13 Maurizio Nagni Exp $\n";
72 kusanagi 1.1 }
73    
74     /**
75     * Initialize the algorithm with a special run. This will initialize the
76     * event reader routines for all packet types.
77     */
78 kusanagi 1.2 void CalibTrk1Reader::Init(PamelaRun *run) {
79 kusanagi 1.6 logger->debug(_T("Initialize"));
80 kusanagi 1.1 SetInputStream(run);
81 kusanagi 1.2 run->WriteSubPacket(this, &calibTrk1, calibTrk1->Class());
82 kusanagi 1.1 }
83    
84     /**
85     * Unpack the CalibTrk event from an input file.
86     */
87     //void CalibTrkReader::RunEvent(int EventNumber, long int length, char *subData) {
88 kusanagi 2.2 void CalibTrk1Reader::RunEvent(int EventNumber, long int dataLength) throw (Exception){
89 kusanagi 1.6 std::stringstream oss;
90 kusanagi 2.2 char subData[dataLength];
91     int ERROR;
92     InputFile->read(subData, sizeof(unsigned char)*dataLength);
93    
94 kusanagi 5.1 char *data = new char[dataLength];
95     memcpy(data, subData, dataLength);
96     int curpos=1;
97    
98 kusanagi 2.2 //Call to the FORTRAN routin that unpack tracker events
99 kusanagi 5.1 trkcalibpkt_( &ERROR,(unsigned char*)data, &dataLength, &curpos);
100 kusanagi 2.6
101     calibTrk1->unpackError = ERROR;
102 kusanagi 2.2 if (ERROR != 0) {
103     oss.str("");
104     oss << "Fortran77 function trkcalibpkt error code = " << ERROR;
105     logger->warn(oss.str().c_str());
106     }
107 kusanagi 1.1
108 kusanagi 2.2 //Store the unpacked data
109 kusanagi 5.1 calibTrk1->good0 = trkcalib_.good0;
110 kusanagi 2.2 memcpy(calibTrk1->DAQmode, trkcalib_.DAQmode, sizeof(calibTrk1->DAQmode));
111     memcpy(calibTrk1->DSPnumber, trkcalib_.DSPnumber, sizeof(calibTrk1->DSPnumber));
112     memcpy(calibTrk1->calibnumber, trkcalib_.calibnumber, sizeof(calibTrk1->calibnumber));
113 kusanagi 5.1 memcpy(calibTrk1->ncalib_event, trkcalib_.ncalib_event, sizeof(calibTrk1->ncalib_event));
114 kusanagi 2.2 memcpy(calibTrk1->ped_l1, trkcalib_.ped_l1, sizeof(calibTrk1->ped_l1));
115     memcpy(calibTrk1->ped_l2, trkcalib_.ped_l2, sizeof(calibTrk1->ped_l2));
116     memcpy(calibTrk1->ped_l3, trkcalib_.ped_l3, sizeof(calibTrk1->ped_l3));
117     memcpy(calibTrk1->sig_l1, trkcalib_.sig_l1, sizeof(calibTrk1->sig_l1));
118     memcpy(calibTrk1->sig_l2, trkcalib_.sig_l2, sizeof(calibTrk1->sig_l2));
119     memcpy(calibTrk1->sig_l3, trkcalib_.sig_l3, sizeof(calibTrk1->sig_l3));
120     memcpy(calibTrk1->nbad_l1, trkcalib_.nbad_l1, sizeof(calibTrk1->nbad_l1));
121     memcpy(calibTrk1->nbad_l2, trkcalib_.nbad_l2, sizeof(calibTrk1->nbad_l2));
122     memcpy(calibTrk1->nbad_l3, trkcalib_.nbad_l3, sizeof(calibTrk1->nbad_l3));
123     memcpy(calibTrk1->cal_flag, trkcalib_.cal_flag, sizeof(calibTrk1->cal_flag));
124 kusanagi 5.1
125 kusanagi 4.2 int tempBad_par[3072][6];
126     float tempPed_par[3072][6];
127     float tempSig_par[3072][6];
128 kusanagi 5.1 int tempCrc_cal[3][6];
129 kusanagi 4.3 memcpy(tempBad_par,trkcalib_.DSPbad_par, sizeof(tempBad_par));
130     memcpy(tempPed_par,trkcalib_.DSPped_par, sizeof(tempPed_par));
131     memcpy(tempSig_par,trkcalib_.DSPsig_par, sizeof(tempSig_par));
132 kusanagi 5.1 memcpy(tempCrc_cal,trkcalib_.crc_cal, sizeof(tempCrc_cal));
133 kusanagi 2.2
134     for (int i = 0; i < 6; i++){
135     for (int j = 0; j < 3072; j++){
136     calibTrk1->DSPbad_par[i][j] = tempBad_par[j][i];
137     calibTrk1->DSPped_par[i][j] = tempPed_par[j][i];
138     calibTrk1->DSPsig_par[i][j] = tempSig_par[j][i];
139 kusanagi 1.4 }
140 kusanagi 5.1 for (int k = 0; k < 3; k++){
141     calibTrk1->crc_cal[i][k] = tempCrc_cal[k][i];
142     }
143 kusanagi 2.2 }
144 kusanagi 1.4 //-----------------------------------------------------------------------------------------
145 kusanagi 4.3
146 kusanagi 1.1 }
147    
148 kusanagi 2.1
  ViewVC Help
Powered by ViewVC 1.1.23