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