EnergySysts.cxx

Go to the documentation of this file.

#include "3FlavorAna/Systs/EnergySysts.h"

#include "StandardRecord/Proxy/SRProxy.h"

#include "TRandom3.h"

namespace ana
{
  //----------------------------------------------------------------------
  const HadEnergyScaleSyst kHadEnergyScaleSyst(.05);
  const MuEnergyScaleSyst kMuEnergyScaleSyst(.02);
  const FDHadEnergyScaleSyst kFDHadEnergyScaleSyst(.05);
  const FDMuEnergyScaleSyst kFDMuEnergyScaleSyst(.01);

  //----------------------------------------------------------------------
  void HadEnergyScaleSyst::
  Shift(double sigma, caf::SRProxy* sr, double& weight) const
  {
    const double hadE = sr->energy.numu.recotrkcchadE;
    const double hadEnew = (1-fScale*sigma)*hadE;
    const double difference = hadEnew - hadE;

    sr->energy.numu.trkccE += difference;
    sr->energy.numu.E += difference;
  }

  //----------------------------------------------------------------------
  void MuEnergyScaleSyst::
  Shift(double sigma, caf::SRProxy* sr, double& weight) const
  {
    const double muE = sr->energy.numu.recomuonE;
    const double muEnew = (1-fScale*sigma)*muE;
    const double difference = muEnew - muE;

    sr->energy.numu.recomuonE = muEnew;
    sr->energy.numu.trkccE += difference;
    sr->energy.numu.E += difference;
  }

  //----------------------------------------------------------------------
  void FDHadEnergyScaleSyst::
  Shift(double sigma, caf::SRProxy* sr, double& weight) const
  {
    if (sr->hdr.det != caf::kFARDET) return;
    HadEnergyScaleSyst::Shift( sigma, sr, weight);
  }

  //----------------------------------------------------------------------
  /// FD muon energy scale systematic
  void FDMuEnergyScaleSyst::
  Shift(double sigma, caf::SRProxy* sr, double& weight) const
  {
    if (sr->hdr.det != caf::kFARDET) return;
    MuEnergyScaleSyst::Shift( sigma, sr, weight);
  }
}