sPhenix code displayed by LXR |
|
|
|
|||
File indexing completed on 2025-08-05 08:12:16
0001 // $Id: $ 0002 0003 /*! 0004 * \file Fun4All_EMCalLikelihood.C 0005 * \brief 0006 * \author Jin Huang <jhuang@bnl.gov> 0007 * \version $Revision: $ 0008 * \date $Date: $ 0009 */ 0010 0011 #include <cassert> 0012 0013 void 0014 Fun4All_EMCalLikelihood(const int nEvents = 1000, 0015 //Fun4All_EMCalLikelihood(const int nEvents = 100000000, // original 0016 TString base_dir = 0017 "../..//sPHENIX_work/production_analysis_cemc2x2/emcstudies/pidstudies/spacal2d/fieldmap/", 0018 TString pid = "e-", TString kine_config = "eta0_8GeV", TString ll_config = 0019 "Edep_Distribution_ll_sample") 0020 { 0021 const TString inputFile = base_dir + "/G4Hits_sPHENIX_" + pid + "_" 0022 + kine_config + "-ALL.root_Ana.root.lst"; 0023 // "/phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/G4Hits_sPHENIX_e-_eta0_2GeV-ALL.root_Ana.root.lst") 0024 0025 TString s_outputFile = inputFile; 0026 s_outputFile += "_EMCalLikelihood.root"; 0027 const char * outputFile = s_outputFile.Data(); 0028 0029 gSystem->Load("libfun4all.so"); 0030 gSystem->Load("libg4vertex.so"); 0031 gSystem->Load("libemcal_ana.so"); 0032 0033 Fun4AllServer *se = Fun4AllServer::instance(); 0034 se->Verbosity(1); 0035 // se->Verbosity(10); 0036 // just if we set some flags somewhere in this macro 0037 recoConsts *rc = recoConsts::instance(); 0038 rc->set_IntFlag("",0); 0039 0040 double center_cemc_iphi = 1000; 0041 double center_cemc_ieta = 1000; 0042 double center_hcalin_iphi = 1000; 0043 double center_hcalin_ieta = 1000; 0044 0045 double width = 0; 0046 double width_emcal_eta = 0; 0047 0048 int charge = 0; 0049 0050 if (pid == "e+" || pid == "pi+" || pid == "kaon+" || pid == "proton") 0051 charge = +1; 0052 else if (pid == "e-" || pid == "pi-" || pid == "kaon-" 0053 || pid == "anti_proton") 0054 charge = -1; 0055 if (charge == 0) 0056 { 0057 cout 0058 << "Fun4All_EMCalLikelihood - Fatal Error - invalid input particle type: " 0059 << pid << endl; 0060 assert(charge != 0); 0061 } 0062 0063 if (base_dir.Contains("spacal") && kine_config.Contains("eta0_")) 0064 { 0065 /////////////////////////////////////////////// 0066 // 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 0067 /////////////////////////////////////////////// 0068 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 0069 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 0070 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 0071 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 0072 0073 // width = 1.4; 0074 // width_emcal_eta = 1.4; 0075 width = 1.6; 0076 width_emcal_eta = 1.6; 0077 } 0078 else if (base_dir.Contains("spacal2d") && charge <0 0079 && kine_config.Contains("eta0.90_2GeV")) 0080 { 0081 0082 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 0083 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 0084 0085 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 0086 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 0087 0088 width = 1.6; 0089 width_emcal_eta = 1.6; 0090 } 0091 else if (base_dir.Contains("spacal2d") && charge <0 0092 && kine_config.Contains("eta0.90_4GeV")) 0093 { 0094 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 0095 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 0096 0097 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 0098 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 0099 0100 width = 1.6; 0101 width_emcal_eta = 1.6; 0102 } 0103 else if (base_dir.Contains("spacal2d") && charge <0 0104 && kine_config.Contains("eta0.90_8GeV")) 0105 { 0106 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 0107 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 0108 0109 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 0110 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 0111 0112 width = 1.6; 0113 width_emcal_eta = 1.6; 0114 } 0115 else if (base_dir.Contains("spacal1d") && charge <0 0116 && kine_config.Contains("eta0.90_2GeV")) 0117 { 0118 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 0119 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 0120 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 0121 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 0122 0123 width = 1.8; 0124 width_emcal_eta = 2.6; 0125 } 0126 else if (base_dir.Contains("spacal1d") && charge <0 0127 && kine_config.Contains("eta0.90_4GeV")) 0128 { 0129 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 0130 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 0131 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 0132 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 0133 0134 width = 1.8; 0135 width_emcal_eta = 2.6; 0136 } 0137 else if (base_dir.Contains("spacal1d") && charge <0 0138 && kine_config.Contains("eta0.90_8GeV")) 0139 { 0140 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 0141 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 0142 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 0143 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 0144 0145 width = 1.8; 0146 width_emcal_eta = 2.6; 0147 } 0148 else if (base_dir.Contains("spacal2d") && charge >0 0149 && kine_config.Contains("eta0.90_2GeV")) 0150 { 0151 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 0152 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 0153 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 0154 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 0155 0156 width = 1.6; 0157 width_emcal_eta = 1.6; 0158 } 0159 else if (base_dir.Contains("spacal2d") && charge >0 0160 && kine_config.Contains("eta0.90_4GeV")) 0161 { 0162 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 0163 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 0164 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 0165 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 0166 0167 0168 width = 1.6; 0169 width_emcal_eta = 1.6; 0170 } 0171 else if (base_dir.Contains("spacal2d") && charge >0 0172 && kine_config.Contains("eta0.90_8GeV")) 0173 { 0174 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 0175 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 0176 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 0177 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 0178 0179 0180 width = 1.6; 0181 width_emcal_eta = 1.6; 0182 } 0183 else if (base_dir.Contains("spacal1d") && charge >0 0184 && kine_config.Contains("eta0.90_2GeV")) 0185 { 0186 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 0187 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 0188 0189 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 0190 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 0191 0192 width = 1.8; 0193 width_emcal_eta = 2.6; 0194 } 0195 else if (base_dir.Contains("spacal1d") && charge >0 0196 && kine_config.Contains("eta0.90_4GeV")) 0197 { 0198 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 0199 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 0200 0201 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 0202 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 0203 0204 width = 1.8; 0205 width_emcal_eta = 2.6; 0206 } 0207 else if (base_dir.Contains("spacal1d") && charge >0 0208 && kine_config.Contains("eta0.90_8GeV")) 0209 { 0210 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 0211 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 0212 0213 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 0214 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 0215 0216 width = 1.8; 0217 width_emcal_eta = 2.6; 0218 } 0219 else 0220 { 0221 cout << "Error !!!!!! Unknown configuraiton"; 0222 0223 return; 0224 } 0225 0226 TH2F * h2_Edep_Distribution_e = NULL; 0227 TH2F * h2_Edep_Distribution_pi = NULL; 0228 0229 // load edep ll input distribution 0230 if (charge <0) 0231 { 0232 TFile * f = 0233 new TFile( 0234 base_dir + "/G4Hits_sPHENIX_e-_" + kine_config 0235 + "-ALL.root_Ana.root.lst_EMCalLikelihood.root_DrawEcal_Likelihood_" 0236 + ll_config + ".root"); 0237 0238 if (f->IsOpen()) 0239 { 0240 f->ls(); 0241 h2_Edep_Distribution_e = (TH2F *) f->GetObjectChecked( 0242 "h2_Edep_Distribution", "TH2F"); 0243 } 0244 0245 // TFile * f = 0246 // new TFile( 0247 // base_dir + "/G4Hits_sPHENIX_pi-_" + kine_config 0248 // + "-ALL.root_Ana.root.lst_EMCalLikelihood.root_DrawEcal_Likelihood_" 0249 // + ll_config + ".root"); 0250 TFile * f = 0251 new TFile( 0252 base_dir + "/G4Hits_sPHENIX_anti_proton_" + kine_config 0253 + "-ALL.root_Ana.root.lst_EMCalLikelihood.root_DrawEcal_Likelihood_" 0254 + ll_config + ".root"); 0255 0256 if (f->IsOpen()) 0257 { 0258 f->ls(); 0259 h2_Edep_Distribution_pi = (TH2F *) f->GetObjectChecked( 0260 "h2_Edep_Distribution", "TH2F"); 0261 } 0262 // assert(h2_Edep_Distribution_e); 0263 // assert(h2_Edep_Distribution_pi); 0264 } 0265 0266 else if (charge > 0) 0267 { 0268 TFile * f = 0269 new TFile( 0270 base_dir + "/G4Hits_sPHENIX_e+_" + kine_config 0271 + "-ALL.root_Ana.root.lst_EMCalLikelihood.root_DrawEcal_Likelihood_" 0272 + ll_config + ".root"); 0273 0274 if (f->IsOpen()) 0275 { 0276 f->ls(); 0277 h2_Edep_Distribution_e = (TH2F *) f->GetObjectChecked( 0278 "h2_Edep_Distribution", "TH2F"); 0279 } 0280 0281 // TFile * f = 0282 // new TFile( 0283 // base_dir + "/G4Hits_sPHENIX_pi-_" + kine_config 0284 // + "-ALL.root_Ana.root.lst_EMCalLikelihood.root_DrawEcal_Likelihood_" 0285 // + ll_config + ".root"); 0286 TFile * f = 0287 new TFile( 0288 base_dir + "/G4Hits_sPHENIX_pi+_" + kine_config 0289 + "-ALL.root_Ana.root.lst_EMCalLikelihood.root_DrawEcal_Likelihood_" 0290 + ll_config + ".root"); 0291 0292 if (f->IsOpen()) 0293 { 0294 f->ls(); 0295 h2_Edep_Distribution_pi = (TH2F *) f->GetObjectChecked( 0296 "h2_Edep_Distribution", "TH2F"); 0297 } 0298 // assert(h2_Edep_Distribution_e); 0299 // assert(h2_Edep_Distribution_pi); 0300 } 0301 0302 EMCalLikelihood * emcal_ana = new EMCalLikelihood( 0303 string(inputFile.Data()) + string("_hist.root")); 0304 emcal_ana->Verbosity(5); 0305 0306 emcal_ana->set_center_cemc_ieta(center_cemc_ieta); 0307 emcal_ana->set_center_cemc_iphi(center_cemc_iphi); 0308 emcal_ana->set_center_hcalin_ieta(center_hcalin_ieta); 0309 emcal_ana->set_center_hcalin_iphi(center_hcalin_iphi); 0310 0311 // ------------------------------------------------------------------- 0312 // No CEMC ganging readout 0313 // ------------------------------------------------------------------- 0314 // emcal_ana->set_width_cemc_ieta(width_emcal_eta); 0315 // emcal_ana->set_width_cemc_iphi(width); 0316 // emcal_ana->set_width_hcalin_ieta(width); 0317 // emcal_ana->set_width_hcalin_iphi(width); 0318 0319 // ------------------------------------------------------------------- 0320 // 2x2 CEMC ganging readout, cluster size x (2x2) 0321 // ------------------------------------------------------------------- 0322 // emcal_ana->do_ganging(2,2); 0323 // emcal_ana->set_width_cemc_ieta(width_emcal_eta*2); 0324 // emcal_ana->set_width_cemc_iphi(width*2); 0325 // emcal_ana->set_width_hcalin_ieta(width*2); 0326 // emcal_ana->set_width_hcalin_iphi(width*2); 0327 0328 // ------------------------------------------------------------------- 0329 // 2x2 CEMC ganging readout, cluster size x (1.4x1.4) 0330 // ------------------------------------------------------------------- 0331 // emcal_ana->do_ganging(2,2); 0332 // emcal_ana->set_width_cemc_ieta(width_emcal_eta*1.4); 0333 // emcal_ana->set_width_cemc_iphi(width*1.4); 0334 // emcal_ana->set_width_hcalin_ieta(width*1.4); 0335 // emcal_ana->set_width_hcalin_iphi(width*1.4); 0336 0337 // ------------------------------------------------------------------- 0338 // 2x2 CEMC ganging readout, cluster size x (1.2x1.2) 0339 // ------------------------------------------------------------------- 0340 emcal_ana->do_ganging(2,2); 0341 emcal_ana->set_width_cemc_ieta(width_emcal_eta*1.2); 0342 emcal_ana->set_width_cemc_iphi(width*1.2); 0343 emcal_ana->set_width_hcalin_ieta(width); 0344 emcal_ana->set_width_hcalin_iphi(width); 0345 0346 if (h2_Edep_Distribution_e) 0347 emcal_ana->set_h2_Edep_Distribution_e(h2_Edep_Distribution_e); 0348 if (h2_Edep_Distribution_pi) 0349 emcal_ana->set_h2_Edep_Distribution_pi(h2_Edep_Distribution_pi); 0350 0351 0352 se->registerSubsystem(emcal_ana); 0353 0354 //-------------- 0355 // IO management 0356 //-------------- 0357 0358 // Hits file 0359 Fun4AllInputManager *hitsin = new Fun4AllDstInputManager("DSTin"); 0360 // hitsin->fileopen(inputFile); 0361 hitsin->AddListFile(inputFile.Data()); 0362 se->registerInputManager(hitsin); 0363 0364 Fun4AllDstOutputManager *out = new Fun4AllDstOutputManager("DSTOUT", 0365 outputFile); 0366 out->AddNode("Sync"); 0367 out->AddNode("UpsilonPair"); 0368 out->AddNode("EMCalTrk"); 0369 out->AddNode("GlobalVertexMap"); 0370 se->registerOutputManager(out); 0371 0372 gSystem->ListLibraries(); 0373 0374 // return; 0375 0376 se->run(nEvents); 0377 // se->dumpHistos(string(inputFile) + string("_hist.root"), "recreate"); 0378 0379 //----- 0380 // Exit 0381 //----- 0382 gSystem->Exec("ps -o sid,ppid,pid,user,comm,vsize,rssize,time"); 0383 0384 se->End(); 0385 0386 std::cout << "All done" << std::endl; 0387 delete se; 0388 gSystem->Exit(0); 0389 } 0390
| [ Source navigation ] | [ Diff markup ] | [ Identifier search ] | [ general search ] |
|
This page was automatically generated by the 2.3.7 LXR engine. The LXR team |
|