AlgFitShowerEMList.cxx

Go to the documentation of this file.

// $Id: AlgFitShowerEMList.cxx,v 1.2 2004/09/29 00:07:29 cbs Exp $
//
// AlgFitShowerEMList.cxx
//

#include <cassert>
#include <vector>
extern "C" {
#include <unistd.h>    // sysconf
#include <sys/times.h> // times()
}

#include "Algorithm/AlgFactory.h"
#include "Algorithm/AlgHandle.h"
#include "Algorithm/AlgConfig.h"
#include "Candidate/CandContext.h"
#include "CandFitShowerEM/AlgFitShowerEMList.h"
#include "CandFitShowerEM/CandFitShowerEM.h"
#include "CandFitShowerEM/CandFitShowerEMHandle.h"
#include "CandFitShowerEM/CandFitShowerEMList.h"
#include "CandFitShowerEM/CandFitShowerEMListHandle.h"
#include "CandShowerEM/CandShowerEMListHandle.h"
#include "Conventions/Mphysical.h"
#include "Conventions/Munits.h"
#include "Conventions/PlaneView.h"
#include "MessageService/MsgService.h"
#include "MinosObjectMap/MomNavigator.h"
#include "Navigation/NavKey.h"
#include "Navigation/NavSet.h"
#include "RecoBase/CandSliceHandle.h"
#include "RecoBase/CandSliceListHandle.h"
#include "RecoBase/LinearFit.h"
#include "UgliGeometry/UgliGeomHandle.h"
#include "Validity/VldContext.h"

ClassImp(AlgFitShowerEMList)


//______________________________________________________________________
CVSID("$Id: AlgFitShowerEMList.cxx,v 1.2 2004/09/29 00:07:29 cbs Exp $");

//______________________________________________________________________
AlgFitShowerEMList::AlgFitShowerEMList()
{

}

//______________________________________________________________________
AlgFitShowerEMList::~AlgFitShowerEMList()
{

}

//______________________________________________________________________
void AlgFitShowerEMList::RunAlg(AlgConfig &ac, CandHandle &ch,CandContext &cx)
{
  MSG("Alg", Msg::kDebug)
    << "Starting AlgFitShowerEMList::RunAlg()" << endl;
  
  assert(cx.GetDataIn());
  
  if (!(cx.GetDataIn()->InheritsFrom("TObjArray"))) {
    return;
  }

  const CandSliceListHandle *slicelist = 0;
  const CandShowerEMListHandle *showerlist = 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("CandShowerEMListHandle")) {
      showerlist = dynamic_cast<CandShowerEMListHandle*>(tobj);
    }
  }
  if (!slicelist || !showerlist) {
    MSG("error",Msg::kError) <<
      "CandSliceListHandle or CandShowerEMListHandle missing\n";
  }
  
  // Create Candcontext
  CandContext cxx(this,cx.GetMom());

  //get FitShowerEM config:
  const char *charFitShowerEMAlgConfig = 0;
  ac.Get("FitShowerEMAlgConfig",charFitShowerEMAlgConfig);

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

  TIter showerItr(showerlist->GetDaughterIterator());

  clock_t dummyt;
  struct tms t1;
  struct tms t2;
  static double ticksPerSecond = sysconf(_SC_CLK_TCK);
  while (CandShowerEMHandle *shower 
         = dynamic_cast<CandShowerEMHandle*>(showerItr())) {
    if(shower->GetShwStatus()!=1) continue;
    TObjArray cxin;
    cxin.Add(shower);
    cxx.SetDataIn(&cxin);
    dummyt = times(&t1);
    CandFitShowerEMHandle fitshower = CandFitShowerEM::MakeCandidate(ah,cxx);
    dummyt = times(&t2);
    Double_t CPUTime = Double_t(t2.tms_utime+t2.tms_stime
                                -t1.tms_utime-t1.tms_stime)/ticksPerSecond;
    fitshower.SetCPUTime(CPUTime);
    MSG("FitShowerEM",Msg::kDebug) << "CPUTime = "<<CPUTime<<endl;
    ch.AddDaughterLink(fitshower);
  }

}

---