goodness.C

Go to the documentation of this file.

goodness()
{
  gStyle->SetOptStat(0);
  
  TCut rsqr1("sqrt((xf-xmc)**2+(yf-ymc)**2+(zf-zmc)**2)<10");
  TCut rsqr2("sqrt((xf-xmc)**2+(yf-ymc)**2+(zf-zmc)**2)<50");
  TCut rsqr3("sqrt((xf-xmc)**2+(yf-ymc)**2+(zf-zmc)**2)>50");

  TCut badf("e/(nhit-3-2*m)>15&&sumvia/nhit<0.5");

  TCanvas* c1 = new TCanvas("c1","c1",800,800);
  c1->Divide(2,2);
  c1->cd(1);
  gPad->SetLogy();
  TH1F* h11 = new TH1F("h11",";E/DoF;Events",70,0.0,70);
  TH1F* h12 = new TH1F("h12",";E/DoF;Events",70,0.0,70);
  TH1F* h13 = new TH1F("h13",";E/DoF;Events",70,0.0,70);
  nt->Project("h11","e/(nhit-3-2*m)",rsqr1);
  nt->Project("h12","e/(nhit-3-2*m)",rsqr2&&!rsqr1);
  nt->Project("h13","e/(nhit-3-2*m)",rsqr3);
  h11->SetMinimum(0.5);
  h11->Draw();
  h12->SetLineColor(kOrange);
  h12->Draw("same");
  h13->SetLineColor(kRed);
  h13->Draw("same");

  c1->cd(2);
  TH2F* h2 = new TH2F("h2",
                      "h2;#Sigma V_{ia} / N_{hit};E/DoF",
                      50,0,1,
                      70,0,70);
  h2->Draw();
  nt->SetMarkerColor(kBlack);
  nt->SetMarkerStyle(20);
  nt->SetMarkerSize(0.6);
  nt->Draw("e/(nhit-3-2*m):sumvia/nhit","","same");

  c1->cd(4);
  gPad->SetLogy(1);
  TH1F* h31 = new TH1F("h31",
                       "h31;#Sigma V_{ia} / N_{hit}",
                       50,0,1);
  TH1F* h32 = new TH1F("h32",
                       "h31;#Sigma V_{ia} / N_{hit}",
                       50,0,1);
  TH1F* h33 = new TH1F("h33",
                       "h31;#Sigma V_{ia} / N_{hit}",
                       50,0,1);
  nt->Project("h31","sumvia/nhit",rsqr1);
  nt->Project("h32","sumvia/nhit",rsqr2&&!rsqr1);
  nt->Project("h33","sumvia/nhit",rsqr3);
  h31->SetMinimum(0.5);
  h31->Draw();
  h32->SetLineColor(kOrange);
  h32->Draw("same");
  h33->SetLineColor(kRed);
  h33->Draw("same");
}
////////////////////////////////////////////////////////////////////////