AlgCaldetDigitList.cxx

Go to the documentation of this file.

// $Id: AlgCaldetDigitList.cxx,v 1.1 2003/09/01 15:51:30 tagg Exp $
//
// AlgCaldetDigitList.cxx
//
// AlgCaldetDigitList is a concrete DigitList Algorithm class.
//
// Author:  G. Irwin 4/2000

#include <cassert>

#include "Algorithm/AlgConfig.h"
#include "Algorithm/AlgFactory.h"
#include "Algorithm/AlgHandle.h"
#include "CandData/CandRecord.h"
#include "CandDigit/AlgCaldetDigitList.h"
#include "CandDigit/CandDigit.h"
#include "CandDigit/CandDigitHandle.h"
#include "CandDigit/CandDigitList.h"
#include "CandDigit/CandDigitListHandle.h"
#include "Candidate/CandContext.h"
#include "Conventions/Munits.h"
#include "MessageService/MsgService.h"
#include "Plex/PlexSEIdAltL.h"
#include "Plex/PlexSEIdAltLItem.h"
#include "RawData/RawDigit.h"
#include "RawData/RawRecord.h"
#include "RawData/RawDigitDataBlock.h"
#include "Validity/VldContext.h"
#include "Validity/VldTimeStamp.h"
#include "Plex/PlexHandle.h"
#include "RawData/RawHeader.h"
#include "CaldetRawDigitWrapper.h"
ClassImp(AlgCaldetDigitList)

//______________________________________________________________________
CVSID("$Id: AlgCaldetDigitList.cxx,v 1.1 2003/09/01 15:51:30 tagg Exp $");

//______________________________________________________________________
AlgCaldetDigitList::AlgCaldetDigitList()
{
}

//______________________________________________________________________
AlgCaldetDigitList::~AlgCaldetDigitList()
{
}

//______________________________________________________________________
void AlgCaldetDigitList::RunAlg(AlgConfig &ac, CandHandle &ch,
                                                        CandContext &cx)
{
   assert(cx.GetDataIn());
   
   // Get VldContext
   const VldContext &vldc = *(cx.GetCandRecord()->GetVldContext());
   
   // Iterate over RawDigitDataBlock or CandDigitList daughters.
   if (cx.GetDataIn()->InheritsFrom("RawRecord")) {
     
     PlexHandle ph(vldc);
     
     // Set up to iterate over input RawDigit's to fill output CandDigitList.
     
     // Get Singleton instance of AlgFactory.
     AlgFactory &af = AlgFactory::GetInstance();
     
     // Get correct Algorithm by name from AlgFactory
     const char *tmpcs = 0;
     const char *digitalgorithm = 0;
     const char *digitalgconfig = 0;

     if (ac.Get("DigitAlgorithm", tmpcs)) digitalgorithm = tmpcs;
     if (ac.Get("DigitAlgConfig", tmpcs)) digitalgconfig = tmpcs;

     // Get an AlgHandle to an Algorithm with requested AlgConfig.
     AlgHandle ah = af.GetAlgHandle(digitalgorithm, digitalgconfig);

     CandContext cxx(this, cx.GetMom());

     // Set CandRecord pointer in CandContext.
     cxx.SetCandRecord(cx.GetCandRecord());

     Bool_t isSparse = true;              // sparsified unless otherwise

     // Iterate over RawRecord's RawDataBlockList to find RawDigitDataBlock.
     RawRecord *rr = (RawRecord *) (cx.GetDataIn());
     TIter rdbit = rr->GetRawBlockIter();
     TObject *tob;
     while ((tob = rdbit())) {
       if (tob->InheritsFrom("RawDigitDataBlock")) {
         RawDigitDataBlock *rdb = (RawDigitDataBlock *) tob;
         TIter rdit = rdb->GetDatumIter();
         RawDigit *rd;
         Int_t rawdigitindex = -1;

         Int_t ndigits = rdb->GetNumberOfDigits();
         Bool_t blkIsSparse = 
            ( rdb->GetNumberOfCrates() == rdb->GetNSparsModeCrates());
         isSparse &= blkIsSparse;
         if (!blkIsSparse) 
           MSG("AlgCaldetDigitList", Msg::kInfo) 
         << "AlgCaldetDigitList finds RawDigitDataBlock is Unsparsified with " 
                               << ndigits << " digits " << endl << endl;


         // This object will contain our data.
         CaldetRawDigitWrapper wrapper;

         // Loop through, looking for a timing fiducial hit.
         rdit.Reset();   
         rawdigitindex = -1;
         while ((rd = (RawDigit *) rdit())) { 
           ++rawdigitindex;
           
           if (ph.GetReadoutType(rd->GetChannel()) == ReadoutType::kTimingFid) {
             wrapper.fFiducialDigit = rd;
           }
         }

         if(wrapper.fFiducialDigit==0) {
           MSG("AlgCaldetDigitList",Msg::kInfo)
             << "Event has no fiducial digit." << endl;
         }

         // Loop through again, getting the actual hits we want.
         // skip all non-scint readouts
         rdit.Reset();   
         rawdigitindex = -1;
         while ((rd = (RawDigit *) rdit())) { 
           ++rawdigitindex;
           
           // Make canddigits out of everything.
           //if (ph.GetReadoutType(rd->GetChannel()) != 
           //    ReadoutType::kScintStrip) continue;

           wrapper.fDigit = rd;

           cxx.SetDataIn(&wrapper);
           CandDigitHandle cdh = CandDigit::MakeCandidate(ah, cxx);
           cdh.SetRawDigitIndex(rawdigitindex);

           ch.AddDaughterLink(cdh, kFALSE);      // Don't check for dups
         }
       }
     }

     // Set absolute time offset in CandDigitList.
     // It's no longer subtracted from CandDigits in AlgDigit.
     CandDigitListHandle &cdlh =
       dynamic_cast<CandDigitListHandle &>(ch);
     cdlh.SetAbsTime(((Double_t) vldc.GetTimeStamp().GetNanoSec()) *
                     Munits::nanosecond);
     
     // Set CandDigitList sparsification flag.
     cdlh.SetIsSparse(isSparse);
   }
   
   else if (cx.GetDataIn()->InheritsFrom("TCollection")) {
     TCollection *cdl = (TCollection *) cx.GetDataIn();
     TIter cdit(cdl);
     CandDigitHandle *cdh;
     while ((cdh = (CandDigitHandle *) cdit())) {
       ch.AddDaughterLink(*cdh, kFALSE);         // Don't check for dups
     }
   }
}

//______________________________________________________________________
void AlgCaldetDigitList::Trace(const char * /* c */) const
{
}

---