Changeset 1600 for trunk/MagicSoft/Mars/mfileio/MCT1ReadPreProc.cc

Timestamp:

11/13/02 17:03:19 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mfileio/MCT1ReadPreProc.cc (modified) (21 diffs)

trunk/MagicSoft/Mars/mfileio/MCT1ReadPreProc.cc

@@ -42 +42 @@
 //
 /////////////////////////////////////////////////////////////////////////////
-
 #include "MCT1ReadPreProc.h"
 
@@ -81 +80 @@
 //
 MCT1ReadPreProc::MCT1ReadPreProc(const char *fname, const char *name,
-                                 const char *title) : fIn(NULL)
-{
-    fName  = name  ? name  : "MCT1ReadPreProc";
+                                 const char *title) : fIn(NULL), fEntries(0)
+{
+    fName  = name  ? name  : "MRead";
     fTitle = title ? title : "Task to loop over events in CT1 ascii file";
 
@@ -113 +112 @@
 void MCT1ReadPreProc::AddFile(const char *txt)
 {
-    TNamed *name = new TNamed(txt, "");
-    fFileNames->AddLast(name);
+    ifstream *store = fIn;
+
+    fIn = new ifstream(gSystem->ExpandPathName(txt));
+
+    if (!(*fIn))
+    {
+        *fLog << warn << "Cannot open file '" << txt << "'... ignored." << endl;
+        fIn = store;
+        return;
+    }
+
+    fEntries += GetNumEvents();
+
+    delete fIn;
+
+    fIn = store;
+
+    fFileNames->AddLast(new TNamed(txt, ""));
 }
 
@@ -123 +138 @@
 // data members
 //
-void MCT1ReadPreProc::ProcessHeader(const struct outputpars &outpars)
-{
-    fNumPixels = outpars.inumpixels; // number of pixels in the camera
-
-    if (fNumPixels>iMAXNUMPIX || fNumPixels==0)
-        fNumPixels=iMAXNUMPIX;
+void MCT1ReadPreProc::ProcessRunHeader(const struct outputpars &outpars)
+{
+    if (outpars.inumpixels != iMAXNUMPIX)
+        *fLog << warn << "WARNING! File doesn't contain " << iMAXNUMPIX << " Pixels... maybe corrupt." << endl;
+
+    fNumEventsInRun = 0;
 
     //
@@ -166 +181 @@
     *fLog << "Source:    RA=" << outpars.dsourcera_hours << "h  DEC=";
     *fLog << outpars.dsourcedec_deg << "deg" << endl;
-
-    *fLog << "Pixels:    " << fNumPixels << endl;
 
     //int     inummuonpixels;   // number of pixels in the muon shield
@@ -262 +275 @@
 
     // FIXME? Is x-y echanged between Mars CT1 geometry and CT1 definition?
-    fSrcPos->SetXY(outpars.fypointcorr_deg/fGeom->GetConvMm2Deg(),
-                   outpars.fxpointcorr_deg/fGeom->GetConvMm2Deg());
+    fSrcPos->SetXY(-outpars.fypointcorr_deg/fGeom->GetConvMm2Deg(),
+                   -outpars.fxpointcorr_deg/fGeom->GetConvMm2Deg());
+    fSrcPos->SetReadyToSave();
 
     /*
@@ -342 +356 @@
 
     fIsMcFile = outpars.bmontecarlo==TRUE;
+
+    fPedest->SetReadyToSave();
 }
 
@@ -348 +364 @@
 // Read CT1 PreProc File Header:
 //
-Bool_t MCT1ReadPreProc::ReadHeader()
+Bool_t MCT1ReadPreProc::ReadRunHeader()
 {
     char cheadertitle[iHEADERTITLELENGTH];
@@ -381 +397 @@
     fIn->read((Byte_t*)&outpars, sizeof(struct outputpars));
 
-    ProcessHeader(outpars);
+    ProcessRunHeader(outpars);
 
     return kTRUE;
 }
 
-Bool_t MCT1ReadPreProc::ReadFooter()
+Bool_t MCT1ReadPreProc::ReadRunFooter()
 {
     char cheadertitle[iHEADERTITLELENGTH];
@@ -400 +416 @@
 
     if (!s.BeginsWith(m(0,p)))
-        return kFALSE;
+    {
+        fIn->seekg(-iHEADERTITLELENGTH, ios::cur);
+        return kFALSE;
+    }
 
     *fLog << inf << cheadertitle << flush;
@@ -423 +442 @@
      */
 
+    fNumFilterEvts += filterres.ifilter_passed_evts;
+    fNumRuns++;
+
+    *fLog << inf << "Read " << fNumEventsInRun << " events from run." << endl;
+
+    if (fNumEventsInRun!=(UInt_t)filterres.ifilter_passed_evts)
+    {
+        *fLog << warn << "WARNING! Number of events in run doesn't match number of read events." << endl;
+        *fLog << "         File might be corrupt." << endl;
+    }
+
     return kTRUE;
 }
@@ -468 +498 @@
     *fLog << inf << "-----------------------------------------------------------------------" << endl;
 
-    //*fLog << "File contains " << GetNumEvents() << " events." << endl;
-
-    Bool_t rc = ReadHeader();
+    *fLog << "File contains " << GetNumEvents() << " events." << endl;
+
+    // WORKAROUND for not working seekg(0) in GetNumEvents
+    fIn->close();
+    fIn->open(gSystem->ExpandPathName(name));
+
+    Bool_t rc = ReadRunHeader();
 
     if (!rc)
-        *fLog << warn << "Unable to read header... skipping file." << endl;
+        *fLog << warn << "Unable to read first run header... skipping file." << endl;
 
     return rc;
@@ -480 +514 @@
 Int_t MCT1ReadPreProc::GetNumEvents()
 {
-    Int_t n = 0;
+    *fLog << inf << "Scanning file for size" << flush;
+
+    const TString m(cEND_EVENTS_TEMPLATE);
+    const Int_t p = m.First('%');
+    const TString test = m(0, p);
+
+    Int_t nevts = 0;
+    Int_t nruns = 0;
 
     while (!fIn->eof())
     {
-        if (fIn->get()!=cEND_EVENTS_TEMPLATE[0])
-            continue;
-
-        char cheadertitle[iHEADERTITLELENGTH];
-        fIn->read(cheadertitle, iHEADERTITLELENGTH-1);
-
-        TString s = cheadertitle;
-        TString m = cEND_EVENTS_TEMPLATE;
-        Int_t p = m.First('%');
-
-        if (!s.BeginsWith(m(1,p-1)))
-            continue;
-
-        int num;
-
-        sscanf(cheadertitle, &cEND_EVENTS_TEMPLATE[1], &num);
-
-        n += num;
-    }
-
-    fIn->seekg(0);
-
-    return n;
+        fIn->seekg(iHEADERTITLELENGTH, ios::cur);
+        fIn->seekg(sizeof(struct outputpars), ios::cur);
+
+        while (1)
+        {
+            if (fIn->peek()==cEND_EVENTS_TEMPLATE[0])
+            {
+                char cheadertitle[iHEADERTITLELENGTH];
+                fIn->read(cheadertitle, iHEADERTITLELENGTH);
+
+                const TString s = cheadertitle;
+                if (s.BeginsWith(test))
+                {
+                    fIn->seekg(sizeof(struct filterresults), ios::cur);
+                    nruns++;
+                    break;
+                }
+
+                fIn->seekg(-iHEADERTITLELENGTH, ios::cur);
+            }
+
+            fIn->seekg(sizeof(struct eventrecord), ios::cur);
+            if (fIn->eof())
+                break;
+            nevts++;
+        }
+        *fLog << "." << flush;
+    }
+
+    *fLog << "done." << endl;
+    *fLog << "Found " << nevts << " events in " << nruns << " runs." << endl;
+
+    return nevts;
 }
 
@@ -542 +593 @@
     //  look for the source position in the camera
     //
-    fSrcPos = (MSrcPosCam*)pList->FindCreateObj("MSrcPosCam");
+    fSrcPos = (MSrcPosCam*)pList->FindCreateObj("MSrcPosCam", "Source");
     if (!fSrcPos)
         return kFALSE;
@@ -566 +617 @@
     if (!fMcTrig)
         return kFALSE;
+
+    fNumFilterEvts = 0;
+    fNumRuns       = 0;
 
     return kTRUE;
@@ -582 +636 @@
      int   isecfrac_200ns;     // fractional part of isecs_since_midday
      short snot_ok_flags;      // the bits in these two bytes are flags for additional information on the event: Everything OK =: all Bits = 0
-     int   ipreproc_alt_arcs;  // "should be" alt according to preproc (arcseconds)
-     int   ipreproc_az_arcs;   // "should be" az according to preproc (arcseconds)
 
      // for ALT-AZ mount telescopes: rotation angle of the field of
@@ -602 +654 @@
      float fhourangle;
      */
-    /*
-     *fLog << event.isecs_since_midday << "s ";
-     *fLog << event.ipreproc_alt_arcs << "s ";
-     *fLog << event.ipreproc_az_arcs << "s ";
-     *fLog << event.ifieldrot_arcs << "s ";
-     *fLog << event.fhourangle << endl;
-     */
 
     //
@@ -614 +659 @@
     // too the list with it's id, number of photons and error
     //
-    fNphot->InitSize(fNumPixels);
+    fNphot->InitSize(iMAXNUMPIX);
 
     // number of photoelectrons measured in each pixel only the
     // actual number of pixels (outputpars.inumpixels) is written out
     // short spixsig_10thphot[iMAXNUMPIX];
-    for (Int_t i=0; i<fNumPixels; i++)
-    {
-        if (event.spixsig_10thphot[i]>0)
-            fNphot->AddPixel(i, 0.1*event.spixsig_10thphot[i],
-                             (*fPedest)[i].GetMeanRms());
-    }
+    for (Int_t i=0; i<iMAXNUMPIX; i++)
+    {
+        if (event.spixsig_10thphot[i]==0)
+            continue;
+
+        fNphot->AddPixel(i, 0.1*event.spixsig_10thphot[i],
+                         (*fPedest)[i].GetMeanRms());
+    }
+    fNphot->SetReadyToSave();
+
+    // int ipreproc_alt_arcs; // "should be" alt according to preproc (arcseconds)
+    // int ipreproc_az_arcs;  // "should be" az according to preproc (arcseconds)
 
     fMcEvt->SetTheta(TMath::Pi()*(0.5-1./180*event.ialt_arcs/3600));   // altitude (arcseconds)
     fMcEvt->SetPhi(TMath::Pi()/180*event.iaz_arcs/3600);               // azimuth (arcseconds)
+
+    fMcEvt->SetReadyToSave();
 
     if (!fIsMcFile)
@@ -636 +689 @@
     fMcEvt->SetImpact(event.imcimpact_m*100);    // simulated impact
 
-    fMcTrig->SetFirstLevel(event.imctriggerflag); // MC data from Dorota get a triggerflag: 1 means triggered, 0 not. */
+    fMcTrig->SetFirstLevel(event.imctriggerflag);    // MC data from Dorota get a triggerflag: 1 means triggered, 0 not. */
+    fMcTrig->SetReadyToSave();
 
     //float fmcsize_phel;   // Simulated SIZE
@@ -655 +709 @@
 
     //
-    // If the first character isn't the first of the footer it must be
-    // an event
-    //
-    if (fIn->peek()!=cEND_EVENTS_TEMPLATE[0])
-        return kTRUE;
-
-    //
-    // Remember the current position if the first character faked a
-    // footer
-    //
-    const Int_t pos = fIn->tellg();
-
-    //
-    // Try reading the footer. If it isn't succefull jump back...
-    // must be an event
-    //
-    if (!ReadFooter())
-    {
-        fIn->seekg(pos, ios::beg);
-        return kTRUE;
-    }
-
-    *fLog << inf << "Footer found." << endl;
-
-    //
-    // No after reading the footer check if we reached the end of the file
-    //
-    if (fIn->eof())
-    {
-        *fLog << "End of file." << endl;
-        return kFALSE;
-    }
-
-    //
-    // If the eof isn't reached a new header must follow. Check for it.
-    //
-    if (fIn->peek()!=cTITLE_TEMPLATE[0])
-    {
-        *fLog << inf << "Error finding new run header in file (possible EOF)... skipping rest of file." << endl;
-        return kFALSE;
-    }
-
-    *fLog << "-----------------------------------------------------------------------" << endl;
-    if (ReadHeader())
-        return kTRUE;
-
-    *fLog << inf << "Error reading run header in file... skipping rest of file." << endl;
-    return kFALSE;
+    // Because we can have 0-event runs in the file we loop as often
+    // as we don't find a new footer-header combination.
+    //
+    while (1)
+    {
+        //
+        // If the first character isn't the first of the footer it must be
+        // an event
+        //
+        if (fIn->peek()!=cEND_EVENTS_TEMPLATE[0])
+            return kTRUE;
+
+        //
+        // Try reading the footer. If it isn't succefull...
+        // must be an event
+        //
+        if (!ReadRunFooter())
+            return kTRUE;
+
+        *fLog << inf << "Footer found." << endl;
+
+        //
+        // No after reading the footer check if we reached the end of the file
+        //
+        if (fIn->eof())
+        {
+            *fLog << "End of file." << endl;
+            return kFALSE;
+        }
+
+        //
+        // If the eof isn't reached a new header must follow. Check for it.
+        //
+        if (fIn->peek()!=cTITLE_TEMPLATE[0])
+        {
+            *fLog << inf << "Error finding new run header in file (possible EOF)... skipping rest of file." << endl;
+            return kFALSE;
+        }
+
+        *fLog << "-----------------------------------------------------------------------" << endl;
+        if (!ReadRunHeader())
+            return kTRUE;
+    }
 }
 
@@ -726 +775 @@
     struct eventrecord event;
 
-    const int size1 = sizeof(event)-sizeof(event.spixsig_10thphot);
-    const int size2 = sizeof(event.spixsig_10thphot[0])*fNumPixels;
-
-    // read the eventrecord with the recorded number of stored pixels
-    fIn->read((Byte_t*)&event, size1+size2);
-
-    if (fIn->eof())
-        return kFALSE;
+    // read the eventrecord from the file
+    fIn->read((Byte_t*)&event, sizeof(struct eventrecord));
 
     ProcessEvent(event);
 
-    fNumEvent++;
+    fNumEvents++;
+    fNumEventsInRun++;
 
     return kTRUE;
 }
 
+Bool_t MCT1ReadPreProc::PostProcess()
+{
+    *fLog << all;
+    *fLog << "Number events passed the filter: " << fNumFilterEvts << endl;
+    *fLog << "Number of Events read from file: " << fNumEvents << endl;
+    *fLog << "Number of Runs read from file:   " << fNumRuns << endl;
+
+    if (fNumEvents!=fNumFilterEvts)
+    {
+        *fLog << warn << "WARNING! Number of events in file doesn't match number of read events..." << endl;
+        *fLog << "         File might be corrupt." << endl;
+    }
+
+    return kTRUE;
+}