TFCS1DFunctionSpline.cxx File Reference

#include "ISF_FastCaloSimEvent/TFCS1DFunctionSpline.h"
#include "ISF_FastCaloSimEvent/TFCS1DFunctionInt32Histogram.h"
#include <algorithm>
#include <iostream>
#include "TMath.h"
#include "TCanvas.h"
#include "TH2F.h"
#include "TRandom.h"
#include "TFile.h"

Include dependency graph for TFCS1DFunctionSpline.cxx:

Go to the source code of this file.

Functions
void TFCS1DFunctionSpline::unit_test	ATLAS_NOT_THREAD_SAFE (TH1 *hist)

Function Documentation

ATLAS_NOT_THREAD_SAFE()

void TFCS1DFunctionSpline::unit_test ATLAS_NOT_THREAD_SAFE

(

TH1 *

hist

)

Definition at line 378 of file TFCS1DFunctionSpline.cxx.

                                                                    {
  ISF_FCS::MLogging logger;
  int nbinsx;
  TH1 *histfine = nullptr;
  if (hist == nullptr) {
    nbinsx = 50;
    double xmin = 1;
    double xpeak = 1.5;
    double sigma = 0.6;
    double xmax = 5;
    hist = new TH1D("test1D", "test1D", nbinsx, xmin, xmax);
    histfine = new TH1D("test1Dfine", "test1Dfine", 10 * nbinsx, xmin, xmax);
    hist->Sumw2();
    for (int i = 1; i <= 100000; ++i) {
      double x = gRandom->Gaus(xpeak, sigma);
      if (x >= xmin && x < xmax) {
        // hist->Fill(TMath::Sqrt(x));
        hist->Fill(x);
        if (histfine)
          histfine->Fill(x, 10);
      }
    }
  }
  if (!histfine)
    histfine = hist;
  TFCS1DFunctionSpline rtof(histfine, 0.01, 2, 15);
  nbinsx = hist->GetNbinsX();
 
  float value[2];
  float rnd[2];
  //cppcheck-suppress uninitvar
  for (rnd[0] = 0; rnd[0] < 0.9999; rnd[0] += 0.25) {
    rtof.rnd_to_fct(value, rnd);
    ATH_MSG_NOCLASS(logger, "rnd0=" << rnd[0] << " -> x=" << value[0]);
  }
 
  TH1 *hist_val = (TH1 *)histfine->Clone("hist_val");
  hist_val->SetTitle("toy simulation");
  hist_val->Reset();
  hist_val->SetLineColor(2);
  TH1 *hist_diff = (TH1 *)hist->Clone("hist_diff");
  hist_diff->SetTitle("difference");
  hist_diff->Reset();
  int nrnd = 5000000;
  double weight = histfine->Integral() / nrnd;
  double weightdiff = hist->Integral() / nrnd;
  hist_val->Sumw2();
  hist_diff->Sumw2();
  for (int i = 0; i < nrnd; ++i) {
    rnd[0] = gRandom->Rndm();
    rtof.rnd_to_fct(value, rnd);
    hist_val->Fill(value[0], weight);
    hist_diff->Fill(value[0], weightdiff);
  }
  hist_diff->Add(hist, -1);
 
  TH1F *hist_pull = new TH1F("pull", "pull", 200, -10, 10);
  for (int ix = 1; ix <= nbinsx; ++ix) {
    float val = hist_diff->GetBinContent(ix);
    float err = hist_diff->GetBinError(ix);
    if (err > 0)
      hist_pull->Fill(val / err);
    // ATH_MSG_NOCLASS(logger,"val="<<val<<" err="<<err);
  }
 
// Screen output in athena won't make sense and would require linking of
// additional libraries
#if defined(__FastCaloSimStandAlone__)
  new TCanvas("input", "Input");
  histfine->SetLineColor(kGray);
  histfine->Draw("hist");
  hist->Draw("same");
  hist_val->Draw("sameshist");
 
  new TCanvas("spline", "spline");
  TFCS1DFunctionInt32Histogram hist_fct(hist);
  int ngr = 101;
  TGraph *gr = new TGraph();
  for (int i = 0; i < ngr; ++i) {
    double r = i * 1.0 / (ngr - 1);
    gr->SetPoint(i, r, hist_fct.rnd_to_fct(r));
  }
  gr->SetMarkerStyle(7);
  gr->Draw("AP");
  TSpline3 *sp = new TSpline3(rtof.spline());
  sp->SetLineColor(2);
  sp->SetMarkerColor(2);
  sp->SetMarkerStyle(2);
  sp->Draw("LPsame");
 
  new TCanvas("difference", "difference");
  hist_diff->Draw();
 
  new TCanvas("pull", "Pull");
  hist_pull->Draw();
#endif
}