macros/leptoquarks_fun4all/G4_HIJetReco.C

0001 void HIJetReco(int verbosity = 0) {
0002
0003   //---------------
0004   // Load libraries
0005   //---------------
0006
0007   gSystem->Load("libjetbackground.so");
0008
0009   //---------------
0010   // Fun4All server
0011   //---------------
0012
0013   Fun4AllServer *se = Fun4AllServer::instance();
0014
0015   RetowerCEMC *rcemc = new RetowerCEMC();
0016   rcemc->Verbosity( verbosity );
0017   se->registerSubsystem( rcemc );
0018
0019   JetReco *towerjetreco = new JetReco();
0020   towerjetreco->add_input(new TowerJetInput(Jet::CEMC_TOWER_RETOWER));
0021   towerjetreco->add_input(new TowerJetInput(Jet::HCALIN_TOWER));
0022   towerjetreco->add_input(new TowerJetInput(Jet::HCALOUT_TOWER));
0023   towerjetreco->add_algo(new FastJetAlgo(Jet::ANTIKT,0.2),"AntiKt_Tower_HIRecoSeedsRaw_r02");
0024   towerjetreco->set_algo_node("ANTIKT");
0025   towerjetreco->set_input_node("TOWER");
0026   towerjetreco->Verbosity( verbosity );
0027   se->registerSubsystem(towerjetreco);
0028
0029   DetermineTowerBackground *dtb = new DetermineTowerBackground();
0030   dtb->SetBackgroundOutputName("TowerBackground_Sub1");
0031   dtb->SetSeedType( 0 );
0032   dtb->Verbosity( verbosity );
0033   se->registerSubsystem( dtb );
0034
0035   CopyAndSubtractJets *casj = new CopyAndSubtractJets();
0036   casj->Verbosity( verbosity );
0037   se->registerSubsystem( casj );
0038
0039   DetermineTowerBackground *dtb2 = new DetermineTowerBackground();
0040   dtb2->SetBackgroundOutputName("TowerBackground_Sub2");
0041   dtb2->SetSeedType( 1 );
0042   dtb2->Verbosity( verbosity );
0043   se->registerSubsystem( dtb2 );
0044
0045   SubtractTowers *st = new SubtractTowers();
0046   st->Verbosity( verbosity );
0047   se->registerSubsystem( st );
0048
0049   JetReco *towerjetreco = new JetReco();
0050   towerjetreco->add_input(new TowerJetInput(Jet::CEMC_TOWER_SUB1));
0051   towerjetreco->add_input(new TowerJetInput(Jet::HCALIN_TOWER_SUB1));
0052   towerjetreco->add_input(new TowerJetInput(Jet::HCALOUT_TOWER_SUB1));
0053   towerjetreco->add_algo(new FastJetAlgoSub(Jet::ANTIKT,0.2,verbosity),"AntiKt_Tower_r02_Sub1");
0054   towerjetreco->add_algo(new FastJetAlgoSub(Jet::ANTIKT,0.3,verbosity),"AntiKt_Tower_r03_Sub1");
0055   towerjetreco->add_algo(new FastJetAlgoSub(Jet::ANTIKT,0.4,verbosity),"AntiKt_Tower_r04_Sub1");
0056   towerjetreco->set_algo_node("ANTIKT");
0057   towerjetreco->set_input_node("TOWER");
0058   towerjetreco->Verbosity( verbosity );
0059   se->registerSubsystem(towerjetreco);
0060
0061   return;
0062
0063 }