RnoGANT.cxx

Go to the documentation of this file.

 
#include "RnoGANT.h"
#include "Riostream.h"
 
ClassImp(RnoGANT); // Class implementation to enable ROOT I/O
 
RnoGANT::RnoGANT(const char* name,const char* title) : NcDevice(name,title)
{
}

RnoGANT::~RnoGANT()
{
}

RnoGANT::RnoGANT(const RnoGANT& q) : NcDevice(q)
{
}

Int_t RnoGANT::GetStation(Int_t id) const
{
 
 Int_t antid=GetUniqueID();
 if (id>0) antid=id;
 
 if (antid<=0) return -1;
 
 Int_t station=-1;
 
 station=antid/1000;
 
 return station;
}

Int_t RnoGANT::GetString(Int_t id) const
{
 
 Int_t antid=GetUniqueID();
 if (id>0) antid=id;
 
 if (antid<=0) return -1;
 
 Int_t string=-1;
 
 if (InheritsFrom("RnoSANT"))
 {
  string=0;
 }
 else
 {
  antid=antid%1000;
  string=antid/100;
 }
 
 return string;
}

Int_t RnoGANT::GetNumber(Int_t id) const
{
 
 Int_t number=-1;
 
 // DAQ channel number was requested
 if (id<0)
 {
  TString name=GetName();
 
  if (!name.Contains("Ch")) return -1;
 
  name.ReplaceAll("Ch","");
  number=name.Atoi();
  return number;
 }
 
 // Antenna number was requested 
 Int_t antid=GetUniqueID();
 if (id>0) antid=id;
 
 Int_t station=GetStation(antid);
 Int_t string=GetString(antid);
 
 if (station<0 || string<0) return -1;
 
 number=antid%100;
 
 return number;
}

Int_t RnoGANT::GetANTId(Int_t station,Int_t string,Int_t number) const
{
 
 Int_t antid=-1;
 
 if (station<11 || station>77 || string<1 || string>3 || number<=0 || number>9) return -1;
 
 if (!(station%10)) return -1; // Station IDs must represent the matrix (col,row) indexing
 
 if ((string==2 || string==3) && number<7) return -1; // Helper string antenna numbers start at 7
 
 antid=1000*station+100*string+number;
 
 return antid;
}

Double_t RnoGANT::GetTimeResidual(NcEvent* evt,NcTrack* t,NcSignal* s,TString name,Int_t mode,Int_t vgroup) const
{
 
 Double_t tres=-99999;
 
 if (!evt || !t || !s) return tres;
 
 RnoGANT* ant=(RnoGANT*)s->GetDevice();
 if (!ant) return tres;
 
 Nc3Vector p=t->Get3Momentum();
 if (!p.HasVector() || !p.GetNorm()) return tres;
 
 const Float_t pi=acos(-1.);
 const Float_t c=0.299792458;         // Light speed in vacuum in meters per ns
 const Float_t npice=1.78829165;      // Phase refractive index (c/v_phase) of ice (IceCube optical 1.31768387)
 const Float_t ngice=1.0251*npice;    // Group refractive index (c/v_group) of ice (IceCube optical 1.35075806)
 const Float_t thetac=acos(1./npice); // Cherenkov angle (in radians)
 
 // Angular reduction of complement of thetac due to v_phase and v_group difference
 Float_t alphac=0;
 if (vgroup) alphac=atan((1.-npice/ngice)/sqrt(npice*npice-1.));
 
 // The track reference point
 NcPosition* r0=t->GetReferencePoint();
 if (!r0) return tres;
 
 // Absolute (UT) time stamp of the track reference point
 NcTimestamp* tt0=r0->GetTimestamp();
 if (!tt0) return tres;
 
 // Time stamp of the track relative to the event time stamp 
 Float_t t0=evt->GetDifference(tt0,"ns");
 
 NcPosition rhit=ant->GetPosition();
 Float_t d=t->GetDistance(rhit);
 Nc3Vector r12=rhit-(*r0);
 Float_t hproj=p.Dot(r12)/p.GetNorm();
 Float_t dist=fabs(hproj)+d/tan(pi/2.-thetac-alphac);
 if (hproj<0) dist=-dist;
 Float_t tgeo=t0+dist/c;               // The predicted geometrical hit time
 Float_t thit=s->GetSignal(name,mode); // Hit time relative to the event time stamp
 
 tres=thit-tgeo;
 return tres;
}

Double_t RnoGANT::GetTimeResidual(NcEvent* evt,NcTrack* t,Int_t j,TString name,Int_t mode,Int_t vgroup) const
{
 
 NcSignal* s=GetHit(j);
 Double_t tres=GetTimeResidual(evt,t,s,name,mode,vgroup);
 return tres;
}

TObject* RnoGANT::Clone(const char* name) const
{
 
 RnoGANT* q=new RnoGANT(*this);
 if (name)
 {
  if (strlen(name)) q->SetName(name);
 }
 return q;
}