| | 1 | |
| | 2 | |
| | 3 | == %$ h \rightarrow A^0A^0 \rightarrow b\overline{b} \tau^+\tau^-$% in VBF at LHC (Les Houches Project) == |
| | 4 | |
| | 5 | === Hypothesis === |
| | 6 | |
| | 7 | We assume SM-like weak vector boson fusion production of %$h$% with a xsec of order 2 to 4 pb depending on the mass. NLO corrections are available but have been shown to be small (typically of order 5 to 10% at most). We assume that %$Br(h\rightarrow A^0 A^0)=1$%. It is a simplifying hypothesis, often justified from the theoretical point of view (%$A^0$% is much heavier than light fermions). Anyway, if it's not true, all signal xsec should simply be scaled by this Br. |
| | 8 | |
| | 9 | We also assume that %$A^0$% is a Higgs like scalar coupling to the mass, so that |
| | 10 | |
| | 11 | |
| | 12 | |
| | 13 | Numerically, this gives %$Br(A^0\rightarrow \tau^+\tau^-)=8\%$% and %$Br(A^0\rightarrow b\overline{b})=92\%$% if one takes %$m_b=3.5$% Gev which is a reasonable value at the scale %$m_{A^0}$%. These BR could be modified if the coupling of %$A^0$% to charm is sizable, but this is not the case with SM like couplings or in large %$\tan(\beta)$% type II models like MSSM. |
| | 14 | |
| | 15 | The total BR for the decay chain %$ h \rightarrow A^0A^0 \rightarrow b\overline{b} \tau^+\tau^-$% is thus 2 * 8% * 92%=15%, and the typical xsec is thus ranging from 300 to 450 fb. |
| | 16 | |
| | 17 | === Benchmark point === |
| | 18 | |
| | 19 | For the LH project, we take %$m_{A^0}$% at 50 Gev and %$m_{h}$% at 120 Gev. The mass of the %$h$% is light enough to have a good production xsec and to avoid to consider its decays into WW or ZZ. The mass of the %$A^0$% is light enough to open the decay %$h\rightarrow A^0 A^0$% but heavy enough to give a good pT and a good separation for b's and tau's. |
| | 20 | |
| | 21 | === Monte Carlo production === |
| | 22 | |
| | 23 | ==== Methodology ==== |
| | 24 | |
| | 25 | Signal and backgrounds are generated separately using MG/ME v4.1. The PDF is CTEQ6L1 and the renormalization and factorization scales are set to 120 Gev. All the following generation cuts are applied. Jet means non b jets and leptons means e, %$\mu$% and %$\tau$%. |
| | 26 | |
| | 27 | * Minimal pT of 20 Gev for the jets and 10 Gev for leptons and b's |
| | 28 | * Maximum rapidity of 5 for jets and 2.5 for leptons and b's |
| | 29 | * Minimal deltaR>0.3 for all objects |
| | 30 | * Maximal invariant mass of 80 Gev for the tau's pair (to avoid the Z peak) |
| | 31 | * Minimal invariant mass of 700 Gev for the two jets (VBF cut) |
| | 32 | * Delta eta >4 for the two jets (VBF cut) |
| | 33 | |
| | 34 | ==== Signal ==== |
| | 35 | |
| | 36 | 100k events with tau's decayed into leptons only: [http://madgraph.phys.ucl.ac.be/MadGraphData/mherquet@fyma.ucl.ac.be/Decay/signal_out.lhe.gz Events] |
| | 37 | |
| | 38 | The cross section '''after tau decay''' into leptons is 9.5fb. There is '''no cut''' on the tau decay products. Plots before tau decay are [http://madgraph.hep.uiuc.edu/MadGraphData/mherquet@fyma.ucl.ac.be/PROC4/Events/run_06_plots.html here]. |
| | 39 | |
| | 40 | If you need to split LHEF files, please use Software.SplitEvents |
| | 41 | |
| | 42 | ==== Considered backgrounds ==== |
| | 43 | |
| | 44 | ===== Irreducible: QCD %$2\tau 2b 2j$% background ===== |
| | 45 | |
| | 46 | 44k events with tau's decayed into leptons only |
| | 47 | [http://madgraph.phys.ucl.ac.be/MadGraphData/mherquet@fyma.ucl.ac.be/Decay/bkg_2tau2b2j_out_out.lhe.gz Events] |
| | 48 | |
| | 49 | The cross section '''after tau decay''' into leptons is 1fb. There is '''no cut''' on the tau decay products. Plots before tau decay are [http://madgraph.phys.ucl.ac.be/MadGraphData/mherquet@fyma.ucl.ac.be/PROC1/Events/run_01_plots.html here]. |
| | 50 | |
| | 51 | ===== Nearly irreducible QCD %$2e 2b 2j$% and %$2mu 2b 2j$% background ===== |
| | 52 | |
| | 53 | Irreducible if no tau tag and no missing Et cut. Cross section is 8.7fb (for each) with cuts on leptons same as cuts on tau's for signal. |
| | 54 | |
| | 55 | A 25k events '''muon''' sample: [http://madgraph.phys.ucl.ac.be/MadGraphData/mherquet@fyma.ucl.ac.be/PROC7/Events/bkg_2mu2b2j_1.lhe.gz Events] |
| | 56 | |
| | 57 | Another 26k events '''muon''' sample |
| | 58 | [http://madgraph.phys.ucl.ac.be/MadGraphData/mherquet@fyma.ucl.ac.be/PROC7/Events/bkg_2mu2b2j_2.lhe.gz Events] |
| | 59 | |
| | 60 | I propose to use one of these two '''muon samples''' as an electron sample (just need to switch the PDG code in the LHEF file since the kinematic is the same). |
| | 61 | |
| | 62 | Some [http://madgraph.phys.ucl.ac.be/MadGraphData/mherquet@fyma.ucl.ac.be/PROC7/Events/run_02_plots.html Plots] for the whole 50k sample. |
| | 63 | |
| | 64 | ===== Reducible tt + 2j with t>Wb and both W going to e, mu or tau and tau going to e or mu ===== |
| | 65 | |
| | 66 | |
| | 67 | A 100k event sample with W's in e or mu (xsec=2300fb): [http://madgraph.phys.ucl.ac.be/MadGraphData/mherquet@fyma.ucl.ac.be/Decay/run_03_unweighted_events_out.lhe.gz Events] |
| | 68 | |
| | 69 | A 50k event sample with W+ in tau's (and tau+ in e or mu) and W- in e or mu (xsec=404fb): |
| | 70 | [http://madgraph.phys.ucl.ac.be/MadGraphData/mherquet@fyma.ucl.ac.be/Decay/run_04_1_unweighted_events_out.lhe.gz Events] |
| | 71 | |
| | 72 | A 50k event sample with W- in tau's (and tau- in e or mu) and W+ in e or mu (xsec=404fb): |
| | 73 | [http://madgraph.phys.ucl.ac.be/MadGraphData/mherquet@fyma.ucl.ac.be/Decay/run_04_2_unweighted_events_out.lhe.gz Events] |
| | 74 | |
| | 75 | A 50k event sample with W- and W+ in tau's (and tau's in e or mu) (xsec=70fb): [http://madgraph.phys.ucl.ac.be/MadGraphData/mherquet@fyma.ucl.ac.be/Decay/run_04_3_unweighted_events_out.lhe.gz Events], another one [http://madgraph.phys.ucl.ac.be/MadGraphData/mherquet@fyma.ucl.ac.be/Decay/run_04_4_unweighted_events_out.lhe.gz Events] |
| | 76 | |
| | 77 | Without any cut on leptonic final state, the total xsec is 3200 fb! But the kinematic is '''very''' different for both leptons pairs and b pairs, large invariant masses and large deltaR are typically preferred. Also, there is a lot (half of it if I'm right) of mismatched pairs e-mu (instead of e-e or mu-mu) which should be easy to remove. The amount of Et is also typically much higher (neutrinos are coming from W's), so a max Et cut should also be rather efficient. So there is maybe a hope... |
| | 78 | |
| | 79 | ==== Not considered backgrounds ==== |
| | 80 | |
| | 81 | ===== QED (i.e. alpha_EW^4) backgrounds ===== |
| | 82 | Should be small (see [http://arxiv.org/pdf/hep-ph/0702119 this paper]) and anyway impossible to simulate completely. |
| | 83 | |
| | 84 | ===== QCD %$2b 4j$% background for jets faking tau's ===== |
| | 85 | |
| | 86 | Not considered here since tau's are decaying into leptons. |
| | 87 | |
| | 88 | ===== 2ta4j without b's (mistagging) ===== |
| | 89 | |
| | 90 | Probably not necessary. CMS TDR I gives a b mistagging probability around %$10^{-2}$% for gluons and light quark jets (for a tagging efficiency of 0.5) and %$10^{-1}$% for c jets. The cross sections of 2ta4j and 2ta2c2j with the strong VBF cuts are probably to small to give sizable contribution after double mistag. |
| | 91 | |
| | 92 | ==== Summary of the situation... ==== |
| | 93 | A 10fb signal with a 1fb irreducible bkg, two 9fb nearly irreducible background which could be efficiently cut by a mass window on leptons and a delta R cut on b's. A huge 3.2 pb tt background with very different kinematic but same final state... |
| | 94 | |
| | 95 | === References === |
| | 96 | |
| | 97 | A good generic reference for VBF is the corresponding section in [http://arxiv.org/abs/hep-ph/0503172 Djouadi]. It gives different references for experimental analysis of %$h\rightarrow \tau^+\tau^-$% in VBF which could be considered as a good starting point (but ATLAS related!). |
| | 98 | |
| | 99 | ==== Les Houches 02 proceeding: ==== |
| | 100 | [http://arxiv.org/abs/hep-ph/0203056 Azuleos] |
| | 101 | |
| | 102 | ==== ATLAS notes ==== |
| | 103 | [http://doc.cern.ch//archive/electronic/cern/others/atlnot/Note/phys/phys-2002-018.pdf phys-2002-018] |
| | 104 | [http://doc.cern.ch//archive/electronic/cern/others/atlnot/Note/phys/phys-2003-004.pdf phys-2002-004] |
| | 105 | |
| | 106 | ==== Parton level Papers: ==== |
| | 107 | |
| | 108 | * [attachment:e014037.pdf PRD59 014037] |
| | 109 | |
| | 110 | * [attachment:e093005.pdf PRD61 093005] |
| | 111 | |
| | 112 | PLB454 297 |
| | 113 | |
| | 114 | ==== CMS ==== |
| | 115 | |
| | 116 | To be checked: Physics TDR and Notes |
| | 117 | |
| | 118 | === Why is it interesting ? === |
| | 119 | |
| | 120 | Theoretical point of view: if %$m_{A^0}>2 m_b$%, most of NMSSM signal with only tau's and mu's decays are suppressed. |
| | 121 | |
| | 122 | Pheno point of view: signal xsec similar to (or even better than) SM one for VBF with h in tau's which has been shown to be a quite interesting channel (4.5 sigmas in each tau's leptonic mode, 6 sigmas for the combined analysis after 30 fb-1). Jet veto cannot be applied straightforward but a "non b jet" veto is an interesting new concept (check if it is new!). The additional double b tag could maybe compensate the less effective background rejection without decreasing too much the signal (25%). |
| | 123 | |
| | 124 | === Final Version of our paper === |
| | 125 | |
| | 126 | * [attachment:VBF_hAA2b2ta.pdf pdf]: Our contribution to Les Houches Proceedings |
| | 127 | |
| | 128 | |
| | 129 | |
| | 130 | |
| | 131 | |
| | 132 | |
| | 133 | |
| | 134 | |
| | 135 | |
| | 136 | |
| | 137 | |
| | 138 | |
| | 139 | |
| | 140 | |
| | 141 | |
| | 142 | |
| | 143 | |
| | 144 | |
| | 145 | |
| | 146 | |
| | 147 | |
| | 148 | |
| | 149 | |
| | 150 | |
| | 151 | |
| | 152 | |
| | 153 | |
| | 154 | |
| | 155 | |
| | 156 | |
| | 157 | |
| | 158 | |
| | 159 | |
| | 160 | |
| | 161 | |
| | 162 | |
| | 163 | |