HTXS Namespace Reference

Higgs Template Cross Section namespace. More...

Namespaces
	Stage0
	Namespace for Stage0 categorization.
	Stage1
	Stage1_1
	Stage1_1_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_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

typedef math::XYZTLorentzVectorD HTXS::TLV

Definition at line 266 of file HiggsTemplateCrossSections.h.

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

Definition at line 267 of file HiggsTemplateCrossSections.h.

Enumeration Type Documentation

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

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

enum HTXS::tH_type

Additional identifier flag for TH production modes.

Enumerator
noTH
THQB
TWH

Definition at line 32 of file HiggsTemplateCrossSections.h.

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

Function Documentation

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

inline

Definition at line 331 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]);
    }

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

inline

Definition at line 358 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);
    }

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

inline

Definition at line 366 of file HiggsTemplateCrossSections.h.

References F(), createfilelist::int, PFRecoTauDiscriminationByIsolation_cfi::offset, particlelevel_cff::stage1_1_cat_pTjet25GeV, particlelevel_cff::stage1_1_cat_pTjet30GeV, particlelevel_cff::stage1_1_fine_cat_pTjet25GeV, and particlelevel_cff::stage1_1_fine_cat_pTjet30GeV.

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

template<class vec4 >

TLV HTXS::MakeTLV

(

vec4 const

p

)

Definition at line 270 of file HiggsTemplateCrossSections.h.

Referenced by MakeTLVs(), and Rivet2Root().

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

template<class Vvec4 >

TLVs HTXS::MakeTLVs ( Vvec4 const &  rivet_jets )

inline

Definition at line 273 of file HiggsTemplateCrossSections.h.

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

Referenced by Rivet2Root().

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

template<class category >

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

inline

Definition at line 309 of file HiggsTemplateCrossSections.h.

References eostools::cat(), HTXS::HiggsClassification::errorCode, HTXS::HiggsClassification::higgs, HTXS::HiggsClassification::jets25, HTXS::HiggsClassification::jets30, MakeTLV(), MakeTLVs(), HTXS::HiggsClassification::p4decay_higgs, HTXS::HiggsClassification::p4decay_V, HTXS::HiggsClassification::prodMode, HTXS::HiggsClassification::stage0_cat, HTXS::HiggsClassification::stage1_1_cat_pTjet25GeV, HTXS::HiggsClassification::stage1_1_cat_pTjet30GeV, HTXS::HiggsClassification::stage1_1_fine_cat_pTjet25GeV, HTXS::HiggsClassification::stage1_1_fine_cat_pTjet30GeV, HTXS::HiggsClassification::stage1_cat_pTjet25GeV, HTXS::HiggsClassification::stage1_cat_pTjet30GeV, and HTXS::HiggsClassification::V.

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;
      return cat;    
    }