d8/d88/EcalTPGParamBuilder_8h_source.html

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

  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