EnergyResolutionVsLumi.h

Go to the documentation of this file.

#ifndef CalibCalorimetry_EnergyResolutionVsLumi_H
#define CalibCalorimetry_EnergyResolutionVsLumi_H
// system include files

#include <vector>
#include <typeinfo>
#include <string>
#include <map>

#include <ctime>
#include <cstdio>

#include <cmath>

#include "SimG4CMS/Calo/interface/EvolutionECAL.h"
#include "DataFormats/DetId/interface/DetId.h"

class EnergyResolutionVsLumi {
public:
  EnergyResolutionVsLumi();

  EnergyResolutionVsLumi(double lumi, double instlumi) {
    m_lumi = lumi;
    m_instlumi = instlumi;
    calcmuTot();
  };

  ~EnergyResolutionVsLumi() = default;

  struct DegradationAtEta {
    double eta;
    double muEM;
    double muHD;
    double ampDropTransparency;
    double ampDropPhotoDetector;
    double ampDropTotal;
    double noiseIncreaseADC;
    double resolutitonConstantTerm;
  };

  DegradationAtEta CalculateDegradation(double eta);
  double Resolution(double eta, double ene);
  //void Decomposition();

  void setLumi(double x) { m_lumi = x; };
  void setInstLumi(double x) { m_instlumi = x; };
  void setLumies(double x, double y) { m_lumi = x, m_instlumi = y, calcmuTot(); };

  double calcmuEM(double eta);
  double calcmuHD(double eta);
  double calcampDropTransparency(double eta);
  double calcampDropPhotoDetector(double eta);
  double calcampDropTotal(double eta);
  double calcnoiseIncreaseADC(double eta);
  double calcnoiseADC(double eta);
  double calcresolutitonConstantTerm(double eta);

  double calcLightCollectionEfficiencyWeighted(DetId id, double z);

  double calcLightCollectionEfficiencyWeighted2(double eta, double z, double mu_ind = -1.0);
  double calcmuTot(double eta);
  void calcmuTot();
  double getmuTot(double eta, int ix, int iy);

private:
  double m_lumi;
  double m_instlumi;
  double mu_eta[10085];
  double vpt_eta[10085];
};

#endif