AlgMergeEvent.cxx

Go to the documentation of this file.

#include <cassert>
#include <cmath>

#include "Algorithm/AlgConfig.h"
#include "Algorithm/AlgFactory.h"
#include "Algorithm/AlgHandle.h"
#include "CandData/CandRecord.h"
#include "MuonRemoval/AlgMergeEvent.h"

#include "CandDigit/CandDigitHandle.h"
#include "CandDigit/CandDigit.h"
#include "CandDigit/CandDigitListHandle.h"

#include "RecoBase/CandTrackListHandle.h"
#include "RecoBase/CandTrackHandle.h"
#include "RecoBase/CandEventListHandle.h"
#include "RecoBase/CandEventHandle.h"
#include "RecoBase/CandStripHandle.h"

#include "Candidate/CandContext.h"
#include "MessageService/MsgService.h"
#include "Validity/VldContext.h"
#include "Validity/VldTimeStamp.h"

#include "RawDigitInfo.h"
#include "TClonesArray.h"

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

#include "PmtMap.h"

CVSID("$Id: AlgMergeEvent.cxx,v 1.3 2006/08/31 23:52:30 cbs Exp $");
 

ClassImp(AlgMergeEvent)

AlgMergeEvent::AlgMergeEvent()
{
}


AlgMergeEvent::~AlgMergeEvent()
{
}


void AlgMergeEvent::RunAlg(AlgConfig &ac, CandHandle &ch, CandContext &cx)
{

  //set up the sub-algs and this alg's config  
  AlgFactory &algfactory = AlgFactory::GetInstance();
  const char *tmpcs = 0;
  const char *conv_algorithm = 0;
  const char *conv_algconfig = 0;
  
  if (ac.Get("RawDigitInfoConvAlgorithm", tmpcs)) conv_algorithm = tmpcs;
  if (ac.Get("RawDigitInfoConvConfig", tmpcs)) conv_algconfig = tmpcs;

  AlgHandle conv_alg = algfactory.GetAlgHandle(conv_algorithm, conv_algconfig);
  CandContext cxx(this, cx.GetMom());
  cxx.SetCandRecord(cx.GetCandRecord());
  const VldContext *vldc = cx.GetCandRecord()->GetVldContext();

  MSG("EvtMrg",Msg::kDebug) << " AlgMergeEvent::RunAlg() " <<endl;
  MSG("EvtMrg",Msg::kDebug) << "   SubAlg         :  " <<conv_algorithm<<endl;
  MSG("EvtMrg",Msg::kDebug) << "   SubAlg Config  :  " <<conv_algconfig<<endl;

  //Get pointers to all the bits we need
  const TObjArray *input = (TObjArray *)(cx.GetDataIn());
  assert(input);
  TClonesArray* input_mc = dynamic_cast<TClonesArray*>(input->At(0));
  CandDigitListHandle* input_digits = dynamic_cast<CandDigitListHandle*>(input->At(1));
  RawRecord* rawrecord = dynamic_cast<RawRecord*>(input->At(2));
  assert(rawrecord);
  RawDigitDataBlock *input_datablock =  NULL;
  if(rawrecord){
    TIter rdbit = rawrecord->GetRawBlockIter();
    TObject *tob;
    while ((tob = rdbit())) {
      if (tob->InheritsFrom("RawDigitDataBlock")) {
        input_datablock = (RawDigitDataBlock *) tob;
      }
    }
  }

  //Build an array of raw data digits that are going to be used on the event
  vector<const RawDigit*> selected_data_digits;
  vector<int> selected_data_digits_index;
  TIter digititer(input_digits->GetDaughterIterator());
  while (CandDigitHandle *digit = dynamic_cast<CandDigitHandle*>(digititer())) {
    const RawDigit* raw_digit = input_datablock->At(digit->GetRawDigitIndex());
    if(raw_digit!=NULL){
      MSG("EvtMrg",Msg::kVerbose) << "Index TDC ADC Plane: "
                                  << digit->GetRawDigitIndex() << " " 
                                  << raw_digit->GetTDC() << " " 
                                  << raw_digit->GetADC() << " " 
                                  << digit->GetPlexSEIdAltL().GetPlane() << endl;
      selected_data_digits.push_back(raw_digit);
      selected_data_digits_index.push_back(digit->GetRawDigitIndex());
    }
    else {
      MSG("EvtMrg",Msg::kError) << "   Failed to find raw digit for reco'd digit " <<endl;
      return;
    }
  }

  //Build the PMT maps
  PmtMap *data_map = new PmtMap();
  data_map->BuildMap(selected_data_digits, selected_data_digits_index,vldc);
  PmtMap* mc_map = new PmtMap();
  mc_map->BuildMap(input_mc);
  
  //Make the merged Detector map
  PmtMap* merged_map = new PmtMap();
  merged_map->Merge(*data_map, *mc_map);

  /*
  //Do some debug print out TO BE REMOVED
  MSG("EvtMrg", Msg::kDebug)  << "***********************************"<<endl;
  MSG("EvtMrg", Msg::kDebug)  << "DATA: " << endl;
  data_map->PrintSummary();
  MSG("EvtMrg", Msg::kDebug)  << "***********************************"<<endl;
  MSG("EvtMrg", Msg::kDebug)  << "MC: " << endl;
  mc_map->PrintSummary();
  MSG("EvtMrg", Msg::kDebug)  << "***********************************"<<endl;
  MSG("EvtMrg", Msg::kDebug)  << "MERGE: " << endl;
  merged_map->PrintSummary();
  MSG("EvtMrg", Msg::kDebug)  << "***********************************"<<endl;
  */

  //delete the work space
  delete data_map; data_map = NULL;
  delete mc_map; mc_map = NULL;
  
  //Iterate over the whole chip set and make the new candidates
  const unsigned int nchannel = merged_map->channel.size();
  for(unsigned int ichannel=0; ichannel<nchannel; ichannel++){
    if(merged_map->channel[ichannel].valid!=0){
      cxx.SetDataIn(&(merged_map->channel[ichannel]));
      CandDigitHandle cdh = CandDigit::MakeCandidate(conv_alg, cxx);
      ch.AddDaughterLink(cdh, kFALSE);    
    }
  }

  //delete the detector map
  delete merged_map; merged_map = NULL;
  
  //global list stuff
  CandDigitListHandle &cdlh =  dynamic_cast<CandDigitListHandle &>(ch);
  cdlh.SetAbsTime(input_digits->GetAbsTime());                   
  cdlh.SetIsSparse(input_digits->GetIsSparse());
}

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

---