DmxChiSqrStat Class Reference

#include <DmxChiSqrStat.h>

Inheritance diagram for DmxChiSqrStat:

List of all members.

Public Member Functions
	DmxChiSqrStat ()
	DmxChiSqrStat (Float_t *SignalW, Float_t *SignalE, Float_t CoG)
virtual	~DmxChiSqrStat ()
Float_t	GetGoodness () const
Float_t	GetTieBreaker () const
Float_t	GetPredictor () const
void	SetPredictor (Float_t predictor)
Private Attributes
Float_t	fChi2
Float_t	fRms
Float_t	fCoG
Float_t	fPredictor

---

Detailed Description

Definition at line 18 of file DmxChiSqrStat.h.

---

Constructor & Destructor Documentation

DmxChiSqrStat::DmxChiSqrStat

(

)

Definition at line 24 of file DmxChiSqrStat.cxx.

{
}

DmxChiSqrStat::DmxChiSqrStat	(	Float_t *	SignalW,
		Float_t *	SignalE,
		Float_t	CoG
	)

Definition at line 29 of file DmxChiSqrStat.cxx.

References fChi2, fPredictor, fRms, and DmxRMSStat::GetGoodness().

                                                                            :
  fChi2(0.),
  fRms(-1.),
  fPredictor(0.)
{

  Float_t sumW = 0.;
  Float_t sumE = 0.;

  //sum up the signal on sides W and E 

  for(Int_t i = 0; i < 192; i++ ){
    sumE += SignalE[i];
    sumW += SignalW[i];
  }
  
  //check to see if we can initialize the predictor as sumE / sumW.  if not
  //set it = 0.
  
  if(sumW > 0.) { fPredictor = sumE / sumW; }
    
  //MSG("DmxHyp", Msg::kInfo) << "predictor\t" << fPredictor << endl;

  //MSG("DmxHyp", Msg::kInfo) << "DmxChiSqrStat calculating goodness of fit" 
  //        << endl << endl;


  //this loop finds the chi^2 statistic for the signal and strip information
  //passed to the constructor.  the statistic is
  //
  //   chi^2 = SUM[ (RW_i - E_i)^2 / (W_i^2 + E_i^2) ]
  //
  //where R is fPredictor and W_i, E_i is the signal for strip i. 

  if( fPredictor > 0. ){
    for(Int_t j = 0; j < 192; j++){
      if(SignalW[j] > 0. || SignalE[j] > 0.){
        fChi2 += (((fPredictor * SignalW[j]) - SignalE[j]) * ((fPredictor * SignalW[j]) - SignalE[j])) / 
          ((SignalW[j] * SignalW[j]) + (SignalE[j] * SignalE[j]));
        //MSG("DmxHyp", Msg::kInfo) <<"strip = " << j << "\tchi^2 = " << fChi2 << "\tpredicter = " 
        //  << fPredictor << "\tW = " << SignalW[j] 
        //  << "\tE = " << SignalE[j] << endl;
      }
    }
  }
  else{ fChi2 = 10000.; }
  //MSG("DmxHyp", Msg::kInfo) << "chi^2 = " << fChi2 <<endl;

  //find the tiebreaker
  DmxRMSStat rms = DmxRMSStat(SignalW, SignalE, CoG);
   
  fRms = rms.GetGoodness();

  return;

}

DmxChiSqrStat::~DmxChiSqrStat

(

)

[virtual]

Definition at line 87 of file DmxChiSqrStat.cxx.

{
  //MSG("Dmx", Msg::kVerbose) << "deleting DmxChiSqrObject" << endl;
}

---

Member Function Documentation

Float_t DmxChiSqrStat::GetGoodness

(

)

const [virtual]

Implements DmxStatistic.

Definition at line 94 of file DmxChiSqrStat.cxx.

References fChi2.

{
  
  return fChi2;
}

Float_t DmxChiSqrStat::GetPredictor

(

)

const

Definition at line 101 of file DmxChiSqrStat.cxx.

References fPredictor.

{
  return fPredictor;
}

Float_t DmxChiSqrStat::GetTieBreaker

(

)

const [virtual]

Implements DmxStatistic.

Definition at line 111 of file DmxChiSqrStat.cxx.

References fRms.

{
  return fRms;
}

void DmxChiSqrStat::SetPredictor

(

Float_t

predictor

)

Definition at line 117 of file DmxChiSqrStat.cxx.

References fPredictor.

{
  fPredictor = predictor;
}

---

Member Data Documentation

Float_t DmxChiSqrStat::fChi2 [private]

Definition at line 35 of file DmxChiSqrStat.h.

Referenced by DmxChiSqrStat(), and GetGoodness().

Float_t DmxChiSqrStat::fCoG [private]

Definition at line 37 of file DmxChiSqrStat.h.

Float_t DmxChiSqrStat::fPredictor [private]

Definition at line 38 of file DmxChiSqrStat.h.

Referenced by DmxChiSqrStat(), GetPredictor(), and SetPredictor().

Float_t DmxChiSqrStat::fRms [private]

Definition at line 36 of file DmxChiSqrStat.h.

Referenced by DmxChiSqrStat(), and GetTieBreaker().

---

The documentation for this class was generated from the following files:

• DmxChiSqrStat.h
• DmxChiSqrStat.cxx

---