WrongSignSysts.h

Go to the documentation of this file.

#pragma once

#include "CAFAna/Core/ISyst.h"

namespace ana
{
  /// Wrong sign scale systematic up and dn by 30%
  class WrongSignScale: public ISyst
  {
  public:
    WrongSignScale() : ISyst("ScaleWS", "Wrong Sign Scale Systematic") {}

    // apply 30% scale to ws
    void Shift(double sigma,
               caf::SRProxy* sr, double& weight) const override;
  };

  extern const WrongSignScale kWrongSignScale;

  //-----------------------------------------------------
  /// Wrong sign scale RHC only systematic by energy bin
  /// only up or dn shift
  class WrongSignEnergyBinRHC: public ISyst
  {
  public:
    WrongSignEnergyBinRHC() : ISyst("WrongSignEnergyBinRHC", "Scale RHC Wrong Sign by Bin (Up or Dn)") {}

    void Shift(double sigma,
               caf::SRProxy* sr, double& weight) const override;
  };

  extern const WrongSignEnergyBinRHC kWrongSignEnergyBinRHC;

  //-----------------------------------------------------
  /// Wrong sign scale systematic by energy bin
  /// only up or dn shift
  class WrongSignEnergyBin: public ISyst
  {
  public:
    WrongSignEnergyBin() : ISyst("WrongSignEnergyBin", "Scale Wrong Sign by Bin") {}

    void Shift(double sigma,
               caf::SRProxy* sr, double& weight) const override;
  };

  extern const WrongSignEnergyBin kWrongSignEnergyBin;

  //-----------------------------------------------------
  class WrongSignScale100: public ISyst
  {
  public:
    WrongSignScale100() : ISyst("ScaleWrongSign100", "Wrong Sign Scale Systematic 100") {}

    void Shift(double sigma,
               caf::SRProxy* sr, double& weight) const override;
  };

  extern const WrongSignScale100 kWrongSignScale100;

  //-------------------------------------------------------------
  /// Wrong sign scale systematic set by the user
  /// scalefhc scales the fhc ws
  /// scalerhc scales the rhc ws
  class MyWrongSignScale: public ISyst
  {
  public:
    MyWrongSignScale(const std::string& shortName,
                     const std::string& latexName,
                     double scaleFHC,
                     double scaleRHC)
      : ISyst(shortName, latexName),
        fScaleFHC(scaleFHC), fScaleRHC(scaleRHC)
    {
    }

    MyWrongSignScale(const double scaleFHC, const double scaleRHC)
      : ISyst("MyWrongSignScale", "My Wrong Sign Scale Systematic"),
        fScaleFHC(scaleFHC),
        fScaleRHC(scaleRHC)
    {
    }


    void Shift(double sigma,
               caf::SRProxy* sr, double& weight) const override;

  private:

    const double fScaleFHC;
    const double fScaleRHC;
  };

  extern const MyWrongSignScale kMyWrongSignScale;
}