Changeset 654 for trunk/MagicSoft/Mars/mraw

Timestamp:

03/01/01 15:56:58 (24 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mraw

Files:

• MRawCrateData.cc (modified) (1 diff)
• MRawCrateData.h (modified) (1 diff)
• MRawEvtData.cc (modified) (7 diffs)
• MRawEvtData.h (modified) (3 diffs)
• MRawEvtHeader.cc (modified) (3 diffs)
• MRawEvtHeader.h (modified) (1 diff)
• MRawFileRead.cc (modified) (5 diffs)
• MRawFileRead.h (modified) (1 diff)
• MRawRunHeader.cc (modified) (2 diffs)
• MRawRunHeader.h (modified) (2 diffs)

trunk/MagicSoft/Mars/mraw/MRawCrateData.cc

@@ -12 +12 @@
 }
 
-void MRawCrateData::ReadEvt(ifstream& fin)
+void MRawCrateData::ReadEvt(istream& fin)
 {
     //

trunk/MagicSoft/Mars/mraw/MRawCrateData.h

@@ -42 +42 @@
     void Print(Option_t *t=NULL);
 
-    void ReadEvt(ifstream& fin);
+    void ReadEvt(istream& fin);
 
     ClassDef(MRawCrateData, 1) //Container to store the Raw CRATE DATA

trunk/MagicSoft/Mars/mraw/MRawEvtData.cc

@@ -216 +216 @@
 }
 
-void MRawEvtData::DeletePixels()
+void MRawEvtData::DeletePixels(Bool_t flag)
 {
     //
@@ -222 +222 @@
     //
     DeleteArrays();
-    InitArrays();
+    InitArrays(flag);
 }
 
@@ -233 +233 @@
 }
 
-void MRawEvtData::InitArrays()
-{
-    fHiGainPixId       = new MArrayS(0); //UShort_t[0];
-    fLoGainPixId       = new MArrayS(0); //new UShort_t[0];
-    fHiGainFadcSamples = new MArrayB(0); //new Byte_t[0];
-    fLoGainFadcSamples = new MArrayB(0); //new Byte_t[0];
+void MRawEvtData::InitArrays(Bool_t flag)
+{
+    // const int npix = !flag ? 0 : fRunHeader->GetNumCrates()*fRunHeader->GetNumPixInCrate();
+
+    fHiGainPixId       = new MArrayS(0);//npix);
+    fLoGainPixId       = new MArrayS(0);//npix); 
+    fHiGainFadcSamples = new MArrayB(0);//npix*fRunHeader->GetNumSamplesHiGain());
+    fLoGainFadcSamples = new MArrayB(0);//npix*fRunHeader->GetNumSamplesLoGain());
 }
 
@@ -248 +250 @@
     //  Add to lo gains if lflag = 1
     //
-
     MArrayS *arrpix = lflag ? fLoGainPixId       : fHiGainPixId;
     MArrayB *arrsam = lflag ? fLoGainFadcSamples : fHiGainFadcSamples;
@@ -286 +287 @@
 }
 
-/*
-void MRawEvtData::AddPixelLo(UShort_t nOfPixel, TArrayC *data, int nr, int pos)
-{
-    //
-    // add the number of the new pixel to the array as last entry
-    //
-    fLoGainPixId->AddAt(nOfPixel, nr);
-
-    //
-    // add the new slices as last entries to array
-    //
-    fLoGainFadcSamples->AddAt((Byte_t*)data->GetArray(), pos,  data->GetSize());
-}
-
-void MRawEvtData::AddPixelHi(UShort_t nOfPixel, TArrayC *data, int nr, int pos)
-{
-    // check whether we got the right number of new samples
-    // if there are no samples already stored: this is the new number of samples
-    //
-    const Byte_t ns = data->GetSize();
-
-    //
-    // add the number of the new pixel to the array as last entry
-    //
-    fHiGainPixId->AddAt(nOfPixel, nr);
-
-    //
-    // add the new slices as last entries to array
-    //
-    fHiGainFadcSamples->AddAt((Byte_t*)data->GetArray(), arrsam->GetSize()-ns, ns);
-}
-*/
-void MRawEvtData::ReadEvt(ifstream &fin)
+void MRawEvtData::ReadEvt(istream &fin)
 {
     //
@@ -327 +296 @@
     const UShort_t nhi = fRunHeader->GetNumSamplesHiGain();
 
-    TArrayC lo(nlo);
-    TArrayC hi(nhi);
-
     const UShort_t npic = fRunHeader->GetNumPixInCrate();
+
+    //
+    // Enhance array by the size which we'll read now
+    //
+    Byte_t lo[nlo];
+    Byte_t hi[nhi];
+
+    const int npixhi = fHiGainPixId->GetSize();
+    const int npixlo = fLoGainPixId->GetSize();
+
+    fHiGainPixId->Set(npixhi+npic);
+    fLoGainPixId->Set(npixlo+npic);
+
+    const int nsamhi = fHiGainFadcSamples->GetSize();
+    const int nsamlo = fLoGainFadcSamples->GetSize();
+
+    fHiGainFadcSamples->Set(nsamhi+nhi*npic);
+    fLoGainFadcSamples->Set(nsamlo+nlo*npic);
 
     for (int i=0; i<npic; i++)
@@ -339 +323 @@
         const UShort_t npix = fRunHeader->GetPixAssignment(i);
 
-        fin.read((Byte_t*)hi.GetArray(), nhi);
-        AddPixel(npix, &hi, kFALSE);
+        //
+        //  This is to fill the data of one pixel to the MRawEvtHeader Class.
+        //  The parameters are the pixelnumber and the FADC_SLICES values of ADCs
+        //  Add to lo gains if lflag = 1
+        //
+        fin.read(hi, nhi);
+        fHiGainPixId->AddAt(npix, npixhi+i);
+        fHiGainFadcSamples->AddAt(hi, nsamhi+i*nhi, nhi);
 
         // FIXME: Not implemented in the raw files yet
         //if (IsLoGainOn(i, j))
         //{
-        fin.read((Byte_t*)lo.GetArray(), nlo);
-        AddPixel(npix, &lo, kTRUE);
+        fin.read(lo, nlo);
+        fLoGainPixId->AddAt(npix, npixlo+i);
+        fLoGainFadcSamples->AddAt(lo, nsamlo+i*nlo, nlo);
         //}
     }

trunk/MagicSoft/Mars/mraw/MRawEvtData.h

@@ -27 +27 @@
     MArrayB *fLoGainFadcSamples;  // list of lo gain samples of all pixels (ordering: see fLoGainPixId)
 
-    void InitArrays();
+    void InitArrays(Bool_t flag=kFALSE);
     void DeleteArrays();
 
@@ -47 +47 @@
     void Draw (Option_t * = NULL);
 
-    void DeletePixels();
+    void DeletePixels(Bool_t flag=kFALSE);
     void AddPixel(UShort_t nOfPixel, TArrayC *data, Bool_t lflag=kFALSE);
 
@@ -54 +54 @@
     Byte_t GetNumPixels() const;
 
-    void ReadEvt(ifstream &fin);
+    void ReadEvt(istream &fin);
 
     ClassDef(MRawEvtData, 1) //Container to store the raw Event Data

trunk/MagicSoft/Mars/mraw/MRawEvtHeader.cc

@@ -169 +169 @@
     *fLog << "Number of Lo Gains On: " << fNumLoGainOn << endl;
 
-    for (int i=0; i<fPixLoGainOn->GetSize(); i++)
+    for (unsigned int i=0; i<fPixLoGainOn->GetSize(); i++)
     {
         for (int j=0; j<8; j++)
@@ -194 +194 @@
 }
 
-int MRawEvtHeader::ReadEvt(ifstream &fin)
+int MRawEvtHeader::ReadEvt(istream &fin)
 {
     //
@@ -220 +220 @@
 
     fNumLoGainOn = 0;
-    for (int i=0; i<fPixLoGainOn->GetSize(); i++)
+    for (unsigned int i=0; i<fPixLoGainOn->GetSize(); i++)
         for (int j=0; j<8; j++)
             if ((*fPixLoGainOn)[i] & (1<<j))

trunk/MagicSoft/Mars/mraw/MRawEvtHeader.h

@@ -46 +46 @@
   UShort_t GetTrigType() const { return fTrigType; }
 
-  int ReadEvt(ifstream& fin);
+  int ReadEvt(istream& fin);
 
   ClassDef(MRawEvtHeader, 1) // Parameter Conatiner for raw EVENT HEADER

trunk/MagicSoft/Mars/mraw/MRawFileRead.cc

@@ -24 +24 @@
 ClassImp(MRawFileRead)
 
+/*
+  ----------- please don't delete ------------
+#define kBUFSZ 1024
+
+class bifstream : public istream, public streambuf
+{
+private:
+    char fBuffer[kBUFSZ]; //!
+    FILE *fd;
+
+    int sync()
+    {
+        memset(fBuffer, 0, kBUFSZ);
+        return 0; 
+    }
+    int underflow()
+    {
+        int sz=fread(fBuffer, 1, kBUFSZ, fd);
+        setg(fBuffer, fBuffer, fBuffer+kBUFSZ);
+
+        return sz==kBUFSZ ? *(unsigned char*)fBuffer : EOF;//EOF;
+    }
+public:
+    bifstream(const char *name) : istream(this)
+    {
+        fd = fopen(name, "rb");
+        setbuf(fBuffer, kBUFSZ);
+    }
+};
+*/
+
 MRawFileRead::MRawFileRead(const char *fname, const char *name, const char *title)
 {
@@ -38 +69 @@
 }
 
-Bool_t MRawFileRead::PreProcess (MParList *pList)
+MRawFileRead::~MRawFileRead()
+{
+    delete fIn;
+}
+
+Bool_t MRawFileRead::PreProcess(MParList *pList)
 {
     //
     // remember the pointer to the parameter list fur further usage
     //
-    pParList = pList;
 
     //
@@ -103 +138 @@
     return kTRUE;
 }
-    
+
 Bool_t MRawFileRead::Process()
 {
@@ -112 +147 @@
     if (!fRawEvtHeader->ReadEvt(*fIn))
         return kFALSE;
-    //fRawEvtHeader->Print();
 
     //
@@ -135 +169 @@
     {
         fRawCrateArray->GetEntry(i)->ReadEvt(*fIn);
-        //fRawCrateArray->GetEntry(i)->Print();
 
         fRawEvtData->ReadEvt(*fIn);
     }
-    //fRawEvtData->Print();
 
     return kTRUE;
-
 }
 

trunk/MagicSoft/Mars/mraw/MRawFileRead.h

@@ -18 +18 @@
 {
 private:
-    MParList *pParList;
+    MRawRunHeader  *fRawRunHeader;  // run header information container to fill from file
+    MRawEvtHeader  *fRawEvtHeader;  // event header information container to fill from file
+    MRawEvtData    *fRawEvtData;    // raw evt header infomation container to fill from file
+    MRawCrateArray *fRawCrateArray; // crate information array container to fill from file
+    MTime          *fRawEvtTime;    // raw evt time information container to fill from file
 
-    MRawRunHeader  *fRawRunHeader;
-    MRawEvtHeader  *fRawEvtHeader;
-    MRawEvtData    *fRawEvtData;
-    MRawCrateArray *fRawCrateArray;
-    MTime          *fRawEvtTime;
-
-    ifstream *fIn; //!
+    ifstream      *fIn;            //! buffered input stream (file to read from)
 
 public:
     MRawFileRead(const char *filename, const char *name=NULL, const char *title=NULL);
+    ~MRawFileRead();
 
     Bool_t PreProcess(MParList *pList);

trunk/MagicSoft/Mars/mraw/MRawRunHeader.cc

@@ -33 +33 @@
 }
 
-void MRawRunHeader::ReadEvt(ifstream& fin)
+void MRawRunHeader::ReadEvt(istream& fin)
 {
     //
@@ -87 +87 @@
     // print run header information on screen
     //
+    *fLog << endl;
     *fLog << "MagicNumber:  0x" << hex << fMagicNumber << " - " << (fMagicNumber==0xc0c0?"OK":"Wrong!") << endl;
     *fLog << "Version:      " << dec << "Format=" << fFormatVersion << "  ";

trunk/MagicSoft/Mars/mraw/MRawRunHeader.h

@@ -11 +11 @@
 #endif
 
+class TBuffer;
 class ifstream;
 class MArrayS;
@@ -77 +78 @@
     void Print(Option_t *t=NULL);
 
-    void ReadEvt(ifstream& fin);
+    void ReadEvt(istream& fin);
 
     ClassDef(MRawRunHeader, 1)  // storage container for general info