HiggsTemplateCrossSections.h

Go to the documentation of this file.

#ifndef SimDataFormats_GeneratorProducts_HiggsTemplateCrossSections_h
#define SimDataFormats_GeneratorProducts_HiggsTemplateCrossSections_h

#include "DataFormats/Math/interface/LorentzVector.h"
#include <vector>

namespace HTXS {

  enum ErrorCode {
    UNDEFINED = -99,
    SUCCESS = 0,                     
    PRODMODE_DEFINED = 1,            
    MOMENTUM_CONSERVATION = 2,       
    HIGGS_IDENTIFICATION = 3,        
    HIGGS_DECAY_IDENTIFICATION = 4,  
    HS_VTX_IDENTIFICATION = 5,       
    VH_IDENTIFICATION = 6,           
    VH_DECAY_IDENTIFICATION = 7,     
    TOP_W_IDENTIFICATION = 8         
  };

  enum HiggsProdMode { UNKNOWN = 0, GGF = 1, VBF = 2, WH = 3, QQ2ZH = 4, GG2ZH = 5, TTH = 6, BBH = 7, TH = 8 };

  enum tH_type { noTH = 0, THQB = 1, TWH = 2 };


  namespace Stage0 {
    enum Category {
      UNKNOWN = 0,
      GG2H_FWDH = 10,
      GG2H = 11,
      VBF_FWDH = 20,
      VBF = 21,
      VH2HQQ_FWDH = 22,
      VH2HQQ = 23,
      QQ2HLNU_FWDH = 30,
      QQ2HLNU = 31,
      QQ2HLL_FWDH = 40,
      QQ2HLL = 41,
      GG2HLL_FWDH = 50,
      GG2HLL = 51,
      TTH_FWDH = 60,
      TTH = 61,
      BBH_FWDH = 70,
      BBH = 71,
      TH_FWDH = 80,
      TH = 81
    };
  }  // namespace Stage0

  namespace Stage1 {
    enum Category {
      UNKNOWN = 0,
      // Gluon fusion
      GG2H_FWDH = 100,
      GG2H_VBFTOPO_JET3VETO = 101,
      GG2H_VBFTOPO_JET3 = 102,
      GG2H_0J = 103,
      GG2H_1J_PTH_0_60 = 104,
      GG2H_1J_PTH_60_120 = 105,
      GG2H_1J_PTH_120_200 = 106,
      GG2H_1J_PTH_GT200 = 107,
      GG2H_GE2J_PTH_0_60 = 108,
      GG2H_GE2J_PTH_60_120 = 109,
      GG2H_GE2J_PTH_120_200 = 110,
      GG2H_GE2J_PTH_GT200 = 111,
      // "VBF"
      QQ2HQQ_FWDH = 200,
      QQ2HQQ_VBFTOPO_JET3VETO = 201,
      QQ2HQQ_VBFTOPO_JET3 = 202,
      QQ2HQQ_VH2JET = 203,
      QQ2HQQ_REST = 204,
      QQ2HQQ_PTJET1_GT200 = 205,
      // qq -> WH
      QQ2HLNU_FWDH = 300,
      QQ2HLNU_PTV_0_150 = 301,
      QQ2HLNU_PTV_150_250_0J = 302,
      QQ2HLNU_PTV_150_250_GE1J = 303,
      QQ2HLNU_PTV_GT250 = 304,
      // qq -> ZH
      QQ2HLL_FWDH = 400,
      QQ2HLL_PTV_0_150 = 401,
      QQ2HLL_PTV_150_250_0J = 402,
      QQ2HLL_PTV_150_250_GE1J = 403,
      QQ2HLL_PTV_GT250 = 404,
      // gg -> ZH
      GG2HLL_FWDH = 500,
      GG2HLL_PTV_0_150 = 501,
      GG2HLL_PTV_GT150_0J = 502,
      GG2HLL_PTV_GT150_GE1J = 503,
      // ttH
      TTH_FWDH = 600,
      TTH = 601,
      // bbH
      BBH_FWDH = 700,
      BBH = 701,
      // tH
      TH_FWDH = 800,
      TH = 801
    };
  }  // namespace Stage1

  namespace Stage1_1 {
    enum Category {
      UNKNOWN = 0,
      // Gluon fusion
      GG2H_FWDH = 100,
      GG2H_PTH_GT200 = 101,
      GG2H_0J_PTH_0_10 = 102,
      GG2H_0J_PTH_GT10 = 103,
      GG2H_1J_PTH_0_60 = 104,
      GG2H_1J_PTH_60_120 = 105,
      GG2H_1J_PTH_120_200 = 106,
      GG2H_GE2J_MJJ_0_350_PTH_0_60 = 107,
      GG2H_GE2J_MJJ_0_350_PTH_60_120 = 108,
      GG2H_GE2J_MJJ_0_350_PTH_120_200 = 109,
      GG2H_MJJ_350_700_PTHJJ_0_25 = 110,
      GG2H_MJJ_350_700_PTHJJ_GT25 = 111,
      GG2H_MJJ_GT700_PTHJJ_0_25 = 112,
      GG2H_MJJ_GT700_PTHJJ_GT25 = 113,
      // "VBF"
      QQ2HQQ_FWDH = 200,
      QQ2HQQ_0J = 201,
      QQ2HQQ_1J = 202,
      QQ2HQQ_MJJ_0_60 = 203,
      QQ2HQQ_MJJ_60_120 = 204,
      QQ2HQQ_MJJ_120_350 = 205,
      QQ2HQQ_MJJ_GT350_PTH_GT200 = 206,
      QQ2HQQ_MJJ_350_700_PTHJJ_0_25 = 207,
      QQ2HQQ_MJJ_350_700_PTHJJ_GT25 = 208,
      QQ2HQQ_MJJ_GT700_PTHJJ_0_25 = 209,
      QQ2HQQ_MJJ_GT700_PTHJJ_GT25 = 210,
      // qq -> WH
      QQ2HLNU_FWDH = 300,
      QQ2HLNU_PTV_0_75 = 301,
      QQ2HLNU_PTV_75_150 = 302,
      QQ2HLNU_PTV_150_250_0J = 303,
      QQ2HLNU_PTV_150_250_GE1J = 304,
      QQ2HLNU_PTV_GT250 = 305,
      // qq -> ZH
      QQ2HLL_FWDH = 400,
      QQ2HLL_PTV_0_75 = 401,
      QQ2HLL_PTV_75_150 = 402,
      QQ2HLL_PTV_150_250_0J = 403,
      QQ2HLL_PTV_150_250_GE1J = 404,
      QQ2HLL_PTV_GT250 = 405,
      // gg -> ZH
      GG2HLL_FWDH = 500,
      GG2HLL_PTV_0_75 = 501,
      GG2HLL_PTV_75_150 = 502,
      GG2HLL_PTV_150_250_0J = 503,
      GG2HLL_PTV_150_250_GE1J = 504,
      GG2HLL_PTV_GT250 = 505,
      // ttH
      TTH_FWDH = 600,
      TTH = 601,
      // bbH
      BBH_FWDH = 700,
      BBH = 701,
      // tH
      TH_FWDH = 800,
      TH = 801
    };
  }  // namespace Stage1_1

  namespace Stage1_1_Fine {
    enum Category {
      UNKNOWN = 0,
      // Gluon fusion
      GG2H_FWDH = 100,
      GG2H_PTH_GT200 = 101,
      GG2H_0J_PTH_0_10 = 102,
      GG2H_0J_PTH_GT10 = 103,
      GG2H_1J_PTH_0_60 = 104,
      GG2H_1J_PTH_60_120 = 105,
      GG2H_1J_PTH_120_200 = 106,
      GG2H_GE2J_MJJ_0_350_PTH_0_60_PTHJJ_0_25 = 107,
      GG2H_GE2J_MJJ_0_350_PTH_60_120_PTHJJ_0_25 = 108,
      GG2H_GE2J_MJJ_0_350_PTH_120_200_PTHJJ_0_25 = 109,
      GG2H_GE2J_MJJ_0_350_PTH_0_60_PTHJJ_GT25 = 110,
      GG2H_GE2J_MJJ_0_350_PTH_60_120_PTHJJ_GT25 = 111,
      GG2H_GE2J_MJJ_0_350_PTH_120_200_PTHJJ_GT25 = 112,
      GG2H_MJJ_350_700_PTHJJ_0_25 = 113,
      GG2H_MJJ_350_700_PTHJJ_GT25 = 114,
      GG2H_MJJ_700_1000_PTHJJ_0_25 = 115,
      GG2H_MJJ_700_1000_PTHJJ_GT25 = 116,
      GG2H_MJJ_1000_1500_PTHJJ_0_25 = 117,
      GG2H_MJJ_1000_1500_PTHJJ_GT25 = 118,
      GG2H_MJJ_GT1500_PTHJJ_0_25 = 119,
      GG2H_MJJ_GT1500_PTHJJ_GT25 = 120,
      // "VBF"
      QQ2HQQ_FWDH = 200,
      QQ2HQQ_0J = 201,
      QQ2HQQ_1J = 202,
      QQ2HQQ_MJJ_0_60_PTHJJ_0_25 = 203,
      QQ2HQQ_MJJ_60_120_PTHJJ_0_25 = 204,
      QQ2HQQ_MJJ_120_350_PTHJJ_0_25 = 205,
      QQ2HQQ_MJJ_0_60_PTHJJ_GT25 = 206,
      QQ2HQQ_MJJ_60_120_PTHJJ_GT25 = 207,
      QQ2HQQ_MJJ_120_350_PTHJJ_GT25 = 208,
      QQ2HQQ_MJJ_350_700_PTHJJ_0_25 = 209,
      QQ2HQQ_MJJ_350_700_PTHJJ_GT25 = 210,
      QQ2HQQ_MJJ_700_1000_PTHJJ_0_25 = 211,
      QQ2HQQ_MJJ_700_1000_PTHJJ_GT25 = 212,
      QQ2HQQ_MJJ_1000_1500_PTHJJ_0_25 = 213,
      QQ2HQQ_MJJ_1000_1500_PTHJJ_GT25 = 214,
      QQ2HQQ_MJJ_GT1500_PTHJJ_0_25 = 215,
      QQ2HQQ_MJJ_GT1500_PTHJJ_GT25 = 216,
      QQ2HQQ_PTH_GT200_MJJ_350_700_PTHJJ_0_25 = 217,
      QQ2HQQ_PTH_GT200_MJJ_350_700_PTHJJ_GT25 = 218,
      QQ2HQQ_PTH_GT200_MJJ_700_1000_PTHJJ_0_25 = 219,
      QQ2HQQ_PTH_GT200_MJJ_700_1000_PTHJJ_GT25 = 220,
      QQ2HQQ_PTH_GT200_MJJ_1000_1500_PTHJJ_0_25 = 221,
      QQ2HQQ_PTH_GT200_MJJ_1000_1500_PTHJJ_GT25 = 222,
      QQ2HQQ_PTH_GT200_MJJ_GT1500_PTHJJ_0_25 = 223,
      QQ2HQQ_PTH_GT200_MJJ_GT1500_PTHJJ_GT25 = 224,
      // qq -> WH
      QQ2HLNU_FWDH = 300,
      QQ2HLNU_PTV_0_75_0J = 301,
      QQ2HLNU_PTV_75_150_0J = 302,
      QQ2HLNU_PTV_150_250_0J = 303,
      QQ2HLNU_PTV_250_400_0J = 304,
      QQ2HLNU_PTV_GT400_0J = 305,
      QQ2HLNU_PTV_0_75_1J = 306,
      QQ2HLNU_PTV_75_150_1J = 307,
      QQ2HLNU_PTV_150_250_1J = 308,
      QQ2HLNU_PTV_250_400_1J = 309,
      QQ2HLNU_PTV_GT400_1J = 310,
      QQ2HLNU_PTV_0_75_GE2J = 311,
      QQ2HLNU_PTV_75_150_GE2J = 312,
      QQ2HLNU_PTV_150_250_GE2J = 313,
      QQ2HLNU_PTV_250_400_GE2J = 314,
      QQ2HLNU_PTV_GT400_GE2J = 315,
      // qq -> ZH
      QQ2HLL_FWDH = 400,
      QQ2HLL_PTV_0_75_0J = 401,
      QQ2HLL_PTV_75_150_0J = 402,
      QQ2HLL_PTV_150_250_0J = 403,
      QQ2HLL_PTV_250_400_0J = 404,
      QQ2HLL_PTV_GT400_0J = 405,
      QQ2HLL_PTV_0_75_1J = 406,
      QQ2HLL_PTV_75_150_1J = 407,
      QQ2HLL_PTV_150_250_1J = 408,
      QQ2HLL_PTV_250_400_1J = 409,
      QQ2HLL_PTV_GT400_1J = 410,
      QQ2HLL_PTV_0_75_GE2J = 411,
      QQ2HLL_PTV_75_150_GE2J = 412,
      QQ2HLL_PTV_150_250_GE2J = 413,
      QQ2HLL_PTV_250_400_GE2J = 414,
      QQ2HLL_PTV_GT400_GE2J = 415,
      // gg -> ZH
      GG2HLL_FWDH = 500,
      GG2HLL_PTV_0_75_0J = 501,
      GG2HLL_PTV_75_150_0J = 502,
      GG2HLL_PTV_150_250_0J = 503,
      GG2HLL_PTV_250_400_0J = 504,
      GG2HLL_PTV_GT400_0J = 505,
      GG2HLL_PTV_0_75_1J = 506,
      GG2HLL_PTV_75_150_1J = 507,
      GG2HLL_PTV_150_250_1J = 508,
      GG2HLL_PTV_250_400_1J = 509,
      GG2HLL_PTV_GT400_1J = 510,
      GG2HLL_PTV_0_75_GE2J = 511,
      GG2HLL_PTV_75_150_GE2J = 512,
      GG2HLL_PTV_150_250_GE2J = 513,
      GG2HLL_PTV_250_400_GE2J = 514,
      GG2HLL_PTV_GT400_GE2J = 515,
      // ttH
      TTH_FWDH = 600,
      TTH = 601,
      // bbH
      BBH_FWDH = 700,
      BBH = 701,
      // tH
      TH_FWDH = 800,
      TH = 801
    };
  }  // namespace Stage1_1_Fine

  //#ifdef ROOT_TLorentzVector
  //typedef TLorentzVector TLV;
  typedef math::XYZTLorentzVectorD TLV;
  typedef std::vector<TLV> TLVs;

  template <class vec4>
  TLV MakeTLV(vec4 const p) {
    return TLV(p.px(), p.py(), p.pz(), p.E());
  }

  template <class Vvec4>
  inline TLVs MakeTLVs(Vvec4 const &rivet_jets) {
    TLVs jets;
    for (auto jet : rivet_jets)
      jets.push_back(MakeTLV(jet));
    return jets;
  }

  // Structure holding information about the current event:
  // Four-momenta and event classification according to the
  // Higgs Template Cross Section
  struct HiggsClassification {
    // Higgs production mode
    HTXS::HiggsProdMode prodMode;
    // The Higgs boson
    TLV higgs;
    // The Higgs boson decay products
    TLV p4decay_higgs;
    // Associated vector bosons
    TLV V;
    // The V-boson decay products
    TLV p4decay_V;
    // Jets are built ignoring Higgs decay products and leptons from V decays
    // jets with pT > 25 GeV and 30 GeV
    TLVs jets25, jets30;
    // Event categorization according to YR4 wrtietup
    // https://cds.cern.ch/record/2138079
    HTXS::Stage0::Category stage0_cat;
    HTXS::Stage1::Category stage1_cat_pTjet25GeV;
    HTXS::Stage1::Category stage1_cat_pTjet30GeV;
    HTXS::Stage1_1::Category stage1_1_cat_pTjet25GeV;
    HTXS::Stage1_1::Category stage1_1_cat_pTjet30GeV;
    HTXS::Stage1_1_Fine::Category stage1_1_fine_cat_pTjet25GeV;
    HTXS::Stage1_1_Fine::Category stage1_1_fine_cat_pTjet30GeV;
    // Error code :: classification was succesful or some error occured
    HTXS::ErrorCode errorCode;
  };

  template <class category>
  inline HTXS::HiggsClassification Rivet2Root(category const &htxs_cat_rivet) {
    HTXS::HiggsClassification cat;
    cat.prodMode = htxs_cat_rivet.prodMode;
    cat.errorCode = htxs_cat_rivet.errorCode;
    cat.higgs = MakeTLV(htxs_cat_rivet.higgs);
    cat.V = MakeTLV(htxs_cat_rivet.V);
    cat.p4decay_higgs = MakeTLV(htxs_cat_rivet.p4decay_higgs);
    cat.p4decay_V = MakeTLV(htxs_cat_rivet.p4decay_V);
    cat.jets25 = MakeTLVs(htxs_cat_rivet.jets25);
    cat.jets30 = MakeTLVs(htxs_cat_rivet.jets30);
    cat.stage0_cat = htxs_cat_rivet.stage0_cat;
    cat.stage1_cat_pTjet25GeV = htxs_cat_rivet.stage1_cat_pTjet25GeV;
    cat.stage1_cat_pTjet30GeV = htxs_cat_rivet.stage1_cat_pTjet30GeV;
    cat.stage1_1_cat_pTjet25GeV = htxs_cat_rivet.stage1_1_cat_pTjet25GeV;
    cat.stage1_1_cat_pTjet30GeV = htxs_cat_rivet.stage1_1_cat_pTjet30GeV;
    cat.stage1_1_fine_cat_pTjet25GeV = htxs_cat_rivet.stage1_1_fine_cat_pTjet25GeV;
    cat.stage1_1_fine_cat_pTjet30GeV = htxs_cat_rivet.stage1_1_fine_cat_pTjet30GeV;
    return cat;
  }

  inline int HTXSstage1_to_HTXSstage1FineIndex(HTXS::Stage1::Category stage1, HiggsProdMode prodMode, tH_type tH) {
    if (stage1 == HTXS::Stage1::Category::UNKNOWN)
      return 0;
    int P = (int)(stage1 / 100);
    int F = (int)(stage1 % 100);
    // 1.a spit tH categories
    if (prodMode == HiggsProdMode::TH) {
      // check that tH splitting is valid for Stage-1 FineIndex
      // else return unknown category
      if (tH == tH_type::noTH)
        return 0;
      // check if forward tH
      int fwdH = F == 0 ? 0 : 1;
      return (49 + 2 * (tH - 1) + fwdH);
    }
    // 1.b QQ2HQQ --> split into VBF, WH, ZH -> HQQ
    // offset vector 1: input is the Higgs prodMode
    // first two indicies are dummies, given that this is only called for prodMode=2,3,4
    std::vector<int> pMode_offset = {0, 0, 13, 19, 25};
    if (P == 2)
      return (F + pMode_offset[prodMode]);
    // 1.c remaining categories
    // offset vector 2: input is the Stage-1 category P
    // third index is dummy, given that this is called for category P=0,1,3,4,5,6,7
    std::vector<int> catP_offset = {0, 1, 0, 31, 36, 41, 45, 47};
    return (F + catP_offset[P]);
  }

  inline int HTXSstage1_to_HTXSstage1FineIndex(const HiggsClassification &stxs,
                                               tH_type tH = noTH,
                                               bool jets_pT25 = false) {
    HTXS::Stage1::Category stage1 = jets_pT25 == false ? stxs.stage1_cat_pTjet30GeV : stxs.stage1_cat_pTjet25GeV;
    return HTXSstage1_to_HTXSstage1FineIndex(stage1, stxs.prodMode, tH);
  }

  inline int HTXSstage1_to_index(HTXS::Stage1::Category stage1) {
    // the Stage-1 categories
    int P = (int)(stage1 / 100);
    int F = (int)(stage1 % 100);
    std::vector<int> offset{0, 1, 13, 19, 24, 29, 33, 35, 37, 39};
    // convert to linear values
    return (F + offset[P]);
  }

  // #endif

}  // namespace HTXS

#ifdef RIVET_Particle_HH
//#ifdef HIGGSTRUTHCLASSIFIER_HIGGSTRUTHCLASSIFIER_CC
//#include "Rivet/Particle.hh"
namespace Rivet {

  struct HiggsClassification {
    HTXS::HiggsProdMode prodMode;
    Rivet::Particle higgs;
    Rivet::Particle V;
    Rivet::FourMomentum p4decay_higgs;
    Rivet::FourMomentum p4decay_V;
    Rivet::Jets jets25, jets30;
    HTXS::Stage0::Category stage0_cat;
    HTXS::Stage1::Category stage1_cat_pTjet25GeV;
    HTXS::Stage1::Category stage1_cat_pTjet30GeV;
    HTXS::Stage1_1::Category stage1_1_cat_pTjet25GeV;
    HTXS::Stage1_1::Category stage1_1_cat_pTjet30GeV;
    HTXS::Stage1_1_Fine::Category stage1_1_fine_cat_pTjet25GeV;
    HTXS::Stage1_1_Fine::Category stage1_1_fine_cat_pTjet30GeV;
    HTXS::ErrorCode errorCode;
  };
}  // namespace Rivet
#endif

#endif