HSCPValidator.h

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    HSCP
// Class:      HSCPValidator
//
//
// Original Author:  Seth Cooper,27 1-024,+41227672342,
//         Created:  Wed Apr 14 14:27:52 CEST 2010
//
//
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/ESHandle.h"

#include "Geometry/RPCGeometry/interface/RPCGeometry.h"
#include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"
#include "FWCore/Common/interface/TriggerResultsByName.h"
#include "DataFormats/HLTReco/interface/TriggerEvent.h"
#include "DataFormats/Common/interface/TriggerResults.h"
#include "DataFormats/HLTReco/interface/TriggerObject.h"
#include "DataFormats/TrackReco/interface/TrackFwd.h"
#include "DataFormats/Common/interface/ValueMap.h"
#include "DataFormats/TrackReco/interface/DeDxData.h"
#include "DataFormats/EcalDigi/interface/EcalDigiCollections.h"
#include "DataFormats/RPCRecHit/interface/RPCRecHitCollection.h"
#include "SimDataFormats/Track/interface/SimTrackContainer.h"
#include "SimDataFormats/CaloHit/interface/PCaloHitContainer.h"

//

#include "TH2F.h"
//
// class declaration
//

class HSCPValidator : public edm::EDAnalyzer {
public:
  explicit HSCPValidator(const edm::ParameterSet&);
  ~HSCPValidator() override;

private:
  void beginJob() override;
  void analyze(const edm::Event&, const edm::EventSetup&) override;
  void endJob() override;
  std::string intToString(int num);
  void makeGenPlots(const edm::Event& iEvent);
  void makeSimTrackPlots(const edm::Event& iEvent);
  void makeSimDigiPlotsECAL(const edm::Event& iEvent);
  void makeSimDigiPlotsRPC(const edm::Event& iEvent);
  void makeHLTPlots(const edm::Event& iEvent);
  void makeRecoPlots(const edm::Event& iEvent);
  bool IncreasedTreshold(const trigger::TriggerEvent& trEv,
                         const edm::InputTag& InputPath,
                         double NewThreshold,
                         double etaCut,
                         int NObjectAboveThreshold,
                         bool averageThreshold);
  // ----------member data ---------------------------
  bool doGenPlots_;
  bool doHLTPlots_;
  bool doSimTrackPlots_;
  bool doSimDigiPlots_;
  bool doRecoPlots_;

  // GEN section
  edm::EDGetTokenT<edm::HepMCProduct> token_;
  edm::EDGetTokenT<edm::SimTrackContainer> simTracksToken_;
  edm::EDGetTokenT<trigger::TriggerEvent> trEvToken_;
  edm::EDGetTokenT<reco::TrackCollection> tkTracksToken_;
  edm::EDGetTokenT<edm::ValueMap<reco::DeDxData> > dEdxTrackToken_;
  edm::EDGetTokenT<RPCRecHitCollection> rpcRecHitsToken_;
  edm::EDGetTokenT<edm::TriggerResults> triggerResultsToken_;
  std::vector<int> particleIds_;
  int particleStatus_;
  std::map<int, int> particleIdsFoundMap_;
  TH1F* particleEtaHist_;
  TH1F* particlePhiHist_;
  TH1F* particlePHist_;
  TH1F* particlePtHist_;
  TH1F* particleMassHist_;
  TH1F* particleStatusHist_;
  TH1F* particleBetaHist_;
  TH1F* particleBetaInverseHist_;
  TH1F* h_genhscp_met;
  TH1F* h_genhscp_met_nohscp;
  TH1F* h_genhscp_scaloret;
  TH1F* h_genhscp_scaloret_nohscp;

  //SIM-Track section
  TH1F* simTrackParticleEtaHist_;
  TH1F* simTrackParticlePhiHist_;
  TH1F* simTrackParticlePHist_;
  TH1F* simTrackParticlePtHist_;
  TH1F* simTrackParticleBetaHist_;

  // SIM-DIGI section
  edm::EDGetTokenT<edm::PCaloHitContainer> ebSimHitToken_;
  edm::EDGetTokenT<edm::PCaloHitContainer> eeSimHitToken_;
  edm::EDGetTokenT<edm::SimTrackContainer> simTrackToken_;
  edm::EDGetTokenT<EBDigiCollection> EBDigiCollectionToken_;
  edm::EDGetTokenT<EEDigiCollection> EEDigiCollectionToken_;
  edm::ESHandle<RPCGeometry> rpcGeo;
  // ECAL
  TH1F* simHitsEcalEnergyHistEB_;
  TH1F* simHitsEcalTimeHistEB_;
  TH1F* simHitsEcalNumHistEB_;
  TH2F* simHitsEcalEnergyVsTimeHistEB_;
  TH1F* simHitsEcalEnergyHistEE_;
  TH1F* simHitsEcalTimeHistEE_;
  TH1F* simHitsEcalNumHistEE_;
  TH2F* simHitsEcalEnergyVsTimeHistEE_;
  TH1F* simHitsEcalDigiMatchEnergyHistEB_;
  TH1F* simHitsEcalDigiMatchTimeHistEB_;
  TH2F* simHitsEcalDigiMatchEnergyVsTimeHistEB_;
  TH1F* simHitsEcalDigiMatchEnergyHistEE_;
  TH1F* simHitsEcalDigiMatchTimeHistEE_;
  TH2F* simHitsEcalDigiMatchEnergyVsTimeHistEE_;
  TH1F* simHitsEcalDigiMatchIEtaHist_;
  TH1F* simHitsEcalDigiMatchIPhiHist_;
  TH1F* digisEcalNumHistEB_;
  TH1F* digisEcalNumHistEE_;
  TH2F* digiOccupancyMapEB_;
  TH2F* digiOccupancyMapEEP_;
  TH2F* digiOccupancyMapEEM_;
  // RPC
  TH1F* residualsRPCRecHitSimDigis_;
  TH1F* efficiencyRPCRecHitSimDigis_;
  TH1F* cluSizeDistribution_;
  TH1F* rpcTimeOfFlightBarrel_[6];
  TH1F* rpcBXBarrel_[6];
  TH1F* rpcTimeOfFlightEndCap_[3];
  TH1F* rpcBXEndCap_[3];
  //HLT
  TH1F* hltmet;
  TH1F* hltjet;
  TH1F* hltmu;
  //RECO
  TH2F* RecoHSCPPtVsGenPt;
  TH2F* dedxVsp;
};