#include <CMSDAS11DijetAnalyzer.h>

Inheritance diagram for CMSDAS11DijetAnalyzer:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	beginJob () override

	CMSDAS11DijetAnalyzer (const edm::ParameterSet &)

void	endJob (void) override

	~CMSDAS11DijetAnalyzer () override

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase &	operator= (EDConsumerBase &&)=default

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static bool	compare_JetPt (const reco::CaloJet &jet1, const reco::CaloJet &jet2)

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Private Attributes
TH1D *	hCorDijetMass

TH1D *	hCorDijetXsec

TH1D *	hDijetDeltaEta

TH1D *	hDijetDeltaPhi

TH1D *	hInnerDijetMass

TH1D *	hJet1EMF

TH1D *	hJet1Eta

TH1D *	hJet1Phi

TH1D *	hJet1Pt

TH1D *	hJet2EMF

TH1D *	hJet2Eta

TH1D *	hJet2Phi

TH1D *	hJet2Pt

TH1D *	hJetCorrPt

TH1D *	hJetEMF

TH1D *	hJetEta

TH1D *	hJetPhi

TH1D *	hJetRawPt

TH1D *	hOuterDijetMass

TH1D *	hRawDijetMass

TH1D *	hVertexZ

double	innerDeltaEta

double	JESbias

std::string	jetCorrections

edm::InputTag	jetSrc

double	outerDeltaEta

edm::InputTag	vertexSrc

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Definition at line 20 of file CMSDAS11DijetAnalyzer.h.

Constructor & Destructor Documentation

◆ CMSDAS11DijetAnalyzer()

CMSDAS11DijetAnalyzer::CMSDAS11DijetAnalyzer ( const edm::ParameterSet & params )

Definition at line 20 of file CMSDAS11DijetAnalyzer.cc.

     : edm::EDAnalyzer(),
       jetSrc(params.getParameter<edm::InputTag>("jetSrc")),
       vertexSrc(params.getParameter<edm::InputTag>("vertexSrc")),
       jetCorrections(params.getParameter<std::string>("jetCorrections")),
       innerDeltaEta(params.getParameter<double>("innerDeltaEta")),
       outerDeltaEta(params.getParameter<double>("outerDeltaEta")),
       JESbias(params.getParameter<double>("JESbias")) {
   // setup file service
   edm::Service<TFileService> fs;
  
   const int NBINS = 36;
   Double_t BOUNDARIES[NBINS] = {220,  244,  270,  296,  325,  354,  386,  419,  453,  489,  526,  565,
                                 606,  649,  693,  740,  788,  838,  890,  944,  1000, 1058, 1118, 1181,
                                 1246, 1313, 1383, 1455, 1530, 1607, 1687, 1770, 1856, 1945, 2037, 2132};
  
   // setup histograms
   hVertexZ = fs->make<TH1D>("hVertexZ", "Z position of the Vertex", 50, -20, 20);
   hJetRawPt = fs->make<TH1D>("hJetRawPt", "Raw Jet Pt", 50, 0, 1000);
   hJetCorrPt = fs->make<TH1D>("hJetCorrPt", "Corrected Jet Pt", 50, 0, 1000);
   hJet1Pt = fs->make<TH1D>("hJet1Pt", "Corrected Jet1 Pt", 50, 0, 1000);
   hJet2Pt = fs->make<TH1D>("hJet2Pt", "Corrected Jet2 Pt", 50, 0, 1000);
  
   hJetEta = fs->make<TH1D>("hJetEta", "Corrected Jet Eta", 10, -5, 5);
   hJet1Eta = fs->make<TH1D>("hJet1Eta", "Corrected Jet1 Eta", 10, -5, 5);
   hJet2Eta = fs->make<TH1D>("hJet2Eta", "Corrected Jet2 Eta", 10, -5, 5);
  
   hJetPhi = fs->make<TH1D>("hJetPhi", "Corrected Jet Phi", 10, -3.1415, 3.1415);
   hJet1Phi = fs->make<TH1D>("hJet1Phi", "Corrected Jet1 Phi", 10, -3.1415, 3.1415);
   hJet2Phi = fs->make<TH1D>("hJet2Phi", "Corrected Jet2 Phi", 10, -3.1415, 3.1415);
  
   hJetEMF = fs->make<TH1D>("hJetEMF", "EM Fraction of Jets", 50, 0, 1);
   hJet1EMF = fs->make<TH1D>("hJet1EMF", "EM Fraction of Jet1", 50, 0, 1);
   hJet2EMF = fs->make<TH1D>("hJet2EMF", "EM Fraction of Jet2", 50, 0, 1);
  
   hCorDijetMass = fs->make<TH1D>("hCorDijetMass", "Corrected Dijet Mass", NBINS - 1, BOUNDARIES);
   hDijetDeltaPhi = fs->make<TH1D>("hDijetDeltaPhi", "Dijet |#Delta #phi|", 50, 0, 3.1415);
   hDijetDeltaEta = fs->make<TH1D>("hDijetDeltaEta", "Dijet |#Delta #eta|", 50, 0, 6);
  
   hInnerDijetMass = fs->make<TH1D>("hInnerDijetMass", "Corrected Inner Dijet Mass", NBINS - 1, BOUNDARIES);
   hOuterDijetMass = fs->make<TH1D>("hOuterDijetMass", "Corrected Outer Dijet Mass", NBINS - 1, BOUNDARIES);
 }

References hCorDijetMass, hDijetDeltaEta, hDijetDeltaPhi, hInnerDijetMass, hJet1EMF, hJet1Eta, hJet1Phi, hJet1Pt, hJet2EMF, hJet2Eta, hJet2Phi, hJet2Pt, hJetCorrPt, hJetEMF, hJetEta, hJetPhi, hJetRawPt, hOuterDijetMass, hVertexZ, TFileService::make(), and NBINS.

◆ ~CMSDAS11DijetAnalyzer()

CMSDAS11DijetAnalyzer::~CMSDAS11DijetAnalyzer ( )

inlineoverride

Definition at line 24 of file CMSDAS11DijetAnalyzer.h.

24 {}

Member Function Documentation

◆ analyze()

void CMSDAS11DijetAnalyzer::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overridevirtual

Implements edm::EDAnalyzer.

Definition at line 65 of file CMSDAS11DijetAnalyzer.cc.

                                                                                        {
   double mWeight;
   edm::Handle<GenEventInfoProduct> hEventInfo;
   iEvent.getByLabel("generator", hEventInfo);
   if (hEventInfo.isValid())
     mWeight = hEventInfo->weight();
   else
     mWeight = 1.;
  
   // Get event ID information
   //int nrun=iEvent.id().run();
   //int nlumi=iEvent.luminosityBlock();
   //int nevent=iEvent.id().event();
  
   // Get Primary Vertex Information
  
   // magic to get the vertices from EDM
   edm::Handle<std::vector<reco::Vertex> > vertices_h;
   iEvent.getByLabel(vertexSrc, vertices_h);
   if (!vertices_h.isValid()) {
     std::cout << "Didja hear the one about the empty vertex collection?\n";
     return;
   }
  
   // require in the event that there is at least one reconstructed vertex
   if (vertices_h->empty())
     return;
  
   // pick the first (i.e. highest sum pt) verte
   const reco::Vertex* theVertex = &(vertices_h->front());
  
   // require that the vertex meets certain criteria
   if (theVertex->ndof() < 5)
     return;
   if (fabs(theVertex->z()) > 24.0)
     return;
   if (fabs(theVertex->position().rho()) > 2.0)
     return;
  
   // Get Jet Information
  
   // magic to get the jets from EDM
   edm::Handle<reco::CaloJetCollection> jets_h;
   iEvent.getByLabel(jetSrc, jets_h);
   if (!jets_h.isValid()) {
     std::cout << "Didja hear the one about the empty jet collection?\n";
     return;
   }
  
   // magic to get the jet energy corrections
   const JetCorrector* corrector = JetCorrector::getJetCorrector(jetCorrections, iSetup);
  
   // collection of selected jets
   std::vector<reco::CaloJet> selectedJets;
  
   // loop over the jet collection
   for (reco::CaloJetCollection::const_iterator j_it = jets_h->begin(); j_it != jets_h->end(); j_it++) {
     reco::CaloJet jet = *j_it;
  
     // calculate and apply the correction
     double scale = corrector->correction(jet.p4());
  
     // Introduce a purposeful bias to the correction, to show what happens
     scale *= JESbias;
  
     // fill the histograms
     hJetRawPt->Fill(jet.pt(), mWeight);
     hJetEta->Fill(jet.eta(), mWeight);
     hJetPhi->Fill(jet.phi(), mWeight);
     hJetEMF->Fill(jet.emEnergyFraction(), mWeight);
  
     // correct the jet energy
     jet.scaleEnergy(scale);
  
     // now fill the correct jet pt after the correction
     hJetCorrPt->Fill(jet.pt(), mWeight);
  
     // put the selected jets into a collection
     selectedJets.push_back(jet);
   }
  
   // require at least two jets to continue
   if (selectedJets.size() < 2)
     return;
  
   //sort by corrected pt (not the same order as raw pt, sometimes)
   sort(selectedJets.begin(), selectedJets.end(), compare_JetPt);
  
   // select high pt, central, non-noise-like jets
   if (selectedJets[0].pt() < 50.0)
     return;
   if (fabs(selectedJets[0].eta()) > 2.5)
     return;
   if (selectedJets[0].emEnergyFraction() < 0.01)
     return;
   if (selectedJets[1].pt() < 50.0)
     return;
   if (fabs(selectedJets[1].eta()) > 2.5)
     return;
   if (selectedJets[1].emEnergyFraction() < 0.01)
     return;
  
   // fill histograms for the jets in our dijets, only
   hJet1Pt->Fill(selectedJets[0].pt(), mWeight);
   hJet1Eta->Fill(selectedJets[0].eta(), mWeight);
   hJet1Phi->Fill(selectedJets[0].phi(), mWeight);
   hJet1EMF->Fill(selectedJets[0].emEnergyFraction(), mWeight);
   hJet2Pt->Fill(selectedJets[1].pt(), mWeight);
   hJet2Eta->Fill(selectedJets[1].eta(), mWeight);
   hJet2Phi->Fill(selectedJets[1].phi(), mWeight);
   hJet2EMF->Fill(selectedJets[1].emEnergyFraction(), mWeight);
  
   //Get the mass of the two leading jets.  Needs their 4-vectors...
   double corMass = (selectedJets[0].p4() + selectedJets[1].p4()).M();
   double deltaEta = fabs(selectedJets[0].eta() - selectedJets[1].eta());
   if (corMass < 489)
     return;
   if (deltaEta > 1.3)
     return;
   hCorDijetMass->Fill(corMass, mWeight);
   hDijetDeltaPhi->Fill(fabs(selectedJets[0].phi() - selectedJets[1].phi()), mWeight);
   hDijetDeltaEta->Fill(deltaEta, mWeight);
  
   //Fill the inner and outer dijet mass spectra, to make the ratio from
   if (deltaEta < innerDeltaEta)
     hInnerDijetMass->Fill(corMass, mWeight);
   else if (deltaEta < outerDeltaEta)
     hOuterDijetMass->Fill(corMass, mWeight);
  
   hVertexZ->Fill(theVertex->z(), mWeight);
   return;
 }

References compare_JetPt(), pfClustersFromHGC3DClusters_cfi::corrector, gather_cfg::cout, spr::deltaEta, PVValHelper::eta, JetCorrector::getJetCorrector(), hCorDijetMass, hDijetDeltaEta, hDijetDeltaPhi, hInnerDijetMass, hJet1EMF, hJet1Eta, hJet1Phi, hJet1Pt, hJet2EMF, hJet2Eta, hJet2Phi, hJet2Pt, hJetCorrPt, hJetEMF, hJetEta, hJetPhi, hJetRawPt, hOuterDijetMass, hVertexZ, iEvent, innerDeltaEta, edm::HandleBase::isValid(), JESbias, metsig::jet, jetCorrections, jetSrc, reco::Vertex::ndof(), outerDeltaEta, phi, reco::Vertex::position(), DiDispStaMuonMonitor_cfi::pt, Scenarios_cff::scale, objectSelection_cff::selectedJets, vertexSrc, GenEventInfoProduct::weight(), and reco::Vertex::z().

◆ beginJob()

void CMSDAS11DijetAnalyzer::beginJob ( void )

inlineoverridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 25 of file CMSDAS11DijetAnalyzer.h.

25 {}

◆ compare_JetPt()

static bool CMSDAS11DijetAnalyzer::compare_JetPt	(	const reco::CaloJet &	jet1,
		const reco::CaloJet &	jet2
	)

inlinestatic

Definition at line 28 of file CMSDAS11DijetAnalyzer.h.

28 { return (jet1.pt() > jet2.pt()); }

References reco::LeafCandidate::pt().

Referenced by analyze().

◆ endJob()

void CMSDAS11DijetAnalyzer::endJob ( void )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 63 of file CMSDAS11DijetAnalyzer.cc.

63 {}

Member Data Documentation

◆ hCorDijetMass

TH1D* CMSDAS11DijetAnalyzer::hCorDijetMass

private

Definition at line 48 of file CMSDAS11DijetAnalyzer.h.

Referenced by analyze(), and CMSDAS11DijetAnalyzer().

◆ hCorDijetXsec

TH1D* CMSDAS11DijetAnalyzer::hCorDijetXsec

private

Definition at line 49 of file CMSDAS11DijetAnalyzer.h.

◆ hDijetDeltaEta

TH1D* CMSDAS11DijetAnalyzer::hDijetDeltaEta

private

Definition at line 59 of file CMSDAS11DijetAnalyzer.h.

Referenced by analyze(), and CMSDAS11DijetAnalyzer().

◆ hDijetDeltaPhi

TH1D* CMSDAS11DijetAnalyzer::hDijetDeltaPhi

private

Definition at line 58 of file CMSDAS11DijetAnalyzer.h.

Referenced by analyze(), and CMSDAS11DijetAnalyzer().

◆ hInnerDijetMass

TH1D* CMSDAS11DijetAnalyzer::hInnerDijetMass

private

Definition at line 61 of file CMSDAS11DijetAnalyzer.h.

Referenced by analyze(), and CMSDAS11DijetAnalyzer().

◆ hJet1EMF

TH1D* CMSDAS11DijetAnalyzer::hJet1EMF

private

Definition at line 53 of file CMSDAS11DijetAnalyzer.h.

Referenced by analyze(), and CMSDAS11DijetAnalyzer().

◆ hJet1Eta

TH1D* CMSDAS11DijetAnalyzer::hJet1Eta

private

Definition at line 51 of file CMSDAS11DijetAnalyzer.h.

Referenced by analyze(), and CMSDAS11DijetAnalyzer().

◆ hJet1Phi

TH1D* CMSDAS11DijetAnalyzer::hJet1Phi

private

Definition at line 52 of file CMSDAS11DijetAnalyzer.h.

Referenced by analyze(), and CMSDAS11DijetAnalyzer().

◆ hJet1Pt

TH1D* CMSDAS11DijetAnalyzer::hJet1Pt

private

Definition at line 50 of file CMSDAS11DijetAnalyzer.h.

Referenced by analyze(), and CMSDAS11DijetAnalyzer().

◆ hJet2EMF

TH1D* CMSDAS11DijetAnalyzer::hJet2EMF

private

Definition at line 57 of file CMSDAS11DijetAnalyzer.h.

Referenced by analyze(), and CMSDAS11DijetAnalyzer().

◆ hJet2Eta

TH1D* CMSDAS11DijetAnalyzer::hJet2Eta

private

Definition at line 55 of file CMSDAS11DijetAnalyzer.h.

Referenced by analyze(), and CMSDAS11DijetAnalyzer().

◆ hJet2Phi

TH1D* CMSDAS11DijetAnalyzer::hJet2Phi

private

Definition at line 56 of file CMSDAS11DijetAnalyzer.h.

Referenced by analyze(), and CMSDAS11DijetAnalyzer().

◆ hJet2Pt

TH1D* CMSDAS11DijetAnalyzer::hJet2Pt

private

Definition at line 54 of file CMSDAS11DijetAnalyzer.h.

Referenced by analyze(), and CMSDAS11DijetAnalyzer().

◆ hJetCorrPt

TH1D* CMSDAS11DijetAnalyzer::hJetCorrPt

private

Definition at line 41 of file CMSDAS11DijetAnalyzer.h.

Referenced by analyze(), and CMSDAS11DijetAnalyzer().

◆ hJetEMF

TH1D* CMSDAS11DijetAnalyzer::hJetEMF

private

Definition at line 45 of file CMSDAS11DijetAnalyzer.h.

Referenced by analyze(), and CMSDAS11DijetAnalyzer().

◆ hJetEta

TH1D* CMSDAS11DijetAnalyzer::hJetEta

private

Definition at line 43 of file CMSDAS11DijetAnalyzer.h.

Referenced by analyze(), and CMSDAS11DijetAnalyzer().

◆ hJetPhi

TH1D* CMSDAS11DijetAnalyzer::hJetPhi

private

Definition at line 44 of file CMSDAS11DijetAnalyzer.h.

Referenced by analyze(), and CMSDAS11DijetAnalyzer().

◆ hJetRawPt

TH1D* CMSDAS11DijetAnalyzer::hJetRawPt

private

Definition at line 42 of file CMSDAS11DijetAnalyzer.h.

Referenced by analyze(), and CMSDAS11DijetAnalyzer().

◆ hOuterDijetMass

TH1D* CMSDAS11DijetAnalyzer::hOuterDijetMass

private

Definition at line 62 of file CMSDAS11DijetAnalyzer.h.

Referenced by analyze(), and CMSDAS11DijetAnalyzer().

◆ hRawDijetMass

TH1D* CMSDAS11DijetAnalyzer::hRawDijetMass

private

Definition at line 47 of file CMSDAS11DijetAnalyzer.h.

◆ hVertexZ

TH1D* CMSDAS11DijetAnalyzer::hVertexZ

private

Definition at line 40 of file CMSDAS11DijetAnalyzer.h.

Referenced by analyze(), and CMSDAS11DijetAnalyzer().

◆ innerDeltaEta

double CMSDAS11DijetAnalyzer::innerDeltaEta

private

Definition at line 35 of file CMSDAS11DijetAnalyzer.h.

Referenced by analyze().

◆ JESbias

double CMSDAS11DijetAnalyzer::JESbias

private

Definition at line 37 of file CMSDAS11DijetAnalyzer.h.

Referenced by analyze().

◆ jetCorrections

std::string CMSDAS11DijetAnalyzer::jetCorrections

private

Definition at line 34 of file CMSDAS11DijetAnalyzer.h.

Referenced by analyze().

◆ jetSrc

edm::InputTag CMSDAS11DijetAnalyzer::jetSrc

private

Definition at line 32 of file CMSDAS11DijetAnalyzer.h.

Referenced by analyze().

◆ outerDeltaEta

double CMSDAS11DijetAnalyzer::outerDeltaEta

private

Definition at line 36 of file CMSDAS11DijetAnalyzer.h.

Referenced by analyze().

◆ vertexSrc

edm::InputTag CMSDAS11DijetAnalyzer::vertexSrc

private

Definition at line 33 of file CMSDAS11DijetAnalyzer.h.

Referenced by analyze().

Public Member Functions

Static Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ CMSDAS11DijetAnalyzer()

◆ ~CMSDAS11DijetAnalyzer()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ compare_JetPt()

◆ endJob()

Member Data Documentation

◆ hCorDijetMass

◆ hCorDijetXsec

◆ hDijetDeltaEta

◆ hDijetDeltaPhi

◆ hInnerDijetMass

◆ hJet1EMF

◆ hJet1Eta

◆ hJet1Phi

◆ hJet1Pt

◆ hJet2EMF

◆ hJet2Eta

◆ hJet2Phi

◆ hJet2Pt

◆ hJetCorrPt

◆ hJetEMF

◆ hJetEta

◆ hJetPhi

◆ hJetRawPt

◆ hOuterDijetMass

◆ hRawDijetMass

◆ hVertexZ

◆ innerDeltaEta

◆ JESbias

◆ jetCorrections

◆ jetSrc

◆ outerDeltaEta

◆ vertexSrc