Skip to content
Snippets Groups Projects
Commit e3946777 authored by Giulia Di Gregorio's avatar Giulia Di Gregorio Committed by Thomas Strebler
Browse files

Update the XS values and change loose definition

parent e1a9b061
No related branches found
Tags 0.14.0
No related merge requests found
......@@ -24,7 +24,8 @@ def btag_decor_cfg(flags):
]
gn2v01_valid_ptag = (
get_valid_ami_tag(split_tags, "p", "p6026") and not flags.Input.isPHYSLITE)
(get_valid_ami_tag(split_tags, "p", "p6026") and not flags.Input.isPHYSLITE)
or get_valid_ami_tag(split_tags, "p", "p6266"))
if gn2v01_valid_ptag:
btag_vars += [
"GN2v01_pb",
......
......@@ -139,7 +139,8 @@ def get_small_R_bjet_branches(
]
gn2v01_valid_ptag = (
get_valid_ami_tag(split_tags, "p", "p6026") and not flags.Input.isPHYSLITE)
(get_valid_ami_tag(split_tags, "p", "p6026") and not flags.Input.isPHYSLITE)
or get_valid_ami_tag(split_tags, "p", "p6266"))
if gn2v01_valid_ptag:
small_R_bjet_branches.variables += [
"GN2v01_pb",
......@@ -150,8 +151,8 @@ def get_small_R_bjet_branches(
if flags.Analysis.small_R_jet.doBtagExpVars:
small_R_bjet_branches.variables += [
"GN2v00_Db",
"GN2v00_pcbtExp"
"GN2v01_Db",
"GN2v01_pcbtExp"
]
return small_R_bjet_branches.get_output_list()
......@@ -46,10 +46,10 @@ namespace Easyjet
// loop over the floatVars_exp and intVars_exp and decorate the jet
if (m_doExp) {
SG::AuxElement::Decorator<float> DbExp("GN2v00_Db");
SG::AuxElement::Decorator<int> pcbtExp("GN2v00_pcbtExp");
SG::AuxElement::Decorator<float> DbExp("GN2v01_Db");
SG::AuxElement::Decorator<int> pcbtExp("GN2v01_pcbtExp");
// evaluate Db score of GN2v00
// evaluate Db score of GN2v01
float Db = evaluate_Db(*jet);
// assign an experimental pcbt score
......@@ -65,24 +65,23 @@ namespace Easyjet
}
float BTaggingDecoratorAlg ::evaluate_Db(const xAOD::Jet& jet) const {
static const SG::AuxElement::ConstAccessor<float> pu("GN2v00_pu");
static const SG::AuxElement::ConstAccessor<float> pc("GN2v00_pc");
static const SG::AuxElement::ConstAccessor<float> pb("GN2v00_pb");
float fc = 0.1; // HARDCODED fc value for GN2
return log( (pb(jet) ) / ( fc * pc(jet) + (1-fc) * pu(jet) ) );
static const SG::AuxElement::ConstAccessor<float> pu("GN2v01_pu");
static const SG::AuxElement::ConstAccessor<float> pc("GN2v01_pc");
static const SG::AuxElement::ConstAccessor<float> pb("GN2v01_pb");
static const SG::AuxElement::ConstAccessor<float> ptau("GN2v01_ptau");
float fc = 0.2; // HARDCODED fc value for GN2v01
float ftau = 0.01; // HARDCODED ftau value for GN2v01
return log( (pb(jet) ) / ( fc * pc(jet) + ftau*ptau(jet) + (1-fc-ftau) * pu(jet) ) );
}
int BTaggingDecoratorAlg ::evaluate_pcbtExp(float Db) const {
int PCBT = 0; // 100-89% -> 0
if (Db > 0.783) PCBT++; // 89-85% -> 1
if (Db > 1.638) PCBT++; // 85-82% -> 2
if (Db > 2.145) PCBT++; // 82-77% -> 3
if (Db > 2.893) PCBT++; // 77-75% -> 4
if (Db > 3.034) PCBT++; // 75-70% -> 5
if (Db > 3.875) PCBT++; // 70-68% -> 6
if (Db > 4.157) PCBT++; // 68-60% -> 7
if (Db > 5.394) PCBT++; // 60-0% -> 8
if (Db > -1.351) PCBT++; // 90-85% -> 1
if (Db > -0.396) PCBT++; // 85-77% -> 2
if (Db > 0.828) PCBT++; // 77-70% -> 3
if (Db > 1.877) PCBT++; // 70-65% -> 4
if (Db > 2.658) PCBT++; // 65-0% -> 5
return PCBT;
}
......
......@@ -2,11 +2,11 @@ ptag=p6266
campaign=v4
#data
easyjet-gridsubmit --data-list ../easyjet/bbyyAnalysis/datasets/PHYSLITE/nominal/data_13TeV.Run2.${ptag}.txt \
--run-config bbyyAnalysis/RunConfig-bbyy-skimming-loose.yaml \
--exec bbyy-ntupler \
--nGBperJob 2 \
--campaign ${campaign}
##easyjet-gridsubmit --data-list ../easyjet/bbyyAnalysis/datasets/PHYSLITE/nominal/data_13TeV.Run2.${ptag}.txt \
## --run-config bbyyAnalysis/RunConfig-bbyy-skimming-loose.yaml \
## --exec bbyy-ntupler \
## --nGBperJob 2 \
## --campaign ${campaign}
#ggF HH(bbyy) SM
easyjet-gridsubmit --mc-list ../easyjet/bbyyAnalysis/datasets/PHYSLITE/nominal/mc20_13TeV.ggFHH_bbyy_SM.${ptag}.txt \
......
......@@ -2,11 +2,11 @@ ptag=p6266
campaign=v4
#data
easyjet-gridsubmit --data-list ../easyjet/bbyyAnalysis/datasets/PHYSLITE/nominal/data_13p6TeV.Run3.p6029.txt \
--run-config bbyyAnalysis/RunConfig-bbyy-skimming-loose.yaml \
--exec bbyy-ntupler \
--nGBperJob 2 \
--campaign ${campaign}
#easyjet-gridsubmit --data-list ../easyjet/bbyyAnalysis/datasets/PHYSLITE/nominal/data_13p6TeV.Run3.${ptag}.txt \
# --run-config bbyyAnalysis/RunConfig-bbyy-skimming-loose.yaml \
# --exec bbyy-ntupler \
# --nGBperJob 2 \
# --campaign ${campaign}
#ggF HH(bbyy) SM
easyjet-gridsubmit --mc-list ../easyjet/bbyyAnalysis/datasets/PHYSLITE/nominal/mc23_13p6TeV.ggFHH_bbyy_SM.${ptag}.txt \
......
dataset_number/I:physics_short/C:crossSection_pb/D:kFactor/D:genFiltEff/D:relUncertUP/D:relUncertDOWN/D
600021 PhPy8EG_PDF4LHC15_HHbbyy_cHHH01d0 8.20992E-05 1.0 1.0 0.0 0.0
600022 PhPy8EG_PDF4LHC15_HHbbyy_cHHH10d0 1.71073E-03 1.0 1.0 0.0 0.0
603559 PhPy8EG_PDF4LHC21_HHbbyy_chhh1p0 8.20992E-05 1.0 1.0 0.0 0.0
600022 PhPy8EG_PDF4LHC15_HHbbyy_cHHH10d0 1.77736E-03 1.0 1.0 0.0 0.0
603696 PhPy8EG_PDF4LHC21_HHbbyy_chhh10p0 8.20992E-05 1.0 1.0 0.0 0.0
603558 PhPy8EG_PDF4LHC21_HHbbyy_chhh0p0 1.82892E-04 1.0 1.0 0.0 0.0
603559 PhPy8EG_PDF4LHC21_HHbbyy_chhh1p0 8.20992E-05 1.0 1.0 0.0 0.0
603560 PhPy8EG_PDF4LHC21_HHbbyy_chhh5p0 2.42279E-04 1.0 1.0 0.0 0.0
603695 PhPy8EG_PDF4LHC21_HHbbyy_chhh2p5 3.64588E-05 1.0 1.0 0.0 0.0
603696 PhPy8EG_PDF4LHC21_HHbbyy_chhh10p0 1.71073E-03 1.0 1.0 0.0 0.0
603697 PhPy8EG_PDF4LHC21_HHbbyy_chhhm1p0 3.40031E-04 1.0 1.0 0.0 0.0
503004 MGPy8EG_hh_bbyy_vbf_novhh_l1cvv1cv1 4.46059E-06 1.0 1.0 0.0 0.0
503007 MGPy8EG_hh_bbyy_vbf_novhh_l1cvv1p5cv1 8.87359E-06 1.0 1.0 0.0 0.0
503011 MGPy8EG_hh_bbyy_vbf_novhh_l2cvv1cv1 4.12215E-06 1.0 1.0 0.0 0.0
......@@ -10,18 +15,25 @@ dataset_number/I:physics_short/C:crossSection_pb/D:kFactor/D:genFiltEff/D:relUnc
503013 MGPy8EG_hh_bbyy_vbf_novhh_l1cvv1cv0p5 2.89740E-05 1.0 1.0 0.0 0.0
508689 MGPy8EG_hh_bbyy_vbf_novhh_l5mcvv1cv0p5 1.71126E-05 1.0 1.0 0.0 0.0
525376 MGPy8EG_hh_bbyy_vbf_novhh_l1cvv1cv1_5fns 4.46059E-06 1.0 1.0 0.0 0.0
525382 MGPy8EG_hh_bbyy_vbf_novhh_l1cvv1p5cv1_5fns 8.87359E-06 1.0 1.0 0.0 0.0
525377 MGPy8EG_hh_bbyy_vbf_novhh_l0cvv1cv1_5fns 1.21840E-05 1.0 1.0 0.0 0.0
525378 MGPy8EG_hh_bbyy_vbf_novhh_l2cvv1cv1_5fns 4.12215E-06 1.0 1.0 0.0 0.0
525379 MGPy8EG_hh_bbyy_vbf_novhh_l10cvv1cv1_5fns 2.59277E-04 1.0 1.0 0.0 0.0
525380 MGPy8EG_hh_bbyy_vbf_novhh_l1cvv0cv1_5fns 7.15730E-05 1.0 1.0 0.0 0.0
525381 MGPy8EG_hh_bbyy_vbf_novhh_l1cvv0p5cv1_5fns 2.54732E-05 1.0 1.0 0.0 0.0
525382 MGPy8EG_hh_bbyy_vbf_novhh_l1cvv1p5cv1_5fns 8.87359E-06 1.0 1.0 0.0 0.0
525383 MGPy8EG_hh_bbyy_vbf_novhh_l1cvv2cv1_5fns 3.83791E-05 1.0 1.0 0.0 0.0
525384 MGPy8EG_hh_bbyy_vbf_novhh_l1cvv3cv1_5fns 1.73006E-04 1.0 1.0 0.0 0.0
525385 MGPy8EG_hh_bbyy_vbf_novhh_l1cvv1cv0p5_5fns 2.89740E-05 1.0 1.0 0.0 0.0
525388 MGPy8EG_hh_bbyy_vbf_novhh_l5mcvv1cv0p5_5fns 1.71126E-05 1.0 1.0 0.0 0.0
525386 MGPy8EG_hh_bbyy_vbf_novhh_l1cvv1cv1p5_5fns 1.74698E-04 1.0 1.0 0.0 0.0
525387 MGPy8EG_hh_bbyy_vbf_novhh_l0cvv0cv1_5fns 9.77390E-05 1.0 1.0 0.0 0.0
525388 MGPy8EG_hh_bbyy_vbf_novhh_lm5cvv1cv0p5_5fns 1.71126E-05 1.0 1.0 0.0 0.0
343981 PowhegPythia8EvtGen_NNLOPS_nnlo_30_ggH125_gamgam 1.10277E-01 1.0 1.0 0.0 0.0
346214 PowhegPy8EG_NNPDF30_AZNLOCTEQ6L1_VBFH125_gamgam 8.58514E-03 1.0 1.0 0.0 0.0
345318 PowhegPythia8EvtGen_NNPDF30_AZNLO_WpH125J_Hyy_Wincl_MINLO 1.90680E-03 1.0 1.0 0.0 0.0
345317 PowhegPythia8EvtGen_NNPDF30_AZNLO_WmH125J_Hyy_Wincl_MINLO 1.20946E-03 1.0 1.0 0.0 0.0
345319 PowhegPythia8EvtGen_NNPDF30_AZNLO_ZH125J_Hyy_Zincl_MINLO 1.72792E-03 1.0 1.0 0.0 0.0
345061 PowhegPythia8EvtGen_NNPDF3_AZNLO_ggZH125_HgamgamZinc 2.78529E-04 1.0 1.0 0.0 0.0
345315 PowhegPythia8EvtGen_A14NNPDF23LO_bbHyy 0.00110390 1.0 1.0 0.0 0.0
345315 PowhegPythia8EvtGen_A14NNPDF23LO_bbHyy 1.10776E-03 1.0 1.0 0.0 0.0
346525 PowhegPythia8EvtGen_A14NNPDF23_NNPDF30ME_ttH125_gamgam 1.15112E-03 1.0 1.0 0.0 0.0
346677 aMcAtNloPythia8EvtGen_tHjb125_4fl_CPalpha_0_gamgam 1.68548E-04 1.0 1.0 0.0 0.0
346759 aMcAtNloPythia8EvtGen_tWH125_CPalpha_0_gamgam 3.44359E-05 1.0 1.0 0.0 0.0
......
dataset_number/I:physics_short/C:crossSection_pb/D:kFactor/D:genFiltEff/D:relUncertUP/D:relUncertDOWN/D
603558 PhPy8EG_PDF4LHC21_HHbbyy_chhh0p0 2.0208826E-04 1.0 1.0 0.0 0.0
603559 PhPy8EG_PDF4LHC21_HHbbyy_chhh1p0 9.1036225E-05 1.0 1.0 0.0 0.0
603560 PhPy8EG_PDF4LHC21_HHbbyy_chhh5p0 2.6573165E-04 1.0 1.0 0.0 0.0
603695 PhPy8EG_PDF4LHC21_HHbbyy_chhh2p5 4.0353135E-05 1.0 1.0 0.0 0.0
603696 PhPy8EG_PDF4LHC21_HHbbyy_chhh10p0 1.8783658E-03 1.0 1.0 0.0 0.0
603697 PhPy8EG_PDF4LHC21_HHbbyy_chhhm1p0 3.7519722E-04 1.0 1.0 0.0 0.0
525376 MGPy8EG_hh_bbyy_vbf_novhh_l1cvv1cv1_5fns 4.9550359E-06 1.0 1.0 0.0 0.0
602421 PhPy8EG_PDF4LHC21_ggH_NNLOPS_gammagamma 0.11857118 1.0 1.0 0.0 0.0
601482 PhPy8EG_PDF4LHC21_VBFH125_gammagamma 0.00925706 1.0 1.0 0.0 0.0
601484 PhPy8EG_PDF4LHC21_WpH125J_Wincl_MINLO_gammagamma 0.002017803 1.0 1.0 0.0 0.0
601483 PhPy8EG_PDF4LHC21_WmH125J_Wincl_MINLO_gammagamma 0.001288679 1.0 1.0 0.0 0.0
601523 PhPy8EG_PDF4LHC21_ZH125J_Zincl_MINLO_gammagamma 0.002142653 1.0 1.0 0.0 0.0
601522 PhPy8EG_PDF4LHC21_ggZH125_Zincl_gammagamma 0.00030872 1.0 1.0 0.0 0.0
601481 PhPy8EG_PDF4LHC21_ttH125_tincl_gammagamma 0.0012939 1.0 1.0 0.0 0.0
602422 PhPy8EG_PDF4LHC21_ttH125_tincl_gammagamma_1file 0.0012939 1.0 1.0 0.0 0.0
601710 PhPy8EG_PDF4LHC21_4FS_A14NNPDF_bbH_gammagamma 0.00119538 1.0 1.0 0.0 0.0
602421 PhPy8EG_PDF4LHC21_ggH_NNLOPS_gammagamma 0.1185621 1.0 1.0 0.0 0.0
601482 PhPy8EG_PDF4LHC21_VBFH125_gammagamma 9.25706E-03 1.0 1.0 0.0 0.0
601484 PhPy8EG_PDF4LHC21_WpH125J_Wincl_MINLO_gammagamma 2.017803E-03 1.0 1.0 0.0 0.0
601483 PhPy8EG_PDF4LHC21_WmH125J_Wincl_MINLO_gammagamma 1.288679E-03 1.0 1.0 0.0 0.0
601523 PhPy8EG_PDF4LHC21_ZH125J_Zincl_MINLO_gammagamma 1.833933E-03 1.0 1.0 0.0 0.0
601522 PhPy8EG_PDF4LHC21_ggZH125_Zincl_gammagamma 3.0872E-04 1.0 1.0 0.0 0.0
601481 PhPy8EG_PDF4LHC21_ttH125_tincl_gammagamma 1.2939E-03 1.0 1.0 0.0 0.0
602422 PhPy8EG_PDF4LHC21_ttH125_tincl_gammagamma_1file 1.2939E-03 1.0 1.0 0.0 0.0
601710 PhPy8EG_PDF4LHC21_4FS_A14NNPDF_bbH_gammagamma 1.19538E-03 1.0 1.0 0.0 0.0
800904 Py8EG_A14NNPDF23LO_XHS_X170_S030_HyySbb 1.0E-03 1.0 1.0 0.0 0.0
800913 Py8EG_A14NNPDF23LO_XHS_X245_S090_HyySbb 1.0E-03 1.0 1.0 0.0 0.0
800915 Py8EG_A14NNPDF23LO_XHS_X205_S030_HyySbb 1.0E-03 1.0 1.0 0.0 0.0
......
......@@ -11,7 +11,7 @@ do_muons: true
do_electrons: true
do_photons: true
do_met: true
do_KinematicFit: false
do_KinematicFit: true
# Toggles for reconstructed objects decorations
small_R_jet:
......
......@@ -7,12 +7,38 @@ bypass: false
selection_name: "loose"
small_R_jet:
doBtagExpVars: true
doBtagExpVars: false
##temporary lines for the resonant analysis
do_resonant_PNN: true
#Low Mass Grid
mX_mS_pairs: [[170, 30], [205, 30], [245, 90]]
#High Mass Grid
mS_values: [30, 170, 400]
mX_values: [750]
#1 bjet Mass Grid
mX_1bjet: [750]
###########
CutList :
- PASS_TRIGGER
- TWO_LOOSE_PHOTONS
- PASS_TRIGGER_MATCHING
- TWO_TIGHTID_ISO_PHOTONS
- DIPHOTON_MASS
- AT_LEAST_TWO_JETS
- PASS_TRIGGER
- TWO_LOOSE_PHOTONS
- PASS_TRIGGER_MATCHING
- TWO_TIGHTID_ISO_PHOTONS
- PASS_RELPT
- DIPHOTON_MASS
- EXACTLY_ZERO_LEPTONS
- AT_LEAST_TWO_JETS
Electron:
forceFullSimConfig: true
Photon:
forceFullSimConfig: true
ttree_output:
- AnalysisMiniTree:
include: bbyy-base-output.yaml
......@@ -7,6 +7,7 @@
collection_options:
small_R_jets:
btag_details: true
btag_info: true
photons:
shower_shapes: false
truth_labels: false
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment