| 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 | } | 
|---|