MixCollectionValidation.h

Go to the documentation of this file.

#ifndef MixCollectionValidation_H
#define MixCollectionValidation_H

// system include files
#include "FWCore/Framework/interface/Frameworkfwd.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/Utilities/interface/InputTag.h"

#include "SimDataFormats/CrossingFrame/interface/CrossingFrame.h"
#include "SimDataFormats/CrossingFrame/interface/MixCollection.h"

#include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"
#include "SimDataFormats/TrackingHit/interface/PSimHit.h"

#include "SimDataFormats/CaloHit/interface/PCaloHitContainer.h"
#include "SimDataFormats/Track/interface/SimTrackContainer.h"

// DQM services for histogram
#include <DQMServices/Core/interface/DQMEDAnalyzer.h>
#include <DQMServices/Core/interface/DQMStore.h>

#include <string>

//
// class declaration
//

class MixCollectionValidation : public DQMEDAnalyzer {
public:
  explicit MixCollectionValidation(const edm::ParameterSet &);
  ~MixCollectionValidation() override;
  void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override;
  void analyze(const edm::Event &, const edm::EventSetup &) override;

private:
  edm::ParameterSet mixObjextsSet_;

  template <class T1, class T2>
  void fillMultiplicity(T1 &theItr_, T2 &theColl_, MonitorElement *theProfile_);

  template <class T1, class T2>
  void fillGenParticleMulti(T1 &theItr_, T2 &theColl_, MonitorElement *theProfile_);

  template <class T1, class T2>
  void fillSimHitTime(T1 &theItr_, T2 &theColl_, MonitorElement *theProfile_);

  template <class T1, class T2>
  void fillCaloHitTime(T1 &theItr_, T2 &theColl_, MonitorElement *theProfile_);

  /* N.B. I see vector<InputTag> as private members of this class, but
     in the corresponding C++ only the first element, if present, is
     used to get products from the event. Hence I did implement a
     single Token for each kind of objects, not a vector of
     Tokens. For all but PSimHitTags_ and PCaloHitTags_, which have a
     corresponding vector of Tokens. */

  edm::EDGetTokenT<CrossingFrame<edm::HepMCProduct>> crossingFrame_Hep_Token_;
  edm::EDGetTokenT<CrossingFrame<SimTrack>> crossingFrame_SimTr_Token_;
  edm::EDGetTokenT<CrossingFrame<SimVertex>> crossingFrame_SimVtx_Token_;
  std::vector<edm::EDGetTokenT<CrossingFrame<PSimHit>>> crossingFrame_PSimHit_Tokens_;
  std::vector<edm::EDGetTokenT<CrossingFrame<PCaloHit>>> crossingFrame_PCaloHit_Tokens_;

  std::string outputFile_;
  int minbunch_;
  int maxbunch_;

  bool verbose_;

  MonitorElement *nrHepMCProductH_;
  MonitorElement *nrSimTrackH_;
  MonitorElement *nrSimVertexH_;

  std::map<std::string, MonitorElement *> SimHitNrmap_;
  std::map<std::string, MonitorElement *> SimHitTimemap_;

  std::map<std::string, MonitorElement *> CaloHitNrmap_;
  std::map<std::string, MonitorElement *> CaloHitTimemap_;

  std::vector<edm::InputTag> HepMCProductTags_;
  std::vector<edm::InputTag> SimTrackTags_;
  std::vector<edm::InputTag> SimVertexTags_;
  std::vector<edm::InputTag> PSimHitTags_;
  std::vector<edm::InputTag> PCaloHitTags_;

  int nbin_;
};

#endif