BPhysicsValidation.h

Go to the documentation of this file.

#ifndef BPhysicsValidation_H
#define BPhysicsValidation_H

/*class BPhysicsValidation
 *  
 *  Class to fill Event Generator dqm monitor elements; works on HepMCProduct
 *
 *
 */
#include <iostream>
#include "TMath.h"
// 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"

//DQM services
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "DQMServices/Core/interface/MonitorElement.h"
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"

#include "DataFormats/HepMCCandidate/interface/GenParticle.h"
#include "DataFormats/HepMCCandidate/interface/GenParticle.h"
#include "DataFormats/HepMCCandidate/interface/GenParticleFwd.h"

#include "Validation/EventGenerator/interface/DQMHelper.h"



class BPhysicsValidation : public DQMEDAnalyzer {
 public:
  explicit BPhysicsValidation(const edm::ParameterSet&);
  virtual ~BPhysicsValidation();

  virtual void bookHistograms(DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &) override;
  virtual void dqmBeginRun(const edm::Run& r, const edm::EventSetup& c) override;
  virtual void analyze(edm::Event const&, edm::EventSetup const&) override;
  
 private:
  
  class ParticleMonitor{
  public:
  ParticleMonitor(std::string name_,const edm::ParameterSet &p_):p(p_),name(name_),pdgid(p.getParameter<int>("pdgid")){};
    ~ParticleMonitor(){};
    
    void Configure(DQMStore::IBooker &i){
      std::string pname=p.getParameter<std::string>("pname");
      double mass_min=p.getParameter<double>("massmin");
      double mass_max=p.getParameter<double>("massmax");
      DQMHelper dqm(&i); i.setCurrentFolder("Generator/BPhysics");
      // Number of analyzed events
      pt   = dqm.book1dHisto(name+"PT", "P_{t} of the "+pname+"s", 100, 0., 100,"P_{t} (GeV)","Number of Events");
      eta  = dqm.book1dHisto(name+"ETA", "#eta of the "+pname+"s", 100, -5., 5.,"#eta","Number of Events");
      phi  = dqm.book1dHisto(name+"PHI", "#phi of the "+pname+"s", 100, 0, 2*TMath::Pi(),"#phi","Number of Events");
      mass = dqm.book1dHisto(name+"MASS", "Mass of the "+pname+"s", 100, mass_min, mass_max,"Mass (GeV)","Number of Events");
    }

    void Fill(const reco::GenParticle* p, double weight){
      if(abs(p->pdgId())==abs(pdgid)){
    pt->Fill(p->pt(),weight);
    eta->Fill(p->eta(),weight);
    phi->Fill(p->phi(),weight);
    mass->Fill(p->mass(),weight);
      }
    }
    int PDGID(){return pdgid;}
    
  private:
    const edm::ParameterSet p;
    std::string name;
    int pdgid;
    MonitorElement *pt, *eta, *phi, *mass; 
  };

  void FillDaughters(const reco::GenParticle* p);
  edm::InputTag genparticleCollection_;
  edm::EDGetTokenT<reco::GenParticleCollection> genparticleCollectionToken_;
  std::string name;
  ParticleMonitor particle;
  std::vector<ParticleMonitor> daughters;
  MonitorElement *Nobj;
};

#endif