MissingMassInput.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
// Class handling the input of the MissingMassCalculator
// author Michael Huebner <michael.huebner@no.spam.cern.ch>
 
#ifndef DITAUMASSTOOLS_MISSINGMASSINPUT_H
#define DITAUMASSTOOLS_MISSINGMASSINPUT_H
 
// EDM include(s):
#include "xAODMissingET/MissingET.h"
 
// Local include(s):
#include "DiTauMassTools/HelperFunctions.h"
 
namespace DiTauMassTools{
  using ROOT::Math::PtEtaPhiMVector;
  using ROOT::Math::XYVector;
 
class MissingMassInput {
 
  public:
    MissingMassInput();
    ~MissingMassInput();
 
    void ClearInput();
    void PrintInputInfo();
 
    void SetBeamEnergy(double val) { m_beamEnergy=val; } // beam energy
    void SetMETresSyst(int val) { m_METresSyst=val; } // MET resolution systematics: +/-1 sigma
    void SetMetScanParamsUE(double sumEt, double phi_scan=0.0, int data_code=0);
    void SetNjet25(int val);
    void SetSumEt(double sumEt);
    void SetVisTauVec(int i, const PtEtaPhiMVector & vec);
    void SetVisTauType(int i, int tautype);
    void SetMetVec(const XYVector & vec);
    void SetNprong(int i, int nprong);
    void SetIsData(int val);
    void SetMHtSigma1(double val) { m_MHtSigma1=val; }
    void SetMHtSigma2(double val) { m_MHtSigma2=val; } 
    void SetMHtGaussFr(double val) { m_MHtGaussFr=val; }
    void SetMETScanScheme(int val) { m_METScanScheme=val; }
    void SetUseDefaults(int val) { m_fUseDefaults=val; }
    void SetUseTailCleanup(bool val) { m_fUseTailCleanup=val; }
    void SetUseVerbose(bool val) { m_fUseVerbose=val; }
    void SetLFVmode(int val) { m_LFVmode=val; }
 
    friend class MissingMassCalculator;
    friend class MissingMassProb;
 
  private:
    XYVector m_MetVec;
    PtEtaPhiMVector m_vistau1;
    PtEtaPhiMVector m_vistau2;
    int m_type_visTau1{}; // 0: l, 1:1-prong, 3:3-prong
    int m_type_visTau2{};  // 0: l, 1:1-prong, 3:3-prong
    int m_Nprong_tau1{};
    int m_Nprong_tau2{};
    int m_dataType{};
    double m_METcovphi{};
    double m_METsigmaP{};
    double m_METsigmaL{};
    double m_SumEt{};
    std::vector<PtEtaPhiMVector> m_jet4vecs;
    int m_Njet25{};
    double m_DelPhiTT{};
    double m_MHtSigma1{}; // sigma of 1st Gaussian in missing Ht resolution
    double m_MHtSigma2{}; // sigma of 2nd Gaussian in missing Ht resolution
    double m_MHtGaussFr{}; // relative fraction of 2nd Gaussian
    double m_HtOffset{}; // HT offset
    double m_beamEnergy{}; // beam energy (Tevatron=980, LHC-1=3500.0) 
    int m_InputReorder{}; // flag for input re-order
    int m_METScanScheme{}; // MET-scan scheme: 0- use JER{}; 1- use simple sumEt & missingHt for Njet=0 events in (lep-had) 
    int m_METresSyst{}; // switch to turn on/off MET resolution systematics
    bool m_fUseDefaults{}; // switch to control defaults: 1== use defaults, 0== don't use defaults (useful for studies) 
    bool m_fUseTailCleanup{}; // switch to apply tail clean-up
    bool m_fUseVerbose{}; // code to turn ON printouts for debugging
    XYVector m_metVec;
    PtEtaPhiMVector m_tlv_tmp;
    double m_inputMEtX{},m_inputMEtY{},m_inputMEtT{};
    double m_MEtX{}, m_MEtY{},m_MEtT{};
    double m_htOffset{};
    int m_LFVmode{}; // flag to determine which LFV decay to be reconstructed: 0=H->e+tau(mu) or 1=H->mu+tau(e) 
    TauTypes::e m_tauTypes{};
};
} // namespace DiTauMassTools
 
#endif