---

Comp2RefKolmogorov Class Reference

Comparison to reference using the Kolmogorov algorithm. More...

#include <DQMServices/Core/interface/QTest.h>

Inheritance diagram for Comp2RefKolmogorov:

List of all members.

Public Member Functions
	Comp2RefKolmogorov (const std::string &name)
float	runTest (const MonitorElement *me)
Static Public Member Functions
static std::string	getAlgoName (void)
Protected Attributes
TH1 *	h
Int_t	ncx1
	# of bins for test & reference histograms
Int_t	ncx2
TH1 *	ref_
Static Protected Attributes
static const Double_t	difprec = 1e-5

---

Detailed Description

Comparison to reference using the Kolmogorov algorithm.

Definition at line 261 of file QTest.h.

---

Constructor & Destructor Documentation

Comp2RefKolmogorov::Comp2RefKolmogorov

(

const std::string &

name

)

[inline]

Definition at line 264 of file QTest.h.

References getAlgoName(), and QCriterion::setAlgoName().

                                            : SimpleTest(name)
  { setAlgoName(getAlgoName()); }

---

Member Function Documentation

static std::string Comp2RefKolmogorov::getAlgoName

(

void

)

[inline, static]

Definition at line 269 of file QTest.h.

Referenced by Comp2RefKolmogorov(), QTestConfigure::EnableComp2RefKolmogorovTest(), QTestConfigure::enableTests(), and QTestParameterNames::QTestParameterNames().

  { return "Comp2RefKolmogorov"; }

float Comp2RefKolmogorov::runTest

(

const MonitorElement *

me

)

[virtual]

Reimplemented from QCriterion.

Definition at line 254 of file QTest.cc.

References GenMuonPlsPt100GeV_cfg::cout, difprec, MonitorElement::DQM_KIND_TH1F, MonitorElement::DQM_KIND_TH1S, MonitorElement::DQM_KIND_TPROFILE, lat::endl(), first, MonitorElement::getRefRootObject(), MonitorElement::getRefTH1F(), MonitorElement::getRefTH1S(), MonitorElement::getRefTProfile(), MonitorElement::getRootObject(), MonitorElement::getTH1F(), MonitorElement::getTH1S(), MonitorElement::getTProfile(), h, MonitorElement::kind(), prof2calltree::last, ncx1, ncx2, ref_, s1, s2, w1, w2, and z.

Referenced by FourVectorHLTClient::endRun().

{
   if (!me) return -1;
   if (!me->getRootObject() || !me->getRefRootObject()) return -1;
   h=0;
   ref_=0;

   //-- TH1F
   if (me->kind()==MonitorElement::DQM_KIND_TH1F){ 
    h = me->getTH1F(); // access Test histo
    ref_ = me->getRefTH1F(); //access Ref histo
   } 
   //-- TH1S
   else if (me->kind()==MonitorElement::DQM_KIND_TH1S){ 
    h = me->getTH1S(); // access Test histo
    ref_ = me->getRefTH1S(); //access Ref histo
   } 
   //-- TProfile
   else if (me->kind()==MonitorElement::DQM_KIND_TPROFILE){
    h = me->getTProfile(); // access Test histo
    ref_ = me->getRefTProfile(); //access Ref histo
    }
   else{ 
    std::cout<< "Comp2RefKolmogorov ERROR: ME does not contain TH1F/TH1S/TProfile" << std::endl; 
    return -1;
   } 
   
   //-- isInvalid ? - Check consistency in number of channels
  ncx1  = h->GetXaxis()->GetNbins(); 
  ncx2   = ref_->GetXaxis()->GetNbins();
  if ( ncx1 !=  ncx2){
  std::cout<<"Comp2RefKolmogorov ERROR: different number of channels! ("
  << ncx1 << ", " << ncx2 << ") " << std::endl;
  return -1;
  } 
  //-- isInvalid ? - Check consistency in channel edges
  Double_t diff1 = TMath::Abs(h->GetXaxis()->GetXmin() - ref_->GetXaxis()->GetXmin());
  Double_t diff2 = TMath::Abs(h->GetXaxis()->GetXmax() - ref_->GetXaxis()->GetXmax());
  if (diff1 > difprec || diff2 > difprec){
  std::cout << "Comp2RefKolmogorov ERROR: histograms with different binning" << std::endl;
  return -1;
  }

   //--  QUALITY TEST itself 
  Bool_t afunc1 = kFALSE; Bool_t afunc2 = kFALSE;
  Double_t sum1 = 0, sum2 = 0;
  Double_t ew1, ew2, w1 = 0, w2 = 0;
  Int_t bin;
  for (bin=1;bin<=ncx1;bin++){
    sum1 += h->GetBinContent(bin);
    sum2 += ref_->GetBinContent(bin);
    ew1   = h->GetBinError(bin);
    ew2   = ref_->GetBinError(bin);
    w1   += ew1*ew1;
    w2   += ew2*ew2;
  }
  if (sum1 == 0){
    std::cout << "Comp2RefKolmogorov ERROR: Test Histogram " << h->GetName() << " integral is zero" << std::endl;
    return -1;
  }
  if (sum2 == 0){
    std::cout << "Comp2RefKolmogorov ERROR: Ref Histogram " << ref_->GetName() << " integral is zero" << std::endl;
    return -1;
  }

  Double_t tsum1 = sum1; Double_t tsum2 = sum2;
  tsum1 += h->GetBinContent(0);
  tsum2 += ref_->GetBinContent(0);
  tsum1 += h->GetBinContent(ncx1+1);
  tsum2 += ref_->GetBinContent(ncx1+1);

  // Check if histograms are weighted.
  // If number of entries = number of channels, probably histograms were
  // not filled via Fill(), but via SetBinContent()
  Double_t ne1 = h->GetEntries();
  Double_t ne2 = ref_->GetEntries();
  // look at first histogram
  Double_t difsum1 = (ne1-tsum1)/tsum1;
  Double_t esum1 = sum1;
  if (difsum1 > difprec && Int_t(ne1) != ncx1)
  {
    if (h->GetSumw2N() == 0)
      std::cout << " Comp2RefKolmogorov WARNING: Weighted events and no Sumw2 for "
                << h->GetName() << std::endl;
    else
      esum1 = sum1*sum1/w1;  //number of equivalent entries
  }
  // look at second histogram
  Double_t difsum2 = (ne2-tsum2)/tsum2;
  Double_t esum2   = sum2;
  if (difsum2 > difprec && Int_t(ne2) != ncx1)
  {
    if (ref_->GetSumw2N() == 0)
      std::cout << " Comp2RefKolmogorov WARNING: Weighted events and no Sumw2 for "
                << ref_->GetName() << std::endl;
    else
      esum2 = sum2*sum2/w2;  //number of equivalent entries
  }

  Double_t s1 = 1/tsum1; Double_t s2 = 1/tsum2;

  // Find largest difference for Kolmogorov Test
  Double_t dfmax =0, rsum1 = 0, rsum2 = 0;

  // use underflow bin
  Int_t first = 0; // 1
  // use overflow bin
  Int_t last  = ncx1+1; // ncx1
  for (bin=first;bin<=last;bin++)
  {
    rsum1 += s1*h->GetBinContent(bin);
    rsum2 += s2*ref_->GetBinContent(bin);
    dfmax = TMath::Max(dfmax,TMath::Abs(rsum1-rsum2));
  }

  // Get Kolmogorov probability
  Double_t z = 0;
  if (afunc1)      z = dfmax*TMath::Sqrt(esum2);
  else if (afunc2) z = dfmax*TMath::Sqrt(esum1);
  else             z = dfmax*TMath::Sqrt(esum1*esum2/(esum1+esum2));

  // This numerical error condition should never occur:
  if (TMath::Abs(rsum1-1) > 0.002)
    std::cout << " Comp2RefKolmogorov WARNING: Numerical problems with histogram "
              << h->GetName() << std::endl;
  if (TMath::Abs(rsum2-1) > 0.002)
    std::cout << " Comp2RefKolmogorov WARNING: Numerical problems with histogram "
              << ref_->GetName() << std::endl;

  return TMath::KolmogorovProb(z);
}

---

Member Data Documentation

const Double_t Comp2RefKolmogorov::difprec = 1e-5 [static, protected]

Definition at line 279 of file QTest.h.

Referenced by runTest().

TH1* Comp2RefKolmogorov::h [protected]

Definition at line 274 of file QTest.h.

Referenced by runTest().

Int_t Comp2RefKolmogorov::ncx1 [protected]

# of bins for test & reference histograms

Definition at line 278 of file QTest.h.

Referenced by runTest().

Int_t Comp2RefKolmogorov::ncx2 [protected]

Definition at line 278 of file QTest.h.

Referenced by runTest().

TH1* Comp2RefKolmogorov::ref_ [protected]

Definition at line 275 of file QTest.h.

Referenced by runTest().

---

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

• DQMServices/Core/interface/QTest.h
• DQMServices/Core/src/QTest.cc

---