AlgDeMuxDigitList.cxx

Go to the documentation of this file.

// $Id: AlgDeMuxDigitList.cxx,v 1.5 2004/02/11 20:37:14 gmieg Exp $
//
// AlgDeMuxDigitList.cxx
//
// Concrete Algorithm for filling CandDeMuxDigitList from CandDigitList
//
// Author:  G. Irwin 9/2002
// AlgConfig parameters (with default values) are:
//
// "DeMuxAlgorithm" ("AlgDeMuxCosmics") - Demuxing Algorithm to call
// "DeMuxAlgConfig" ("default") - Demuxing Algorithm's AlgConfig "name"
// "NormalizeWeights" (1) - Normalize weights to 1. if non-zero
// "TrimHyps" (1) - Drop "0" weights if neg., or keep top N (=TrimHyps)

#include <cassert>

#include "Algorithm/AlgConfig.h"
#include "Algorithm/AlgFactory.h"
#include "Algorithm/AlgHandle.h"
#include "CandDigit/AlgDeMuxDigitList.h"
#include "CandDigit/CandDeMuxDigit.h"
#include "CandDigit/CandDeMuxDigitHandle.h"
#include "CandDigit/CandDeMuxDigitList.h"
#include "CandDigit/CandDeMuxDigitListHandle.h"
#include "CandDigit/CandDigit.h"
#include "CandDigit/CandDigitHandle.h"
#include "CandDigit/CandDigitList.h"
#include "CandDigit/CandDigitListHandle.h"
#include "Candidate/CandContext.h"
#include "MessageService/MsgService.h"

ClassImp(AlgDeMuxDigitList)

//______________________________________________________________________
CVSID("$Id: AlgDeMuxDigitList.cxx,v 1.5 2004/02/11 20:37:14 gmieg Exp $");

//______________________________________________________________________
AlgDeMuxDigitList::AlgDeMuxDigitList()
{
}

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

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

// Check for input CandDigitListHandle in CandContext
   assert(cx.GetDataIn());

// Obtain CandContext's CandDigitList as input
   const CandDigitListHandle *cdlh =
             dynamic_cast<const CandDigitListHandle *> (cx.GetDataIn());
   assert(cdlh);   // Check for input CandDigitListHandle in CandContext

// Set CandRecord pointer in CandDeMuxDigitList to that of CandDigitList
   ch.SetCandRecord((cdlh->GetCandRecord())); //Sub-algorithms need this

// Copy member variables from orig CandDigitList to CandDeMuxDigitList
   CandDeMuxDigitListHandle &cddlh =
                           dynamic_cast<CandDeMuxDigitListHandle &>(ch);
   cddlh.SetAbsTime(cdlh->GetAbsTime());
   cddlh.SetIsSparse(cdlh->GetIsSparse());

// Iterate over CandContext's CandDigitList to find CandDigitHandles.
   TIter cdhiter = cdlh->GetDaughterIterator();
   TObject *tob;
   while ((tob = cdhiter())) {
     CandDigitHandle *cdh = dynamic_cast<CandDigitHandle *>(tob);
     if (cdh == 0) continue;   // Skip daughter if not a CandDigitHandle

// Effect copy from CandDigit source to a new CandDeMuxDigit
     CandDeMuxDigitHandle cddh =
                            CandDeMuxDigitHandle::DupCandFromBase(*cdh);
     ch.AddDaughterLink(cddh, kFALSE);           // Don't check for dups
   }

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

// Get correct Algorithm by name from AlgFactory
   const char *tmpcs = 0;
   const char *demuxalgorithm = 0;
   const char *demuxalgconfig = 0;

   if (ac.Get("DeMuxAlgorithm", tmpcs)) demuxalgorithm = tmpcs;
   if (ac.Get("DeMuxAlgConfig", tmpcs)) demuxalgconfig = tmpcs;
   ah = af.GetAlgHandle(demuxalgorithm, demuxalgconfig);

   CandContext cxx(this, cx.GetMom());
   
   ah.RunAlg(ch, cxx);

// Normalize and/or de-sparsify PlexSEIdAltL, if requested
   if (cddlh.GetDeMuxDigitListFlagWord() == 0) {    // Demuxing is valid
     Int_t tmpi = 0;
     Int_t normwts = 0;  Int_t trimhyps = 0;
     if (ac.Get("NormalizeWeights", tmpi)) normwts = tmpi;
     if (ac.Get("TrimHyps", tmpi)) trimhyps = tmpi;
     if (normwts || trimhyps)  {     // Need to loop over DeMuxDigitList
       TIter cddhiter = ch.GetDaughterIterator();
       while ((tob = cddhiter())) {
         CandDeMuxDigitHandle *cddh =
                              dynamic_cast<CandDeMuxDigitHandle *>(tob);
         if (cddh == 0) continue;      // Require a CandDeMuxDigitHandle
         if (cddh->GetPlexSEIdAltL().GetDemuxVetoFlag() != 0) continue;
         PlexSEIdAltL &seidw = cddh->GetPlexSEIdAltLWritable();
         if (normwts) seidw.NormalizeWeights(); // Norm wts to 1 if ne 0
         if (trimhyps<0) seidw.DropZeroWeights();   // Drop 0 wts if neg
         else if (trimhyps>0) seidw.KeepTopWeights(trimhyps);  // Keep N
       }
     }
   }
}

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

---