PhotonJet/example_macros/G4_TruthJets.C

0001
0002 void Jet_Reco(int verbosity = 0) {
0003
0004   gSystem->Load("libg4jets.so");
0005
0006   Fun4AllServer *se = Fun4AllServer::instance();
0007
0008   // truth particle level jets
0009   JetReco *truthjetreco = new JetReco();
0010   truthjetreco->add_input(new TruthJetInput(Jet::PARTICLE));
0011   truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.2),"AntiKt_Truth_r02");
0012   truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.3),"AntiKt_Truth_r03");
0013   truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.4),"AntiKt_Truth_r04");
0014   truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.5),"AntiKt_Truth_r05");
0015   truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.6),"AntiKt_Truth_r06");
0016   truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.7),"AntiKt_Truth_r07");
0017   truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.8),"AntiKt_Truth_r08");
0018   truthjetreco->set_algo_node("ANTIKT");
0019   truthjetreco->set_input_node("TRUTH");
0020   truthjetreco->Verbosity(verbosity);
0021   se->registerSubsystem(truthjetreco);
0022
0023
0024   return;
0025 }
0026
0027 void Jet_Eval(std::string outfilename = "g4jets_eval.root",
0028           int verbosity = 0) {
0029
0030   gSystem->Load("libg4eval.so");
0031
0032   Fun4AllServer *se = Fun4AllServer::instance();
0033
0034   JetEvaluator* eval = new JetEvaluator("JETEVALUATOR",
0035                     "AntiKt_Tower_r03",
0036                     "AntiKt_Truth_r03",
0037                     outfilename);
0038   eval->Verbosity(verbosity);
0039   se->registerSubsystem(eval);
0040
0041   return;
0042 }