ShowerCalo.cxx

Go to the documentation of this file.

#include "ShowerCalo.h"

#include "TString.h"
#include "larreco/ShowerFinder/ShowerReco3D/ShowerRecoException.h"
#include <string>

namespace showerreco {

  std::pair<double,double> ShowerCalo::ShowerEnergy(double energy,
                                                    ::showerreco::energy::RecoType_t type,
                                                    bool mc)
  {
    double result = -1;
    double result_err = -1;

    switch(type) {

    case ::showerreco::energy::kDEFAULT:

      if(!mc) throw ShowerRecoException(Form("DATA correction not yet supported!"));

      result = ::showerreco::energy::DEFAULT_ECorr * energy + 10;
      result_err = ::showerreco::energy::DEFAULT_ECorrError * result;
      break;

    default:

      throw ShowerRecoException(Form("RecoType_t: %d not supported!",type));
    }

    return std::make_pair(result,result_err);
  }

}