LaserClusterQA/macros/Fun4All_laserQA.C

0001 /*
0002  * This macro shows a minimum working example of running the tracking
0003  * hit unpackers with some basic seeding algorithms to try to put together
0004  * tracks. There are some analysis modules run at the end which package
0005  * hits, clusters, and clusters on tracks into trees for analysis.
0006  */
0007
0008 #include <fun4all/Fun4AllUtils.h>
0009 #include <GlobalVariables.C>
0010 //#include <G4_ActsGeom.C>
0011 //#include <G4_Global.C>
0012 //#include <G4_Magnet.C>
0013 //#include <G4_Mbd.C>
0014 //#include <Trkr_Clustering.C>
0015 //#include <Trkr_LaserClustering.C>
0016 //#include <Trkr_RecoInit.C>
0017 //#include <Trkr_TpcReadoutInit.C>
0018
0019 #include <ffamodules/CDBInterface.h>
0020 #include <fun4all/Fun4AllDstInputManager.h>
0021 #include <fun4all/Fun4AllDstOutputManager.h>
0022 #include <fun4all/Fun4AllInputManager.h>
0023 #include <fun4all/Fun4AllOutputManager.h>
0024 #include <fun4all/Fun4AllRunNodeInputManager.h>
0025 #include <fun4all/Fun4AllServer.h>
0026
0027 #include <phool/recoConsts.h>
0028
0029 //#include <tpc/LaserEventIdentifier.h>
0030
0031 #include <laserclusterqa/laserQA.h>
0032
0033 #include <stdio.h>
0034 R__LOAD_LIBRARY(libfun4all.so)
0035 R__LOAD_LIBRARY(libffamodules.so)
0036 R__LOAD_LIBRARY(libtpc.so)
0037 R__LOAD_LIBRARY(libtpccalib.so)
0038 R__LOAD_LIBRARY(libLaserClusterQA.so)
0039 void Fun4All_laserQA(
0040     const int nEvents = 10,
0041     const std::string filelist = "/sphenix/lustre01/sphnxpro/commissioning/slurp/tpcbeam/run_00041900_00042000/",
0042     const std::string outdir = "/sphenig/tg/tg01/jets/bkimelman/",
0043     const std::string prod = "ana441_2024p007",
0044     const std::string type = "run2pp")
0045
0046 {
0047
0048   auto se = Fun4AllServer::instance();
0049   se->Verbosity(2);
0050
0051   //TRACKING::tpc_zero_supp = true;
0052   //TRACKING::tpc_baseline_corr = true;
0053   //Enable::MVTX_APPLYMISALIGNMENT = true;
0054   //ACTSGEOM::mvtx_applymisalignment = Enable::MVTX_APPLYMISALIGNMENT;
0055
0056   std::ifstream ifs(filelist);
0057   std::string filepath;
0058
0059   auto rc = recoConsts::instance();
0060
0061   Enable::CDB = true;
0062   rc->set_StringFlag("CDB_GLOBALTAG", "ProdA_2024");
0063
0064   int i=0;
0065   int runnumber = 0;
0066   int segment = 0;
0067   while(std::getline(ifs,filepath))
0068   {
0069     std::cout << "Adding DST with filepath: " << filepath << std::endl;
0070     if(i==0)
0071     {
0072       std::pair<int, int> runseg = Fun4AllUtils::GetRunSegment(filepath);
0073       runnumber = runseg.first;
0074       segment = runseg.second;
0075       rc->set_IntFlag("RUNNUMBER",runnumber);
0076       rc->set_uint64Flag("TIMESTAMP",runnumber);
0077     }
0078
0079     //std::string inputname = "InputManager" + std::to_string(i);
0080     //auto hitsin = new Fun4AllDstInputManager(inputname);
0081     //hitsin->fileopen(filepath);
0082     //hitsin->AddFile(filepath);
0083     //se->registerInputManager(hitsin);
0084     i++;
0085   }
0086
0087   auto hitsin = new Fun4AllDstInputManager("InputManager");
0088   hitsin->AddListFile(filelist);
0089   se->registerInputManager(hitsin);
0090
0091   /*
0092   std::string geofile = CDBInterface::instance()->getUrl("Tracking_Geometry");
0093   Fun4AllRunNodeInputManager *ingeo = new Fun4AllRunNodeInputManager("GeoIn");
0094   ingeo->AddFile(geofile);
0095   se->registerInputManager(ingeo);
0096   */
0097
0098   //  G4TPC::tpc_drift_velocity_reco = (8.0 / 1000) * 107.0 / 105.0;
0099   //G4TPC::ENABLE_MODULE_EDGE_CORRECTIONS = true;
0100   //to turn on the default static corrections, enable the two lines below
0101   //G4TPC::ENABLE_STATIC_CORRECTIONS = true;
0102
0103   //if(runnumber == 53098)
0104   //{
0105   //G4TPC::ENABLE_STATIC_CORRECTIONS = false;
0106   //}
0107
0108   //G4TPC::USE_PHI_AS_RAD_STATIC_CORRECTIONS=false;
0109
0110
0111   //G4TPC::ENABLE_AVERAGE_CORRECTIONS = false;
0112   //G4TPC::average_correction_filename = std::string(Form("/sphenix/tg/tg01/jets/bkimelman/BenProduction/Feb25_2025/Laminations_run2pp_ana466_2024p012_v001-%08d.root",runnumber));
0113   //G4TPC::USE_PHI_AS_RAD_AVERAGE_CORRECTIONS=false;
0114   //G4TPC::average_correction_interpolate = false;
0115
0116   //G4TPC::DISTORTIONS_USE_PHI_AS_RADIANS = false;
0117   //G4TPC::ENABLE_CENTRAL_MEMBRANE_CLUSTERING = true;
0118
0119   //TRACKING::reco_tpc_maxtime_sample = 1023;
0120   //TRACKING::tpc_zero_supp = true;
0121
0122   //TrackingInit();
0123
0124
0125   laserQA *lcQA = new laserQA("laserClusterQA");
0126   lcQA->setOutfile(Form("%s/laserHitQA_%s_%s-%08d.root",outdir.c_str(),type.c_str(),prod.c_str(),runnumber));
0127   se->registerSubsystem(lcQA);
0128
0129   se->run(nEvents);
0130   //se->run(-1);
0131   se->End();
0132   se->PrintTimer();
0133
0134   delete se;
0135   std::cout << "Finished" << std::endl;
0136   gSystem->Exit(0);
0137 }