EcalTPGParamBuilder.h

Go to the documentation of this file.

#ifndef CalibCalorimetry_EcalTPGTools_EcalTPGParamBuilder_h
#define CalibCalorimetry_EcalTPGTools_EcalTPGParamBuilder_h

//Author: Pascal Paganini - LLR
//Date: 2006/07/10 15:58:06 $

// system include files
#include <memory>

// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Utilities/interface/ESGetToken.h"
#include "FWCore/Framework/interface/ESHandle.h"

#include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
#include "Geometry/CaloTopology/interface/EcalTrigTowerConstituentsMap.h"
#include "Geometry/EcalMapping/interface/EcalElectronicsMapping.h"
#include "Geometry/EcalMapping/interface/EcalMappingRcd.h"
#include "Geometry/Records/interface/CaloGeometryRecord.h"

#include "CondFormats/DataRecord/interface/EcalADCToGeVConstantRcd.h"
#include "CondFormats/DataRecord/interface/EcalIntercalibConstantsRcd.h"
#include "CondFormats/DataRecord/interface/EcalGainRatiosRcd.h"
#include "CondFormats/DataRecord/interface/EcalLaserAlphasRcd.h"
#include "CondFormats/DataRecord/interface/EcalLaserAPDPNRatiosRcd.h"
#include "CondFormats/DataRecord/interface/EcalPedestalsRcd.h"
#include "CondFormats/EcalObjects/interface/EcalADCToGeVConstant.h"
#include "CondFormats/EcalObjects/interface/EcalIntercalibConstants.h"
#include "CondFormats/EcalObjects/interface/EcalGainRatios.h"
#include "CondFormats/EcalObjects/interface/EcalLaserAPDPNRatios.h"
#include "CondFormats/EcalObjects/interface/EcalPedestals.h"
#include "OnlineDB/EcalCondDB/interface/all_monitoring_types.h"
#include "OnlineDB/EcalCondDB/interface/all_fe_config_types.h"

//modif-alex-27-july-2015
#include "CondFormats/EcalObjects/interface/EcalLaserAlphas.h"

#include "SimCalorimetry/EcalSimAlgos/interface/EcalShapeBase.h"
#include "SimCalorimetry/EcalSimAlgos/interface/EBShape.h"
#include "SimCalorimetry/EcalSimAlgos/interface/EEShape.h"

#include <TH1F.h>

#include <vector>
#include <string>
#include <map>
#include <iostream>

class EcalTPGDBApp;

class coeffStruc {
public:
  coeffStruc() {}
  double calibCoeff_;
  double gainRatio_[3];
  int pedestals_[3];
};

class linStruc {
public:
  linStruc() {}
  int pedestal_[3];
  int mult_[3];
  int shift_[3];
};

class EcalTPGParamBuilder : public edm::one::EDAnalyzer<> {
public:
  explicit EcalTPGParamBuilder(edm::ParameterSet const& pSet);
  ~EcalTPGParamBuilder() override;
  void analyze(const edm::Event& evt, const edm::EventSetup& evtSetup) override;
  void beginJob() override;

private:
  bool checkIfOK(EcalPedestals::Item item);
  bool computeLinearizerParam(
      double theta, double gainRatio, double calibCoeff, std::string subdet, int& mult, int& shift);
  void create_header();
  int uncodeWeight(double weight, int complement2 = 7);
  double uncodeWeight(int iweight, int complement2 = 7);
  std::vector<unsigned int> computeWeights(EcalShapeBase& shape, TH1F* histo);
  void computeLUT(int* lut, std::string det = "EB");
  //void getCoeff(coeffStruc & coeff, const EcalIntercalibConstantMap & calibMap, uint rawId) ; //modif-alex-27-july-2015 uncomment to go back
  void getCoeff(coeffStruc& coeff, const EcalGainRatioMap& gainMap, uint rawId);
  void getCoeff(coeffStruc& coeff, const EcalPedestalsMap& pedMap, uint rawId);
  void getCoeff(coeffStruc& coeff, const std::map<EcalLogicID, MonPedestalsDat>& pedMap, const EcalLogicID& logicId);

  //modif-alex-27-july-2015
  void getCoeff(coeffStruc& coeff,
                const EcalIntercalibConstantMap& calibMap,
                const EcalLaserAlphaMap& laserAlphaMap,
                uint rawId,
                std::string& ss);

  void computeFineGrainEBParameters(uint& lowRatio, uint& highRatio, uint& lowThreshold, uint& highThreshold, uint& lut);
  void computeFineGrainEEParameters(uint& threshold, uint& lut_strip, uint& lut_tower);
  int getEtaSlice(int tccId, int towerInTCC);
  bool realignBaseline(linStruc& lin, float forceBase12, int cmsswid);
  int getGCTRegionPhi(int ttphi);
  int getGCTRegionEta(int tteta);
  std::string getDet(int tcc);
  std::pair<std::string, int> getCrate(int tcc);

  edm::ESGetToken<CaloSubdetectorGeometry, EcalEndcapGeometryRecord> theEndcapGeometryToken_;
  edm::ESGetToken<CaloSubdetectorGeometry, EcalBarrelGeometryRecord> theBarrelGeometryToken_;
  const CaloSubdetectorGeometry* theEndcapGeometry_;
  const CaloSubdetectorGeometry* theBarrelGeometry_;
  edm::ESGetToken<EcalTrigTowerConstituentsMap, IdealGeometryRecord> eTTmapToken_;
  const EcalTrigTowerConstituentsMap* eTTmap_;
  edm::ESGetToken<EcalElectronicsMapping, EcalMappingRcd> ecalmappingToken_;
  const EcalElectronicsMapping* theMapping_;

  edm::ESGetToken<EcalLaserAlphas, EcalLaserAlphasRcd> ecalLaserAlphasToken_;
  edm::ESGetToken<EcalLaserAPDPNRatios, EcalLaserAPDPNRatiosRcd> ecalLaserAPDPNRatiosToken_;
  edm::ESGetToken<EcalPedestals, EcalPedestalsRcd> ecalPedestalsToken_;
  edm::ESGetToken<EcalIntercalibConstants, EcalIntercalibConstantsRcd> ecalIntercalibConstantsToken_;
  edm::ESGetToken<EcalGainRatios, EcalGainRatiosRcd> ecalGainRatiosToken_;
  edm::ESGetToken<EcalADCToGeVConstant, EcalADCToGeVConstantRcd> ecalADCToGeVConstantToken_;

  bool useTransverseEnergy_;
  double xtal_LSB_EB_, xtal_LSB_EE_;
  double Et_sat_EB_, Et_sat_EE_;
  unsigned int sliding_;
  unsigned int sampleMax_;
  double weight_timeShift_;
  bool weight_unbias_recovery_;
  bool weight_even_computeFromShape_;
  std::string weight_even_idMapFile_;
  std::string weight_even_weightGroupFile_;
  std::string weight_odd_idMapFile_;
  std::string weight_odd_weightGroupFile_;
  static constexpr uint NWEIGHTS = 5;  // It can be up to 6, but in Run 1,2,3 only 5 weights have been used

  bool weight_useDoubleWeights_;
  bool TPmode_EnableEBOddFilter_;
  bool TPmode_EnableEEOddFilter_;
  bool TPmode_EnableEBOddPeakFinder_;
  bool TPmode_EnableEEOddPeakFinder_;
  bool TPmode_DisableEBEvenPeakFinder_;
  bool TPmode_DisableEEEvenPeakFinder_;
  unsigned int TPmode_FenixEBStripOutput_;
  unsigned int TPmode_FenixEEStripOutput_;
  unsigned int TPmode_FenixEBStripInfobit2_;
  unsigned int TPmode_FenixEEStripInfobit2_;
  unsigned int TPmode_FenixEBTcpOutput_;
  unsigned int TPmode_FenixEBTcpInfobit1_;
  unsigned int TPmode_FenixEETcpOutput_;
  unsigned int TPmode_FenixEETcpInfobit1_;

  unsigned int nSample_;
  unsigned int complement2_;
  std::string LUT_option_;
  double LUT_threshold_EB_, LUT_threshold_EE_;
  double LUT_stochastic_EB_, LUT_noise_EB_, LUT_constant_EB_;
  double LUT_stochastic_EE_, LUT_noise_EE_, LUT_constant_EE_;
  double TTF_lowThreshold_EB_, TTF_highThreshold_EB_;
  double TTF_lowThreshold_EE_, TTF_highThreshold_EE_;
  double FG_lowThreshold_EB_, FG_highThreshold_EB_, FG_lowRatio_EB_, FG_highRatio_EB_;
  unsigned int FG_lut_EB_;
  double FG_Threshold_EE_;
  unsigned int FG_lut_strip_EE_, FG_lut_tower_EE_;
  int forcedPedestalValue_;
  bool forceEtaSlice_;
  unsigned int SFGVB_Threshold_, SFGVB_lut_, pedestal_offset_;
  int SFGVB_SpikeKillingThreshold_;  //modif-alex 01/21/11
  bool useInterCalibration_, H2_;

  //modif-alex-30/01/2012
  std::string Transparency_Corr_;
  bool useTransparencyCorr_;

  //modif-alex-02/02/11
  std::string TimingDelays_EB_;
  std::string TimingDelays_EE_;
  std::string TimingPhases_EB_;
  std::string TimingPhases_EE_;
  std::map<int, std::vector<int> > delays_EB_;
  std::map<int, std::vector<int> > phases_EB_;
  std::map<int, std::vector<int> > delays_EE_;
  std::map<int, std::vector<int> > phases_EE_;

  //modif-alex 30/01/2012
  std::map<int, double> Transparency_Correction_;

  EBShape shapeEB_;
  EEShape shapeEE_;

  std::ofstream* out_file_;
  std::ofstream* geomFile_;
  EcalTPGDBApp* db_;
  bool writeToDB_;
  bool writeToFiles_;
  unsigned int DBrunNb_;
  bool DBEE_;

  int ped_conf_id_;
  int lin_conf_id_;
  int lut_conf_id_;
  int fgr_conf_id_;
  int sli_conf_id_;
  int wei_conf_id_;
  int wei2_conf_id_;
  int spi_conf_id_;  //modif-alex 21/01.11
  int del_conf_id_;  //modif-alex 21/01.11
  int bxt_conf_id_;
  int btt_conf_id_;
  int bst_conf_id_;
  int coke_conf_id_;
  std::string tag_;
  int version_;
  int m_write_ped;
  int m_write_lin;
  int m_write_lut;
  int m_write_wei;
  int m_write_wei2;
  int m_write_fgr;
  int m_write_sli;
  int m_write_spi;  //modif-alex 21/01/11
  int m_write_del;  //modif-alex 21/01/11
  int m_write_bxt;
  int m_write_btt;
  int m_write_bst;

  Int_t* ntupleInts_;
  Char_t ntupleDet_[10];
  Char_t ntupleCrate_[10];
};
#endif