#include <memory>
#include <string>
#include <cstdlib>
#include <sstream>
#include <fstream>
#include <iostream>
#include "FWCore/FWLite/interface/FWLiteEnabler.h"
#include "DataFormats/PatCandidates/interface/Jet.h"
#include "TopQuarkAnalysis/Examples/bin/NiceStyle.cc"
#include "TopQuarkAnalysis/Examples/interface/RootSystem.h"
#include "TopQuarkAnalysis/Examples/interface/RootHistograms.h"
#include "TopQuarkAnalysis/Examples/interface/RootPostScript.h"

Functions
int	main (int argc, char *argv[])

Function Documentation

◆ main()

int main	(	int	argc,
		char *	argv[]
	)

===============================================================================================================================================================================================

variant2: for each run define phi-averaged A for normalization channel (Dref,16) and then, divide Rijk on it, i.e. get RRijk

eta=27

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=19

eta=17

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=20

eta=19

eta=18

eta=27 L1=1

eta=25 L1=1

eta=23 L1=1

eta=22 L1=1

eta=21 L1=1

eta=29 L1=1

eta=26 L1=1

eta=24 L1=1

eta=20 L1=1

eta=19 L1=1

eta=18 L1=1

eta=17 L1=1

eta=28 L7=1

eta=27 L7=1

eta=25 L7=1

eta=23 L7=1

eta=22 L7=1

eta=21 L7=1

eta=26 L7=1

eta=24 L7=1

eta=20 L7=1

eta=19 L7=1

eta=18 L7=1

eta=17 L7=1

eta=27

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=19

eta=17

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=20

eta=19

eta=18

eta=27 L1=1

eta=25 L1=1

eta=23 L1=1

eta=22 L1=1

eta=21 L1=1

eta=26 L1=1

eta=24 L1=1

eta=20 L1=1

eta=19 L1=1

eta=18 L1=1

eta=17 L1=1

eta=28 L7=1

eta=27 L7=1

eta=25 L7=1

eta=23 L7=1

eta=22 L7=1

eta=21 L7=1

eta=26 L7=1

eta=24 L7=1

eta=20 L7=1

eta=19 L7=1

eta=18 L7=1

eta=17 L7=1

eta=27

eta=28

errA with average Amplitudes

Summed Amplitude Plots:

RBX:

errA with average Amplitudes

Summed Amplitude Plots:

RBX:

errA with average Amplitudes

Summed Amplitude Plots:

RBX:

Prepare maps of good/bad channels:

CALO JETS

PF JETS

Definition at line 16 of file TopJetFWLiteAnalyzer.cc.

                                  {
   if (argc < 3) {
     // -------------------------------------------------
     std::cerr << "ERROR:: "
               << "Wrong number of arguments! Please specify:" << std::endl
               << "        * filepath" << std::endl
               << "        * process name" << std::endl;
     // -------------------------------------------------
     return -1;
   }
  
   // load framework libraries
   gSystem->Load("libFWCoreFWLite");
   FWLiteEnabler::enable();
  
   // set nice style for histograms
   setNiceStyle();
  
   // define some histograms
   TH1I* noJets = new TH1I("noJets", "N_{Jets}", 10, 0, 10);
   TH1F* ptJets = new TH1F("ptJets", "pt_{Jets}", 100, 0., 300.);
   TH1F* enJets = new TH1F("enJets", "energy_{Jets}", 100, 0., 300.);
   TH1F* etaJets = new TH1F("etaJets", "eta_{Jets}", 100, -3., 3.);
   TH1F* phiJets = new TH1F("phiJets", "phi_{Jets}", 100, -5., 5.);
  
   // -------------------------------------------------
   std::cout << "open  file: " << argv[1] << std::endl;
   // -------------------------------------------------
   TFile* inFile = TFile::Open(argv[1]);
   TTree* events_ = nullptr;
   if (inFile)
     inFile->GetObject("Events", events_);
   if (events_ == nullptr) {
     // -------------------------------------------------
     std::cerr << "ERROR:: "
               << "Unable to retrieve TTree Events!" << std::endl
               << "        Eighter wrong file name or the the tree doesn't exists" << std::endl;
     // -------------------------------------------------
     return -1;
   }
  
   // acess branch of elecs
   char jetsName[50];
   sprintf(jetsName, "patJets_selectedPatJets__%s.obj", argv[2]);
   TBranch* jets_ = events_->GetBranch(jetsName);
   assert(jets_ != nullptr);
  
   // loop over events and fill histograms
   std::vector<pat::Jet> jets;
   int nevt = events_->GetEntries();
  
   // -------------------------------------------------
   std::cout << "start looping " << nevt << " events..." << std::endl;
   // -------------------------------------------------
   for (int evt = 0; evt < nevt; ++evt) {
     // set branch address
     jets_->SetAddress(&jets);
     // get event
     jets_->GetEntry(evt);
     events_->GetEntry(evt, 0);
  
     // -------------------------------------------------
     if (evt > 0 && !(evt % 10))
       std::cout << "  processing event: " << evt << std::endl;
     // -------------------------------------------------
  
     // fill histograms
     noJets->Fill(jets.size());
     for (unsigned idx = 0; idx < jets.size(); ++idx) {
       // fill histograms
       ptJets->Fill(jets[idx].pt());
       enJets->Fill(jets[idx].energy());
       etaJets->Fill(jets[idx].eta());
       phiJets->Fill(jets[idx].phi());
     }
   }
   // -------------------------------------------------
   std::cout << "close file" << std::endl;
   // -------------------------------------------------
   inFile->Close();
  
   // save histograms to file
   TFile outFile("analyzeJets.root", "recreate");
   outFile.mkdir("analyzeJet");
   outFile.cd("analyzeJet");
   noJets->Write();
   ptJets->Write();
   enJets->Write();
   etaJets->Write();
   phiJets->Write();
   outFile.Close();
  
   // free allocated space
   delete noJets;
   delete ptJets;
   delete enJets;
   delete etaJets;
   delete phiJets;
  
   return 0;
 }

References dir2webdir::argc, cmsBatch::argv, cms::cuda::assert(), beam_dqm_sourceclient-live_cfg::cerr, gather_cfg::cout, FWLiteEnabler::enable(), HCALHighEnergyHPDFilter_cfi::energy, PVValHelper::eta, b2gHadronicHLTEventValidation_cfi::etaJets, training_settings::idx, singleTopDQM_cfi::jets, nevt, L1TdeCSCTF_cfi::outFile, PVValHelper::phi, DiDispStaMuonMonitor_cfi::pt, and setNiceStyle().

Functions

Function Documentation

◆ main()

CALO JETS

PF JETS