RnoStation.cxx

Go to the documentation of this file.

 
#include "RnoStation.h"
#include "Riostream.h"
 
ClassImp(RnoStation); // Class implementation to enable ROOT I/O
 
RnoStation::RnoStation(const char* name,const char* title) : NcDetectorUnit(name,title)
{
 
 fCanvas=0;
}

RnoStation::~RnoStation()
{
 
 if (fCanvas)
 {
  if (gROOT->GetListOfCanvases()->FindObject(ClassName())) delete fCanvas;
  fCanvas=0;
 }
}

RnoStation::RnoStation(const RnoStation& q) : NcDetectorUnit(q)
{
 
 fCanvas=0;
}

TGraph* RnoStation::DisplaySampling(Int_t ich,Int_t j)
{
 
 if (fCanvas)
 {
  if (gROOT->GetListOfCanvases()->FindObject(ClassName())) delete fCanvas;
  fCanvas=0;
 }
 
 if (ich<0 || j<1) return 0;
 
 fCanvas=new TCanvas();
 TVirtualPad* pad=fCanvas->cd();
 gROOT->SetSelectedPad(pad);
 
 TString devname="Ch";
 devname+=ich;
 
 Int_t id=GetUniqueID(); // The station ID
 TString staname="Station";
 staname+=id;
 
 TString title="";
 
 title=staname;
 title+=" ";
 title+=devname;
 fCanvas->SetName(ClassName());
 fCanvas->SetTitle(title);
 fCanvas->SetGrid();
 
 GetSamplingGraph(ich,j).DrawClone("AL");
 
 // Obtain the pointer to the displayed graph
 TGraph* gr=(TGraph*)fCanvas->FindObject("NcSample");
 
 return gr;
}

TCanvas* RnoStation::DisplaySamplings(Int_t j)
{
 
 if (fCanvas)
 {
  if (gROOT->GetListOfCanvases()->FindObject(ClassName())) delete fCanvas;
  fCanvas=0;
 }
 
 if (j<1) return 0;
 
 fCanvas=new TCanvas();
 
 Int_t id=GetUniqueID(); // The station ID
 TString staname="Station";
 staname+=id;
 
 TString title="";
 TString devname="";
 
 title=staname;
 fCanvas->SetName(ClassName());
 fCanvas->SetTitle(title);
 fCanvas->Divide(4,6);
 
 TVirtualPad* pad=0;
 for (Int_t jch=0; jch<24; jch++)
 {
  pad=fCanvas->cd(jch+1);
 
  if (!pad) continue;
 
  gROOT->SetSelectedPad(pad);
  pad->SetGrid();
 
  GetSamplingGraph(jch,j).DrawClone("AL");
 } // End of loop over the channels
 
 return fCanvas;
}

TGraph RnoStation::GetSamplingGraph(Int_t ich,Int_t j)
{
 
 TGraph gr;
 
 if (ich<0 || j<1) return gr;
 
 TString devname="Ch";
 devname+=ich;
 
 Int_t id=GetUniqueID(); // The station ID
 TString staname="Station";
 staname+=id;
 
 NcDevice* daq=GetDevice("DAQ",kFALSE);
 Double_t fsample=0;
 if (daq) fsample=daq->GetSignal("Sampling-rate");
 
 TString rate="";
 rate.Form(" (DAQ: %-.3g Samples/sec)",fsample);
 
 TString title=staname;
 title+=" ";
 title+=devname;
 title+=rate;
 
 RnoGANT* ant=(RnoGANT*)GetDevice(devname,kTRUE);
 NcSample* sx=ant->GetSample(1);
 
 if (!sx) return gr;
 
 Int_t ndim=sx->GetDimension();
 if (j>ndim) return gr;
 
 gr=sx->GetGraph(j);
 
 TString xname="Sample";
 TString yname=sx->GetVariableName(j);
 
 gr.SetTitle(title);
 gr.GetXaxis()->SetTitle(xname);
 gr.GetYaxis()->SetTitle(yname);
 
 return gr;
}

TH1F RnoStation::GetSamplingDFT(Int_t ich,TString sel,Int_t j)
{
 
 TH1F his;
 
 TGraph gr=GetSamplingGraph(ich,j);
 
 NcDevice* daq=GetDevice("DAQ",kFALSE);
 Double_t fsample=0;
 if (daq) fsample=daq->GetSignal("Sampling-rate");
 
 NcDSP q;
 q.Load(&gr,fsample);
 q.Fourier("R2C",&his,sel);
 
 TString title=gr.GetTitle();
 
 his.SetTitle(title);
 his.SetStats(0);
 
 return his;
}

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