FakeDecomp.h

Go to the documentation of this file.

//********************************************************
//Docdb:30405
//Not a real Decomp for now. As I hate to change the 
//following Pred/Extrap frame for a flexible Type.
//********************************************************


#pragma once

#include "CAFAna/Decomp/IDecomp.h"

#include "CAFAna/Core/Loaders.h"

namespace ana
{

  namespace DecompType{
    enum Type{kNue, kAntiNue, kNumu, kAntiNumu, kNC, kAntiNC, kNCTotal, kNum};
  }

  /// Uses MC for \f$ \nu_\mu \f$ CC and NC components, assigns remainder of data to \f$ \nu_e \f$ CC
  class FakeDecomp: public IDecomp
  {
  public:
    
    FakeDecomp(SpectrumLoaderBase& loaderMC,
               SpectrumLoaderBase& loaderData,
               const std::string& label,
               const Binning& bins,
               const Var& var,
               const Cut& cut,
               const SystShifts& shiftMC = kNoShift,
               const SystShifts& shiftData = kNoShift,
               const Var& wei = kUnweighted);

    FakeDecomp(SpectrumLoaderBase& loaderMC,
               SpectrumLoaderBase& loaderData,
               const HistAxis& axis,
               const Cut& cut,
               const SystShifts& shiftMC = kNoShift,
               const SystShifts& shiftData = kNoShift,
               const Var& wei = kUnweighted);

    FakeDecomp(Loaders& loaders,
               const HistAxis& axis,
               const Cut& cut,
               const SystShifts& shiftMC = kNoShift,
               const SystShifts& shiftData = kNoShift,
               const Var& wei = kUnweighted);

    Spectrum NumuComponent()     const override; 
    Spectrum AntiNumuComponent() const override; 
    Spectrum NCTotalComponent()       const override; 
    Spectrum NCComponent()            const override; 
    Spectrum NCAntiComponent()        const override;
    Spectrum NueComponent()      const override; 
    Spectrum AntiNueComponent()  const override;

    virtual ~FakeDecomp();

    void Decompose () const;

    Spectrum Data_Component()                const override {return fData;}
    virtual Spectrum Pre_MCTotal()           const{return fTotal;}
    virtual Spectrum MCTotal()               const{return NumuComponent()+AntiNumuComponent()+NCComponent()+NCAntiComponent()+NueComponent()+AntiNueComponent();}
    virtual Spectrum Pre_NumuComponent()     const{return fNumu;}
    virtual Spectrum Pre_AntiNumuComponent() const{return fAntiNumu;}
    virtual Spectrum Pre_NCTotalComponent()       const{return fNCTot;}
    virtual Spectrum Pre_NCComponent()            const{return fNC;}
    virtual Spectrum Pre_NCAntiComponent()        const{return fNCAnti;}
    virtual Spectrum Pre_NueComponent()      const{return fNue;}
    virtual Spectrum Pre_AntiNueComponent()  const{return fAntiNue;}

    void SaveTo(TDirectory* dir, const std::string& name) const override;
    void SetType(DecompType::Type type)  const;

    static std::unique_ptr<FakeDecomp> LoadFrom(TDirectory* dir, const std::string& name);

  protected:
    FakeDecomp()
      : fData(Spectrum::Uninitialized()),
        fTotal(Spectrum::Uninitialized()),
        fNCTot  (Spectrum::Uninitialized()),
        fNC     (Spectrum::Uninitialized()),
        fNCAnti (Spectrum::Uninitialized()),
        fNue (Spectrum::Uninitialized()),
        fAntiNue(Spectrum::Uninitialized()),
        fNumu(Spectrum::Uninitialized()),
        fAntiNumu(Spectrum::Uninitialized()),
        fIsDecomposed(0)
    {};

    Spectrum fData;
    Spectrum fTotal;

    Spectrum fNCTot;
    Spectrum fNC;
    Spectrum fNCAnti;
    
    Spectrum fNue;
    Spectrum fAntiNue;

    Spectrum fNumu;
    Spectrum fAntiNumu;

// 1 ***********
    mutable DecompType::Type fType = DecompType::Type::kNCTotal;
    mutable bool fIsSet = false;
// *************

// 2 ***********
//    mutable DecompType::Type fType = DecompType::Type::kAntiNue;
//    mutable bool fIsSet = true;
//**************
    mutable bool fIsDecomposed = false;
    mutable std::vector<Spectrum*> fDiff;

//    mutable Spectrum fSpecCache;

  };
}