StripCands.cxx

Go to the documentation of this file.

#include <cassert>

#include "MuonRemoval/StripCands.h"
#include "MessageService/MsgService.h"
#include "MinosObjectMap/MomNavigator.h"
#include "JobControl/JobCModuleRegistry.h" // For JOBMODULE macro

#include <CandData/CandRecord.h>
#include <CandData/CandHeader.h>
#include <CandDigit/CandDigitListHandle.h>
#include <RecoBase/CandTrackListHandle.h>
#include <RecoBase/CandFitTrackHandle.h>
#include <RecoBase/CandTrackHandle.h>
#include <RecoBase/CandShowerHandle.h>
#include <RecoBase/CandShowerListHandle.h>
#include <RecoBase/CandStripHandle.h>
#include <RecoBase/CandStripListHandle.h>
#include <RecoBase/CandEventListHandle.h>
#include <RecoBase/CandEventHandle.h>
#include "StandardNtuple/NtpStRecord.h"

#include "RawData/RawRecord.h"
#include "RawData/RawDigit.h"
#include "RawData/RawDigitDataBlock.h"

JOBMODULE(StripCands, "StripCands",
          "strips cands");
CVSID("$Id: StripCands.cxx,v 1.5 2007/02/13 01:11:18 cbs Exp $");
//......................................................................

StripCands::StripCands() : 
  fKeepDigitList("")
{
}
//......................................................................

StripCands::~StripCands()
{
}

//......................................................................

void StripCands::BeginJob()
{
}

//......................................................................

void StripCands::EndJob()
{
}

//......................................................................

JobCResult StripCands::Ana(const MomNavigator* mom)
{
 CandRecord* candrec = dynamic_cast<CandRecord*>(mom->GetFragment("CandRecord"));
  if(!candrec){
    MSG("StripCands",Msg::kDebug) << " Unable to find CandRecord in mom " << endl;
    return JobCResult::kFailed;
  }

  const CandHeader* header  = dynamic_cast<const CandHeader*>(candrec->GetHeader());

  TIter iter(candrec->GetCandHandleIter());
  TObject* tobj = NULL;
  MSG("StripCands",Msg::kInfo)  << "************* New Mom **************" <<endl;
  if(header){
    MSG("StripCands",Msg::kInfo)  << "  Run:  " << header->GetRun()
                               <<  "    " << header->GetSnarl()<<endl;
  }

  RawRecord *rr = dynamic_cast<RawRecord *>(mom->GetFragment("RawRecord"));
  if (rr == 0) {
    MSG("StripCands", Msg::kWarning) << "No RawRecord in MOM." << endl;
    return JobCResult::kFailed;
  }
  const RawDigitDataBlock *rddb=dynamic_cast<const RawDigitDataBlock *>
    (rr->FindRawBlock("RawDigitDataBlock"));

  while((tobj = iter()) ){
    const CandHandle* cand = dynamic_cast<const CandHandle*>(tobj);
    assert(cand);
    MSG("StripCands",Msg::kInfo)  << cand->ClassName()<<":"
                               <<((cand->GetName())?cand->GetName():" ? ")
                               << "  N "  << cand->GetNDaughters() << endl;
    if(strcmp("CandDeMuxDigitListHandle", tobj->ClassName()) !=0 &&
       strcmp("CandDigitListHandle", tobj->ClassName()) !=0 &&
       strcmp("CandStripListHandle", tobj->ClassName()) !=0){
      TIter iter2(cand->GetDaughterIterator());
      TObject* tobj2 = 0;
      while((tobj2 = iter2())){
        const CandHandle* cand2= dynamic_cast<const CandHandle*>(tobj2);
        assert(cand2);
        MSG("StripCands",Msg::kInfo)  << "  " << cand2->ClassName()<<":"
                                   <<((cand2->GetName()!=NULL)?cand2->GetName():" ")<<endl;
      }
    }
    else if(MsgService::Instance()->IsActive("StripCands",Msg::kVerbose) && 
            strcmp("CandDigitListHandle", tobj->ClassName())==0){      
      TIter iter2(cand->GetDaughterIterator());
      TObject* tobj2 = 0;
      while((tobj2 = iter2())){
        const CandDigitHandle* cand2 = dynamic_cast<const CandDigitHandle*>(tobj2);
        const RawDigit* raw_digit = rddb->At(cand2->GetRawDigitIndex());
        MSG("StripCands",Msg::kVerbose) << cand2->GetRawDigitIndex() << " " 
                                        << raw_digit->GetADC() << " " 
                                        << raw_digit->GetTDC() << endl;
      }
    }
  }
  return JobCResult::kPassed; // kNoDecision, kFailed, etc.
}

//......................................................................

JobCResult StripCands::Reco(MomNavigator* mom)
{
  
 CandRecord* candrec = dynamic_cast<CandRecord*>(mom->GetFragment("CandRecord"));
  if(!candrec){
    MSG("StripCands",Msg::kDebug)  << " Unable to find CandRecord in mom " << endl;
    return JobCResult::kFailed;
  }

  const CandHeader* header  = dynamic_cast<const CandHeader*>(candrec->GetHeader());
  if(header){
    MSG("StripCands",Msg::kDebug)  << "Run:  " << header->GetRun()<<endl;
  }

  TIter iter(candrec->GetCandHandleIter());
  TObject* tobj = NULL;
  MSG("StripCands",Msg::kDebug)  << " New Mom : " <<endl;
  while( (tobj = iter()) ){
    CandHandle* cand = dynamic_cast<CandHandle*>(tobj);
    if( strcmp("candrmmulist", tobj->GetName()) ==0 ) { 
      MSG("StripCands",Msg::kDebug) << " Found and keeping: " 
                                 << tobj->ClassName() 
                                 <<":"<<tobj->GetName()<<endl;
      continue;
    }
    else if( strcmp(fKeepDigitList.c_str(), tobj->GetName()) !=0 ) {
      MSG("StripCands",Msg::kDebug) << " REMOVING: " <<  tobj->ClassName() 
                                 <<":"<<tobj->GetName()<<endl;
      candrec->RemoveCandHandle(cand);
    }
    else{
      MSG("StripCands",Msg::kDebug) <<" KEEPING: " <<  tobj->ClassName()
                                 <<":"<<tobj->GetName()<<endl;
      cand->SetName("canddigitlist");
      MSG("StripCands",Msg::kDebug) << "   Changed name " 
                                 << cand->GetName() <<endl; 
    }
  }

  /* NOT DOING THIS ANYMORE - CAN ALL BE CONFIGURED VIA MACRO
  //rename any Ntp objects:
  NtpStRecord *ntpst = dynamic_cast<NtpStRecord*> (mom->GetFragment("NtpStRecord"));
  if(ntpst) {
    MSG("StripCands",Msg::kDebug) << " RENAMING: " <<  ntpst->ClassName() 
                                  << ":" << ntpst->GetName() << endl;
    ntpst->SetName("NtpStRecordOld");
  }
  */
  
  return JobCResult::kPassed; // kNoDecision, kFailed, etc.
}

//......................................................................

const Registry& StripCands::DefaultConfig() const
{
  static Registry r; // Default configuration for module

  // Set name of config
  std::string name = this->GetName();
  name += ".config.default";
  r.SetName(name.c_str());

  // Set values in configuration
  r.UnLockValues();
  r.Set("keepdigitlist","mergedigitlist");
  r.LockValues();

  return r;
}

//......................................................................

void StripCands::Config(const Registry& r)
{
  fKeepDigitList = r.GetCharString("keepdigitlist");
}

---