10 gStyle->SetPalette(1,0);
11 gStyle->SetOptStat(kTRUE);
12 gStyle->SetOptStat(11111111);
13 gStyle->SetOptFit(kTRUE);
14 gStyle->SetPadRightMargin(0.15);
15 gStyle->SetPadLeftMargin(0.15);
16 gStyle->SetPadBottomMargin(0.15);
26 Double_t g2931[nchan];
27 Double_t g3079[nchan];
28 Double_t g3179[nchan];
29 Double_t g3180[nchan];
30 Double_t g3182[nchan];
31 Double_t g2931_r20_60[nchan];
32 Double_t g3079_r20_60[nchan];
33 Double_t g3179_r20_60[nchan];
34 Double_t g3180_r20_60[nchan];
35 Double_t g3182_r20_60[nchan];
37 Double_t g3079_3182_r20[nchan];
38 Double_t g3079_3182_r20_60[nchan];
39 Double_t g2931_3182_r20[nchan];
40 Double_t g2931_3182_r20_60[nchan];
44 sprintf(
string,
"GainFactorIter1_2931-3182_r20-60.root");
45 printf (
"Histogram input filename=%s\n",
string);
46 TFile *in =
new TFile(
string,
"read");
48 TH1D *h_unordered2931_3182_r20_60 = (TH1D*)in->Get(
"h_unordered");
49 TH2D *h_yvsx2931_3182_r20_60 = (TH2D*)in->Get(
"h_yvsx");
53 outgain.open (
"plot_gain_factors2.list");
59 TString
infile =
"GainFactorIter1_2931-3182_r20-60.txt";
60 cout <<
"Opening file: " << infile.Data() << endl;
61 gainfile.open (infile.Data());
63 cout <<
"ERROR: Failed to open data file= " << infile.Data() << endl;
68 while (line.ReadLine(gainfile)){
70 TObjArray *tokens = line.Tokenize(
" ");
71 Float_t gain = (((TObjString*)tokens->At(0))->GetString()).Atof();
72 TString s = ((TObjString*)tokens->At(1))->GetString().Remove(0,8);
73 Int_t channel = s.Atoi();
74 g2931_3182_r20_60[channel] = gain;
75 outgain <<
" R3182_r20_60, channel=" << channel <<
" gain=" << gain << endl;
82 TCanvas *c10 =
new TCanvas(
"c10",
"c10 plot_gain_factors2",200,10,1000,500);
93 sprintf(
string,
"Electron Gain Factors Run 2931-3182, Jul27");
94 h_yvsx2931_3182_r20_60->SetTitle(
string);
95 h_yvsx2931_3182_r20_60->GetZaxis()->SetRangeUser(zmin,zmax);
96 h_yvsx2931_3182_r20_60->GetYaxis()->SetTitle(
"y");
97 h_yvsx2931_3182_r20_60->GetXaxis()->SetTitle(
"x");
98 h_yvsx2931_3182_r20_60->Draw(
"colz");
104 h_unordered2931_3182_r20_60->SetTitle(
string);;
105 h_unordered2931_3182_r20_60->GetYaxis()->SetRangeUser(ymin,ymax);
106 h_unordered2931_3182_r20_60->GetYaxis()->SetTitle(
"Counts");
107 h_unordered2931_3182_r20_60->GetXaxis()->SetTitle(
"Gain Factor");
108 h_unordered2931_3182_r20_60->Fit(
"gaus");
109 h_unordered2931_3182_r20_60->Draw(
"");
121 sprintf(
string,
"plot_gain_factors2.pdf");
sprintf(text,"Post KinFit Cut")
void plot_gain_factors2(void)
printf("string=%s", string)