CandSliceHandle.cxx

Go to the documentation of this file.

// $Id: CandSliceHandle.cxx,v 1.12 2004/08/19 16:15:01 musser Exp $
//
// CandSliceHandle.cxx
//
// Begin_Html<img src="../../pedestrians.gif" align=center>
// <a href="../source_warning.html">Warning for beginners</a>.<br> 
//
// CandSliceHandle is the specialized access handle to CandSlice.
//
// Each concrete CandHandle must define a DupHandle function.
//
// Author:  R. Lee
//
// Also see <a href="../../root_crib/index.html">The ROOT Crib</a> and 
// <a href="../CandSlice.html"> CandSlice Classes</a> (part of
// <a href="../index.html">The MINOS Class User Guide</a>)End_Html

#include <cassert>
#include <iostream>
#include <map>

#include "RecoBase/CandSliceHandle.h"
#include "RecoBase/CandSlice.h"
#include "MessageService/MsgService.h"

ClassImp(CandSliceHandle)

//______________________________________________________________________
CVSID("$Id: CandSliceHandle.cxx,v 1.12 2004/08/19 16:15:01 musser Exp $");

//______________________________________________________________________
CandSliceHandle::CandSliceHandle()
{
}

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

//______________________________________________________________________
CandSliceHandle::CandSliceHandle(CandSlice *cd) :
CandHandle(cd)
{
}

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

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


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

Int_t CandSliceHandle::GetNStrip(PlaneView::PlaneView_t planeview_t) const
{
  if (planeview_t==PlaneView::kUnknown) {
    return GetNDaughters();
  }
  Int_t n=0;
  TIter stripItr(GetDaughterIterator());
  CandStripHandle *strip;
  while ((strip = dynamic_cast<CandStripHandle*>(stripItr()))) {
    PlaneView::PlaneView_t planeview = strip->GetPlaneView();
    if (planeview == planeview_t) {
      n++;
    }
  }
  return n;
}

Int_t CandSliceHandle::GetNDigit(StripEnd::StripEnd_t stripend_t) const
{
  Int_t n=0;
  TIter stripItr(GetDaughterIterator());
  CandStripHandle *strip;
  while ((strip = dynamic_cast<CandStripHandle*>(stripItr()))) {
    n += strip->GetNDigit(stripend_t);
  }
  return n;
}


Int_t CandSliceHandle::GetNPlane(PlaneView::PlaneView_t planeview_t) const
{
  TIter stripItr(GetDaughterIterator());
  map<Int_t,Bool_t> occupancy;
  CandStripHandle *strip;
  Int_t minplane = -1;
  Int_t maxplane = -1;
  while ((strip = dynamic_cast<CandStripHandle*>(stripItr()))) {
    PlaneView::PlaneView_t planeview = strip->GetPlaneView();
    if (planeview!=PlaneView::kA && planeview!=PlaneView::kB &&
        (planeview_t==PlaneView::kUnknown || planeview==planeview_t)) {
      occupancy[strip->GetPlane()] = kTRUE;
      if (minplane<0 || strip->GetPlane()<minplane) minplane = strip->GetPlane();
      if (maxplane<0 || strip->GetPlane()>maxplane) maxplane = strip->GetPlane();
    }
  }
  Int_t nplane = 0;
  map<Int_t,Bool_t>::iterator iter;
  for (Int_t iplane=minplane; iplane<=maxplane; iplane++) {
    iter = occupancy.find(iplane);
    if (iter!=occupancy.end()) {
      nplane++;
    }
  }

  return nplane;
}


Int_t CandSliceHandle::GetBegPlane(PlaneView::PlaneView_t planeview_t) const
{
  TIter stripItr(GetDaughterIterator());
  CandStripHandle *strip;
  Bool_t first(1);
  Int_t plane=0;
  while ((strip = dynamic_cast<CandStripHandle*>(stripItr()))) {
    PlaneView::PlaneView_t planeview = strip->GetPlaneView();
    if (planeview!=PlaneView::kA && planeview!=PlaneView::kB &&
        (planeview_t==PlaneView::kUnknown || planeview_t==planeview) &&
        (first || strip->GetPlane()<plane)) {
      first = 0;
      plane = strip->GetPlane();
    }
  }
  return plane;
}

Int_t CandSliceHandle::GetEndPlane(PlaneView::PlaneView_t planeview_t) const
{
  TIter stripItr(GetDaughterIterator());
  CandStripHandle *strip;
  Bool_t first(1);
  Int_t plane=0;
  while ((strip = dynamic_cast<CandStripHandle*>(stripItr()))) {
    PlaneView::PlaneView_t planeview = strip->GetPlaneView();
    if (planeview!=PlaneView::kA && planeview!=PlaneView::kB &&
        (planeview_t==PlaneView::kUnknown || planeview_t==planeview) &&
        (first || strip->GetPlane()>plane)) {
      first = 0;
      plane = strip->GetPlane();
    }
  }
  return plane;
}


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


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


Double_t CandSliceHandle::GetCharge(CalDigitType::CalDigitType_t caltype_t) const 
{
  TIter stripItr(GetDaughterIterator());
  Double_t charge = 0.;
  while (CandStripHandle *strip = dynamic_cast<CandStripHandle*>(stripItr())) {
    charge += strip->GetCharge(caltype_t);
  }
  return charge;
}

//______________________________________________________________________
Double_t CandSliceHandle::GetCorrTime() const
{
 Double_t corrtime=0.;
  Double_t totph=0.;
  TIter stripItr(GetDaughterIterator());
  while (CandStripHandle *striphandle = dynamic_cast<CandStripHandle*>(stripItr())) {
     corrtime += striphandle->GetCharge()*striphandle->GetCorrBegTime();
     totph += striphandle->GetCharge();    
  }
  if (totph>0.) corrtime /= totph;
  return corrtime;
}


//______________________________________________________________________
Double_t CandSliceHandle::GetCorrBegTime() const
{
 Double_t corrtime=0.;
  TIter stripItr(GetDaughterIterator());
  while (CandStripHandle *striphandle = dynamic_cast<CandStripHandle*>(stripItr())) {
    if (!corrtime || striphandle->GetCorrBegTime()<corrtime) {
      corrtime = striphandle->GetCorrBegTime();
     }
  }
  return corrtime;
}



//______________________________________________________________________
Double_t CandSliceHandle::GetCorrEndTime() const
{
 Double_t corrtime=0.;
 TIter stripItr(GetDaughterIterator());
 while (CandStripHandle *striphandle = dynamic_cast<CandStripHandle*>(stripItr())) {
   if (!corrtime || striphandle->GetCorrBegTime()>corrtime) {
     corrtime = striphandle->GetCorrBegTime();
   }
 }
 return corrtime;
}





XXXITRIMP(CandSliceHandle)

---