AlgFilterChannels.cxx

Go to the documentation of this file.

// $Id: AlgFilterChannels.cxx,v 1.2 2008/11/18 17:09:29 rhatcher Exp $
//
// AlgFilterChannels.cxx
//
// Concrete Algorithm for filtering CandDigit's from CandChannels
//
// Author:  G. Irwin 2/2003
// AlgConfig parameters (with default values) are:
//
// "FilterStrategy" - Filtering strategy to use

#include <cassert>

#include "CalDetDST/AlgFilterChannels.h"
#include "CalDetDST/ChannelsToFilter.h"
#include "DatabaseInterface/DbiResultPtr.h"

#include "Algorithm/AlgConfig.h"
#include "Algorithm/AlgFactory.h"
#include "Algorithm/AlgHandle.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(AlgFilterChannels)

//______________________________________________________________________
CVSID("$Id: AlgFilterChannels.cxx,v 1.2 2008/11/18 17:09:29 rhatcher Exp $");

//______________________________________________________________________
AlgFilterChannels::AlgFilterChannels()
{
}

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

//______________________________________________________________________
void AlgFilterChannels::RunAlg(AlgConfig & /* ac */, CandHandle &ch,
                                                        CandContext &cx)
{
   MSG("AlgFilterChannels", Msg::kDebug)
                     << "Starting AlgFilterChannels::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 new CandDigitList to that of CandDigitList
   ch.SetCandRecord((cdlh->GetCandRecord())); //Sub-algorithms need this

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

// open a connection to the db

/*
   // this bit of code doesn't work!
   // apparently, cx.GetCandRecord() returns a bad pointer
   // as cx.SetCandRecord() never gets called by the module
   // arghh!

   const VldContext* vcp = cx.GetCandRecord()->GetVldContext();
   if(!vcp ){
      MSG("AlgFilterChannels", Msg::kFatal)
         << "Couldn't get a VldContext?!" << endl;
      return;
   }
*/
   // try to get a vld context from the digitlisthandle
   const VldContext* vcp = cdlh->GetVldContext();
   if(!vcp ){
      MSG("AlgFilterChannels", Msg::kFatal)
         << "Couldn't get a VldContext?!" << endl;
      return;
   }
   
 
   DbiResultPtr<ChannelsToFilter> rptr(*vcp);   

// 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
     MSG("AlgFilterChannels", Msg::kVerbose)
        << "Considering: "<<cdh->GetChannelId()
        <<" --> "<<cdh->GetPlexSEIdAltL()<< endl;

     // very simple, just ask the result pointer if this is a 
     // channel in its list
     bool cutit=false;
     MSG("AlgFilterChannels", Msg::kVerbose)
        << "Filtering out "<<rptr.GetNumRows()<<" channels."<<endl;

     for(unsigned int i=0; i<rptr.GetNumRows(); i++){
        const ChannelsToFilter* ctof = rptr.GetRow(i);
        MSG("AlgFilterChannels", Msg::kVerbose)
           <<"Is it this one? :  "<<(*ctof)<<endl;
           if( ctof->IsSameChannel(cdh) ) {
              cutit=true;
              MSG("AlgFilterChannels", Msg::kDebug)
                 <<"Cut this "<< (*ctof) <<endl;
              break;
           }
     }
     if(cutit) continue; // don't add this digit to the new list!

// Effect copy of wanted digits from orig CandDigitList to daughter list
     CandDigitHandle cddh(*cdh);
     ch.AddDaughterLink(cddh, kFALSE);           // Don't check for dups
   }
}

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

---