macros/micro_analysis/micro_DST_analysis_new_formatted.C

0001 #include <anneutral/AnNeutralMeson_micro_dst_formatted.h>
0002 #include <fun4all/Fun4AllServer.h>
0003 #include <phool/recoConsts.h>
0004
0005 R__LOAD_LIBRARY(libffamodules.so)
0006 R__LOAD_LIBRARY(libfun4all.so)
0007 R__LOAD_LIBRARY(libanneutral.so)
0008
0009 void micro_DST_analysis_new_formatted(const int runnumber = 47289)
0010 {
0011   std::string inputfilename = "/sphenix/tg/tg01/coldqcd/vmahaut/AnNeutralMeson/combine_out_run_ana450_2024p009/root_common/preliminary/OUTDIPHOTON_PRELIMINARY_COMMON_" + std::to_string(runnumber) + ".root";
0012   //std::string inputfilename = "/sphenix/tg/tg01/coldqcd/vmahaut/AnNeutralMeson/combine_out_run_ana450_2024p009/root_ordered_triggered/OUTDIPHOTON_ORDERED_TRIGGERED_" + std::to_string(runnumber) + ".root";
0013   std::string outputfilename = "/sphenix/u/virgilemahaut/work/analysis/AnNeutralMeson/macros/comparison_micro_DST_analysis/analysis_per_run_allruns_new_formatted/analysis_" + std::to_string(runnumber) + ".root";
0014   std::string outputfiletreename = "/sphenix/u/virgilemahaut/work/analysis/AnNeutralMeson/macros/comparison_micro_DST_analysis/analysis_per_run_allruns_new_formatted/diphoton_minimal_" + std::to_string(runnumber) + ".root";
0015   //std::string inputfilename = "../combine_out_run_48746_05_1000/OUTDIPHOTON_" + std::to_string(runnumber) + ".root";
0016   //std::string outputfilename = "analysis_per_run_48746_asymmetry/analysis_" + std::to_string(runnumber) + ".root";
0017
0018   Fun4AllServer *se = Fun4AllServer::instance();
0019   se->Verbosity(1000);
0020
0021   recoConsts *rc = recoConsts::instance();
0022   rc->set_StringFlag("CDB_GLOBALTAG", "ProdA_2024");
0023   rc->set_uint64Flag("TIMESTAMP", runnumber);
0024
0025   AnNeutralMeson_micro_dst_formatted *AN = new AnNeutralMeson_micro_dst_formatted("AnNeutralMeson_micro_formatted", runnumber, inputfilename, outputfilename,outputfiletreename);
0026   std::vector<float> chi2_cuts = {1000}; // {1000, 6, 4, 2};
0027   //std::vector<float> chi2_cuts = {1000, 6, 4, 2};
0028   std::vector<float> ecore_cuts = {1.0}; // {0.5, 1.0, 1.5};
0029   //std::vector<float> ecore_cuts = {0.5, 1.0, 1.5};
0030   float alpha_cut = 0.7;
0031   float pT_cut = 1.0;
0032   AN->set_chi2cut(chi2_cuts);
0033   AN->set_ecorecut(ecore_cuts);
0034   AN->set_alphacut(alpha_cut);
0035   AN->set_ptcut(pT_cut);
0036   AN->set_photon_trigger_bit_requirement(true);
0037   AN->set_photon_trigger_efficiency_matching_requirement(true);
0038   AN->set_photon_trigger_emulator_matching_requirement(false);
0039   AN->set_event_mixing(false);
0040   se->registerSubsystem(AN);
0041
0042   se->run(1);
0043   se->End();
0044   delete se;
0045
0046   gSystem->Exit(0);
0047 }