dd/d5a/MissingMassInput_8h_source.html

/*

  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