MichelTaggingSyst.h

Go to the documentation of this file.

#pragma once

#include "CAFAna/Core/ISyst.h"

#include "TRandom3.h"

namespace ana
{
  // 2020 version -- easiest way to just copy the old version to avoid naming conflicts and not destroying the old pieces of code
  class MichelTaggingSyst2020 : public ISyst
  {
  public:
    MichelTaggingSyst2020() : ISyst("michel_tagging2020", "Michel Electrons Tagging Uncertainty") {}

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

  private:
    // Jeremy's trick for random numbers
    int Hash(caf::SRProxy* sr) const;

    // 2020 first preliminary numbers, waiting for full validation
    double fProbUP = .08;
    double fProbDO = .05;

    mutable TRandom3 fMRnd;
  };
  
  extern const MichelTaggingSyst2020 kMichelTaggingSyst2020;
  
  //----------------------------------------------------------------------------------------
  
  
  // 2018 version
  class MichelTaggingSyst2018 : public ISyst
  {
  public:
    MichelTaggingSyst2018() : ISyst("michel_tagging2018", "Michel Electrons Tagging Uncertainty") {}

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

  private:
    // Jeremy's trick for random numbers
    int Hash(caf::SRProxy* sr) const;
    
    // made up numbers to hack +-5% shift of tagging efficiency
    // this needs to be treated carefully and tested for every new/old datasets
    // it can change quite a lot
    // DO NOT USE FOR ANY OTHER DATASETS THAN 2018 (prod4) FHC
    double fProbUP = .36;
    double fProbDO = .08;

    mutable TRandom3 fMRnd;
  };

  extern const MichelTaggingSyst2018 kMichelTaggingSyst2018;
}