1TString
old_str =
"../RungEangle/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");
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;
51 TTree *tree_wireg = (TTree*)fconst1.Get(
"wiregcalib");
55 cout<<
"wiregcalib const OK!!!"<<endl;
59 Double_t ggs0[43], ggs1[43], ggs2[43], ggs3[43];
60 TTree *tree_ggs = (TTree*)fconst2.Get(
"ggscalib");
65 tree_ggs->GetEntry(0);
67 cout<<
"ggscalib const OK!!!"<<endl;
71 Double_t ddg0[43], ddg1[43], ddg2[43], ddg3[43];
72 TTree *tree_ddg = (TTree*)fconst3.Get(
"ddgcalib");
79 cout<<
"ddgcalib const OK!!!"<<endl;
83 Double_t entra0[43], entra1[43], entra2[43], entra3[43];
84 TTree *tree_entra = (TTree*)fconst4.Get(
"entracalib");
93 tree_entra->GetEntry(0);
98 int entry[1];
double denangle[100];
99 TTree *tree_entra1 = (TTree*)fconst4_1.Get(
"entracalib");
102 tree_entra1 ->GetEntry(0);
104 cout<<
"entrance angle const OK!!!"<<endl;
108 Double_t zdep0[43], zdep1[43], zdep2[43], zdep3[43];
109 TTree *tree_zdep = (TTree*)fconst5.Get(
"zdepcalib");
114 tree_zdep->GetEntry(0);
116 cout<<
"zdepcalib const OK!!!"<<endl;
121 TTree *tree_layerg = (TTree*)fconst6.Get(
"layergcalib");
123 tree_layerg->GetEntry(0);
125 cout<<
"layergcalib const OK!!!"<<endl;
131 TTree *tree_costheta = (TTree*)fconst7.Get(
"ggscalib");
132 double gcostheta[80];
133 tree_costheta ->SetBranchAddress(
"costheta",gcostheta);
134 tree_costheta->GetEntry(0);
136 cout<<
"costhetacalib const OK!!!"<<endl;
140 TTree *tree_t0=(TTree*)fconst8.Get(
"gaincalib");
143 tree_t0->SetBranchAddress(
"t0", gt0);
144 tree_t0->SetBranchAddress(
"dedx",gdedx);
145 for(
int i=0; i<tree_t0->GetEntries(); i++) tree_t0->GetEntry(i);
147 cout<<
"t0calib const OK!!!"<<endl;
153 double Iner_gain[
n], Iner_chi[
n], Iner_hits[
n];
154 double Out_gain[
n], Out_chi[
n], Out_hits[
n];
155 double Id_doca[
n], Ip_eangle[
n];
157 TTree *tree_docasin = (TTree*)fconst9.Get(
"ddgcalib");
169 cout<<
"ddgcalib const OK!!!"<<endl;
171 TFile fconst10(
"hadron.root");
176 TTree *tree_hadron = (TTree*)fconst10.Get(
"ggscalib");
179 tree_hadron->GetEntry(0);
181 cout<<
"hadron const OK!!!"<<endl;
185 TTree* gain =
new TTree(
"gaincalib",
"gaincalib");
186 gain ->
Branch(
"gain", gainpar,
"gain[1]/D");
187 gain ->
Branch(
"t0", gt0,
"t0[35]/D");
188 gain ->
Branch(
"dedx", gdedx,
"dedx[35]/D");
190 TTree* resol =
new TTree(
"resolcalib",
"resolcalib");
191 resol ->
Branch(
"resol", resolpar,
"resol[1]/D");
193 TTree* wireg =
new TTree(
"wiregcalib",
"wiregcalib");
194 wireg ->
Branch(
"wireg", gwire,
"wireg[6796]/D");
195 wireg-> SetBasketSize(
"wireg",2000000);
197 TTree* ggs =
new TTree(
"ggscalib",
"ggscalib");
198 ggs ->
Branch(
"ggs0", ggs0,
"ggs[43]/D");
199 ggs ->
Branch(
"ggs1", ggs1,
"ggs[43]/D");
200 ggs ->
Branch(
"ggs2", ggs2,
"ggs[43]/D");
201 ggs ->
Branch(
"ggs3", ggs3,
"ggs[43]/D");
202 ggs ->
Branch(
"costheta", gcostheta,
"costheta[80]/D");
203 ggs ->
Branch(
"hadron",ghadron,
"hadron[5]/D");
204 ggs ->
Branch(
"hadronNo",&hadronNo,
"hadronNo/I");
206 TTree* ddg =
new TTree(
"ddgcalib",
"ddgcalib");
207 ddg ->
Branch(
"ddg0", ddg0,
"ddg[43]/D");
208 ddg ->
Branch(
"ddg1", ddg1,
"ddg[43]/D");
209 ddg ->
Branch(
"ddg2", ddg2,
"ddg[43]/D");
210 ddg ->
Branch(
"ddg3", ddg3,
"ddg[43]/D");
211 ddg ->
Branch(
"Iner_gain", Iner_gain,
"Iner_gain[1600]/D");
212 ddg ->
Branch(
"Out_gain", Out_gain,
"Out_gain[1600]/D");
213 ddg ->
Branch(
"Iner_chi", Iner_chi,
"Iner_chi[1600]/D");
214 ddg ->
Branch(
"Out_chi", Out_chi,
"Out_chi[1600]/D");
215 ddg ->
Branch(
"Iner_hits", Iner_hits,
"Iner_hits[1600]/D");
216 ddg ->
Branch(
"Out_hits", Out_hits,
"Out_hits[1600]/D");
217 ddg ->
Branch(
"Id_doca", Id_doca,
"Id_doca[1600]/D");
218 ddg ->
Branch(
"Ip_eangle", Ip_eangle,
"Ip_eangle[1600]/D");
220 TTree* entra =
new TTree(
"entracalib",
"entracalib");
221 entra ->
Branch(
"entra0", entra0,
"entra[43]/D");
222 entra ->
Branch(
"entra1", entra1,
"entra[43]/D");
223 entra ->
Branch(
"entra2", entra2,
"entra[43]/D");
224 entra ->
Branch(
"entra3", entra3,
"entra[43]/D");
225 entra ->
Branch(
"1denangle_entry", entry,
"1denangle_entry[1]/I");
226 entra ->
Branch(
"1denangle", denangle,
"1denangle[100]/D");
228 TTree* zdep =
new TTree(
"zdepcalib",
"zdepcalib");
229 zdep ->
Branch(
"zdep0", zdep0,
"zdep[43]/D");
230 zdep ->
Branch(
"zdep1", zdep1,
"zdep[43]/D");
231 zdep ->
Branch(
"zdep2", zdep2,
"zdep[43]/D");
232 zdep ->
Branch(
"zdep3", zdep3,
"zdep[43]/D");
234 TTree* layerg =
new TTree(
"layergcalib",
"layergcalib");
235 layerg ->
Branch(
"layerg", glayer,
"layerg[43]/D");
258 cout<<
"All const OK!!!"<<endl;
curve Branch("CurveSize",&CurveSize,"CurveSize/I")
data SetBranchAddress("time",&time)