CandClusterHandle.cxx

Go to the documentation of this file.

// $Id: CandClusterHandle.cxx,v 1.20 2005/11/01 20:12:24 gmieg Exp $
//
// CandClusterHandle.cxx
//
// CandClusterHandle is the specialized access handle to CandCluster.
//
// Each concrete CandHandle must define a DupHandle function.
//
// Author:  R. Lee

#include <cassert>
#include <iostream>

#include "MessageService/MsgService.h"
#include "Navigation/NavKey.h"
#include "Navigation/NavSet.h"
#include "Plex/PlexStripEndId.h"
#include "RecoBase/CandCluster.h"
#include "RecoBase/CandClusterHandle.h"
#include "RecoBase/CandSliceHandle.h"
#include "RecoBase/CandStripHandle.h"

ClassImp(CandClusterHandle)

//______________________________________________________________________
CVSID("$Id: CandClusterHandle.cxx,v 1.20 2005/11/01 20:12:24 gmieg Exp $");

//______________________________________________________________________
CandClusterHandle::CandClusterHandle()
{
}

//______________________________________________________________________
CandClusterHandle::CandClusterHandle(const CandClusterHandle &cdh) :
  CandHandle(cdh)
{
}

//______________________________________________________________________
CandClusterHandle::CandClusterHandle(CandCluster *cd) :
  CandHandle(cd)
{
}

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

//______________________________________________________________________
CandClusterHandle *CandClusterHandle::DupHandle() const
{
  return (new CandClusterHandle(*this));
}

//______________________________________________________________________
void CandClusterHandle::Trace(const char *c) const
{
  MSG("Cand", Msg::kDebug)
  << "**********Begin CandClusterHandle::Trace(\"" << c << "\")" << endl
         << "Information from CandClusterHandle's CandHandle: " << endl;
  CandHandle::Trace(c);
  MSG("Cand", Msg::kDebug)
   << "**********End CandClusterHandle::Trace(\"" << c << "\")" << endl;
}

//______________________________________________________________________
void CandClusterHandle::SetCandSlice(const CandSliceHandle *slice)
{
  if (slice) {                 // fCandSlice is now an owned CandHandle*
    CandSliceHandle *ch = slice->DupHandle();
    delete (dynamic_cast<CandCluster *>
                                      (GetOwnedCandBase()))->fCandSlice;
    dynamic_cast<CandCluster *>(GetOwnedCandBase())->fCandSlice = ch;
  }
}

const CandSliceHandle *CandClusterHandle::GetCandSlice() const
{
  return dynamic_cast<const CandCluster *>(GetCandBase())->fCandSlice;
}

//______________________________________________________________________

void CandClusterHandle::IsTrackLike(Bool_t bvar)
{
  dynamic_cast<CandCluster *>(GetOwnedCandBase())->fTrackLike = bvar;
}

Bool_t CandClusterHandle::IsTrackLike() const
{
  return dynamic_cast<const CandCluster *>(GetCandBase())->fTrackLike;
}

//______________________________________________________________________

void CandClusterHandle::IsShowerLike(Bool_t bvar)
{
  dynamic_cast<CandCluster *>(GetOwnedCandBase())->fShowerLike = bvar;

}

Bool_t CandClusterHandle::IsShowerLike() const
{
  return dynamic_cast<const CandCluster *>(GetCandBase())->fShowerLike;
}

//______________________________________________________________________

PlaneView::PlaneView_t CandClusterHandle::GetPlaneView() const 
{
  CandStripHandleItr cshItr(GetDaughterIterator());
  if (CandStripHandle *csh = cshItr()) {
    return csh->GetPlaneView();
  }
  else {
    return PlaneView::kUnknown;
  }
}

//______________________________________________________________________

Int_t CandClusterHandle::GetNStrip() const 
{
  return GetNDaughters();
}

Int_t CandClusterHandle::GetNPlane() const
{
  Int_t begplane = GetBegPlane();
  if (begplane==0) {
    return 0;
  }
  Int_t endplane = GetEndPlane();
  if (begplane==endplane) {
    return 1;
  }
  Int_t ncount = 2;
  CandStripHandleItr stripItr(GetDaughterIterator());
  CandStripHandleKeyFunc *stripKf = stripItr.CreateKeyFunc();
  stripKf->SetFun(CandStripHandle::KeyFromPlane);
  stripItr.GetSet()->AdoptSortKeyFunc(stripKf);
  stripKf = 0;
  while (stripItr.IsValid() && stripItr.Ptr()->GetPlane()<=begplane) {
    stripItr.Next();
  }
  Int_t oldplane = stripItr.Ptr()->GetPlane();
  while (CandStripHandle *strip =
                           dynamic_cast<CandStripHandle*>(stripItr())) {
    if (strip->GetPlaneView()!=PlaneView::kA &&
        strip->GetPlaneView()!=PlaneView::kB) {
      if (strip->GetPlane()!=oldplane) {
        ncount++;
      }
      oldplane = strip->GetPlane();
    }
  }
  return ncount;
}

//______________________________________________________________________

Int_t CandClusterHandle::GetBegPlane() const 
{
  TIter cshItr(GetDaughterIterator());
  Int_t plane = 0;
  Bool_t first = 1;
  while (CandStripHandle *csh =
                             dynamic_cast<CandStripHandle*>(cshItr())) {
    if (first || csh->GetPlane()<plane) {
      plane = csh->GetPlane();
      first = 0;
    }
  }
  return plane;
}

Int_t CandClusterHandle::GetEndPlane() const 
{
  TIter cshItr(GetDaughterIterator());
  Int_t plane = 0;
  Bool_t first = 1;
  while (CandStripHandle *csh =
                             dynamic_cast<CandStripHandle*>(cshItr())) {
    if (first || csh->GetPlane()>plane) {
      plane = csh->GetPlane();
      first = 0;
    }
  }
  return plane;
}

//______________________________________________________________________

Double_t CandClusterHandle::GetBegTime() const 
{
  TIter cshItr(GetDaughterIterator());
  Double_t time = 0.;
  Bool_t first = 1;
  while (CandStripHandle *csh =
                             dynamic_cast<CandStripHandle*>(cshItr())) {
    if (first || csh->GetBegTime()<time) {
      time = csh->GetBegTime();
      first = 0;
    }
  }
  return time;
}

Double_t CandClusterHandle::GetEndTime() const 
{
  TIter cshItr(GetDaughterIterator());
  Double_t time = 0.;
  Bool_t first = 1;
  while (CandStripHandle *csh =
                             dynamic_cast<CandStripHandle*>(cshItr())) {
    if (first || csh->GetEndTime()>time) {
      time = csh->GetEndTime();
      first = 0;
    }
  }
  return time;
}

//______________________________________________________________________

Double_t CandClusterHandle::GetCharge() const 
{
  TIter cshItr(GetDaughterIterator());
  Double_t charge = 0.;
  while (CandStripHandle *csh =
                             dynamic_cast<CandStripHandle*>(cshItr())) {
    charge += csh->GetCharge();
  }
  return charge;
}

//______________________________________________________________________

Double_t CandClusterHandle::GetTPos() const 
{
  TIter cshItr(GetDaughterIterator());
  Double_t charge = 0.;
  Double_t tpos = 0.;
  while (CandStripHandle *csh =
                             dynamic_cast<CandStripHandle*>(cshItr())) {
    charge += csh->GetCharge();
    tpos += csh->GetCharge()*csh->GetTPos();
  }
  tpos /= charge;
  return tpos;
}


//______________________________________________________________________
// Navigation Helpers
//______________________________________________________________________

NavKey CandClusterHandle::KeyFromObject(const CandClusterHandle *cluster)
{
  if (cluster->IsShowerLike()) {
    return 1;
  }
  if (cluster->IsTrackLike()) {
    return 2;
  }
  return 0;
}


NavKey CandClusterHandle::KeyFromView(const CandClusterHandle *cluster)
{
  if (const CandStripHandle *strip =
       dynamic_cast<const CandStripHandle *>(cluster->FindDaughter())) {
    return static_cast<Int_t>(const_cast<CandStripHandle *>
      (strip)->GetPlaneView());
  }
  return static_cast<Int_t>(PlaneView::kUnknown);
}

NavKey CandClusterHandle::KeyFromSlice(const CandClusterHandle *cluster)
{
  if (cluster->GetCandSlice()) {
    return static_cast<Int_t>(cluster->GetCandSlice()->GetUidInt());
  } 
  return 0;

}


XXXITRIMP(CandClusterHandle)

---