macros/nano_analysis/nano_DST_analysis_efficiency.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_efficiency(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_efficiency_30/diphoton_minimal_";
0013   std::string outputfiletemplate = "analysis_complete_ana509_01312026_backward_efficiency_30_pt_g3_0mrad/analysis_";
0014
0015   Fun4AllServer *se = Fun4AllServer::instance();
0016   //se->Verbosity(0);
0017
0018   recoConsts *rc = recoConsts::instance();
0019   rc->set_StringFlag("CDB_GLOBALTAG", "ProdA_2024"); // Required for the spin DB
0020   rc->set_uint64Flag("TIMESTAMP", 48746); // I am not using the TIMESTAMP anyway.
0021
0022   AnNeutralMeson_nano *AN = new AnNeutralMeson_nano("AnNeutralMeson_nano", inputlistname, inputfiletemplate, outputfiletemplate);
0023   AN->set_store_tree(true);
0024   //AN->set_forward_cut(true);
0025   AN->set_backward_cut(true);
0026   AN->set_ptcut(3.0, 1000.0);
0027   se->registerSubsystem(AN);
0028   se->run(1);
0029   se->End();
0030   delete se;
0031
0032   gSystem->Exit(0);
0033 }