PFJetTester.h

Go to the documentation of this file.

#ifndef ValidationRecoJetsPFJetTester_h
#define ValidationRecoJetsPFJetTester_h

// Producer for validation histograms for PFJet objects
// F. Ratnikov, Sept. 7, 2006
// Modified by Chiyoung.Jeong Feb 2, 2010
// $Id: PFJetTester.h,v 1.5 2009/12/18 20:45:14 wmtan Exp $

#include <string>

#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "DQMServices/Core/interface/MonitorElement.h"

namespace reco {
  class PFJet;
  class GenJet;
}

class MonitorElement;

class PFJetTester : public edm::EDAnalyzer {
public:

  PFJetTester (const edm::ParameterSet&);
  ~PFJetTester();

  virtual void analyze(const edm::Event&, const edm::EventSetup&);
  virtual void beginJob() ;
  virtual void endJob() ;
 
private:
  
  void fillMatchHists (const reco::GenJet& fGenJet, const reco::PFJet& fPFJet);

  edm::InputTag mInputCollection;
  edm::InputTag mInputGenCollection;
  std::string mOutputFile;
  edm::InputTag inputMETLabel_;
  std::string METType_;
  std::string inputGenMETLabel_;
  std::string inputCaloMETLabel_;

  // count number of events
  MonitorElement* numberofevents;

  // Generic Jet Parameters
  MonitorElement* mEta;
  MonitorElement* mEtaFineBin;
  MonitorElement* mEtaFineBin1p;
  MonitorElement* mEtaFineBin2p;
  MonitorElement* mEtaFineBin3p;
  MonitorElement* mEtaFineBin1m;
  MonitorElement* mEtaFineBin2m;
  MonitorElement* mEtaFineBin3m;
  MonitorElement* mPhi;
  MonitorElement* mPhiFineBin;
  MonitorElement* mE;
  MonitorElement* mE_80;
  MonitorElement* mE_3000;
  MonitorElement* mP;
  MonitorElement* mP_80;
  MonitorElement* mP_3000;
  MonitorElement* mPt;
  MonitorElement* mPt_80;
  MonitorElement* mPt_3000;
  MonitorElement* mMass;
  MonitorElement* mMass_80;
  MonitorElement* mMass_3000;
  MonitorElement* mConstituents;
  MonitorElement* mConstituents_80;
  MonitorElement* mConstituents_3000;
  MonitorElement* mHadTiming;
  MonitorElement* mEmTiming;

  // Leading Jet Parameters
  MonitorElement* mEtaFirst;
  MonitorElement* mPhiFirst;
  MonitorElement* mEFirst;
  MonitorElement* mEFirst_80;
  MonitorElement* mEFirst_3000;
  MonitorElement* mPtFirst;
  MonitorElement* mPtFirst_80;
  MonitorElement* mPtFirst_3000;

  MonitorElement* mNJetsEtaC;
  MonitorElement* mNJetsEtaF;

  MonitorElement* mNJets1;
  MonitorElement* mNJets2;

  // DiJet Parameters
  MonitorElement* mMjj;
  MonitorElement* mMjj_3000;

  // PFJet specific
  MonitorElement* mChargedEmEnergy_80;
  MonitorElement* mChargedHadronEnergy_80;
  MonitorElement* mNeutralEmEnergy_80;
  MonitorElement* mNeutralHadronEnergy_80;

  MonitorElement* mChargedEmEnergy_3000;
  MonitorElement* mChargedHadronEnergy_3000;
  MonitorElement* mNeutralEmEnergy_3000;
  MonitorElement* mNeutralHadronEnergy_3000;

  MonitorElement* mChargedEmEnergyFraction;
  MonitorElement* mChargedHadronEnergyFraction;
  MonitorElement* mNeutralEmEnergyFraction;
  MonitorElement* mNeutralHadronEnergyFraction;

  //  MonitorElement* mMaxEInEmTowers;
  //  MonitorElement* mMaxEInHadTowers;
  //  MonitorElement* mHadEnergyInHO;
  //  MonitorElement* mHadEnergyInHB;
  //  MonitorElement* mHadEnergyInHF;
  //  MonitorElement* mHadEnergyInHE;
  //  MonitorElement* mHadEnergyInHO_80;
  //  MonitorElement* mHadEnergyInHB_80;
  //  MonitorElement* mHadEnergyInHE_80;
  //  MonitorElement* mHadEnergyInHO_3000;
  //  MonitorElement* mHadEnergyInHB_3000;
  //  MonitorElement* mHadEnergyInHE_3000;
  //  MonitorElement* mEmEnergyInEB;
  //  MonitorElement* mEmEnergyInEE;
  //  MonitorElement* mEmEnergyInHF;
  //  MonitorElement* mEmEnergyInEB_80;
  //  MonitorElement* mEmEnergyInEE_80;
  //  MonitorElement* mEmEnergyInEB_3000;
  //  MonitorElement* mEmEnergyInEE_3000;
  //  MonitorElement* mEnergyFractionHadronic;
  //  MonitorElement* mEnergyFractionEm;
  //  MonitorElement* mHFTotal;
  //  MonitorElement* mHFTotal_80;
  //  MonitorElement* mHFTotal_3000;
  //  MonitorElement* mHFLong;
  //  MonitorElement* mHFLong_80;
  //  MonitorElement* mHFLong_3000;
  //  MonitorElement* mHFShort;
  //  MonitorElement* mHFShort_80;
  //  MonitorElement* mHFShort_3000;
  //  MonitorElement* mN90;

  // pthat
  MonitorElement* mPthat_80;
  MonitorElement* mPthat_3000;

  // GenJet Generic Jet Parameters
  MonitorElement* mGenEta;
  MonitorElement* mGenPhi;
  MonitorElement* mGenPt;
  MonitorElement* mGenPt_80;
  MonitorElement* mGenPt_3000;

  // GenJet Leading Jet Parameters
  MonitorElement* mGenEtaFirst;
  MonitorElement* mGenPhiFirst;

  // PFJet<->GenJet matching
  MonitorElement* mAllGenJetsPt;
  MonitorElement* mMatchedGenJetsPt;
  MonitorElement* mAllGenJetsEta;
  MonitorElement* mMatchedGenJetsEta;
  MonitorElement* mGenJetMatchEnergyFraction;
  MonitorElement* mReverseMatchEnergyFraction;
  MonitorElement* mRMatch;
  MonitorElement* mDeltaEta;
  MonitorElement* mDeltaPhi;
  MonitorElement* mEScale;
  MonitorElement* mlinEScale;  //new
  MonitorElement* mDeltaE;

  MonitorElement* mEScale_pt10;  //new
  MonitorElement* mEScaleFineBin;  //new

  MonitorElement* mpTScaleB_s;
  MonitorElement* mpTScaleE_s;
  MonitorElement* mpTScaleF_s;
  MonitorElement* mpTScaleB_d;
  MonitorElement* mpTScaleE_d;
  MonitorElement* mpTScaleF_d;

  MonitorElement* mpTScale_60_120_s;
  MonitorElement* mpTScale_200_300_s;
  MonitorElement* mpTScale_600_900_s;
  MonitorElement* mpTScale_2700_3500_s;

  MonitorElement* mpTScale_60_120_d;
  MonitorElement* mpTScale_200_300_d;
  MonitorElement* mpTScale_600_900_d;
  MonitorElement* mpTScale_2700_3500_d;

  MonitorElement* mpTScale1DB_60_120;
  MonitorElement* mpTScale1DE_60_120;
  MonitorElement* mpTScale1DF_60_120;
  MonitorElement* mpTScale1DB_200_300;
  MonitorElement* mpTScale1DE_200_300;
  MonitorElement* mpTScale1DF_200_300;
  MonitorElement* mpTScale1DB_600_900;
  MonitorElement* mpTScale1DE_600_900;
  MonitorElement* mpTScale1DF_600_900;
  MonitorElement* mpTScale1DB_2700_3500;
  MonitorElement* mpTScale1DE_2700_3500;
  MonitorElement* mpTScale1DF_2700_3500;
  MonitorElement* mpTScale1D_60_120;
  MonitorElement* mpTScale1D_200_300;
  MonitorElement* mpTScale1D_600_900;
  MonitorElement* mpTScale1D_2700_3500;

  MonitorElement* mDelEta;
  MonitorElement* mDelPhi;
  MonitorElement* mDelPt;

  // Matching parameters
  double mMatchGenPtThreshold;
  double mGenEnergyFractionThreshold;
  double mReverseEnergyFractionThreshold;
  double mRThreshold;

  // Switch on/off unimportant histogram
  std::string  mTurnOnEverything;

  // Energy Profiles
  MonitorElement* mHadEnergyProfile;
  MonitorElement* mEmEnergyProfile;
  MonitorElement* mJetEnergyProfile;
  //  MonitorElement* mHadJetEnergyProfile;
  //  MonitorElement* mEMJetEnergyProfile;

  // CaloMET
  MonitorElement* mCaloMEx;
  MonitorElement* mCaloMEx_3000;
  MonitorElement* mCaloMEy;
  MonitorElement* mCaloMEy_3000;
  MonitorElement* mCaloMET;
  MonitorElement* mCaloMET_3000;
  MonitorElement* mCaloMETPhi;
  MonitorElement* mCaloSumET;
  MonitorElement* mCaloSumET_3000;
  MonitorElement* mCaloMETSig;
  MonitorElement* mCaloMETSig_3000;

  // RecHits
  MonitorElement* mHBEne;
  MonitorElement* mHBTime;
  MonitorElement* mHEEne;
  MonitorElement* mHETime;
  MonitorElement* mHOEne;
  MonitorElement* mHOTime;
  MonitorElement* mHFEne;
  MonitorElement* mHFTime;
  MonitorElement* mEBEne;
  MonitorElement* mEBTime;
  MonitorElement* mEEEne;
  MonitorElement* mEETime;


};
#endif