TEfficiencyCal.cxx

Go to the documentation of this file.

#include "TEfficiencyCal.h"
 
TEfficiencyCal::TEfficiencyCal(const TEfficiencyCal& copy) : TCal(copy)
{
   copy.Copy(*this);
}
 
void TEfficiencyCal::Copy(TObject& obj) const
{
   static_cast<TEfficiencyCal&>(obj).fScaleFactor = fScaleFactor;
   TCal::Copy(obj);
}
 
void TEfficiencyCal::ScaleGraph(Double_t scaleFactor)
{
   if(GetN() == 0) {
      Error("ScaleGraph", "Graph does not exist");
      return;
   }
 
   for(int i = 0; i < GetN(); i++) {
      GetY()[i] *= scaleFactor;
      GetEY()[i] *= scaleFactor;
   }
 
   fScaleFactor = scaleFactor;
}
 
void TEfficiencyCal::AddPoint(TPeak* peak)
{
   if(peak == nullptr) {
      Error("AddPoint", "Peak is empty");
      return;
   }
   AddPoint(peak->GetCentroid(), peak->GetArea(), peak->GetCentroidErr(), peak->GetAreaErr());
}
 
void TEfficiencyCal::AddPoint(Double_t energy, Double_t area, Double_t dEnergy, Double_t dArea)
{
   if(GetNucleus() == nullptr) {
      Error("AddPoint", "No nucleus set");
      return;
   }
   // Will eventually write a method that doesn't need a nucleus
 
   Double_t intensity = 1.0;   // nuc;
 
   Double_t efficiency  = area / intensity;
   Double_t dEfficiency = dArea / intensity;
   // dEfficiency = efficiency*TMath::Sqrt(TMath::Power(dEfficiency/efficiency,2.0) + TMath::Power(dArea/area,2.0));
 
   SetPoint(GetN(), energy, efficiency);
   SetPointError(GetN() - 1, dEnergy, dEfficiency);
 
   Sort();   // This keeps the points in order of energy;
}
 
void TEfficiencyCal::Print(Option_t*) const
{
   TCal::Print();
}
 
void TEfficiencyCal::Clear(Option_t*)
{
   fScaleFactor = 1.0;
   TCal::Clear();
}