PhiSymmetryCalibration_step2_SM.h

Go to the documentation of this file.

#include "Calibration/EcalCalibAlgos/interface/EcalGeomPhiSymHelper.h"
#include "CondFormats/EcalObjects/interface/EcalIntercalibConstants.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ProducerBase.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/ESHandle.h"
 
class TH1F;
class TH2F;
 
class PhiSymmetryCalibration_step2_SM : public edm::EDAnalyzer {
public:
  PhiSymmetryCalibration_step2_SM(const edm::ParameterSet& iConfig);
  ~PhiSymmetryCalibration_step2_SM() override;
 
  void beginJob() override;
  void endJob() override;
 
  void analyze(const edm::Event&, const edm::EventSetup&) override;
 
  void fillHistos();
  void fillConstantsHistos();
  void setupResidHistos();
  void outResidHistos();
 
  void setUp(const edm::EventSetup& setup);
 
  void readEtSums();
 
private:
  // Transverse energy sum arrays
  double etsum_barl_[kBarlRings][kBarlWedges][kSides];
 
  // energy sum for SM calibration
  // number of crystals in a "supercrystal"
  static const int nscx = 20;
  double etsum_barl_SM_[kBarlRings][int(kBarlWedges / nscx)][kSides];
  int nBads_barl_SM_[kBarlRings][int(kBarlWedges / nscx)][kSides];
  double etsumMean_barl_SM_[kBarlRings];
  double epsilon_M_barl_SM_[kBarlRings][int(kBarlWedges / nscx)][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 esum_barl_[kBarlRings][kBarlWedges][kSides];
  double esum_endc_[kEndcWedgesX][kEndcWedgesX][kSides];
 
  double esumMean_barl_[kBarlRings];
  double esumMean_endc_[kEndcEtaRings];
 
  double k_barl_[kBarlRings];
  double k_endc_[kEndcEtaRings];
 
  // calibration const not corrected for k
  float rawconst_barl[kBarlRings][kBarlWedges][kSides];
  float rawconst_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];
 
  EcalGeomPhiSymHelper e_;
 
  std::vector<DetId> barrelCells;
  std::vector<DetId> endcapCells;
 
  bool firstpass_;
  int statusThreshold_;
 
  bool reiteration_;
  std::string oldcalibfile_;
 
  EcalIntercalibConstants oldCalibs_;
 
  EcalIntercalibConstants newCalibs_;
 
  EcalIntercalibConstants miscalib_;
 
  bool have_initial_miscalib_;
  std::string initialmiscalibfile_;
 
  std::vector<TH1F*> miscal_resid_barl_histos;
  std::vector<TH2F*> correl_barl_histos;
 
  std::vector<TH1F*> miscal_resid_endc_histos;
  std::vector<TH2F*> correl_endc_histos;
};