15 #include <TDirectory.h>
17 #include <TPaveLabel.h>
18 #include <TPaletteAxis.h>
75 Double_t fitval_pp = par[0] + par[1]*x[0];
82 TFile *
df =
new TFile(input_filename);
83 std::ofstream st_time_res;
84 st_time_res.open (
"st_time_res.txt", std::ofstream::out);
90 for (
unsigned int j = 0; j < 30; j++)
93 PT_can[j] =
new TCanvas( Form(
"PT_can_%i",j+1), Form(
"PT_can_%i",j+1), 800, 450);
97 char*
total = Form(
"h2_CorrectedTime_z_%i",j+1);
98 TH2I*
h2_total = (TH2I*) df->Get(total);
99 char* pj = Form(
"py_%i",j+1);
100 TH1D* pytotal = (TH1D*)pj;
102 cout <<
"==================================================" << endl;
103 cout <<
"Processing Channel " << j+1 << endl;
104 cout <<
"==================================================" << endl;
112 gStyle->SetOptStat(0);
113 gStyle->SetErrorX(0);
116 h2_total->Draw(
"colz");
118 h2_total->GetYaxis()->SetTitle(
"ST Corrected Time (ns)");
121 gStyle->SetOptFit(11);
122 gStyle->SetErrorX(0);
126 pytotal = h2_total->ProjectionY();
127 t_max[j]= pytotal->GetMaximum();
128 t_min[j]= pytotal->GetMinimum();
129 for (
unsigned int i = 0; i < 200; i++)
130 {
bins[i]= pytotal->GetBinContent(i);
133 pytotal->SetBinContent(i,0);
135 pytotal->GetXaxis()->SetRangeUser(-1.5,1.5);
136 pytotal->GetXaxis()->SetTitle(
"ST Corrected Time (ns)");
139 pytotal->Fit(
"gaus");
150 in.open(Form(
"st_time_res.txt"));
151 float sector,ser,tall,taller;
152 TNtuple *ntuple =
new TNtuple(
"ntuple",
"data from text file",
"sector:ser:tall:taller");
153 ntuple->ReadFile(
"st_time_res.txt");
157 Time_can =
new TCanvas(
"Time_can",
"Time_can", 800, 600);
159 gStyle->SetOptFit(011);
160 ntuple->Draw(
"sector:tall:taller:ser");
161 TGraphErrors *
gr=
new TGraphErrors(ntuple->GetSelectedRows(),ntuple->GetV1(),ntuple->GetV2(),ntuple->GetV4(),ntuple->GetV3());
163 gr->SetTitle(
"Time Resolution After Applying Propagation Time Corrections");
164 gr->GetXaxis()->SetTitle(
"Sector Number");
165 gr->GetYaxis()->SetTitle(
"Time Resolution (ps)");
166 gr->SetMarkerStyle(21);
167 gr->SetMarkerSize(3);
168 gr->SetMarkerColor(4);
172 TLine *line =
new TLine(0,350,33,350);
173 line->SetLineWidth(4);
174 line->SetLineColor(28);
175 line->SetLineStyle(9);
Double_t tns_min[NCHANNELS]
Double_t t_bs_fit_err[NCHANNELS][3]
Double_t t_bs_fit[NCHANNELS][3]
TF1 * t_vs_z_fit_chan[NCHANNELS]
Double_t tbs_min[NCHANNELS]
Double_t t_total_fit[NCHANNELS][3]
void st_time_resolution(char const *input_filename)
Double_t t_bn_fit[NCHANNELS][3]
Double_t t_min[NCHANNELS]
Double_t tbs_max[NCHANNELS]
Double_t t_ns_fit_err[NCHANNELS][3]
Double_t tall_max[NCHANNELS]
Double_t tss_min[NCHANNELS]
Double_t t_bn_fit_err[NCHANNELS][3]
Double_t t_total_fit_err[NCHANNELS][3]
Double_t t_all_fit[NCHANNELS][3]
Double_t t_ns_fit[NCHANNELS][3]
Double_t t_ss_fit_err[NCHANNELS][3]
Double_t tns_max[NCHANNELS]
Double_t tbn_max[NCHANNELS]
Double_t t_ss_fit[NCHANNELS][3]
Double_t t_all_fit_err[NCHANNELS][3]
Double_t tall_min[NCHANNELS]
Double_t t_max[NCHANNELS]
Double_t tss_max[NCHANNELS]
Double_t fitf_pp(Double_t *x, Double_t *par)
Double_t tbn_min[NCHANNELS]