BasicHepMCValidation.h

Go to the documentation of this file.

#ifndef BASICHEPMCVALIDATION_H
#define BASICHEPMCVALIDATION_H

/*class BasicHepMCValidation
 *  
 *  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"

//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 "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"

#include "SimGeneral/HepPDTRecord/interface/ParticleDataTable.h"

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

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

        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:
    WeightManager wmanager_;
    edm::InputTag hepmcCollection_;

    edm::ESHandle<HepPDT::ParticleDataTable> fPDGTable ;
    

    MonitorElement* nEvt;
  
    MonitorElement *uNumber, *dNumber, *sNumber, *cNumber, *bNumber, *tNumber;
    MonitorElement *ubarNumber, *dbarNumber, *sbarNumber, *cbarNumber, *bbarNumber, *tbarNumber;
    //
    MonitorElement *eminusNumber, *nueNumber, *muminusNumber, *numuNumber, *tauminusNumber, *nutauNumber;
    MonitorElement *eplusNumber, *nuebarNumber, *muplusNumber, *numubarNumber, *tauplusNumber, *nutaubarNumber;
    //
    MonitorElement *gluNumber, *WplusNumber,*WminusNumber, *ZNumber, *gammaNumber;
    //
    MonitorElement *piplusNumber, *piminusNumber, *pizeroNumber, *KplusNumber, *KminusNumber, *KlzeroNumber, *KszeroNumber;
    MonitorElement *pNumber, *pbarNumber, *nNumber, *nbarNumber, *l0Number, *l0barNumber;
    //
    MonitorElement *DplusNumber, *DminusNumber, *DzeroNumber, *BplusNumber, *BminusNumber, *BzeroNumber, *BszeroNumber;
    //
    MonitorElement *otherPtclNumber;
    
    MonitorElement *uMomentum, *dMomentum, *sMomentum, *cMomentum, *bMomentum, *tMomentum;
    MonitorElement *ubarMomentum, *dbarMomentum, *sbarMomentum, *cbarMomentum, *bbarMomentum, *tbarMomentum;
    //
    MonitorElement *eminusMomentum, *nueMomentum, *muminusMomentum, *numuMomentum, *tauminusMomentum, *nutauMomentum;
    MonitorElement *eplusMomentum, *nuebarMomentum, *muplusMomentum, *numubarMomentum, *tauplusMomentum, *nutaubarMomentum;
    //
    MonitorElement *gluMomentum, *WplusMomentum,*WminusMomentum, *ZMomentum, *gammaMomentum;
    //
    MonitorElement *piplusMomentum, *piminusMomentum, *pizeroMomentum, *KplusMomentum, *KminusMomentum, *KlzeroMomentum,  *KszeroMomentum;
    //
    MonitorElement *pMomentum, *pbarMomentum, *nMomentum, *nbarMomentum, *l0Momentum, *l0barMomentum;
    //
    MonitorElement *DplusMomentum, *DminusMomentum, *DzeroMomentum,  *BplusMomentum, *BminusMomentum, *BzeroMomentum, *BszeroMomentum;
    //
    MonitorElement *otherPtclMomentum;

    MonitorElement *genPtclNumber; 
    MonitorElement *genVrtxNumber;
    MonitorElement *unknownPDTNumber;
    MonitorElement *outVrtxPtclNumber;
    MonitorElement *genPtclStatus;
    //
    MonitorElement *stablePtclNumber;
    MonitorElement *stableChaNumber;
    MonitorElement *stablePtclPhi;
    MonitorElement *stablePtclEta;
    MonitorElement *stablePtclCharge;
    MonitorElement *stablePtclp;
    MonitorElement *stablePtclpT;
    MonitorElement *partonNumber;
    MonitorElement *partonpT;
    MonitorElement *outVrtxStablePtclNumber;
    //
    MonitorElement *vrtxZ;
    MonitorElement *vrtxRadius;
    //
    MonitorElement *Bjorken_x;

    MonitorElement *status1ShortLived;

    MonitorElement *DeltaEcms;
    MonitorElement *DeltaPx;
    MonitorElement *DeltaPy;
    MonitorElement *DeltaPz;

    edm::EDGetTokenT<edm::HepMCProduct> hepmcCollectionToken_;

};

#endif