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 = "../../sPHENIX_work//test_production/Upsilon/spacal2d/fieldon/SimALL_PythiaUpsilon.root_Ana.root" 0028 0029 // const TString infile = "../../sPHENIX_work/production_analysis//emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_8GeV-ALL.root_Ana.root.lst_EMCalLikelihood.root" 0030 const TString infile = "../../sPHENIX_work/production_analysis_cemc2x2//emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_8GeV-ALL.root_Ana.root.lst_EMCalLikelihood.root" 0031 0032 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_8GeV-ALL.root_Ana.root" 0033 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_8GeV-ALL.root_Ana.root" 0034 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_4GeV-ALL.root_Ana.root" 0035 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_4GeV-ALL.root_Ana.root" 0036 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_2GeV-ALL.root_Ana.root" 0037 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_2GeV-ALL.root_Ana.root" 0038 0039 // const TString infile = "../../sPHENIX_work/production_analysis/embedding/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_8GeV-ALL.root_Ana.root" 0040 // const TString infile = "../../sPHENIX_work/production_analysis/embedding/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_8GeV-ALL.root_Ana.root" 0041 // const TString infile = "../../sPHENIX_work/production_analysis/embedding/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_4GeV-ALL.root_Ana.root" 0042 // const TString infile = "../../sPHENIX_work/production_analysis/embedding/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_4GeV-ALL.root_Ana.root" 0043 // const TString infile = "../../sPHENIX_work/production_analysis/embedding/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_2GeV-ALL.root_Ana.root" 0044 // const TString infile = "../../sPHENIX_work/production_analysis/embedding/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_2GeV-ALL.root_Ana.root" 0045 0046 // const TString infile = "../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root" 0047 // const TString infile = "../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_8GeV-ALL.root_Ana.root" 0048 // const TString infile = "../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root" 0049 // const TString infile = "../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root" 0050 // const TString infile = "../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root" 0051 // const TString infile = "../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root" 0052 // 0053 // const TString infile = "../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root" 0054 // const TString infile = "../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_8GeV-ALL.root_Ana.root" 0055 // const TString infile = "../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root" 0056 // const TString infile = "../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root" 0057 // const TString infile = "../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root" 0058 // const TString infile = "../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root" 0059 // 0060 // const TString infile = "../../sPHENIX_work/production_analysis/embedding/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_8GeV-ALL.root_Ana.root" 0061 // const TString infile = "../../sPHENIX_work/production_analysis/embedding/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root" 0062 // const TString infile = "../../sPHENIX_work/production_analysis/embedding/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root" 0063 // const TString infile = "../../sPHENIX_work/production_analysis/embedding/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root" 0064 // const TString infile = "../../sPHENIX_work/production_analysis/embedding/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root" 0065 // const TString infile = "../../sPHENIX_work/production_analysis/embedding/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root" 0066 0067 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_8GeV-ALL.root_Ana.root" 0068 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_8GeV-ALL.root_Ana.root" 0069 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_4GeV-ALL.root_Ana.root" 0070 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_4GeV-ALL.root_Ana.root" 0071 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_2GeV-ALL.root_Ana.root" 0072 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_2GeV-ALL.root_Ana.root" 0073 0074 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_8GeV-ALL.root_Ana.root" 0075 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_8GeV-ALL.root_Ana.root" 0076 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_4GeV-ALL.root_Ana.root" 0077 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_4GeV-ALL.root_Ana.root" 0078 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_2GeV-ALL.root_Ana.root" 0079 // const TString infile = "../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_2GeV-ALL.root_Ana.root" 0080 0081 // const TString infile = "G4Hits_sPHENIX_e-_eta0_8GeV-0000.root_Ana.root" // 0082 // const TString infile = "G4Hits_sPHENIX_pi-_eta0_8GeV-0000.root_Ana.root"// 0083 ) 0084 { 0085 SetOKStyle(); 0086 gStyle->SetOptStat(0); 0087 gStyle->SetOptFit(1111); 0088 TVirtualFitter::SetDefaultFitter("Minuit2"); 0089 0090 gSystem->Load("libg4eval.so"); 0091 gSystem->Load("libemcal_ana.so"); 0092 gSystem->Load("libg4vertex.so"); 0093 0094 if (!_file0) 0095 { 0096 TString chian_str = infile; 0097 chian_str.ReplaceAll("ALL", "*"); 0098 chian_str.ReplaceAll("+", "\\+"); 0099 0100 TChain * t = new TChain("T"); 0101 const int n = t->Add(chian_str); 0102 0103 cout << "Loaded " << n << " root files with " << chian_str << endl; 0104 assert(n>0); 0105 0106 T = t; 0107 0108 _file0 = new TFile; 0109 _file0->SetName(infile); 0110 0111 fstream flst(infile + ".lst", ios_base::out); 0112 0113 TObjArray *fileElements = t->GetListOfFiles(); 0114 TIter next(fileElements); 0115 TChainElement *chEl = 0; 0116 while ((chEl = (TChainElement*) next())) 0117 { 0118 flst << chEl->GetTitle() << endl; 0119 } 0120 flst.close(); 0121 0122 cout << "Saved file list to " << infile + ".lst" << endl; 0123 } 0124 0125 assert(_file0); 0126 0127 T->SetAlias("UpsilonPair_trk_gpt", 0128 "1*sqrt(DST.UpsilonPair.trk.gpx**2 + DST.UpsilonPair.trk.gpy**2)"); 0129 T->SetAlias("UpsilonPair_trk_pt", 0130 "1*sqrt(DST.UpsilonPair.trk.px**2 + DST.UpsilonPair.trk.py**2)"); 0131 0132 T->SetAlias("EMCalTrk_pt", "1*sqrt(DST.EMCalTrk.px**2 + DST.EMCalTrk.py**2)"); 0133 T->SetAlias("EMCalTrk_gpt", 0134 "1*sqrt(DST.EMCalTrk.gpx**2 + DST.EMCalTrk.gpy**2)"); 0135 0136 const TCut event_sel = "1*1"; 0137 cuts = "_all_event"; 0138 0139 cout << "Build event selection of " << (const char *) event_sel << endl; 0140 0141 T->Draw(">>EventList", event_sel); 0142 TEventList * elist = gDirectory->GetObjectChecked("EventList", "TEventList"); 0143 cout << elist->GetN() << " / " << T->GetEntriesFast() << " events selected" 0144 << endl; 0145 0146 T->SetEventList(elist); 0147 //// 0148 // UpsilonPair_Checks(infile); 0149 0150 // TrackProjection_Checks(infile, "abs(EMCalTrk_pt/EMCalTrk_gpt - 1)<0.05"); 0151 // Edep_Checks(infile, 1.4, 1.4, "abs(EMCalTrk_pt/EMCalTrk_gpt - 1)<0.05"); 0152 0153 TrackProjection_Checks(infile, "abs(EMCalTrk_pt/EMCalTrk_gpt - 1)<0.05 && DST.EMCalTrk.cemce>2"); 0154 Edep_Checks(infile, 1.4, 1.4, "abs(EMCalTrk_pt/EMCalTrk_gpt - 1)<0.05 && DST.EMCalTrk.cemce>2"); 0155 } 0156 0157 void 0158 MakeRadiusCut(TString infile, const double R_CEMC, const double R_HCALIN, 0159 TCut good_track_cut) 0160 { 0161 TCanvas *c1 = new TCanvas("MakeRadiusCut" + cuts, "MakeRadiusCut" + cuts, 0162 1800, 900); 0163 c1->Divide(4, 2); 0164 int idx = 1; 0165 TPad * p; 0166 0167 p = (TPad *) c1->cd(idx++); 0168 c1->Update(); 0169 p->SetLogy(); 0170 0171 T->Draw("DST.EMCalTrk.cemc_iphi>>hcemc_iphi(130,-6.5,6.5)", 0172 TString("DST.EMCalTrk.cemc_energy * (") 0173 + TString(good_track_cut.GetTitle()) + ")"); 0174 hcemc_iphi->SetTitle( 0175 ";CEMC Azimuthal Distance (Tower Width);Energy Distribution"); 0176 hcemc_iphi->Fit("gaus", "M"); 0177 TF1* f = (TF1*) (hcemc_iphi->GetListOfFunctions()->At(0)); 0178 assert(f); 0179 T->SetAlias("EMCalTrk_cor_cemc_iphi", 0180 Form("DST.EMCalTrk.cemc_iphi - %f", f->GetParameter(1))); 0181 0182 const double center_cemc_iphi = f->GetParameter(1); 0183 0184 p = (TPad *) c1->cd(idx++); 0185 c1->Update(); 0186 p->SetLogy(); 0187 0188 T->Draw("abs(EMCalTrk_cor_cemc_iphi)>>hEMCalTrk_cor_cemc_iphi(130,0,6.5)", 0189 TString("DST.EMCalTrk.cemc_energy * (") 0190 + TString(good_track_cut.GetTitle()) + ")"); 0191 0192 p->DrawFrame(0, 1e-3, 6.5, 1, 0193 ";CEMC Azimuthal Distance (Tower Width);Energy Leakage Ratio"); 0194 TGraphErrors * ge = Distribution2Efficiency(hEMCalTrk_cor_cemc_iphi); 0195 ge->SetLineColor(kBlue + 2); 0196 ge->SetMarkerColor(kBlue + 21); 0197 ge->SetMarkerColor(kFullCircle); 0198 ge->SetLineWidth(3); 0199 ge->Draw("lp"); 0200 0201 p = (TPad *) c1->cd(idx++); 0202 c1->Update(); 0203 p->SetLogy(); 0204 0205 T->Draw("DST.EMCalTrk.cemc_ieta>>hcemc_ieta(130,-6.5,6.5)", 0206 TString("DST.EMCalTrk.cemc_energy * (") 0207 + TString(good_track_cut.GetTitle()) + ")"); 0208 hcemc_ieta->SetTitle( 0209 ";CEMC Polar Distance (Tower Width);Energy Distribution"); 0210 hcemc_ieta->Fit("gaus", "M"); 0211 TF1* f = (TF1*) (hcemc_ieta->GetListOfFunctions()->At(0)); 0212 assert(f); 0213 T->SetAlias("EMCalTrk_cor_cemc_ieta", 0214 Form("DST.EMCalTrk.cemc_ieta - %f", f->GetParameter(1))); 0215 0216 const double center_cemc_ieta = f->GetParameter(1); 0217 0218 p = (TPad *) c1->cd(idx++); 0219 c1->Update(); 0220 p->SetLogy(); 0221 0222 T->Draw("abs(EMCalTrk_cor_cemc_ieta)>>hEMCalTrk_cor_cemc_ieta(130,0,6.5)", 0223 TString("DST.EMCalTrk.cemc_energy * (") 0224 + TString(good_track_cut.GetTitle()) + ")"); 0225 0226 p->DrawFrame(0, 1e-3, 6.5, 1, 0227 ";CEMC Polar Distance (Tower Width);Energy Leakage Ratio"); 0228 TGraphErrors * ge = Distribution2Efficiency(hEMCalTrk_cor_cemc_ieta); 0229 ge->SetLineColor(kBlue + 2); 0230 ge->SetMarkerColor(kBlue + 21); 0231 ge->SetMarkerColor(kFullCircle); 0232 ge->SetLineWidth(3); 0233 ge->Draw("lp"); 0234 0235 p = (TPad *) c1->cd(idx++); 0236 c1->Update(); 0237 p->SetLogy(); 0238 0239 T->Draw("DST.EMCalTrk.hcalin_iphi>>hhcalin_iphi(130,-6.5,6.5)", 0240 TString("DST.EMCalTrk.hcalin_energy * (") 0241 + TString(good_track_cut.GetTitle()) + ")"); 0242 hhcalin_iphi->SetTitle( 0243 ";HCal_{IN} Azimuthal Distance (Tower Width);Energy Distribution"); 0244 hhcalin_iphi->Fit("gaus", "M"); 0245 TF1* f = (TF1*) (hhcalin_iphi->GetListOfFunctions()->At(0)); 0246 assert(f); 0247 T->SetAlias("EMCalTrk_cor_hcalin_iphi", 0248 Form("DST.EMCalTrk.hcalin_iphi - %f", f->GetParameter(1))); 0249 0250 const double center_hcalin_iphi = f->GetParameter(1); 0251 0252 p = (TPad *) c1->cd(idx++); 0253 c1->Update(); 0254 p->SetLogy(); 0255 0256 T->Draw("abs(EMCalTrk_cor_hcalin_iphi)>>hEMCalTrk_cor_hcalin_iphi(130,0,6.5)", 0257 TString("DST.EMCalTrk.hcalin_energy * (") 0258 + TString(good_track_cut.GetTitle()) + ")"); 0259 0260 p->DrawFrame(0, 1e-3, 6.5, 1, 0261 ";HCal_{IN} Azimuthal Distance (Tower Width);Energy Leakage Ratio"); 0262 TGraphErrors * ge = Distribution2Efficiency(hEMCalTrk_cor_hcalin_iphi); 0263 ge->SetLineColor(kBlue + 2); 0264 ge->SetMarkerColor(kBlue + 21); 0265 ge->SetMarkerColor(kFullCircle); 0266 ge->SetLineWidth(3); 0267 ge->Draw("lp"); 0268 0269 p = (TPad *) c1->cd(idx++); 0270 c1->Update(); 0271 p->SetLogy(); 0272 0273 T->Draw("DST.EMCalTrk.hcalin_ieta>>hhcalin_ieta(130,-6.5,6.5)", 0274 TString("DST.EMCalTrk.hcalin_energy * (") 0275 + TString(good_track_cut.GetTitle()) + ")"); 0276 hhcalin_ieta->SetTitle( 0277 ";HCal_{IN} Polar Distance (Tower Width);Energy Distribution"); 0278 hhcalin_ieta->Fit("gaus", "M"); 0279 TF1* f = (TF1*) (hhcalin_ieta->GetListOfFunctions()->At(0)); 0280 assert(f); 0281 T->SetAlias("EMCalTrk_cor_hcalin_ieta", 0282 Form("DST.EMCalTrk.hcalin_ieta - %f", f->GetParameter(1))); 0283 0284 const double center_hcalin_ieta = f->GetParameter(1); 0285 0286 p = (TPad *) c1->cd(idx++); 0287 c1->Update(); 0288 p->SetLogy(); 0289 0290 T->Draw("abs(EMCalTrk_cor_hcalin_ieta)>>hEMCalTrk_cor_hcalin_ieta(130,0,6.5)", 0291 TString("DST.EMCalTrk.hcalin_energy * (") 0292 + TString(good_track_cut.GetTitle()) + ")"); 0293 0294 p->DrawFrame(0, 1e-3, 6.5, 1, 0295 ";HCal_{IN} Polar Distance (Tower Width);Energy Leakage Ratio"); 0296 TGraphErrors * ge = Distribution2Efficiency(hEMCalTrk_cor_hcalin_ieta); 0297 ge->SetLineColor(kBlue + 2); 0298 ge->SetMarkerColor(kBlue + 21); 0299 ge->SetMarkerColor(kFullCircle); 0300 ge->SetLineWidth(3); 0301 ge->Draw("lp"); 0302 0303 T->SetAlias("EMCalTrk_cemc_e", 0304 Form( 0305 "1*Sum$( DST.EMCalTrk.cemc_energy * ( sqrt(EMCalTrk_cor_cemc_iphi*EMCalTrk_cor_cemc_iphi + EMCalTrk_cor_cemc_ieta*EMCalTrk_cor_cemc_ieta) < %f ) )", 0306 R_CEMC)); 0307 T->SetAlias("EMCalTrk_hcalin_e", 0308 Form( 0309 "1*Sum$( DST.EMCalTrk.hcalin_energy * ( sqrt(EMCalTrk_cor_hcalin_iphi*EMCalTrk_cor_hcalin_iphi + EMCalTrk_cor_hcalin_ieta*EMCalTrk_cor_hcalin_ieta) < %f ) )", 0310 R_CEMC)); 0311 0312 T->SetAlias("EMCalTrk_cemc_ntower", 0313 Form( 0314 "1*Sum$(( sqrt(EMCalTrk_cor_cemc_iphi*EMCalTrk_cor_cemc_iphi + EMCalTrk_cor_cemc_ieta*EMCalTrk_cor_cemc_ieta) < %f ) )", 0315 R_HCALIN)); 0316 T->SetAlias("EMCalTrk_hcalin_ntower", 0317 Form( 0318 "1*Sum$( ( sqrt(EMCalTrk_cor_hcalin_iphi*EMCalTrk_cor_hcalin_iphi + EMCalTrk_cor_hcalin_ieta*EMCalTrk_cor_hcalin_ieta) < %f ) )", 0319 R_HCALIN)); 0320 0321 SaveCanvas(c1, 0322 TString(_file0->GetName()) + TString("_DrawEcal_pDST_") 0323 + Form("RCEMC%.1f_RCEMC%.1f_", R_CEMC, R_HCALIN) 0324 + TString(c1->GetName()), kFALSE); 0325 0326 cout << "///////////////////////////////////////////////" << endl; 0327 cout << "// Projection center based on " << _file0->GetName() << endl; 0328 cout << "///////////////////////////////////////////////" << endl; 0329 0330 cout << "const double center_cemc_iphi = " << center_cemc_iphi << "; // " 0331 << _file0->GetName() << endl; 0332 cout << "const double center_cemc_ieta = " << center_cemc_ieta << "; // " 0333 << _file0->GetName() << endl; 0334 cout << "const double center_hcalin_iphi = " << center_hcalin_iphi << "; // " 0335 << _file0->GetName() << endl; 0336 cout << "const double center_hcalin_ieta = " << center_hcalin_ieta << "; // " 0337 << _file0->GetName() << endl; 0338 0339 } 0340 0341 void 0342 Edep_Checks(TString infile, const double R_CEMC, const double R_HCALIN, 0343 TCut good_track_cut) 0344 { 0345 MakeRadiusCut(infile, R_CEMC, R_HCALIN, good_track_cut); 0346 0347 double N_Event = T->GetEntries(); 0348 0349 TCanvas *c1 = new TCanvas("Edep_Checks" + cuts, "Edep_Checks" + cuts, 1900, 0350 950); 0351 c1->Divide(2, 2); 0352 int idx = 1; 0353 TPad * p; 0354 0355 p = (TPad *) c1->cd(idx++); 0356 c1->Update(); 0357 p->SetLogy(); 0358 T->Draw("EMCalTrk_cemc_ntower>>hEMCalTrk_cemc_ntower(16,-.5,15.5)", 0359 good_track_cut); 0360 hEMCalTrk_cemc_ntower->SetTitle( 0361 Form("CEMC Cluster Size (R = %.1f);Cluster Size (Towers);Probability", 0362 R_CEMC)); 0363 hEMCalTrk_cemc_ntower->Scale(1. / N_Event); 0364 0365 p = (TPad *) c1->cd(idx++); 0366 c1->Update(); 0367 p->SetLogy(); 0368 T->Draw("EMCalTrk_hcalin_ntower>>hEMCalTrk_hcalin_ntower(16,-.5,15.5)", 0369 good_track_cut); 0370 hEMCalTrk_hcalin_ntower->SetTitle( 0371 Form( 0372 "HCal_{in} Cluster Size (R = %.1f);Cluster Size (Towers);Probability", 0373 R_HCALIN)); 0374 hEMCalTrk_hcalin_ntower->Scale(1. / N_Event); 0375 0376 p = (TPad *) c1->cd(idx++); 0377 c1->Update(); 0378 p->SetLogy(); 0379 T->Draw("EMCalTrk_cemc_e>>hEMCalTrk_cemc_e(240,-.0,12)", good_track_cut); 0380 hEMCalTrk_cemc_e->SetTitle( 0381 Form( 0382 "CEMC Cluster Energy (R = %.1f);Cluster Energy (/bin);Probability/bin", 0383 R_CEMC)); 0384 hEMCalTrk_cemc_e->Scale(1. / N_Event); 0385 0386 p = (TPad *) c1->cd(idx++); 0387 c1->Update(); 0388 p->SetLogy(); 0389 T->Draw("EMCalTrk_hcalin_e>>hEMCalTrk_hcalin_e(240,-.0,12)", good_track_cut); 0390 hEMCalTrk_hcalin_e->SetTitle( 0391 Form( 0392 "HCal_{in} Cluster Energy (R = %.1f);Cluster Energy (GeV);Probability/bin", 0393 R_HCALIN)); 0394 hEMCalTrk_hcalin_e->Scale(1. / N_Event); 0395 0396 SaveCanvas(c1, 0397 TString(_file0->GetName()) + TString("_DrawEcal_pDST_") 0398 + Form("RCEMC%.1f_RCEMC%.1f_", R_CEMC, R_HCALIN) 0399 + TString(c1->GetName()), kFALSE); 0400 0401 TCanvas *c1 = new TCanvas("Edep_Checks_2D" + cuts, "Edep_Checks_2D" + cuts, 0402 900, 900); 0403 // c1->Divide(2, 2); 0404 // int idx = 1; 0405 // TPad * p; 0406 0407 p = (TPad *) c1->cd(idx++); 0408 c1->Update(); 0409 p->SetLogz(); 0410 T->Draw( 0411 "EMCalTrk_hcalin_e:EMCalTrk_cemc_e>>h2_EMCalTrk_hcalin_e_EMCalTrk_cemc_e(240,-.0,12, 240,-.0,12)", 0412 good_track_cut, "colz"); 0413 h2_EMCalTrk_hcalin_e_EMCalTrk_cemc_e->SetTitle( 0414 Form( 0415 "Energy distribution;CEMC Cluster Energy (R = %.1f) in GeV;HCal_{in} Cluster Energy (R = %.1f) in GeV", 0416 R_CEMC, R_HCALIN)); 0417 h2_EMCalTrk_hcalin_e_EMCalTrk_cemc_e->Scale(1. / N_Event); 0418 h2_EMCalTrk_hcalin_e_EMCalTrk_cemc_e->GetZaxis()->SetRangeUser(1e-6, 1); 0419 0420 SaveCanvas(c1, 0421 TString(_file0->GetName()) + TString("_DrawEcal_pDST_") 0422 + Form("RCEMC%.1f_RCEMC%.1f_", R_CEMC, R_HCALIN) 0423 + TString(c1->GetName()), kFALSE); 0424 0425 } 0426 0427 void 0428 TrackProjection_Checks(TString infile, TCut good_track_cut) 0429 { 0430 0431 double N_Event = T->GetEntries(); 0432 0433 TCanvas *c1 = new TCanvas("TrackProjection_Checks_Raw" + cuts, 0434 "TrackProjection_Checks_Raw" + cuts, 1800, 900); 0435 c1->Divide(2, 1); 0436 int idx = 1; 0437 TPad * p; 0438 0439 p = (TPad *) c1->cd(idx++); 0440 c1->Update(); 0441 p->SetLogz(); 0442 T->Draw( 0443 "DST.EMCalTrk.cemc_iphi:DST.EMCalTrk.cemc_ieta>>hcemc(130,-6.5,6.5,130,-6.5,6.5)", 0444 "DST.EMCalTrk.cemc_energy", "colz"); 0445 hcemc->SetTitle( 0446 "CEMC Tower Energy Distribution;Polar Distance (Tower Width);Azimuthal Distance (Tower Width)"); 0447 hcemc->Scale(1. / N_Event); 0448 hcemc->GetZaxis()->SetRangeUser(1e-5, 30); 0449 0450 p = (TPad *) c1->cd(idx++); 0451 c1->Update(); 0452 p->SetLogz(); 0453 T->Draw( 0454 "DST.EMCalTrk.hcalin_iphi:DST.EMCalTrk.hcalin_ieta>>hcalin(130,-6.5,6.5,130,-6.5,6.5)", 0455 "DST.EMCalTrk.hcalin_energy", "colz"); 0456 hcalin->SetTitle( 0457 "HCal_{in} Tower Energy Distribution;Polar Distance (Tower Width);Azimuthal Distance (Tower Width)"); 0458 hcalin->Scale(1. / N_Event); 0459 hcalin->GetZaxis()->SetRangeUser(1e-5, 30); 0460 0461 SaveCanvas(c1, 0462 TString(_file0->GetName()) + TString("_DrawEcal_pDST_") 0463 + TString(c1->GetName()), kFALSE); 0464 0465 TCanvas *c1 = new TCanvas("TrackProjection_Checks_Tracking" + cuts, 0466 "TrackProjection_Checks_Tracking" + cuts, 1800, 900); 0467 c1->Divide(2, 1); 0468 int idx = 1; 0469 TPad * p; 0470 0471 p = (TPad *) c1->cd(idx++); 0472 c1->Update(); 0473 p->SetLogy(); 0474 T->Draw("EMCalTrk_pt/EMCalTrk_gpt>>htracking_pt(100,0,2)", "", ""); 0475 htracking_pt->SetTitle( 0476 "Tracking reco precision;p_{T, reco}/p_{T, truth};Probability per bin"); 0477 htracking_pt->Scale(1. / N_Event); 0478 T->Draw("EMCalTrk_pt/EMCalTrk_gpt>>htracking_pt_cut(100,0,2)", good_track_cut, 0479 "same"); 0480 htracking_pt_cut->Scale(1. / N_Event); 0481 htracking_pt_cut->SetLineColor(kRed); 0482 0483 p = (TPad *) c1->cd(idx++); 0484 c1->Update(); 0485 p->SetLogy(); 0486 0487 T->Draw("EMCalTrk_pt/EMCalTrk_gpt>0>>h_have_trak(3,-.5,2.5)", "", ""); 0488 h_have_trak->Scale(1. / N_Event); 0489 T->Draw("2>>h_have_trak_cut(3,-.5,2.5)", good_track_cut, "same"); 0490 h_have_trak_cut->Scale(1. / N_Event); 0491 h_have_trak_cut->SetLineColor(kRed); 0492 0493 SaveCanvas(c1, 0494 TString(_file0->GetName()) + TString("_DrawEcal_pDST_") 0495 + TString(c1->GetName()), kFALSE); 0496 0497 TCanvas *c1 = new TCanvas("TrackProjection_Checks_Cut" + cuts, 0498 "TrackProjection_Checks_Cut" + cuts, 1800, 900); 0499 c1->Divide(2, 1); 0500 int idx = 1; 0501 TPad * p; 0502 0503 p = (TPad *) c1->cd(idx++); 0504 c1->Update(); 0505 p->SetLogz(); 0506 T->Draw( 0507 "DST.EMCalTrk.cemc_iphi:DST.EMCalTrk.cemc_ieta>>hcemcc(130,-6.5,6.5,130,-6.5,6.5)", 0508 TString("DST.EMCalTrk.cemc_energy * (") 0509 + TString(good_track_cut.GetTitle()) + ")", "colz"); 0510 hcemcc->SetTitle( 0511 "CEMC Tower Energy Distribution;Polar Distance (Tower Width);Azimuthal Distance (Tower Width)"); 0512 hcemcc->Scale(1. / N_Event); 0513 hcemcc->GetZaxis()->SetRangeUser(1e-5, 30); 0514 0515 p = (TPad *) c1->cd(idx++); 0516 c1->Update(); 0517 p->SetLogz(); 0518 T->Draw( 0519 "DST.EMCalTrk.hcalin_iphi:DST.EMCalTrk.hcalin_ieta>>hcalinc(130,-6.5,6.5,130,-6.5,6.5)", 0520 TString("DST.EMCalTrk.hcalin_energy * (") 0521 + TString(good_track_cut.GetTitle()) + ")", "colz"); 0522 hcalinc->SetTitle( 0523 "HCal_{in} Tower Energy Distribution;Polar Distance (Tower Width);Azimuthal Distance (Tower Width)"); 0524 hcalinc->Scale(1. / N_Event); 0525 hcalinc->GetZaxis()->SetRangeUser(1e-5, 30); 0526 0527 SaveCanvas(c1, 0528 TString(_file0->GetName()) + TString("_DrawEcal_pDST_") 0529 + TString(c1->GetName()), kFALSE); 0530 0531 TCanvas *c1 = new TCanvas("TrackProjection_Checks_Proj" + cuts, 0532 "TrackProjection_Checks_Proj" + cuts, 1800, 900); 0533 c1->Divide(2, 2); 0534 int idx = 1; 0535 TPad * p; 0536 0537 p = (TPad *) c1->cd(idx++); 0538 c1->Update(); 0539 p->SetLogy(); 0540 0541 TH1 * proj = hcemc->ProjectionX(); 0542 proj->SetLineColor(kBlack); 0543 proj->Draw(""); 0544 TH1 * proj = hcemcc->ProjectionX(); 0545 proj->SetLineColor(kRed); 0546 proj->Draw("same"); 0547 0548 p = (TPad *) c1->cd(idx++); 0549 c1->Update(); 0550 p->SetLogy(); 0551 0552 TH1 * proj = hcalin->ProjectionX(); 0553 proj->SetLineColor(kBlack); 0554 proj->Draw(""); 0555 TH1 * proj = hcalinc->ProjectionX(); 0556 proj->SetLineColor(kRed); 0557 proj->Draw("same"); 0558 0559 p = (TPad *) c1->cd(idx++); 0560 c1->Update(); 0561 p->SetLogy(); 0562 0563 TH1 * proj = hcemc->ProjectionY(); 0564 proj->SetLineColor(kBlack); 0565 proj->Draw(""); 0566 TH1 * proj = hcemcc->ProjectionY(); 0567 proj->SetLineColor(kRed); 0568 proj->Draw("same"); 0569 0570 p = (TPad *) c1->cd(idx++); 0571 c1->Update(); 0572 p->SetLogy(); 0573 0574 TH1 * proj = hcalin->ProjectionY(); 0575 proj->SetLineColor(kBlack); 0576 proj->Draw(""); 0577 TH1 * proj = hcalinc->ProjectionY(); 0578 proj->SetLineColor(kRed); 0579 proj->Draw("same"); 0580 0581 SaveCanvas(c1, 0582 TString(_file0->GetName()) + TString("_DrawEcal_pDST_") 0583 + TString(c1->GetName()), kFALSE); 0584 } 0585 0586 void 0587 UpsilonPair_Checks(TString infile) 0588 { 0589 0590 TCanvas *c1 = new TCanvas("UpsilonPair_Checks" + cuts, 0591 "UpsilonPair_Checks" + cuts, 1800, 900); 0592 c1->Divide(3, 2); 0593 int idx = 1; 0594 TPad * p; 0595 0596 p = (TPad *) c1->cd(idx++); 0597 c1->Update(); 0598 T->Draw("DST.UpsilonPair.gmass>>gmass(100,0,10)"); 0599 T->Draw("DST.UpsilonPair.gmass>>gmassc(100,0,10)", 0600 "DST.UpsilonPair.good_upsilon", "same"); 0601 gmassc->SetLineColor(kRed); 0602 0603 p = (TPad *) c1->cd(idx++); 0604 c1->Update(); 0605 T->Draw("DST.UpsilonPair.mass>>mass(100,0,10)"); 0606 T->Draw("DST.UpsilonPair.mass>>massc(100,0,10)", 0607 "DST.UpsilonPair.good_upsilon", "same"); 0608 massc->SetLineColor(kRed); 0609 0610 p = (TPad *) c1->cd(idx++); 0611 c1->Update(); 0612 T->Draw("UpsilonPair_trk_gpt>>UpsilonPair_trk_gpt(100,0,10)"); 0613 0614 p = (TPad *) c1->cd(idx++); 0615 c1->Update(); 0616 T->Draw("UpsilonPair_trk_pt:UpsilonPair_trk_gpt>>pt2(100,0,10,100,0,10)", "", 0617 "colz"); 0618 0619 p = (TPad *) c1->cd(idx++); 0620 c1->Update(); 0621 T->Draw("UpsilonPair_trk_pt:UpsilonPair_trk_gpt>>pt2c(100,0,10,100,0,10)", 0622 "DST.UpsilonPair.good_upsilon", "colz"); 0623 0624 SaveCanvas(c1, 0625 TString(_file0->GetName()) + TString("_DrawEcal_pDST_") 0626 + TString(c1->GetName()), kFALSE); 0627 } 0628 0629 TGraphErrors * 0630 Distribution2Efficiency(TH1F * hCEMC3_Max) 0631 { 0632 double threshold[10000] = 0633 { 0 }; 0634 double eff[10000] = 0635 { 0 }; 0636 double eff_err[10000] = 0637 { 0 }; 0638 0639 assert(hCEMC3_Max->GetNbinsX()<10000); 0640 0641 const double n = hCEMC3_Max->GetSum(); 0642 double pass = 0; 0643 int cnt = 0; 0644 for (int i = hCEMC3_Max->GetNbinsX(); i >= 1; i--) 0645 { 0646 pass += hCEMC3_Max->GetBinContent(i); 0647 0648 const double pp = pass / n; 0649 // const double z = 1.96; 0650 const double z = 1.; 0651 0652 const double A = z * sqrt(1. / n * pp * (1 - pp) + z * z / 4 / n / n); 0653 const double B = 1 / (1 + z * z / n); 0654 0655 threshold[cnt] = hCEMC3_Max->GetBinCenter(i); 0656 eff[cnt] = (pp + z * z / 2 / n) * B; 0657 eff_err[cnt] = A * B; 0658 0659 cout << threshold[cnt] << ": " << "CL " << eff[cnt] << "+/-" 0660 << eff_err[cnt] << endl; 0661 cnt++; 0662 } 0663 TGraphErrors * ge = new TGraphErrors(cnt, threshold, eff, NULL, eff_err); 0664 0665 return ge; 0666 } 0667 0668 void 0669 TrackProjection_Checks_Comparison() 0670 { 0671 0672 SetOKStyle(); 0673 gStyle->SetOptStat(0); 0674 gStyle->SetOptFit(1111); 0675 TVirtualFitter::SetDefaultFitter("Minuit2"); 0676 0677 gSystem->Load("libg4eval.so"); 0678 gSystem->Load("libemcal_ana.so"); 0679 gSystem->Load("libg4vertex.so"); 0680 0681 0682 TFile * f = 0683 new TFile( 0684 "../../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"); 0685 0686 assert(f->IsOpen()); 0687 TH2F * hcemcc_2d = (TH2F *) f->GetObjectChecked("hcemcc", "TH2F"); 0688 assert(hcemcc_2d); 0689 0690 TFile * f = 0691 new TFile( 0692 "../../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"); 0693 0694 assert(f->IsOpen()); 0695 TH2F * hcemcc_1d = (TH2F *) f->GetObjectChecked("hcemcc", "TH2F"); 0696 assert(hcemcc_1d); 0697 0698 TCanvas *c1 = new TCanvas("TrackProjection_Checks_Comparison" , 0699 "TrackProjection_Checks_Comparison" , 1800, 1000); 0700 c1->Divide(2, 1); 0701 int idx = 1; 0702 TPad * p; 0703 0704 p = (TPad *) c1->cd(idx++); 0705 c1->Update(); 0706 p->SetLogz(); 0707 0708 hcemcc_2d -> Draw("colz"); 0709 hcemcc_2d->SetTitle(";Polar Distance (Tower Width);Azimuthal Distance (Tower Width)"); 0710 0711 TLatex *text = new TLatex(0,7,"2-D Projective SPACAL, 0.9 < #eta_{e^{-}} < 1.0, E_{e^{-}} = 8 GeV"); 0712 text->SetTextAlign(22); 0713 text->SetTextSize(0.044); 0714 text->Draw(); 0715 0716 0717 p = (TPad *) c1->cd(idx++); 0718 c1->Update(); 0719 p->SetLogz(); 0720 0721 hcemcc_1d -> Draw("colz"); 0722 hcemcc_1d->SetTitle(";Polar Distance (Tower Width);Azimuthal Distance (Tower Width)"); 0723 0724 TLatex *text = new TLatex(0,7,"1-D Projective SPACAL, 0.9 < #eta_{e^{-}} < 1.0, E_{e^{-}} = 8 GeV"); 0725 text->SetTextAlign(22); 0726 text->SetTextSize(0.044); 0727 text->Draw(); 0728 0729 0730 SaveCanvas(c1, 0731 "../../sPHENIX_work/production_analysis/" + TString("DrawEcal_pDST_") 0732 + TString(c1->GetName())); 0733 } 0734 0735 void 0736 Logs() 0737 { 0738 0739 /////////////////////////////////////////////// 0740 // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0_8GeV-ALL.root_Ana.root 0741 /////////////////////////////////////////////// 0742 const double center_cemc_iphi = 0.0202159; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0_8GeV-ALL.root_Ana.root 0743 const double center_cemc_ieta = 0.00453938; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0_8GeV-ALL.root_Ana.root 0744 const double center_hcalin_iphi = 0.488432; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0_8GeV-ALL.root_Ana.root 0745 const double center_hcalin_ieta = -0.00936002; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0_8GeV-ALL.root_Ana.root 0746 0747 /////////////////////////////////////////////// 0748 // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0_8GeV-ALL.root_Ana.root 0749 /////////////////////////////////////////////// 0750 const double center_cemc_iphi = 0.0402477; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0_8GeV-ALL.root_Ana.root 0751 const double center_cemc_ieta = 0.00574989; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0_8GeV-ALL.root_Ana.root 0752 const double center_hcalin_iphi = 0.331081; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0_8GeV-ALL.root_Ana.root 0753 const double center_hcalin_ieta = -0.000527382; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0_8GeV-ALL.root_Ana.root 0754 0755 0756 /////////////////////////////////////////////// 0757 // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_8GeV-ALL.root_Ana.root 0758 /////////////////////////////////////////////// 0759 const double center_cemc_iphi = 0.0649585; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_8GeV-ALL.root_Ana.root 0760 const double center_cemc_ieta = 0.00465546; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_8GeV-ALL.root_Ana.root 0761 const double center_hcalin_iphi = 0.552238; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_8GeV-ALL.root_Ana.root 0762 const double center_hcalin_ieta = -0.0102162; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_8GeV-ALL.root_Ana.root 0763 0764 /////////////////////////////////////////////// 0765 // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_8GeV-ALL.root_Ana.root 0766 /////////////////////////////////////////////// 0767 const double center_cemc_iphi = 0.044375; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_8GeV-ALL.root_Ana.root 0768 const double center_cemc_ieta = 0.00680547; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_8GeV-ALL.root_Ana.root 0769 const double center_hcalin_iphi = 0.336741; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_8GeV-ALL.root_Ana.root 0770 const double center_hcalin_ieta = 0.00108616; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_8GeV-ALL.root_Ana.root 0771 0772 /////////////////////////////////////////////// 0773 // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_4GeV-ALL.root_Ana.root 0774 /////////////////////////////////////////////// 0775 const double center_cemc_iphi = 0.0952638; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_4GeV-ALL.root_Ana.root 0776 const double center_cemc_ieta = 0.00397331; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_4GeV-ALL.root_Ana.root 0777 const double center_hcalin_iphi = 0.592338; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_4GeV-ALL.root_Ana.root 0778 const double center_hcalin_ieta = -0.0118107; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_4GeV-ALL.root_Ana.root 0779 0780 /////////////////////////////////////////////// 0781 // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_4GeV-ALL.root_Ana.root 0782 /////////////////////////////////////////////// 0783 const double center_cemc_iphi = 0.0577477; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_4GeV-ALL.root_Ana.root 0784 const double center_cemc_ieta = 0.00609392; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_4GeV-ALL.root_Ana.root 0785 const double center_hcalin_iphi = 0.354109; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_4GeV-ALL.root_Ana.root 0786 const double center_hcalin_ieta = -0.000273002; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_4GeV-ALL.root_Ana.root 0787 0788 /////////////////////////////////////////////// 0789 // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_2GeV-ALL.root_Ana.root 0790 /////////////////////////////////////////////// 0791 const double center_cemc_iphi = 0.170699; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_2GeV-ALL.root_Ana.root 0792 const double center_cemc_ieta = 0.00345903; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_2GeV-ALL.root_Ana.root 0793 const double center_hcalin_iphi = 0.668343; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_2GeV-ALL.root_Ana.root 0794 const double center_hcalin_ieta = -0.0130487; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_2GeV-ALL.root_Ana.root 0795 0796 /////////////////////////////////////////////// 0797 // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_2GeV-ALL.root_Ana.root 0798 /////////////////////////////////////////////// 0799 const double center_cemc_iphi = 0.0898377; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_2GeV-ALL.root_Ana.root 0800 const double center_cemc_ieta = 0.00445598; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_2GeV-ALL.root_Ana.root 0801 const double center_hcalin_iphi = 0.376608; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_2GeV-ALL.root_Ana.root 0802 const double center_hcalin_ieta = -0.0074984; // ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0_2GeV-ALL.root_Ana.root 0803 0804 /////////////////////////////////////////////// 0805 // Projection center based on ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root 0806 /////////////////////////////////////////////// 0807 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 0808 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 0809 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 0810 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 0811 /////////////////////////////////////////////// 0812 // Projection center based on ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root 0813 /////////////////////////////////////////////// 0814 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 0815 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 0816 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 0817 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 0818 /////////////////////////////////////////////// 0819 // Projection center based on ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root 0820 /////////////////////////////////////////////// 0821 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 0822 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 0823 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 0824 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 0825 /////////////////////////////////////////////// 0826 // Projection center based on ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root 0827 /////////////////////////////////////////////// 0828 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 0829 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 0830 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 0831 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 0832 /////////////////////////////////////////////// 0833 // Projection center based on ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root 0834 /////////////////////////////////////////////// 0835 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 0836 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 0837 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 0838 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 0839 /////////////////////////////////////////////// 0840 // Projection center based on ../../sPHENIX_work/production_analysis/single_particle/spacal1d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root 0841 /////////////////////////////////////////////// 0842 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 0843 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 0844 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 0845 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 0846 0847 /////////////////////////////////////////////// 0848 // Projection center based on ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_4GeV-ALL.root_Ana.root 0849 /////////////////////////////////////////////// 0850 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 0851 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 0852 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 0853 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 0854 0855 /////////////////////////////////////////////// 0856 // Projection center based on ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_4GeV-ALL.root_Ana.root 0857 /////////////////////////////////////////////// 0858 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 0859 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 0860 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 0861 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 0862 0863 /////////////////////////////////////////////// 0864 // Projection center based on ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_2GeV-ALL.root_Ana.root 0865 /////////////////////////////////////////////// 0866 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 0867 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 0868 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 0869 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 0870 0871 /////////////////////////////////////////////// 0872 // Projection center based on ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_2GeV-ALL.root_Ana.root 0873 /////////////////////////////////////////////// 0874 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 0875 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 0876 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 0877 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 0878 0879 /////////////////////////////////////////////// 0880 // Projection center based on ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0.90_8GeV-ALL.root_Ana.root 0881 /////////////////////////////////////////////// 0882 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 0883 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 0884 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 0885 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 0886 0887 /////////////////////////////////////////////// 0888 // Projection center based on ../../sPHENIX_work/production_analysis/single_particle/spacal2d/fieldmap/G4Hits_sPHENIX_pi-_eta0.90_8GeV-ALL.root_Ana.root 0889 /////////////////////////////////////////////// 0890 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 0891 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 0892 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 0893 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 0894 0895 0896 0897 0898 /////////////////////////////////////////////// 0899 // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_2GeV-ALL.root_Ana.root 0900 /////////////////////////////////////////////// 0901 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 0902 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 0903 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 0904 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 0905 0906 /////////////////////////////////////////////// 0907 // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_2GeV-ALL.root_Ana.root 0908 /////////////////////////////////////////////// 0909 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 0910 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 0911 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 0912 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 0913 0914 /////////////////////////////////////////////// 0915 // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_4GeV-ALL.root_Ana.root 0916 /////////////////////////////////////////////// 0917 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 0918 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 0919 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 0920 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 0921 0922 /////////////////////////////////////////////// 0923 // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_4GeV-ALL.root_Ana.root 0924 /////////////////////////////////////////////// 0925 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 0926 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 0927 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 0928 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 0929 0930 0931 /////////////////////////////////////////////// 0932 // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_8GeV-ALL.root_Ana.root 0933 /////////////////////////////////////////////// 0934 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 0935 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 0936 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 0937 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 0938 0939 /////////////////////////////////////////////// 0940 // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_8GeV-ALL.root_Ana.root 0941 /////////////////////////////////////////////// 0942 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 0943 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 0944 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 0945 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 0946 0947 0948 0949 0950 /////////////////////////////////////////////// 0951 // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_8GeV-ALL.root_Ana.root 0952 /////////////////////////////////////////////// 0953 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 0954 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 0955 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 0956 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 0957 /////////////////////////////////////////////// 0958 // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_8GeV-ALL.root_Ana.root 0959 /////////////////////////////////////////////// 0960 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 0961 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 0962 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 0963 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 0964 0965 /////////////////////////////////////////////// 0966 // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_4GeV-ALL.root_Ana.root 0967 /////////////////////////////////////////////// 0968 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 0969 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 0970 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 0971 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 0972 0973 /////////////////////////////////////////////// 0974 // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_4GeV-ALL.root_Ana.root 0975 /////////////////////////////////////////////// 0976 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 0977 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 0978 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 0979 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 0980 0981 /////////////////////////////////////////////// 0982 // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_e+_eta0.90_2GeV-ALL.root_Ana.root 0983 /////////////////////////////////////////////// 0984 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 0985 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 0986 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 0987 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 0988 0989 /////////////////////////////////////////////// 0990 // Projection center based on ../../sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal1d/fieldmap/G4Hits_sPHENIX_pi+_eta0.90_2GeV-ALL.root_Ana.root 0991 /////////////////////////////////////////////// 0992 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 0993 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 0994 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 0995 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 0996 0997 0998 0999 }
| [ Source navigation ] | [ Diff markup ] | [ Identifier search ] | [ general search ] |
|
This page was automatically generated by the 2.3.7 LXR engine. The LXR team |
|