CandStripHandle.cxx

Go to the documentation of this file.

// $Id: CandStripHandle.cxx,v 1.33 2007/11/08 13:50:05 musser Exp $
//
// CandStripHandle.cxx
//
// Begin_Html<img src="../../pedestrians.gif" align=center>
// <a href="../source_warning.html">Warning for beginners</a>.<br> 
//
// CandStripHandle is the specialized access handle to CandStrip.
//
// 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="../CandStrip.html"> CandStrip Classes</a> (part of
// <a href="../index.html">The MINOS Class User Guide</a>)End_Html

#include <cassert>
#include <iostream>

#include "CandDigit/CandDigitHandle.h"
#include "Conventions/CalDigitType.h"
#include "Conventions/StripEnd.h"
#include "MessageService/MsgService.h"
#include "Navigation/NavKey.h"
#include "Navigation/NavSet.h"
#include "Plex/PlexStripEndId.h"
#include "RecoBase/CandStripHandle.h"
#include "RecoBase/CandStrip.h"
#include "Validity/VldContext.h"
#include "Conventions/Mphysical.h"
#include "Conventions/Munits.h"
#include "UgliGeometry/UgliGeomHandle.h"
#include "UgliGeometry/UgliScintPlnHandle.h"
#include "Plex/PlexPlaneId.h"
ClassImp(CandStripHandle)

//______________________________________________________________________
CVSID("$Id: CandStripHandle.cxx,v 1.33 2007/11/08 13:50:05 musser Exp $");

//______________________________________________________________________
CandStripHandle::CandStripHandle()
{
}

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

//______________________________________________________________________
CandStripHandle::CandStripHandle(CandStrip *cd) :
  CandHandle(cd)
{
}

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

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


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


//______________________________________________________________________


Int_t CandStripHandle::GetPlane() const
{
  return GetStripEndId().GetPlane();
}

//______________________________________________________________________


Int_t CandStripHandle::GetStrip() const
{
  return  GetStripEndId().GetStrip();
}

//______________________________________________________________________

void CandStripHandle::SetZPos(Double_t dvar)
{
  dynamic_cast<CandStrip *>(GetOwnedCandBase())->fZPos = dvar;
}
Double_t CandStripHandle::GetZPos() const
{
  return   dynamic_cast<const CandStrip *>(GetCandBase())->fZPos;
}

//______________________________________________________________________

void CandStripHandle::SetTPos(Double_t dvar)
{
  dynamic_cast<CandStrip *>(GetOwnedCandBase())->fTPos = dvar;
}
Double_t CandStripHandle::GetTPos() const
{
  return dynamic_cast<const CandStrip *>(GetCandBase())->fTPos;
}
Double_t CandStripHandle::GetTPos(Double_t orthopos) const
{
  UgliGeomHandle ugh(*GetVldContext());
  PlexPlaneId plnid(GetVldContext()->GetDetector(),GetPlane(),false);
  UgliScintPlnHandle scintpln = ugh.GetScintPlnHandle(plnid); 
  UgliStripHandle ush = scintpln.GetStripHandle(GetStripEndId());
  return ush.GetTPos(orthopos);
}


//______________________________________________________________________
void CandStripHandle::SetCalTimeType(CalTimeType::CalTimeType_t caltimetype)
{
  dynamic_cast<CandStrip *>(GetOwnedCandBase())->fCalTimeType = caltimetype;
}

CalTimeType::CalTimeType_t CandStripHandle::GetCalTimeType() const
{
  return dynamic_cast<const CandStrip *>(GetCandBase())->fCalTimeType;
}
//______________________________________________________________________

Double_t CandStripHandle::GetTime(StripEnd::StripEnd_t stripend_t) const {
  Int_t i=0;
  Double_t time=0.;
  Double_t ph=0.;
  CalTimeType::CalTimeType_t caltimetype = GetCalTimeType();
  TIter cdhItr(GetDaughterIterator());
  while (CandDigitHandle *cdh = dynamic_cast<CandDigitHandle*>(cdhItr())) {
    Double_t digittime = cdh->GetSubtractedTime(caltimetype);
    if (stripend_t==StripEnd::kWhole ||
        cdh->GetPlexSEIdAltL().GetEnd()==stripend_t) {
      time += digittime*cdh->GetCharge();
      ph += cdh->GetCharge();
      i++;
    }
  }
  if (ph>0.) time /= ph;
  return time;
}

Double_t CandStripHandle::GetBegTime(StripEnd::StripEnd_t stripend_t) const {
  Int_t i=0;
  Double_t time=0.;
  CalTimeType::CalTimeType_t caltimetype = GetCalTimeType();
  TIter cdhItr(GetDaughterIterator());
  while (CandDigitHandle *cdh = dynamic_cast<CandDigitHandle*>(cdhItr())) {
    Double_t digittime = cdh->GetSubtractedTime(caltimetype);
    if ((!i || digittime<time) && (stripend_t==StripEnd::kWhole ||
         cdh->GetPlexSEIdAltL().GetEnd()==stripend_t)) {
      time = digittime;
      i++;
    }
  }
  return time;
}

Double_t CandStripHandle::GetCorrBegTime() const {
  return GetTime()-GetZPos()/Mphysical::c_light;
}

Double_t CandStripHandle::GetEndTime(StripEnd::StripEnd_t stripend_t) const {
  Int_t i=0;
  Double_t time=0.;
  CalTimeType::CalTimeType_t caltimetype = GetCalTimeType();
  TIter cdhItr(GetDaughterIterator());
  while (CandDigitHandle *cdh = dynamic_cast<CandDigitHandle*>(cdhItr())) {
     if ((!i || cdh->GetSubtractedTime(caltimetype)>time) && (stripend_t==StripEnd::kWhole ||
         cdh->GetPlexSEIdAltL().GetEnd()==stripend_t)) {
       time = cdh->GetSubtractedTime(caltimetype);
       i++;
     }
  }
  return time;
}

//______________________________________________________________________

Double_t CandStripHandle::GetCharge(StripEnd::StripEnd_t stripend_t, CalDigitType::CalDigitType_t caltype_t) const {
  Double_t charge=0.;
  TIter cdhItr(GetDaughterIterator());
  while (CandDigitHandle *cdh = dynamic_cast<CandDigitHandle*>(cdhItr())) {
     if (stripend_t==StripEnd::kWhole ||
         cdh->GetPlexSEIdAltL().GetEnd()==stripend_t) {
       charge += cdh->GetCharge(caltype_t);
     }
  }
  return charge;
}
Double_t CandStripHandle::GetCharge(CalDigitType::CalDigitType_t caltype_t, StripEnd::StripEnd_t stripend_t) const {
  return GetCharge(stripend_t,caltype_t);
}

Double_t CandStripHandle::GetCharge() const
{
  return GetCharge(StripEnd::kWhole,CalDigitType::kPE);
}

Double_t CandStripHandle::GetCharge(StripEnd::StripEnd_t stripend_t) const
{
  return GetCharge(stripend_t,CalDigitType::kPE);
}

Double_t CandStripHandle::GetCharge(CalDigitType::CalDigitType_t caltype_t) const
{
  return GetCharge(StripEnd::kWhole,caltype_t);
}

//______________________________________________________________________
Int_t CandStripHandle::GetNDigit(StripEnd::StripEnd_t stripend_t) const {
  if (stripend_t == StripEnd::kWhole) {
    return GetNDaughters();
  }
  Int_t n=0;
  TIter cdhItr(GetDaughterIterator());
  while (CandDigitHandle *cdh = dynamic_cast<CandDigitHandle*>(cdhItr())) {
    if (cdh->GetPlexSEIdAltL().GetEnd() == stripend_t) {
      n++;
    }
  }
  return n;
}

//______________________________________________________________________

UShort_t CandStripHandle::GetDemuxVetoFlag(StripEnd::StripEnd_t stripend_t) const
{
  UShort_t demuxveto(1);
  TIter digitItr(GetDaughterIterator());
  while (CandDigitHandle *digit = dynamic_cast<CandDigitHandle*>(digitItr())) {
    if (stripend_t==StripEnd::kWhole || digit->GetPlexSEIdAltL().GetBestSEId().GetEnd()==stripend_t) {
      if (digit->GetPlexSEIdAltL().GetDemuxVetoFlag()==0) demuxveto = 0;
    }
  }
  return demuxveto;
}

//______________________________________________________________________

PlaneView::PlaneView_t CandStripHandle::GetPlaneView() const 
{
  return GetStripEndId().GetPlaneView();
}

//______________________________________________________________________

void CandStripHandle::SetBestStripEndId(PlexStripEndId stripendid)
{
  CandStrip *candstrip = dynamic_cast<CandStrip *>(GetOwnedCandBase());
  candstrip->fBestStripEnd = stripendid;
  candstrip->fBestStripEnd.SetEnd(StripEnd::kWhole);
}

PlexStripEndId CandStripHandle::GetStripEndId(StripEnd::StripEnd_t stripend_t)
const {
  const CandStrip *candstrip = dynamic_cast<const CandStrip *>(GetCandBase());
  PlexStripEndId seid = candstrip->fBestStripEnd;
  seid.SetEnd(stripend_t);
  return seid;
}


//______________________________________________________________________
//  Navigation Helpers
//______________________________________________________________________
//______________________________________________________________________
NavKey CandStripHandle::KeyFromPlane(const CandStripHandle *strip)
{
  return const_cast<CandStripHandle *>(strip)->GetPlane();
}

//______________________________________________________________________
NavKey CandStripHandle::KeyFromPlaneStripTime(const CandStripHandle *strip)
{

   Int_t iplane = strip->GetPlane();
   Int_t istrip = strip->GetStrip();
   Float_t time = strip->GetBegTime();
   Int_t itime = static_cast<Int_t>((time*.1/18.7)+.5);

   Int_t navkey = (iplane<<22) | (istrip<<14) | itime;

   return navkey;

}

//______________________________________________________________________
NavKey CandStripHandle::KeyFromView(const CandStripHandle *strip)
{
  return static_cast<Int_t>(strip->GetStripEndId().GetPlaneView());
}

//______________________________________________________________________
NavKey CandStripHandle::StripSRKeyFromPSEId(const CandStripHandle *csh)
{ 
        return const_cast<CandStripHandle *>(csh)->GetStripEndId().GetEncoded(); 
}

//______________________________________________________________________
NavKey CandStripHandle::StripSRKeyFromBegTime(const CandStripHandle *csh)
{   
        Double_t time = (Double_t)csh->GetBegTime();
    return time;
}


XXXITRIMP(CandStripHandle)

---