| | 1 | |
| | 2 | |
| | 3 | == Jet Radiation in SUSY events == |
| | 4 | |
| | 5 | |
| | 6 | |
| | 7 | === People === |
| | 8 | 1. Johan Alwall |
| | 9 | 1. Simon de Visscher |
| | 10 | 1. Fabio Maltoni |
| | 11 | |
| | 12 | === Aim of the project === |
| | 13 | We study QCD radiation in SUSY events, such as gluino-gluino, stop-stop, stop-gluino processes. |
| | 14 | We use the MLM matching and verify/improve the results of |
| | 15 | [http://arxiv.org/abs/hep-ph/0510144 Tilman, David and Peter], and in particular we investigate |
| | 16 | whether the extra hard jets in SUSY events are well described by the parton-shower, as it is usually claimed, |
| | 17 | or if matching gives a better description. We also study the matching in heavy tops events as a reference and |
| | 18 | the SM background needed in these kind of analysis. The SM samples will be used also for other applications. |
| | 19 | |
| | 20 | === Tools === |
| | 21 | * Madgraph/MadEvent v4 : SUSY and matching |
| | 22 | * [http://cp3wks05.fynu.ucl.ac.be/twiki/bin/view/Software/MadMatchChecker MadMatchChecker] |
| | 23 | |
| | 24 | === Work === |
| | 25 | We have identified two phases: |
| | 26 | 1. MC samples generation and validation |
| | 27 | 2. Physics results and paper writeup |
| | 28 | |
| | 29 | ==== MC samples generation and validation ==== |
| | 30 | The strategy here is to first generate samples at the after-showering level to perform the validation. |
| | 31 | Once all samples are validated, we'll do the full Pythia-PGS simulation. |
| | 32 | |
| | 33 | ===== Generation ===== |
| | 34 | ====== Standard Model Samples (Simon) ====== |
| | 35 | To get a faster generation we: |
| | 36 | * define a new model (smzerobmass), where the bmass is set to "Zero" in the particles.dat. This optimizes the flavor sum |
| | 37 | * We don't decay the Z, W+,W-, and tops. We'll let these particles be decayed by Pythia. In doing this we loose spin correlations, which, however are irrelevant for the multi-jet samples for this type of studies. |
| | 38 | * Not to generate the single particle samples for Z and W. These have a negligible impact on the multi-jet |
| | 39 | configurations we are interested in. |
| | 40 | * In Pythia, we consider only decays of W into leptons (e, mu and tau) and also into neutrinos for Z. |
| | 41 | |
| | 42 | ====== SUSY and SM Samples ====== |
| | 43 | |
| | 44 | ||'''Process'''||'''Webpage'''||'''section(s)'''||'''LHE (on disks)'''||'''Pythia'''||'''MC'''||'''PGS'''|| |
| | 45 | ||1. '''gogo at 607 GeV (sps1a)'''|| n/a|| all ||ok || ok||ok || || |
| | 46 | ||1a. '''SC 1'''||n/a||4||ok||ok|||| n/a || |
| | 47 | ||1b. '''SC 2'''||n/a||4||ok|||||| n/a || |
| | 48 | ||1c. '''SC 3'''||n/a||4||ok|||||| n/a || |
| | 49 | ||1d. '''SC 4'''||n/a||4||ok|||||| n/a || |
| | 50 | ||1e. '''QCUT x 2'''||n/a||3.1||ok|||||| n/a || |
| | 51 | ||2. '''gogo at 607 GeV with xqcut and QCUT / 2'''|| [http://madgraph.phys.ucl.ac.be/MadGraphData/simon.devisscher@uclouvain.be/PROC40/HTML/crossx.html here] || 3.1 || || || || n/a || |
| | 52 | ||3. '''gogo at 607 GeV with scales divided by 2'''||[http://madgraph.phys.ucl.ac.be/MadGraphData/simon.devisscher@uclouvain.be/PROC40/HTML/crossx.html here] ||3 ||ok || || || n/a || |
| | 53 | ||4. '''gogo at 607 GeV with scales multiplied by 2'''|| [http://madgraph.phys.ucl.ac.be/MadGraphData/simon.devisscher@uclouvain.be/PROC40/HTML/crossx.html here]||3 || ok|| || || n/a || |
| | 54 | ||5. '''gogo at 300 GeV'''|| ||4 ||ok || ok|| ok|| n/a || |
| | 55 | ||6. '''gogo at 1200 GeV'''|| ||4 || ok||ok || ok|| n/a || |
| | 56 | ||7. '''sqsq at sps1a'''||n/a || 4,5|| ok|| || || || |
| | 57 | ||8. '''sqgo at sps1a'''|| n/a|| 4,5||ok|| || || || |
| | 58 | ||9. '''ulul at 607 GeV'''||[http://madgraph.hep.uiuc.edu/MadGraphData/devisscher@fynu.ucl.ac.be/PROC30/HTML/crossx.html here] ||3,4 ||ok ||ok || || n/a || |
| | 59 | ||10. '''ulul at 607 GeV with scales divided by 2'''|| [http://madgraph.hep.uiuc.edu/MadGraphData/johan@alwall.net/PROC60/HTML/crossx.html here]||3 ||ok ||ok ||ok || n/a || |
| | 60 | ||11. '''ulul at 607 GeV with scales multiplied by 2'''|| [http://madgraph.hep.uiuc.edu/MadGraphData/johan@alwall.net/PROC60/HTML/crossx.html here] ||3 ||ok ||ok || ok|| n/a || |
| | 61 | ||12. '''ululbar at 607 GeV'''|| [http://madgraph.hep.uiuc.edu/MadGraphData/devisscher@fynu.ucl.ac.be/PROC31/HTML/crossx.html here] ||3,4 ||ok ||ok ||ok || n/a || |
| | 62 | ||13. '''ululbar at 607 GeV with scales divided by 2'''||[http://madgraph.hep.uiuc.edu/MadGraphData/devisscher@fynu.ucl.ac.be/PROC37/HTML/crossx.html here] ||3 ||ok ||ok ||ok || n/a || |
| | 63 | ||14. '''ululbar at 607 GeV with scales multiplied by 2'''||[http://madgraph.hep.uiuc.edu/MadGraphData/devisscher@fynu.ucl.ac.be/PROC37/HTML/crossx.html here] ||3 || ok|| ok|| ok|| n/a || |
| | 64 | ||15. '''ttbar'''|| n/a|| 3,5|| ok||ok || ok|| || |
| | 65 | ||15a. '''QCUT x 2'''|| n/a|| 3,5|| ok|| || || n/a || |
| | 66 | ||16. '''ttbar with xqcut and QCUT / 2'''|| n/a|| 3.1 || || || || n/a || |
| | 67 | ||17. '''ttbar with scales divided by 2'''|| [http://madgraph.phys.ucl.ac.be/MadGraphData/johan@alwall.net/PROC169/HTML/crossx.html here]||3 ||ok||ok || ok|| n/a || |
| | 68 | ||18. '''ttbar with scales multiplied by 2'''||[http://madgraph.phys.ucl.ac.be/MadGraphData/johan@alwall.net/PROC169/HTML/crossx.html here] ||3 ||ok ||ok ||ok || n/a || |
| | 69 | ||19. '''W+4 jets biaised'''|| n/a||5 ||ok || || || || |
| | 70 | ||20. '''W+3 jets biaised'''|| [http://madgraph.phys.ucl.ac.be/MadGraphData/simon.devisscher@uclouvain.be/PROC37/HTML/crossx.html here]|| 5||ok|| || || || |
| | 71 | ||21. '''Z+4 jets biaised'''|| n/a||5 ||ok|| || || || |
| | 72 | |
| | 73 | |
| | 74 | ==== Physics results and paper writeup ==== |
| | 75 | |
| | 76 | ===== Plots for the paper ===== |
| | 77 | |
| | 78 | Explanation of notation: |
| | 79 | * 4x2 runs means 4 unmatched, 4 matched runs with the different shower params |
| | 80 | * 4 params means with the 4 different shower params |
| | 81 | * decay jet #s means number of jets from decays |
| | 82 | * jet numbers means plot of number of jets above (say) 50 GeV in the events |
| | 83 | |
| | 84 | ||'''Figure'''||'''section'''||'''Process(es)'''||'''Description'''||'''Done (plan)'''||'''plot'''|| |
| | 85 | || 1 || 2 || n/a || Feynman diagrams illustrating resonance double counting (3 plots) ||done|||| |
| | 86 | || 2 || 3.1 || 1,1e,2,15,15a,16 || Standard validation plots: diff. jet rates 0,1,2 with QCUT, x2, /2 for pT- and Q2-ordered showers, for gogo and ttbar (4x3x3 plots)|| ? |||| |
| | 87 | || 3a || 3.2 || 15 (4x2 runs),16,17 || pT of 1st and 2nd extra jet in ttbar, unmatched and matched dep on shower params, scale dep. for matched (4 plots) || 6-7/08 ||[attachment:Pt_ttbar.eps Pt_ttbar.eps]|| |
| | 88 | || 3b || 3.2 || 1 (4x2 runs),3,4 || pT of 1st and 2nd extra jet in gogo, unmatched and matched dep on shower params, scale dep. for matched (4 plots) || 8/08? ||[attachment:Pt_gogo.eps Pt_gogo.eps]|| |
| | 89 | || 4(7) || 3.3 || 1,9,12 || HTred and pTs, matched for gogo, ulul and ulul* 600 GeV (3 plots) || 11/08 ||[attachment:Htred.eps Htred.eps]|| |
| | 90 | || 5(6) || 3.3 || 1,5,6 || pT of 1st-2nd extra jet and HTred, matched for gogo at 300, 600, 1200 GeV (2x3 plots) || 11/08 ||[attachment:pt_300_600_1200.eps pt_300_600_1200.eps]|| |
| | 91 | || 6a(8) || 4 || 1a (4x2 runs) || HT(n),n=3,4,5(?) matched for gogo with decay jet #s, ETmiss unmatched and matched for sc. 1 (4 plots) || 6-7/08 ||[attachment:htcontrib-mockup.eps htcontrib-mockup.eps]|| |
| | 92 | || 6b(9) || 4 || 1b (4x2 runs) || HT(n),n=3,4,5(?) matched for gogo with decay jet #s, ETmiss unmatched and matched for sc. 2 (4 plots) || ||[attachment:htcontrib-sc2.eps htcontrib-sc2.eps]|| |
| | 93 | || 7(10) || 4.1 || 1a (4 unm runs),12 || HT(n),n=2,3,4 and jet numbers for unmatched gogo sc. 1 (4 params) and matched ululbar sc. 3 (2x4 plots) || 12/08 |||| |
| | 94 | || 8(11) || 4.2 || 1d (4x2 runs) || ETmiss unmatched and matched (1 plot) || 12/08 |||| |
| | 95 | |
| | 96 | For sec. 5 we should plot the most important quantities used in typical experimental analyses, such as ETmiss, HT(4), and what else? All after PGS. |
| | 97 | |
| | 98 | |
| | 99 | |
| | 100 | -- Main.JohanAlwall - 05 Aug 2008 |
| | 101 | |
| | 102 | * [attachment:Pt_ulul.eps Pt_ulul.eps]: Pt of 1rst and 2nd jet for ulul |
| | 103 | |
| | 104 | * [attachment:Pt_ululbar.eps Pt_ululbar.eps]: Pt of 1rst and 2nd jet for ululbar |
| | 105 | |
| | 106 | * [attachment:Htred.eps Htred.eps]: Ht red for gogo, ulul and ululbar, as well as Pt of 1rst and second, all with and without matching |
| | 107 | |
| | 108 | * [attachment:htcontrib-mockup.eps htcontrib-mockup.eps]: contrib of jets to Ht scenario 1 |
| | 109 | |
| | 110 | * [attachment:Pt_ttbar.eps Pt_ttbar.eps]: pt for ttbar |
| | 111 | |
| | 112 | |
| | 113 | * [attachment:pt_300_600_1200.eps pt_300_600_1200.eps]: characteristic of jets for in gogo for go mass=300, 600 and 1200 GeV |
| | 114 | |
| | 115 | * [attachment:Pt_gogo.eps Pt_gogo.eps]: pt of the 1rst and 2nd jet for gogo (with ME and PS) |
| | 116 | |
| | 117 | * [attachment:Pt_ulul.eps Pt_ulul.eps]: Pt of two leading jets for ulul |
| | 118 | |
| | 119 | * [attachment:Pt_ululbar.eps Pt_ululbar.eps]: pt of leading jets for ululbar |
| | 120 | |
| | 121 | * [attachment:Htred.eps Htred.eps]: Htred and Pt for gogo, ulul and ululbar, Pythia vs ME |
| | 122 | |
| | 123 | * [attachment:Pt_ttbar.eps Pt_ttbar.eps]: Pt_ttbar.eps |
| | 124 | |
| | 125 | * [attachment:Pt_ulul.eps Pt_ulul.eps]: Pt_ulul.eps |
| | 126 | |
| | 127 | * [attachment:Pt_ululbar.eps Pt_ululbar.eps]: Pt_ululbar.eps |
| | 128 | |
| | 129 | * [attachment:Pt_ulul.eps Pt_ulul.eps]: Pt_ulul.eps |
| | 130 | |
| | 131 | * [attachment:pt_300_600_1200.eps pt_300_600_1200.eps]: pt_300_600_1200.eps |
| | 132 | |
| | 133 | * [attachment:htcontrib-sc2.eps htcontrib-sc2.eps]: preliminary plot showing the contribution of decay jets to Ht and MET for scenario2 |
| | 134 | |
| | 135 | * [attachment:Htred.eps Htred.eps]: Htred.eps |
| | 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 | |
| | 164 | |
| | 165 | |
| | 166 | |
| | 167 | |
| | 168 | |
| | 169 | |
| | 170 | |