KVQ3.cpp

//Created by KVClassFactory on Tue Jan 13 15:11:11 2015
//Author: ,,,
 
#include "KVQ3.h"
#include "KVDBParameterList.h"
 
ClassImp(KVQ3)
 
// BEGIN_HTML <!--
/* -->
<h2>KVQ3</h2>
<h4>digitized charge signal</h4>
<!-- */
// --> END_HTML
 
 
void KVQ3::init()
{
   SetDefaultValues();
   fChannel = kQ3;
   SetType("Q3");
   fFPGAOutputNumbers = 2;
   LoadPSAParameters();
 
}
 
 
 
 
KVQ3::KVQ3()
{
   init();
}
 
 
 
 
 
KVQ3::KVQ3(const char* name) : KVFAZIASignal(name, "Charge")
{
   init();
}
 
 
 
 
void KVQ3::SetDefaultValues()
{
   SetChannelWidth(10);
   SetBaseLineLength(250);
}
 
 
 
 
void KVQ3::LoadPSAParameters()
{
 
   Double_t val = GetPSAParameter("BaseLineLength");
   SetBaseLineLength(val);
 
   val = GetPSAParameter("ChannelWidth");
   SetChannelWidth(val);
 
   val = GetPSAParameter("TauRC");
   SetTauRC(val);
 
   Double_t rise = GetPSAParameter("ShaperRiseTime");
   val = GetPSAParameter("ShaperFlatTop");
   SetTrapShaperParameters(rise, val);
 
   val = GetPSAParameter("PZCorrection");
   SetPoleZeroCorrection((val == 1));
 
   rise = GetPSAParameter("FastShaperRiseTime");
   val = GetPSAParameter("FastShaperFlatTop");
   SetFastTrapShaperParameters(rise, val);
 
   val = GetPSAParameter("MinimumAmplitude");
   SetAmplitudeTriggerValue(val);
 
}
 
 
 
 
void KVQ3::UpdatePSAParameter(KVNameValueList* par)
{
   for (Int_t ii = 0; ii < par->GetNpar(); ii += 1) {
      TString nameat(par->GetNameAt(ii));
      if (nameat == "BaseLineLength") SetBaseLineLength(par->GetDoubleValue(ii));
      else if (nameat == "ChannelWidth")    SetChannelWidth(par->GetDoubleValue(ii));
      else if (nameat == "ShaperRiseTime") SetShaperRiseTime(par->GetDoubleValue(ii));
      else if (nameat == "ShaperFlatTop")   SetShaperFlatTop(par->GetDoubleValue(ii));
      else if (nameat == "FastShaperRiseTime") SetFastShaperRiseTime(par->GetDoubleValue(ii));
      else if (nameat == "FastShaperFlatTop")  SetFastShaperFlatTop(par->GetDoubleValue(ii));
      else if (nameat == "TauRC")        SetTauRC(par->GetDoubleValue(ii));
      else if (nameat == "MinimumAmplitude")   SetAmplitudeTriggerValue(par->GetDoubleValue(ii));
      else if (nameat == "InterpolatedChannelWidth")    SetInterpolatedChannelWidth(par->GetDoubleValue(ii));
      else if (nameat == "Interpolation")   SetInterpolation((par->GetDoubleValue(ii) == 1));
      else if (nameat == "PZCorrection")    SetPoleZeroCorrection((par->GetDoubleValue(ii) == 1));
      else {
         if (nameat == "Detector" || nameat == "Signal" || nameat == "RunRange") {
 
         }
         else {
            Warning("UpdatePSAParameter", "Not supported PSA parameter : %d %s\n", ii, nameat.Data());
         }
      }
   }
}
 
 
 
 
 
void KVQ3::TreateSignal()
{
   if (PSAHasBeenComputed()) return;
 
   if (!IsLongEnough()) return;
   if (!TestWidth())
      ChangeChannelWidth(GetChannelWidth());
 
   RemoveBaseLine();
 
   FIR_ApplyTrapezoidal(fFastTrapRiseTime, fFastTrapFlatTop);
   ComputeAmplitude();
   fFastAmplitude = GetAmplitude();
 
 
 
   SetADCData();
 
   FIR_ApplyTrapezoidal(fTrapRiseTime, fTrapFlatTop);
   ComputeAmplitude();
   fAmplitude = GetAmplitude();// - (sQ3fast[fMtot-1])*0.5;
 
   ComputeRiseTime();
 
   fPSAIsDone = kTRUE;
 
}
 
 
 
 
void KVQ3::GetPSAResults(KVNameValueList&val) const
{
   if (!fPSAIsDone) return;
 
   val.SetValue("BaseLine", fBaseLine);
   val.SetValue("SigmaBaseLine", fSigmaBase);
   val.SetValue("Amplitude", fAmplitude);
   val.SetValue("FastAmplitude", fFastAmplitude);
   val.SetValue("RiseTime", fRiseTime);
   val.SetValue("RawAmplitude", GetRawAmplitude());
}