1TString
old_str =
"../Eangle1D/DedxConst_test.root";
5 Int_t runno, evtfrom, evtto;
6 Double_t runmean, rungain, runresol;
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 runbyrun ->
Branch(
"evtfrom", &evtfrom,
"evtfrom/I");
15 runbyrun ->
Branch(
"evtto", &evtto,
"evtto/I");
18 TTree *tree_rung = (TTree*)fconst0.Get(
"runcalib");
26 long NRUN = tree_rung -> GetEntries();
27 cout<<
"Run Number == "<<NRUN<<endl;
28 for(
long runid =0; runid<NRUN; runid++)
30 tree_rung->GetEntry(runid);
34 TTree *tree_resg = (TTree*)fconst0.Get(
"resolcalib");
39 TTree *tree_gaing = (TTree*)fconst0.Get(
"gaincalib");
45 cout<<
"runbyrun const OK!!!"<<endl;
53 TTree *tree_wireg = (TTree*)fconst1.Get(
"wiregcalib");
57 cout<<
"wiregcalib const OK!!!"<<endl;
61 Double_t ggs0[43], ggs1[43], ggs2[43], ggs3[43];
62 TTree *tree_ggs = (TTree*)fconst2.Get(
"ggscalib");
67 tree_ggs->GetEntry(0);
69 cout<<
"ggscalib const OK!!!"<<endl;
73 Double_t ddg0[43], ddg1[43], ddg2[43], ddg3[43];
74 TTree *tree_ddg = (TTree*)fconst3.Get(
"ddgcalib");
81 cout<<
"ddgcalib const OK!!!"<<endl;
85 Double_t entra0[43], entra1[43], entra2[43], entra3[43];
86 TTree *tree_entra = (TTree*)fconst4.Get(
"entracalib");
95 tree_entra->GetEntry(0);
101 int entry[1];
double denangle[100];
102 TTree *tree_entra1 = (TTree*)fconst4_1.Get(
"entracalib");
105 tree_entra1 ->GetEntry(0);
107 cout<<
"entrance angle const OK!!!"<<endl;
111 Double_t zdep0[43], zdep1[43], zdep2[43], zdep3[43];
112 TTree *tree_zdep = (TTree*)fconst5.Get(
"zdepcalib");
117 tree_zdep->GetEntry(0);
119 cout<<
"zdepcalib const OK!!!"<<endl;
124 TTree *tree_layerg = (TTree*)fconst6.Get(
"layergcalib");
126 tree_layerg->GetEntry(0);
128 cout<<
"layergcalib const OK!!!"<<endl;
133 TTree *tree_costheta = (TTree*)fconst7.Get(
"ggscalib");
135 double gcostheta[80];
136 tree_costheta ->SetBranchAddress(
"costheta",gcostheta);
137 tree_costheta->GetEntry(0);
139 cout<<
"costhetacalib const OK!!!"<<endl;
142 TFile fconst8(
"t0.root");
144 TTree *tree_t0=(TTree*)fconst8.Get(
"gaincalib");
147 tree_t0->SetBranchAddress(
"t0",gt0);
148 tree_t0->SetBranchAddress(
"dedx",gdedx);
149 for(
int i=0; i<tree_t0->GetEntries(); i++) tree_t0->GetEntry(i);
151 cout<<
"t0calib const OK!!!"<<endl;
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];
161 TTree *tree_docasin = (TTree*)fconst9.Get(
"ddgcalib");
173 cout<<
"ddgcalib const OK!!!"<<endl;
180 TTree *tree_hadron = (TTree*)fconst10.Get(
"ggscalib");
183 tree_hadron->GetEntry(0);
185 cout<<
"hadron const OK!!!"<<endl;
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");
194 TTree* resol =
new TTree(
"resolcalib",
"resolcalib");
195 resol ->
Branch(
"resol", resolpar,
"resol[1]/D");
197 TTree* wireg =
new TTree(
"wiregcalib",
"wiregcalib");
198 wireg ->
Branch(
"wireg", gwire,
"wireg[6796]/D");
199 wireg-> SetBasketSize(
"wireg",2000000);
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");
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");
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");
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");
238 TTree* layerg =
new TTree(
"layergcalib",
"layergcalib");
239 layerg ->
Branch(
"layerg", glayer,
"layerg[43]/D");
262 cout<<
"All const OK!!!"<<endl;
curve Branch("CurveSize",&CurveSize,"CurveSize/I")
data SetBranchAddress("time",&time)