JetRTrack/macros/Fun4All_RTrack.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 #include <g4jets/FastJetAlgo.h>
0017 #include <g4jets/JetReco.h>
0018 #include <g4jets/TowerJetInput.h>
0019 #include <g4jets/TruthJetInput.h>
0020
0021 #include <g4centrality/PHG4CentralityReco.h>
0022
0023 #include <jetbackground/CopyAndSubtractJets.h>
0024 #include <jetbackground/DetermineTowerBackground.h>
0025 #include <jetbackground/FastJetAlgoSub.h>
0026 #include <jetbackground/RetowerCEMC.h>
0027 #include <jetbackground/SubtractTowers.h>
0028 #include <jetbackground/SubtractTowersCS.h>
0029
0030
0031 #include <G4_Bbc.C>
0032 #include <G4_CaloTrigger.C>
0033 #include <G4_Centrality.C>
0034 #include <G4_DSTReader.C>
0035 #include <G4_Global.C>
0036 #include <G4_Input.C>
0037 #include <G4_KFParticle.C>
0038 #include <G4_ParticleFlow.C>
0039 #include <G4_Production.C>
0040 #include <G4_TopoClusterReco.C>
0041 #include <G4_Tracking.C>
0042 #include <G4_User.C>
0043 #include <QA.C>
0044 #include "HIJetReco.C"
0045
0046 // #include <G4_CEmc_Spacal.C>
0047 #include <jetrtrack/jetrtrack.h>
0048
0049 R__LOAD_LIBRARY(libfun4all.so)
0050 R__LOAD_LIBRARY(libg4jets.so)
0051 R__LOAD_LIBRARY(libjetbackground.so)
0052 R__LOAD_LIBRARY(libjetrtrack.so)
0053 R__LOAD_LIBRARY(libg4centrality.so)
0054 R__LOAD_LIBRARY(libg4dst.so)
0055
0056 #endif
0057
0058
0059 void Fun4All_RTrack(const char *filetruthjet = "DST_TRACKS_pythia8_Jet30-0000000062-00000.root",
0060                      const char *filecalo = "DST_CALO_CLUSTER_pythia8_Jet30-0000000062-00000.root",
0061             const char *filetruthjets = "DST_TRUTH_JET_pythia8_Jet30-0000000062-00000.root",
0062             const char *filetruth = "DST_TRUTH_pythia8_Jet30-0000000062-00000.root",
0063             const char * oufilename = "tree_test.root"
0064             )
0065 {
0066
0067
0068   Fun4AllServer *se = Fun4AllServer::instance();
0069   int verbosity = 0;
0070
0071   se->Verbosity(verbosity);
0072   recoConsts *rc = recoConsts::instance();
0073
0074
0075
0076
0077   PHG4CentralityReco *cent = new PHG4CentralityReco();
0078   cent->Verbosity(0);
0079   cent->GetCalibrationParameters().ReadFromFile("centrality", "xml", 0, 0, string(getenv("CALIBRATIONROOT")) + string("/Centrality/"));
0080   se->registerSubsystem( cent );
0081
0082   HIJetReco(); // Performs the officially reccomended jet reconstruction
0083
0084
0085   jetrtrack *myJetVal = new jetrtrack("analysis");
0086   myJetVal->SetOutputFileName(oufilename);
0087   se->registerSubsystem(myJetVal);
0088
0089
0090
0091   Fun4AllInputManager *intrue = new Fun4AllDstInputManager("DSTcalocluster");
0092   intrue->AddFile(filetruthjet);
0093   se->registerInputManager(intrue);
0094
0095   Fun4AllInputManager *in2 = new Fun4AllDstInputManager("DSTcalo");
0096   in2->AddFile(filecalo);
0097   se->registerInputManager(in2);
0098
0099   Fun4AllInputManager *in3 = new Fun4AllDstInputManager("DSTtruthjets");
0100   in3->AddFile(filetruthjets);
0101   se->registerInputManager(in3);
0102
0103   Fun4AllInputManager *in4 = new Fun4AllDstInputManager("DSTtrackseeds");
0104   in4->AddFile(filetruth);
0105   se->registerInputManager(in4);
0106
0107
0108   se->run();
0109   se->End();
0110
0111   gSystem->Exit(0);
0112   return 0;
0113
0114 }