Analysis Software
Documentation for sPHENIX simulation software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4_HIJetReco.C
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4_HIJetReco.C
1 void HIJetReco(int verbosity = 0) {
2 
3  //---------------
4  // Load libraries
5  //---------------
6 
7  gSystem->Load("libjetbackground.so");
8 
9  //---------------
10  // Fun4All server
11  //---------------
12 
13  Fun4AllServer *se = Fun4AllServer::instance();
14 
15  RetowerCEMC *rcemc = new RetowerCEMC();
16  rcemc->Verbosity( verbosity );
17  se->registerSubsystem( rcemc );
18 
19  JetReco *towerjetreco = new JetReco();
20  towerjetreco->add_input(new TowerJetInput(Jet::CEMC_TOWER_RETOWER));
21  towerjetreco->add_input(new TowerJetInput(Jet::HCALIN_TOWER));
22  towerjetreco->add_input(new TowerJetInput(Jet::HCALOUT_TOWER));
23  towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.2),"AntiKt_Tower_HIRecoSeedsRaw_r02");
24  towerjetreco->set_algo_node("ANTIKT");
25  towerjetreco->set_input_node("TOWER");
26  towerjetreco->Verbosity( verbosity );
27  se->registerSubsystem(towerjetreco);
28 
29  DetermineTowerBackground *dtb = new DetermineTowerBackground();
30  dtb->SetBackgroundOutputName("TowerBackground_Sub1");
31  dtb->SetSeedType( 0 );
32  dtb->Verbosity( verbosity );
33  se->registerSubsystem( dtb );
34 
35  CopyAndSubtractJets *casj = new CopyAndSubtractJets();
36  casj->Verbosity( verbosity );
37  se->registerSubsystem( casj );
38 
39  DetermineTowerBackground *dtb2 = new DetermineTowerBackground();
40  dtb2->SetBackgroundOutputName("TowerBackground_Sub2");
41  dtb2->SetSeedType( 1 );
42  dtb2->Verbosity( verbosity );
43  se->registerSubsystem( dtb2 );
44 
45  SubtractTowers *st = new SubtractTowers();
46  st->Verbosity( verbosity );
47  se->registerSubsystem( st );
48 
49  JetReco *towerjetreco = new JetReco();
50  towerjetreco->add_input(new TowerJetInput(Jet::CEMC_TOWER_SUB1));
51  towerjetreco->add_input(new TowerJetInput(Jet::HCALIN_TOWER_SUB1));
52  towerjetreco->add_input(new TowerJetInput(Jet::HCALOUT_TOWER_SUB1));
53  towerjetreco->add_algo(new FastJetAlgoSub(Jet::ANTIKT,0.2,verbosity),"AntiKt_Tower_r02_Sub1");
54  towerjetreco->add_algo(new FastJetAlgoSub(Jet::ANTIKT,0.3,verbosity),"AntiKt_Tower_r03_Sub1");
55  towerjetreco->add_algo(new FastJetAlgoSub(Jet::ANTIKT,0.4,verbosity),"AntiKt_Tower_r04_Sub1");
56  towerjetreco->set_algo_node("ANTIKT");
57  towerjetreco->set_input_node("TOWER");
58  towerjetreco->Verbosity( verbosity );
59  se->registerSubsystem(towerjetreco);
60 
61  return;
62 
63 }