How to perform b-tagging #297
Description
Question
Dear madanalysis5 team,
For process generate p p > h+ h, (h+> W+ h, h > b b~), h > b b~ , I want to study its background, so I define three categories:
- 4b0j (four b-jets, no light jets)
- 3b1j (three b-jets, one light jet)
- 2b2j (two b-jets, two light jets) . I used MG5 + Pythia8 + Delphes, then I used the following script. ./bin/ma5 -R /home/cz/science/HEPtool_install/testmg5/run_bcats.ma5
import /home/cz/science/HEPtool_install/mg5amcnlo/myexample1/Events/run_01/tag_1_delphes_events.root as SIG
set main.stacking_method = normalize2one
（pb）
set main.lumi = 300
set SIG.xsection = 4.3764e-08
========== ==========
reject (j) PT < 20
reject (j) abs(ETA) > 2.5
reject (bjets) PT < 20
reject (bjets) abs(ETA) > 2.5
reject (e) PT < 10
reject (e) abs(ETA) > 2.5
reject (mu) PT < 10
reject (mu) abs(ETA) > 2.5
========== PC1 ==========
select N(l) < 2 or DR(l, l) > 0.4
select DR(j, j) > 0.4
select MET > 5
select N(j) == 4
========== ==========
define_region fourb0j
define_region threeb1j
define_region twob2j
select (fourb0j) N(bjets) == 4
select (threeb1j) N(bjets) == 3
select (twob2j) N(bjets) == 2
plot N(j)
plot N(bjets)
plot MET 40 0 500
plot THT 40 0 500
plot PT(j[1]) 40 0 500 [logY]
plot ETA(j[1]) 40 -10 10 [logY]
plot MT_MET(j[1]) 40 0 500 [logY]
plot PT(j[2]) 40 0 500 [logY]
plot ETA(j[2]) 40 -10 10 [logY]
plot MT_MET(j[2]) 40 0 500 [logY]
plot M(j[1] j[2]) 40 0 500 [logY]
plot DELTAR(j[1],j[2]) 40 0 10 [logY]
submit Hplus_PC1_bcats .
However, there are many errors that I cannot resolve. I’m using my own model, and in R-mode I can’t redefine particles, so I don’t know how to fix this and need help.