macros/micro_analysis/micro_DST_analysis.C

0001 #include <anneutral/AnNeutralMeson_micro.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(const int runnumber = 47289)
0010 {
0011   std::string inputfilename = "/sphenix/tg/tg01/coldqcd/vmahaut/AnNeutralMeson/combine_out_run_ana450_2024p009/OUTDIPHOTON_" + std::to_string(runnumber) + ".bin";
0012   std::string outputfilename = "/sphenix/u/virgilemahaut/work/analysis/AnNeutralMeson/macros/comparison_micro_DST_analysis/analysis_per_run_allruns/analysis_" + std::to_string(runnumber) + ".root";
0013   std::string outputfiletreename = "/sphenix/u/virgilemahaut/work/analysis/AnNeutralMeson/macros/comparison_micro_DST_analysis/analysis_per_run_allruns/diphoton_minimal_" + std::to_string(runnumber) + ".root";
0014   //std::string inputfilename = "../combine_out_run_48746_05_1000/OUTDIPHOTON_" + std::to_string(runnumber) + ".root";
0015   //std::string outputfilename = "analysis_per_run_48746_asymmetry/analysis_" + std::to_string(runnumber) + ".root";
0016
0017   Fun4AllServer *se = Fun4AllServer::instance();
0018   se->Verbosity(1000);
0019
0020   recoConsts *rc = recoConsts::instance();
0021   rc->set_StringFlag("CDB_GLOBALTAG", "ProdA_2024");
0022   rc->set_uint64Flag("TIMESTAMP", runnumber);
0023
0024   AnNeutralMeson_micro *AN = new AnNeutralMeson_micro("AnNeutralMeson_micro", runnumber, inputfilename, outputfilename,outputfiletreename);
0025   std::vector<float> chi2_cuts = {1000}; // {1000, 6, 4, 2};
0026   //std::vector<float> chi2_cuts = {1000, 6, 4, 2};
0027   std::vector<float> ecore_cuts = {1.0}; // {0.5, 1.0, 1.5};
0028   //std::vector<float> ecore_cuts = {0.5, 1.0, 1.5};
0029   float alpha_cut = 0.7;
0030   float pT_cut = 1.0;
0031   AN->set_chi2cut(chi2_cuts);
0032   AN->set_ecorecut(ecore_cuts);
0033   AN->set_alphacut(alpha_cut);
0034   AN->set_ptcut(pT_cut);
0035   AN->set_triggerselection(true);
0036   se->registerSubsystem(AN);
0037
0038   se->run(1);
0039   se->End();
0040   delete se;
0041
0042   gSystem->Exit(0);
0043 }