test_shiftHistFit.cxx File Reference

#include "dunecore/DuneCommon/Utility/shiftHistFit.h"
#include "dunecore/DuneCommon/Utility/gausTF1.h"
#include "dunecore/DuneCommon/Utility/TPadManipulator.h"
#include <string>
#include <iostream>
#include <sstream>
#include <iomanip>
#include <vector>
#include <cassert>
#include "TH1F.h"
#include "TF1.h"
#include "TRandom.h"

Go to the source code of this file.

Typedefs
using	Index = unsigned int

Functions
int	test_shiftHistFit (double xshift)
int	main (int narg, const char *argc[])

Typedef Documentation

using Index = unsigned int

Definition at line 30 of file test_shiftHistFit.cxx.

Function Documentation

int main	(	int	narg,
		const char *	argc[]
	)

Definition at line 88 of file test_shiftHistFit.cxx.

                                       {
  double xshift = 1000.0;
  if ( narg > 1 ) {
    istringstream ssin(argc[1]);
    ssin >> xshift;
  }
  return test_shiftHistFit(xshift);
}

int test_shiftHistFit

(

double

xshift

)

Definition at line 34 of file test_shiftHistFit.cxx.

                                     {
  const string myname = "test_shiftHistFit: ";
  cout << myname << "Starting test" << endl;
  cout << myname << "Shift is " << xshift << endl;
#ifdef NDEBUG
  cout << myname << "NDEBUG must be off." << endl;
  abort();
#endif
  string line = "-----------------------------";
  string scfg;

  cout << myname << line << endl;
  cout << "Testing shifHistFit." << endl;

  cout << myname << line << endl;
  cout << "Build histogram" << endl;
  // Truth parameters before shift.
  double x0 = 50.0;
  double sig0 = 10.0;
  // Histogram range.
  double hx1 = x0 + xshift - 50.0;
  double hx2 = x0 + xshift + 50.0;
  // Fit range.
  double fx1 = x0 + xshift - 25.0;
  double fx2 = x0 + xshift + 25.0;
  // Fill histgram with data distributed around shifted truth.
  TH1* ph = new TH1F("hgau", "My gaus", 100, hx1, hx2);
  for ( int ient=0; ient<1000; ++ient ) {
    ph->Fill(gRandom->Gaus(x0 + xshift, sig0));
  }

  cout << myname << line << endl;
  cout << "Build function." << endl;
  TF1* pf = gausTF1(10.0, x0 + xshift + 0.3*sig0, 0.8*sig0, "mygaus");
  pf->SetRange(fx1, fx2);
  pf->Print();

  //ph->Fit(pf);
  shiftHistFit(ph, pf, "R", "Mean", xshift);
  pf->Print();

  TPadManipulator man;
  man.add(ph, "hist");
  man.add(pf, "same");
  man.addAxis();
  string fnout = "test_shiftHistFit.png";
  cout << myname << "Printing " << fnout << endl;
  man.print(fnout);

  return 0;
}