1TString
old_str =
"../Eangle1DSec/DedxConst_test.root";
6 Int_t runno, evtfrom, evtto;
7 Double_t runmean, rungain, runresol;
9 TFile fdedx(
"DedxConst_test.root",
"recreate");
10 TTree* runbyrun =
new TTree(
"runcalib",
"runcalib");
11 runbyrun ->
Branch(
"runno", &runno,
"runno/I");
12 runbyrun ->
Branch(
"runmean", &runmean,
"runmean/D");
13 runbyrun ->
Branch(
"rungain", &rungain,
"rungain/D");
14 runbyrun ->
Branch(
"runresol", &runresol,
"runresol/D");
15 runbyrun ->
Branch(
"evtfrom", &evtfrom,
"evtfrom/I");
16 runbyrun ->
Branch(
"evtto", &evtto,
"evtto/I");
18 TFile fconst0(
"rungain.root");
19 TTree *tree_rung = (TTree*)fconst0.Get(
"runcalib");
27 long NRUN = tree_rung -> GetEntries();
28 cout<<
"Run Number == "<<NRUN<<endl;
29 for(
long runid =0; runid<NRUN; runid++)
31 tree_rung->GetEntry(runid);
35 TTree *tree_resg = (TTree*)fconst0.Get(
"resolcalib");
40 TTree *tree_gaing = (TTree*)fconst0.Get(
"gaincalib");
46 cout<<
"runbyrun const OK!!!"<<endl;
54 TTree *tree_wireg = (TTree*)fconst1.Get(
"wiregcalib");
58 cout<<
"wiregcalib const OK!!!"<<endl;
62 Double_t ggs0[43], ggs1[43], ggs2[43], ggs3[43];
63 TTree *tree_ggs = (TTree*)fconst2.Get(
"ggscalib");
68 tree_ggs->GetEntry(0);
70 cout<<
"ggscalib const OK!!!"<<endl;
74 Double_t ddg0[43], ddg1[43], ddg2[43], ddg3[43];
75 TTree *tree_ddg = (TTree*)fconst3.Get(
"ddgcalib");
82 cout<<
"ddgcalib const OK!!!"<<endl;
86 Double_t entra0[43], entra1[43], entra2[43], entra3[43];
87 TTree *tree_entra = (TTree*)fconst4.Get(
"entracalib");
96 tree_entra->GetEntry(0);
102 int entry[1];
double denangle[100];
103 TTree *tree_entra1 = (TTree*)fconst4_1.Get(
"entracalib");
106 tree_entra1 ->GetEntry(0);
108 cout<<
"entrance angle const OK!!!"<<endl;
112 Double_t zdep0[43], zdep1[43], zdep2[43], zdep3[43];
113 TTree *tree_zdep = (TTree*)fconst5.Get(
"zdepcalib");
118 tree_zdep->GetEntry(0);
120 cout<<
"zdepcalib const OK!!!"<<endl;
125 TTree *tree_layerg = (TTree*)fconst6.Get(
"layergcalib");
127 tree_layerg->GetEntry(0);
129 cout<<
"layergcalib const OK!!!"<<endl;
135 TTree *tree_costheta = (TTree*)fconst7.Get(
"ggscalib");
136 double gcostheta[80];
137 tree_costheta ->SetBranchAddress(
"costheta",gcostheta);
138 tree_costheta->GetEntry(0);
140 cout<<
"costhetacalib const OK!!!"<<endl;
145 TTree *tree_t0=(TTree*)fconst8.Get(
"gaincalib");
149 tree_t0->SetBranchAddress(
"t0",gt0);
150 tree_t0->SetBranchAddress(
"dedx",gdedx);
151 for(
int i=0; i<tree_t0->GetEntries(); i++)
153 tree_t0->GetEntry(i);
158 cout<<
"t0calib const OK!!!"<<endl;
164 double Iner_gain[
n], Iner_chi[
n], Iner_hits[
n];
165 double Out_gain[
n], Out_chi[
n], Out_hits[
n];
166 double Id_doca[
n], Ip_eangle[
n];
168 TTree *tree_docasin = (TTree*)fconst9.Get(
"ddgcalib");
180 cout<<
"ddgcalib const OK!!!"<<endl;
187 TTree *tree_hadron = (TTree*)fconst10.Get(
"ggscalib");
190 tree_hadron->GetEntry(0);
192 cout<<
"hadron const OK!!!"<<endl;
196 TTree* gain =
new TTree(
"gaincalib",
"gaincalib");
197 gain ->
Branch(
"gain", gainpar,
"gain[1]/D");
198 gain ->
Branch(
"t0", gt0,
"t0[35]/D");
199 gain ->
Branch(
"dedx", gdedx,
"dedx[35]/D");
201 TTree* resol =
new TTree(
"resolcalib",
"resolcalib");
202 resol ->
Branch(
"resol", resolpar,
"resol[1]/D");
204 TTree* wireg =
new TTree(
"wiregcalib",
"wiregcalib");
205 wireg ->
Branch(
"wireg", gwire,
"wireg[6796]/D");
206 wireg-> SetBasketSize(
"wireg",2000000);
208 TTree* ggs =
new TTree(
"ggscalib",
"ggscalib");
209 ggs ->
Branch(
"ggs0", ggs0,
"ggs[43]/D");
210 ggs ->
Branch(
"ggs1", ggs1,
"ggs[43]/D");
211 ggs ->
Branch(
"ggs2", ggs2,
"ggs[43]/D");
212 ggs ->
Branch(
"ggs3", ggs3,
"ggs[43]/D");
213 ggs ->
Branch(
"costheta", gcostheta,
"costheta[80]/D");
214 ggs ->
Branch(
"hadron",ghadron,
"hadron[5]/D");
215 ggs ->
Branch(
"hadronNo",&hadronNo,
"hadronNo/I");
217 TTree* ddg =
new TTree(
"ddgcalib",
"ddgcalib");
218 ddg ->
Branch(
"ddg0", ddg0,
"ddg[43]/D");
219 ddg ->
Branch(
"ddg1", ddg1,
"ddg[43]/D");
220 ddg ->
Branch(
"ddg2", ddg2,
"ddg[43]/D");
221 ddg ->
Branch(
"ddg3", ddg3,
"ddg[43]/D");
222 ddg ->
Branch(
"Iner_gain", Iner_gain,
"Iner_gain[1600]/D");
223 ddg ->
Branch(
"Out_gain", Out_gain,
"Out_gain[1600]/D");
224 ddg ->
Branch(
"Iner_chi", Iner_chi,
"Iner_chi[1600]/D");
225 ddg ->
Branch(
"Out_chi", Out_chi,
"Out_chi[1600]/D");
226 ddg ->
Branch(
"Iner_hits", Iner_hits,
"Iner_hits[1600]/D");
227 ddg ->
Branch(
"Out_hits", Out_hits,
"Out_hits[1600]/D");
228 ddg ->
Branch(
"Id_doca", Id_doca,
"Id_doca[1600]/D");
229 ddg ->
Branch(
"Ip_eangle", Ip_eangle,
"Ip_eangle[1600]/D");
231 TTree* entra =
new TTree(
"entracalib",
"entracalib");
232 entra ->
Branch(
"entra0", entra0,
"entra[43]/D");
233 entra ->
Branch(
"entra1", entra1,
"entra[43]/D");
234 entra ->
Branch(
"entra2", entra2,
"entra[43]/D");
235 entra ->
Branch(
"entra3", entra3,
"entra[43]/D");
236 entra ->
Branch(
"1denangle_entry", entry,
"1denangle_entry[1]/I");
237 entra ->
Branch(
"1denangle", denangle,
"1denangle[100]/D");
239 TTree* zdep =
new TTree(
"zdepcalib",
"zdepcalib");
240 zdep ->
Branch(
"zdep0", zdep0,
"zdep[43]/D");
241 zdep ->
Branch(
"zdep1", zdep1,
"zdep[43]/D");
242 zdep ->
Branch(
"zdep2", zdep2,
"zdep[43]/D");
243 zdep ->
Branch(
"zdep3", zdep3,
"zdep[43]/D");
245 TTree* layerg =
new TTree(
"layergcalib",
"layergcalib");
246 layerg ->
Branch(
"layerg", glayer,
"layerg[43]/D");
269 cout<<
"All const OK!!!"<<endl;
curve Branch("CurveSize",&CurveSize,"CurveSize/I")
data SetBranchAddress("time",&time)