MBUEandQCDValidation.h

Go to the documentation of this file.

#ifndef MBUEandQCDVALIDATION_H
#define MBUEandQCDVALIDATION_H
 
/*class MBUEandQCDValidation
 *  
 *  Class to fill Event Generator dqm monitor elements; works on HepMCProduct
 *
 *
 */
 
// framework & common header files
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/Run.h"
 
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/InputTag.h"
 
#include "FWCore/MessageLogger/interface/MessageLogger.h"
 
//DQM services
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
 
#include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"
#include "DataFormats/JetReco/interface/GenJetCollection.h"
 
#include "SimGeneral/HepPDTRecord/interface/ParticleDataTable.h"
#include "Validation/EventGenerator/interface/CaloCellManager.h"
#include "Validation/EventGenerator/interface/WeightManager.h"
 
#include <vector>
 
class MBUEandQCDValidation : public DQMEDAnalyzer {
public:
  explicit MBUEandQCDValidation(const edm::ParameterSet&);
  ~MBUEandQCDValidation() override;
 
  void analyze(const edm::Event&, const edm::EventSetup&) override;
  void bookHistograms(DQMStore::IBooker& i, edm::Run const&, edm::EventSetup const&) override;
  void dqmBeginRun(const edm::Run& r, const edm::EventSetup& c) override;
 
private:
  WeightManager wmanager_;
  edm::InputTag hepmcCollection_;
  edm::InputTag genchjetCollection_;
  edm::InputTag genjetCollection_;
 
  unsigned int verbosity_;
 
  edm::ESHandle<HepPDT::ParticleDataTable> fPDGTable;
 
  std::vector<const HepMC::GenParticle*> hepmcGPCollection;
  std::vector<double> hepmcCharge;
 
  CaloCellManager* theCalo;
 
  unsigned int getHFbin(double eta);
 
  bool isCharged(unsigned int i);
  bool isNeutral(unsigned int i);
  bool isNeutrino(unsigned int i);
 
  std::vector<double> eneInCell;
 
  MonitorElement* nEvt;
 
  MonitorElement* nNoFwdTrig;
  MonitorElement* nSaFwdTrig;
 
  MonitorElement* nbquark;
  MonitorElement* ncandbquark;
  MonitorElement* ncnobquark;
 
  MonitorElement* nEvt1;
  MonitorElement* dNchdpt1;
  MonitorElement* dNchdeta1;
 
  //QCD-10-001 analysis
  MonitorElement* nEvt2;
  MonitorElement* leadTrackpt;
  MonitorElement* leadTracketa;
  MonitorElement* dNchdeta2;
  MonitorElement* dNchdpt2;
  MonitorElement* nCha;
  MonitorElement* dNchdSpt;
  MonitorElement* dNchdphi;
  MonitorElement* dSptdphi;
  MonitorElement* nChaDenLpt;
  MonitorElement* sptDenLpt;
 
  //Charged jets
  MonitorElement* nChj;
  MonitorElement* dNchjdeta;
  MonitorElement* dNchjdpt;
  MonitorElement* leadChjpt;
  MonitorElement* leadChjeta;
  MonitorElement* pt1pt2optotch;
 
  //Identified particles multiplicities
  MonitorElement* nPPbar;
  MonitorElement* nKpm;
  MonitorElement* nK0s;
  MonitorElement* nL0;
  MonitorElement* nNNbar;
  MonitorElement* nGamma;
  MonitorElement* nXim;
  MonitorElement* nOmega;
 
  //Identified particles momentum specturm
  MonitorElement* pPPbar;
  MonitorElement* pKpm;
  MonitorElement* pK0s;
  MonitorElement* pL0;
  MonitorElement* pNNbar;
  MonitorElement* pGamma;
  MonitorElement* pXim;
  MonitorElement* pOmega;
 
  MonitorElement* elePt;
  MonitorElement* muoPt;
 
  //Jets no neutrino
  MonitorElement* nDijet;
  MonitorElement* nj;
  MonitorElement* dNjdeta;
  MonitorElement* dNjdpt;
  MonitorElement* pt1pt2optot;
  MonitorElement* pt1pt2balance;
  MonitorElement* pt1pt2Dphi;
  MonitorElement* pt1pt2InvM;
  MonitorElement* pt3Frac;
  MonitorElement* sumJEt;
  MonitorElement* missEtosumJEt;
  MonitorElement* sumPt;
  MonitorElement* sumChPt;
 
  //Forward energy flow
  MonitorElement* nHFflow;
  MonitorElement* dEdetaHFmb;
  MonitorElement* dEdetaHFdj;
 
  MonitorElement* nHFSD;
  MonitorElement* EmpzHFm;
  MonitorElement* ntHFm;
  MonitorElement* eneHFmSel;
 
  // Jet Multiplicity Analysis
  MonitorElement* _JM25njets;
  MonitorElement* _JM25ht;
  MonitorElement* _JM25pt1;
  MonitorElement* _JM25pt2;
  MonitorElement* _JM25pt3;
  MonitorElement* _JM25pt4;
  MonitorElement* _JM80njets;
  MonitorElement* _JM80ht;
  MonitorElement* _JM80pt1;
  MonitorElement* _JM80pt2;
  MonitorElement* _JM80pt3;
  MonitorElement* _JM80pt4;
 
  //differential jet rates
  MonitorElement *djr10, *djr21, *djr32, *djr43;
 
  // SumET hiostograms
  MonitorElement* _sumEt;
  MonitorElement* _sumEt1;
  MonitorElement* _sumEt2;
  MonitorElement* _sumEt3;
  MonitorElement* _sumEt4;
  MonitorElement* _sumEt5;
 
  const static unsigned int nphiBin;
  const static unsigned int initSize;
 
  edm::EDGetTokenT<edm::HepMCProduct> hepmcCollectionToken_;
  edm::EDGetTokenT<reco::GenJetCollection> genchjetCollectionToken_;
  edm::EDGetTokenT<reco::GenJetCollection> genjetCollectionToken_;
};
 
#endif