ClusterAtNu.cxx

Go to the documentation of this file.

#include "ClusterAtNu.h"

ClassImp(ClusterAtNu)

//CVSID("$Id: ClusterAtNu.cxx,v 1.8 2005/02/04 21:07:54 blake Exp $");

ClusterAtNu::ClusterAtNu() :
  fPlane(-1),fBegStrip(-1), fEndStrip(-1), 
  fBegTime(0.0), fEndTime(0.0), 
  fZPos(0.0), fTPos(0.0), fCharge(0.0),
  fTrkFlag(0), fShwFlag(0),
  fTrkPlnFlag(0), fShwPlnFlag(0), fClrPlnFlag(0)
{

}

ClusterAtNu::ClusterAtNu(HitAtNu* hit) :
  fPlane(-1),fBegStrip(-1), fEndStrip(-1), 
  fBegTime(0.0), fEndTime(0.0), 
  fZPos(0.0), fTPos(0.0), fCharge(0.0),
  fTrkFlag(0), fShwFlag(0),
  fTrkPlnFlag(0), fShwPlnFlag(0), fClrPlnFlag(0)
{
  this->AddHit(hit);
}

ClusterAtNu::ClusterAtNu(const ClusterAtNu& rhs) :
  ObjAtNu(rhs),
  fPlane(rhs.fPlane),
  fBegStrip(rhs.fBegStrip),
  fEndStrip(rhs.fEndStrip),
  fBegTime(rhs.fBegTime),
  fEndTime(rhs.fEndTime),
  fZPos(rhs.fZPos),
  fTPos(rhs.fTPos),
  fCharge(rhs.fCharge),
  fTrkFlag(rhs.fTrkFlag),
  fShwFlag(rhs.fShwFlag),
  fTrkPlnFlag(rhs.fTrkPlnFlag),
  fShwPlnFlag(rhs.fShwPlnFlag),
  fClrPlnFlag(rhs.fClrPlnFlag)
{

}

ClusterAtNu::~ClusterAtNu()
{

}

Int_t ClusterAtNu::IsTrkAssoc(ClusterAtNu* clrm, ClusterAtNu* clrp) const
{
  Double_t win=99.9;
  Int_t assoc=0; 
  Double_t min=0.20,max=0.80;
  if( (fEndTime-clrm->GetBegTime()>-win && clrm->GetEndTime()-fBegTime>-win)
   && (fEndTime-clrp->GetBegTime()>-win && clrp->GetEndTime()-fBegTime>-win) ){

    Double_t clrmGetBegStrip=clrm->GetBegStrip();
    Double_t clrmGetEndStrip=clrm->GetEndStrip();
    if(fPlane-clrm->GetPlane()>2){
      Double_t km=2.0/(fPlane-clrm->GetPlane());
      clrmGetBegStrip=fBegStrip+km*(clrmGetBegStrip-fBegStrip);
      clrmGetEndStrip=fEndStrip+km*(clrmGetEndStrip-fEndStrip);
    }

    Double_t clrpGetBegStrip=clrp->GetBegStrip();
    Double_t clrpGetEndStrip=clrp->GetEndStrip();
    if(clrp->GetPlane()-fPlane>2){
      Double_t kp=2.0/(clrp->GetPlane()-fPlane);
      clrpGetBegStrip=fBegStrip+kp*(clrpGetBegStrip-fBegStrip);
      clrpGetEndStrip=fEndStrip+kp*(clrpGetEndStrip-fEndStrip);
    }

    Double_t k0 = 0.5*(clrp->GetPlane()-clrm->GetPlane()-4.0);
    min = min + 0.1*k0; max = max - 0.1*k0; 

    if(fEndStrip-clrm->GetBegStrip()>-1.1 && clrp->GetEndStrip()-fBegStrip>-1.1){
      if(fBegStrip-clrm->GetEndStrip()>-0.1 || clrp->GetBegStrip()-fEndStrip>-0.1){
        if( ( (clrp->GetBegStrip()-fEndStrip)-(fBegStrip-clrm->GetEndStrip())<2.1
          &&  (clrp->GetBegStrip()-fEndStrip)-(fBegStrip-clrm->GetEndStrip())>-2.1 )
         || ( min*clrmGetEndStrip+max*clrpGetEndStrip>fBegStrip-0.5
          &&  max*clrmGetBegStrip+min*clrpGetBegStrip<fEndStrip+0.5 ) )
          if(assoc<2) assoc=2; 
      }
      if(clrm->GetEndStrip()-fBegStrip>-1.1 && fEndStrip-clrp->GetBegStrip()>-1.1){
        if(assoc<1) assoc=1;
      }
    }
    if(fBegStrip-clrm->GetEndStrip()<1.1 && clrp->GetBegStrip()-fEndStrip<1.1){
      if(fEndStrip-clrm->GetBegStrip()<0.1 || clrp->GetEndStrip()-fBegStrip<0.1){
        if( ( (clrp->GetEndStrip()-fBegStrip)-(fEndStrip-clrm->GetBegStrip())>-2.1
          &&  (clrp->GetEndStrip()-fBegStrip)-(fEndStrip-clrm->GetBegStrip())<2.1 )
         || ( min*clrpGetEndStrip+max*clrmGetEndStrip>fBegStrip-0.5
          &&  max*clrpGetBegStrip+min*clrmGetBegStrip<fEndStrip+0.5 ) )
         if(assoc<2) assoc=2;
      }
      if(clrm->GetBegStrip()-fEndStrip<1.1 && fBegStrip-clrp->GetEndStrip()<1.1){
        if(assoc<1) assoc=1;
      }
    }

  }
  return assoc;
}

Int_t ClusterAtNu::IsShwAssoc(ClusterAtNu* clr) const
{
  Double_t win = 99.9;
  Int_t assoc = 0;
  if(fEndTime-clr->GetBegTime()>-win || clr->GetEndTime()-fBegTime>-win){
    if( clr->GetPlane()-fPlane<5 && clr->GetPlane()-fPlane>-5
     && clr->GetEndStrip()-fBegStrip>-6 && fEndStrip-clr->GetBegStrip()>-6 ){
       if( ( clr->GetPlane()-fPlane<3 && clr->GetPlane()-fPlane>-3
          && clr->GetEndStrip()-fBegStrip>-1 && fEndStrip-clr->GetBegStrip()>-1 )
        || ( clr->GetPlane()-fPlane==0 
          && clr->GetEndStrip()-fBegStrip>-3 && fEndStrip-clr->GetBegStrip()>-3 ) ) 
         assoc=2; else assoc=1;
    }
  }
  return assoc;
}

Int_t ClusterAtNu::IsDiffuseShwAssoc(ClusterAtNu* clr) const
{
  Double_t win = 99.9;
  Int_t assoc = 0;
  if(fEndTime-clr->GetBegTime()>-win || clr->GetEndTime()-fBegTime>-win){
    if( clr->GetPlane()-fPlane<9 && clr->GetPlane()-fPlane>-9
     && clr->GetEndStrip()-fBegStrip>-21 && fEndStrip-clr->GetBegStrip()>-21 ){
       if( ( clr->GetPlane()-fPlane<5 && clr->GetPlane()-fPlane>-5
        && clr->GetEndStrip()-fBegStrip>-11 && fEndStrip-clr->GetBegStrip()>-11 ) )
         assoc=2; else assoc=1;
    }
  }
  return assoc;
}

Int_t ClusterAtNu::IsHitAssoc(HitAtNu* hit) const
{
  Double_t win = 9999.9;
  if( hit->GetPlane()==fPlane 
   && hit->GetStrip()>fBegStrip-2 && hit->GetStrip()<fEndStrip+2
   && hit->GetTime()>fBegTime-win && hit->GetTime()<fEndTime+win )
    return 1; else return 0;
}

Int_t ClusterAtNu::GetPlane() const
{
  return fPlane;
}

Int_t ClusterAtNu::GetBegStrip() const
{
  return fBegStrip;
}
  
Int_t ClusterAtNu::GetEndStrip() const
{
  return fEndStrip;
}

Double_t ClusterAtNu::GetZPos()
{
  return fZPos;
}

Double_t ClusterAtNu::GetTPos()
{
  return fTPos;
}

Double_t ClusterAtNu::GetCharge()
{
  return fCharge;
}

Int_t ClusterAtNu::GetDigits()
{
  Int_t i,digits=0;
  for(i=0;i<1+this->GetHitLast();i++){
    HitAtNu* hit = (HitAtNu*)(this->GetHitAt(i));
    digits+=hit->GetDigits();
  }
  return digits;
}

Double_t ClusterAtNu::GetBegTime() const
{
  return fBegTime;
}

Double_t ClusterAtNu::GetEndTime() const
{
  return fEndTime;
}


Int_t ClusterAtNu::GetTrkFlag() const
{
  return fTrkFlag;
}
 
Int_t ClusterAtNu::GetShwFlag() const
{
  return fShwFlag;
}
 
Int_t ClusterAtNu::GetTrkPlnFlag() const
{
  return fTrkPlnFlag;
}
  
Int_t ClusterAtNu::GetShwPlnFlag() const
{
  return fShwPlnFlag;
}

Int_t ClusterAtNu::GetClrPlnFlag() const
{
  return fClrPlnFlag;
}

void ClusterAtNu::SetTrkFlag(Int_t flag)
{
  fTrkFlag=flag;
}
  
void ClusterAtNu::SetShwFlag(Int_t flag)
{
  fShwFlag=flag;
}

void ClusterAtNu::SetTrkPlnFlag(Int_t flag)
{
  fTrkPlnFlag=flag;
}

void ClusterAtNu::SetShwPlnFlag(Int_t flag)
{
  fShwPlnFlag=flag;
}

void ClusterAtNu::SetClrPlnFlag(Int_t flag)
{
  fClrPlnFlag=flag;
}

void ClusterAtNu::AddHit(HitAtNu* hit)
{
  if(this->AddHitToList(hit)){
    if(this->GetHitLast()==0){
      fPlane=hit->GetPlane();
      fBegStrip=hit->GetStrip();
      fEndStrip=hit->GetStrip();
      fBegTime=hit->GetTime();
      fEndTime=hit->GetTime();
      fZPos=hit->GetZPos();
    }
    else{
      if(hit->GetStrip()<fBegStrip) fBegStrip=hit->GetStrip();
      if(hit->GetStrip()>fEndStrip) fEndStrip=hit->GetStrip();
      if(hit->GetTime()<fBegTime) fBegTime=hit->GetTime();
      if(hit->GetTime()>fEndTime) fEndTime=hit->GetTime();
    }
    fTPos = (fTPos*fCharge+hit->GetTPos()*hit->GetCharge())/(fCharge+hit->GetCharge());
    fCharge += hit->GetCharge();
    return;
  }
  else return;
}

---