macros/leptoquarks_fun4all/G4_Jets.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   // tower jets
0024   JetReco *towerjetreco = new JetReco();
0025   towerjetreco->add_input(new TowerJetInput(Jet::CEMC_TOWER));
0026   towerjetreco->add_input(new TowerJetInput(Jet::HCALIN_TOWER));
0027   towerjetreco->add_input(new TowerJetInput(Jet::HCALOUT_TOWER));
0028   towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.2),"AntiKt_Tower_r02");
0029   towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.3),"AntiKt_Tower_r03");
0030   towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.4),"AntiKt_Tower_r04");
0031   towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.5),"AntiKt_Tower_r05");
0032   towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.6),"AntiKt_Tower_r06");
0033   towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.7),"AntiKt_Tower_r07");
0034   towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.8),"AntiKt_Tower_r08");
0035   towerjetreco->set_algo_node("ANTIKT");
0036   towerjetreco->set_input_node("TOWER");
0037   towerjetreco->Verbosity(verbosity);
0038   se->registerSubsystem(towerjetreco);
0039
0040   // cluster jets
0041   JetReco *clusterjetreco = new JetReco();
0042   clusterjetreco->add_input(new ClusterJetInput(Jet::CEMC_CLUSTER));
0043   clusterjetreco->add_input(new ClusterJetInput(Jet::HCALIN_CLUSTER));
0044   clusterjetreco->add_input(new ClusterJetInput(Jet::HCALOUT_CLUSTER));
0045   clusterjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.2),"AntiKt_Cluster_r02");
0046   clusterjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.3),"AntiKt_Cluster_r03");
0047   clusterjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.4),"AntiKt_Cluster_r04");
0048   clusterjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.5),"AntiKt_Cluster_r05");
0049   clusterjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.6),"AntiKt_Cluster_r06");
0050   clusterjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.7),"AntiKt_Cluster_r07");
0051   clusterjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.8),"AntiKt_Cluster_r08");
0052   clusterjetreco->set_algo_node("ANTIKT");
0053   clusterjetreco->set_input_node("CLUSTER");
0054   clusterjetreco->Verbosity(verbosity);
0055   se->registerSubsystem(clusterjetreco);
0056
0057   // track jets
0058   JetReco *trackjetreco = new JetReco();
0059   trackjetreco->add_input(new TrackJetInput(Jet::TRACK));
0060   trackjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.2),"AntiKt_Track_r02");
0061   trackjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.3),"AntiKt_Track_r03");
0062   trackjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.4),"AntiKt_Track_r04");
0063   trackjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.5),"AntiKt_Track_r05");
0064   trackjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.6),"AntiKt_Track_r06");
0065   trackjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.7),"AntiKt_Track_r07");
0066   trackjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.8),"AntiKt_Track_r08");
0067   trackjetreco->set_algo_node("ANTIKT");
0068   trackjetreco->set_input_node("TRACK");
0069   trackjetreco->Verbosity(verbosity);
0070   se->registerSubsystem(trackjetreco);
0071
0072   return;
0073 }
0074
0075 void Jet_Eval(std::string outfilename = "g4jets_eval.root",
0076           int verbosity = 0) {
0077
0078   gSystem->Load("libg4eval.so");
0079
0080   Fun4AllServer *se = Fun4AllServer::instance();
0081
0082   JetEvaluator* eval = new JetEvaluator("JETEVALUATOR",
0083                     "AntiKt_Tower_r03",
0084                     "AntiKt_Truth_r03",
0085                     outfilename);
0086   eval->Verbosity(verbosity);
0087   se->registerSubsystem(eval);
0088
0089   return;
0090 }