Back to home page

sPhenix code displayed by LXR

 
 

    


File indexing completed on 2025-08-05 08:12:16

0001 // $Id: $
0002 
0003 /*!
0004  * \file Draw_PHG4DSTReader.C
0005  * \brief
0006  * \author Jin Huang <jhuang@bnl.gov>
0007  * \version $Revision:   $
0008  * \date $Date: $
0009  */
0010 
0011 #include <cmath>
0012 #include <TFile.h>
0013 #include <TString.h>
0014 #include <TLine.h>
0015 #include <TTree.h>
0016 #include <cassert>
0017 #include "SaveCanvas.C"
0018 #include "SetOKStyle.C"
0019 using namespace std;
0020 
0021 TFile * _file0 = NULL;
0022 TTree * T = NULL;
0023 TString cuts = "";
0024 
0025 void
0026 DrawEcal_pDST( //
0027             const TString infile = "G4Hits_sPHENIX_gamma_eta0_8GeV-0000.root_Ana.root"//
0028     )
0029 {
0030   SetOKStyle();
0031   gStyle->SetOptStat(0);
0032   gStyle->SetOptFit(1111);
0033   TVirtualFitter::SetDefaultFitter("Minuit2");
0034 
0035   gSystem->Load("libg4eval.so");
0036   gSystem->Load("libemcal_ana.so");
0037   gSystem->Load("libg4vertex.so");
0038 
0039   if (!_file0)
0040     {
0041       TString chian_str = infile;
0042       chian_str.ReplaceAll("ALL", "*");
0043       chian_str.ReplaceAll("+", "\\+");
0044 
0045       TChain * t = new TChain("T");
0046       const int n = t->Add(chian_str);
0047 
0048       cout << "Loaded " << n << " root files with " << chian_str << endl;
0049       assert(n>0);
0050 
0051       T = t;
0052 
0053       _file0 = new TFile;
0054       _file0->SetName(infile);
0055 
0056       fstream flst(infile + ".lst", ios_base::out);
0057 
0058       TObjArray *fileElements = t->GetListOfFiles();
0059       TIter next(fileElements);
0060       TChainElement *chEl = 0;
0061       while ((chEl = (TChainElement*) next()))
0062         {
0063           flst << chEl->GetTitle() << endl;
0064         }
0065       flst.close();
0066 
0067       cout << "Saved file list to " << infile + ".lst" << endl;
0068     }
0069 
0070   assert(_file0);
0071 
0072   T->SetAlias("UpsilonPair_trk_gpt",
0073       "1*sqrt(DST.UpsilonPair.trk.gpx**2 + DST.UpsilonPair.trk.gpy**2)");
0074   T->SetAlias("UpsilonPair_trk_pt",
0075       "1*sqrt(DST.UpsilonPair.trk.px**2 + DST.UpsilonPair.trk.py**2)");
0076 
0077   T->SetAlias("MCPhoton_pt", "1*sqrt(DST.MCPhoton.px**2 + DST.MCPhoton.py**2)");
0078   T->SetAlias("MCPhoton_gpt",
0079       "1*sqrt(DST.MCPhoton.gpx**2 + DST.MCPhoton.gpy**2)");
0080 
0081   const TCut event_sel = "1*1";
0082   cuts = "_all_event";
0083 
0084   cout << "Build event selection of " << (const char *) event_sel << endl;
0085 
0086   T->Draw(">>EventList", event_sel);
0087   TEventList * elist = gDirectory->GetObjectChecked("EventList", "TEventList");
0088   cout << elist->GetN() << " / " << T->GetEntriesFast() << " events selected"
0089       << endl;
0090 
0091   T->SetEventList(elist);
0092 ////
0093 //  UpsilonPair_Checks(infile);
0094 
0095   TrackProjection_Checks(infile, " ");
0096 //  Edep_Checks(infile, 1.4, 1.4, " ");
0097 
0098 }
0099 
0100 void
0101 MakeRadiusCut(TString infile, const double R_CEMC, const double R_HCALIN,
0102     TCut good_track_cut)
0103 {
0104   TCanvas *c1 = new TCanvas("MakeRadiusCut" + cuts, "MakeRadiusCut" + cuts,
0105       1800, 900);
0106   c1->Divide(4, 2);
0107   int idx = 1;
0108   TPad * p;
0109 
0110   p = (TPad *) c1->cd(idx++);
0111   c1->Update();
0112   p->SetLogy();
0113 
0114   T->Draw("DST.MCPhoton.cemc_iphi>>hcemc_iphi(130,-6.5,6.5)",
0115       TString("DST.MCPhoton.cemc_energy * (")
0116           + TString(good_track_cut.GetTitle()) + ")");
0117   hcemc_iphi->SetTitle(
0118       ";CEMC Azimuthal Distance (Tower Width);Energy Distribution");
0119   hcemc_iphi->Fit("gaus", "M");
0120   TF1* f = (TF1*) (hcemc_iphi->GetListOfFunctions()->At(0));
0121   assert(f);
0122   T->SetAlias("MCPhoton_cor_cemc_iphi",
0123       Form("DST.MCPhoton.cemc_iphi - %f", f->GetParameter(1)));
0124 
0125   const double center_cemc_iphi = f->GetParameter(1);
0126 
0127   p = (TPad *) c1->cd(idx++);
0128   c1->Update();
0129   p->SetLogy();
0130 
0131   T->Draw("abs(MCPhoton_cor_cemc_iphi)>>hMCPhoton_cor_cemc_iphi(130,0,6.5)",
0132       TString("DST.MCPhoton.cemc_energy * (")
0133           + TString(good_track_cut.GetTitle()) + ")");
0134 
0135   p->DrawFrame(0, 1e-3, 6.5, 1,
0136       ";CEMC Azimuthal Distance (Tower Width);Energy Leakage Ratio");
0137   TGraphErrors * ge = Distribution2Efficiency(hMCPhoton_cor_cemc_iphi);
0138   ge->SetLineColor(kBlue + 2);
0139   ge->SetMarkerColor(kBlue + 21);
0140   ge->SetMarkerColor(kFullCircle);
0141   ge->SetLineWidth(3);
0142   ge->Draw("lp");
0143 
0144   p = (TPad *) c1->cd(idx++);
0145   c1->Update();
0146   p->SetLogy();
0147 
0148   T->Draw("DST.MCPhoton.cemc_ieta>>hcemc_ieta(130,-6.5,6.5)",
0149       TString("DST.MCPhoton.cemc_energy * (")
0150           + TString(good_track_cut.GetTitle()) + ")");
0151   hcemc_ieta->SetTitle(
0152       ";CEMC Polar Distance (Tower Width);Energy Distribution");
0153   hcemc_ieta->Fit("gaus", "M");
0154   TF1* f = (TF1*) (hcemc_ieta->GetListOfFunctions()->At(0));
0155   assert(f);
0156   T->SetAlias("MCPhoton_cor_cemc_ieta",
0157       Form("DST.MCPhoton.cemc_ieta - %f", f->GetParameter(1)));
0158 
0159   const double center_cemc_ieta = f->GetParameter(1);
0160 
0161   p = (TPad *) c1->cd(idx++);
0162   c1->Update();
0163   p->SetLogy();
0164 
0165   T->Draw("abs(MCPhoton_cor_cemc_ieta)>>hMCPhoton_cor_cemc_ieta(130,0,6.5)",
0166       TString("DST.MCPhoton.cemc_energy * (")
0167           + TString(good_track_cut.GetTitle()) + ")");
0168 
0169   p->DrawFrame(0, 1e-3, 6.5, 1,
0170       ";CEMC Polar Distance (Tower Width);Energy Leakage Ratio");
0171   TGraphErrors * ge = Distribution2Efficiency(hMCPhoton_cor_cemc_ieta);
0172   ge->SetLineColor(kBlue + 2);
0173   ge->SetMarkerColor(kBlue + 21);
0174   ge->SetMarkerColor(kFullCircle);
0175   ge->SetLineWidth(3);
0176   ge->Draw("lp");
0177 
0178   p = (TPad *) c1->cd(idx++);
0179   c1->Update();
0180   p->SetLogy();
0181 
0182   T->Draw("DST.MCPhoton.hcalin_iphi>>hhcalin_iphi(130,-6.5,6.5)",
0183       TString("DST.MCPhoton.hcalin_energy * (")
0184           + TString(good_track_cut.GetTitle()) + ")");
0185   hhcalin_iphi->SetTitle(
0186       ";HCal_{IN} Azimuthal Distance (Tower Width);Energy Distribution");
0187   hhcalin_iphi->Fit("gaus", "M");
0188   TF1* f = (TF1*) (hhcalin_iphi->GetListOfFunctions()->At(0));
0189   assert(f);
0190   T->SetAlias("MCPhoton_cor_hcalin_iphi",
0191       Form("DST.MCPhoton.hcalin_iphi - %f", f->GetParameter(1)));
0192 
0193   const double center_hcalin_iphi = f->GetParameter(1);
0194 
0195   p = (TPad *) c1->cd(idx++);
0196   c1->Update();
0197   p->SetLogy();
0198 
0199   T->Draw("abs(MCPhoton_cor_hcalin_iphi)>>hMCPhoton_cor_hcalin_iphi(130,0,6.5)",
0200       TString("DST.MCPhoton.hcalin_energy * (")
0201           + TString(good_track_cut.GetTitle()) + ")");
0202 
0203   p->DrawFrame(0, 1e-3, 6.5, 1,
0204       ";HCal_{IN} Azimuthal Distance (Tower Width);Energy Leakage Ratio");
0205   TGraphErrors * ge = Distribution2Efficiency(hMCPhoton_cor_hcalin_iphi);
0206   ge->SetLineColor(kBlue + 2);
0207   ge->SetMarkerColor(kBlue + 21);
0208   ge->SetMarkerColor(kFullCircle);
0209   ge->SetLineWidth(3);
0210   ge->Draw("lp");
0211 
0212   p = (TPad *) c1->cd(idx++);
0213   c1->Update();
0214   p->SetLogy();
0215 
0216   T->Draw("DST.MCPhoton.hcalin_ieta>>hhcalin_ieta(130,-6.5,6.5)",
0217       TString("DST.MCPhoton.hcalin_energy * (")
0218           + TString(good_track_cut.GetTitle()) + ")");
0219   hhcalin_ieta->SetTitle(
0220       ";HCal_{IN} Polar Distance (Tower Width);Energy Distribution");
0221   hhcalin_ieta->Fit("gaus", "M");
0222   TF1* f = (TF1*) (hhcalin_ieta->GetListOfFunctions()->At(0));
0223   assert(f);
0224   T->SetAlias("MCPhoton_cor_hcalin_ieta",
0225       Form("DST.MCPhoton.hcalin_ieta - %f", f->GetParameter(1)));
0226 
0227   const double center_hcalin_ieta = f->GetParameter(1);
0228 
0229   p = (TPad *) c1->cd(idx++);
0230   c1->Update();
0231   p->SetLogy();
0232 
0233   T->Draw("abs(MCPhoton_cor_hcalin_ieta)>>hMCPhoton_cor_hcalin_ieta(130,0,6.5)",
0234       TString("DST.MCPhoton.hcalin_energy * (")
0235           + TString(good_track_cut.GetTitle()) + ")");
0236 
0237   p->DrawFrame(0, 1e-3, 6.5, 1,
0238       ";HCal_{IN} Polar Distance (Tower Width);Energy Leakage Ratio");
0239   TGraphErrors * ge = Distribution2Efficiency(hMCPhoton_cor_hcalin_ieta);
0240   ge->SetLineColor(kBlue + 2);
0241   ge->SetMarkerColor(kBlue + 21);
0242   ge->SetMarkerColor(kFullCircle);
0243   ge->SetLineWidth(3);
0244   ge->Draw("lp");
0245 
0246   T->SetAlias("MCPhoton_cemc_e",
0247       Form(
0248           "1*Sum$( DST.MCPhoton.cemc_energy * ( sqrt(MCPhoton_cor_cemc_iphi*MCPhoton_cor_cemc_iphi + MCPhoton_cor_cemc_ieta*MCPhoton_cor_cemc_ieta) < %f )   )",
0249           R_CEMC));
0250   T->SetAlias("MCPhoton_hcalin_e",
0251       Form(
0252           "1*Sum$( DST.MCPhoton.hcalin_energy * ( sqrt(MCPhoton_cor_hcalin_iphi*MCPhoton_cor_hcalin_iphi + MCPhoton_cor_hcalin_ieta*MCPhoton_cor_hcalin_ieta) < %f )   )",
0253           R_CEMC));
0254 
0255   T->SetAlias("MCPhoton_cemc_ntower",
0256       Form(
0257           "1*Sum$(( sqrt(MCPhoton_cor_cemc_iphi*MCPhoton_cor_cemc_iphi + MCPhoton_cor_cemc_ieta*MCPhoton_cor_cemc_ieta) < %f )   )",
0258           R_HCALIN));
0259   T->SetAlias("MCPhoton_hcalin_ntower",
0260       Form(
0261           "1*Sum$( ( sqrt(MCPhoton_cor_hcalin_iphi*MCPhoton_cor_hcalin_iphi + MCPhoton_cor_hcalin_ieta*MCPhoton_cor_hcalin_ieta) < %f )   )",
0262           R_HCALIN));
0263 
0264   SaveCanvas(c1,
0265       TString(_file0->GetName()) + TString("_DrawEcal_pDST_")
0266           + Form("RCEMC%.1f_RCEMC%.1f_", R_CEMC, R_HCALIN)
0267           + TString(c1->GetName()), kFALSE);
0268 
0269   cout << "///////////////////////////////////////////////" << endl;
0270   cout << "// Projection center based on " << _file0->GetName() << endl;
0271   cout << "///////////////////////////////////////////////" << endl;
0272 
0273   cout << "const double center_cemc_iphi = " << center_cemc_iphi << "; // "
0274       << _file0->GetName() << endl;
0275   cout << "const double center_cemc_ieta = " << center_cemc_ieta << "; // "
0276       << _file0->GetName() << endl;
0277   cout << "const double center_hcalin_iphi = " << center_hcalin_iphi << "; // "
0278       << _file0->GetName() << endl;
0279   cout << "const double center_hcalin_ieta = " << center_hcalin_ieta << "; // "
0280       << _file0->GetName() << endl;
0281 
0282 }
0283 
0284 void
0285 Edep_Checks(TString infile, const double R_CEMC, const double R_HCALIN,
0286     TCut good_track_cut)
0287 {
0288   MakeRadiusCut(infile, R_CEMC, R_HCALIN, good_track_cut);
0289 
0290   double N_Event = T->GetEntries();
0291 
0292   TCanvas *c1 = new TCanvas("Edep_Checks" + cuts, "Edep_Checks" + cuts, 1900,
0293       950);
0294   c1->Divide(2, 2);
0295   int idx = 1;
0296   TPad * p;
0297 
0298   p = (TPad *) c1->cd(idx++);
0299   c1->Update();
0300   p->SetLogy();
0301   T->Draw("MCPhoton_cemc_ntower>>hMCPhoton_cemc_ntower(16,-.5,15.5)",
0302       good_track_cut);
0303   hMCPhoton_cemc_ntower->SetTitle(
0304       Form("CEMC Cluster Size (R = %.1f);Cluster Size (Towers);Probability",
0305           R_CEMC));
0306   hMCPhoton_cemc_ntower->Scale(1. / N_Event);
0307 
0308   p = (TPad *) c1->cd(idx++);
0309   c1->Update();
0310   p->SetLogy();
0311   T->Draw("MCPhoton_hcalin_ntower>>hMCPhoton_hcalin_ntower(16,-.5,15.5)",
0312       good_track_cut);
0313   hMCPhoton_hcalin_ntower->SetTitle(
0314       Form(
0315           "HCal_{in} Cluster Size (R = %.1f);Cluster Size (Towers);Probability",
0316           R_HCALIN));
0317   hMCPhoton_hcalin_ntower->Scale(1. / N_Event);
0318 
0319   p = (TPad *) c1->cd(idx++);
0320   c1->Update();
0321   p->SetLogy();
0322   T->Draw("MCPhoton_cemc_e>>hMCPhoton_cemc_e(240,-.0,12)", good_track_cut);
0323   hMCPhoton_cemc_e->SetTitle(
0324       Form(
0325           "CEMC Cluster Energy (R = %.1f);Cluster Energy (/bin);Probability/bin",
0326           R_CEMC));
0327   hMCPhoton_cemc_e->Scale(1. / N_Event);
0328 
0329   p = (TPad *) c1->cd(idx++);
0330   c1->Update();
0331   p->SetLogy();
0332   T->Draw("MCPhoton_hcalin_e>>hMCPhoton_hcalin_e(240,-.0,12)", good_track_cut);
0333   hMCPhoton_hcalin_e->SetTitle(
0334       Form(
0335           "HCal_{in} Cluster Energy (R = %.1f);Cluster Energy (GeV);Probability/bin",
0336           R_HCALIN));
0337   hMCPhoton_hcalin_e->Scale(1. / N_Event);
0338 
0339   SaveCanvas(c1,
0340       TString(_file0->GetName()) + TString("_DrawEcal_pDST_")
0341           + Form("RCEMC%.1f_RCEMC%.1f_", R_CEMC, R_HCALIN)
0342           + TString(c1->GetName()), kFALSE);
0343 
0344   TCanvas *c1 = new TCanvas("Edep_Checks_2D" + cuts, "Edep_Checks_2D" + cuts,
0345       900, 900);
0346 //  c1->Divide(2, 2);
0347 //  int idx = 1;
0348 //  TPad * p;
0349 
0350   p = (TPad *) c1->cd(idx++);
0351   c1->Update();
0352   p->SetLogz();
0353   T->Draw(
0354       "MCPhoton_hcalin_e:MCPhoton_cemc_e>>h2_MCPhoton_hcalin_e_MCPhoton_cemc_e(240,-.0,12, 240,-.0,12)",
0355       good_track_cut, "colz");
0356   h2_MCPhoton_hcalin_e_MCPhoton_cemc_e->SetTitle(
0357       Form(
0358           "Energy distribution;CEMC Cluster Energy (R = %.1f) in GeV;HCal_{in} Cluster Energy (R = %.1f) in GeV",
0359           R_CEMC, R_HCALIN));
0360   h2_MCPhoton_hcalin_e_MCPhoton_cemc_e->Scale(1. / N_Event);
0361   h2_MCPhoton_hcalin_e_MCPhoton_cemc_e->GetZaxis()->SetRangeUser(1e-6, 1);
0362 
0363   SaveCanvas(c1,
0364       TString(_file0->GetName()) + TString("_DrawEcal_pDST_")
0365           + Form("RCEMC%.1f_RCEMC%.1f_", R_CEMC, R_HCALIN)
0366           + TString(c1->GetName()), kFALSE);
0367 
0368 }
0369 
0370 void
0371 TrackProjection_Checks(TString infile, TCut good_track_cut)
0372 {
0373 
0374   double N_Event = T->GetEntries();
0375 
0376   TCanvas *c1 = new TCanvas("TrackProjection_Checks_Raw" + cuts,
0377       "TrackProjection_Checks_Raw" + cuts, 1800, 900);
0378   c1->Divide(2, 1);
0379   int idx = 1;
0380   TPad * p;
0381 
0382   p = (TPad *) c1->cd(idx++);
0383   c1->Update();
0384   p->SetLogz();
0385   T->Draw(
0386       "DST.MCPhoton.cemc_iphi:DST.MCPhoton.cemc_ieta>>hcemc(130,-6.5,6.5,130,-6.5,6.5)",
0387       "DST.MCPhoton.cemc_energy", "colz");
0388   hcemc->SetTitle(
0389       "CEMC Tower Energy Distribution;Polar Distance (Tower Width);Azimuthal Distance (Tower Width)");
0390   hcemc->Scale(1. / N_Event);
0391   hcemc->GetZaxis()->SetRangeUser(1e-5, 30);
0392 
0393   p = (TPad *) c1->cd(idx++);
0394   c1->Update();
0395   p->SetLogz();
0396   T->Draw(
0397       "DST.MCPhoton.hcalin_iphi:DST.MCPhoton.hcalin_ieta>>hcalin(130,-6.5,6.5,130,-6.5,6.5)",
0398       "DST.MCPhoton.hcalin_energy", "colz");
0399   hcalin->SetTitle(
0400       "HCal_{in} Tower Energy Distribution;Polar Distance (Tower Width);Azimuthal Distance (Tower Width)");
0401   hcalin->Scale(1. / N_Event);
0402   hcalin->GetZaxis()->SetRangeUser(1e-5, 30);
0403 
0404   SaveCanvas(c1,
0405       TString(_file0->GetName()) + TString("_DrawEcal_pDST_")
0406           + TString(c1->GetName()), kFALSE);
0407 
0408 
0409 
0410   TCanvas *c1 = new TCanvas("TrackProjection_Checks_Proj" + cuts,
0411       "TrackProjection_Checks_Proj" + cuts, 1800, 900);
0412   c1->Divide(2, 2);
0413   int idx = 1;
0414   TPad * p;
0415 
0416   p = (TPad *) c1->cd(idx++);
0417   c1->Update();
0418   p->SetLogy();
0419 
0420   TH1 * proj = hcemc->ProjectionX();
0421   proj->SetLineColor(kBlack);
0422   proj->Draw("");
0423 //  TH1 * proj = hcemcc->ProjectionX();
0424 //  proj->SetLineColor(kRed);
0425 //  proj->Draw("same");
0426 
0427   p = (TPad *) c1->cd(idx++);
0428   c1->Update();
0429   p->SetLogy();
0430 
0431   TH1 * proj = hcalin->ProjectionX();
0432   proj->SetLineColor(kBlack);
0433   proj->Draw("");
0434 //  TH1 * proj = hcalinc->ProjectionX();
0435 //  proj->SetLineColor(kRed);
0436 //  proj->Draw("same");
0437 
0438   p = (TPad *) c1->cd(idx++);
0439   c1->Update();
0440   p->SetLogy();
0441 
0442   TH1 * proj = hcemc->ProjectionY();
0443   proj->SetLineColor(kBlack);
0444   proj->Draw("");
0445 //  TH1 * proj = hcemcc->ProjectionY();
0446 //  proj->SetLineColor(kRed);
0447 //  proj->Draw("same");
0448 
0449   p = (TPad *) c1->cd(idx++);
0450   c1->Update();
0451   p->SetLogy();
0452 
0453   TH1 * proj = hcalin->ProjectionY();
0454   proj->SetLineColor(kBlack);
0455   proj->Draw("");
0456 //  TH1 * proj = hcalinc->ProjectionY();
0457 //  proj->SetLineColor(kRed);
0458 //  proj->Draw("same");
0459 
0460   SaveCanvas(c1,
0461       TString(_file0->GetName()) + TString("_DrawEcal_pDST_")
0462           + TString(c1->GetName()), kFALSE);
0463 }
0464 
0465 void
0466 UpsilonPair_Checks(TString infile)
0467 {
0468 
0469   TCanvas *c1 = new TCanvas("UpsilonPair_Checks" + cuts,
0470       "UpsilonPair_Checks" + cuts, 1800, 900);
0471   c1->Divide(3, 2);
0472   int idx = 1;
0473   TPad * p;
0474 
0475   p = (TPad *) c1->cd(idx++);
0476   c1->Update();
0477   T->Draw("DST.UpsilonPair.gmass>>gmass(100,0,10)");
0478   T->Draw("DST.UpsilonPair.gmass>>gmassc(100,0,10)",
0479       "DST.UpsilonPair.good_upsilon", "same");
0480   gmassc->SetLineColor(kRed);
0481 
0482   p = (TPad *) c1->cd(idx++);
0483   c1->Update();
0484   T->Draw("DST.UpsilonPair.mass>>mass(100,0,10)");
0485   T->Draw("DST.UpsilonPair.mass>>massc(100,0,10)",
0486       "DST.UpsilonPair.good_upsilon", "same");
0487   massc->SetLineColor(kRed);
0488 
0489   p = (TPad *) c1->cd(idx++);
0490   c1->Update();
0491   T->Draw("UpsilonPair_trk_gpt>>UpsilonPair_trk_gpt(100,0,10)");
0492 
0493   p = (TPad *) c1->cd(idx++);
0494   c1->Update();
0495   T->Draw("UpsilonPair_trk_pt:UpsilonPair_trk_gpt>>pt2(100,0,10,100,0,10)", "",
0496       "colz");
0497 
0498   p = (TPad *) c1->cd(idx++);
0499   c1->Update();
0500   T->Draw("UpsilonPair_trk_pt:UpsilonPair_trk_gpt>>pt2c(100,0,10,100,0,10)",
0501       "DST.UpsilonPair.good_upsilon", "colz");
0502 
0503   SaveCanvas(c1,
0504       TString(_file0->GetName()) + TString("_DrawEcal_pDST_")
0505           + TString(c1->GetName()), kFALSE);
0506 }
0507 
0508 TGraphErrors *
0509 Distribution2Efficiency(TH1F * hCEMC3_Max)
0510 {
0511   double threshold[10000] =
0512     { 0 };
0513   double eff[10000] =
0514     { 0 };
0515   double eff_err[10000] =
0516     { 0 };
0517 
0518   assert(hCEMC3_Max->GetNbinsX()<10000);
0519 
0520   const double n = hCEMC3_Max->GetSum();
0521   double pass = 0;
0522   int cnt = 0;
0523   for (int i = hCEMC3_Max->GetNbinsX(); i >= 1; i--)
0524     {
0525       pass += hCEMC3_Max->GetBinContent(i);
0526 
0527       const double pp = pass / n;
0528 //      const double z = 1.96;
0529       const double z = 1.;
0530 
0531       const double A = z * sqrt(1. / n * pp * (1 - pp) + z * z / 4 / n / n);
0532       const double B = 1 / (1 + z * z / n);
0533 
0534       threshold[cnt] = hCEMC3_Max->GetBinCenter(i);
0535       eff[cnt] = (pp + z * z / 2 / n) * B;
0536       eff_err[cnt] = A * B;
0537 
0538       cout << threshold[cnt] << ": " << "CL " << eff[cnt] << "+/-"
0539           << eff_err[cnt] << endl;
0540       cnt++;
0541     }
0542   TGraphErrors * ge = new TGraphErrors(cnt, threshold, eff, NULL, eff_err);
0543 
0544   return ge;
0545 }
0546 
0547 void
0548 TrackProjection_Checks_Comparison()
0549 {
0550 
0551   SetOKStyle();
0552   gStyle->SetOptStat(0);
0553   gStyle->SetOptFit(1111);
0554   TVirtualFitter::SetDefaultFitter("Minuit2");
0555 
0556   gSystem->Load("libg4eval.so");
0557   gSystem->Load("libemcal_ana.so");
0558   gSystem->Load("libg4vertex.so");
0559 
0560 
0561   TFile * f =
0562       new TFile(
0563           "/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");
0564 
0565   assert(f->IsOpen());
0566   TH2F * hcemcc_2d = (TH2F *) f->GetObjectChecked("hcemcc", "TH2F");
0567   assert(hcemcc_2d);
0568 
0569   TFile * f =
0570       new TFile(
0571           "/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");
0572 
0573   assert(f->IsOpen());
0574   TH2F * hcemcc_1d = (TH2F *) f->GetObjectChecked("hcemcc", "TH2F");
0575   assert(hcemcc_1d);
0576 
0577   TCanvas *c1 = new TCanvas("TrackProjection_Checks_Comparison" ,
0578       "TrackProjection_Checks_Comparison" , 1800, 1000);
0579   c1->Divide(2, 1);
0580   int idx = 1;
0581   TPad * p;
0582 
0583   p = (TPad *) c1->cd(idx++);
0584   c1->Update();
0585   p->SetLogz();
0586 
0587   hcemcc_2d -> Draw("colz");
0588   hcemcc_2d->SetTitle(";Polar Distance (Tower Width);Azimuthal Distance (Tower Width)");
0589 
0590   TLatex *text = new TLatex(0,7,"2-D Projective SPACAL, 0.9 < #eta_{e^{-}} < 1.0, E_{e^{-}} = 8 GeV");
0591   text->SetTextAlign(22);
0592   text->SetTextSize(0.044);
0593   text->Draw();
0594 
0595 
0596   p = (TPad *) c1->cd(idx++);
0597   c1->Update();
0598   p->SetLogz();
0599 
0600   hcemcc_1d -> Draw("colz");
0601   hcemcc_1d->SetTitle(";Polar Distance (Tower Width);Azimuthal Distance (Tower Width)");
0602 
0603   TLatex *text = new TLatex(0,7,"1-D Projective SPACAL, 0.9 < #eta_{e^{-}} < 1.0, E_{e^{-}} = 8 GeV");
0604   text->SetTextAlign(22);
0605   text->SetTextSize(0.044);
0606   text->Draw();
0607 
0608 
0609   SaveCanvas(c1,
0610       "/phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/" + TString("DrawEcal_pDST_")
0611           + TString(c1->GetName()));
0612 }
0613 
0614 void
0615 Logs()
0616 {
0617 
0618   ///////////////////////////////////////////////
0619   // 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
0620   ///////////////////////////////////////////////
0621   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
0622   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
0623   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
0624   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
0625 
0626   ///////////////////////////////////////////////
0627   // 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
0628   ///////////////////////////////////////////////
0629   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
0630   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
0631   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
0632   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
0633 
0634 
0635   ///////////////////////////////////////////////
0636   // 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
0637   ///////////////////////////////////////////////
0638   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
0639   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
0640   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
0641   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
0642 
0643   ///////////////////////////////////////////////
0644   // 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
0645   ///////////////////////////////////////////////
0646   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
0647   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
0648   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
0649   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
0650 
0651   ///////////////////////////////////////////////
0652   // 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
0653   ///////////////////////////////////////////////
0654   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
0655   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
0656   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
0657   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
0658 
0659   ///////////////////////////////////////////////
0660   // 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
0661   ///////////////////////////////////////////////
0662   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
0663   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
0664   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
0665   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
0666 
0667   ///////////////////////////////////////////////
0668   // 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
0669   ///////////////////////////////////////////////
0670   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
0671   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
0672   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
0673   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
0674 
0675   ///////////////////////////////////////////////
0676   // 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
0677   ///////////////////////////////////////////////
0678   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
0679   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
0680   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
0681   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
0682 
0683   ///////////////////////////////////////////////
0684   // 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
0685   ///////////////////////////////////////////////
0686   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
0687   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
0688   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
0689   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
0690   ///////////////////////////////////////////////
0691   // 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
0692   ///////////////////////////////////////////////
0693   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
0694   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
0695   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
0696   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
0697   ///////////////////////////////////////////////
0698   // 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
0699   ///////////////////////////////////////////////
0700   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
0701   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
0702   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
0703   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
0704   ///////////////////////////////////////////////
0705   // 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
0706   ///////////////////////////////////////////////
0707   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
0708   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
0709   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
0710   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
0711   ///////////////////////////////////////////////
0712   // 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
0713   ///////////////////////////////////////////////
0714   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
0715   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
0716   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
0717   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
0718   ///////////////////////////////////////////////
0719   // 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
0720   ///////////////////////////////////////////////
0721   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
0722   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
0723   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
0724   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
0725 
0726   ///////////////////////////////////////////////
0727   // 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
0728   ///////////////////////////////////////////////
0729   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
0730   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
0731   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
0732   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
0733 
0734   ///////////////////////////////////////////////
0735   // 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
0736   ///////////////////////////////////////////////
0737   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
0738   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
0739   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
0740   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
0741 
0742   ///////////////////////////////////////////////
0743   // 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
0744   ///////////////////////////////////////////////
0745   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
0746   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
0747   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
0748   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
0749 
0750   ///////////////////////////////////////////////
0751   // 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
0752   ///////////////////////////////////////////////
0753   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
0754   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
0755   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
0756   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
0757 
0758   ///////////////////////////////////////////////
0759   // 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
0760   ///////////////////////////////////////////////
0761   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
0762   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
0763   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
0764   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
0765 
0766   ///////////////////////////////////////////////
0767   // 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
0768   ///////////////////////////////////////////////
0769   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
0770   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
0771   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
0772   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
0773 
0774 
0775 
0776 
0777   ///////////////////////////////////////////////
0778   // 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
0779   ///////////////////////////////////////////////
0780   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
0781   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
0782   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
0783   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
0784 
0785   ///////////////////////////////////////////////
0786   // 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
0787   ///////////////////////////////////////////////
0788   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
0789   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
0790   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
0791   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
0792 
0793   ///////////////////////////////////////////////
0794   // 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
0795   ///////////////////////////////////////////////
0796   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
0797   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
0798   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
0799   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
0800 
0801   ///////////////////////////////////////////////
0802   // 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
0803   ///////////////////////////////////////////////
0804   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
0805   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
0806   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
0807   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
0808 
0809 
0810   ///////////////////////////////////////////////
0811   // 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
0812   ///////////////////////////////////////////////
0813   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
0814   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
0815   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
0816   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
0817 
0818   ///////////////////////////////////////////////
0819   // 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
0820   ///////////////////////////////////////////////
0821   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
0822   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
0823   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
0824   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
0825 
0826 
0827 
0828 
0829   ///////////////////////////////////////////////
0830   // 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
0831   ///////////////////////////////////////////////
0832   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
0833   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
0834   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
0835   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
0836   ///////////////////////////////////////////////
0837   // 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
0838   ///////////////////////////////////////////////
0839   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
0840   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
0841   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
0842   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
0843 
0844   ///////////////////////////////////////////////
0845   // 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
0846   ///////////////////////////////////////////////
0847   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
0848   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
0849   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
0850   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
0851 
0852   ///////////////////////////////////////////////
0853   // 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
0854   ///////////////////////////////////////////////
0855   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
0856   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
0857   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
0858   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
0859 
0860   ///////////////////////////////////////////////
0861   // 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
0862   ///////////////////////////////////////////////
0863   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
0864   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
0865   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
0866   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
0867 
0868   ///////////////////////////////////////////////
0869   // 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
0870   ///////////////////////////////////////////////
0871   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
0872   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
0873   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
0874   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
0875 
0876 
0877 
0878 }