RawLiTpmtDigitsBlock.cxx

Go to the documentation of this file.

// $Id: RawLiTpmtDigitsBlock.cxx,v 1.5 2005/08/26 19:04:53 rhatcher Exp $
// 
// RawLiTpmtDigitsBlock 
// 
// RawLiTpmtDigitsBlock holds LightInjection TriggerPMT hits
//
// Author:  R. Hatcher 2003.06.30
//

#include "RawData/RawLiTpmtDigitsBlock.h"

enum EDebugDigitDataBlock {
   dbg_DumpHeaderOnUnpack = 0x0001,
   dbg_DumpCrateOnUnpack  = 0x0002,
   dbg_DumpDigitOnUnpack  = 0x0004
};
UInt_t RawLiTpmtDigitsBlock::fgDebugFlags = 0;

#include "RawData/RawCrateStatus.h"
#include "RawData/RawVaMCDigit.h"
#include "RawData/RawQieMCDigit.h"

#include "MessageService/MsgStream.h"
#include "MessageService/MsgFormat.h"

#include "MessageService/MsgService.h"
CVSID("$Id: RawLiTpmtDigitsBlock.cxx,v 1.5 2005/08/26 19:04:53 rhatcher Exp $");

#include "RawData/RawBlockRegistry.h"
REGISTERRAWBLOCK(RawLiTpmtDigitsBlock,kMdBlockLiTpmtDigits,0);

ClassImp(RawLiTpmtDigitsBlock)

enum ELiTpmtDigitsBlkPos {
   indx_size       =  0,
   indx_checksum   =  1,
   indx_blockid    =  2,
   indx_run        =  3,
   indx_subrun     =  4,
   indx_runtype    =  indx_subrun,
   indx_loadsec    =  5,
   indx_loadnsec   =  6,
   indx_tf_v1      =  7,
   indx_crateid_v1 =  8,
   indx_crateid    =  7,
   indx_tf         =  8,
   indx_ndigits    =  9,
   zzzz_last       = 10
};

//_____________________________________________________________________________
RawLiTpmtDigitsBlock::RawLiTpmtDigitsBlock()
  : fRawCrateStatus(0), fRawDigits(0)
{
   // Default constructor
}

//_____________________________________________________________________________
RawLiTpmtDigitsBlock::RawLiTpmtDigitsBlock(const Int_t *block)
  : RawDataBlock(block), fRawCrateStatus(0), fRawDigits(0)
{
   //  stored block format is:
   //---------------------
   //  0   # words in block           
   //  1   checksum                   
   //  2   Block Id                   
   //-----                                                      
   //  3   Run Number                 
   //  4   SubRun | RunType           
   //  5   Time (sec)                 
   //  6   Time (ns)                  
   //  7   CrateId (+ special bits)   TimeFrame # (v1)
   //  8   TimeFrame #                CrateId     (v1)
   //  8   # of digits to follow      
   //  ... packed two 32-bit words per digit

}

//_____________________________________________________________________________
RawLiTpmtDigitsBlock::~RawLiTpmtDigitsBlock()
{
   // dtor
   if ( fRawDigits ) {
     fRawDigits->Delete();
     delete fRawDigits;
     fRawDigits = 0;
   }
   if ( fRawCrateStatus ) {
     delete fRawCrateStatus;
     fRawCrateStatus = 0;
   }

}

//____________________________________________________________________________
RawLiTpmtDigitsBlock& RawLiTpmtDigitsBlock::operator=(const RawLiTpmtDigitsBlock& rhs)
{
   // deep copy assignment 
  if (this != &rhs) {
    RawDataBlock::operator=(rhs);
    if (fRawCrateStatus) { delete fRawCrateStatus; fRawCrateStatus = 0; }
    if (fRawDigits) { 
      fRawDigits->Delete(); 
      delete fRawDigits; 
      fRawDigits = 0;
    }
  }
  return *this;
}

//_____________________________________________________________________________
std::ostream& RawLiTpmtDigitsBlock::FormatToOStream(std::ostream& os, 
                                                    Option_t *option) const
{
   RawDataBlock::FormatToOStream(os,option);
   if (option[0] == 'X') return os;
   
   // additional block specific formatted output is done here

   os << " TimeStamp " << GetTimeStamp().AsString("c") << endl;
   os << " Run " << GetRun() 
      << " SubRun " << GetSubRun() 
      << " RunType " << GetRunType() 
      << " TimeFrame " << GetTimeFrameNum()
      << " with " << GetNumberOfDigits() << " digits "
      << endl; 
   os << " " << *GetCrateStatus() << endl;

   if (strchr(option,'q') != 0) return os;

   TIter iter = GetDatumIter();
   TObject *tobj;
   RawDigit *udigit;
   while ( ( tobj = iter.Next() ) ) {
     udigit = dynamic_cast<RawDigit *>(tobj);
     if (udigit) os << " " << (*udigit) << endl;
   }
   os << endl;
   return os;
}

//_____________________________________________________________________________
Int_t RawLiTpmtDigitsBlock::GetRun() const
{
   // get the run number
   if (fSize > indx_run) return fRawBlock[indx_run];
   return -1;
}

//_____________________________________________________________________________
Short_t RawLiTpmtDigitsBlock::GetSubRun() const
{
   // get the subrun number
   if (fSize > indx_subrun) return (fRawBlock[indx_subrun]>>16)&0xffff;
   return -1;
}

//_____________________________________________________________________________
Short_t RawLiTpmtDigitsBlock::GetRunType() const
{
   // get the run type
   if (fSize > indx_runtype) return fRawBlock[indx_runtype]&0xffff;
   return -1;
}

//_____________________________________________________________________________
VldTimeStamp RawLiTpmtDigitsBlock::GetTimeStamp() const
{
   // get the start time
   if (fSize >= zzzz_last) 
      return VldTimeStamp(fRawBlock[indx_loadsec],fRawBlock[indx_loadnsec]);

   return VldTimeStamp((time_t)0,(Int_t)0);
}

//_____________________________________________________________________________
VldContext RawLiTpmtDigitsBlock::GetVldContext() const
{
   // get a VldContext

   RawBlockId rbid = GetBlockId();
   return VldContext(rbid.GetDetector(),rbid.GetSimFlag(),GetTimeStamp());
}

//_____________________________________________________________________________
Int_t RawLiTpmtDigitsBlock::GetTimeFrameNum() const
{
   // get the time frame #
   int indx_tf_v = indx_tf;
   if (GetMinorId() <= 1) indx_tf_v = indx_tf_v1;
   if (fSize >= indx_tf_v) return fRawBlock[indx_tf_v];
   return -1;
}

//_____________________________________________________________________________
const RawDigitCrateStatus* RawLiTpmtDigitsBlock::GetCrateStatus() const
{
   // get the crate id --> status
   if ( ! fRawDigits ) FillRawDigits();
   return fRawCrateStatus;

}

//_____________________________________________________________________________
Int_t RawLiTpmtDigitsBlock::GetNumberOfDigits() const
{
   // get the number of digits
   if ( ! fRawDigits ) FillRawDigits();
   return fRawDigits->GetEntriesFast();  // known not to have gaps
}

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

   if ( ! fRawDigits ) FillRawDigits();
   return (RawDigit*) fRawDigits->At(idx);
}

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

   if ( ! fRawDigits ) FillRawDigits();
   return TIter(fRawDigits,dir);
}

//_____________________________________________________________________________
Int_t RawLiTpmtDigitsBlock::IndexOf(RawDigit *rawdigit) const
{
   // create/fill the TObjArray of RawDigits (if it doesn't exist)
   // return the index for this object (if contained in the list)

   if ( ! fRawDigits ) FillRawDigits();
   return fRawDigits->IndexOf(rawdigit);
}

//_____________________________________________________________________________
void RawLiTpmtDigitsBlock::FillRawDigits() const
{
   // create the TObjArray of RawDigit (if it doesn't exist)
   
   if ( fRawDigits ) return;   // already filled
   if ( fSize <= 0 || fRawBlock == 0) {
      MSG("RawData",Msg::kWarning)
         << "RawLiTpmtDigitsBlock::FillRawDigits empty block? "
         << fSize << " " << fRawBlock << endl;
      return;
   }

   fRawDigits          = new TObjArray();
   RawDigit  *rawdigit = 0;
   const Int_t    *p   = fRawBlock;
   const Int_t    *end = fRawBlock + fSize; // 1 beyond end

   int minorId = GetMinorId();
   RawBlockId rbid = GetBlockId();
   // skip #words, blockid and checksum, etc
   p += ( (minorId<=1) ? indx_crateid_v1 : indx_crateid ); 

   if (fgDebugFlags&dbg_DumpHeaderOnUnpack) {
      MSG("RawData",Msg::kInfo) 
         << "RawLiTpmtDigitsBlock::FillRawDigits "
         << " fSize " << fSize 
         << " id " << rbid.AsString()
         << endl;
   }

   // unpack crate info
   // version >1 has TF word between crateId and # of entries
   Bool_t skipTFword = ( minorId > 1 );  
   fRawCrateStatus = 
     new RawDigitCrateStatus(rbid,p,GetTimeStamp(),skipTFword);
   Int_t npairs = fRawCrateStatus->GetEntries(); 
   if (fgDebugFlags&dbg_DumpCrateOnUnpack) fRawCrateStatus->Print();

   while (p<end) {

      for (Int_t ipair = 0; ipair < npairs; ipair++) {
         // assume that the readout type is consistent for whole crate
         // and that the status block knows what it is
         ElecType::Elec_t etype = fRawCrateStatus->GetElecType();

         // pointer is advanced by RawDigit ctor
         // knowing how many words it needed to eat
         switch (rbid.GetSimFlag()) {
         case SimFlag::kData:
         case SimFlag::kDaqFakeData:
            switch (etype) {
            case ElecType::kVA:
               rawdigit = new RawVaDigit(p,fRawCrateStatus);
               break;
            case ElecType::kQIE:
               rawdigit = new RawQieDigit(p,fRawCrateStatus);
               break;
            default:
               rawdigit = new RawDigit(p,fRawCrateStatus);
               break;
            }
            break;
         case SimFlag::kMC:
         case SimFlag::kReroot:
            switch (etype) {
            case ElecType::kVA:
               rawdigit = new RawVaMCDigit(p,fRawCrateStatus);
               break;
            case ElecType::kQIE:
               rawdigit = new RawQieMCDigit(p,fRawCrateStatus);
               break;
            case ElecType::kReroot:
              {
                if (rbid.GetDetector() == Detector::kNear) {
                  static int nmsg = 10;
                  if (nmsg) {
                    MSG("RawData",Msg::kWarning)
                      << "RawLiTpmtDigitsBlock::FillRawDigits SimFlag "
                      << SimFlag::AsString(rbid.GetSimFlag())
                      << " ElecType " << ElecType::AsString(etype)
                      << " in Detector " 
                      << Detector::AsString(rbid.GetDetector())
                      << endl
                      << " treat as kQIE for RawQieMCDigit creation"
                      << endl;
                    if (nmsg) {
                      if (--nmsg == 0)
                        MSG("RawData",Msg::kWarning)
                          << " ... last warning of this type" << endl;
                    }
                  }
                  rawdigit = new RawQieMCDigit(p,fRawCrateStatus);
                  break;
                }
                else 
                  MSG("RawData",Msg::kError) 
                    << "RawLiTpmtDigitsBlock::FillRawDigits SimFlag "
                    << SimFlag::AsString(rbid.GetSimFlag())
                    << " ElecType " << ElecType::AsString(etype)
                    << endl << "  attempt to move past problem " << endl;
                p += 3;
                break;
              }
            default:
               MSG("RawData",Msg::kError) 
                  << "RawLiTpmtDigitsBlock::FillRawDigits SimFlag "
                  << SimFlag::AsString(rbid.GetSimFlag())
                  << " ElecType " << ElecType::AsString(etype)
                  << endl << "  attempt to move past problem " << endl;
               p += 3;
               break;
            }
            break;
         case SimFlag::kUnknown:
               MSG("RawData",Msg::kError) 
                  << "RawLiTpmtDigitsBlock::FillRawDigits SimFlag "
                  << SimFlag::AsString(rbid.GetSimFlag())
                  << endl << "  attempt to move past problem " << endl;
               p += 2;
               break;
         }

         fRawDigits->Add(rawdigit);

         if (fgDebugFlags&dbg_DumpDigitOnUnpack) {
           rawdigit->Print(); cout << endl;
         }
         if (p>end) {
            MSG("RawData",Msg::kWarning)
               << "RawLiTpmtDigitsBlock::FillRawDigits crate info seems trashed"
               << endl << " SimFlag " << SimFlag::AsString(rbid.GetSimFlag())
               << " ElecType " << ElecType::AsString(etype)
               << endl;
            break;
         }
      }
   }

}

//_____________________________________________________________________________

---