Changes between Version 2 and Version 3 of AntwerpCapitaSelecta08
- Timestamp:
- Apr 17, 2012, 11:13:07 AM (12 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
AntwerpCapitaSelecta08
v2 v3 13 13 The biggest difficulty in generating events with multiple jets is to prevent double counting. For example, an event with three jets in the final state could have come from a three parton configuration of which all three showered into different jets, or a 2-parton event of which a third jet is generated by the parton shower (or even an 1 (or zero) parton event with two (or three) jets generated by the shower). The otherway around is also possible. A two jet final state could also have come from a three parton state, of which two parton are belonging to the same jet, and this could therefore lead to double counting with a two-parton sample of which both the partons are showered to a two jet final state. 14 14 15 There are basically several schemes to avoid this double counting. MadGraph uses a method based on the MLM scheme (also a version of the CKKW method is implemented for more details see [wiki: Software.Matching this page]). It works as follows.15 There are basically several schemes to avoid this double counting. MadGraph uses a method based on the MLM scheme (also a version of the CKKW method is implemented for more details see [wiki:Matching this page]). It works as follows. 16 16 17 Multi-parton events are produced by first generating all the diagrams belonging to the +0, 1, 2, 3, ... jet configurations. Any number of additional particles can be included, the more the 'better' the result will be. In practice, up to 3 jets is the maximum due to the large number of diagrams contributing to these high multiplicity states. For the $p\bar{p} \to Z/\gamma^* \to \mu^+\mu^- +0,1,2,3\textrm{ jets}$ sample, the proc_card can be found [attachment: ../../Physics/YETI08/proc_card.dat here]. Note that all the diagrams with multiplicities up to three jets will be generated (a couple of thousand diagrams) and this takes in general 1 or 2 hours or so.17 Multi-parton events are produced by first generating all the diagrams belonging to the +0, 1, 2, 3, ... jet configurations. Any number of additional particles can be included, the more the 'better' the result will be. In practice, up to 3 jets is the maximum due to the large number of diagrams contributing to these high multiplicity states. For the $p\bar{p} \to Z/\gamma^* \to \mu^+\mu^- +0,1,2,3\textrm{ jets}$ sample, the proc_card can be found [attachment:proc_card.dat here]. Note that all the diagrams with multiplicities up to three jets will be generated (a couple of thousand diagrams) and this takes in general 1 or 2 hours or so. 18 18 19 19 The second step is, as usual, to generate events by integrating the matrix elements over the phase space. Later-on we will remove events where the 'distance' between two partons is very small, so to improve the effficiency a cut has to be put on the 'kT'-distance between two particles that could form a jet togheter. This 'kT'-distance between particle i and j is defined as $k_T^{ij}=\sqrt{2\textrm{ min}[p_T^i,p_T^j]^2(\cosh(\eta^i-\eta^j)-\cos(\phi^i-\phi^j))^2}$. A cut on this variable can be set with the 'xqcut' parameter in the run_card. Altough this is a non-phyical cut, and the final result should be independent of its value, chosing it in an appropriate way improves the efficiency quite a lot. Note that for consistency it has to be smaller than other scales in the process, in particular it should be smaller than the jet measure cut-off (described later). At this stage also the minimum transverse momentum (and possibly invariant mass) of the jets has to be set equal to the xqcut to make sure that there are no gaps in the phase space and increase the efficiency. However, the dR-distance can be set (almost) to zero (e.g. 0.001), because the xqcut already takes care of this. … … 21 21 The third step is to shower the events. Now we have to make sure that we remove events that 'change the number of jets' to avoid the couple counting. A showered event is clustered to jets using a k_T jet algorithm with a jet measure cut-off (generally speaking 1.5-2 as large as the xqcut). This cut-off makes sure that no new jets are generated from a given parton level multiplicity: if an event generates an extra jet, it is removed from the event sample (except for the highest multiplicity samples, so that also events with 4, 5, 6 or even more jets are generated from the 3 jet sample). Because no new jets are generated for the low multiplicity samples, all double counting is removed. And because the minimum transverse momentum of a jet was equal to the xqcut, the full phase space is covered, and we are not removing too many events. 22 22 23 For more technical details also have a look at this page about [wiki: Software.Matching Matching]. Also [wiki:Software.MatchChecker this page] has a lot of information about matching/merging and testing your event samples.23 For more technical details also have a look at this page about [wiki:Matching Matching]. Also [wiki:MatchChecker this page] has a lot of information about matching/merging and testing your event samples. 24 24 25 25 26 26 ==== $p\bar{p} \to Z/\gamma^* \to \mu^+\mu^- +$ (no matching) ==== 27 [attachment: ../../Physics/YETI08/run_01_banner.txt run_01_banner.txt]: Z banner28 * at hadron level with jets above $p_t >10$ GeV [attachment: ../../Physics/YETI08/run_01_pythia_events.lhe.gz event file], [attachment:../../Physics/YETI08/run_01_pythia_events_plots.ps plots].29 * at detector level with jets above $p_t >10$ GeV [attachment: ../../Physics/YETI08/run_01_pgs_events.lhco.gz event file], [attachment:../../Physics/YETI08/run_01_pgs_events_plots.ps plots].27 [attachment:run_01_banner.txt run_01_banner.txt]: Z banner 28 * at hadron level with jets above $p_t >10$ GeV [attachment:run_01_pythia_events.lhe.gz event file], [attachment:run_01_pythia_events_plots.ps plots]. 29 * at detector level with jets above $p_t >10$ GeV [attachment:run_01_pgs_events.lhco.gz event file], [attachment:run_01_pgs_events_plots.ps plots]. 30 30 31 31 32 32 ==== $p\bar{p} \to Z/\gamma^* \to \mu^+\mu^- +0,1,2,3\textrm{ jets}$ matched ==== 33 [attachment: ../../Physics/YETI08/run_01_1_banner.txt run_01_1_banner.txt]: Z+0,1,2,3 jets banner34 * at hadron level with jets above $p_t >10$ GeV [attachment: ../../Physics/YETI08/run_01_1_pythia_events.lhe.gz event file], [attachment:../../Physics/YETI08/run_01_1_pythia_events_plots.ps plots].35 * at detector level with jets above $p_t >10$ GeV [attachment: ../../Physics/YETI08/run_01_1_pgs_events.lhco.gz event file], [attachment:../../Physics/YETI08/run_01_1_pgs_events_plots.ps plots].33 [attachment:run_01_1_banner.txt run_01_1_banner.txt]: Z+0,1,2,3 jets banner 34 * at hadron level with jets above $p_t >10$ GeV [attachment:run_01_1_pythia_events.lhe.gz event file], [attachment:run_01_1_pythia_events_plots.ps plots]. 35 * at detector level with jets above $p_t >10$ GeV [attachment:run_01_1_pgs_events.lhco.gz event file], [attachment:run_01_1_pgs_events_plots.ps plots]. 36 36 37 37 38 38 ==== Example difference with and without matching ==== 39 39 40 [attachment: ../../Physics/YETI08/plots.ps Here] is the (reconstructed) transverse momentum of the Z boson plotted for the two event samples. From this plot it is clear that the event sample generated without the matching clearly underestimates the number of event for the high pt region. This is exactly as expected, because the pt of the Z for the sample without matching is purely generated from the parton shower. As we know the shower only generates partons in the soft and collinear region, hence the high pt region is underestimated. Using the matching ensures that at the parton level higher multiplicities are already included, so the matrix elements itself can already generate a non-zero pt for the Z boson, therefore it is much larger.40 [attachment:plots.ps Here] is the (reconstructed) transverse momentum of the Z boson plotted for the two event samples. From this plot it is clear that the event sample generated without the matching clearly underestimates the number of event for the high pt region. This is exactly as expected, because the pt of the Z for the sample without matching is purely generated from the parton shower. As we know the shower only generates partons in the soft and collinear region, hence the high pt region is underestimated. Using the matching ensures that at the parton level higher multiplicities are already included, so the matrix elements itself can already generate a non-zero pt for the Z boson, therefore it is much larger. 41 41 42 42 43 43 === Back-up Files === 44 Events for p\bar{p}>Zh>bb~l+l- at [attachment: ../../Physics/YETI08/unweighted_events.lhe.gz parton level] ([attachment:../../Physics/YETI08/unweighted_events_plots.ps plots]),45 [attachment: ../../Physics/YETI08/pythia_events.lhe.gz hadron level] ([attachment:../../Physics/YETI08/pythia_events_plots.ps plots])46 and [attachment: ../../Physics/YETI08/pgs_events.lhco.gz detector level] ([attachment:../../Physics/YETI08/pgs_events_plots.ps plots]). Generated with the standard set of cuts (i.e. pt(l)>10GeV, ||eta(l)||<2.5 and dR(ll)>0.4, with no cuts on the b quarks (at parton level))44 Events for p\bar{p}>Zh>bb~l+l- at [attachment:unweighted_events.lhe.gz parton level] ([attachment:unweighted_events_plots.ps plots]), 45 [attachment:pythia_events.lhe.gz hadron level] ([attachment:pythia_events_plots.ps plots]) 46 and [attachment:pgs_events.lhco.gz detector level] ([attachment:pgs_events_plots.ps plots]). Generated with the standard set of cuts (i.e. pt(l)>10GeV, ||eta(l)||<2.5 and dR(ll)>0.4, with no cuts on the b quarks (at parton level)) 47 47 48 [attachment: ../../Physics/YETI08/s1_higgs.tar.gz Here] can the model for the "spin-1 Higgs" be found and some parton level [attachment:../../Physics/YETI08/spin-1_higgs_events.lhe.gz events] and the [attachment:../../Physics/YETI08/spin-1_higgs_plots.ps plots].48 [attachment:s1_higgs.tar.gz Here] can the model for the "spin-1 Higgs" be found and some parton level [attachment:spin-1_higgs_events.lhe.gz events] and the [attachment:spin-1_higgs_plots.ps plots]. 49 49 50 50 … … 53 53 54 54 55
