12 TH1F*
h_means =
new TH1F(
"means_dt",
"Mean time difference for each PSC counter;Counter (1-8 left, 9-16 right);#Deltat (PSC - RF) [ns]",16,1.0,17.0);
13 TH1F*
h_sigmas =
new TH1F(
"sigmas_dt",
"Sigmas for time difference for each PSC counter;Counter (1-8 left, 9-16);Sigma [ns]",16,1.0,17.0);
17 h_means->SetMinimum(-0.3);
18 h_means->SetMaximum(0.3);
19 h_means->GetYaxis()->SetTitleOffset(1.8);
20 h_sigmas->SetMinimum(0.08);
21 h_sigmas->SetMaximum(0.22);
22 h_sigmas->GetYaxis()->SetTitleOffset(1.8);
26 h_dt_vs_pp = (TH2I*)gROOT->FindObject(Form(
"h_dt_vs_pp_l_%i",i+1));
27 h_proj = h_dt_vs_pp->ProjectionY();
28 TFitResultPtr fitptr = h_proj->Fit(
"gaus",
"sqn");
29 double fitmean = fitptr->Parameters()[1];
30 double fitsig = fitptr->Parameters()[2];
32 fitptr = h_proj->Fit(
"gaus",
"sqn",
"",fitmean-2*fitsig,fitmean+2*fitsig);
33 fitmean = fitptr->Parameters()[1];
34 double fitmeanerr = fitptr->ParError(1);
35 fitsig = fitptr->Parameters()[2];
36 double fitsigerr = fitptr->ParError(2);
38 h_means->Fill(i+1,fitmean);
39 h_means->SetBinError(i+1,fitmeanerr);
40 h_sigmas->Fill(i+1,fitsig);
41 h_sigmas->SetBinError(i+1,fitsigerr);
45 h_dt_vs_pp = (TH2I*)gROOT->FindObject(Form(
"h_dt_vs_pp_r_%i",i+1));
46 h_proj = h_dt_vs_pp->ProjectionY();
47 fitptr = h_proj->Fit(
"gaus",
"sqn");
48 fitmean = fitptr->Parameters()[1];
49 fitsig = fitptr->Parameters()[2];
51 fitptr = h_proj->Fit(
"gaus",
"sqn",
"",fitmean-2*fitsig,fitmean+2*fitsig);
52 fitmean = fitptr->Parameters()[1];
53 fitmeanerr = fitptr->ParError(1);
54 fitsig = fitptr->Parameters()[2];
55 fitsigerr = fitptr->ParError(2);
57 h_means->Fill(i+NMODULES+1,fitmean);
58 h_means->SetBinError(i+NMODULES+1,fitmeanerr);
59 h_sigmas->Fill(i+NMODULES+1,fitsig);
60 h_sigmas->SetBinError(i+NMODULES+1,fitsigerr);
64 if (TVirtualPad::Pad() == NULL)
65 c =
new TCanvas(
"PSC_TW",
"PSC_TW",1200,600);
67 c = gPad->GetCanvas();
70 gStyle->SetOptStat(0);
73 TPad*
pad = (TPad*)c->GetPad(1);
74 pad->SetLeftMargin(0.15);
79 pad = (TPad*)c->GetPad(2);
80 pad->SetLeftMargin(0.15);