PatMcMatchFWLiteAnalyzer.cc File Reference

#include <memory>
#include <string>
#include <vector>
#include <sstream>
#include <fstream>
#include <iostream>
#include <TH2F.h>
#include <TROOT.h>
#include <TFile.h>
#include <TSystem.h>
#include "DataFormats/FWLite/interface/Handle.h"
#include "DataFormats/PatCandidates/interface/Muon.h"
#include "DataFormats/PatCandidates/interface/Jet.h"
#include "FWCore/FWLite/interface/FWLiteEnabler.h"

Go to the source code of this file.

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

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

Summed Amplitude Plots:





Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

RBX:

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

Summed Amplitude Plots:





Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

RBX:

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

Summed Amplitude Plots:





Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

RBX:

Prepare maps of good/bad channels:

Prepare maps of good/bad channels:

Prepare maps of good/bad channels:

Prepare maps of good/bad channels:

---

CALO JETS

---

PF JETS

Definition at line 19 of file PatMcMatchFWLiteAnalyzer.cc.

References gather_cfg::cout, FWLiteEnabler::enable(), edmPickEvents::event, mps_fire::i, iEvent, testHGCalDigi_cfg::inFile, PDWG_EXODelayedJetMET_cff::jets, DiMuonV_cfg::muons, L1TdeCSCTF_cfi::outFile, and reco::LeafCandidate::pt().

{
  // ----------------------------------------------------------------------
  // First Part: 
  //
  //  * enable FWLite 
  //  * book the histograms of interest 
  //  * open the input file
  // ----------------------------------------------------------------------

  // load framework libraries
  gSystem->Load( "libFWCoreFWLite" );
  FWLiteEnabler::enable();
  
  // book a set of histograms
  TH2F* muonPt_  = new TH2F( "muonPt", "Muon Pt", 60, 0., 300., 60, 0., 300. ); // 2-D histo for muon Pt
  muonPt_->SetXTitle( "gen." );
  muonPt_->SetYTitle( "reco." );
  TH2F* jetPt_   = new TH2F( "jetPt", "Jet Pt", 100, 0., 500., 100, 0., 500. ); // 2-D histo for jet Pt
  jetPt_->SetXTitle( "gen." );
  jetPt_->SetYTitle( "reco." );
  
  // open input file (can be located on castor)
  TFile* inFile = TFile::Open( "file:edmPatMcMatch.root" ); // Adapt the input file name

  // ----------------------------------------------------------------------
  // Second Part: 
  //
  //  * loop the events in the input file 
  //  * receive the collections of interest via fwlite::Handle
  //  * fill the histograms
  //  * after the loop close the input file
  // ----------------------------------------------------------------------

  // loop the events
  unsigned int iEvent=0;
  fwlite::Event event(inFile);
  for(event.toBegin(); !event.atEnd(); ++event, ++iEvent){
    // break loop after end of file is reached 
    // or after 1000 events have been processed
    if( iEvent==1000 ) break;
    
    // simple event counter
    if(iEvent>0 && iEvent%25==0){ // Reduce print-out
      std::cout << "  processing event: " << iEvent << std::endl;
    }

    // fwlite::Handle to the muon collection
    fwlite::Handle<std::vector<pat::Muon> > muons; // Access the muon collection
    muons.getByLabel(event, "cleanLayer1Muons");
    // fwlite::Handle to the muon collection
    fwlite::Handle<std::vector<pat::Jet> > jets; // Access the jet collection
    jets.getByLabel(event, "cleanLayer1Jets");
    
    // loop muon collection and fill histograms
    for(unsigned i=0; i<muons->size(); ++i){
      const reco::GenParticle * genMuon = (*muons)[i].genParticle(); // Get the matched generator muon
      if ( genMuon ) {                                               // Check for valid reference
        muonPt_->Fill( genMuon->pt(), (*muons)[i].pt() );            // Fill 2-D histo
      }
    }
    // loop jet collection and fill histograms
    for(unsigned i=0; i<jets->size(); ++i){
      const reco::GenJet * genJet = (*jets)[i].genJet(); // Get the matched generator jet
      if ( genJet ) {                                    // Check for valid reference
        jetPt_->Fill( genJet->pt(), (*jets)[i].pt() );   // Fill 2-D histo
      }
    }
  }  
  // close input file
  inFile->Close();

  // ----------------------------------------------------------------------
  // Third Part: 
  //
  //  * open the output file 
  //  * write the histograms to the output file
  //  * close the output file
  // ----------------------------------------------------------------------
  
  //open output file
  TFile outFile( "rootPatMcMatch.root", "recreate" ); // Adapt the output file name
  outFile.mkdir("analyzeMcMatchPat");                 // Adapt output file according to modifications
  outFile.cd("analyzeMcMatchPat");
  muonPt_->Write( );
  jetPt_->Write( );
  outFile.Close();
  
  // ----------------------------------------------------------------------
  // Fourth Part: 
  //
  //  * never forgett to free the memory of the histograms
  // ----------------------------------------------------------------------

  // free allocated space
  delete muonPt_; // Delete the muon histo
  delete jetPt_;  // Delete the jet histo
  
  // that's it!
  return 0;
}