EcalPreshowerSimHitsValidation.h

Go to the documentation of this file.

#ifndef EcalPreshowerSimHitsValidation_H
#define EcalPreshowerSimHitsValidation_H

/*
 * \file EcalPreshowerSimHitsValidation.h
 *
 * \author C.Rovelli
 *
*/

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/ConsumesCollector.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/ServiceRegistry/interface/Service.h"

#include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"
#include "SimDataFormats/CaloHit/interface/PCaloHit.h"
#include "SimDataFormats/CaloHit/interface/PCaloHitContainer.h"
#include "SimDataFormats/ValidationFormats/interface/PValidationFormats.h"

#include <iostream>
#include <fstream>
#include <vector>
#include <map>
#include "DQMServices/Core/interface/MonitorElement.h"


class EcalPreshowerSimHitsValidation: public edm::EDAnalyzer{

    typedef std::map<uint32_t,float,std::less<uint32_t> >  MapType;

public:

EcalPreshowerSimHitsValidation(const edm::ParameterSet& ps);

~EcalPreshowerSimHitsValidation() override;

protected:

void analyze(const edm::Event& e, const edm::EventSetup& c) override;

// BeginJob
void beginJob() override;

// EndJob
void endJob(void) override;

private:

 std::string HepMCLabel; 
 std::string g4InfoLabel; 
 std::string EEHitsCollection;
 std::string ESHitsCollection;

 edm::EDGetTokenT<edm::HepMCProduct> HepMCToken;
 edm::EDGetTokenT<edm::PCaloHitContainer> EEHitsToken;
 edm::EDGetTokenT<edm::PCaloHitContainer> ESHitsToken;
 
 bool verbose_;
 
 DQMStore* dbe_;
 
 std::string outputFile_;

 MonitorElement* menESHits1zp_;
 MonitorElement* menESHits2zp_;

 MonitorElement* menESHits1zm_;
 MonitorElement* menESHits2zm_;

 MonitorElement* meEShitLog10Energy_;
 MonitorElement* meEShitLog10EnergyNorm_;

 MonitorElement* meESEnergyHits1zp_;
 MonitorElement* meESEnergyHits2zp_;

 MonitorElement* meESEnergyHits1zm_;
 MonitorElement* meESEnergyHits2zm_;

 MonitorElement* meE1alphaE2zp_;
 MonitorElement* meE1alphaE2zm_;

 MonitorElement* meEEoverESzp_;
 MonitorElement* meEEoverESzm_;

 MonitorElement* me2eszpOver1eszp_; 
 MonitorElement* me2eszmOver1eszm_; 

};

#endif