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;

};