FirstRung/check_two_sigma_violate.cxx

Go to the documentation of this file.

void check_two_sigma_violate(TString file="rungain.root"){
  TFile *f = new TFile(file);
  TTree *t = (TTree*)f->Get("runcalib");
  TH1F *h = new TH1F("h","h", 100, -0.5, 1.5);
  t->Draw("rungain>>h");
  h->Fit("gaus");
  double mean =  h->GetFunction("gaus")->GetParameter(1);
  double sigma =  h->GetFunction("gaus")->GetParameter(2);
  TString cut = Form("rungain<%f-2*%f || rungain>%f+2*%f",mean, sigma,mean, sigma);
  cout <<"cut: " <<  cut << endl;
  t->Scan("runno:runmean:rungain:runresol", cut);
}