Changeset 9212 for trunk/MagicSoft/Mars/mraw

Timestamp:

01/06/09 13:07:00 (16 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mraw

Files:

• MRawEvtData.cc (modified) (3 diffs)
• MRawEvtData.h (modified) (3 diffs)
• MRawRunHeader.cc (modified) (6 diffs)
• MRawRunHeader.h (modified) (4 diffs)

trunk/MagicSoft/Mars/mraw/MRawEvtData.cc

@@ -91 +91 @@
 #include "MArrayS.h"
 #include "MArrayB.h"
+#include "MArrayI.h"
 #include "MGeomCam.h"
 
@@ -396 +397 @@
 //  The parameters are the pixelnumber and the FADC_SLICES values of ADCs
 //
-void MRawEvtData::AddPixel(UShort_t nOfPixel, TArrayC *data)
+void MRawEvtData::AddPixel(UShort_t nOfPixel, const TArrayC &data)
 {
     const Int_t n = fRunHeader->GetNumSamples();
-    if (data->GetSize()!=n)
+    if (data.GetSize()!=n)
     {
         *fLog << err << "RawEvtData::AddPixel: Error, number of samples in ";
-        *fLog << "TArrayC " << data->GetSize() << " doesn't match current number " << n << endl;
+        *fLog << "TArrayC " << data.GetSize() << " doesn't match current number " << n << endl;
         return;
     }
@@ -424 +425 @@
     // add the new slices as last entries to array
     //
-    fHiGainFadcSamples->AddAt((Byte_t*)data->GetArray(), fHiGainFadcSamples->GetSize()-n, n);
-}
+    fHiGainFadcSamples->AddAt((Byte_t*)data.GetArray(), fHiGainFadcSamples->GetSize()-n, n);
+}
+
+// --------------------------------------------------------------------------
+//
+// Add the contents of the MArrayI to the fHiGainFadcSamples
+// One Integer is added to one sample in the array.
+//
+void MRawEvtData::Set(const MArrayI &data)
+{
+    const UInt_t n = fRunHeader->GetNumSamples()*fConnectedPixels;
+
+    Byte_t *dest = fHiGainFadcSamples->GetArray();
+
+    UInt_t *src  = reinterpret_cast<UInt_t*>(data.GetArray());
+    UInt_t *end  = reinterpret_cast<UInt_t*>(data.GetArray()) + n;
+
+    switch (fNumBytesPerSample)
+    {
+    case 1:
+        {
+            Byte_t *ptr = reinterpret_cast<Byte_t*>(dest);
+            while (src<end)
+                *ptr++ = *src++;
+        }
+        return;
+
+    case 2:
+        {
+            UShort_t *ptr = reinterpret_cast<UShort_t*>(dest);
+            while (src<end)
+                *ptr++ = *src++;
+        }
+        return;
+
+    case 4:
+        memcpy(dest, data.GetArray(), n*4);
+        return;
+    }
+}
+
 /*
 void MRawEvtData::AddPixel(UShort_t nOfPixel, TArrayC *data, Bool_t lflag)

trunk/MagicSoft/Mars/mraw/MRawEvtData.h

@@ -19 +19 @@
 class MArrayS;
 class MArrayB;
+class MArrayI;
 
 class MRawEvtData : public MParContainer, public MCamEvent
@@ -51 +52 @@
         case 1: return reinterpret_cast<const Byte_t*>(ptr)[n];
         case 2: return reinterpret_cast<const UShort_t*>(ptr)[n];
-        case 3: return reinterpret_cast<const UInt_t*>(ptr)[n];
+        case 4: return reinterpret_cast<const UInt_t*>(ptr)[n];
         }
         return 0;
@@ -74 +75 @@
 
     void ResetPixels(UShort_t npix, UShort_t maxid);
-    void AddPixel(UShort_t nOfPixel, TArrayC *data);
+    void AddPixel(UShort_t nOfPixel, const TArrayC &data);
+    void Set(const MArrayI &data);
 
     UShort_t GetNumHiGainSamples() const;

trunk/MagicSoft/Mars/mraw/MRawRunHeader.cc

@@ -181 +181 @@
     // the default. Do not mix files with and without a value if the
     // files with the value do not match the default!
-    fFormatVersion=0;
+    fFormatVersion=11;
     fSoftVersion=0;
     fTelescopeNumber=1;
@@ -351 +351 @@
     }
     return kTRUE;
+}
+
+void MRawRunHeader::FixRunNumbers()
+{
+    if (fFormatVersion<11 || fTelescopeNumber!=1)
+        return;
+
+    // Map 1:1001349 to 47:1001395
+    if (fRunNumber<48 &&
+        fRunStart.GetMjd()>54674.5 && fRunStart.GetMjd()<56476.5)
+    {
+        fRunNumber += 1001348;
+        *fLog << warn << "Format >V10: Wrong run number increased by 1001348 to " << fRunNumber << "." << endl;
+    }
+
+    // Map 1001916:1001922 to 1002349:1002355
+    if (fRunNumber>1001915 && fRunNumber<1001923 &&
+        fRunStart.GetMjd()>54710.5 && fRunStart.GetMjd()<54711.5)
+    {
+        fRunNumber += 433;
+        *fLog << warn << "Format >V10: Wrong run number increased by 434 to " << fRunNumber << "." << endl;
+    }
+
+    // Map 10000342:1000343 to 10000351:1000351
+    if (fRunNumber>10000342 && fRunNumber<10000352 &&
+        fRunStart.GetMjd()>54254.5 && fRunStart.GetMjd()<54255.5)
+    {
+        fRunNumber -= 9000000;
+        *fLog << warn << "Format >V10: Wrong run number decreased by 9000000 to " << fRunNumber << "." << endl;
+    }
+
+    if (fRunNumber==1000382 &&
+        fRunStart.GetMjd()>54256.5 && fRunStart.GetMjd()<54257.5 &&
+        !strcmp(fSourceName, "GRB080605-2347"))
+    {
+        fFileNumber += 99;
+        *fLog << warn << "Format >V10: Ambiguous file number increased by 99 to " << fFileNumber << "." << endl;
+    }
 }
 
@@ -492 +530 @@
 Bool_t MRawRunHeader::Fixes()
 {
+    FixRunNumbers();
+
     if (fFormatVersion>8 && fRunNumber>326152 && fTelescopeNumber==1)
     {
@@ -497 +537 @@
         fNumEventsRead--;
         *fLog << inf << "Format >V8 and Run>M1:326152: Stored number of events decreased by 1." << endl;
-    }
-
-    if (fFormatVersion>10 && fRunNumber<48 &&
-        fRunStart.GetMjd()>54674.5 && fRunStart.GetMjd()<56476.5)
-    {
-        fRunNumber += 1001348;
-        *fLog << warn << "Format >V10: Wrong run number increased by 1001348 to " << fRunNumber << "." << endl;
     }
 
@@ -760 +793 @@
     return IsConsistent();
 }
-
+/*
+Bool_t MRawRunHeader::WriteEvt(ostream& out) const
+{
+    //
+    // write one RUN HEADER from the input stream
+    //
+    const UInt_t n = fNumCrates*fNumPixInCrate;
+
+    const UShort_t magicnumber     = kMagicNumber;
+    const UInt_t   formatversion   = 11;
+    const UInt_t   headersizerun   = (97+n)*4; //???
+
+    // FIXME: Write fixed number (c0c0)
+    out.write((char*)&magicnumber,   2);     // Total=2
+    out.write((char*)&formatversion, 2);     // Total=4
+    out.write((char*)&formatversion, 4);
+    out.write((char*)&headersizerun, 4);
+
+    TArrayC h(headersizerun-12);
+
+    // ----- convert -----
+    Byte_t  *Char  = reinterpret_cast<Byte_t* >(h.GetArray());
+    UInt_t  *Int   = reinterpret_cast<UInt_t* >(h.GetArray());
+    //const Float_t *Float = reinterpret_cast<Float_t*>(h.GetArray());
+
+    // ----- Start interpretation -----
+
+    Int[0] = 0; // fHeaderSizeEvt;
+    Int[1] = 0; // fHeaderSizeCrate;
+    Int[2] = 0; // fSoftVersion;
+    Int[3] = fFadcType;
+    Int[4] = fCameraVersion;
+    Int[5] = fTelescopeNumber;
+    Int[5] = fRunType;
+    Int[6] = fRunNumber;
+    Int[7] = fFileNumber;
+
+    memcpy(Char+ 36, fProjectName,    100);  // 25
+    memcpy(Char+136, fSourceName,      80);  // 20
+    memcpy(Char+216, fObservationMode, 60);  // 15
+
+    //F32       fSourceRA     = Float[69];
+    //F32       fSourceDEC    = Float[70];
+    //F32       fTelescopeRA  = Float[71];
+    //F32       fTelescopeDEC = Float[72];
+
+    memcpy(Char+232, fSourceEpochChar, 4);
+
+    Int[74] = fSourceEpochDate;
+    Int[75] = fNumCrates;
+    Int[76] = fNumPixInCrate;
+    Int[77] = fNumSamplesHiGain;
+
+    //fNumSamplesRemovedHead = Int[78];
+    //fNumSamplesRemovedTail = Int[79];
+
+    Int[80] = fNumEvents;
+    Int[81] = fNumEvents; //fNumEventsRead;
+    Int[82] = fNumBytesPerSample;
+    Int[83] = fSamplingFrequency;
+    Int[84] = fFadcResolution;
+
+    fRunStart.GetBinary(Int+85);
+    fRunStop.GetBinary(Int+91);
+
+    // ----- 388 bytes so far -----
+
+    //const UInt_t n = fNumCrates*fNumPixInCrate;
+    //if (fHeaderSizeRun<(97+n)*4)
+    //{
+    //    *fLog << err << "ERROR - Event header too small to contain pix assignment." << endl;
+    //    return kFALSE;
+    //}
+
+    // ----- Pixel Assignment -----
+    for (UInt_t i=0; i<n; i++)
+        Int[97+i] = (*fPixAssignment)[i];
+
+    out.write(h.GetArray(), h.GetSize());
+
+    return out;
+}
+*/
 // --------------------------------------------------------------------------
 //
@@ -990 +1105 @@
     return 0;
 }
+
+// --------------------------------------------------------------------------
+//
+// Monte Carlo Interface
+//
+// This is a (prelimiary) way to setup an existing FADC system.
+//
+//  1: Siegen FADCs
+//  2: MUX FADCs
+//
+void MRawRunHeader::InitFadcType(UShort_t type)
+{
+    switch (type)
+    {
+    case 1:
+        fNumSamplesLoGain  =   15;
+        fNumSamplesHiGain  =   15;
+        fNumBytesPerSample =    1;    // number of bytes per sample
+        fSamplingFrequency =  300;    // Sampling Frequency [MHz]
+        fFadcResolution    =    8;    // number of significant bits
+        break;
+    case 2:
+        fNumSamplesLoGain  =    0;
+        fNumSamplesHiGain  =   50;
+        fNumBytesPerSample =    2;    // number of bytes per sample
+        fSamplingFrequency = 2000;    // Sampling Frequency [MHz]
+        fFadcResolution    =   12;    // number of significant bits
+        break;
+    }
+
+    fFadcType = type;
+}
+
+// --------------------------------------------------------------------------
+//
+// Monte Carlo Interface
+//
+//  Init a camera
+//
+void MRawRunHeader::InitCamera(UShort_t type)
+{
+    switch (type)
+    {
+    case 1:
+        fNumCrates     =   1;
+        fNumPixInCrate = 577;
+        break;
+    }
+
+    fCameraVersion = type;
+
+    const Int_t n = fNumCrates*fNumPixInCrate;
+
+    fPixAssignment->Set(n);
+
+    for (int i=0; i<n; i++)
+        (*fPixAssignment)[i] = i;
+}
+
+// --------------------------------------------------------------------------
+//
+// Monte Carlo Interface
+//
+//  Set run-type, telescope number, run-number and file-number
+//
+void MRawRunHeader::SetRunInfo(UShort_t type, UShort_t tel, UInt_t run, UInt_t file)
+{
+    fRunType         = type;
+    fTelescopeNumber = tel;
+    fRunNumber       = run;
+    fFileNumber      = file;
+}
+
+// --------------------------------------------------------------------------
+//
+// Monte Carlo Interface
+//
+//  Set source-name, epoch (default J) and date (default 2000)
+//
+void MRawRunHeader::SetSourceInfo(const char src[80], char epoch, UShort_t date)
+{
+    memcpy(fSourceName, src, 80);
+
+    fSourceEpochChar[0] = epoch;     // epoch char of the source
+    fSourceEpochChar[1] = 0;         // epoch char of the source
+
+    fSourceEpochDate = date;         // epoch date of the source
+}
+
+// --------------------------------------------------------------------------
+//
+// Monte Carlo Interface
+//
+//  Set run-start and -stop time
+//
+void MRawRunHeader::SetRunTime(const MTime &start, const MTime &end)
+{
+    fRunStart = start;
+    fRunStop  = end;
+}
+
+// --------------------------------------------------------------------------
+//
+// Monte Carlo Interface
+//
+//  Set project name and observation mode
+//
+void MRawRunHeader::SetObservation(const char mode[60], const char proj[100])
+{
+    memcpy(fProjectName,     proj, 100);
+    memcpy(fObservationMode, mode,  60);
+}
+
+// --------------------------------------------------------------------------
+//
+// Monte Carlo Interface
+//
+//  Set number of events in file.
+//
+void MRawRunHeader::SetNumEvents(UInt_t num)
+{
+    fNumEvents     = num;
+    fNumEventsRead = num;
+}

trunk/MagicSoft/Mars/mraw/MRawRunHeader.h

@@ -72 +72 @@
 
     Bool_t SwapAssignment(Short_t id0, Short_t id1);
+    void   FixRunNumbers();
     Bool_t FixAssignment();
     Bool_t Fixes();
@@ -83 +84 @@
 
     // This is to be used in the MC chain only!
+/*
     void SetMagicNumber(UShort_t a)       { fMagicNumber=a; }
     void SetFormatVersion(UShort_t a)     { fFormatVersion=a; }
@@ -96 +98 @@
     void SetNumCrates(UShort_t a)         { fNumCrates=a; }
     void SetNumPixInCrate(UShort_t a)     { fNumPixInCrate=a; }
-    void SetRunTime(Float_t start, Float_t stop)
-    { fRunStart.SetMjd(start); fRunStop.SetMjd(stop); }
+*/
+    void InitFadcType(UShort_t type);
+    void InitCamera(UShort_t type);
+    void SetRunInfo(UShort_t type, UShort_t tel, UInt_t run, UInt_t file=0);
+    void SetSourceInfo(const char src[80], char epoch='J', UShort_t date=2000);
+    void SetRunTime(const MTime &start, const MTime &end);
+    void SetRunTimeMjd(Float_t start, Float_t stop) { fRunStart.SetMjd(start); fRunStop.SetMjd(stop); }
+    void SetObservation(const char mode[60], const char proj[100]);
+    void SetNumEvents(UInt_t num);
 
     // This is to get the numbers...
@@ -153 +162 @@
 
     Bool_t ReadEvt(istream& fin);
+    //Bool_t WriteEvt(ostream& fout) const;
 
     ClassDef(MRawRunHeader, 10) // storage container for general info