Analysis Software
Documentation for sPHENIX simulation software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4_FwdJets.C
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4_FwdJets.C
1 
2 void Jet_FwdReco(int verbosity = 0) {
3 
4  gSystem->Load("libg4jets.so");
5 
7 
8  // truth particle level jets
9  JetReco *truthjetreco = new JetReco();
10  truthjetreco->add_input(new TruthJetInput(Jet::PARTICLE));
11  //truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.2),"AntiKt_Truth_r02");
12  truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.3),"AntiKt_Truth_r03");
13  // truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.4),"AntiKt_Truth_r04");
14  truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.5),"AntiKt_Truth_r05");
15  // truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.6),"AntiKt_Truth_r06");
16  truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.7),"AntiKt_Truth_r07");
17  // truthjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.8),"AntiKt_Truth_r08");
18  truthjetreco->set_algo_node("ANTIKT");
19  truthjetreco->set_input_node("TRUTH");
20  truthjetreco->Verbosity(verbosity);
21  se->registerSubsystem(truthjetreco);
22 
23  // tower jets
24  JetReco *towerjetreco = new JetReco();
25  towerjetreco->add_input(new TowerJetInput(Jet::FEMC_TOWER));
26  towerjetreco->add_input(new TowerJetInput(Jet::FHCAL_TOWER));
27  towerjetreco->add_input(new TowerJetInput(Jet::CEMC_TOWER));
28  towerjetreco->add_input(new TowerJetInput(Jet::HCALIN_TOWER));
29  towerjetreco->add_input(new TowerJetInput(Jet::HCALOUT_TOWER));
30  // towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.2),"AntiKt_Tower_r02");
31  towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.3),"AntiKt_Tower_r03");
32  // towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.4),"AntiKt_Tower_r04");
33  towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.5),"AntiKt_Tower_r05");
34  // towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.6),"AntiKt_Tower_r06");
35  towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.7),"AntiKt_Tower_r07");
36  // towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.8),"AntiKt_Tower_r08");
37  towerjetreco->set_algo_node("ANTIKT");
38  towerjetreco->set_input_node("TOWER");
39  towerjetreco->Verbosity(verbosity);
40  se->registerSubsystem(towerjetreco);
41 
42  return;
43 }
44 
45 void Jet_FwdEval(std::string outfilename = "g4fwdjets_eval.root",
46  int verbosity = 0) {
47 
48  gSystem->Load("libg4eval.so");
49 
51 
52  JetEvaluator* eval = new JetEvaluator("JETEVALUATOR",
53  "AntiKt_Tower_r05",
54  "AntiKt_Truth_r05",
55  outfilename);
56  eval->Verbosity(verbosity);
57  se->registerSubsystem(eval);
58 
59  return;
60 }