AlgCluster3DList.cxx

Go to the documentation of this file.

//
// AlgCluster3DList.cxx (based on R.Lee's  AlgCluster3DList)
//
// AlgCluster3DList is a concrete Cluster3DList Algorithm class.
//
// Author:  K.Grzelak1@physics.ox.ac.uk
//

#include <cassert>

#include "Algorithm/AlgFactory.h"
#include "Algorithm/AlgHandle.h"
#include "Algorithm/AlgConfig.h"
#include "Candidate/CandContext.h"
#include "Candidate/CandHandle.h"
#include "Cluster3D/AlgCluster3DList.h"
#include "Cluster3D/CandCluster3D.h"
#include "Cluster3D/CandCluster3D.h"
#include "Cluster3D/CandCluster3DHandle.h"
#include "Cluster3D/CandCluster3DList.h"
#include "Cluster3D/CandCluster3DListHandle.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/CandSliceHandle.h"
#include "RecoBase/CandSliceListHandle.h"
#include "UgliGeometry/UgliGeomHandle.h"
#include "Validity/VldContext.h"
#include "Cluster3D/StripKG.h"
#include "Cluster3D/StripListKG.h"
#include "Cluster3D/CellKG.h"
#include "Cluster3D/CellListKG.h"
#include "Cluster3D/ClusterKG.h"
#include "Cluster3D/ClusterListKG.h"

ClassImp(AlgCluster3DList)

CVSID("$Id:");

//______________________________________________________________________
void AlgCluster3DList::RunAlg(AlgConfig &ac, CandHandle &ch, CandContext &cx)
//
//  Purpose:  From here, the methods which are the heart of the algorithm,
//            are called 
//
//  Arguments: 
//              
//  Return:  n/a
//
//  Contact:   K.Grzelak1@physics.ox.ac.uk
//
//
//  Program Notes:-
//  =============
//
{
   MSG("Cluster3D", Msg::kDebug)
                           << "Starting AlgCluster3DList::RunAlg()" << endl;

  assert(cx.GetDataIn());

// Check for CandSliceListHandle input
  if (!(cx.GetDataIn()->InheritsFrom("CandSliceListHandle"))) {
    MSG("Cluster3D", Msg::kError)
        << "AlgCluster3DList::RunAlg() No CandSliceListHandle input" << endl;
    return;
  }

  const CandSliceListHandle *cslh =
        dynamic_cast<const CandSliceListHandle*>(cx.GetDataIn());

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

  // Algorithm parameters

// Get singleton instance of AlgFactory
  AlgFactory &af = AlgFactory::GetInstance();

// Get an AlgHandle to AlgCluster3D with default AlgConfig
  AlgHandle ah = af.GetAlgHandle("AlgCluster3D","default");

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

  UgliGeomHandle ugh0(vldc);
  UgliGeomHandle *ugh = &ugh0;

  Int_t islice=0;

  CandSliceHandleItr sliceItr(cslh->GetDaughterIterator());
  // Different way to iterate over the structure
  //while (CandSliceHandle *slice = sliceItr()) {
  while ( sliceItr.IsValid() ) {

      MSG("Cluster3D", Msg::kDebug)
        << "Slice number = " << islice++ << endl;

      CandSliceHandle *slice = sliceItr.Ptr();    

      StripListKG stripList;
      stripList.CreateStrips(ac,slice,ugh);

      vector<StripKG> stps;   
      stripList.GetStrips(stps);

      MSG("Cluster3D", Msg::kDebug)
        << "Number of strips =  " << stps.size() << endl;
 
      CellListKG cellList;

      cellList.CreateCells(ac,stripList);        

      vector<CellKG> cells;   
      cellList.GetCells(cells);

      MSG("Cluster3D", Msg::kDebug)
        << "Number of cells = " << cells.size() << endl;

      ClusterListKG clusterList;

      clusterList.CreateClusters(ac,stripList,cellList);

      vector<ClusterKG> clsts;   
      clusterList.GetClusters(clsts);

      MSG("Cluster3D", Msg::kDebug)
        << "Number of clusters = " << clsts.size() << endl;

      for(UInt_t i=0;i<clsts.size();++i){
        TObjArray striparray;
        MSG("Cluster3D", Msg::kDebug)
          << "Number of strips in " << i << " cluster= " 
          << clsts[i].stpindexes.size() << endl;
        for(UInt_t j=0;j<clsts[i].stpindexes.size();++j){
          CandStripHandle *strip = 
              stripList.GetStripHandle(clsts[i].stpindexes[j]);
          striparray.Add(strip);
        }
        cxx.SetDataIn(&striparray);
        CandCluster3DHandle clusterhandle = 
                 CandCluster3D::MakeCandidate(ah,cxx);
        ch.AddDaughterLink(clusterhandle);         
      }

      // something that inherits from TObject and is used 
      // to create clusters - in this case we should use strips
      // as an argument to SetDataIn

      sliceItr.Next();
  }

        /*
// don't cluster in muon spectrometer
        if (strip->GetPlaneView()==view &&
            (vldc.GetDetector()!=DetectorType::kNear ||
            strip->GetPlane()<=121)) {
        }
        */

}

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

---