4{
5 Int_t runno;
6 Double_t runmean, rungain, runresol;
7
8 TFile fdedx("DedxConst_test.root","recreate");
9 TTree* runbyrun = new TTree("runcalib", "runcalib");
10 runbyrun ->
Branch(
"runno", &runno,
"runno/I");
11 runbyrun ->
Branch(
"runmean", &runmean,
"runmean/D");
12 runbyrun ->
Branch(
"rungain", &rungain,
"rungain/D");
13 runbyrun ->
Branch(
"runresol", &runresol,
"runresol/D");
14
15 TFile fconst0("rungain.root");
16 TTree *tree_rung = (TTree*)fconst0.Get("runcalib");
21
22 long NRUN = tree_rung -> GetEntries();
23 cout<<"Run Number == "<<NRUN<<endl;
24 for(long runid =0; runid<NRUN; runid++)
25 {
26 tree_rung->GetEntry(runid);
27 runbyrun ->Fill();
28 }
29
30 TTree *tree_resg = (TTree*)fconst0.Get("resolcalib");
31 Double_t resolpar[1];
34
35 TTree *tree_gaing = (TTree*)fconst0.Get("gaincalib");
36 Double_t gainpar[1];
39
40 fconst0.Close();
41 cout<<"runbyrun const OK!!!"<<endl;
42
43
45
46
47
48 Double_t gwire[6796];
49 TTree *tree_wireg = (TTree*)fconst1.Get("wiregcalib");
52 fconst1.Close();
53 cout<<"wiregcalib const OK!!!"<<endl;
54
55
57 Double_t ggs0[43], ggs1[43], ggs2[43], ggs3[43];
58 TTree *tree_ggs = (TTree*)fconst2.Get("ggscalib");
63 tree_ggs->GetEntry(0);
64 fconst2.Close();
65 cout<<"ggscalib const OK!!!"<<endl;
66
67
69 Double_t ddg0[43], ddg1[43], ddg2[43], ddg3[43];
70 TTree *tree_ddg = (TTree*)fconst3.Get("ddgcalib");
76 fconst3.Close();
77 cout<<"ddgcalib const OK!!!"<<endl;
78
79
81 Double_t entra0[43], entra1[43], entra2[43], entra3[43];
82 TTree *tree_entra = (TTree*)fconst4.Get("entracalib");
87
88
89
90
91 tree_entra->GetEntry(0);
92 fconst4.Close();
93
94
96 int entry[1]; double denangle[100];
97 TTree *tree_entra1 = (TTree*)fconst4_1.Get("entracalib");
100 tree_entra1 ->GetEntry(0);
101 fconst4_1.Close();
102 cout<<"entrance angle const OK!!!"<<endl;
103
104
106 Double_t zdep0[43], zdep1[43], zdep2[43], zdep3[43];
107 TTree *tree_zdep = (TTree*)fconst5.Get("zdepcalib");
112 tree_zdep->GetEntry(0);
113 fconst5.Close();
114 cout<<"zdepcalib const OK!!!"<<endl;
115
116
118 Double_t glayer[43];
119 TTree *tree_layerg = (TTree*)fconst6.Get("layergcalib");
121 tree_layerg->GetEntry(0);
122 fconst6.Close();
123 cout<<"layergcalib const OK!!!"<<endl;
124
125
126
127
129 TTree *tree_costheta = (TTree*)fconst7.Get("ggscalib");
130 double gcostheta[80];
131 tree_costheta ->SetBranchAddress("costheta",gcostheta);
132 tree_costheta->GetEntry(0);
133 fconst7.Close();
134 cout<<"costhetacalib const OK!!!"<<endl;
135
136
138 TTree *tree_t0=(TTree*)fconst8.Get("gaincalib");
139
140 double gt0[35];
141 double gdedx[35];
142 tree_t0->SetBranchAddress("t0",gt0);
143 tree_t0->SetBranchAddress("dedx",gdedx);
144 for(int i=0; i<tree_t0->GetEntries(); i++)
145 {
146 tree_t0->GetEntry(i);
147
148
149 }
150 fconst8.Close();
151 cout<<"t0calib const OK!!!"<<endl;
152
153
154
157 double Iner_gain[
n], Iner_chi[
n], Iner_hits[
n];
158 double Out_gain[
n], Out_chi[
n], Out_hits[
n];
159 double Id_doca[
n], Ip_eangle[
n];
160
161 TTree *tree_docasin = (TTree*)fconst9.Get("ddgcalib");
171 fconst9.Close();
172
173 cout<<"ddgcalib const OK!!!"<<endl;
174
176
177
178 int hadronNo;
179 double ghadron[5];
180 TTree *tree_hadron = (TTree*)fconst10.Get("ggscalib");
183 tree_hadron->GetEntry(0);
184 fconst10.Close();
185 cout<<"hadron const OK!!!"<<endl;
186
187
188 fdedx.cd();
189 TTree* gain = new TTree("gaincalib", "gaincalib");
190 gain ->
Branch(
"gain", gainpar,
"gain[1]/D");
191 gain ->
Branch(
"t0", gt0,
"t0[35]/D");
192 gain ->
Branch(
"dedx", gdedx,
"dedx[35]/D");
193
194 TTree* resol = new TTree("resolcalib", "resolcalib");
195 resol ->
Branch(
"resol", resolpar,
"resol[1]/D");
196
197 TTree* wireg = new TTree("wiregcalib", "wiregcalib");
198 wireg ->
Branch(
"wireg", gwire,
"wireg[6796]/D");
199 wireg-> SetBasketSize("wireg",2000000);
200
201 TTree* ggs = new TTree("ggscalib", "ggscalib");
202 ggs ->
Branch(
"ggs0", ggs0,
"ggs[43]/D");
203 ggs ->
Branch(
"ggs1", ggs1,
"ggs[43]/D");
204 ggs ->
Branch(
"ggs2", ggs2,
"ggs[43]/D");
205 ggs ->
Branch(
"ggs3", ggs3,
"ggs[43]/D");
206 ggs ->
Branch(
"costheta", gcostheta,
"costheta[80]/D");
207 ggs ->
Branch(
"hadron",ghadron,
"hadron[5]/D");
208 ggs ->
Branch(
"hadronNo",&hadronNo,
"hadronNo/I");
209
210 TTree* ddg = new TTree("ddgcalib", "ddgcalib");
211 ddg ->
Branch(
"ddg0", ddg0,
"ddg[43]/D");
212 ddg ->
Branch(
"ddg1", ddg1,
"ddg[43]/D");
213 ddg ->
Branch(
"ddg2", ddg2,
"ddg[43]/D");
214 ddg ->
Branch(
"ddg3", ddg3,
"ddg[43]/D");
215 ddg ->
Branch(
"Iner_gain", Iner_gain,
"Iner_gain[1600]/D");
216 ddg ->
Branch(
"Out_gain", Out_gain,
"Out_gain[1600]/D");
217 ddg ->
Branch(
"Iner_chi", Iner_chi,
"Iner_chi[1600]/D");
218 ddg ->
Branch(
"Out_chi", Out_chi,
"Out_chi[1600]/D");
219 ddg ->
Branch(
"Iner_hits", Iner_hits,
"Iner_hits[1600]/D");
220 ddg ->
Branch(
"Out_hits", Out_hits,
"Out_hits[1600]/D");
221 ddg ->
Branch(
"Id_doca", Id_doca,
"Id_doca[1600]/D");
222 ddg ->
Branch(
"Ip_eangle", Ip_eangle,
"Ip_eangle[1600]/D");
223
224 TTree* entra = new TTree("entracalib", "entracalib");
225 entra ->
Branch(
"entra0", entra0,
"entra[43]/D");
226 entra ->
Branch(
"entra1", entra1,
"entra[43]/D");
227 entra ->
Branch(
"entra2", entra2,
"entra[43]/D");
228 entra ->
Branch(
"entra3", entra3,
"entra[43]/D");
229 entra ->
Branch(
"1denangle_entry", entry,
"1denangle_entry[1]/I");
230 entra ->
Branch(
"1denangle", denangle,
"1denangle[100]/D");
231
232 TTree* zdep = new TTree("zdepcalib", "zdepcalib");
233 zdep ->
Branch(
"zdep0", zdep0,
"zdep[43]/D");
234 zdep ->
Branch(
"zdep1", zdep1,
"zdep[43]/D");
235 zdep ->
Branch(
"zdep2", zdep2,
"zdep[43]/D");
236 zdep ->
Branch(
"zdep3", zdep3,
"zdep[43]/D");
237
238 TTree* layerg = new TTree("layergcalib", "layergcalib");
239 layerg ->
Branch(
"layerg", glayer,
"layerg[43]/D");
240
241
242 gain->Fill();
243 resol->Fill();
250
251
252
253
254
255
256
257
258
259
260 fdedx.Write();
261 fdedx.Close();
262 cout<<"All const OK!!!"<<endl;
263}
curve Branch("CurveSize",&CurveSize,"CurveSize/I")
data SetBranchAddress("time",&time)