KVFAZIASignal.cpp

//Created by KVClassFactory on Mon Dec 22 15:46:46 2014
//Author: ,,,
 
#include "KVFAZIASignal.h"
#include "KVString.h"
#include "TMath.h"
#include "KVDigitalFilter.h"
#include "KVEnv.h"
#include "KVDataSet.h"
#include "TMatrixD.h"
#include "TMatrixF.h"
#include "TClass.h"
 
#define LOG2 (double)6.93147180559945286e-01
# define M_PI     3.14159265358979323846  /* pi */
 
ClassImp(KVFAZIASignal)
 
// BEGIN_HTML <!--
/* -->
<h2>KVFAZIASignal</h2>
<h4>Base class for FAZIA signal processing</h4>
<!-- */
// --> END_HTML
//
// SIGNAL NAME
// ===========
// KVFAZIASignal::GetName() returns the name of the signal as defined in FAZIA
// raw data, e.g. "T1-Q3-B001-QH1"
//
// SIGNAL TYPE
// ===========
// KVFAZIASignal::GetType() returns one of: "QH1", "QL1", "Q2", "Q3", "I1", "I2"
 
 
 
void KVFAZIASignal::init()
{
   fChannel = kUNKDT;
   ResetIndexes();
   SetDefaultValues();
}
 
 
 
 
void KVFAZIASignal::ResetIndexes()
{
   fIndex = 0;
   fDetName = "";
   fFPGAOutputNumbers = 0;
}
 
 
 
 
KVFAZIASignal::KVFAZIASignal()
{
   // Default constructor
   init();
}
 
 
 
 
 
KVFAZIASignal::KVFAZIASignal(const char* name, const char* title) : KVSignal(name,title)
{
   init();
}
 
 
 
 
 
KVFAZIASignal::KVFAZIASignal(const TString& name, const TString& title) : KVSignal(name,title)
{
   init();
}
 
 
 
 
 
 
 
void KVFAZIASignal::TreateOldSignalName()
{
   TString stit = GetTitle();
   stit.ToUpper();
 
   KVString tmp = GetName();
   KVString part = "";
   if (tmp.BeginsWith("B")) {
      tmp.Begin("-");
      part = tmp.Next();
      part.ReplaceAll("B", "");
      Int_t bb = part.Atoi();
      part = tmp.Next();
      part.ReplaceAll("Q", "");
      Int_t qq = part.Atoi();
      part = tmp.Next();
      part.ReplaceAll("T", "");
      Int_t tt = part.Atoi();
      fType = tmp.Next();
 
      fIndex = 100 * bb + 10 * qq + tt;
      fDetName.Form("%s-%d", stit.Data(), fIndex);
   }
   else if (tmp.BeginsWith("RUTH")) {
      //Old FAZIA telescope denomination
      //    Rutherford telescope case for FAZIASYM experiment
      tmp.Begin("-");
      part = tmp.Next();
      fType = tmp.Next();
 
      TString stit = GetTitle();
      stit.ToUpper();
 
      fDetName.Form("%s-RUTH", stit.Data());
   }
}
 
 
 
 
void KVFAZIASignal::DeduceFromName()
{
   //methods used to identify to from which detector/telescope/quartet/block
   //the signal is associated
   //it is assumed that the name of the signal is defined as it is in the raw data
   //generated by fazia_reader which convert raw acquisition file in raw ROOT files
 
   ResetIndexes();
   KVString tmp = GetName();
   KVString part = "";
   if (tmp.BeginsWith("B") || tmp.BeginsWith("RUTH")) {
      //Old FAZIA telescope denomination
      //Info("DeduceFromName","Old format %s",GetName());
      TreateOldSignalName();
   }
   else {
 
      if (tmp.GetNValues("-") == 2) {
         //new general denomination
         //in KaliVeda :
         //    signal name : [signal_type]-[index]
         //    if  index is digit (number), it is the telescope number : 100*b+10*q+t
         //    if index is a string of character, we kept it as it is
         //          for example RUTH for the FAZIASYM dataset
         TString stit = GetTitle();
         stit.ToUpper();
         //new format
         //Info("DeduceFromName", "New format %s", GetName());
         tmp.Begin("-");
         fType = tmp.Next();
         KVString ss = tmp.Next();
 
         if (ss.IsDigit()) {
            fIndex = ss.Atoi();
            fDetName.Form("%s-%d", stit.Data(), fIndex);
         }
         else {
            fDetName.Form("%s-%s", stit.Data(), ss.Data());
         }
      }
      else {
         Warning("DeduceFromName", "Unkown format %s", GetName());
      }
   }
 
 
}
 
 
 
 
Double_t KVFAZIASignal::GetPSAParameter(const Char_t* parname)
{
   //DeduceFromName has to be called before
 
   Double_t lval = -1;
   KVString spar;
   spar.Form("%s.%s", GetType(), parname);
   if (gDataSet)  lval = gDataSet->GetDataSetEnv(spar.Data(), 0.0);
   else           lval = gEnv->GetValue(spar.Data(), 0.0);
   return lval;
}
 
 
 
 
void KVFAZIASignal::Print(Option_t*) const
{
   Info("Print", "\nName: %s - Title: %s", GetName(), GetTitle());
   if (fDetName != "") {
      printf("\tAssociated to the detector %s\n", fDetName.Data());
   }
   print_psa_parameters();
}