macros/nano_analysis/nano_DST_analysis_mbd_4_phenix.C

0001 #include <anneutral/AnNeutralMeson_nano.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 nano_DST_analysis_mbd_4_phenix(const int seednumber = 0)
0010 {
0011   std::string inputlistname = "inputruns.txt";
0012   std::string inputfiletemplate = "/sphenix/u/virgilemahaut/work/analysis/AnNeutralMeson/macros/comparison_micro_DST_analysis/nano_analysis/trees_complete_ana509_01312026_MBD/diphoton_minimal_";
0013   std::string outputfolder = "analysis_complete_ana509_01312026_phenix_MBD_pt_l4_0mrad/";
0014   gSystem->Exec(("mkdir -p " + outputfolder).c_str());
0015   std::string outputfiletemplate = outputfolder + "/analysis_";
0016   std::string outputbunchtemplate = outputfolder + "/fast_analysis_";
0017
0018   Fun4AllServer *se = Fun4AllServer::instance();
0019
0020   recoConsts *rc = recoConsts::instance();
0021   rc->set_StringFlag("CDB_GLOBALTAG", "ProdA_2024"); // Required for the spin DB
0022   rc->set_uint64Flag("TIMESTAMP", 48746); // I am not using the TIMESTAMP anyway.
0023
0024   AnNeutralMeson_nano *AN = new AnNeutralMeson_nano("AnNeutralMeson_nano", inputlistname, inputfiletemplate, outputfiletemplate);
0025   AN->set_store_bunch_yields(true, outputbunchtemplate);
0026   AN->set_trigger_mbd(true);
0027   AN->set_ptcut(1.0, 4.0);
0028   AN->set_phenix_cut(true);
0029   se->registerSubsystem(AN);
0030   se->run(1);
0031   se->End();
0032   delete se;
0033
0034   gSystem->Exit(0);
0035 }