EcalTPGParamBuilder.h

Go to the documentation of this file.

#ifndef ECALTPGPARAMBUILDER_H
#define 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/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/Framework/interface/ESHandle.h"

#include "Geometry/CaloTopology/interface/EcalTrigTowerConstituentsMap.h"

#include "CondFormats/EcalObjects/interface/EcalIntercalibConstants.h"
#include "CondFormats/EcalObjects/interface/EcalGainRatios.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 <TH1F.h>

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

class CaloSubdetectorGeometry;
class EcalElectronicsMapping;
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::EDAnalyzer {
public:
  explicit EcalTPGParamBuilder(edm::ParameterSet const& pSet);
  ~EcalTPGParamBuilder() override;
  void analyze(const edm::Event& evt, const edm::EventSetup& evtSetup) override;
  void beginJob() override;
  bool checkIfOK(EcalPedestals::Item item);

private:
  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 getGCTRegionPhi(int ttphi);
  int getGCTRegionEta(int tteta);
  std::string getDet(int tcc);
  std::pair<std::string, int> getCrate(int tcc);

  const CaloSubdetectorGeometry* theEndcapGeometry_;
  const CaloSubdetectorGeometry* theBarrelGeometry_;
  edm::ESHandle<EcalTrigTowerConstituentsMap> eTTmap_;
  const EcalElectronicsMapping* theMapping_;

  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_;
  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_;

  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 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_;
  std::string tag_;
  int version_;
  int m_write_ped;
  int m_write_lin;
  int m_write_lut;
  int m_write_wei;
  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];

  bool useDBShape_;
};
#endif