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 = "../../sPHENIX_work//test_production/Upsilon/spacal2d/fieldon/SimALL_PythiaUpsilon.root_Ana.root"
28 
29 // const TString infile = "../../sPHENIX_work/production_analysis//emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_8GeV-ALL.root_Ana.root.lst_EMCalLikelihood.root"
30  const TString infile = "../../sPHENIX_work/production_analysis_cemc2x2//emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_8GeV-ALL.root_Ana.root.lst_EMCalLikelihood.root"
31 
32 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_8GeV-ALL.root_Ana.root"
33 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_8GeV-ALL.root_Ana.root"
34 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_4GeV-ALL.root_Ana.root"
35 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_4GeV-ALL.root_Ana.root"
36 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_2GeV-ALL.root_Ana.root"
37 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_2GeV-ALL.root_Ana.root"
38 
39 // const TString infile = "../../sPHENIX_work/production_analysis/embedding/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_8GeV-ALL.root_Ana.root"
40 // const TString infile = "../../sPHENIX_work/production_analysis/embedding/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_8GeV-ALL.root_Ana.root"
41 // const TString infile = "../../sPHENIX_work/production_analysis/embedding/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_4GeV-ALL.root_Ana.root"
42 // const TString infile = "../../sPHENIX_work/production_analysis/embedding/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_4GeV-ALL.root_Ana.root"
43 // const TString infile = "../../sPHENIX_work/production_analysis/embedding/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_2GeV-ALL.root_Ana.root"
44 // const TString infile = "../../sPHENIX_work/production_analysis/embedding/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_2GeV-ALL.root_Ana.root"
45 
46 // const TString infile = "../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root"
47 // const TString infile = "../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_8GeV-ALL.root_Ana.root"
48 // const TString infile = "../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root"
49 // const TString infile = "../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root"
50 // const TString infile = "../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root"
51 // const TString infile = "../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root"
52  //
53  // const TString infile = "../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root"
54  // const TString infile = "../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_8GeV-ALL.root_Ana.root"
55  // const TString infile = "../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root"
56  // const TString infile = "../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root"
57  // const TString infile = "../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root"
58  // const TString infile = "../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root"
59  //
60 // const TString infile = "../../sPHENIX_work/production_analysis/embedding/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_8GeV-ALL.root_Ana.root"
61 // const TString infile = "../../sPHENIX_work/production_analysis/embedding/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root"
62 // const TString infile = "../../sPHENIX_work/production_analysis/embedding/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root"
63 // const TString infile = "../../sPHENIX_work/production_analysis/embedding/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root"
64 // const TString infile = "../../sPHENIX_work/production_analysis/embedding/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root"
65 // const TString infile = "../../sPHENIX_work/production_analysis/embedding/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root"
66 
67  // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_8GeV-ALL.root_Ana.root"
68  // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_8GeV-ALL.root_Ana.root"
69  // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_4GeV-ALL.root_Ana.root"
70  // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_4GeV-ALL.root_Ana.root"
71  // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_2GeV-ALL.root_Ana.root"
72  // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_2GeV-ALL.root_Ana.root"
73 
74 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_8GeV-ALL.root_Ana.root"
75 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_8GeV-ALL.root_Ana.root"
76 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_4GeV-ALL.root_Ana.root"
77 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_4GeV-ALL.root_Ana.root"
78 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_2GeV-ALL.root_Ana.root"
79 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_2GeV-ALL.root_Ana.root"
80 
81 // const TString infile = "G4Hits_sPHENIX_e-_eta0_8GeV-0000.root_Ana.root" //
82 // const TString infile = "G4Hits_sPHENIX_pi-_eta0_8GeV-0000.root_Ana.root"//
83  )
84 {
85  SetOKStyle();
86  gStyle->SetOptStat(0);
87  gStyle->SetOptFit(1111);
88  TVirtualFitter::SetDefaultFitter("Minuit2");
89 
90  gSystem->Load("libg4eval.so");
91  gSystem->Load("libemcal_ana.so");
92  gSystem->Load("libg4vertex.so");
93 
94  if (!_file0)
95  {
96  TString chian_str = infile;
97  chian_str.ReplaceAll("ALL", "*");
98  chian_str.ReplaceAll("+", "\\+");
99 
100  TChain * t = new TChain("T");
101  const int n = t->Add(chian_str);
102 
103  cout << "Loaded " << n << " root files with " << chian_str << endl;
104  assert(n>0);
105 
106  T = t;
107 
108  _file0 = new TFile;
109  _file0->SetName(infile);
110 
111  fstream flst(infile + ".lst", ios_base::out);
112 
113  TObjArray *fileElements = t->GetListOfFiles();
114  TIter next(fileElements);
115  TChainElement *chEl = 0;
116  while ((chEl = (TChainElement*) next()))
117  {
118  flst << chEl->GetTitle() << endl;
119  }
120  flst.close();
121 
122  cout << "Saved file list to " << infile + ".lst" << endl;
123  }
124 
125  assert(_file0);
126 
127  T->SetAlias("UpsilonPair_trk_gpt",
128  "1*sqrt(DST.UpsilonPair.trk.gpx**2 + DST.UpsilonPair.trk.gpy**2)");
129  T->SetAlias("UpsilonPair_trk_pt",
130  "1*sqrt(DST.UpsilonPair.trk.px**2 + DST.UpsilonPair.trk.py**2)");
131 
132  T->SetAlias("EMCalTrk_pt", "1*sqrt(DST.EMCalTrk.px**2 + DST.EMCalTrk.py**2)");
133  T->SetAlias("EMCalTrk_gpt",
134  "1*sqrt(DST.EMCalTrk.gpx**2 + DST.EMCalTrk.gpy**2)");
135 
136  const TCut event_sel = "1*1";
137  cuts = "_all_event";
138 
139  cout << "Build event selection of " << (const char *) event_sel << endl;
140 
141  T->Draw(">>EventList", event_sel);
142  TEventList * elist = gDirectory->GetObjectChecked("EventList", "TEventList");
143  cout << elist->GetN() << " / " << T->GetEntriesFast() << " events selected"
144  << endl;
145 
146  T->SetEventList(elist);
148 // UpsilonPair_Checks(infile);
149 
150 // TrackProjection_Checks(infile, "abs(EMCalTrk_pt/EMCalTrk_gpt - 1)<0.05");
151 // Edep_Checks(infile, 1.4, 1.4, "abs(EMCalTrk_pt/EMCalTrk_gpt - 1)<0.05");
152 
153  TrackProjection_Checks(infile, "abs(EMCalTrk_pt/EMCalTrk_gpt - 1)<0.05 && DST.EMCalTrk.cemce>2");
154  Edep_Checks(infile, 1.4, 1.4, "abs(EMCalTrk_pt/EMCalTrk_gpt - 1)<0.05 && DST.EMCalTrk.cemce>2");
155 }
156 
157 void
158 MakeRadiusCut(TString infile, const double R_CEMC, const double R_HCALIN,
159  TCut good_track_cut)
160 {
161  TCanvas *c1 = new TCanvas("MakeRadiusCut" + cuts, "MakeRadiusCut" + cuts,
162  1800, 900);
163  c1->Divide(4, 2);
164  int idx = 1;
165  TPad * p;
166 
167  p = (TPad *) c1->cd(idx++);
168  c1->Update();
169  p->SetLogy();
170 
171  T->Draw("DST.EMCalTrk.cemc_iphi>>hcemc_iphi(130,-6.5,6.5)",
172  TString("DST.EMCalTrk.cemc_energy * (")
173  + TString(good_track_cut.GetTitle()) + ")");
174  hcemc_iphi->SetTitle(
175  ";CEMC Azimuthal Distance (Tower Width);Energy Distribution");
176  hcemc_iphi->Fit("gaus", "M");
177  TF1* f = (TF1*) (hcemc_iphi->GetListOfFunctions()->At(0));
178  assert(f);
179  T->SetAlias("EMCalTrk_cor_cemc_iphi",
180  Form("DST.EMCalTrk.cemc_iphi - %f", f->GetParameter(1)));
181 
182  const double center_cemc_iphi = f->GetParameter(1);
183 
184  p = (TPad *) c1->cd(idx++);
185  c1->Update();
186  p->SetLogy();
187 
188  T->Draw("abs(EMCalTrk_cor_cemc_iphi)>>hEMCalTrk_cor_cemc_iphi(130,0,6.5)",
189  TString("DST.EMCalTrk.cemc_energy * (")
190  + TString(good_track_cut.GetTitle()) + ")");
191 
192  p->DrawFrame(0, 1e-3, 6.5, 1,
193  ";CEMC Azimuthal Distance (Tower Width);Energy Leakage Ratio");
194  TGraphErrors * ge = Distribution2Efficiency(hEMCalTrk_cor_cemc_iphi);
195  ge->SetLineColor(kBlue + 2);
196  ge->SetMarkerColor(kBlue + 21);
197  ge->SetMarkerColor(kFullCircle);
198  ge->SetLineWidth(3);
199  ge->Draw("lp");
200 
201  p = (TPad *) c1->cd(idx++);
202  c1->Update();
203  p->SetLogy();
204 
205  T->Draw("DST.EMCalTrk.cemc_ieta>>hcemc_ieta(130,-6.5,6.5)",
206  TString("DST.EMCalTrk.cemc_energy * (")
207  + TString(good_track_cut.GetTitle()) + ")");
208  hcemc_ieta->SetTitle(
209  ";CEMC Polar Distance (Tower Width);Energy Distribution");
210  hcemc_ieta->Fit("gaus", "M");
211  TF1* f = (TF1*) (hcemc_ieta->GetListOfFunctions()->At(0));
212  assert(f);
213  T->SetAlias("EMCalTrk_cor_cemc_ieta",
214  Form("DST.EMCalTrk.cemc_ieta - %f", f->GetParameter(1)));
215 
216  const double center_cemc_ieta = f->GetParameter(1);
217 
218  p = (TPad *) c1->cd(idx++);
219  c1->Update();
220  p->SetLogy();
221 
222  T->Draw("abs(EMCalTrk_cor_cemc_ieta)>>hEMCalTrk_cor_cemc_ieta(130,0,6.5)",
223  TString("DST.EMCalTrk.cemc_energy * (")
224  + TString(good_track_cut.GetTitle()) + ")");
225 
226  p->DrawFrame(0, 1e-3, 6.5, 1,
227  ";CEMC Polar Distance (Tower Width);Energy Leakage Ratio");
228  TGraphErrors * ge = Distribution2Efficiency(hEMCalTrk_cor_cemc_ieta);
229  ge->SetLineColor(kBlue + 2);
230  ge->SetMarkerColor(kBlue + 21);
231  ge->SetMarkerColor(kFullCircle);
232  ge->SetLineWidth(3);
233  ge->Draw("lp");
234 
235  p = (TPad *) c1->cd(idx++);
236  c1->Update();
237  p->SetLogy();
238 
239  T->Draw("DST.EMCalTrk.hcalin_iphi>>hhcalin_iphi(130,-6.5,6.5)",
240  TString("DST.EMCalTrk.hcalin_energy * (")
241  + TString(good_track_cut.GetTitle()) + ")");
242  hhcalin_iphi->SetTitle(
243  ";HCal_{IN} Azimuthal Distance (Tower Width);Energy Distribution");
244  hhcalin_iphi->Fit("gaus", "M");
245  TF1* f = (TF1*) (hhcalin_iphi->GetListOfFunctions()->At(0));
246  assert(f);
247  T->SetAlias("EMCalTrk_cor_hcalin_iphi",
248  Form("DST.EMCalTrk.hcalin_iphi - %f", f->GetParameter(1)));
249 
250  const double center_hcalin_iphi = f->GetParameter(1);
251 
252  p = (TPad *) c1->cd(idx++);
253  c1->Update();
254  p->SetLogy();
255 
256  T->Draw("abs(EMCalTrk_cor_hcalin_iphi)>>hEMCalTrk_cor_hcalin_iphi(130,0,6.5)",
257  TString("DST.EMCalTrk.hcalin_energy * (")
258  + TString(good_track_cut.GetTitle()) + ")");
259 
260  p->DrawFrame(0, 1e-3, 6.5, 1,
261  ";HCal_{IN} Azimuthal Distance (Tower Width);Energy Leakage Ratio");
262  TGraphErrors * ge = Distribution2Efficiency(hEMCalTrk_cor_hcalin_iphi);
263  ge->SetLineColor(kBlue + 2);
264  ge->SetMarkerColor(kBlue + 21);
265  ge->SetMarkerColor(kFullCircle);
266  ge->SetLineWidth(3);
267  ge->Draw("lp");
268 
269  p = (TPad *) c1->cd(idx++);
270  c1->Update();
271  p->SetLogy();
272 
273  T->Draw("DST.EMCalTrk.hcalin_ieta>>hhcalin_ieta(130,-6.5,6.5)",
274  TString("DST.EMCalTrk.hcalin_energy * (")
275  + TString(good_track_cut.GetTitle()) + ")");
276  hhcalin_ieta->SetTitle(
277  ";HCal_{IN} Polar Distance (Tower Width);Energy Distribution");
278  hhcalin_ieta->Fit("gaus", "M");
279  TF1* f = (TF1*) (hhcalin_ieta->GetListOfFunctions()->At(0));
280  assert(f);
281  T->SetAlias("EMCalTrk_cor_hcalin_ieta",
282  Form("DST.EMCalTrk.hcalin_ieta - %f", f->GetParameter(1)));
283 
284  const double center_hcalin_ieta = f->GetParameter(1);
285 
286  p = (TPad *) c1->cd(idx++);
287  c1->Update();
288  p->SetLogy();
289 
290  T->Draw("abs(EMCalTrk_cor_hcalin_ieta)>>hEMCalTrk_cor_hcalin_ieta(130,0,6.5)",
291  TString("DST.EMCalTrk.hcalin_energy * (")
292  + TString(good_track_cut.GetTitle()) + ")");
293 
294  p->DrawFrame(0, 1e-3, 6.5, 1,
295  ";HCal_{IN} Polar Distance (Tower Width);Energy Leakage Ratio");
296  TGraphErrors * ge = Distribution2Efficiency(hEMCalTrk_cor_hcalin_ieta);
297  ge->SetLineColor(kBlue + 2);
298  ge->SetMarkerColor(kBlue + 21);
299  ge->SetMarkerColor(kFullCircle);
300  ge->SetLineWidth(3);
301  ge->Draw("lp");
302 
303  T->SetAlias("EMCalTrk_cemc_e",
304  Form(
305  "1*Sum$( DST.EMCalTrk.cemc_energy * ( sqrt(EMCalTrk_cor_cemc_iphi*EMCalTrk_cor_cemc_iphi + EMCalTrk_cor_cemc_ieta*EMCalTrk_cor_cemc_ieta) < %f ) )",
306  R_CEMC));
307  T->SetAlias("EMCalTrk_hcalin_e",
308  Form(
309  "1*Sum$( DST.EMCalTrk.hcalin_energy * ( sqrt(EMCalTrk_cor_hcalin_iphi*EMCalTrk_cor_hcalin_iphi + EMCalTrk_cor_hcalin_ieta*EMCalTrk_cor_hcalin_ieta) < %f ) )",
310  R_CEMC));
311 
312  T->SetAlias("EMCalTrk_cemc_ntower",
313  Form(
314  "1*Sum$(( sqrt(EMCalTrk_cor_cemc_iphi*EMCalTrk_cor_cemc_iphi + EMCalTrk_cor_cemc_ieta*EMCalTrk_cor_cemc_ieta) < %f ) )",
315  R_HCALIN));
316  T->SetAlias("EMCalTrk_hcalin_ntower",
317  Form(
318  "1*Sum$( ( sqrt(EMCalTrk_cor_hcalin_iphi*EMCalTrk_cor_hcalin_iphi + EMCalTrk_cor_hcalin_ieta*EMCalTrk_cor_hcalin_ieta) < %f ) )",
319  R_HCALIN));
320 
321  SaveCanvas(c1,
322  TString(_file0->GetName()) + TString("_DrawEcal_pDST_")
323  + Form("RCEMC%.1f_RCEMC%.1f_", R_CEMC, R_HCALIN)
324  + TString(c1->GetName()), kFALSE);
325 
326  cout << "///////////////////////////////////////////////" << endl;
327  cout << "// Projection center based on " << _file0->GetName() << endl;
328  cout << "///////////////////////////////////////////////" << endl;
329 
330  cout << "const double center_cemc_iphi = " << center_cemc_iphi << "; // "
331  << _file0->GetName() << endl;
332  cout << "const double center_cemc_ieta = " << center_cemc_ieta << "; // "
333  << _file0->GetName() << endl;
334  cout << "const double center_hcalin_iphi = " << center_hcalin_iphi << "; // "
335  << _file0->GetName() << endl;
336  cout << "const double center_hcalin_ieta = " << center_hcalin_ieta << "; // "
337  << _file0->GetName() << endl;
338 
339 }
340 
341 void
342 Edep_Checks(TString infile, const double R_CEMC, const double R_HCALIN,
343  TCut good_track_cut)
344 {
345  MakeRadiusCut(infile, R_CEMC, R_HCALIN, good_track_cut);
346 
347  double N_Event = T->GetEntries();
348 
349  TCanvas *c1 = new TCanvas("Edep_Checks" + cuts, "Edep_Checks" + cuts, 1900,
350  950);
351  c1->Divide(2, 2);
352  int idx = 1;
353  TPad * p;
354 
355  p = (TPad *) c1->cd(idx++);
356  c1->Update();
357  p->SetLogy();
358  T->Draw("EMCalTrk_cemc_ntower>>hEMCalTrk_cemc_ntower(16,-.5,15.5)",
359  good_track_cut);
360  hEMCalTrk_cemc_ntower->SetTitle(
361  Form("CEMC Cluster Size (R = %.1f);Cluster Size (Towers);Probability",
362  R_CEMC));
363  hEMCalTrk_cemc_ntower->Scale(1. / N_Event);
364 
365  p = (TPad *) c1->cd(idx++);
366  c1->Update();
367  p->SetLogy();
368  T->Draw("EMCalTrk_hcalin_ntower>>hEMCalTrk_hcalin_ntower(16,-.5,15.5)",
369  good_track_cut);
370  hEMCalTrk_hcalin_ntower->SetTitle(
371  Form(
372  "HCal_{in} Cluster Size (R = %.1f);Cluster Size (Towers);Probability",
373  R_HCALIN));
374  hEMCalTrk_hcalin_ntower->Scale(1. / N_Event);
375 
376  p = (TPad *) c1->cd(idx++);
377  c1->Update();
378  p->SetLogy();
379  T->Draw("EMCalTrk_cemc_e>>hEMCalTrk_cemc_e(240,-.0,12)", good_track_cut);
380  hEMCalTrk_cemc_e->SetTitle(
381  Form(
382  "CEMC Cluster Energy (R = %.1f);Cluster Energy (/bin);Probability/bin",
383  R_CEMC));
384  hEMCalTrk_cemc_e->Scale(1. / N_Event);
385 
386  p = (TPad *) c1->cd(idx++);
387  c1->Update();
388  p->SetLogy();
389  T->Draw("EMCalTrk_hcalin_e>>hEMCalTrk_hcalin_e(240,-.0,12)", good_track_cut);
390  hEMCalTrk_hcalin_e->SetTitle(
391  Form(
392  "HCal_{in} Cluster Energy (R = %.1f);Cluster Energy (GeV);Probability/bin",
393  R_HCALIN));
394  hEMCalTrk_hcalin_e->Scale(1. / N_Event);
395 
396  SaveCanvas(c1,
397  TString(_file0->GetName()) + TString("_DrawEcal_pDST_")
398  + Form("RCEMC%.1f_RCEMC%.1f_", R_CEMC, R_HCALIN)
399  + TString(c1->GetName()), kFALSE);
400 
401  TCanvas *c1 = new TCanvas("Edep_Checks_2D" + cuts, "Edep_Checks_2D" + cuts,
402  900, 900);
403 // c1->Divide(2, 2);
404 // int idx = 1;
405 // TPad * p;
406 
407  p = (TPad *) c1->cd(idx++);
408  c1->Update();
409  p->SetLogz();
410  T->Draw(
411  "EMCalTrk_hcalin_e:EMCalTrk_cemc_e>>h2_EMCalTrk_hcalin_e_EMCalTrk_cemc_e(240,-.0,12, 240,-.0,12)",
412  good_track_cut, "colz");
413  h2_EMCalTrk_hcalin_e_EMCalTrk_cemc_e->SetTitle(
414  Form(
415  "Energy distribution;CEMC Cluster Energy (R = %.1f) in GeV;HCal_{in} Cluster Energy (R = %.1f) in GeV",
416  R_CEMC, R_HCALIN));
417  h2_EMCalTrk_hcalin_e_EMCalTrk_cemc_e->Scale(1. / N_Event);
418  h2_EMCalTrk_hcalin_e_EMCalTrk_cemc_e->GetZaxis()->SetRangeUser(1e-6, 1);
419 
420  SaveCanvas(c1,
421  TString(_file0->GetName()) + TString("_DrawEcal_pDST_")
422  + Form("RCEMC%.1f_RCEMC%.1f_", R_CEMC, R_HCALIN)
423  + TString(c1->GetName()), kFALSE);
424 
425 }
426 
427 void
428 TrackProjection_Checks(TString infile, TCut good_track_cut)
429 {
430 
431  double N_Event = T->GetEntries();
432 
433  TCanvas *c1 = new TCanvas("TrackProjection_Checks_Raw" + cuts,
434  "TrackProjection_Checks_Raw" + cuts, 1800, 900);
435  c1->Divide(2, 1);
436  int idx = 1;
437  TPad * p;
438 
439  p = (TPad *) c1->cd(idx++);
440  c1->Update();
441  p->SetLogz();
442  T->Draw(
443  "DST.EMCalTrk.cemc_iphi:DST.EMCalTrk.cemc_ieta>>hcemc(130,-6.5,6.5,130,-6.5,6.5)",
444  "DST.EMCalTrk.cemc_energy", "colz");
445  hcemc->SetTitle(
446  "CEMC Tower Energy Distribution;Polar Distance (Tower Width);Azimuthal Distance (Tower Width)");
447  hcemc->Scale(1. / N_Event);
448  hcemc->GetZaxis()->SetRangeUser(1e-5, 30);
449 
450  p = (TPad *) c1->cd(idx++);
451  c1->Update();
452  p->SetLogz();
453  T->Draw(
454  "DST.EMCalTrk.hcalin_iphi:DST.EMCalTrk.hcalin_ieta>>hcalin(130,-6.5,6.5,130,-6.5,6.5)",
455  "DST.EMCalTrk.hcalin_energy", "colz");
456  hcalin->SetTitle(
457  "HCal_{in} Tower Energy Distribution;Polar Distance (Tower Width);Azimuthal Distance (Tower Width)");
458  hcalin->Scale(1. / N_Event);
459  hcalin->GetZaxis()->SetRangeUser(1e-5, 30);
460 
461  SaveCanvas(c1,
462  TString(_file0->GetName()) + TString("_DrawEcal_pDST_")
463  + TString(c1->GetName()), kFALSE);
464 
465  TCanvas *c1 = new TCanvas("TrackProjection_Checks_Tracking" + cuts,
466  "TrackProjection_Checks_Tracking" + cuts, 1800, 900);
467  c1->Divide(2, 1);
468  int idx = 1;
469  TPad * p;
470 
471  p = (TPad *) c1->cd(idx++);
472  c1->Update();
473  p->SetLogy();
474  T->Draw("EMCalTrk_pt/EMCalTrk_gpt>>htracking_pt(100,0,2)", "", "");
475  htracking_pt->SetTitle(
476  "Tracking reco precision;p_{T, reco}/p_{T, truth};Probability per bin");
477  htracking_pt->Scale(1. / N_Event);
478  T->Draw("EMCalTrk_pt/EMCalTrk_gpt>>htracking_pt_cut(100,0,2)", good_track_cut,
479  "same");
480  htracking_pt_cut->Scale(1. / N_Event);
481  htracking_pt_cut->SetLineColor(kRed);
482 
483  p = (TPad *) c1->cd(idx++);
484  c1->Update();
485  p->SetLogy();
486 
487  T->Draw("EMCalTrk_pt/EMCalTrk_gpt>0>>h_have_trak(3,-.5,2.5)", "", "");
488  h_have_trak->Scale(1. / N_Event);
489  T->Draw("2>>h_have_trak_cut(3,-.5,2.5)", good_track_cut, "same");
490  h_have_trak_cut->Scale(1. / N_Event);
491  h_have_trak_cut->SetLineColor(kRed);
492 
493  SaveCanvas(c1,
494  TString(_file0->GetName()) + TString("_DrawEcal_pDST_")
495  + TString(c1->GetName()), kFALSE);
496 
497  TCanvas *c1 = new TCanvas("TrackProjection_Checks_Cut" + cuts,
498  "TrackProjection_Checks_Cut" + cuts, 1800, 900);
499  c1->Divide(2, 1);
500  int idx = 1;
501  TPad * p;
502 
503  p = (TPad *) c1->cd(idx++);
504  c1->Update();
505  p->SetLogz();
506  T->Draw(
507  "DST.EMCalTrk.cemc_iphi:DST.EMCalTrk.cemc_ieta>>hcemcc(130,-6.5,6.5,130,-6.5,6.5)",
508  TString("DST.EMCalTrk.cemc_energy * (")
509  + TString(good_track_cut.GetTitle()) + ")", "colz");
510  hcemcc->SetTitle(
511  "CEMC Tower Energy Distribution;Polar Distance (Tower Width);Azimuthal Distance (Tower Width)");
512  hcemcc->Scale(1. / N_Event);
513  hcemcc->GetZaxis()->SetRangeUser(1e-5, 30);
514 
515  p = (TPad *) c1->cd(idx++);
516  c1->Update();
517  p->SetLogz();
518  T->Draw(
519  "DST.EMCalTrk.hcalin_iphi:DST.EMCalTrk.hcalin_ieta>>hcalinc(130,-6.5,6.5,130,-6.5,6.5)",
520  TString("DST.EMCalTrk.hcalin_energy * (")
521  + TString(good_track_cut.GetTitle()) + ")", "colz");
522  hcalinc->SetTitle(
523  "HCal_{in} Tower Energy Distribution;Polar Distance (Tower Width);Azimuthal Distance (Tower Width)");
524  hcalinc->Scale(1. / N_Event);
525  hcalinc->GetZaxis()->SetRangeUser(1e-5, 30);
526 
527  SaveCanvas(c1,
528  TString(_file0->GetName()) + TString("_DrawEcal_pDST_")
529  + TString(c1->GetName()), kFALSE);
530 
531  TCanvas *c1 = new TCanvas("TrackProjection_Checks_Proj" + cuts,
532  "TrackProjection_Checks_Proj" + cuts, 1800, 900);
533  c1->Divide(2, 2);
534  int idx = 1;
535  TPad * p;
536 
537  p = (TPad *) c1->cd(idx++);
538  c1->Update();
539  p->SetLogy();
540 
541  TH1 * proj = hcemc->ProjectionX();
542  proj->SetLineColor(kBlack);
543  proj->Draw("");
544  TH1 * proj = hcemcc->ProjectionX();
545  proj->SetLineColor(kRed);
546  proj->Draw("same");
547 
548  p = (TPad *) c1->cd(idx++);
549  c1->Update();
550  p->SetLogy();
551 
552  TH1 * proj = hcalin->ProjectionX();
553  proj->SetLineColor(kBlack);
554  proj->Draw("");
555  TH1 * proj = hcalinc->ProjectionX();
556  proj->SetLineColor(kRed);
557  proj->Draw("same");
558 
559  p = (TPad *) c1->cd(idx++);
560  c1->Update();
561  p->SetLogy();
562 
563  TH1 * proj = hcemc->ProjectionY();
564  proj->SetLineColor(kBlack);
565  proj->Draw("");
566  TH1 * proj = hcemcc->ProjectionY();
567  proj->SetLineColor(kRed);
568  proj->Draw("same");
569 
570  p = (TPad *) c1->cd(idx++);
571  c1->Update();
572  p->SetLogy();
573 
574  TH1 * proj = hcalin->ProjectionY();
575  proj->SetLineColor(kBlack);
576  proj->Draw("");
577  TH1 * proj = hcalinc->ProjectionY();
578  proj->SetLineColor(kRed);
579  proj->Draw("same");
580 
581  SaveCanvas(c1,
582  TString(_file0->GetName()) + TString("_DrawEcal_pDST_")
583  + TString(c1->GetName()), kFALSE);
584 }
585 
586 void
587 UpsilonPair_Checks(TString infile)
588 {
589 
590  TCanvas *c1 = new TCanvas("UpsilonPair_Checks" + cuts,
591  "UpsilonPair_Checks" + cuts, 1800, 900);
592  c1->Divide(3, 2);
593  int idx = 1;
594  TPad * p;
595 
596  p = (TPad *) c1->cd(idx++);
597  c1->Update();
598  T->Draw("DST.UpsilonPair.gmass>>gmass(100,0,10)");
599  T->Draw("DST.UpsilonPair.gmass>>gmassc(100,0,10)",
600  "DST.UpsilonPair.good_upsilon", "same");
601  gmassc->SetLineColor(kRed);
602 
603  p = (TPad *) c1->cd(idx++);
604  c1->Update();
605  T->Draw("DST.UpsilonPair.mass>>mass(100,0,10)");
606  T->Draw("DST.UpsilonPair.mass>>massc(100,0,10)",
607  "DST.UpsilonPair.good_upsilon", "same");
608  massc->SetLineColor(kRed);
609 
610  p = (TPad *) c1->cd(idx++);
611  c1->Update();
612  T->Draw("UpsilonPair_trk_gpt>>UpsilonPair_trk_gpt(100,0,10)");
613 
614  p = (TPad *) c1->cd(idx++);
615  c1->Update();
616  T->Draw("UpsilonPair_trk_pt:UpsilonPair_trk_gpt>>pt2(100,0,10,100,0,10)", "",
617  "colz");
618 
619  p = (TPad *) c1->cd(idx++);
620  c1->Update();
621  T->Draw("UpsilonPair_trk_pt:UpsilonPair_trk_gpt>>pt2c(100,0,10,100,0,10)",
622  "DST.UpsilonPair.good_upsilon", "colz");
623 
624  SaveCanvas(c1,
625  TString(_file0->GetName()) + TString("_DrawEcal_pDST_")
626  + TString(c1->GetName()), kFALSE);
627 }
628 
629 TGraphErrors *
630 Distribution2Efficiency(TH1F * hCEMC3_Max)
631 {
632  double threshold[10000] =
633  { 0 };
634  double eff[10000] =
635  { 0 };
636  double eff_err[10000] =
637  { 0 };
638 
639  assert(hCEMC3_Max->GetNbinsX()<10000);
640 
641  const double n = hCEMC3_Max->GetSum();
642  double pass = 0;
643  int cnt = 0;
644  for (int i = hCEMC3_Max->GetNbinsX(); i >= 1; i--)
645  {
646  pass += hCEMC3_Max->GetBinContent(i);
647 
648  const double pp = pass / n;
649 // const double z = 1.96;
650  const double z = 1.;
651 
652  const double A = z * sqrt(1. / n * pp * (1 - pp) + z * z / 4 / n / n);
653  const double B = 1 / (1 + z * z / n);
654 
655  threshold[cnt] = hCEMC3_Max->GetBinCenter(i);
656  eff[cnt] = (pp + z * z / 2 / n) * B;
657  eff_err[cnt] = A * B;
658 
659  cout << threshold[cnt] << ": " << "CL " << eff[cnt] << "+/-"
660  << eff_err[cnt] << endl;
661  cnt++;
662  }
663  TGraphErrors * ge = new TGraphErrors(cnt, threshold, eff, NULL, eff_err);
664 
665  return ge;
666 }
667 
668 void
669 TrackProjection_Checks_Comparison()
670 {
671 
672  SetOKStyle();
673  gStyle->SetOptStat(0);
674  gStyle->SetOptFit(1111);
675  TVirtualFitter::SetDefaultFitter("Minuit2");
676 
677  gSystem->Load("libg4eval.so");
678  gSystem->Load("libemcal_ana.so");
679  gSystem->Load("libg4vertex.so");
680 
681 
682  TFile * f =
683  new TFile(
684  "../../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");
685 
686  assert(f->IsOpen());
687  TH2F * hcemcc_2d = (TH2F *) f->GetObjectChecked("hcemcc", "TH2F");
688  assert(hcemcc_2d);
689 
690  TFile * f =
691  new TFile(
692  "../../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");
693 
694  assert(f->IsOpen());
695  TH2F * hcemcc_1d = (TH2F *) f->GetObjectChecked("hcemcc", "TH2F");
696  assert(hcemcc_1d);
697 
698  TCanvas *c1 = new TCanvas("TrackProjection_Checks_Comparison" ,
699  "TrackProjection_Checks_Comparison" , 1800, 1000);
700  c1->Divide(2, 1);
701  int idx = 1;
702  TPad * p;
703 
704  p = (TPad *) c1->cd(idx++);
705  c1->Update();
706  p->SetLogz();
707 
708  hcemcc_2d -> Draw("colz");
709  hcemcc_2d->SetTitle(";Polar Distance (Tower Width);Azimuthal Distance (Tower Width)");
710 
711  TLatex *text = new TLatex(0,7,"2-D Projective SPACAL, 0.9 < #eta_{e^{-}} < 1.0, E_{e^{-}} = 8 GeV");
712  text->SetTextAlign(22);
713  text->SetTextSize(0.044);
714  text->Draw();
715 
716 
717  p = (TPad *) c1->cd(idx++);
718  c1->Update();
719  p->SetLogz();
720 
721  hcemcc_1d -> Draw("colz");
722  hcemcc_1d->SetTitle(";Polar Distance (Tower Width);Azimuthal Distance (Tower Width)");
723 
724  TLatex *text = new TLatex(0,7,"1-D Projective SPACAL, 0.9 < #eta_{e^{-}} < 1.0, E_{e^{-}} = 8 GeV");
725  text->SetTextAlign(22);
726  text->SetTextSize(0.044);
727  text->Draw();
728 
729 
730  SaveCanvas(c1,
731  "../../sPHENIX_work/production_analysis/" + TString("DrawEcal_pDST_")
732  + TString(c1->GetName()));
733 }
734 
735 void
736 Logs()
737 {
738 
740  // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0_8GeV-ALL.root_Ana.root
742  const double center_cemc_iphi = 0.0202159; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0_8GeV-ALL.root_Ana.root
743  const double center_cemc_ieta = 0.00453938; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0_8GeV-ALL.root_Ana.root
744  const double center_hcalin_iphi = 0.488432; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0_8GeV-ALL.root_Ana.root
745  const double center_hcalin_ieta = -0.00936002; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0_8GeV-ALL.root_Ana.root
746 
748  // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0_8GeV-ALL.root_Ana.root
750  const double center_cemc_iphi = 0.0402477; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0_8GeV-ALL.root_Ana.root
751  const double center_cemc_ieta = 0.00574989; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0_8GeV-ALL.root_Ana.root
752  const double center_hcalin_iphi = 0.331081; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0_8GeV-ALL.root_Ana.root
753  const double center_hcalin_ieta = -0.000527382; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0_8GeV-ALL.root_Ana.root
754 
755 
757  // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_8GeV-ALL.root_Ana.root
759  const double center_cemc_iphi = 0.0649585; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_8GeV-ALL.root_Ana.root
760  const double center_cemc_ieta = 0.00465546; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_8GeV-ALL.root_Ana.root
761  const double center_hcalin_iphi = 0.552238; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_8GeV-ALL.root_Ana.root
762  const double center_hcalin_ieta = -0.0102162; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_8GeV-ALL.root_Ana.root
763 
765  // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_8GeV-ALL.root_Ana.root
767  const double center_cemc_iphi = 0.044375; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_8GeV-ALL.root_Ana.root
768  const double center_cemc_ieta = 0.00680547; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_8GeV-ALL.root_Ana.root
769  const double center_hcalin_iphi = 0.336741; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_8GeV-ALL.root_Ana.root
770  const double center_hcalin_ieta = 0.00108616; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_8GeV-ALL.root_Ana.root
771 
773  // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_4GeV-ALL.root_Ana.root
775  const double center_cemc_iphi = 0.0952638; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_4GeV-ALL.root_Ana.root
776  const double center_cemc_ieta = 0.00397331; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_4GeV-ALL.root_Ana.root
777  const double center_hcalin_iphi = 0.592338; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_4GeV-ALL.root_Ana.root
778  const double center_hcalin_ieta = -0.0118107; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_4GeV-ALL.root_Ana.root
779 
781  // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_4GeV-ALL.root_Ana.root
783  const double center_cemc_iphi = 0.0577477; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_4GeV-ALL.root_Ana.root
784  const double center_cemc_ieta = 0.00609392; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_4GeV-ALL.root_Ana.root
785  const double center_hcalin_iphi = 0.354109; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_4GeV-ALL.root_Ana.root
786  const double center_hcalin_ieta = -0.000273002; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_4GeV-ALL.root_Ana.root
787 
789  // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_2GeV-ALL.root_Ana.root
791  const double center_cemc_iphi = 0.170699; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_2GeV-ALL.root_Ana.root
792  const double center_cemc_ieta = 0.00345903; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_2GeV-ALL.root_Ana.root
793  const double center_hcalin_iphi = 0.668343; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_2GeV-ALL.root_Ana.root
794  const double center_hcalin_ieta = -0.0130487; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_2GeV-ALL.root_Ana.root
795 
797  // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_2GeV-ALL.root_Ana.root
799  const double center_cemc_iphi = 0.0898377; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_2GeV-ALL.root_Ana.root
800  const double center_cemc_ieta = 0.00445598; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_2GeV-ALL.root_Ana.root
801  const double center_hcalin_iphi = 0.376608; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_2GeV-ALL.root_Ana.root
802  const double center_hcalin_ieta = -0.0074984; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_2GeV-ALL.root_Ana.root
803 
805  // Projection center based on ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
807  const double center_cemc_iphi = 0.122055; // ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
808  const double center_cemc_ieta = 1.52588; // ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
809  const double center_hcalin_iphi = 0.644177; // ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
810  const double center_hcalin_ieta = -0.799621; // ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
812  // Projection center based on ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
814  const double center_cemc_iphi = 0.122055; // ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
815  const double center_cemc_ieta = 1.52588; // ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
816  const double center_hcalin_iphi = 0.644177; // ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
817  const double center_hcalin_ieta = -0.799621; // ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
819  // Projection center based on ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root
821  const double center_cemc_iphi = 0.251162; // ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root
822  const double center_cemc_ieta = 1.35067; // ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root
823  const double center_hcalin_iphi = 0.761301; // ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root
824  const double center_hcalin_ieta = -0.874644; // ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root
826  // Projection center based on ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root
828  const double center_cemc_iphi = 0.177891; // ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root
829  const double center_cemc_ieta = 2.62604; // ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root
830  const double center_hcalin_iphi = 0.393433; // ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root
831  const double center_hcalin_ieta = -0.116632; // ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root
833  // Projection center based on ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root
835  const double center_cemc_iphi = 0.53236; // ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root
836  const double center_cemc_ieta = 1.15917; // ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root
837  const double center_hcalin_iphi = 0.899527; // ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root
838  const double center_hcalin_ieta = -0.889783; // ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root
840  // Projection center based on ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root
842  const double center_cemc_iphi = 0.378081; // ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root
843  const double center_cemc_ieta = 2.47524; // ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root
844  const double center_hcalin_iphi = 0.445272; // ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root
845  const double center_hcalin_ieta = -0.180098; // ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root
846 
848  // Projection center based on ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root
850  const double center_cemc_iphi = 0.129035; // ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root
851  const double center_cemc_ieta = 0.0222246; // ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root
852  const double center_hcalin_iphi = 0.701734; // ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root
853  const double center_hcalin_ieta = -0.474402; // ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root
854 
856  // Projection center based on ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root
858  const double center_cemc_iphi = 0.0736277; // ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root
859  const double center_cemc_ieta = 0.0556152; // ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root
860  const double center_hcalin_iphi = 0.385379; // ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root
861  const double center_hcalin_ieta = -0.0596968; // ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root
862 
864  // Projection center based on ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root
866  const double center_cemc_iphi = 0.249273; // ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root
867  const double center_cemc_ieta = 0.0146006; // ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root
868  const double center_hcalin_iphi = 0.843407; // ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root
869  const double center_hcalin_ieta = -0.592856; // ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root
870 
872  // Projection center based on ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root
874  const double center_cemc_iphi = 0.12442; // ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root
875  const double center_cemc_ieta = 0.0514677; // ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root
876  const double center_hcalin_iphi = 0.447927; // ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root
877  const double center_hcalin_ieta = -0.101503; // ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root
878 
880  // Projection center based on ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
882  const double center_cemc_iphi = 0.0829824; // ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
883  const double center_cemc_ieta = 0.0275653; // ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
884  const double center_hcalin_iphi = 0.635325; // ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
885  const double center_hcalin_ieta = -0.408261; // ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root
886 
888  // Projection center based on ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_8GeV-ALL.root_Ana.root
890  const double center_cemc_iphi = 0.054392; // ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_8GeV-ALL.root_Ana.root
891  const double center_cemc_ieta = 0.0852583; // ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_8GeV-ALL.root_Ana.root
892  const double center_hcalin_iphi = 0.348779; // ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_8GeV-ALL.root_Ana.root
893  const double center_hcalin_ieta = -0.0348952; // ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_8GeV-ALL.root_Ana.root
894 
895 
896 
897 
899  // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_2GeV-ALL.root_Ana.root
901  const double center_cemc_iphi = -0.16237; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_2GeV-ALL.root_Ana.root
902  const double center_cemc_ieta = 0.00993428; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_2GeV-ALL.root_Ana.root
903  const double center_hcalin_iphi = 0.312564; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_2GeV-ALL.root_Ana.root
904  const double center_hcalin_ieta = -0.568171; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_2GeV-ALL.root_Ana.root
905 
907  // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_2GeV-ALL.root_Ana.root
909  const double center_cemc_iphi = -0.039828; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_2GeV-ALL.root_Ana.root
910  const double center_cemc_ieta = 0.0489088; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_2GeV-ALL.root_Ana.root
911  const double center_hcalin_iphi = 0.338006; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_2GeV-ALL.root_Ana.root
912  const double center_hcalin_ieta = -0.0939571; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_2GeV-ALL.root_Ana.root
913 
915  // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_4GeV-ALL.root_Ana.root
917  const double center_cemc_iphi = -0.0455698; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_4GeV-ALL.root_Ana.root
918  const double center_cemc_ieta = 0.0181539; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_4GeV-ALL.root_Ana.root
919  const double center_hcalin_iphi = 0.444123; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_4GeV-ALL.root_Ana.root
920  const double center_hcalin_ieta = -0.475483; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_4GeV-ALL.root_Ana.root
921 
923  // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_4GeV-ALL.root_Ana.root
925  const double center_cemc_iphi = 0.0114778; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_4GeV-ALL.root_Ana.root
926  const double center_cemc_ieta = 0.0580504; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_4GeV-ALL.root_Ana.root
927  const double center_hcalin_iphi = 0.340425; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_4GeV-ALL.root_Ana.root
928  const double center_hcalin_ieta = -0.0670486; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_4GeV-ALL.root_Ana.root
929 
930 
932  // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_8GeV-ALL.root_Ana.root
934  const double center_cemc_iphi = 0.00278605; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_8GeV-ALL.root_Ana.root
935  const double center_cemc_ieta = 0.0237387; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_8GeV-ALL.root_Ana.root
936  const double center_hcalin_iphi = 0.504268; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_8GeV-ALL.root_Ana.root
937  const double center_hcalin_ieta = -0.401039; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_8GeV-ALL.root_Ana.root
938 
940  // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_8GeV-ALL.root_Ana.root
942  const double center_cemc_iphi = 0.0313807; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_8GeV-ALL.root_Ana.root
943  const double center_cemc_ieta = 0.0800593; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_8GeV-ALL.root_Ana.root
944  const double center_hcalin_iphi = 0.336144; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_8GeV-ALL.root_Ana.root
945  const double center_hcalin_ieta = -0.0484305; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_8GeV-ALL.root_Ana.root
946 
947 
948 
949 
951  // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_8GeV-ALL.root_Ana.root
953  const double center_cemc_iphi = -0.122248; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_8GeV-ALL.root_Ana.root
954  const double center_cemc_ieta = 1.52555; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_8GeV-ALL.root_Ana.root
955  const double center_hcalin_iphi = 0.500768; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_8GeV-ALL.root_Ana.root
956  const double center_hcalin_ieta = -0.802395; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_8GeV-ALL.root_Ana.root
958  // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_8GeV-ALL.root_Ana.root
960  const double center_cemc_iphi = -0.0819904; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_8GeV-ALL.root_Ana.root
961  const double center_cemc_ieta = 2.83299; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_8GeV-ALL.root_Ana.root
962  const double center_hcalin_iphi = 0.33193; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_8GeV-ALL.root_Ana.root
963  const double center_hcalin_ieta = -0.0749197; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_8GeV-ALL.root_Ana.root
964 
966  // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_4GeV-ALL.root_Ana.root
968  const double center_cemc_iphi = -0.252243; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_4GeV-ALL.root_Ana.root
969  const double center_cemc_ieta = 1.34503; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_4GeV-ALL.root_Ana.root
970  const double center_hcalin_iphi = 0.411387; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_4GeV-ALL.root_Ana.root
971  const double center_hcalin_ieta = -0.860032; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_4GeV-ALL.root_Ana.root
972 
974  // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_4GeV-ALL.root_Ana.root
976  const double center_cemc_iphi = -0.171931; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_4GeV-ALL.root_Ana.root
977  const double center_cemc_ieta = 2.65013; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_4GeV-ALL.root_Ana.root
978  const double center_hcalin_iphi = 0.334954; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_4GeV-ALL.root_Ana.root
979  const double center_hcalin_ieta = -0.0992204; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_4GeV-ALL.root_Ana.root
980 
982  // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_2GeV-ALL.root_Ana.root
984  const double center_cemc_iphi = -0.530751; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_2GeV-ALL.root_Ana.root
985  const double center_cemc_ieta = 1.15259; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_2GeV-ALL.root_Ana.root
986  const double center_hcalin_iphi = 0.228558; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_2GeV-ALL.root_Ana.root
987  const double center_hcalin_ieta = -0.923731; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_2GeV-ALL.root_Ana.root
988 
990  // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_2GeV-ALL.root_Ana.root
992  const double center_cemc_iphi = -0.371817; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_2GeV-ALL.root_Ana.root
993  const double center_cemc_ieta = 2.46165; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_2GeV-ALL.root_Ana.root
994  const double center_hcalin_iphi = 0.317738; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_2GeV-ALL.root_Ana.root
995  const double center_hcalin_ieta = -0.119789; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_2GeV-ALL.root_Ana.root
996 
997 
998 
999 }