macros/leptoquarks_fun4all/G4_FwdJets.C

0001
0002 void Jet_FwdReco(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   // tower jets
0024   JetReco *towerjetreco = new JetReco();
0025   towerjetreco->add_input(new TowerJetInput(Jet::FEMC_TOWER));
0026   towerjetreco->add_input(new TowerJetInput(Jet::FHCAL_TOWER));
0027   towerjetreco->add_input(new TowerJetInput(Jet::CEMC_TOWER));
0028   towerjetreco->add_input(new TowerJetInput(Jet::HCALIN_TOWER));
0029   towerjetreco->add_input(new TowerJetInput(Jet::HCALOUT_TOWER));
0030   // towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.2),"AntiKt_Tower_r02");
0031   towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.3),"AntiKt_Tower_r03");
0032   // towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.4),"AntiKt_Tower_r04");
0033    towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.5),"AntiKt_Tower_r05");
0034   // towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.6),"AntiKt_Tower_r06");
0035    towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.7),"AntiKt_Tower_r07");
0036   // towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.8),"AntiKt_Tower_r08");
0037   towerjetreco->set_algo_node("ANTIKT");
0038   towerjetreco->set_input_node("TOWER");
0039   towerjetreco->Verbosity(verbosity);
0040   se->registerSubsystem(towerjetreco);
0041
0042   return;
0043 }
0044
0045 void Jet_FwdEval(std::string outfilename = "g4fwdjets_eval.root",
0046           int verbosity = 0) {
0047
0048   gSystem->Load("libg4eval.so");
0049
0050   Fun4AllServer *se = Fun4AllServer::instance();
0051
0052   JetEvaluator* eval = new JetEvaluator("JETEVALUATOR",
0053                     "AntiKt_Tower_r05",
0054                     "AntiKt_Truth_r05",
0055                     outfilename);
0056   eval->Verbosity(verbosity);
0057   se->registerSubsystem(eval);
0058
0059   return;
0060 }