PhiSymmetryCalibration.h

Go to the documentation of this file.

#ifndef Calibration_EcalCalibAlgos_PhiSymmetryCalibration_h
#define Calibration_EcalCalibAlgos_PhiSymmetryCalibration_h
 
//
// Package:    Calibration/EcalCalibAlgos
// Class:      PhiSymmetryCalibration
//
//
// Description: performs phi-symmetry calibration
//
//
// Original Author:  David Futyan
//
// Responsible    :  Stefano Argiro, Valentina Sola
//
 
#include <vector>
 
#include "Calibration/EcalCalibAlgos/interface/EcalGeomPhiSymHelper.h"
 
// Framework
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ProducerBase.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/ESHandle.h"
 
#include "CondFormats/EcalObjects/interface/EcalIntercalibConstants.h"
#include "DataFormats/DetId/interface/DetId.h"
 
class TH1F;
 
class PhiSymmetryCalibration : public edm::EDAnalyzer {
public:
  PhiSymmetryCalibration(const edm::ParameterSet& iConfig);
 
  ~PhiSymmetryCalibration() override;
 
  void beginJob() override;
  void endRun(edm::Run const&, const edm::EventSetup&) override;
  void endLuminosityBlock(edm::LuminosityBlock const&, edm::EventSetup const&) override;
 
  void endJob() override;
 
  void analyze(const edm::Event&, const edm::EventSetup&) override;
 
  void setUp(const edm::EventSetup& setup);
 
private:
  // private member functions
 
  void getKfactors();
 
  // private data members
 
  static const int kNMiscalBinsEB = 21;
  static const float kMiscalRangeEB;
 
  static const int kNMiscalBinsEE = 41;
  static const float kMiscalRangeEE;
 
  EcalGeomPhiSymHelper e_;
 
  // Transverse energy sum arrays
  double etsum_barl_[kBarlRings][kBarlWedges][kSides];
  double etsum_endc_[kEndcWedgesX][kEndcWedgesX][kSides];
  double etsum_endc_uncorr[kEndcWedgesX][kEndcWedgesX][kSides];
  double etsumMean_barl_[kBarlRings];
  double etsumMean_endc_[kEndcEtaRings];
 
  unsigned int nhits_barl_[kBarlRings][kBarlWedges][kSides];
  unsigned int nhits_endc_[kEndcWedgesX][kEndcWedgesX][kSides];
 
  double etsum_barl_miscal_[kNMiscalBinsEB][kBarlRings];
  double etsum_endc_miscal_[kNMiscalBinsEE][kEndcEtaRings];
 
  double esumMean_barl_[kBarlRings];
  double esumMean_endc_[kEndcEtaRings];
 
  // factors to convert from ET sum deviation to miscalibration
  double k_barl_[kBarlRings];
  double k_endc_[kEndcEtaRings];
  double miscalEB_[kNMiscalBinsEB];
  double miscalEE_[kNMiscalBinsEE];
 
  std::vector<DetId> barrelCells;
  std::vector<DetId> endcapCells;
 
  // input calibration constants
  double oldCalibs_barl[kBarlRings][kBarlWedges][kSides];
  double oldCalibs_endc[kEndcWedgesX][kEndcWedgesY][kSides];
 
  // new calibration constants
  double newCalibs_barl[kBarlRings][kBarlWedges][kSides];
  double newCalibs_endc[kEndcWedgesX][kEndcWedgesX][kSides];
 
  // calibration constants not multiplied by old ones
  float epsilon_M_barl[kBarlRings][kBarlWedges][kSides];
  float epsilon_M_endc[kEndcWedgesX][kEndcWedgesY][kSides];
 
  // calibration const not corrected for k
  float rawconst_barl[kBarlRings][kBarlWedges][kSides];
  float rawconst_endc[kEndcWedgesX][kEndcWedgesX][kSides];
 
  // steering parameters
 
  std::string ecalHitsProducer_;
  std::string barrelHits_;
  std::string endcapHits_;
 
  // energy cut in the barrel
  double eCut_barl_;
 
  // parametrized energy cut EE : e_cut = ap + eta_ring*b
  double ap_;
  double b_;
 
  int eventSet_;
  int statusThreshold_;
 
  static const int kMaxEndciPhi = 360;
 
  float phi_endc[kMaxEndciPhi][kEndcEtaRings];
 
  bool reiteration_;
  std::string oldcalibfile_;  //searched for in Calibration/EcalCalibAlgos/data
 
  EcalIntercalibConstants oldCalibs_;
 
  bool isfirstpass_;
 
  // Et and E spectra
  std::vector<TH1F*> et_spectrum_b_histos;  //kBarlEtaRings
  std::vector<TH1F*> e_spectrum_b_histos;
  std::vector<TH1F*> et_spectrum_e_histos;  //kEndcEtaRings
  std::vector<TH1F*> e_spectrum_e_histos;
 
  bool spectra;
  int nevents_;
  int eventsinrun_;
  int eventsinlb_;
};
 
#endif