User Tools

Site Tools


2017:groups:tools:contur

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
2017:groups:tools:contur [2017/09/20 23:09]
jonathan.butterworth
2017:groups:tools:contur [2017/10/02 10:44] (current)
jonathan.butterworth
Line 4: Line 4:
   * [[https://​contur.hepforge.org/​|Contur at HepForge]]   * [[https://​contur.hepforge.org/​|Contur at HepForge]]
   * [[http://​contur.hepforge.org/​svn/​|Contur sources]] (CLTestSingle works. Many improvements in progress.)   * [[http://​contur.hepforge.org/​svn/​|Contur sources]] (CLTestSingle works. Many improvements in progress.)
-  * [[contur#2HDM Studies|2HDM Studies]]+  * [[2HDM Studies]]
   * [[2017:​groups:​higgs:​fxsnp:​susycascades|FXS:​ Higgs production from SUSY cascades]]   * [[2017:​groups:​higgs:​fxsnp:​susycascades|FXS:​ Higgs production from SUSY cascades]]
-  * [[contur#​Light Scalar Particles|Light Scalar Particles]] +  * [[Contur for Light Scalar Particles]] 
-  * [[Explanation of Contur plot format]]+  * [[http://​contur.hepforge.org/​share/​plot-format.html|Explanation of Contur plot format]]
 -------------------------------------------------------------------------------- --------------------------------------------------------------------------------
- 
- 
- 
-==== 2HDM Studies ==== 
- 
-//Kristin Lohwasser, Ken Lane, Jon Butterworth,​ David Grellscheid,​ Lukas Pritchett, ...// 
- 
-=== First look with default parameters === 
- 
-Made in Les Houches, Masses as below, Beta = 0.0813, MixH = 0.0883 
- 
-  * [[http://​www.hep.ucl.ac.uk/​~jmb/​Work/​2hdm/​wplots8/​ATLAS_2016_I1426515/​index.html|8 TeV WW plots]] 
- 
-=== Proposed parameter scan from Ken === 
- 
-Start with: 
- 
-  * m(h_1) = 125 GeV 
-  * m(h_2) = 30 GeV 
-  * m(h_3) = 30 GeV 
-  * beta = arcsin(v2/​246 GeV) = 0.05084 for v2 = 12.5 GeV (and = 0.08139 20 GeV, if you feel ambitious) 
-  * mixh = beta - alpha for alpha = -0.435 x 10^{-2} ​ (and -0.695 x 10{-2} for v2 = 20 GeV; ditto) 
- 
-The signal processes of interest are: 
- 
-$$\begin{aligned} 
-\bar{q} q \rightarrow \gamma, Z \rightarrow h^+ h^- \rightarrow \mu^+ \mu^- + E_T^{\rm missing} 
- \\ \bar{d} u \rightarrow W^+ \rightarrow h^+ h_2 \ {\rm AND}\ h^+ h_3 \rightarrow \mu^+ \mu^+ \mu^- + E_T^{\rm missing}, \\ {\rm (plus\ charge\ conjugate\ )} \\  \bar{q} q \rightarrow Z \rightarrow h_2 h_3 -> \bar{b} b \mu^+ \mu^-\end{aligned}$$ 
- 
-For these, you want to scan over $h^+$ mass from 95 GeV to 250 GeV. 
- 
-=== Update from Jon 18 Sept === 
- 
-Here's a [[http://​www.hep.ucl.ac.uk/​~jmb/​Work/​LH/​2HDM/​8_WEAK/​|first look]], for inclusive $h^+$ or $h^-$ production, 8 TeV, parameters as above and beta = 0.0508, mixh = 0.0547 and $m(h^{+/-}) = 100$ GeV. 
- 
-Total cross section 1.8 pb. 
- 
-Note that the ATLAS 4l result still has a lot of exclusion power, even though the $h_2 h_3$ production channel is switched off (with this switched on, the total cross section is ~346 pb). From looking at the Herwig logs, I think this contribution may come from events like $q\bar{q} \rightarrow h_2 h^+ \rightarrow h_2 h_3 \mu \nu$ with subsequent $h_2, h_3$ decays to muons. 
- 
-However, you can see that the [[http://​www.hep.ucl.ac.uk/​~jmb/​Work/​LH/​2HDM/​8_WEAK/​ATLAS_2016_I1426515/​index.html|ATLAS WW measurement]] also has some exclusion - about 92%. I don't understand why the [[http://​www.hep.ucl.ac.uk/​~jmb/​Work/​LH/​2HDM/​8_WEAK/​CMS_2017_I1467451/​index.html|CMS H->WW]] one doesn'​t have an exclusion anywhere in the same league. Either the Higgs selection cuts suppress it or there is a bug/feature in the rivet routine (either ATLAS or CMS!). I will look into it. 
- 
- 
--------------------------------------------------------------------------------- 
- 
-==== Light Scalar Particles ==== 
- 
-//Susan Shotkin Gascon-Shotkin,​ Linda Finco, Sijing Zhang, Grégory Moreau, Peter Richardson, Sylvain Fichet, David Grellscheid,​ Jon Butterworth,​ David Yallup...// 
- 
-Joint with [[2017:​groups:​higgs:​lsp|Higgs group]]. 
- 
-=== Herwig status === 
- 
-<​del>​Herwig not dealing with UFO file properly. PR aiming to fix the 3D vertex.</​del> ​ 
-This is now fixed. Peter has implemented the relevant vertices in Herwig/​ThePEG and its UFO interface, and validated the cross section against MadGraph. This is not yet in 7.1.1 but is in the Herwig 7 trunk in the 
-repository, hopefully in next release. The exception is the 5-point vertices which are used in the CP-even model: 
- 
- 
-   ​Warning:​ Lorentz structure VVVVS ( g g g g phiNP ) in V_49 is not supported. 
-   ​Warning:​ Lorentz structure VVSSS ( W- W+ H H phiNP ) in V_118 is not supported. 
-   ​Warning:​ Lorentz structure VVVVS ( a a W- W+ phiNP ) in V_122 is not supported. 
-   ​Warning:​ Lorentz structure VVVVS ( W- W- W+ W+ phiNP ) in V_126 is not supported. 
-   ​Warning:​ Lorentz structure VVVVS ( a W- W+ Z phiNP ) in V_150 is not supported. 
-   ​Warning:​ Lorentz structure VVSSS ( Z Z H H phiNP ) in V_156 is not supported. 
-   ​Warning:​ Lorentz structure VVVVS ( W- W+ Z Z phiNP ) in V_160 is not supported. 
- 
- 
-These are currently ignored - hope they aren't a major effect (they don't look like they should be...) 
- 
- 
-=== Some (preliminary!) results with Contur === 
- 
-Suggested parameter ranges from Sylvain: 
- 
- 
-   ​CP-even scalar: 
-   mphi in [10, 100] GeV 
-   f0B in [1, 10] TeV 
-   f0W in [1, 10] TeV 
-   f0H in [1, 10] TeV 
- 
-   ​CP-odd scalar: 
-   mphi in [10, 100] GeV 
-   f0B in [1, 10] TeV 
-   f0W in [1, 10] TeV 
- 
-   For each case, all the other f's should be set to very large values. ​ 
- 
-So, here are some scans using these ranges (f0 value labelled $\Lambda$ on the vertical axis, $M_\phi$ on the horizontal. ​ 
-[[2017:​groups:​tools:​colorbarkey.png]] 
- 
- 
- 
- 
-== Older Stuff == 
- 
-Not really sure of what sensible parameter ranges would be at the moment, but here is an example run for CP-Odd scalar with some 8 TeV data. I (Jon) set $m_\phi = 80$ GeV, and setting the mass scales for the effective couplings at 50 or 100 TeV as follows: $f0B, f0W, f0H, f0gam, f0Z = 50$ TeV, $f0G, f0u, f0d, f0l = 100$ TeV. Total cross section is $4.4 pb^{-1}$. 
- 
-The Contur "​exclusion"​ for these parameters with the 8 TeV measurements currently in Rivet is at 77% cl. The exclusion dominantly comes from the ATLAS diphoton cross sections: [[http://​inspirehep.net/​record/​1591327|arXiv:​1704.03839]]. 
-The comparison plots are here for [[http://​www.hep.ucl.ac.uk/​~jmb/​Work/​LH/​NSCPO/​7_HAD/​index.html|7 TeV]] [[http://​www.hep.ucl.ac.uk/​~jmb/​Work/​LH/​NSCPO/​8_HAD/​index.html|8 TeV]]. You can see that the most significant exclusion is from the mass peak, though many other differential cross sections show a significant effect (they are not all combined, since they are the same events - only the most significant contributes to the overall c.l.) 
- 
-I also made a run where I increased the couplings to weak bosons by setting $f0B, f0W, f0H, f0Z = 5$ TeV, just to see what happens. I also turned off the $\phi$+jet process, to focus on the electroweak-type signatures. The exclusion is now at 97% c.l., dominated by the ATLAS $W/Z + \gamma$ [[http://​inspirehep.net/​record/​1448301|measurement]]. The results are here for [[http://​www.hep.ucl.ac.uk/​~jmb/​Work/​LH/​NSCPO/​7_WEAK/​index.html|7 TeV]] and [[http://​www.hep.ucl.ac.uk/​~jmb/​Work/​LH/​NSCPO/​8_WEAK/​index.html|8 TeV]] (8 TeV is by far the most interesting). 
- 
-(Modulo BR changes, the exclusion from the diphotons would presumably still be similar to the above, since the production cross section is driven by $f0G$ which didn't change. This could be combined with the 97%, as the events are independent.) 
- 
-Next step? Agree some sensible parameter ranges to scan, presumably mass and some subset of couplings; I can combine 7, 8 and 13 TeV limits, of course. 
- 
- 
  
  
  
  
2017/groups/tools/contur.1505941794.txt.gz · Last modified: 2017/09/20 23:09 by jonathan.butterworth