Shift_emcal_g4hits/macros/Fun4All_GeoExtract.C

0001 #pragma once
0002 #if ROOT_VERSION_CODE >= ROOT_VERSION(6,00,0)
0003
0004 #include <fun4all/SubsysReco.h>
0005 #include <fun4all/Fun4AllServer.h>
0006 #include <fun4all/Fun4AllInputManager.h>
0007 #include <fun4all/Fun4AllDstInputManager.h>
0008
0009 #include <fun4all/Fun4AllDstOutputManager.h>
0010 #include <fun4all/Fun4AllOutputManager.h>
0011 #include <fun4all/Fun4AllServer.h>
0012
0013 #include <phool/PHRandomSeed.h>
0014 #include <phool/recoConsts.h>
0015
0016
0017
0018 R__LOAD_LIBRARY(libfun4all.so)
0019
0020
0021 #endif
0022
0023
0024 void Fun4All_GeoExtract(const char *filegeometry = "G4Hits_sHijing_pAu_0_10fm-0000000006-00000.root",
0025             const char * oufilename = "tree_test.root"
0026             )
0027 {
0028   gSystem->Load("libg4dst.so");
0029
0030   Fun4AllServer *se = Fun4AllServer::instance();
0031   int verbosity = 0;
0032
0033   se->Verbosity(verbosity);
0034   recoConsts *rc = recoConsts::instance();
0035
0036  //===============
0037   // conditions DB flags
0038   //===============
0039
0040   rc->set_StringFlag("CDB_GLOBALTAG","MDC2");
0041   rc->set_uint64Flag("TIMESTAMP",6);
0042
0043
0044
0045
0046   Fun4AllInputManager *intrue = new Fun4AllDstInputManager("DSTgeo");
0047   intrue->AddFile(filegeometry);
0048   se->registerInputManager(intrue);
0049
0050
0051   Fun4AllDstOutputManager *out = new Fun4AllDstOutputManager("DSTOUT", "updated_geo.root");
0052   out->AddRunNode("CYLINDERGEOM_CEMC");
0053   out->AddRunNode("RunHeader");
0054   out->AddNode("just_pick_an_invalid_name");
0055   se->registerOutputManager(out);
0056
0057
0058   se->run(1);
0059   se->End();
0060
0061   gSystem->Exit(0);
0062   return 0;
0063
0064 }