AlgMCDigitList.cxx

Go to the documentation of this file.

// $Id: AlgMCDigitList.cxx,v 1.12 2003/03/29 01:26:13 gmieg Exp $
//
// AlgMCDigitList.cxx
//
// AlgMCDigitList is a concrete MCDigitList 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/AlgMCDigitList.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"

ClassImp(AlgMCDigitList)

//______________________________________________________________________
CVSID("$Id: AlgMCDigitList.cxx,v 1.12 2003/03/29 01:26:13 gmieg Exp $");

//______________________________________________________________________
AlgMCDigitList::AlgMCDigitList()
{
}

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

//______________________________________________________________________
void AlgMCDigitList::RunAlg(AlgConfig &ac, CandHandle &ch,
                                                        CandContext &cx)
{
   MSG("AlgMCDigitList", Msg::kDebug)
                         << "Starting AlgMCDigitList::RunAlg()" << endl;

   assert(cx.GetDataIn());

// Get VldContext
   const VldContext &vldc = *(cx.GetCandRecord()->GetVldContext());

// Iterate over RawDigitDataBlock or CandDigitList daughters.
   if (cx.GetDataIn()->InheritsFrom("RawRecord")) {

// Set up to iterate over input RawDigit's to fill output CandDigitList.

// Get Singleton instance of AlgFactory.
     MSG("AlgMCDigitList", Msg::kDebug)
               << "AlgFactory &af = AlgFactory::GetInstance();" << endl;
     AlgFactory &af = AlgFactory::GetInstance();

// Get correct Algorithm by name from AlgFactory
     const char *tmpcs = 0;
     const char *mcdigitalgorithm = 0;
     const char *mcdigitalgconfig = 0;

     if (ac.Get("MCDigitAlgorithm", tmpcs)) mcdigitalgorithm = tmpcs;
     if (ac.Get("MCDigitAlgConfig", tmpcs)) mcdigitalgconfig = tmpcs;

// Get an AlgHandle to an Algorithm with requested AlgConfig.
     MSG("AlgMCDigitList", Msg::kDebug)
         << "AlgHandle ah = af.GetAlgHandle("
         << mcdigitalgorithm << ", "
         << mcdigitalgconfig << ");" << endl;
     AlgHandle ah = af.GetAlgHandle(mcdigitalgorithm, mcdigitalgconfig);

     MSG("AlgMCDigitList", Msg::kDebug)
                               << "Create CandContext instance." << endl
                                    << "CandContext cxx(this);" << endl;
     CandContext cxx(this, cx.GetMom());

// Set CandRecord pointer in CandContext.
     MSG("AlgMCDigitList", Msg::kDebug)
                    << "Set CandRecord pointer in CandContext." << endl;
     cxx.SetCandRecord(cx.GetCandRecord());

// 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;
         while ((rd = (RawDigit *) rdit())) { 
           ++rawdigitindex;

           MSG("AlgMCDigitList", Msg::kDebug)
             << "Set RawDigit object as CandContext input data." << endl
                                        << "cxx.SetDataIn(rd);" << endl;
           cxx.SetDataIn(rd);

           MSG("AlgMCDigitList", Msg::kDebug)
                  << "Get handle to CandDigit." << endl
                  << "cdh = CandDigit::MakeCandidate(ah, cxx);" << endl;

           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);
   }

   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 AlgMCDigitList::Trace(const char * /* c */) const
{
}

---