DmxRMSStat.cxx

Go to the documentation of this file.

//$Id: DmxRMSStat.cxx,v 1.19 2003/05/16 14:02:44 brebel Exp $
//
//DmxRMSStat.cxx
//
//DmxRMSStat is a Statistic class returning an RMS statistic for the
//demuxing of signal/strip possibilities
//
//Author:  B. Rebel 6/2000

#include "TMath.h"
#include "DeMux/DmxRMSStat.h"
#include "MessageService/MsgService.h"

ClassImp(DmxRMSStat)

//______________________________________________________________________
//CVSID("$Id: DmxRMSStat.cxx,v 1.19 2003/05/16 14:02:44 brebel Exp $");

//------------------------------------------------------------------
DmxRMSStat::DmxRMSStat()
{
}

//------------------------------------------------------------------
DmxRMSStat::DmxRMSStat(Float_t *SignalW, Float_t *SignalE, Float_t CoG) :
  fRms(0.)
{
  Float_t numerator = 0.;
  Float_t denominator = 0.;
  Float_t totalHits = 0.;

  for(Int_t i = 0; i < 192; i++){
    if( SignalW[i] > 0.){ 
      totalHits += 1.; 
      numerator += SignalW[i] * ((1. * i) - CoG) * ((1. * i) - CoG);
      denominator += SignalW[i];
    }
    if( SignalE[i] > 0.){ 
      totalHits += 1.; 
      numerator += SignalE[i] * ((1. * i) - CoG) * ((1. * i) - CoG);
      denominator += SignalE[i];
    }
  }
    
  //find the RMS and return it

  if(denominator != 0. && totalHits > 1.){
    Float_t rmssqr = (numerator / denominator) * (totalHits/(totalHits - 1.));
    fRms = TMath::Sqrt(rmssqr);
    //MSG("Dmx", Msg::kVerbose)  << "rms = " << rms << endl;
  }
  
  return;
}

//------------------------------------------------------------------
DmxRMSStat::~DmxRMSStat()
{
  //  MSG("Dmx", Msg::kVerbose) << "deleting DmxRMSStat object" << endl;
}

//------------------------------------------------------------------
//the rms of a reconstruction is defined by
//
//rms^2 = ((SUM[S * (x - x_avg)^2]) / (SUM[signal])) * (N / (N -1))
//
//where S is the ADC signal, x the strip number, x_avg the cog, and
//N the number of hits.

Float_t DmxRMSStat::GetGoodness() const
{
  
  return fRms;
}

//------------------------------------------------------------------

---