HTXS Namespace Reference

Higgs Template Cross Section namespace. More...

Namespaces
	Stage0
	Namespace for Stage0 categorization.
	Stage1
	Stage1_1
	Stage1_1_Fine
	Stage1_2
	Stage1_2_Fine

Classes
struct	HiggsClassification

Typedefs
typedef math::XYZTLorentzVectorD	TLV
typedef std::vector< TLV >	TLVs

Enumerations
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 }
	Error code: whether the classification was successful or failed. More...
enum	HiggsProdMode {   UNKNOWN = 0, GGF = 1, VBF = 2, WH = 3,   QQ2ZH = 4, GG2ZH = 5, TTH = 6, BBH = 7,   TH = 8 }
	Higgs production modes, corresponding to input sample. More...
enum	tH_type { noTH = 0, THQB = 1, TWH = 2 }
	Additional identifier flag for TH production modes. More...

Functions
int	HTXSstage1_2_Fine_to_HTXSstage1_2_Fine_FineIndex (HTXS::Stage1_2_Fine::Category Stage1_2_Fine, HiggsProdMode prodMode, tH_type tH)
int	HTXSstage1_2_Fine_to_HTXSstage1_2_Fine_FineIndex (const HiggsClassification &stxs, tH_type tH=noTH, bool jets_pT25=false)
int	HTXSstage1_2_Fine_to_index (HTXS::Stage1_2_Fine::Category Stage1_2_Fine)
int	HTXSstage1_2_to_HTXSstage1_2_FineIndex (HTXS::Stage1_2::Category stage1_2, HiggsProdMode prodMode, tH_type tH)
int	HTXSstage1_2_to_HTXSstage1_2_FineIndex (const HiggsClassification &stxs, tH_type tH=noTH, bool jets_pT25=false)
int	HTXSstage1_2_to_index (HTXS::Stage1_2::Category stage1_2)
int	HTXSstage1_to_HTXSstage1FineIndex (HTXS::Stage1::Category stage1, HiggsProdMode prodMode, tH_type tH)
int	HTXSstage1_to_HTXSstage1FineIndex (const HiggsClassification &stxs, tH_type tH=noTH, bool jets_pT25=false)
int	HTXSstage1_to_index (HTXS::Stage1::Category stage1)
template<class vec4 >
TLV	MakeTLV (vec4 const p)
template<class Vvec4 >
TLVs	MakeTLVs (Vvec4 const &rivet_jets)
template<class category >
HTXS::HiggsClassification	Rivet2Root (category const &htxs_cat_rivet)

Detailed Description

Higgs Template Cross Section namespace.

Typedef Documentation

TLV

typedef math::XYZTLorentzVectorD HTXS::TLV

Definition at line 495 of file HiggsTemplateCrossSections.h.

TLVs

typedef std::vector<TLV> HTXS::TLVs

Definition at line 496 of file HiggsTemplateCrossSections.h.

Enumeration Type Documentation

ErrorCode

enum HTXS::ErrorCode

Error code: whether the classification was successful or failed.

Enumerator
UNDEFINED
SUCCESS	successful classification
PRODMODE_DEFINED	production mode not defined
MOMENTUM_CONSERVATION	failed momentum conservation
HIGGS_IDENTIFICATION	failed to identify Higgs boson
HIGGS_DECAY_IDENTIFICATION	failed to identify Higgs boson decay products
HS_VTX_IDENTIFICATION	failed to identify hard scatter vertex
VH_IDENTIFICATION	failed to identify associated vector boson
VH_DECAY_IDENTIFICATION	failed to identify associated vector boson decay products
TOP_W_IDENTIFICATION	failed to identify top decay

Definition at line 11 of file HiggsTemplateCrossSections.h.

                 {
    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         
  };

HiggsProdMode

enum HTXS::HiggsProdMode

Higgs production modes, corresponding to input sample.

Enumerator
UNKNOWN
GGF
VBF
WH
QQ2ZH
GG2ZH
TTH
BBH
TH

Definition at line 25 of file HiggsTemplateCrossSections.h.

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

tH_type

enum HTXS::tH_type

Additional identifier flag for TH production modes.

Enumerator
noTH
THQB
TWH

Definition at line 28 of file HiggsTemplateCrossSections.h.

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

Function Documentation

HTXSstage1_2_Fine_to_HTXSstage1_2_Fine_FineIndex() [1/2]

int HTXS::HTXSstage1_2_Fine_to_HTXSstage1_2_Fine_FineIndex ( HTXS::Stage1_2_Fine::Category  Stage1_2_Fine, HiggsProdMode  prodMode, tH_type  tH  )

inline

Definition at line 679 of file HiggsTemplateCrossSections.h.

References F(), GG2ZH, createfilelist::int, noTH, TH, and UNKNOWN.

Referenced by HTXSstage1_2_Fine_to_HTXSstage1_2_Fine_FineIndex().

                                                                          {
    if (Stage1_2_Fine == HTXS::Stage1_2_Fine::Category::UNKNOWN)
      return 0;
    int P = (int)(Stage1_2_Fine / 100);
    int F = (int)(Stage1_2_Fine % 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 (189 + 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, 57, 82, 107};
    if (P == 2)
      return (F + pMode_offset[prodMode]);
    // 1.c GG2ZH split into gg->ZH-had and gg->ZH-lep
    if (prodMode == HiggsProdMode::GG2ZH && P == 1)
      return F + 29;
    // 1.d 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, 132, 148, 164, 180, 187};
    return (F + catP_offset[P]);
  }

HTXSstage1_2_Fine_to_HTXSstage1_2_Fine_FineIndex() [2/2]

int HTXS::HTXSstage1_2_Fine_to_HTXSstage1_2_Fine_FineIndex ( const HiggsClassification &  stxs, tH_type  tH = noTH, bool  jets_pT25 = false  )

inline

Definition at line 712 of file HiggsTemplateCrossSections.h.

References HTXSstage1_2_Fine_to_HTXSstage1_2_Fine_FineIndex(), HTXS::HiggsClassification::prodMode, HTXS::HiggsClassification::stage1_2_fine_cat_pTjet25GeV, and HTXS::HiggsClassification::stage1_2_fine_cat_pTjet30GeV.

                                                                                      {
    HTXS::Stage1_2_Fine::Category Stage1_2_Fine =
        jets_pT25 == false ? stxs.stage1_2_fine_cat_pTjet30GeV : stxs.stage1_2_fine_cat_pTjet25GeV;
    return HTXSstage1_2_Fine_to_HTXSstage1_2_Fine_FineIndex(Stage1_2_Fine, stxs.prodMode, tH);
  }

HTXSstage1_2_Fine_to_index()

int HTXS::HTXSstage1_2_Fine_to_index ( HTXS::Stage1_2_Fine::Category  Stage1_2_Fine )

inline

Definition at line 720 of file HiggsTemplateCrossSections.h.

References F(), createfilelist::int, and HLT_IsoTrack_cff::offset.

                                                                                 {
    // the Stage-1_2_Fine categories
    int P = (int)(Stage1_2_Fine / 100);
    int F = (int)(Stage1_2_Fine % 100);
    std::vector<int> offset{0, 1, 29, 54, 70, 86, 102, 109, 111, 113};
    // convert to linear values
    return (F + offset[P]);
  }

HTXSstage1_2_to_HTXSstage1_2_FineIndex() [1/2]

int HTXS::HTXSstage1_2_to_HTXSstage1_2_FineIndex ( HTXS::Stage1_2::Category  stage1_2, HiggsProdMode  prodMode, tH_type  tH  )

inline

Definition at line 627 of file HiggsTemplateCrossSections.h.

References F(), GG2ZH, createfilelist::int, noTH, TH, and UNKNOWN.

Referenced by HTXSstage1_2_to_HTXSstage1_2_FineIndex().

                                                                {
    if (stage1_2 == HTXS::Stage1_2::Category::UNKNOWN)
      return 0;
    int P = (int)(stage1_2 / 100);
    int F = (int)(stage1_2 % 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 (94 + 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, 35, 46, 57};
    if (P == 2)
      return (F + pMode_offset[prodMode]);
    // 1.c GG2ZH split into gg->ZH-had and gg->ZH-lep
    if (prodMode == HiggsProdMode::GG2ZH && P == 1)
      return F + 18;
    // 1.d 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, 68, 74, 80, 86, 92};
    return (F + catP_offset[P]);
  }

HTXSstage1_2_to_HTXSstage1_2_FineIndex() [2/2]

int HTXS::HTXSstage1_2_to_HTXSstage1_2_FineIndex ( const HiggsClassification &  stxs, tH_type  tH = noTH, bool  jets_pT25 = false  )

inline

Definition at line 660 of file HiggsTemplateCrossSections.h.

References HTXSstage1_2_to_HTXSstage1_2_FineIndex(), HTXS::HiggsClassification::prodMode, HTXS::HiggsClassification::stage1_2_cat_pTjet25GeV, and HTXS::HiggsClassification::stage1_2_cat_pTjet30GeV.

                                                                            {
    HTXS::Stage1_2::Category stage1_2 =
        jets_pT25 == false ? stxs.stage1_2_cat_pTjet30GeV : stxs.stage1_2_cat_pTjet25GeV;
    return HTXSstage1_2_to_HTXSstage1_2_FineIndex(stage1_2, stxs.prodMode, tH);
  }

HTXSstage1_2_to_index()

int HTXS::HTXSstage1_2_to_index ( HTXS::Stage1_2::Category  stage1_2 )

inline

Definition at line 668 of file HiggsTemplateCrossSections.h.

References F(), createfilelist::int, and HLT_IsoTrack_cff::offset.

                                                                  {
    // the Stage-1 categories
    int P = (int)(stage1_2 / 100);
    int F = (int)(stage1_2 % 100);
    //std::vector<int> offset{0,1,13,19,24,29,33,35,37,39};
    std::vector<int> offset{0, 1, 18, 29, 35, 41, 47, 53, 55, 57};
    // convert to linear values
    return (F + offset[P]);
  }

HTXSstage1_to_HTXSstage1FineIndex() [1/2]

int HTXS::HTXSstage1_to_HTXSstage1FineIndex ( HTXS::Stage1::Category  stage1, HiggsProdMode  prodMode, tH_type  tH  )

inline

Definition at line 583 of file HiggsTemplateCrossSections.h.

References F(), createfilelist::int, noTH, TH, and UNKNOWN.

Referenced by HTXSstage1_to_HTXSstage1FineIndex().

                                                                                                              {
    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]);
  }

HTXSstage1_to_HTXSstage1FineIndex() [2/2]

int HTXS::HTXSstage1_to_HTXSstage1FineIndex ( const HiggsClassification &  stxs, tH_type  tH = noTH, bool  jets_pT25 = false  )

inline

Definition at line 611 of file HiggsTemplateCrossSections.h.

References HTXSstage1_to_HTXSstage1FineIndex(), HTXS::HiggsClassification::prodMode, HTXS::HiggsClassification::stage1_cat_pTjet25GeV, and HTXS::HiggsClassification::stage1_cat_pTjet30GeV.

                                                                       {
    HTXS::Stage1::Category stage1 = jets_pT25 == false ? stxs.stage1_cat_pTjet30GeV : stxs.stage1_cat_pTjet25GeV;
    return HTXSstage1_to_HTXSstage1FineIndex(stage1, stxs.prodMode, tH);
  }

HTXSstage1_to_index()

int HTXS::HTXSstage1_to_index ( HTXS::Stage1::Category  stage1 )

inline

Definition at line 618 of file HiggsTemplateCrossSections.h.

References F(), createfilelist::int, and HLT_IsoTrack_cff::offset.

                                                            {
    // 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]);
  }

MakeTLV()

template<class vec4 >

TLV HTXS::MakeTLV

(

vec4 const

p

)

Definition at line 499 of file HiggsTemplateCrossSections.h.

References AlCaHLTBitMon_ParallelJobs::p.

Referenced by MakeTLVs(), and Rivet2Root().

                            {
    return TLV(p.px(), p.py(), p.pz(), p.E());
  }

MakeTLVs()

template<class Vvec4 >

TLVs HTXS::MakeTLVs ( Vvec4 const &  rivet_jets )

inline

Definition at line 504 of file HiggsTemplateCrossSections.h.

References metsig::jet, PDWG_EXODelayedJetMET_cff::jets, and MakeTLV().

Referenced by Rivet2Root().

                                                {
    TLVs jets;
    for (const auto &jet : rivet_jets)
      jets.push_back(MakeTLV(jet));
    return jets;
  }

Rivet2Root()

template<class category >

HTXS::HiggsClassification HTXS::Rivet2Root ( category const &  htxs_cat_rivet )

inline

Definition at line 553 of file HiggsTemplateCrossSections.h.

References eostools::cat(), MakeTLV(), and MakeTLVs().

Referenced by HTXSRivetProducer::produce().

                                                                            {
    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;
    cat.stage1_2_cat_pTjet25GeV = htxs_cat_rivet.stage1_2_cat_pTjet25GeV;
    cat.stage1_2_cat_pTjet30GeV = htxs_cat_rivet.stage1_2_cat_pTjet30GeV;
    cat.stage1_2_fine_cat_pTjet25GeV = htxs_cat_rivet.stage1_2_fine_cat_pTjet25GeV;
    cat.stage1_2_fine_cat_pTjet30GeV = htxs_cat_rivet.stage1_2_fine_cat_pTjet30GeV;
    cat.isZ2vvDecay = htxs_cat_rivet.isZ2vvDecay;
    // Temporary fix: add variables to perform STXS 1.3 classification with nanoAOD on-the-fly
    cat.V_pt = htxs_cat_rivet.V_pt;
    cat.Mjj = htxs_cat_rivet.Mjj;
    cat.ptHjj = htxs_cat_rivet.ptHjj;
    cat.dPhijj = htxs_cat_rivet.dPhijj;
    return cat;
  }