RawLIAdcSummaryBlock.cxx

Go to the documentation of this file.

// $Id: RawLIAdcSummaryBlock.cxx,v 1.16 2004/09/02 21:21:23 rhatcher Exp $
// 
// RawLIAdcSummaryBlock 
// 
// RawLIAdcSummaryBlock holds a block of raw light injection summaries
//  and knows how to unpack it into individual RawLIAdcSummary's
//
// Author:  R. Hatcher 2001.05.10
//

#include "RawData/RawLIAdcSummaryBlock.h"

enum EDebugLIAdcSummaryBlock {
   dbg_DumpHeaderOnUnpack = 0x0001,
   dbg_DumpCrateOnUnpack  = 0x0002,
   dbg_DumpAdcSumOnUnpack = 0x0004
};

UInt_t RawLIAdcSummaryBlock::fgDebugFlags = 0;
//  dbg_DumpCrateOnUnpack | dbg_DumpAdcSumOnUnpack;

#include "RawData/RawLIHeader.h"
#include "RawData/RawLIAdcSummary.h"

#include "MessageService/MsgStream.h"
#include "MessageService/MsgService.h"
CVSID("$Id: RawLIAdcSummaryBlock.cxx,v 1.16 2004/09/02 21:21:23 rhatcher Exp $");

#include "RawData/RawBlockRegistry.h"
REGISTERRAWBLOCK(RawLIAdcSummaryBlock,kMdBlockLiSummaryAdc,0);

ClassImp(RawLIAdcSummaryBlock)

//_____________________________________________________________________________
RawLIAdcSummaryBlock::RawLIAdcSummaryBlock() :
  RawDataBlock(), fRawLIHeader(0), fRawCrateStatuses(0), fRawAdcSums(0)
{
   // Default constructor
}

//_____________________________________________________________________________
RawLIAdcSummaryBlock::RawLIAdcSummaryBlock(const Int_t *block)
  : RawDataBlock(block), 
    fRawLIHeader(0), fRawCrateStatuses(0), fRawAdcSums(0)
{
   //  stored block format is:
   //---------------------
   //  0   # words in block
   //  1   checksum
   //  2   Block Id
   //-----
   //  3   CalibType   | CalibPoint#
   //  4   PulserBox # | LED #
   //  5   PulseHeight | PulseWidth
   //  6   # pulses    | flashing period
   //  7   T0 (1st ROP) sec
   //  8   T0 (1st ROP) nsec
   //  9   spare field 0 (v0,v1) : # pulses used in this summary (v2)
   // 10   spare field 1 (v0,v1) : spare field 0  (v2)
   //-----
   // 11   Crate # i (+extra)
   // 12   # of crate entries (N_i)
   //-----
   // 13   ChAdd 1_i  | # of entries
   // 14   Mean*4     | RMS*4
   // 15   ChAdd 2_i  | # of entries
   // 16   Mean*4     | RMS*4
   //  ...
   //      ChAdd N_i  | # of entries
   //      Mean*4     | RMS*4
   //-----
   //  x   Crate # j (+extra)
   //      # of crate entries (N_j)
   //-----
   //      ChAdd 1_j  | # of entries
   //      Mean*4     | RMS*4
   //      ChAdd 2_j  | # of entries
   //      Mean*4     | RMS*4
   //  ...
   //      ChAdd N_j  | # of entries
   //      Mean*4     | RMS*4

   // version 0:
   // 1st 32-bit word: no of entries (16); Chadd (16)
   // 2nd 32-bit word: RMS (16); Mean (16)
   // Mean and RMS are both scaled by a factor of 4, so that they have a
   // two-binal resolution.

   // version 1:
   // 1st 32-bit word: Chadd (16); no of entries (16)
   // 2nd 32-bit word: Mean (16);  RMS (16)
   // Mean and RMS are both scaled by a factor of 4, so that they have a
   // two-binal resolution.

   // version 2:
   // 1st 32-bit word: Chadd (16); no of entries (16)
   // 2nd 32-bit word: Mean (17);  RMS (15)
   // Mean is scaled by a factor of 4 and RMS scaled by 2.
   // Mean is a 2's complement # (ie. when unpacking high order bit
   // of the field must be sign extended).

   // We're now up to v4 - look in RawLIAdcSummary.cxx

}

//_____________________________________________________________________________
RawLIAdcSummaryBlock::~RawLIAdcSummaryBlock()
{

   if (fRawLIHeader) {
      delete fRawLIHeader;
      fRawLIHeader = 0;
   }

   if (fRawAdcSums) {
      fRawAdcSums->Delete();
      delete fRawAdcSums;
      fRawAdcSums = 0;
   }
   if (fRawCrateStatuses) {
      fRawCrateStatuses->Delete();
      delete fRawCrateStatuses;
      fRawCrateStatuses = 0;
   }
   return;
}


//_____________________________________________________________________________
RawLIAdcSummaryBlock& RawLIAdcSummaryBlock::operator=(const RawLIAdcSummaryBlock& rhs)
{

   // deep copy assignment
   if (this != &rhs) {
     RawDataBlock::operator=(rhs);

     if (fRawLIHeader) {
       delete fRawLIHeader;
       fRawLIHeader = 0;
     }
     
     if (fRawAdcSums) {
       fRawAdcSums->Delete();
       delete fRawAdcSums;
       fRawAdcSums = 0;
     }
     if (fRawCrateStatuses) {
       fRawCrateStatuses->Delete();
       delete fRawCrateStatuses;
       fRawCrateStatuses = 0;
     }
   }
   return *this;
}

//_____________________________________________________________________________
Int_t RawLIAdcSummaryBlock::GetNumberOfSummaries(void) const
{
   // create/fill the TObjArray of RawAdcSummarys (if it doesn't exist)
   // return the number of entries

   if ( ! fRawAdcSums ) FillRawAdcSums();
   return fRawAdcSums->GetEntriesFast(); // known not to have gaps
}

//_____________________________________________________________________________
const RawLIAdcSummary* RawLIAdcSummaryBlock::At(Int_t idx) const
{
   // create/fill the TObjArray of RawAdcSummarys (if it doesn't exist)
   // return pointer to the i-th element

   if ( ! fRawAdcSums ) FillRawAdcSums();
   return (RawLIAdcSummary*) fRawAdcSums->At(idx);
}

//_____________________________________________________________________________
TIter RawLIAdcSummaryBlock::GetDatumIter(Bool_t dir) const
{
   // create/fill the TObjArray of RawAdcSummarys (if it doesn't exist)
   // return an iterator to look over them

   if ( ! fRawAdcSums ) FillRawAdcSums();
   return TIter(fRawAdcSums,dir);
}

//_____________________________________________________________________________
Int_t RawLIAdcSummaryBlock::IndexOf(RawLIAdcSummary *rawAdcSum) const
{
   // create/fill the TObjArray of RawAdcSummarys (if it doesn't exist)
   // return the index for this object (if contained in the list)

   if ( ! fRawAdcSums ) FillRawAdcSums();
   return fRawAdcSums->IndexOf(rawAdcSum);
}

//_____________________________________________________________________________
void RawLIAdcSummaryBlock::FillRawAdcSums() const
{
   // create the TObjArray of RawAdcSum (if it doesn't exist)

   if ( fRawAdcSums ) return;   // already filled
   if ( fSize <= 0 || fRawBlock == 0) {
      MSG("RawData",Msg::kWarning)
         << "RawLIAdcSummaryBlock::FillRawAdcSums empty block? "
         << fSize << " " << fRawBlock << endl;
      // build a dummy header at least
      fRawLIHeader = new RawLIHeader;
      return;
   }

   fRawCrateStatuses = new TObjArray();
   fRawAdcSums       = new TObjArray();
   RawCrateStatus  *cratestatus = 0;
   RawLIAdcSummary *rawAdcSum   = 0;
   const Int_t    *p   = fRawBlock;
   const Int_t    *end = fRawBlock + fSize; // 1 beyond end

   RawBlockId rbid = GetBlockId();
   Int_t minor = rbid.GetMinor();

   p += 3; // skip #words, blockid and checksum

   // eat up the common header info
   fRawLIHeader = new RawLIHeader(rbid,p);

   if (fgDebugFlags&dbg_DumpHeaderOnUnpack) fRawLIHeader->Print();

   while (p<end) {
      // start a new crate (pointer advanced by RawCrateStatus ctor
      cratestatus = new RawCrateStatus(rbid,p);
      fRawCrateStatuses->Add(cratestatus);
      Int_t npairs = cratestatus->GetEntries(); 

      if (fgDebugFlags&dbg_DumpCrateOnUnpack) cratestatus->Print();

      for (Int_t ipair = 0; ipair < npairs; ipair++) {
         rawAdcSum = new RawLIAdcSummary(p,cratestatus,minor);
         fRawAdcSums->Add(rawAdcSum);

         if (fgDebugFlags&dbg_DumpAdcSumOnUnpack) {
           std::cout << " i=" << ipair << " ";
           rawAdcSum->Print();
         }
         if (p>end) {
            MSG("RawData",Msg::kWarning)
               << "RawLIAdcSummaryBlock::FillRawAdcSums crate info seems trashed"
               << endl << " SimFlag " << SimFlag::AsString(rbid.GetSimFlag())
               << endl;
            break;
         }
      }
   }

}

//_____________________________________________________________________________
std::ostream& RawLIAdcSummaryBlock::FormatToOStream(std::ostream& os, 
                                                    Option_t *option) const
{
   RawDataBlock::FormatToOStream(os,option);
   if (option[0] == 'X') return os;

   TIter iter = GetDatumIter();

   const RawLIHeader* header = GetHeader();
   os << (*header);

   os << " " << "Channel         " 
      << "entries"
      << " " << "      Mean" << "          "
      << " " << "       RMS" << "          "
      << endl;

   TObject *tobj;
   RawLIAdcSummary *usum;
   while ( ( tobj = iter.Next() ) ) {
      usum = dynamic_cast<RawLIAdcSummary *>(tobj);
      if (usum)  os << (*usum);
   }
   os << endl;

   return os;   

}

//_____________________________________________________________________________

---