#include <CMSDAS11DijetTestAnalyzer.h>

Inheritance diagram for CMSDAS11DijetTestAnalyzer:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob ()
	CMSDAS11DijetTestAnalyzer (const edm::ParameterSet &)
virtual void	endJob (void)
virtual	~CMSDAS11DijetTestAnalyzer ()
Static Public Member Functions
static bool	compare_JetPt (const reco::CaloJet &jet1, const reco::CaloJet &jet2)
Private Attributes
TH1D *	hCorDijetMass
TH1D *	hCorDijetXsec
TH1D *	hDijetDeltaEta
TH1D *	hDijetDeltaPhi
TH2D *	hDijetDeltaPhiNJets
TH2D *	hDijetEta1Eta2
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

Detailed Description

Definition at line 20 of file CMSDAS11DijetTestAnalyzer.h.

Constructor & Destructor Documentation

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

Definition at line 22 of file CMSDAS11DijetTestAnalyzer.cc.

References hCorDijetMass, hDijetDeltaEta, hDijetDeltaPhi, hDijetDeltaPhiNJets, hDijetEta1Eta2, hInnerDijetMass, hJet1EMF, hJet1Eta, hJet1Phi, hJet1Pt, hJet2EMF, hJet2Eta, hJet2Phi, hJet2Pt, hJetCorrPt, hJetEMF, hJetEta, hJetPhi, hJetRawPt, hOuterDijetMass, hVertexZ, and NBINS.

                                                                                  :
  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",   50,-5,5);
  hJet1Eta      = fs->make<TH1D>("hJet1Eta","Corrected Jet1 Eta",50,-5,5);
  hJet2Eta      = fs->make<TH1D>("hJet2Eta","Corrected Jet2 Eta",50,-5,5);

  hJetPhi      = fs->make<TH1D>("hJetPhi","Corrected Jet Phi",   50,-3.1415,3.1415);
  hJet1Phi      = fs->make<TH1D>("hJet1Phi","Corrected Jet1 Phi",50,-3.1415,3.1415);
  hJet2Phi      = fs->make<TH1D>("hJet2Phi","Corrected Jet2 Phi",50,-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,1.3);
  hDijetDeltaPhiNJets= fs->make<TH2D>("hDijetDeltaPhiNJets","Dijet |#Delta #phi| v. the number of jets",50,0,3.1415,7,0.5,7.5);
  hDijetEta1Eta2=fs->make<TH2D>("hDijetEta1Eta2","Eta 1 versus Eta 2 of dijet events",50,-5,5,50,-5,5);

  hInnerDijetMass = fs->make<TH1D>("hInnerDijetMass","Corrected Inner Dijet Mass",NBINS-1,BOUNDARIES);
  hOuterDijetMass = fs->make<TH1D>("hOuterDijetMass","Corrected Outer Dijet Mass",NBINS-1,BOUNDARIES);
}

virtual CMSDAS11DijetTestAnalyzer::~CMSDAS11DijetTestAnalyzer ( ) [inline, virtual]

Definition at line 24 of file CMSDAS11DijetTestAnalyzer.h.

{}

Member Function Documentation

void CMSDAS11DijetTestAnalyzer::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)		`[virtual]`

Implements edm::EDAnalyzer.

Definition at line 72 of file CMSDAS11DijetTestAnalyzer.cc.

References gather_cfg::cout, edm::Event::getByLabel(), JetCorrector::getJetCorrector(), hVertexZ, edm::HandleBase::isValid(), metsig::jet, jetCorrections, jetSrc, reco::Vertex::ndof(), reco::Vertex::position(), JetCorrector::vectorialCorrection(), vertexSrc, and reco::Vertex::z().

{
  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->size()<=0) 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);
  corrector->vectorialCorrection();
  // 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;

    // put the selected jets into a collection
    selectedJets.push_back(jet);
  }

  hVertexZ->Fill(theVertex->z(), mWeight);
  return;
}

virtual void CMSDAS11DijetTestAnalyzer::beginJob ( void ) [inline, virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 25 of file CMSDAS11DijetTestAnalyzer.h.

{}

static bool CMSDAS11DijetTestAnalyzer::compare_JetPt	(	const reco::CaloJet &	jet1,
		const reco::CaloJet &	jet2
	)		`[inline, static]`

Definition at line 28 of file CMSDAS11DijetTestAnalyzer.h.

References reco::LeafCandidate::pt().

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

void CMSDAS11DijetTestAnalyzer::endJob ( void ) [virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 69 of file CMSDAS11DijetTestAnalyzer.cc.

{
}

Member Data Documentation

TH1D* CMSDAS11DijetTestAnalyzer::hCorDijetMass [private]

Definition at line 50 of file CMSDAS11DijetTestAnalyzer.h.

Referenced by CMSDAS11DijetTestAnalyzer().

TH1D* CMSDAS11DijetTestAnalyzer::hCorDijetXsec [private]

Definition at line 51 of file CMSDAS11DijetTestAnalyzer.h.

TH1D* CMSDAS11DijetTestAnalyzer::hDijetDeltaEta [private]

Definition at line 61 of file CMSDAS11DijetTestAnalyzer.h.

Referenced by CMSDAS11DijetTestAnalyzer().

TH1D* CMSDAS11DijetTestAnalyzer::hDijetDeltaPhi [private]

Definition at line 60 of file CMSDAS11DijetTestAnalyzer.h.

Referenced by CMSDAS11DijetTestAnalyzer().

TH2D* CMSDAS11DijetTestAnalyzer::hDijetDeltaPhiNJets [private]

Definition at line 62 of file CMSDAS11DijetTestAnalyzer.h.

Referenced by CMSDAS11DijetTestAnalyzer().

TH2D* CMSDAS11DijetTestAnalyzer::hDijetEta1Eta2 [private]

Definition at line 63 of file CMSDAS11DijetTestAnalyzer.h.

Referenced by CMSDAS11DijetTestAnalyzer().

TH1D* CMSDAS11DijetTestAnalyzer::hInnerDijetMass [private]

Definition at line 65 of file CMSDAS11DijetTestAnalyzer.h.

Referenced by CMSDAS11DijetTestAnalyzer().

TH1D* CMSDAS11DijetTestAnalyzer::hJet1EMF [private]

Definition at line 55 of file CMSDAS11DijetTestAnalyzer.h.

Referenced by CMSDAS11DijetTestAnalyzer().

TH1D* CMSDAS11DijetTestAnalyzer::hJet1Eta [private]

Definition at line 53 of file CMSDAS11DijetTestAnalyzer.h.

Referenced by CMSDAS11DijetTestAnalyzer().

TH1D* CMSDAS11DijetTestAnalyzer::hJet1Phi [private]

Definition at line 54 of file CMSDAS11DijetTestAnalyzer.h.

Referenced by CMSDAS11DijetTestAnalyzer().

TH1D* CMSDAS11DijetTestAnalyzer::hJet1Pt [private]

Definition at line 52 of file CMSDAS11DijetTestAnalyzer.h.

Referenced by CMSDAS11DijetTestAnalyzer().

TH1D* CMSDAS11DijetTestAnalyzer::hJet2EMF [private]

Definition at line 59 of file CMSDAS11DijetTestAnalyzer.h.

Referenced by CMSDAS11DijetTestAnalyzer().

TH1D* CMSDAS11DijetTestAnalyzer::hJet2Eta [private]

Definition at line 57 of file CMSDAS11DijetTestAnalyzer.h.

Referenced by CMSDAS11DijetTestAnalyzer().

TH1D* CMSDAS11DijetTestAnalyzer::hJet2Phi [private]

Definition at line 58 of file CMSDAS11DijetTestAnalyzer.h.

Referenced by CMSDAS11DijetTestAnalyzer().

TH1D* CMSDAS11DijetTestAnalyzer::hJet2Pt [private]

Definition at line 56 of file CMSDAS11DijetTestAnalyzer.h.

Referenced by CMSDAS11DijetTestAnalyzer().

TH1D* CMSDAS11DijetTestAnalyzer::hJetCorrPt [private]

Definition at line 43 of file CMSDAS11DijetTestAnalyzer.h.

Referenced by CMSDAS11DijetTestAnalyzer().

TH1D* CMSDAS11DijetTestAnalyzer::hJetEMF [private]

Definition at line 47 of file CMSDAS11DijetTestAnalyzer.h.

Referenced by CMSDAS11DijetTestAnalyzer().

TH1D* CMSDAS11DijetTestAnalyzer::hJetEta [private]

Definition at line 45 of file CMSDAS11DijetTestAnalyzer.h.

Referenced by CMSDAS11DijetTestAnalyzer().

TH1D* CMSDAS11DijetTestAnalyzer::hJetPhi [private]

Definition at line 46 of file CMSDAS11DijetTestAnalyzer.h.

Referenced by CMSDAS11DijetTestAnalyzer().

TH1D* CMSDAS11DijetTestAnalyzer::hJetRawPt [private]

Definition at line 44 of file CMSDAS11DijetTestAnalyzer.h.

Referenced by CMSDAS11DijetTestAnalyzer().

TH1D* CMSDAS11DijetTestAnalyzer::hOuterDijetMass [private]

Definition at line 66 of file CMSDAS11DijetTestAnalyzer.h.

Referenced by CMSDAS11DijetTestAnalyzer().

TH1D* CMSDAS11DijetTestAnalyzer::hRawDijetMass [private]

Definition at line 49 of file CMSDAS11DijetTestAnalyzer.h.

TH1D* CMSDAS11DijetTestAnalyzer::hVertexZ [private]

Definition at line 42 of file CMSDAS11DijetTestAnalyzer.h.

Referenced by analyze(), and CMSDAS11DijetTestAnalyzer().

double CMSDAS11DijetTestAnalyzer::innerDeltaEta [private]

Definition at line 37 of file CMSDAS11DijetTestAnalyzer.h.

double CMSDAS11DijetTestAnalyzer::JESbias [private]

Definition at line 39 of file CMSDAS11DijetTestAnalyzer.h.

std::string CMSDAS11DijetTestAnalyzer::jetCorrections [private]

Definition at line 36 of file CMSDAS11DijetTestAnalyzer.h.

Referenced by analyze().

edm::InputTag CMSDAS11DijetTestAnalyzer::jetSrc [private]

Definition at line 34 of file CMSDAS11DijetTestAnalyzer.h.

Referenced by analyze().

double CMSDAS11DijetTestAnalyzer::outerDeltaEta [private]

Definition at line 38 of file CMSDAS11DijetTestAnalyzer.h.

edm::InputTag CMSDAS11DijetTestAnalyzer::vertexSrc [private]

Definition at line 35 of file CMSDAS11DijetTestAnalyzer.h.

Referenced by analyze().

CMSDAS11DijetTestAnalyzer Class Reference

Public Member Functions

Static Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation