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();
  };

  virtual ~EnergyResolutionVsLumi();
      
  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