| 1 | #ifndef MARS_MRawRunHeader
|
|---|
| 2 | #define MARS_MRawRunHeader
|
|---|
| 3 | ///////////////////////////////////////////////////////////////////////
|
|---|
| 4 | // //
|
|---|
| 5 | // MRunHeader //
|
|---|
| 6 | // //
|
|---|
| 7 | ///////////////////////////////////////////////////////////////////////
|
|---|
| 8 |
|
|---|
| 9 | #ifndef MARS_MTime
|
|---|
| 10 | #include "MTime.h"
|
|---|
| 11 | #endif
|
|---|
| 12 |
|
|---|
| 13 | class TBuffer;
|
|---|
| 14 | class MArrayS;
|
|---|
| 15 |
|
|---|
| 16 | class MRawRunHeader : public MParContainer
|
|---|
| 17 | {
|
|---|
| 18 | public:
|
|---|
| 19 | //
|
|---|
| 20 | // enum for the Run Type. Monte Carlo Runs have
|
|---|
| 21 | // to have a value greater than 255 (>0xff)
|
|---|
| 22 | //
|
|---|
| 23 | enum {
|
|---|
| 24 | kRTData = 0x0000,
|
|---|
| 25 | kRTPedestal = 0x0001,
|
|---|
| 26 | kRTCalibration = 0x0002,
|
|---|
| 27 | kRTDominoCal = 0x0003,
|
|---|
| 28 | kRTLinearity = 0x0004,
|
|---|
| 29 | kRTPointRun = 0x0007,
|
|---|
| 30 | kRTMonteCarlo = 0x0100,
|
|---|
| 31 | kRTFact = 0xfac7,
|
|---|
| 32 | kRTNone = 0xffff
|
|---|
| 33 | };
|
|---|
| 34 |
|
|---|
| 35 | //
|
|---|
| 36 | // Magic number to detect the magic file type
|
|---|
| 37 | //
|
|---|
| 38 | static const UShort_t kMagicNumber;
|
|---|
| 39 | static const Byte_t kMaxFormatVersion;
|
|---|
| 40 |
|
|---|
| 41 | private:
|
|---|
| 42 | /* ---- Run Header Informations ---- */
|
|---|
| 43 | UShort_t fMagicNumber; // File type identifier
|
|---|
| 44 |
|
|---|
| 45 | UInt_t fHeaderSizeRun; // Size of run header
|
|---|
| 46 | UInt_t fHeaderSizeEvt; // Size of evt header
|
|---|
| 47 | UInt_t fHeaderSizeCrate; // Size of crate header
|
|---|
| 48 |
|
|---|
| 49 | UShort_t fFormatVersion; // File format version
|
|---|
| 50 | UShort_t fSoftVersion; // DAQ software version
|
|---|
| 51 | UShort_t fFadcType; // FADC type (1=Siegen, 2=MUX)
|
|---|
| 52 | UShort_t fCameraVersion; // Camera Version (1=MAGIC I)
|
|---|
| 53 | UShort_t fTelescopeNumber; // Telescope number (1=Magic I)
|
|---|
| 54 | UShort_t fRunType; // Run Type
|
|---|
| 55 | UInt_t fRunNumber; // Run number
|
|---|
| 56 | UInt_t fFileNumber; // File number
|
|---|
| 57 | Char_t fProjectName[101]; // Project name
|
|---|
| 58 | Char_t fSourceName[81]; // Source name
|
|---|
| 59 | Char_t fObservationMode[61]; // observation mode
|
|---|
| 60 | Char_t fSourceEpochChar[4]; // epoch char of the source
|
|---|
| 61 | UShort_t fSourceEpochDate; // epoch date of the source
|
|---|
| 62 | UShort_t fNumCrates; // number of electronic boards
|
|---|
| 63 | UShort_t fNumPixInCrate; // number of pixels in crate
|
|---|
| 64 | UShort_t fNumSamplesLoGain; // number of logain samples stored
|
|---|
| 65 | UShort_t fNumSamplesHiGain; // number of higain samples stored
|
|---|
| 66 | UShort_t fNumBytesPerSample; // number of bytes per sample
|
|---|
| 67 | Bool_t fIsSigned;
|
|---|
| 68 | UInt_t fNumEvents; // number of events stored
|
|---|
| 69 | UInt_t fNumEventsRead; // number of events read by the electronics
|
|---|
| 70 | UShort_t fSamplingFrequency; // Sampling Frequency [MHz]
|
|---|
| 71 | Byte_t fFadcResolution; // number of significant bits
|
|---|
| 72 | MTime fRunStart; // time of run start
|
|---|
| 73 | MTime fRunStop; // time of run stop
|
|---|
| 74 | MArrayS *fPixAssignment; //-> pixel assignment table
|
|---|
| 75 |
|
|---|
| 76 | Bool_t SwapAssignment(Short_t id0, Short_t id1);
|
|---|
| 77 | void FixRunNumbers();
|
|---|
| 78 | Bool_t FixAssignment();
|
|---|
| 79 | Bool_t Fixes();
|
|---|
| 80 | Bool_t IsConsistent() const;
|
|---|
| 81 |
|
|---|
| 82 | Bool_t ReadEvtOld(std::istream& fin);
|
|---|
| 83 |
|
|---|
| 84 | // MParContainer
|
|---|
| 85 | Int_t ReadEnv(const TEnv &env, TString prefix, Bool_t print);
|
|---|
| 86 |
|
|---|
| 87 | public:
|
|---|
| 88 | MRawRunHeader(const char *name=NULL, const char *title=NULL);
|
|---|
| 89 | MRawRunHeader(const MRawRunHeader &h);
|
|---|
| 90 | ~MRawRunHeader();
|
|---|
| 91 |
|
|---|
| 92 | // This is to be used in the MC chain only!
|
|---|
| 93 | /*
|
|---|
| 94 | void SetMagicNumber(UShort_t a) { fMagicNumber=a; }
|
|---|
| 95 | void SetFormatVersion(UShort_t a) { fFormatVersion=a; }
|
|---|
| 96 | void SetSoftVersion(UShort_t a) { fSoftVersion=a; }
|
|---|
| 97 | void SetRunType(UShort_t a) { fRunType=a; }
|
|---|
| 98 | void SetRunNumber(UInt_t a) { fRunNumber=a; }
|
|---|
| 99 | void SetNumEvents(UInt_t a) { fNumEvents=a; }
|
|---|
| 100 | void SetNumSamples(UShort_t low, UShort_t high)
|
|---|
| 101 | {
|
|---|
| 102 | fNumSamplesLoGain=low;
|
|---|
| 103 | fNumSamplesHiGain=high;
|
|---|
| 104 | }
|
|---|
| 105 | void SetNumCrates(UShort_t a) { fNumCrates=a; }
|
|---|
| 106 | void SetNumPixInCrate(UShort_t a) { fNumPixInCrate=a; }
|
|---|
| 107 | */
|
|---|
| 108 | void InitFadcType(UShort_t type);
|
|---|
| 109 | void InitCamera(UShort_t type, UShort_t pix=0);
|
|---|
| 110 | void InitPixels(UShort_t pix) { InitCamera((UShort_t)-1, pix); }
|
|---|
| 111 | void InitFact(UShort_t num, UShort_t pix, UShort_t samples, UShort_t *map=0);
|
|---|
| 112 | void SetRunType(UShort_t type) { fRunType=type; }
|
|---|
| 113 | void SetRunInfo(UShort_t tel, UInt_t run, UInt_t file=0);
|
|---|
| 114 | void SetFileNumber(UInt_t file) { fFileNumber=file; }
|
|---|
| 115 | void SetSourceInfo(const TString src, char epoch='J', UShort_t date=2000);
|
|---|
| 116 | void SetRunTime(const MTime &start, const MTime &end);
|
|---|
| 117 | void SetRunTimeMjd(Float_t start, Float_t stop) { fRunStart.SetMjd(start); fRunStop.SetMjd(stop); }
|
|---|
| 118 | void SetObservation(const TString mode, const TString proj);
|
|---|
| 119 | void SetNumEvents(UInt_t num);
|
|---|
| 120 | void SetValidMagicNumber() { fMagicNumber=kMagicNumber; }
|
|---|
| 121 | void SetFormat(UShort_t file, UShort_t soft) { fFormatVersion=file; fSoftVersion=soft; }
|
|---|
| 122 |
|
|---|
| 123 | // This is to get the numbers...
|
|---|
| 124 | UShort_t GetMagicNumber() const { return fMagicNumber; }
|
|---|
| 125 | UInt_t GetHeaderSizeEvt() const { return fHeaderSizeEvt; }
|
|---|
| 126 | UInt_t GetHeaderSizeCrate() const { return fHeaderSizeCrate; }
|
|---|
| 127 | UShort_t GetFormatVersion() const { return fFormatVersion; }
|
|---|
| 128 | UShort_t GetSoftVersion() const { return fSoftVersion; }
|
|---|
| 129 | UShort_t GetCameraVersion() const { return fCameraVersion; }
|
|---|
| 130 | UInt_t GetRunNumber() const { return fRunNumber; }
|
|---|
| 131 | UInt_t GetFileNumber() const { return fFileNumber; }
|
|---|
| 132 | UInt_t GetTypeID() const { return (fRunNumber/1000000)%100; }
|
|---|
| 133 | UInt_t GetFileID() const { return fRunNumber>1000000?(fRunNumber%1000000)*1000+(fFileNumber%1000):fRunNumber; }
|
|---|
| 134 | TString GetStringID() const;
|
|---|
| 135 | UShort_t GetTelescopeNumber() const { return fTelescopeNumber; }
|
|---|
| 136 | UShort_t GetRunType() const { return fRunType; }
|
|---|
| 137 | const Char_t *GetRunTypeStr() const;
|
|---|
| 138 | Char_t GetRunTypeChar() const;
|
|---|
| 139 | const Char_t *GetProjectName() const { return fProjectName; }
|
|---|
| 140 | const Char_t *GetSourceName() const { return fSourceName; }
|
|---|
| 141 | const Char_t *GetSourceEpocheChar() const { return fSourceEpochChar; }
|
|---|
| 142 | const Char_t *GetObservationMode() const { return fObservationMode; }
|
|---|
| 143 | UShort_t GetSourceEpocheDate() const { return fSourceEpochDate; }
|
|---|
| 144 | UShort_t GetNumCrates() const { return fNumCrates; }
|
|---|
| 145 | UShort_t GetNumPixInCrate() const { return fNumPixInCrate; }
|
|---|
| 146 | UShort_t GetNumSamplesLoGain() const { return fNumSamplesLoGain; }
|
|---|
| 147 | UShort_t GetNumSamplesHiGain() const { return fNumSamplesHiGain; }
|
|---|
| 148 | UInt_t GetNumSamples() const { return fNumSamplesHiGain+fNumSamplesLoGain; }
|
|---|
| 149 | UShort_t GetNumBytesPerSample() const { return fNumBytesPerSample; }
|
|---|
| 150 | Bool_t IsSigned() const { return fIsSigned; }
|
|---|
| 151 | UInt_t GetNumEvents() const { return fNumEvents; }
|
|---|
| 152 | UInt_t GetNumEventsRead() const { return fNumEventsRead; }
|
|---|
| 153 | UShort_t GetFreqSampling() const { return fSamplingFrequency; }
|
|---|
| 154 | const MTime &GetRunStart() const { return fRunStart; }
|
|---|
| 155 | const MTime &GetRunEnd() const { return fRunStop; }
|
|---|
| 156 | Double_t GetRunLength() const { return !fRunStart || !fRunStop ? 0 : fRunStop-fRunStart; }
|
|---|
| 157 | Short_t GetPixAssignment(UShort_t i) const;
|
|---|
| 158 | UShort_t GetMaxPixId() const;
|
|---|
| 159 | UShort_t GetMinPixId() const;
|
|---|
| 160 | UShort_t GetNumConnectedPixels() const;
|
|---|
| 161 | UShort_t GetNumNormalPixels() const;
|
|---|
| 162 | UShort_t GetNumSpecialPixels() const;
|
|---|
| 163 | UInt_t GetScale() const { switch (fNumBytesPerSample) { case 1: return 0x1; case 2: return 0x100; case 4: return 0x1000000; } return 0; }
|
|---|
| 164 | UInt_t GetMax() const { switch (fNumBytesPerSample) { case 1: return 0xff; case 2: return 0xffff; case 4: return 0xffffffff; } return 0; }
|
|---|
| 165 |
|
|---|
| 166 | const MArrayS &GetPixAssignment() const { return *fPixAssignment; }
|
|---|
| 167 |
|
|---|
| 168 | UInt_t GetNumSamplesPerCrate() const
|
|---|
| 169 | {
|
|---|
| 170 | return fNumPixInCrate*(fNumSamplesLoGain+fNumSamplesHiGain);
|
|---|
| 171 | }
|
|---|
| 172 |
|
|---|
| 173 | UShort_t GetNumPixel() const;
|
|---|
| 174 | Int_t GetNumTotalBytes() const;
|
|---|
| 175 |
|
|---|
| 176 | Bool_t IsValid() const { return fMagicNumber==0xc0c0 || fMagicNumber==0xc0c1; }
|
|---|
| 177 | Bool_t IsMonteCarloRun() const { return fRunType>0x00ff; }
|
|---|
| 178 | Bool_t IsDataRun() const { return (fRunType&0xff)==kRTData; }
|
|---|
| 179 | Bool_t IsPedestalRun() const { return (fRunType&0xff)==kRTPedestal; }
|
|---|
| 180 | Bool_t IsCalibrationRun() const { return (fRunType&0xff)==kRTCalibration; }
|
|---|
| 181 | Bool_t IsPointRun() const { return (fRunType&0xff)==kRTPointRun; }
|
|---|
| 182 |
|
|---|
| 183 | Bool_t IsValidRun() const { return fRunType!=kRTNone; }
|
|---|
| 184 |
|
|---|
| 185 | Bool_t IsConsistent(const MRawRunHeader &h) const;
|
|---|
| 186 |
|
|---|
| 187 | void Print(Option_t *t=NULL) const;
|
|---|
| 188 |
|
|---|
| 189 | Bool_t ReadEvt(std::istream& fin);
|
|---|
| 190 | //Bool_t WriteEvt(ostream& fout) const;
|
|---|
| 191 |
|
|---|
| 192 | ClassDef(MRawRunHeader, 12) // storage container for general info
|
|---|
| 193 | };
|
|---|
| 194 | #endif
|
|---|