ElectronEnergyCalibrator.h

Go to the documentation of this file.

#ifndef ElectronEnergyCalibrator_H
#define ElectronEnergyCalibrator_H

#include "EgammaAnalysis/ElectronTools/interface/SimpleElectron.h"
#include <string>
#include <vector>
#include <fstream>
#include <sstream>
#include <iostream>

namespace edm {
  class StreamID;
}

struct correctionValues
{
    double nRunMin;
    double nRunMax;
    double corrCat0;
    double corrCat1;
    double corrCat2;
    double corrCat3;
    double corrCat4;
    double corrCat5;
    double corrCat6;
    double corrCat7;
};

struct linearityCorrectionValues
{
    double ptMin;
    double ptMax;
    double corrCat0;
    double corrCat1;
    double corrCat2;
    double corrCat3;
    double corrCat4;
    double corrCat5;
};

class ElectronEnergyCalibrator
{
 public:
 ElectronEnergyCalibrator(const std::string pathData, 
              const std::string pathLinData,
              const std::string dataset, 
              int correctionsType, 
              bool applyLinearityCorrection, 
              double lumiRatio, 
              bool isMC, 
              bool updateEnergyErrors, 
              bool verbose, 
              bool synchronization
              ) : 
    pathData_(pathData), 
    pathLinData_(pathLinData), 
    dataset_(dataset), 
    correctionsType_(correctionsType), 
    applyLinearityCorrection_(applyLinearityCorrection),
    lumiRatio_(lumiRatio), 
    isMC_(isMC), 
    updateEnergyErrors_(updateEnergyErrors), 
    verbose_(verbose), 
    synchronization_(synchronization) {
      init();
    }

    void calibrate(SimpleElectron &electron, edm::StreamID const&);
    void correctLinearity(SimpleElectron &electron);
                  
 private:
    void init();
    void splitString( const std::string &fullstr, 
              std::vector<std::string> &elements, 
              const std::string &delimiter
              );
    double stringToDouble(const std::string &str);
      
    double newEnergy_ ;
    double newEnergyError_ ;
    
    std::string pathData_;
    std::string pathLinData_;
    std::string dataset_;
    int correctionsType_;
    bool applyLinearityCorrection_;
    double lumiRatio_;
    bool isMC_;
    bool updateEnergyErrors_;
    bool verbose_;
    bool synchronization_;
    
    correctionValues corrValArray[100];
    correctionValues corrValMC;
    linearityCorrectionValues linCorrValArray[100];
    int nCorrValRaw, nLinCorrValRaw;
};

#endif