RawLITimingSummary.cxx

Go to the documentation of this file.

// $Id: RawLITimingSummary.cxx,v 1.16 2005/08/26 19:04:53 rhatcher Exp $
// 
// RawLITimingSummary
// 
// RawLITimingSummary is the unpacked Timing summary for one channel
//
// Author:  R. Hatcher 2000.04.19
//

#include "RawData/RawLITimingSummary.h"

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

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

ClassImp(RawLITimingSummary)

//_____________________________________________________________________________
std::ostream& operator<<(std::ostream& os, const RawLITimingSummary& s)
{ return s.FormatToOStream(os); }
 
//_____________________________________________________________________________
RawLITimingSummary::RawLITimingSummary()
   : fCrateStatus(0), fRawChannelId(0), fEntries(-1), fIntMean(-1), fIntRms(-1)
{
   // Default constructor

   // MSG("RawData",Msg::kWarning) << "RawLITimingSummary default ctor" << endl;
}

//_____________________________________________________________________________
RawLITimingSummary::RawLITimingSummary(const Int_t *&p, 
                                       const RawCrateStatus* cstat,
                                       Int_t minorVersion)
   : fCrateStatus(cstat), fRawChannelId(0), fEntries(-1), 
     fIntMean(-1), fIntRms(-1)
{
   // this is where the everything is actually unpacked

   // version 0,1,2:
   // 1st 32-bit word: Chadd (16); no of entries (16)   (v0 U/L swapped)
   // 2nd 32-bit word: Mean 
   // 3rd 32-bit word: RMS
   // Mean and RMS are both scaled by a factor of 4

   const Int_t maskShort = 0xffff;

   switch (minorVersion) {
   case 0:
   {
      // note this version had upper/lower swapped from v1, v2...
      Int_t chadd =   p[0]         & maskShort;
      BuildRawChannelId(cstat,chadd);

      fEntries   =  ( p[0] >> 16 ) & maskShort;

      fIntMean   =    p[1];
      fIntRms    =    p[2];

      // advance pointer by amount it eats up
      p += 3;
      break;
   }
   default:
   {
     MAXMSG("RawData",Msg::kWarning,10)
       << "RawLITimingSummary::ctor with minor id "
       << cstat->GetRawBlockId().GetMinor()
       << " unknown layout" << endl;
     // don't break, fall through to latest
   }
   case 1:
   case 2:
   case 3:
   case 4:
   case 5:
   {
      Int_t chadd = ( p[0] >> 16) & maskShort;
      BuildRawChannelId(cstat,chadd);
      
      fEntries   =    p[0]         & maskShort;

      fIntMean   =    p[1];
      fIntRms    =    p[2];

      // advance pointer by amount it eats up
      p += 3;
      break;
   }  
   } // end minor version number switch

}

//_____________________________________________________________________________
RawLITimingSummary::~RawLITimingSummary()
{
   // delete all owned sub-objects

}

//_____________________________________________________________________________
Float_t RawLITimingSummary::GetMean() const
{
   // return unscaled mean
   
   switch (fCrateStatus->GetRawBlockId().GetMinor()) {
   default:
     {
      MAXMSG("RawData",Msg::kWarning,100)
         << "RawLITimingSummary::GetMean() with minor id "
         << fCrateStatus->GetRawBlockId().GetMinor()
         << " assumes default scale factor" << endl;
     }
      // don't break, fall through to latest
   case 5:
     return fIntMean * 0.0625;
     break;
   case 0:
   case 1:
   case 2:
   case 3:
   case 4:
      return fIntMean * 0.25;
      break;
   }
}

//_____________________________________________________________________________
Float_t RawLITimingSummary::GetRms() const
{
   // return unscaled RMS
   
   switch (fCrateStatus->GetRawBlockId().GetMinor()) {
   default:
     {
      MAXMSG("RawData",Msg::kWarning,100)
         << "RawLITimingSummary::GetRms() with minor id "
         << fCrateStatus->GetRawBlockId().GetMinor()
         << " assumes default scale factor" << endl;
     }
      // don't break, fall through to latest
   case 5:
     return fIntRms * 0.0625;
     break;
   case 0:
   case 1:
   case 2:
   case 3:
   case 4:
      return fIntRms * 0.25;
      break;
   }
}

//_____________________________________________________________________________
void RawLITimingSummary::Print(Option_t *option) const
{
   FormatToOStream(cout,option);
}

//_____________________________________________________________________________
std::ostream& RawLITimingSummary::FormatToOStream(std::ostream& os,
                                                  Option_t * /* option */) const
{
   MsgFormat ffmt("%9.2f");
   MsgFormat ifmt("%7d");
   
   os << " "  << fRawChannelId
      << "  " << setw(5)  << fEntries 
      << "  " << ffmt(GetMean())
      << " (" << ifmt(fIntMean) << ")"
      << "  " << ffmt(GetRms())
      << " (" << ifmt(fIntRms) << ")"
      << endl;

   return os;
}

//_____________________________________________________________________________
void RawLITimingSummary::BuildRawChannelId(const RawCrateStatus *cstat, Int_t chadd)
{
   Detector::Detector_t detector = 
                            cstat->GetRawBlockId().GetDetector();
   ElecType::Elec_t etype = cstat->GetElecType();
   Bool_t         pedMode = cstat->GetPedMode();
   Bool_t       sparsMode = cstat->GetSparsMode();
   Bool_t      commonMode = cstat->GetCommonMode();
   Int_t            crate = cstat->GetCrate();

   fRawChannelId = RawChannelId(detector,etype,crate,chadd);
   fRawChannelId.SetModeBits(commonMode,sparsMode,pedMode);
}
//_____________________________________________________________________________

---