RnoDetector.cxx

Go to the documentation of this file.

 
#include "RnoDetector.h"
#include "Riostream.h"
 
ClassImp(RnoDetector); // Class implementation to enable ROOT I/O
 
RnoDetector::RnoDetector(const char* name,const char* title) : NcDetector(name,title)
{
 
 TString sname=name;
 TString stitle=title;
 
 if (sname=="") sname="RNO-G";
 if (stitle=="") stitle="The Radio Neutrino Observatory in Greenland";
 
 SetNameTitle(sname,stitle);
}

RnoDetector::~RnoDetector()
{
}

RnoDetector::RnoDetector(const RnoDetector& d) : NcDetector(d)
{
}

RnoStation* RnoDetector::GetStation(Int_t id,Bool_t create)
{
 
 // Only allow the "matrix element" labeled station IDs
 if (id<0 || !id%10)
 {
  cout << " *" << ClassName() << "::GetStation* Invalid station ID : " << id << endl;
  return 0;
 }
 
 RnoStation* station=0;
 
 TString name="Station";
 name+=id;
 
 station=(RnoStation*)GetDetectorUnit(name,kTRUE);
 
 if (!create) return station;
 
 // Creation of a new station
 CreateStation(id);
 
 station=(RnoStation*)GetDetectorUnit(name,kTRUE);
 
 return station;
}

void RnoDetector::CreateStation(Int_t id)
{
 
 RnoStation station("Station","RNO-G detector station");
 
 RnoString string("String","RNO-G string");
 
//@@@@ NcDetectorUnit trench("Trench","RNO-G surface trench");
 
 RnoULPDA ulpda("Channel","LPDA-up antenna");
 RnoDLPDA dlpda("Channel","LPDA-down antenna");
 RnoVPOL vpol("Channel","Vpol antenna");
 RnoHPOL hpol("Channel","Hpol antenna");
 
 Int_t nstr=3; // Number of strings per station
 Int_t ntre=3; // Number of surface trenches per station
 
 Int_t staid=id; // The unique station ID
 Int_t strid=0;  // The unique string ID
 Int_t antid=0;  // The unique antenna ID
 TString name="";
 Float_t depths[9]={-40,-60,-80,-92,-93,-94,-95,-96,-97}; // All antenna depths
 
 Int_t icha=0; // The RADIANT channel number
 
 name="Station";
 name+=staid;
 station.SetName(name);
 station.SetUniqueID(staid);
 
 // Loop over the strings of this station
 for (Int_t istr=1; istr<=nstr; istr++)
 {
  strid=1000*staid+100*istr;
  string.Reset();
  string.SetUniqueID(strid);
 
  name="String";
  name+=istr;
  string.SetName(name);
 
  // Power string
  if (istr==1) // Power string
  {
   for (Int_t ia=1; ia<=9; ia++)
   {
    antid=strid+ia;
    vpol.Reset();
    vpol.SetUniqueID(antid);
    vpol.SetPosition(0,-20,depths[ia-1],"car");
 
    hpol.Reset();
    hpol.SetUniqueID(antid);
    hpol.SetPosition(0,-20,depths[ia-1],"car");
 
    if (ia<4 || ia>5) // VPOL setting
    {
     icha=8-ia;
     if (ia>5) icha++;
     name="Ch";
     name+=icha;
     vpol.SetName(name);
     string.AddDevice(vpol);
    }
    if (ia==4 || ia==5) // HPOL setting
    {
     if (ia==4) icha=4;
     if (ia==5) icha=8;
     name="Ch";
     name+=icha;
     hpol.SetName(name);
     string.AddDevice(hpol);
    }
   }
  }
 
  // Helper string
  if (istr>1)
  {
   for (Int_t ia=7; ia<=9; ia++)
   {
    antid=strid+ia;
    vpol.Reset();
    vpol.SetUniqueID(antid);
    vpol.SetPosition(45-istr*15,10,depths[ia-1],"car");
 
    hpol.Reset();
    hpol.SetUniqueID(antid);
    hpol.SetPosition(45-istr*15,10,depths[ia-1],"car");
 
    icha=18-ia; // Helper string 1
    if (istr==3) icha=14+ia; // Helper string 2
 
    if (ia==7)
    {
     name="Ch";
     name+=icha;
     hpol.SetName(name);
     string.AddDevice(hpol);
    }
    else
    {
     name="Ch";
     name+=icha;
     vpol.SetName(name);
     string.AddDevice(vpol);
    }
   }
  }
  station.AddDevice(string);
 } // End of string loop
 
 // Loop over the surface trenches with surface antennas of this station
 icha=12;
 Int_t iadd=0;
 for (Int_t itre=1; itre<=ntre; itre++)
 {
  for (Int_t ia=1; ia<=3; ia++)
  {
   ulpda.Reset();
   ulpda.SetUniqueID(ia+iadd+1000*staid);
   ulpda.SetPosition(ia-2,itre,0,"car");
 
   dlpda.Reset();
   dlpda.SetUniqueID(ia+iadd+1000*staid);
   dlpda.SetPosition(ia,itre,0,"car");
 
   if (ia==1 || ia==3)
   {
    name="Ch";
    name+=icha;
    dlpda.SetName(name);
    station.AddDevice(dlpda);
   }
   else
   {
    name="Ch";
    name+=icha;
    ulpda.SetName(name);
    station.AddDevice(ulpda);
   }
   icha++;
  }
  iadd+=3;
 } // End of trench loop
 
 AddDevice(station);
}

TGraph* RnoDetector::DisplaySampling(Int_t ista,Int_t ich,Int_t j)
{
 
 RnoStation* sta=GetStation(ista,kFALSE);
 
 if (!sta)
 {
  cout << " *" << ClassName() << "::DisplayWaveform* Could not find Station" << ista << endl;
  return 0;
 }
 
 TGraph* gr=sta->DisplaySampling(ich,j);
 return gr;
}

TCanvas* RnoDetector::DisplaySamplings(Int_t ista,Int_t j)
{
 
 RnoStation* sta=GetStation(ista,kFALSE);
 
 if (!sta)
 {
  cout << " *" << ClassName() << "::DisplayWaveform* Could not find Station" << ista << endl;
  return 0;
 }
 
 TCanvas* c=sta->DisplaySamplings(j);
 return c;
}

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