Short description of the tool
The aim is to build an event generator for quarkonium physics. Production rates are computed within the NRQCD theory (Non-Relativistic QCD, Bodwin, Braaten and Lepage, see hep-ph/9407339), where cross sections are expanded in $\alpha_s$ and $\normalsize v$, the relative velocity in the heavy quarkonium state. As a result, cross sections read
where the creation of a heavy-quark pair is described by the short-distance coefficients $\hat \sigma(Q\bar Q(n))$, whereas the non-perturbative evolution of this heavy-quark pair into a quarkonium state is encoded into the long-distance matrix elements $<O{\mathcal Q}(n)>$. The label $\normalsize n $ stands for the intermediate state of the heavy quark pair, which is usually expressed in the spectroscopic notation:
with $\normalsize S $ the spin of the heavy quark pair, $ \normalsize L $ its orbital angular momentum, $ \normalsize J $ its total angular momentum, and $ \normalsize c $ its color state. MadOnia automatically generates any tree-level amplitude for an arbitrary transition $\normalsize n$ (up to P-wave state) required for the computation of the short-distance coefficients. The amplitude generator is interfaced with MadEvent to produce unweighted events. These events are written in the same format in the lhe format and can be passed through Pythia for the showering and hadronization.
How to use the code
Edit the proc_card.dat and generate the process
- Enter the processus name in the proc_card.dat in the following format:
pp>jcc~[3P18to443]
This corresponds to the production of a $\normalsize J/\psi$ (PID=443) via an intermediate state $\normalsize 2S+1=3$, $\normalsize L=P$, $\normalsize J=1$, $\normalsize c=8$, i.e. a spin-1, P-wave, color-octet intermediate state. As another example:
pp>jbb~[1S01to551]
corresponds to the production of an $\normalsize \eta_b$ via a color-singlet transition.
Pay attention
the heavy-quark pair must be placed in that order, and at the end! For example,
pp>jc~c[3P18to443]
or
pp>cc~[3P18to443]j
will not work!
For the model, select sm_onium.
- Then type
./bin/newprocess
to generate the process
Edit the param_card.dat
If you edit the file param_card.dat, you will see that in addition to the sm parameters, you have a block called LDME. The parameters in this block correspond to the long distance matrix element, expressed in the BBL normalization (see hep-ph/9407339). You can specify the value of the LDME's for each transition. The transition $\normalsize n$ is indicated at the end of the line as a comment. For example, the line
3. 1.16E+00 # 3S11
sets
Note that only one LDME parameter is assigned to each transition ${2S+1}L_J{[c]}$. For example, for all the following processes
pp >j cc~[3S18to443] color-octet $\normalsize J/\psi$ production pp > jcc~[3S18to10441] color-octet $\normalsize \chi_{c0}$ production pp > jbb~[3S18to100553] color-octet $\normalsize \Upsilon(2S)$ production
the same LDME parameter will be used in the param_card.dat.
For the rest…
Just proceed as for a usual generation of events:
./bin/generate_events to generate events
Note that you can put cuts on the quarkonium state via the run_card.dat.
Once your events have been generated, you can use the tool MadAnalysis to draw distributions of events.
-- Main.PierreArtoisenet - 02 Mar 2009
Attachments (22)
- latexc146e812252376c83a83eecffb9998be.png (184 bytes ) - added by 13 years ago.
- latex96912292ebc2970303e369e9e2bc2638.png (222 bytes ) - added by 13 years ago.
- latex6833c9c1ddd9495c522dcfd7842f8ad1.png (178 bytes ) - added by 13 years ago.
- latex67f74d1a9d8d6769156e5ee0a69854d0.png (513 bytes ) - added by 13 years ago.
- latexc3864a509f5d79e4576f8fa92464c94f.png (975 bytes ) - added by 13 years ago.
- latex2f9aff114224d6e7bc97e052794793cc.png (531 bytes ) - added by 13 years ago.
- latexf44dea91124b061e1a66a0f7d216ee8b.png (255 bytes ) - added by 13 years ago.
- latex9a4e91c2238e56b7b726197769de7128.png (255 bytes ) - added by 13 years ago.
- latex413ef4593d80b91b5ff141f3759a03b2.png (1.1 KB ) - added by 13 years ago.
- latexa67fdebc01c6f1746f7f96855ab06567.png (308 bytes ) - added by 13 years ago.
- latex86b02a29140c0051e6b9b8492de4e189.png (192 bytes ) - added by 13 years ago.
- latexb2d2652febddb9da786a06d06569a0fd.png (416 bytes ) - added by 13 years ago.
- latex66834e4a8ff37fe4e158b33aa45b57dc.png (202 bytes ) - added by 13 years ago.
- latex2f4da44b6b814b31526c5ad697867023.png (202 bytes ) - added by 13 years ago.
- latexad644ce59195c369d4d6cd3e4c613244.png (515 bytes ) - added by 13 years ago.
- latex408f5382750626982c06632a4c5b8707.png (396 bytes ) - added by 13 years ago.
- latexffca13fb64e59050c69ee1bafe1fa7dc.png (199 bytes ) - added by 13 years ago.
- latex2d55ab09f9e80793b83aee0ebb65a9ac.png (263 bytes ) - added by 13 years ago.
- latexe20aeb8fb6ebb2bea761aea8c6f79b90.png (233 bytes ) - added by 13 years ago.
- latex5a705fa7a54641e68ae2bafc76baed16.png (437 bytes ) - added by 13 years ago.
- latexe9472dbcae11e6b3603715f03bf06bf2.png (271 bytes ) - added by 13 years ago.
- latex0577fb38ea5ad2c3c800a6a0295ca2a8.png (232 bytes ) - added by 13 years ago.
Download all attachments as: .zip