CMS 3D CMS Logo

Namespaces | Classes | Typedefs | Enumerations | Functions
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< TLVTLVs
 

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

Definition at line 495 of file HiggsTemplateCrossSections.h.

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

Definition at line 496 of file HiggsTemplateCrossSections.h.

Enumeration Type Documentation

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.

11  {
12  UNDEFINED = -99,
13  SUCCESS = 0,
14  PRODMODE_DEFINED = 1,
19  VH_IDENTIFICATION = 6,
22  };
failed to identify associated vector boson
successful classification
failed to identify Higgs boson decay products
failed momentum conservation
failed to identify associated vector boson decay products
failed to identify Higgs boson
failed to identify hard scatter vertex
production mode not defined
failed to identify top decay

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.

Additional identifier flag for TH production modes.

Enumerator
noTH 
THQB 
TWH 

Definition at line 28 of file HiggsTemplateCrossSections.h.

Function Documentation

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 671 of file HiggsTemplateCrossSections.h.

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

Referenced by HTXSstage1_2_Fine_to_HTXSstage1_2_Fine_FineIndex().

673  {
674  if (Stage1_2_Fine == HTXS::Stage1_2_Fine::Category::UNKNOWN)
675  return 0;
676  int P = (int)(Stage1_2_Fine / 100);
677  int F = (int)(Stage1_2_Fine % 100);
678  // 1.a spit tH categories
679  if (prodMode == HiggsProdMode::TH) {
680  // check that tH splitting is valid for Stage-1 FineIndex
681  // else return unknown category
682  if (tH == tH_type::noTH)
683  return 0;
684  // check if forward tH
685  int fwdH = F == 0 ? 0 : 1;
686  return (189 + 2 * (tH - 1) + fwdH);
687  }
688  // 1.b QQ2HQQ --> split into VBF, WH, ZH -> HQQ
689  // offset vector 1: input is the Higgs prodMode
690  // first two indicies are dummies, given that this is only called for prodMode=2,3,4
691  std::vector<int> pMode_offset = {0, 0, 57, 82, 107};
692  if (P == 2)
693  return (F + pMode_offset[prodMode]);
694  // 1.c GG2ZH split into gg->ZH-had and gg->ZH-lep
695  if (prodMode == HiggsProdMode::GG2ZH && P == 1)
696  return F + 29;
697  // 1.d remaining categories
698  // offset vector 2: input is the Stage-1 category P
699  // third index is dummy, given that this is called for category P=0,1,3,4,5,6,7
700  std::vector<int> catP_offset = {0, 1, 0, 132, 148, 164, 180, 187};
701  return (F + catP_offset[P]);
702  }
std::pair< OmniClusterRef, TrackingParticleRef > P
static uInt32 F(BLOWFISH_CTX *ctx, uInt32 x)
Definition: blowfish.cc:281
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 704 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.

706  {
707  HTXS::Stage1_2_Fine::Category Stage1_2_Fine =
708  jets_pT25 == false ? stxs.stage1_2_fine_cat_pTjet30GeV : stxs.stage1_2_fine_cat_pTjet25GeV;
709  return HTXSstage1_2_Fine_to_HTXSstage1_2_Fine_FineIndex(Stage1_2_Fine, stxs.prodMode, tH);
710  }
int HTXSstage1_2_Fine_to_HTXSstage1_2_Fine_FineIndex(const HiggsClassification &stxs, tH_type tH=noTH, bool jets_pT25=false)
int HTXS::HTXSstage1_2_Fine_to_index ( HTXS::Stage1_2_Fine::Category  Stage1_2_Fine)
inline

Definition at line 712 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, particlelevel_cff::stage1_1_fine_cat_pTjet30GeV, particlelevel_cff::stage1_2_cat_pTjet25GeV, particlelevel_cff::stage1_2_cat_pTjet30GeV, particlelevel_cff::stage1_2_fine_cat_pTjet25GeV, and particlelevel_cff::stage1_2_fine_cat_pTjet30GeV.

712  {
713  // the Stage-1_2_Fine categories
714  int P = (int)(Stage1_2_Fine / 100);
715  int F = (int)(Stage1_2_Fine % 100);
716  std::vector<int> offset{0, 1, 29, 54, 70, 86, 102, 109, 111, 113};
717  // convert to linear values
718  return (F + offset[P]);
719  }
std::pair< OmniClusterRef, TrackingParticleRef > P
static uInt32 F(BLOWFISH_CTX *ctx, uInt32 x)
Definition: blowfish.cc:281
int HTXS::HTXSstage1_2_to_HTXSstage1_2_FineIndex ( HTXS::Stage1_2::Category  stage1_2,
HiggsProdMode  prodMode,
tH_type  tH 
)
inline

Definition at line 619 of file HiggsTemplateCrossSections.h.

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

Referenced by HTXSstage1_2_to_HTXSstage1_2_FineIndex().

621  {
622  if (stage1_2 == HTXS::Stage1_2::Category::UNKNOWN)
623  return 0;
624  int P = (int)(stage1_2 / 100);
625  int F = (int)(stage1_2 % 100);
626  // 1.a spit tH categories
627  if (prodMode == HiggsProdMode::TH) {
628  // check that tH splitting is valid for Stage-1 FineIndex
629  // else return unknown category
630  if (tH == tH_type::noTH)
631  return 0;
632  // check if forward tH
633  int fwdH = F == 0 ? 0 : 1;
634  return (94 + 2 * (tH - 1) + fwdH);
635  }
636  // 1.b QQ2HQQ --> split into VBF, WH, ZH -> HQQ
637  // offset vector 1: input is the Higgs prodMode
638  // first two indicies are dummies, given that this is only called for prodMode=2,3,4
639  std::vector<int> pMode_offset = {0, 0, 35, 46, 57};
640  if (P == 2)
641  return (F + pMode_offset[prodMode]);
642  // 1.c GG2ZH split into gg->ZH-had and gg->ZH-lep
643  if (prodMode == HiggsProdMode::GG2ZH && P == 1)
644  return F + 18;
645  // 1.d remaining categories
646  // offset vector 2: input is the Stage-1 category P
647  // third index is dummy, given that this is called for category P=0,1,3,4,5,6,7
648  std::vector<int> catP_offset = {0, 1, 0, 68, 74, 80, 86, 92};
649  return (F + catP_offset[P]);
650  }
std::pair< OmniClusterRef, TrackingParticleRef > P
static uInt32 F(BLOWFISH_CTX *ctx, uInt32 x)
Definition: blowfish.cc:281
int HTXS::HTXSstage1_2_to_HTXSstage1_2_FineIndex ( const HiggsClassification stxs,
tH_type  tH = noTH,
bool  jets_pT25 = false 
)
inline

Definition at line 652 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.

654  {
655  HTXS::Stage1_2::Category stage1_2 =
656  jets_pT25 == false ? stxs.stage1_2_cat_pTjet30GeV : stxs.stage1_2_cat_pTjet25GeV;
657  return HTXSstage1_2_to_HTXSstage1_2_FineIndex(stage1_2, stxs.prodMode, tH);
658  }
int HTXSstage1_2_to_HTXSstage1_2_FineIndex(const HiggsClassification &stxs, tH_type tH=noTH, bool jets_pT25=false)
int HTXS::HTXSstage1_2_to_index ( HTXS::Stage1_2::Category  stage1_2)
inline

Definition at line 660 of file HiggsTemplateCrossSections.h.

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

660  {
661  // the Stage-1 categories
662  int P = (int)(stage1_2 / 100);
663  int F = (int)(stage1_2 % 100);
664  //std::vector<int> offset{0,1,13,19,24,29,33,35,37,39};
665  std::vector<int> offset{0, 1, 18, 29, 35, 41, 47, 53, 55, 57};
666  // convert to linear values
667  return (F + offset[P]);
668  }
std::pair< OmniClusterRef, TrackingParticleRef > P
static uInt32 F(BLOWFISH_CTX *ctx, uInt32 x)
Definition: blowfish.cc:281
int HTXS::HTXSstage1_to_HTXSstage1FineIndex ( HTXS::Stage1::Category  stage1,
HiggsProdMode  prodMode,
tH_type  tH 
)
inline

Definition at line 573 of file HiggsTemplateCrossSections.h.

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

Referenced by HTXSstage1_to_HTXSstage1FineIndex().

573  {
574  if (stage1 == HTXS::Stage1::Category::UNKNOWN)
575  return 0;
576  int P = (int)(stage1 / 100);
577  int F = (int)(stage1 % 100);
578  // 1.a spit tH categories
579  if (prodMode == HiggsProdMode::TH) {
580  // check that tH splitting is valid for Stage-1 FineIndex
581  // else return unknown category
582  if (tH == tH_type::noTH)
583  return 0;
584  // check if forward tH
585  int fwdH = F == 0 ? 0 : 1;
586  return (49 + 2 * (tH - 1) + fwdH);
587  }
588  // 1.b QQ2HQQ --> split into VBF, WH, ZH -> HQQ
589  // offset vector 1: input is the Higgs prodMode
590  // first two indicies are dummies, given that this is only called for prodMode=2,3,4
591  std::vector<int> pMode_offset = {0, 0, 13, 19, 25};
592  if (P == 2)
593  return (F + pMode_offset[prodMode]);
594  // 1.c remaining categories
595  // offset vector 2: input is the Stage-1 category P
596  // third index is dummy, given that this is called for category P=0,1,3,4,5,6,7
597  std::vector<int> catP_offset = {0, 1, 0, 31, 36, 41, 45, 47};
598  return (F + catP_offset[P]);
599  }
std::pair< OmniClusterRef, TrackingParticleRef > P
static uInt32 F(BLOWFISH_CTX *ctx, uInt32 x)
Definition: blowfish.cc:281
int HTXS::HTXSstage1_to_HTXSstage1FineIndex ( const HiggsClassification stxs,
tH_type  tH = noTH,
bool  jets_pT25 = false 
)
inline

Definition at line 601 of file HiggsTemplateCrossSections.h.

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

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

Definition at line 608 of file HiggsTemplateCrossSections.h.

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

608  {
609  // the Stage-1 categories
610  int P = (int)(stage1 / 100);
611  int F = (int)(stage1 % 100);
612  std::vector<int> offset{0, 1, 13, 19, 24, 29, 33, 35, 37, 39};
613  // convert to linear values
614  return (F + offset[P]);
615  }
std::pair< OmniClusterRef, TrackingParticleRef > P
static uInt32 F(BLOWFISH_CTX *ctx, uInt32 x)
Definition: blowfish.cc:281
template<class vec4 >
TLV HTXS::MakeTLV ( vec4 const  p)

Definition at line 499 of file HiggsTemplateCrossSections.h.

Referenced by MakeTLVs(), and Rivet2Root().

499  {
500  return TLV(p.px(), p.py(), p.pz(), p.E());
501  }
math::XYZTLorentzVectorD TLV
template<class Vvec4 >
TLVs HTXS::MakeTLVs ( Vvec4 const &  rivet_jets)
inline

Definition at line 504 of file HiggsTemplateCrossSections.h.

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

Referenced by Rivet2Root().

504  {
505  TLVs jets;
506  for (auto jet : rivet_jets)
507  jets.push_back(MakeTLV(jet));
508  return jets;
509  }
std::vector< TLV > TLVs
vector< PseudoJet > jets
TLV MakeTLV(vec4 const p)
template<class category >
HTXS::HiggsClassification HTXS::Rivet2Root ( category const &  htxs_cat_rivet)
inline

Definition at line 548 of file HiggsTemplateCrossSections.h.

References eostools::cat(), HTXS::HiggsClassification::errorCode, HTXS::HiggsClassification::higgs, HTXS::HiggsClassification::isZ2vvDecay, 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_2_cat_pTjet25GeV, HTXS::HiggsClassification::stage1_2_cat_pTjet30GeV, HTXS::HiggsClassification::stage1_2_fine_cat_pTjet25GeV, HTXS::HiggsClassification::stage1_2_fine_cat_pTjet30GeV, HTXS::HiggsClassification::stage1_cat_pTjet25GeV, HTXS::HiggsClassification::stage1_cat_pTjet30GeV, and HTXS::HiggsClassification::V.

Referenced by HTXSRivetProducer::produce().

548  {
550  cat.prodMode = htxs_cat_rivet.prodMode;
551  cat.errorCode = htxs_cat_rivet.errorCode;
552  cat.higgs = MakeTLV(htxs_cat_rivet.higgs);
553  cat.V = MakeTLV(htxs_cat_rivet.V);
554  cat.p4decay_higgs = MakeTLV(htxs_cat_rivet.p4decay_higgs);
555  cat.p4decay_V = MakeTLV(htxs_cat_rivet.p4decay_V);
556  cat.jets25 = MakeTLVs(htxs_cat_rivet.jets25);
557  cat.jets30 = MakeTLVs(htxs_cat_rivet.jets30);
558  cat.stage0_cat = htxs_cat_rivet.stage0_cat;
559  cat.stage1_cat_pTjet25GeV = htxs_cat_rivet.stage1_cat_pTjet25GeV;
560  cat.stage1_cat_pTjet30GeV = htxs_cat_rivet.stage1_cat_pTjet30GeV;
561  cat.stage1_1_cat_pTjet25GeV = htxs_cat_rivet.stage1_1_cat_pTjet25GeV;
562  cat.stage1_1_cat_pTjet30GeV = htxs_cat_rivet.stage1_1_cat_pTjet30GeV;
563  cat.stage1_1_fine_cat_pTjet25GeV = htxs_cat_rivet.stage1_1_fine_cat_pTjet25GeV;
564  cat.stage1_1_fine_cat_pTjet30GeV = htxs_cat_rivet.stage1_1_fine_cat_pTjet30GeV;
565  cat.stage1_2_cat_pTjet25GeV = htxs_cat_rivet.stage1_2_cat_pTjet25GeV;
566  cat.stage1_2_cat_pTjet30GeV = htxs_cat_rivet.stage1_2_cat_pTjet30GeV;
567  cat.stage1_2_fine_cat_pTjet25GeV = htxs_cat_rivet.stage1_2_fine_cat_pTjet25GeV;
568  cat.stage1_2_fine_cat_pTjet30GeV = htxs_cat_rivet.stage1_2_fine_cat_pTjet30GeV;
569  cat.isZ2vvDecay = htxs_cat_rivet.isZ2vvDecay;
570  return cat;
571  }
HTXS::Stage1_2_Fine::Category stage1_2_fine_cat_pTjet25GeV
HTXS::Stage1::Category stage1_cat_pTjet25GeV
HTXS::Stage1_2::Category stage1_2_cat_pTjet30GeV
HTXS::Stage1_1_Fine::Category stage1_1_fine_cat_pTjet30GeV
def cat(path)
Definition: eostools.py:401
HTXS::Stage1_1_Fine::Category stage1_1_fine_cat_pTjet25GeV
TLV MakeTLV(vec4 const p)
HTXS::Stage1::Category stage1_cat_pTjet30GeV
HTXS::Stage1_1::Category stage1_1_cat_pTjet25GeV
HTXS::Stage1_2::Category stage1_2_cat_pTjet25GeV
HTXS::Stage1_1::Category stage1_1_cat_pTjet30GeV
HTXS::Stage1_2_Fine::Category stage1_2_fine_cat_pTjet30GeV
TLVs MakeTLVs(Vvec4 const &rivet_jets)