void Jet_FwdReco(int verbosity = 0) { gSystem->Load("libg4jets.so"); Fun4AllServer *se = Fun4AllServer::instance(); // truth particle level jets JetReco *truthjetreco = new JetReco(); truthjetreco->add_input(new TruthJetInput(Jet::PARTICLE)); //truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.2),"AntiKt_Truth_r02"); truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.3),"AntiKt_Truth_r03"); // truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.4),"AntiKt_Truth_r04"); // truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.5),"AntiKt_Truth_r05"); // truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.6),"AntiKt_Truth_r06"); // truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.7),"AntiKt_Truth_r07"); // truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.8),"AntiKt_Truth_r08"); truthjetreco->set_algo_node("ANTIKT"); truthjetreco->set_input_node("TRUTH"); truthjetreco->Verbosity(verbosity); se->registerSubsystem(truthjetreco); // tower jets JetReco *towerjetreco = new JetReco(); towerjetreco->add_input(new TowerJetInput(Jet::FEMC_TOWER)); towerjetreco->add_input(new TowerJetInput(Jet::FHCAL_TOWER)); // towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.2),"AntiKt_Tower_r02"); towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.3),"AntiKt_Tower_r03"); // towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.4),"AntiKt_Tower_r04"); // towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.5),"AntiKt_Tower_r05"); // towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.6),"AntiKt_Tower_r06"); // towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.7),"AntiKt_Tower_r07"); // towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.8),"AntiKt_Tower_r08"); towerjetreco->set_algo_node("ANTIKT"); towerjetreco->set_input_node("TOWER"); towerjetreco->Verbosity(verbosity); se->registerSubsystem(towerjetreco); return; }
void Jet_Reco(int verbosity = 0) { gSystem->Load("libg4jets.so"); Fun4AllServer *se = Fun4AllServer::instance(); // truth particle level jets JetReco *truthjetreco = new JetReco(); truthjetreco->add_input(new TruthJetInput(Jet::PARTICLE)); truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.2),"AntiKt_Truth_r0.2"); truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.3),"AntiKt_Truth_r0.3"); truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.4),"AntiKt_Truth_r0.4"); truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.5),"AntiKt_Truth_r0.5"); truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.6),"AntiKt_Truth_r0.6"); truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.7),"AntiKt_Truth_r0.7"); truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.8),"AntiKt_Truth_r0.8"); truthjetreco->set_algo_node("ANTIKT"); truthjetreco->set_input_node("TRUTH"); truthjetreco->Verbosity(verbosity); se->registerSubsystem(truthjetreco); // tower jets JetReco *towerjetreco = new JetReco(); towerjetreco->add_input(new TowerJetInput(Jet::CEMC_TOWER)); towerjetreco->add_input(new TowerJetInput(Jet::HCALIN_TOWER)); towerjetreco->add_input(new TowerJetInput(Jet::HCALOUT_TOWER)); towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.2),"AntiKt_Tower_r0.2"); towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.3),"AntiKt_Tower_r0.3"); towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.4),"AntiKt_Tower_r0.4"); towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.5),"AntiKt_Tower_r0.5"); towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.6),"AntiKt_Tower_r0.6"); towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.7),"AntiKt_Tower_r0.7"); towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.8),"AntiKt_Tower_r0.8"); towerjetreco->set_algo_node("ANTIKT"); towerjetreco->set_input_node("TOWER"); towerjetreco->Verbosity(verbosity); se->registerSubsystem(towerjetreco); // cluster jets JetReco *clusterjetreco = new JetReco(); clusterjetreco->add_input(new ClusterJetInput(Jet::CEMC_CLUSTER)); clusterjetreco->add_input(new ClusterJetInput(Jet::HCALIN_CLUSTER)); clusterjetreco->add_input(new ClusterJetInput(Jet::HCALOUT_CLUSTER)); clusterjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.2),"AntiKt_Cluster_r0.2"); clusterjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.3),"AntiKt_Cluster_r0.3"); clusterjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.4),"AntiKt_Cluster_r0.4"); clusterjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.5),"AntiKt_Cluster_r0.5"); clusterjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.6),"AntiKt_Cluster_r0.6"); clusterjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.7),"AntiKt_Cluster_r0.7"); clusterjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.8),"AntiKt_Cluster_r0.8"); clusterjetreco->set_algo_node("ANTIKT"); clusterjetreco->set_input_node("CLUSTER"); clusterjetreco->Verbosity(verbosity); se->registerSubsystem(clusterjetreco); // track jets JetReco *trackjetreco = new JetReco(); trackjetreco->add_input(new TrackJetInput(Jet::TRACK)); trackjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.2),"AntiKt_Track_r0.2"); trackjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.3),"AntiKt_Track_r0.3"); trackjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.4),"AntiKt_Track_r0.4"); trackjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.5),"AntiKt_Track_r0.5"); trackjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.6),"AntiKt_Track_r0.6"); trackjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.7),"AntiKt_Track_r0.7"); trackjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.8),"AntiKt_Track_r0.8"); trackjetreco->set_algo_node("ANTIKT"); trackjetreco->set_input_node("TRACK"); trackjetreco->Verbosity(verbosity); se->registerSubsystem(trackjetreco); return; }