1 | /* ======================================================================== *\
|
---|
2 | !
|
---|
3 | ! *
|
---|
4 | ! * This file is part of MARS, the MAGIC Analysis and Reconstruction
|
---|
5 | ! * Software. It is distributed to you in the hope that it can be a useful
|
---|
6 | ! * and timesaving tool in analysing Data of imaging Cerenkov telescopes.
|
---|
7 | ! * It is distributed WITHOUT ANY WARRANTY.
|
---|
8 | ! *
|
---|
9 | ! * Permission to use, copy, modify and distribute this software and its
|
---|
10 | ! * documentation for any purpose is hereby granted without fee,
|
---|
11 | ! * provided that the above copyright notice appear in all copies and
|
---|
12 | ! * that both that copyright notice and this permission notice appear
|
---|
13 | ! * in supporting documentation. It is provided "as is" without express
|
---|
14 | ! * or implied warranty.
|
---|
15 | ! *
|
---|
16 | !
|
---|
17 | !
|
---|
18 | ! Author(s): Thomas Bretz 12/2000 <mailto:tbretz@astro.uni-wuerzburg.de>
|
---|
19 | !
|
---|
20 | ! Copyright: MAGIC Software Development, 2000-2004
|
---|
21 | !
|
---|
22 | !
|
---|
23 | \* ======================================================================== */
|
---|
24 |
|
---|
25 | /////////////////////////////////////////////////////////////////////////////
|
---|
26 | //
|
---|
27 | // MRawCrateData
|
---|
28 | //
|
---|
29 | // This container stores the information about one crate. A list of this
|
---|
30 | // informations can be find at MRawCrateArray
|
---|
31 | //
|
---|
32 | /////////////////////////////////////////////////////////////////////////////
|
---|
33 | #include "MRawCrateData.h"
|
---|
34 |
|
---|
35 | #include <iostream>
|
---|
36 | #include <iomanip>
|
---|
37 |
|
---|
38 | #include <fstream>
|
---|
39 |
|
---|
40 | #include "MLog.h"
|
---|
41 | #include "MLogManip.h"
|
---|
42 |
|
---|
43 | ClassImp(MRawCrateData);
|
---|
44 |
|
---|
45 | using namespace std;
|
---|
46 |
|
---|
47 | MRawCrateData::MRawCrateData() : fDAQCrateNumber(0), fFADCEvtNumber(0), fFADCClockTick(0), fABFlags(0)
|
---|
48 | {
|
---|
49 | }
|
---|
50 |
|
---|
51 | // --------------------------------------------------------------------------
|
---|
52 | //
|
---|
53 | // read the information from a binary input stream about the CRATE DATA,
|
---|
54 | // like specified in a TDAS note
|
---|
55 | //
|
---|
56 | void MRawCrateData::ReadEvt(istream& fin, UShort_t ver)
|
---|
57 | {
|
---|
58 | if (ver<7)
|
---|
59 | {
|
---|
60 | fin.read((char*)&fDAQCrateNumber, 2);
|
---|
61 | fin.read((char*)&fFADCEvtNumber, 4);
|
---|
62 | fin.read((char*)&fFADCClockTick, 4);
|
---|
63 | if (ver>1)
|
---|
64 | fin.read((char*)&fABFlags, 1);
|
---|
65 | }
|
---|
66 | else
|
---|
67 | {
|
---|
68 | // U8 CrateNumber; // 0-4
|
---|
69 | // U8 BoardNumber; // 0-1
|
---|
70 | // U8 ChannelNumber; // 0-3
|
---|
71 | fDAQCrateNumber = 0;
|
---|
72 | fin.read((char*)&fDAQCrateNumber, 1); // U8: CrateNumber 0-4
|
---|
73 |
|
---|
74 | Byte_t dummyb;
|
---|
75 | fin.read((char*)&dummyb, 1); // U8 Board number 0-1
|
---|
76 | fin.read((char*)&dummyb, 1); // U8 Channel number 0-3
|
---|
77 |
|
---|
78 | fin.read((char*)&fFADCEvtNumber, 4); // U32 CrateEvtNumber
|
---|
79 |
|
---|
80 | // Clock count. The Clock is synchronized with the 10 MHz external clock,
|
---|
81 | // which is feed to for all FADC channels.
|
---|
82 | // The units are [psec], which is obviously much smaller than the real accuracy.
|
---|
83 | // The ClockTick should be identical for all channels of the same board.
|
---|
84 | // Still keep it for debugging purposes.
|
---|
85 | // Like in the 300MHz system this number is extremely useful to check the
|
---|
86 | // integrity of the data.
|
---|
87 | UInt_t dummyi;
|
---|
88 | fin.read((char*)&dummyi, 4); // U32 FadcClockTickHi; // high significant bits
|
---|
89 | fin.read((char*)&dummyi, 4); // U32 FadcClockTickLo; // low significant bits
|
---|
90 |
|
---|
91 | // Trigger Time Interpolation in [psec] (originally it is a double
|
---|
92 | // in Acqiris software). Again this number should be identical for
|
---|
93 | // all channels in the same board. It is not clear at the moment
|
---|
94 | // if this number will be useful in the end, but I propose to keep
|
---|
95 | // it. The data volume is increase by <0.1%
|
---|
96 | UInt_t dummys;
|
---|
97 | fin.read((char*)&dummys, 2); // U16 TrigTimeInterpol;
|
---|
98 | }
|
---|
99 | }
|
---|
100 |
|
---|
101 | void MRawCrateData::SkipEvt(istream &fin, UShort_t ver)
|
---|
102 | {
|
---|
103 | fin.seekg(ver>1?11:10);
|
---|
104 | }
|
---|
105 |
|
---|
106 | // --------------------------------------------------------------------------
|
---|
107 | //
|
---|
108 | // print all stored information to gLog
|
---|
109 | //
|
---|
110 | void MRawCrateData::Print(Option_t *t) const
|
---|
111 | {
|
---|
112 | *fLog << all;
|
---|
113 | *fLog << "Crate #" << dec << fDAQCrateNumber << ": ";
|
---|
114 | *fLog << "FADCEventNr = " << fFADCEvtNumber << " ";
|
---|
115 | *fLog << "FADCClockTick = " << fFADCClockTick << " (20MHz) ";
|
---|
116 | *fLog << "ABFlags = 0x" << Form("%02x", fABFlags) << endl;
|
---|
117 | }
|
---|