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;

};