RungWire/check_two_sigma_violate.cxx File Reference

Go to the source code of this file.

Functions
void	check_two_sigma_violate (TString file="rungain.root")

Function Documentation

check_two_sigma_violate()

void check_two_sigma_violate

(

TString

file = "rungain.root"

)

Definition at line 1 of file RungWire/check_two_sigma_violate.cxx.

                                                         {
  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);
}