AlgShowerEMList.cxx

Go to the documentation of this file.

// $Id: AlgShowerEMList.cxx,v 1.4 2004/09/28 23:47:27 cbs Exp $
//
// AlgShowerEMList.cxx
//
// AlgShowerEMList is a concrete ShowerEMList Algorithm class.

#include <cassert>

#include "Algorithm/AlgFactory.h"
#include "Algorithm/AlgHandle.h"
#include "Algorithm/AlgConfig.h"
#include "Candidate/CandContext.h"
#include "CandShowerEM/AlgShowerEMList.h"
#include "CandShowerEM/CandShowerEM.h"
#include "CandShowerEM/CandShowerEM.h"
#include "CandShowerEM/CandShowerEMHandle.h"
#include "CandShowerEM/CandShowerEMList.h"
#include "CandShowerEM/CandShowerEMListHandle.h"
#include "Conventions/PlaneView.h"
#include "MessageService/MsgService.h"
#include "MinosObjectMap/MomNavigator.h"
#include "Navigation/NavKey.h"
#include "Navigation/NavSet.h"
#include "RecoBase/CandClusterHandle.h"
#include "RecoBase/CandClusterListHandle.h"
#include "RecoBase/CandShowerHandle.h"
#include "RecoBase/CandSliceHandle.h"
#include "RecoBase/CandSliceListHandle.h"
#include "UgliGeometry/UgliGeomHandle.h"
#include "Validity/VldContext.h"

ClassImp(AlgShowerEMList)

//______________________________________________________________________
CVSID("$Id: AlgShowerEMList.cxx,v 1.4 2004/09/28 23:47:27 cbs Exp $");

//______________________________________________________________________
AlgShowerEMList::AlgShowerEMList()
{
}

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

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

  assert(cx.GetDataIn());
  if (!(cx.GetDataIn()->InheritsFrom("TObjArray"))) {
    return;
  }

  const CandSliceListHandle *slicelist = 0;
  const CandClusterListHandle *clusterlist = 0;
  const TObjArray *cxin = dynamic_cast<const TObjArray *>(cx.GetDataIn());
  for (Int_t i=0; i<=cxin->GetLast(); i++) {
    TObject *tobj = cxin->At(i);
    if (tobj->InheritsFrom("CandSliceListHandle")) {
      slicelist = dynamic_cast<CandSliceListHandle*>(tobj);
    }
    if (tobj->InheritsFrom("CandClusterListHandle")) {
      clusterlist = dynamic_cast<CandClusterListHandle*>(tobj);
    }
  }
  if (!slicelist || !clusterlist) {
    MSG("error",Msg::kError) <<
      "CandSliceListHandle or CandClusterListHandle missing\n";
  }

  // Create Candcontext
  CandContext cxx(this,cx.GetMom());

  //get config for AlgShowerEM
  const char *charShowerEMAlgConfig = 0;
  ac.Get("ShowerEMAlgConfig",charShowerEMAlgConfig);

  // Get singleton instance of AlgFactory
  AlgFactory &af = AlgFactory::GetInstance();  
  AlgHandle ah = af.GetAlgHandle("AlgShowerEM",charShowerEMAlgConfig);

  const CandRecord *candrec = cx.GetCandRecord();
  assert(candrec);
  const VldContext *vldcptr = candrec->GetVldContext();
  assert(vldcptr);
  VldContext vldc = *vldcptr;

  UgliGeomHandle ugh(vldc);
  
  CandSliceHandleItr sliceItr(slicelist->GetDaughterIterator());
  while (CandSliceHandle *slice = sliceItr()) {
    CandClusterHandleItr clusterItr(clusterlist->GetDaughterIterator());
    TObjArray newshower;
    while (CandClusterHandle *cluster = clusterItr()){
      if (*cluster->GetCandSlice()==*slice) {
        newshower.Add(cluster);
      }
    }
    cxx.SetDataIn(&newshower);
    MSG("ShowerEM",Msg::kDebug) << "forming shower\n";
    CandShowerEMHandle showerhandle = CandShowerEM::MakeCandidate(ah,cxx);
    if(showerhandle.GetNStrip()!=0) {
      showerhandle.SetCandSlice(slice);
      ch.AddDaughterLink(showerhandle);    
    }
    else MSG("ShowerEM",Msg::kWarning) << "Removing empty shower\n";
  }
}

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

---