IceGOM.cxx

Go to the documentation of this file.

 
#include "IceGOM.h"
#include "Riostream.h"
 
ClassImp(IceGOM); // Class implementation to enable ROOT I/O
 
IceGOM::IceGOM() : NcDevice()
{
}

IceGOM::~IceGOM()
{
}

IceGOM::IceGOM(const IceGOM& m) : NcDevice(m)
{
}

Int_t IceGOM::GetString(Int_t id) const
{
 
 Int_t omid=GetUniqueID();
 if (id>0) omid=id;
 
 if (omid<=0) return 0;
 
 Int_t string=0;
 if (InheritsFrom("IceAOM"))
 {
  if (omid<=20) return -1;
  if (omid>=21 && omid<=40) return -2;
  if (omid>=41 && omid<=60) return -3;
  if (omid>=61 && omid<=86) return -4;
  if (omid>=87 && omid<=122) return -5;
  if (omid>=123 && omid<=158) return -6;
  if (omid>=159 && omid<=194) return -7;
  if (omid>=195 && omid<=230) return -8;
  if (omid>=231 && omid<=266) return -9;
  if (omid>=267 && omid<=302) return -10;
  if (omid>=303 && omid<=344) return -11;
  if (omid>=345 && omid<=386) return -12;
  if (omid>=387 && omid<=428) return -13;
  if (omid>=429 && omid<=470) return -14;
  if (omid>=471 && omid<=512) return -15;
  if (omid>=513 && omid<=554) return -16;
  if (omid>=555 && omid<=596) return -17;
  if (omid>=597 && omid<=638) return -18;
  if (omid>=639 && omid<=680) return -19;
 
  // OM681 is a special case on string 18
  if (omid==681) return -18;
 }
 else
 {
  string=omid/100;
 }
 
 return string;
}

Int_t IceGOM::GetLevel(Int_t id) const
{
 
 Int_t omid=GetUniqueID();
 if (id>0) omid=id;
 
 if (omid<=0) return 0;
 
 Int_t level=0;
 if (InheritsFrom("IceAOM"))
 {
  if (omid<=20) return omid;
  if (omid>=21 && omid<=40) return (omid-20);
  if (omid>=41 && omid<=60) return (omid-40);
  if (omid>=61 && omid<=86) return (omid-60);
  if (omid>=87 && omid<=122) return (omid-86);
  if (omid>=123 && omid<=158) return (omid-122);
  if (omid>=159 && omid<=194) return (omid-158);
  if (omid>=195 && omid<=230) return (omid-194);
  if (omid>=231 && omid<=266) return (omid-230);
  if (omid>=267 && omid<=302) return (omid-266);
  if (omid>=303 && omid<=344) return (omid-302);
  if (omid>=345 && omid<=386) return (omid-344);
  if (omid>=387 && omid<=428) return (omid-386);
  if (omid>=429 && omid<=470) return (omid-428);
  if (omid>=471 && omid<=512) return (omid-470);
  if (omid>=513 && omid<=554) return (omid-512);
  if (omid>=555 && omid<=596) return (omid-554);
  if (omid>=597 && omid<=638) return (omid-596);
  if (omid>=639 && omid<=680) return (omid-638);
 
  // OM681 is physically the 4th module on string 18
  // but the database convention is to regard it as
  // a module at the bottom of string 18
  if (omid==681) return 43;
 }
 else
 {
  level=omid%100;
 }
 
 return level;
}

Int_t IceGOM::GetOMId(Int_t string,Int_t level) const
{
 
 Int_t omid=0;
 Int_t s=0,l=0;
 if (InheritsFrom("IceAOM"))
 {
  for (Int_t i=1; i<=681; i++)
  {
   s=GetString(i);
   l=GetLevel(i);
   if (abs(s)==abs(string) && l==level) return i;
  }
 }
 else
 {
  omid=100*string+level;
 }
 
 return omid;
}

Double_t IceGOM::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;
 
 IceGOM* om=(IceGOM*)s->GetDevice();
 if (!om) return tres;
 
 // No implementation (yet) for IceTOP signals
 if (om->InheritsFrom("IceTDOM")) 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.31768387;      // Phase refractive index (c/v_phase) of ice
 const Float_t ngice=1.35075806;      // Group refractive index (c/v_group) of ice
 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=om->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("LE",7); // Hit time relative to the event time stamp
 
 tres=thit-tgeo;
 return tres;
}

Double_t IceGOM::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* IceGOM::Clone(const char* name) const
{
 
 IceGOM* m=new IceGOM(*this);
 if (name)
 {
  if (strlen(name)) m->SetName(name);
 }
 return m;
}