14 #include <TDirectory.h>
16 #include <TPaveLabel.h>
17 #include <TPaletteAxis.h>
89 void SC_PTC(
char const* input_filename)
92 TFile *
df =
new TFile(input_filename);
94 SC_PTC.open (
"SC_PTC.txt", std::ofstream::out);
98 for (
unsigned int j = 0; j <
NCHANNELS; j++)
101 PTC[j] =
new TCanvas( Form(
"PTC_%i",j+1), Form(
"PTC_%i",j+1), 800, 450);
103 PTC_BS[j] =
new TCanvas( Form(
"PTC_BS_%i",j+1), Form(
"PTC_BS_%i",j+1), 800, 450);
105 PTC_NS[j] =
new TCanvas( Form(
"PTC_NS_%i",j+1), Form(
"PTC_NS_%i",j+1), 800, 450);
108 TopDirectory = (TDirectory*) df->FindObjectAny(
"ST_Propagation_Time");
111 char* All = Form(
"h2_PpropTime_z_%i",j+1);
114 cout <<
"==================================================" << endl;
115 cout <<
"Processing Channel " << j+1 << endl;
116 cout <<
"==================================================" << endl;
121 for (Int_t i = 10; i < 40; i+=5)
124 cout <<
" k = " << k << endl;
126 bin1[j][i] =
h2All->GetXaxis()->FindBin(i);
127 bin2[j][i] =
h2All->GetXaxis()->FindBin(i+5);
131 double LL_fit = (
t_max[j][i] -100)*0.1 - 0.5 ;
132 double UL_fit = (
t_max[j][i] -100)*0.1 + 0.5 ;
134 h1_ss_Bin[j][i]->Fit(
"gaus",
"",
"",LL_fit,UL_fit);
135 TF1 *gaus =
h1_ss_Bin[j][i]->GetFunction(
"gaus");
146 for (Int_t i = 40; i < 43; i+=1)
156 double LL_fit = (
t_max_BS[j][i] -100)*0.1 - 0.5 ;
157 double UL_fit = (
t_max_BS[j][i] -100)*0.1 + 0.5 ;
159 h1_bs_Bin[j][i]->Fit(
"gaus",
"",
"",LL_fit,UL_fit);
160 TF1 *gaus_BS =
h1_bs_Bin[j][i]->GetFunction(
"gaus");
173 for (Int_t i = 43; i < 61; i+=2)
182 double LL_fit = (
t_max_NS[j][i] -100)*0.1 - 0.5 ;
183 double UL_fit = (
t_max_NS[j][i] -100)*0.1 + 0.5 ;
185 h1_ns_Bin[j][i]->Fit(
"gaus",
"",
"",LL_fit,UL_fit);
186 TF1 *gaus_NS =
h1_ns_Bin[j][i]->GetFunction(
"gaus");
Double_t t_z_ss_fit_intercept[NCHANNELS][3]
TH1D * h1_ns_Bin[30][100]
Double_t t_z_ss_fit_intercept_err[NCHANNELS][3]
Double_t t_z_bs_fit_slope_err[NCHANNELS][3]
void SC_PTC(char const *input_filename)
Double_t t_ns_fit_sigma_err[NCHANNELS][300][3]
Double_t t_max_NS[NCHANNELS][300]
Double_t t_z_bs_fit_intercept[NCHANNELS][3]
Double_t t_ss_fit_sigma_err[NCHANNELS][300][3]
TDirectory * TopDirectory
Double_t t_ss_fit_mean[NCHANNELS][300][3]
Double_t t_bs_fit_mean[NCHANNELS][300][3]
Double_t t_z_ns_fit_slope_err[NCHANNELS][3]
Double_t global[NCHANNELS]
Double_t t_bs_fit_sigma[NCHANNELS][300][3]
Double_t t_max_BS[NCHANNELS][300]
Double_t t_z_bs_fit_slope[NCHANNELS][3]
TH1D * h1_ss_Bin[30][100]
Double_t t_ss_fit_mean_err[NCHANNELS][300][3]
Double_t t_z_ss_fit_slope[NCHANNELS][3]
Double_t t_z_ns_fit_intercept[NCHANNELS][3]
Double_t t_bs_fit_sigma_err[NCHANNELS][300][3]
Double_t t_z_ns_fit_intercept_err[NCHANNELS][3]
Double_t t_bs_fit_mean_err[NCHANNELS][300][3]
Double_t t_z_ns_fit_slope[NCHANNELS][3]
Double_t t_ss_fit_sigma[NCHANNELS][300][3]
Double_t max_y_axis_value[NCHANNELS]
Double_t t_ns_fit_mean_err[NCHANNELS][300][3]
int sss[NCHANNELS][NOXbins][NOYbins]
Double_t N_y_bin[NCHANNELS]
Double_t t_z_bs_fit_intercept_err[NCHANNELS][3]
Double_t t_z_ss_fit_slope_err[NCHANNELS][3]
TH1D * h1_bs_Bin[30][100]
Double_t t_ns_fit_mean[NCHANNELS][300][3]
Double_t N_x_bin[NCHANNELS]
Double_t t_max[NCHANNELS]
Double_t ss_max[NCHANNELS]
Double_t t_ns_fit_sigma[NCHANNELS][300][3]
Double_t max_y_axis_bin[NCHANNELS]