HitCamAtNu.cxx

Go to the documentation of this file.

// Package: AtNuReco
//
// HitCamAtNu.cxx
//
// chapman@hep.phy.cam.ac.uk
// blake@hep.phy.cam.ac.uk
#include "AtNuReco/HitCamAtNu.h"
#include "Validity/VldContext.h"

#include <cmath>
#include <cstdlib>

ClassImp(HitCamAtNu)
  
  //CVSID("$Id: HitCamAtNu.cxx,v 1.2 2009/02/28 21:46:10 gmieg Exp $");


HitCamAtNu::HitCamAtNu(CandStripHandle* strip) :
   fStrip(0), fUid(0), 
   fPlaneView(-1), fStripNum(-1), fTime(-999.), fTrackFlag(0), 
   fShowerFlag(0), fCharge(0.), fPlane(-1), fTPos(-1), fZPos(-1)
{
  if(strip) {
    fStrip=strip;
    fTime=1.0e9*strip->GetTime();
    fStripNum = strip->GetStrip();
    fCharge = strip->GetCharge();
    fPlane = strip->GetPlane();
    fTPos = strip->GetTPos();
    fZPos = strip->GetZPos();
    
    if( strip->GetPlaneView()==PlaneView::kU
        || strip->GetPlaneView()==PlaneView::kX
        || strip->GetPlaneView()==PlaneView::kA ) {fPlaneView=0;}
    if( strip->GetPlaneView()==PlaneView::kV
        || strip->GetPlaneView()==PlaneView::kY
        || strip->GetPlaneView()==PlaneView::kB ) {fPlaneView=1;}
  }
  
}


HitCamAtNu::~HitCamAtNu()
{
}


int HitCamAtNu::IsShwAssoc(HitCamAtNu* hit) const
{
  double win=99.9;
  if(fStrip!=0)
    {
      if( fabs(hit->GetTime()-this->GetTime())<win )
        {
          if( abs(hit->GetPlane()-this->GetPlane())<3 && abs(hit->GetStrip()-this->GetStrip())<4 )
            {
              return 2;
            }
          else if( abs(hit->GetPlane()-this->GetPlane())<5 && abs(hit->GetStrip()-this->GetStrip())<6 )
            {
              return 1;
            }
        }
    }
  return 0;
}


int HitCamAtNu::IsDiffuseShwAssoc(HitCamAtNu* hit) const
{
  double win=99.9;
  if(fStrip!=0)
    {
      if( fabs(hit->GetTime()-this->GetTime())<win )
        {
          if( abs(hit->GetPlane()-this->GetPlane())<5 && abs(hit->GetStrip()-this->GetStrip())<11 )
            {
              return 2;
            }
          else if( abs(hit->GetPlane()-this->GetPlane())<9 && abs(hit->GetStrip()-this->GetStrip())<21 )
            {
              return 1;
            }
        }
    }
  return 0;
}

---