Analysis Software
Documentation for sPHENIX simulation software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
DrawEcal_pDST.C
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file DrawEcal_pDST.C
1 // $Id: $
2 
11 #include <cmath>
12 #include <TFile.h>
13 #include <TString.h>
14 #include <TLine.h>
15 #include <TTree.h>
16 #include <cassert>
17 #include "SaveCanvas.C"
18 #include "SetOKStyle.C"
19 using namespace std;
20 
21 TFile * _file0 = NULL;
22 TTree * T = NULL;
23 TString cuts = "";
24 
25 void
26 DrawEcal_pDST( //
27  const TString infile = "G4Hits_sPHENIX_gamma_eta0_8GeV-0000.root_Ana.root"//
28  )
29 {
30  SetOKStyle();
31  gStyle->SetOptStat(0);
32  gStyle->SetOptFit(1111);
33  TVirtualFitter::SetDefaultFitter("Minuit2");
34 
35  gSystem->Load("libg4eval.so");
36  gSystem->Load("libemcal_ana.so");
37  gSystem->Load("libg4vertex.so");
38 
39  if (!_file0)
40  {
41  TString chian_str = infile;
42  chian_str.ReplaceAll("ALL", "*");
43  chian_str.ReplaceAll("+", "\\+");
44 
45  TChain * t = new TChain("T");
46  const int n = t->Add(chian_str);
47 
48  cout << "Loaded " << n << " root files with " << chian_str << endl;
49  assert(n>0);
50 
51  T = t;
52 
53  _file0 = new TFile;
54  _file0->SetName(infile);
55 
56  fstream flst(infile + ".lst", ios_base::out);
57 
58  TObjArray *fileElements = t->GetListOfFiles();
59  TIter next(fileElements);
60  TChainElement *chEl = 0;
61  while ((chEl = (TChainElement*) next()))
62  {
63  flst << chEl->GetTitle() << endl;
64  }
65  flst.close();
66 
67  cout << "Saved file list to " << infile + ".lst" << endl;
68  }
69 
70  assert(_file0);
71 
72  T->SetAlias("UpsilonPair_trk_gpt",
73  "1*sqrt(DST.UpsilonPair.trk.gpx**2 + DST.UpsilonPair.trk.gpy**2)");
74  T->SetAlias("UpsilonPair_trk_pt",
75  "1*sqrt(DST.UpsilonPair.trk.px**2 + DST.UpsilonPair.trk.py**2)");
76 
77  T->SetAlias("MCPhoton_pt", "1*sqrt(DST.MCPhoton.px**2 + DST.MCPhoton.py**2)");
78  T->SetAlias("MCPhoton_gpt",
79  "1*sqrt(DST.MCPhoton.gpx**2 + DST.MCPhoton.gpy**2)");
80 
81  const TCut event_sel = "1*1";
82  cuts = "_all_event";
83 
84  cout << "Build event selection of " << (const char *) event_sel << endl;
85 
86  T->Draw(">>EventList", event_sel);
87  TEventList * elist = gDirectory->GetObjectChecked("EventList", "TEventList");
88  cout << elist->GetN() << " / " << T->GetEntriesFast() << " events selected"
89  << endl;
90 
91  T->SetEventList(elist);
93 // UpsilonPair_Checks(infile);
94 
95  TrackProjection_Checks(infile, " ");
96 // Edep_Checks(infile, 1.4, 1.4, " ");
97 
98 }
99 
100 void
101 MakeRadiusCut(TString infile, const double R_CEMC, const double R_HCALIN,
102  TCut good_track_cut)
103 {
104  TCanvas *c1 = new TCanvas("MakeRadiusCut" + cuts, "MakeRadiusCut" + cuts,
105  1800, 900);
106  c1->Divide(4, 2);
107  int idx = 1;
108  TPad * p;
109 
110  p = (TPad *) c1->cd(idx++);
111  c1->Update();
112  p->SetLogy();
113 
114  T->Draw("DST.MCPhoton.cemc_iphi>>hcemc_iphi(130,-6.5,6.5)",
115  TString("DST.MCPhoton.cemc_energy * (")
116  + TString(good_track_cut.GetTitle()) + ")");
117  hcemc_iphi->SetTitle(
118  ";CEMC Azimuthal Distance (Tower Width);Energy Distribution");
119  hcemc_iphi->Fit("gaus", "M");
120  TF1* f = (TF1*) (hcemc_iphi->GetListOfFunctions()->At(0));
121  assert(f);
122  T->SetAlias("MCPhoton_cor_cemc_iphi",
123  Form("DST.MCPhoton.cemc_iphi - %f", f->GetParameter(1)));
124 
125  const double center_cemc_iphi = f->GetParameter(1);
126 
127  p = (TPad *) c1->cd(idx++);
128  c1->Update();
129  p->SetLogy();
130 
131  T->Draw("abs(MCPhoton_cor_cemc_iphi)>>hMCPhoton_cor_cemc_iphi(130,0,6.5)",
132  TString("DST.MCPhoton.cemc_energy * (")
133  + TString(good_track_cut.GetTitle()) + ")");
134 
135  p->DrawFrame(0, 1e-3, 6.5, 1,
136  ";CEMC Azimuthal Distance (Tower Width);Energy Leakage Ratio");
137  TGraphErrors * ge = Distribution2Efficiency(hMCPhoton_cor_cemc_iphi);
138  ge->SetLineColor(kBlue + 2);
139  ge->SetMarkerColor(kBlue + 21);
140  ge->SetMarkerColor(kFullCircle);
141  ge->SetLineWidth(3);
142  ge->Draw("lp");
143 
144  p = (TPad *) c1->cd(idx++);
145  c1->Update();
146  p->SetLogy();
147 
148  T->Draw("DST.MCPhoton.cemc_ieta>>hcemc_ieta(130,-6.5,6.5)",
149  TString("DST.MCPhoton.cemc_energy * (")
150  + TString(good_track_cut.GetTitle()) + ")");
151  hcemc_ieta->SetTitle(
152  ";CEMC Polar Distance (Tower Width);Energy Distribution");
153  hcemc_ieta->Fit("gaus", "M");
154  TF1* f = (TF1*) (hcemc_ieta->GetListOfFunctions()->At(0));
155  assert(f);
156  T->SetAlias("MCPhoton_cor_cemc_ieta",
157  Form("DST.MCPhoton.cemc_ieta - %f", f->GetParameter(1)));
158 
159  const double center_cemc_ieta = f->GetParameter(1);
160 
161  p = (TPad *) c1->cd(idx++);
162  c1->Update();
163  p->SetLogy();
164 
165  T->Draw("abs(MCPhoton_cor_cemc_ieta)>>hMCPhoton_cor_cemc_ieta(130,0,6.5)",
166  TString("DST.MCPhoton.cemc_energy * (")
167  + TString(good_track_cut.GetTitle()) + ")");
168 
169  p->DrawFrame(0, 1e-3, 6.5, 1,
170  ";CEMC Polar Distance (Tower Width);Energy Leakage Ratio");
171  TGraphErrors * ge = Distribution2Efficiency(hMCPhoton_cor_cemc_ieta);
172  ge->SetLineColor(kBlue + 2);
173  ge->SetMarkerColor(kBlue + 21);
174  ge->SetMarkerColor(kFullCircle);
175  ge->SetLineWidth(3);
176  ge->Draw("lp");
177 
178  p = (TPad *) c1->cd(idx++);
179  c1->Update();
180  p->SetLogy();
181 
182  T->Draw("DST.MCPhoton.hcalin_iphi>>hhcalin_iphi(130,-6.5,6.5)",
183  TString("DST.MCPhoton.hcalin_energy * (")
184  + TString(good_track_cut.GetTitle()) + ")");
185  hhcalin_iphi->SetTitle(
186  ";HCal_{IN} Azimuthal Distance (Tower Width);Energy Distribution");
187  hhcalin_iphi->Fit("gaus", "M");
188  TF1* f = (TF1*) (hhcalin_iphi->GetListOfFunctions()->At(0));
189  assert(f);
190  T->SetAlias("MCPhoton_cor_hcalin_iphi",
191  Form("DST.MCPhoton.hcalin_iphi - %f", f->GetParameter(1)));
192 
193  const double center_hcalin_iphi = f->GetParameter(1);
194 
195  p = (TPad *) c1->cd(idx++);
196  c1->Update();
197  p->SetLogy();
198 
199  T->Draw("abs(MCPhoton_cor_hcalin_iphi)>>hMCPhoton_cor_hcalin_iphi(130,0,6.5)",
200  TString("DST.MCPhoton.hcalin_energy * (")
201  + TString(good_track_cut.GetTitle()) + ")");
202 
203  p->DrawFrame(0, 1e-3, 6.5, 1,
204  ";HCal_{IN} Azimuthal Distance (Tower Width);Energy Leakage Ratio");
205  TGraphErrors * ge = Distribution2Efficiency(hMCPhoton_cor_hcalin_iphi);
206  ge->SetLineColor(kBlue + 2);
207  ge->SetMarkerColor(kBlue + 21);
208  ge->SetMarkerColor(kFullCircle);
209  ge->SetLineWidth(3);
210  ge->Draw("lp");
211 
212  p = (TPad *) c1->cd(idx++);
213  c1->Update();
214  p->SetLogy();
215 
216  T->Draw("DST.MCPhoton.hcalin_ieta>>hhcalin_ieta(130,-6.5,6.5)",
217  TString("DST.MCPhoton.hcalin_energy * (")
218  + TString(good_track_cut.GetTitle()) + ")");
219  hhcalin_ieta->SetTitle(
220  ";HCal_{IN} Polar Distance (Tower Width);Energy Distribution");
221  hhcalin_ieta->Fit("gaus", "M");
222  TF1* f = (TF1*) (hhcalin_ieta->GetListOfFunctions()->At(0));
223  assert(f);
224  T->SetAlias("MCPhoton_cor_hcalin_ieta",
225  Form("DST.MCPhoton.hcalin_ieta - %f", f->GetParameter(1)));
226 
227  const double center_hcalin_ieta = f->GetParameter(1);
228 
229  p = (TPad *) c1->cd(idx++);
230  c1->Update();
231  p->SetLogy();
232 
233  T->Draw("abs(MCPhoton_cor_hcalin_ieta)>>hMCPhoton_cor_hcalin_ieta(130,0,6.5)",
234  TString("DST.MCPhoton.hcalin_energy * (")
235  + TString(good_track_cut.GetTitle()) + ")");
236 
237  p->DrawFrame(0, 1e-3, 6.5, 1,
238  ";HCal_{IN} Polar Distance (Tower Width);Energy Leakage Ratio");
239  TGraphErrors * ge = Distribution2Efficiency(hMCPhoton_cor_hcalin_ieta);
240  ge->SetLineColor(kBlue + 2);
241  ge->SetMarkerColor(kBlue + 21);
242  ge->SetMarkerColor(kFullCircle);
243  ge->SetLineWidth(3);
244  ge->Draw("lp");
245 
246  T->SetAlias("MCPhoton_cemc_e",
247  Form(
248  "1*Sum$( DST.MCPhoton.cemc_energy * ( sqrt(MCPhoton_cor_cemc_iphi*MCPhoton_cor_cemc_iphi + MCPhoton_cor_cemc_ieta*MCPhoton_cor_cemc_ieta) < %f ) )",
249  R_CEMC));
250  T->SetAlias("MCPhoton_hcalin_e",
251  Form(
252  "1*Sum$( DST.MCPhoton.hcalin_energy * ( sqrt(MCPhoton_cor_hcalin_iphi*MCPhoton_cor_hcalin_iphi + MCPhoton_cor_hcalin_ieta*MCPhoton_cor_hcalin_ieta) < %f ) )",
253  R_CEMC));
254 
255  T->SetAlias("MCPhoton_cemc_ntower",
256  Form(
257  "1*Sum$(( sqrt(MCPhoton_cor_cemc_iphi*MCPhoton_cor_cemc_iphi + MCPhoton_cor_cemc_ieta*MCPhoton_cor_cemc_ieta) < %f ) )",
258  R_HCALIN));
259  T->SetAlias("MCPhoton_hcalin_ntower",
260  Form(
261  "1*Sum$( ( sqrt(MCPhoton_cor_hcalin_iphi*MCPhoton_cor_hcalin_iphi + MCPhoton_cor_hcalin_ieta*MCPhoton_cor_hcalin_ieta) < %f ) )",
262  R_HCALIN));
263 
264  SaveCanvas(c1,
265  TString(_file0->GetName()) + TString("_DrawEcal_pDST_")
266  + Form("RCEMC%.1f_RCEMC%.1f_", R_CEMC, R_HCALIN)
267  + TString(c1->GetName()), kFALSE);
268 
269  cout << "///////////////////////////////////////////////" << endl;
270  cout << "// Projection center based on " << _file0->GetName() << endl;
271  cout << "///////////////////////////////////////////////" << endl;
272 
273  cout << "const double center_cemc_iphi = " << center_cemc_iphi << "; // "
274  << _file0->GetName() << endl;
275  cout << "const double center_cemc_ieta = " << center_cemc_ieta << "; // "
276  << _file0->GetName() << endl;
277  cout << "const double center_hcalin_iphi = " << center_hcalin_iphi << "; // "
278  << _file0->GetName() << endl;
279  cout << "const double center_hcalin_ieta = " << center_hcalin_ieta << "; // "
280  << _file0->GetName() << endl;
281 
282 }
283 
284 void
285 Edep_Checks(TString infile, const double R_CEMC, const double R_HCALIN,
286  TCut good_track_cut)
287 {
288  MakeRadiusCut(infile, R_CEMC, R_HCALIN, good_track_cut);
289 
290  double N_Event = T->GetEntries();
291 
292  TCanvas *c1 = new TCanvas("Edep_Checks" + cuts, "Edep_Checks" + cuts, 1900,
293  950);
294  c1->Divide(2, 2);
295  int idx = 1;
296  TPad * p;
297 
298  p = (TPad *) c1->cd(idx++);
299  c1->Update();
300  p->SetLogy();
301  T->Draw("MCPhoton_cemc_ntower>>hMCPhoton_cemc_ntower(16,-.5,15.5)",
302  good_track_cut);
303  hMCPhoton_cemc_ntower->SetTitle(
304  Form("CEMC Cluster Size (R = %.1f);Cluster Size (Towers);Probability",
305  R_CEMC));
306  hMCPhoton_cemc_ntower->Scale(1. / N_Event);
307 
308  p = (TPad *) c1->cd(idx++);
309  c1->Update();
310  p->SetLogy();
311  T->Draw("MCPhoton_hcalin_ntower>>hMCPhoton_hcalin_ntower(16,-.5,15.5)",
312  good_track_cut);
313  hMCPhoton_hcalin_ntower->SetTitle(
314  Form(
315  "HCal_{in} Cluster Size (R = %.1f);Cluster Size (Towers);Probability",
316  R_HCALIN));
317  hMCPhoton_hcalin_ntower->Scale(1. / N_Event);
318 
319  p = (TPad *) c1->cd(idx++);
320  c1->Update();
321  p->SetLogy();
322  T->Draw("MCPhoton_cemc_e>>hMCPhoton_cemc_e(240,-.0,12)", good_track_cut);
323  hMCPhoton_cemc_e->SetTitle(
324  Form(
325  "CEMC Cluster Energy (R = %.1f);Cluster Energy (/bin);Probability/bin",
326  R_CEMC));
327  hMCPhoton_cemc_e->Scale(1. / N_Event);
328 
329  p = (TPad *) c1->cd(idx++);
330  c1->Update();
331  p->SetLogy();
332  T->Draw("MCPhoton_hcalin_e>>hMCPhoton_hcalin_e(240,-.0,12)", good_track_cut);
333  hMCPhoton_hcalin_e->SetTitle(
334  Form(
335  "HCal_{in} Cluster Energy (R = %.1f);Cluster Energy (GeV);Probability/bin",
336  R_HCALIN));
337  hMCPhoton_hcalin_e->Scale(1. / N_Event);
338 
339  SaveCanvas(c1,
340  TString(_file0->GetName()) + TString("_DrawEcal_pDST_")
341  + Form("RCEMC%.1f_RCEMC%.1f_", R_CEMC, R_HCALIN)
342  + TString(c1->GetName()), kFALSE);
343 
344  TCanvas *c1 = new TCanvas("Edep_Checks_2D" + cuts, "Edep_Checks_2D" + cuts,
345  900, 900);
346 // c1->Divide(2, 2);
347 // int idx = 1;
348 // TPad * p;
349 
350  p = (TPad *) c1->cd(idx++);
351  c1->Update();
352  p->SetLogz();
353  T->Draw(
354  "MCPhoton_hcalin_e:MCPhoton_cemc_e>>h2_MCPhoton_hcalin_e_MCPhoton_cemc_e(240,-.0,12, 240,-.0,12)",
355  good_track_cut, "colz");
356  h2_MCPhoton_hcalin_e_MCPhoton_cemc_e->SetTitle(
357  Form(
358  "Energy distribution;CEMC Cluster Energy (R = %.1f) in GeV;HCal_{in} Cluster Energy (R = %.1f) in GeV",
359  R_CEMC, R_HCALIN));
360  h2_MCPhoton_hcalin_e_MCPhoton_cemc_e->Scale(1. / N_Event);
361  h2_MCPhoton_hcalin_e_MCPhoton_cemc_e->GetZaxis()->SetRangeUser(1e-6, 1);
362 
363  SaveCanvas(c1,
364  TString(_file0->GetName()) + TString("_DrawEcal_pDST_")
365  + Form("RCEMC%.1f_RCEMC%.1f_", R_CEMC, R_HCALIN)
366  + TString(c1->GetName()), kFALSE);
367 
368 }
369 
370 void
371 TrackProjection_Checks(TString infile, TCut good_track_cut)
372 {
373 
374  double N_Event = T->GetEntries();
375 
376  TCanvas *c1 = new TCanvas("TrackProjection_Checks_Raw" + cuts,
377  "TrackProjection_Checks_Raw" + cuts, 1800, 900);
378  c1->Divide(2, 1);
379  int idx = 1;
380  TPad * p;
381 
382  p = (TPad *) c1->cd(idx++);
383  c1->Update();
384  p->SetLogz();
385  T->Draw(
386  "DST.MCPhoton.cemc_iphi:DST.MCPhoton.cemc_ieta>>hcemc(130,-6.5,6.5,130,-6.5,6.5)",
387  "DST.MCPhoton.cemc_energy", "colz");
388  hcemc->SetTitle(
389  "CEMC Tower Energy Distribution;Polar Distance (Tower Width);Azimuthal Distance (Tower Width)");
390  hcemc->Scale(1. / N_Event);
391  hcemc->GetZaxis()->SetRangeUser(1e-5, 30);
392 
393  p = (TPad *) c1->cd(idx++);
394  c1->Update();
395  p->SetLogz();
396  T->Draw(
397  "DST.MCPhoton.hcalin_iphi:DST.MCPhoton.hcalin_ieta>>hcalin(130,-6.5,6.5,130,-6.5,6.5)",
398  "DST.MCPhoton.hcalin_energy", "colz");
399  hcalin->SetTitle(
400  "HCal_{in} Tower Energy Distribution;Polar Distance (Tower Width);Azimuthal Distance (Tower Width)");
401  hcalin->Scale(1. / N_Event);
402  hcalin->GetZaxis()->SetRangeUser(1e-5, 30);
403 
404  SaveCanvas(c1,
405  TString(_file0->GetName()) + TString("_DrawEcal_pDST_")
406  + TString(c1->GetName()), kFALSE);
407 
408 
409 
410  TCanvas *c1 = new TCanvas("TrackProjection_Checks_Proj" + cuts,
411  "TrackProjection_Checks_Proj" + cuts, 1800, 900);
412  c1->Divide(2, 2);
413  int idx = 1;
414  TPad * p;
415 
416  p = (TPad *) c1->cd(idx++);
417  c1->Update();
418  p->SetLogy();
419 
420  TH1 * proj = hcemc->ProjectionX();
421  proj->SetLineColor(kBlack);
422  proj->Draw("");
423 // TH1 * proj = hcemcc->ProjectionX();
424 // proj->SetLineColor(kRed);
425 // proj->Draw("same");
426 
427  p = (TPad *) c1->cd(idx++);
428  c1->Update();
429  p->SetLogy();
430 
431  TH1 * proj = hcalin->ProjectionX();
432  proj->SetLineColor(kBlack);
433  proj->Draw("");
434 // TH1 * proj = hcalinc->ProjectionX();
435 // proj->SetLineColor(kRed);
436 // proj->Draw("same");
437 
438  p = (TPad *) c1->cd(idx++);
439  c1->Update();
440  p->SetLogy();
441 
442  TH1 * proj = hcemc->ProjectionY();
443  proj->SetLineColor(kBlack);
444  proj->Draw("");
445 // TH1 * proj = hcemcc->ProjectionY();
446 // proj->SetLineColor(kRed);
447 // proj->Draw("same");
448 
449  p = (TPad *) c1->cd(idx++);
450  c1->Update();
451  p->SetLogy();
452 
453  TH1 * proj = hcalin->ProjectionY();
454  proj->SetLineColor(kBlack);
455  proj->Draw("");
456 // TH1 * proj = hcalinc->ProjectionY();
457 // proj->SetLineColor(kRed);
458 // proj->Draw("same");
459 
460  SaveCanvas(c1,
461  TString(_file0->GetName()) + TString("_DrawEcal_pDST_")
462  + TString(c1->GetName()), kFALSE);
463 }
464 
465 void
466 UpsilonPair_Checks(TString infile)
467 {
468 
469  TCanvas *c1 = new TCanvas("UpsilonPair_Checks" + cuts,
470  "UpsilonPair_Checks" + cuts, 1800, 900);
471  c1->Divide(3, 2);
472  int idx = 1;
473  TPad * p;
474 
475  p = (TPad *) c1->cd(idx++);
476  c1->Update();
477  T->Draw("DST.UpsilonPair.gmass>>gmass(100,0,10)");
478  T->Draw("DST.UpsilonPair.gmass>>gmassc(100,0,10)",
479  "DST.UpsilonPair.good_upsilon", "same");
480  gmassc->SetLineColor(kRed);
481 
482  p = (TPad *) c1->cd(idx++);
483  c1->Update();
484  T->Draw("DST.UpsilonPair.mass>>mass(100,0,10)");
485  T->Draw("DST.UpsilonPair.mass>>massc(100,0,10)",
486  "DST.UpsilonPair.good_upsilon", "same");
487  massc->SetLineColor(kRed);
488 
489  p = (TPad *) c1->cd(idx++);
490  c1->Update();
491  T->Draw("UpsilonPair_trk_gpt>>UpsilonPair_trk_gpt(100,0,10)");
492 
493  p = (TPad *) c1->cd(idx++);
494  c1->Update();
495  T->Draw("UpsilonPair_trk_pt:UpsilonPair_trk_gpt>>pt2(100,0,10,100,0,10)", "",
496  "colz");
497 
498  p = (TPad *) c1->cd(idx++);
499  c1->Update();
500  T->Draw("UpsilonPair_trk_pt:UpsilonPair_trk_gpt>>pt2c(100,0,10,100,0,10)",
501  "DST.UpsilonPair.good_upsilon", "colz");
502 
503  SaveCanvas(c1,
504  TString(_file0->GetName()) + TString("_DrawEcal_pDST_")
505  + TString(c1->GetName()), kFALSE);
506 }
507 
508 TGraphErrors *
509 Distribution2Efficiency(TH1F * hCEMC3_Max)
510 {
511  double threshold[10000] =
512  { 0 };
513  double eff[10000] =
514  { 0 };
515  double eff_err[10000] =
516  { 0 };
517 
518  assert(hCEMC3_Max->GetNbinsX()<10000);
519 
520  const double n = hCEMC3_Max->GetSum();
521  double pass = 0;
522  int cnt = 0;
523  for (int i = hCEMC3_Max->GetNbinsX(); i >= 1; i--)
524  {
525  pass += hCEMC3_Max->GetBinContent(i);
526 
527  const double pp = pass / n;
528 // const double z = 1.96;
529  const double z = 1.;
530 
531  const double A = z * sqrt(1. / n * pp * (1 - pp) + z * z / 4 / n / n);
532  const double B = 1 / (1 + z * z / n);
533 
534  threshold[cnt] = hCEMC3_Max->GetBinCenter(i);
535  eff[cnt] = (pp + z * z / 2 / n) * B;
536  eff_err[cnt] = A * B;
537 
538  cout << threshold[cnt] << ": " << "CL " << eff[cnt] << "+/-"
539  << eff_err[cnt] << endl;
540  cnt++;
541  }
542  TGraphErrors * ge = new TGraphErrors(cnt, threshold, eff, NULL, eff_err);
543 
544  return ge;
545 }
546 
547 void
548 TrackProjection_Checks_Comparison()
549 {
550 
551  SetOKStyle();
552  gStyle->SetOptStat(0);
553  gStyle->SetOptFit(1111);
554  TVirtualFitter::SetDefaultFitter("Minuit2");
555 
556  gSystem->Load("libg4eval.so");
557  gSystem->Load("libemcal_ana.so");
558  gSystem->Load("libg4vertex.so");
559 
560 
561  TFile * f =
562  new TFile(
563  "/phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root_DrawEcal_pDST_TrackProjection_Checks_Cut_all_event.root");
564 
565  assert(f->IsOpen());
566  TH2F * hcemcc_2d = (TH2F *) f->GetObjectChecked("hcemcc", "TH2F");
567  assert(hcemcc_2d);
568 
569  TFile * f =
570  new TFile(
571  "/phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root_DrawEcal_pDST_TrackProjection_Checks_Cut_all_event.root");
572 
573  assert(f->IsOpen());
574  TH2F * hcemcc_1d = (TH2F *) f->GetObjectChecked("hcemcc", "TH2F");
575  assert(hcemcc_1d);
576 
577  TCanvas *c1 = new TCanvas("TrackProjection_Checks_Comparison" ,
578  "TrackProjection_Checks_Comparison" , 1800, 1000);
579  c1->Divide(2, 1);
580  int idx = 1;
581  TPad * p;
582 
583  p = (TPad *) c1->cd(idx++);
584  c1->Update();
585  p->SetLogz();
586 
587  hcemcc_2d -> Draw("colz");
588  hcemcc_2d->SetTitle(";Polar Distance (Tower Width);Azimuthal Distance (Tower Width)");
589 
590  TLatex *text = new TLatex(0,7,"2-D Projective SPACAL, 0.9 < #eta_{e^{-}} < 1.0, E_{e^{-}} = 8 GeV");
591  text->SetTextAlign(22);
592  text->SetTextSize(0.044);
593  text->Draw();
594 
595 
596  p = (TPad *) c1->cd(idx++);
597  c1->Update();
598  p->SetLogz();
599 
600  hcemcc_1d -> Draw("colz");
601  hcemcc_1d->SetTitle(";Polar Distance (Tower Width);Azimuthal Distance (Tower Width)");
602 
603  TLatex *text = new TLatex(0,7,"1-D Projective SPACAL, 0.9 < #eta_{e^{-}} < 1.0, E_{e^{-}} = 8 GeV");
604  text->SetTextAlign(22);
605  text->SetTextSize(0.044);
606  text->Draw();
607 
608 
609  SaveCanvas(c1,
610  "/phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/" + TString("DrawEcal_pDST_")
611  + TString(c1->GetName()));
612 }
613 
614 void
615 Logs()
616 {
617 
619  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0_8GeV-ALL.root_Ana.root
621  const double center_cemc_iphi = 0.0202159; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0_8GeV-ALL.root_Ana.root
622  const double center_cemc_ieta = 0.00453938; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0_8GeV-ALL.root_Ana.root
623  const double center_hcalin_iphi = 0.488432; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0_8GeV-ALL.root_Ana.root
624  const double center_hcalin_ieta = -0.00936002; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0_8GeV-ALL.root_Ana.root
625 
627  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0_8GeV-ALL.root_Ana.root
629  const double center_cemc_iphi = 0.0402477; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0_8GeV-ALL.root_Ana.root
630  const double center_cemc_ieta = 0.00574989; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0_8GeV-ALL.root_Ana.root
631  const double center_hcalin_iphi = 0.331081; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0_8GeV-ALL.root_Ana.root
632  const double center_hcalin_ieta = -0.000527382; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0_8GeV-ALL.root_Ana.root
633 
634 
636  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_8GeV-ALL.root_Ana.root
638  const double center_cemc_iphi = 0.0649585; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_8GeV-ALL.root_Ana.root
639  const double center_cemc_ieta = 0.00465546; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_8GeV-ALL.root_Ana.root
640  const double center_hcalin_iphi = 0.552238; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_8GeV-ALL.root_Ana.root
641  const double center_hcalin_ieta = -0.0102162; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_8GeV-ALL.root_Ana.root
642 
644  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_8GeV-ALL.root_Ana.root
646  const double center_cemc_iphi = 0.044375; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_8GeV-ALL.root_Ana.root
647  const double center_cemc_ieta = 0.00680547; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_8GeV-ALL.root_Ana.root
648  const double center_hcalin_iphi = 0.336741; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_8GeV-ALL.root_Ana.root
649  const double center_hcalin_ieta = 0.00108616; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_8GeV-ALL.root_Ana.root
650 
652  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_4GeV-ALL.root_Ana.root
654  const double center_cemc_iphi = 0.0952638; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_4GeV-ALL.root_Ana.root
655  const double center_cemc_ieta = 0.00397331; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_4GeV-ALL.root_Ana.root
656  const double center_hcalin_iphi = 0.592338; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_4GeV-ALL.root_Ana.root
657  const double center_hcalin_ieta = -0.0118107; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_4GeV-ALL.root_Ana.root
658 
660  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_4GeV-ALL.root_Ana.root
662  const double center_cemc_iphi = 0.0577477; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_4GeV-ALL.root_Ana.root
663  const double center_cemc_ieta = 0.00609392; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_4GeV-ALL.root_Ana.root
664  const double center_hcalin_iphi = 0.354109; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_4GeV-ALL.root_Ana.root
665  const double center_hcalin_ieta = -0.000273002; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_4GeV-ALL.root_Ana.root
666 
668  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_2GeV-ALL.root_Ana.root
670  const double center_cemc_iphi = 0.170699; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_2GeV-ALL.root_Ana.root
671  const double center_cemc_ieta = 0.00345903; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_2GeV-ALL.root_Ana.root
672  const double center_hcalin_iphi = 0.668343; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_2GeV-ALL.root_Ana.root
673  const double center_hcalin_ieta = -0.0130487; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_2GeV-ALL.root_Ana.root
674 
676  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_2GeV-ALL.root_Ana.root
678  const double center_cemc_iphi = 0.0898377; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_2GeV-ALL.root_Ana.root
679  const double center_cemc_ieta = 0.00445598; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_2GeV-ALL.root_Ana.root
680  const double center_hcalin_iphi = 0.376608; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_2GeV-ALL.root_Ana.root
681  const double center_hcalin_ieta = -0.0074984; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_2GeV-ALL.root_Ana.root
682 
684  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
686  const double center_cemc_iphi = 0.122055; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
687  const double center_cemc_ieta = 1.52588; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
688  const double center_hcalin_iphi = 0.644177; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
689  const double center_hcalin_ieta = -0.799621; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
691  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
693  const double center_cemc_iphi = 0.122055; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
694  const double center_cemc_ieta = 1.52588; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
695  const double center_hcalin_iphi = 0.644177; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
696  const double center_hcalin_ieta = -0.799621; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
698  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root
700  const double center_cemc_iphi = 0.251162; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root
701  const double center_cemc_ieta = 1.35067; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root
702  const double center_hcalin_iphi = 0.761301; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root
703  const double center_hcalin_ieta = -0.874644; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root
705  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root
707  const double center_cemc_iphi = 0.177891; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root
708  const double center_cemc_ieta = 2.62604; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root
709  const double center_hcalin_iphi = 0.393433; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root
710  const double center_hcalin_ieta = -0.116632; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root
712  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root
714  const double center_cemc_iphi = 0.53236; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root
715  const double center_cemc_ieta = 1.15917; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root
716  const double center_hcalin_iphi = 0.899527; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root
717  const double center_hcalin_ieta = -0.889783; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root
719  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root
721  const double center_cemc_iphi = 0.378081; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root
722  const double center_cemc_ieta = 2.47524; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root
723  const double center_hcalin_iphi = 0.445272; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root
724  const double center_hcalin_ieta = -0.180098; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root
725 
727  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root
729  const double center_cemc_iphi = 0.129035; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root
730  const double center_cemc_ieta = 0.0222246; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root
731  const double center_hcalin_iphi = 0.701734; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root
732  const double center_hcalin_ieta = -0.474402; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root
733 
735  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root
737  const double center_cemc_iphi = 0.0736277; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root
738  const double center_cemc_ieta = 0.0556152; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root
739  const double center_hcalin_iphi = 0.385379; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root
740  const double center_hcalin_ieta = -0.0596968; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root
741 
743  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root
745  const double center_cemc_iphi = 0.249273; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root
746  const double center_cemc_ieta = 0.0146006; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root
747  const double center_hcalin_iphi = 0.843407; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root
748  const double center_hcalin_ieta = -0.592856; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root
749 
751  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root
753  const double center_cemc_iphi = 0.12442; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root
754  const double center_cemc_ieta = 0.0514677; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root
755  const double center_hcalin_iphi = 0.447927; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root
756  const double center_hcalin_ieta = -0.101503; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root
757 
759  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
761  const double center_cemc_iphi = 0.0829824; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
762  const double center_cemc_ieta = 0.0275653; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
763  const double center_hcalin_iphi = 0.635325; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
764  const double center_hcalin_ieta = -0.408261; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
765 
767  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_8GeV-ALL.root_Ana.root
769  const double center_cemc_iphi = 0.054392; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_8GeV-ALL.root_Ana.root
770  const double center_cemc_ieta = 0.0852583; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_8GeV-ALL.root_Ana.root
771  const double center_hcalin_iphi = 0.348779; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_8GeV-ALL.root_Ana.root
772  const double center_hcalin_ieta = -0.0348952; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_8GeV-ALL.root_Ana.root
773 
774 
775 
776 
778  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_2GeV-ALL.root_Ana.root
780  const double center_cemc_iphi = -0.16237; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_2GeV-ALL.root_Ana.root
781  const double center_cemc_ieta = 0.00993428; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_2GeV-ALL.root_Ana.root
782  const double center_hcalin_iphi = 0.312564; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_2GeV-ALL.root_Ana.root
783  const double center_hcalin_ieta = -0.568171; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_2GeV-ALL.root_Ana.root
784 
786  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_2GeV-ALL.root_Ana.root
788  const double center_cemc_iphi = -0.039828; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_2GeV-ALL.root_Ana.root
789  const double center_cemc_ieta = 0.0489088; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_2GeV-ALL.root_Ana.root
790  const double center_hcalin_iphi = 0.338006; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_2GeV-ALL.root_Ana.root
791  const double center_hcalin_ieta = -0.0939571; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_2GeV-ALL.root_Ana.root
792 
794  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_4GeV-ALL.root_Ana.root
796  const double center_cemc_iphi = -0.0455698; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_4GeV-ALL.root_Ana.root
797  const double center_cemc_ieta = 0.0181539; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_4GeV-ALL.root_Ana.root
798  const double center_hcalin_iphi = 0.444123; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_4GeV-ALL.root_Ana.root
799  const double center_hcalin_ieta = -0.475483; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_4GeV-ALL.root_Ana.root
800 
802  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_4GeV-ALL.root_Ana.root
804  const double center_cemc_iphi = 0.0114778; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_4GeV-ALL.root_Ana.root
805  const double center_cemc_ieta = 0.0580504; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_4GeV-ALL.root_Ana.root
806  const double center_hcalin_iphi = 0.340425; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_4GeV-ALL.root_Ana.root
807  const double center_hcalin_ieta = -0.0670486; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_4GeV-ALL.root_Ana.root
808 
809 
811  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_8GeV-ALL.root_Ana.root
813  const double center_cemc_iphi = 0.00278605; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_8GeV-ALL.root_Ana.root
814  const double center_cemc_ieta = 0.0237387; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_8GeV-ALL.root_Ana.root
815  const double center_hcalin_iphi = 0.504268; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_8GeV-ALL.root_Ana.root
816  const double center_hcalin_ieta = -0.401039; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_8GeV-ALL.root_Ana.root
817 
819  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_8GeV-ALL.root_Ana.root
821  const double center_cemc_iphi = 0.0313807; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_8GeV-ALL.root_Ana.root
822  const double center_cemc_ieta = 0.0800593; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_8GeV-ALL.root_Ana.root
823  const double center_hcalin_iphi = 0.336144; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_8GeV-ALL.root_Ana.root
824  const double center_hcalin_ieta = -0.0484305; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_8GeV-ALL.root_Ana.root
825 
826 
827 
828 
830  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_8GeV-ALL.root_Ana.root
832  const double center_cemc_iphi = -0.122248; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_8GeV-ALL.root_Ana.root
833  const double center_cemc_ieta = 1.52555; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_8GeV-ALL.root_Ana.root
834  const double center_hcalin_iphi = 0.500768; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_8GeV-ALL.root_Ana.root
835  const double center_hcalin_ieta = -0.802395; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_8GeV-ALL.root_Ana.root
837  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_8GeV-ALL.root_Ana.root
839  const double center_cemc_iphi = -0.0819904; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_8GeV-ALL.root_Ana.root
840  const double center_cemc_ieta = 2.83299; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_8GeV-ALL.root_Ana.root
841  const double center_hcalin_iphi = 0.33193; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_8GeV-ALL.root_Ana.root
842  const double center_hcalin_ieta = -0.0749197; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_8GeV-ALL.root_Ana.root
843 
845  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_4GeV-ALL.root_Ana.root
847  const double center_cemc_iphi = -0.252243; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_4GeV-ALL.root_Ana.root
848  const double center_cemc_ieta = 1.34503; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_4GeV-ALL.root_Ana.root
849  const double center_hcalin_iphi = 0.411387; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_4GeV-ALL.root_Ana.root
850  const double center_hcalin_ieta = -0.860032; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_4GeV-ALL.root_Ana.root
851 
853  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_4GeV-ALL.root_Ana.root
855  const double center_cemc_iphi = -0.171931; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_4GeV-ALL.root_Ana.root
856  const double center_cemc_ieta = 2.65013; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_4GeV-ALL.root_Ana.root
857  const double center_hcalin_iphi = 0.334954; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_4GeV-ALL.root_Ana.root
858  const double center_hcalin_ieta = -0.0992204; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_4GeV-ALL.root_Ana.root
859 
861  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_2GeV-ALL.root_Ana.root
863  const double center_cemc_iphi = -0.530751; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_2GeV-ALL.root_Ana.root
864  const double center_cemc_ieta = 1.15259; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_2GeV-ALL.root_Ana.root
865  const double center_hcalin_iphi = 0.228558; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_2GeV-ALL.root_Ana.root
866  const double center_hcalin_ieta = -0.923731; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_2GeV-ALL.root_Ana.root
867 
869  // Projection center based on /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_2GeV-ALL.root_Ana.root
871  const double center_cemc_iphi = -0.371817; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_2GeV-ALL.root_Ana.root
872  const double center_cemc_ieta = 2.46165; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_2GeV-ALL.root_Ana.root
873  const double center_hcalin_iphi = 0.317738; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_2GeV-ALL.root_Ana.root
874  const double center_hcalin_ieta = -0.119789; // /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_2GeV-ALL.root_Ana.root
875 
876 
877 
878 }