MEtop - a generator for single top production via FCNC interactions

We present a generator for single top quark production via flavour-changing neutral currents. The MEtop event generator allows for Next-to-Leading-Order direct top production $pp \to t$ and Leading-Order production of several other single top processes. A few packages with definite sets of dimension six operators are available. We discuss how to improve the bounds on the effective operators and how well new physics can be probed with each set of independent dimension six operators.


Introduction
The Large Hadron Collider (LHC) at CERN has resumed operation with its center-ofmass energy increased to 8 TeV. The LHC top factory will allow us to scrutinise the heaviest of all known quarks with unprecedented precision. Flavour physics is on the top of the agenda as one of most interesting research topic that can be address at this collider, through the study of flavour changing neutral currents (FCNC) in top quark production and decay. In fact, a wide variety of models shows a strong dependence in the measurable FCNC quantities: for instance, top quark FCNC branching ratios can vary from extremely small in the Standard Model (SM) to measurable values at the LHC in a variety of the SM extensions [1]. Therefore, the large number of top quarks being produced provide a means to search for deviation from the SM, however small they are. It is clear that the simplest way to search for new FCNC physics is to look for the rare top decays as for example in t → qγ where q = u, c is the up-quark or a c-quark. Limits on the BR(t → qγ) were set indirectly at LEP [2] and HERA [3,4] and directly at the Tevatron [5] (see [6] for references and details). Presently, the best bound on the photon FCNC current is the one from HERA while the best experimental bounds on BR(t → qZ) were obtained at the Tevatron [7] and at the LHC [8,9] . Finally, the best bound on the strong FCNC current tqg was recently obtained in direct top production at the LHC [10].
Our goal is to provide a tool to measure FCNC related to the top quark at the production level. More evolved top FCNC searches can only be performed if a dedicated generator for top FCNC studies is available. So far only the PROTOS generator [11] and the TopReX generator [12] can be used to generate events for direct top production pp(gu + gc) → t studies as well as the top FCNC decays. Therefore, we considered that it was both necessary and timely to make available a generator for top FCNC physics that included a larger set of FCNC operators together with a larger set of physical processes at the production level. This is the purpose of the MEtop event generator.
MEtop main process is direct top production, pp(gu + gc) → t, but calculated at Nextto-Leading-Order (NLO). NLO direct top was implemented by adopting an effective NLO approximation as described in [13]. Besides direct top, MEtop can be used to generate events at LO for all FCNC processes with a top and a gluon or any quark other then the top in the final state. We plan to include other processes like pp → tV , with V = γ, Z, W in the near future. From the theoretical point of view we will adopt the effective operator formalism as described in [14]. We use a set of dimension six effective operators always involving at least one top quark. The set of operators is classified in three different groups: strong, electroweak and four-fermion (4F) operators. MEtop comes with several choices of packages where different sets of effective operators are available.
The plan of the paper is as follows. The next section describes the complete set of FCNC operators needed for top physics. In section 3 we provide a detailed description of the physical processes available in MEtop. Section 4 is devoted to the implementation in MEtop of the effective NLO approximation for direct top production. In section 5 we compare the contributions of the different classes of operators to single top production. Our conclusions are drawn in section 6. There are also three appendixes dealing with more technical issues.

The FCNC operators for top physics
The effective operator formalism assumes that some general theory which has the SM as its low energy limit can be written as a series in Λ with operators obeying the SM symmetries, where L SM is the SM lagrangian and L (5) and L (6) contain all the dimension five and six operators respectively. This formalism provides a model-independent parametrisation of physics beyond the SM. This lagrangian contains only SM fields and therefore any new particles and any new interactions are hidden in the effective operators. The term L (5) vanishes if baryon and lepton number conservation is imposed. The complete set of dimension six operators is quite vast. In order to simplify the discussion we classify the operators in three categories: strong FCNC operators [15,16], the ones generating a vertex of the formtug, where g is a gluon and u is an up-quark; electroweak FCNC operators [17,18], which are the ones giving rise to a vertex with one top quark, an up-quark and one electroweak gauge boson, and finally four fermion (4F) operators which are Fermi interactions with one top quark and three other quarks.
When writing all allowed dimension six operators obeying the required symmetries of the SM lagrangian, one readily understands that not all operators are independent [14]. They are related by the equations of motion and also by Fierz transformations. Therefore, the total number of operators can be reduced to a minimum set of independent operators. Moreover, this set can be further reduced when only specific processes are studied, like in our case, where all operators have at least one top quark in the interaction. A minimal set of operators for top quark physics was discussed in [15,16,19,20] and here we will just present this minimal set according to our classification.
We will start by considering the non 4F operators. As previously stated, we divide these operators in two classes: strong FCNC operators, when the gluonic tensor is involved, and electroweak FCNC operators when electroweak gauge bosons are present in the interaction. We assume that O ij and O ji are independent operators and the hermitian conjugate of all the operators are included in the final lagrangian.
Following the notation of [14], the independent operators contributing to the strong FCNC vertices can be written as where G a µν is the gluonic field tensor, u i R stands for a right-handed quark singlet and q i L represents the left-handed quark doublet. FCNC occurs because one of the indices is always equal to 3 while the other is either 1 or 2, that is, there is always one (and one only) top-quark present in the operator; the remaining fermion field in the interaction is either a u or a c-quark. These operators will give rise to the FCNC vertices of the form g tū i (with u i = u, c) and the corresponding hermitian conjugate interaction with an independent coefficient. Operator in equation (2) also appears in the literature as a dimension 5 operator. In that case, the corresponding FCNC lagrangian is written as where κ u is real, g s is the strong coupling and f u and h u are complex numbers with |f u | 2 + |h u | 2 = 1 (see appendix A for a detailed discussion relating the forms of the strong FCNC operators).
In the electroweak sector we now have to look for FCNC vertices of the type V tū i (with u i = u, c and V = Z, γ). The minimal set of operators that give rise to the above interactions can be written as where B µν and W I µν are the U(1) Y and SU(2) L field tensors, respectively. As was shown in [19], for all the operators in (5) and (6), O ij and O ji are not independent. In fact, by writing the combinations O i+j and O i−j and using the equations of motions, it can be shown that only one of these combinations is independent. This means that the number of independent operators in (5) and (6) is reduced to three (for each light flavour). The above discussion leads us to the conclusion that the minimal number of non-4F operators needed to study top FCNC physics is 9 for each light flavour (u and c).
The equations of motion used to reduce the number of operators relate the operators from the strong and electroweak sectors with the 4F operators. The number of independent 4F operators depends on the process considered. Following [21] we found the minimal number of independent 4F operators needed for t plus quark production, pp → tq, which are shown in appendix B. The final lagrangian for the study of single top production via FCNC currents can then be written as where L 4f u and L 4f c are described in appendix B. One should note that in order to keep a manageable number of 4F operators we only consider initial states with up-quarks in the hadron colliders case. There are reasons for considering a reduced set of 4F operators, namely the ones that have in the initial state either uu and uū. First, and assuming that all 4F coupling constants are of the same order, these initial states provide the largest contribution for the cross section. Second, our main goal is to provide a means to distinguish between operators by analysing different distributions and this can only be done if the number of operators is not too large. The addition of further 4F operators will be done in the future if found necessary. The minimal number of 4F operators in the case of FCNC pp → tt was recently considered in [22]. Before ending this section we will briefly discuss the bounds on the coupling constants α i . In the effective operator approach, the lagrangian has the SM symmetries. Therefore, physics of the top quark is related with B physics. In reference [23] this relation was explored in order to constrain the electroweak FCNC operators 1 . The most constrained operators are obviously the ones built with quark doublets only while the less constrained are the ones built with quark singlets only. Consequently, for the first generation, bounds on operators with doublets only, α (3,ij) φq /Λ 2 and α (1,ij) φq /Λ 2 , are of the order 0.01 TeV −2 . For operators α ut uW φ /Λ 2 and α ut uBφ /Λ 2 the bounds are of the order 0.3 TeV −2 while for α tu uW φ /Λ 2 and α tu uBφ /Λ 2 we have 1 TeV −2 . Finally, regarding operators with singlets only, like α uφ /Λ 2 , reference [23] obtained a bound of the order 3 TeV −2 . The bounds for operators relating the second and third generation are of the same order of magnitude.
As stated in the introduction, there are new direct bounds from the LHC that lowered the limit on BR(t → qZ) to 0.34% [9]. A new indirect bound from HERA [3] is also available BR(t → qγ) < 0.5%. Also, a combined study on B physics and Tevatron data on top quark production cross section places an indirect bound on the sum of the FCNC branching ratios forcing them to be below the percent level [25]. All these new bounds do not imply any dramatic changes on the bounds in the electroweak sector.
Regarding the strong FCNC operators the most recent search is the one from the ATLAS collaboration [10] in direct top production at the LHC. The upper limits obtained at 95 % CL for the strong couplings are κ u /Λ < 6.9 × 10 −3 TeV −1 and κ c /Λ < 1.6 × 10 −2 TeV −1 which in turn can be translated into strong branching ratio bounds BR(t → ug) < 5.7 × 10 −5 and BR(t → cg) < 2.7 × 10 −4 .
Contrary to the the dimension six FCNC operators from the strong and electroweak sector, there are no useful bounds on the four fermion operators involving two top quarks and this is even more so if the top is right-handed. Therefore, the LHC can place constraints on these operators.

Physical processes
MEtop generates events according to the Von Neumann algorithm (see [26] for details). The amplitudes for each process were generated with CalcHEP [27], and the Feynman rules for the effective operators were derived with LanHEP [28]. Integrations are performed using the CUBA library [29], configured to use VEGAS algorithm [30] Generation of events for hadron colliders need the linking with the LHAPDF package [31]. The events are written in the standard format of Les Houches event file [32] (.LHE). Whenever possible our results were checked via a completely different path. First, the Feynman rules were generated by the implementation of the effective operators in UFO [33]. Then, cross sections calculation and event generation was performed using MadGraph 5 [34]. We always found an excellent agreement with MEtop.
The following hard processes are already included in MEtop: direct top production at LO and NLO both for pp and for pp colliders, which at the parton level amounts to the processes gq → t and gu → gt where q = u, c; pp(p) → qt, with all possible parton contributions in the initial state taken into account and q is now any quark other then the top-quark. The corresponding conjugate processes, with an anti-top in the final state, are also included, and the processes can be generated independently. There are several available packages in MEtop, each containing a different set of operators. The list of operators in each package is presented in appendix C.
In figure 1 we present direct top production together with the top+gluon processes. The two have in common the fact that only the strong FCNC couplings contribute to the process. Only one diagram contributes to direct top channel while several diagrams with gluon exchange contribute to top+gluon. MEtop allows for generation of events at LO and at NLO. As we have already mentioned, currently, FCNC direct top production events can be generated at LO with the PROTOS generator. In figure 2 the diagrams for top + quark production are shown. Both q 1 , q 2 and q run through all quarks other then the top-quark, that is u, d, c, s, b and the respective anti-particles. Again, conjugate processes are also included. In the diagram on the right, only strong FCNC operators are present. On the left diagram all operators can take part, including the 4F ones, contributing to LO single top production at the parton level. A detailed description of each process will be presented in appendix C.

NLO approximation to direct top production
In an event generator, initial (ISR) and final (FSR) state radiation is included trough a shower mechanism. The shower mechanism assumes collinear factorization, that is, that the real radiation process can be approximated by a branching mechanism, where the first QCD radiation is emitted by one of the legs of the Born configuration. In equation 9 we present the relation between the transition amplitude for the case of q → qg splitting where M i is the amplitude and dΦ i is the phase space for the ith body processes and P q,qg (z) is the Altarelli-Parisi splitting function. This approximation breaks down in the hard P T region where the matrix element M n+1 should be used. A factorization prescription or matching scheme is then used to merge these two regions in a smooth and optimised way (see for example the merging approaches CKKW [35] and MLM [36]).
In the previous section we have presented in figure 1 the parton level contributions to direct top production, together with the parton level contribution to the hard process gq → gt. The later process contributes to the inclusive direct top production. The gq → gt process has soft and collinear divergences, and this problem can only be solved by including the NLO corrections. Furthermore, the FCNC direct top production cross section was calculated in [37] and a considerable enhancement of about 40 % was found relative to the leading order cross section for the LHC with √ s = 14 TeV. Therefore it is desirable to have a NLO generator for direct top production at the LHC. It was also shown in [37] that the NLO QCD corrections vastly reduces the dependence of the total cross section in the renormalization and factorization scales which in turn increases the confidence in the predictions.  In order to calculate the inclusive NLO FCNC direct top cross section one has to consider the Born diagram, the virtual contributions and the real emission diagrams. A sketch of the diagrams from each of the above mentioned contribution are shown in figure 3. It is well known that at NLO, the integration in the full phase-space gives rise to infrared divergence from the virtual-born interference part and these divergences are cured by the addition of the real emission process. Although easily fixed for a total cross section analytical calculation, this problem is not straightforward to deal with at an event generator level. In fact, there is no way to extract the infrared infinities as it is usually done with dimensional regularization in a theoretical calculation. There are methods to deal with these infinities such as Phase Space Slicing [38,39,40] and Subtraction Methods [41,42]. In this work we will adopt an effective NLO approximation [13] to simulate direct top events at the NLO level. In this approach a merging scheme between 2 → 1 and 2 → 2 events is performed, where each process will separately populate two distinct but joint regions of the phase-space. A resolution parameter must then be defined, which in the present case is the cut in transverse momentum of the top quark applied to the real radiation process. This P cut T will then play a role of a matching variable, P match T . The phase-space region for small P T will be described by the 2 → 1 process and the subsequent parton shower (PS) mechanism, whereas the hard P T region will be described by the 2 → 2 process. One must then just make sure that the transition is done in a smooth way. The virtual corrections are included via a K-factor applied to the cross section of the 2 → 1 process. We assume this to be a good approximation because the kinematics of the Born and Virtual configurations of the direct top process should be very similar. The events will then be generated according to the following relation where σ LO is the tree-level direct top contribution, σ Real is the real radiation part, K is the K-factor and P P S T and P match T are the transverse momentum of the first PS emission and the integration cut of the real radiation process, respectively. Once the direct top events are produced, they will be radiated through a radiator like the one in PYTHIA [43]. In order to avoid double counting, the matching must ensure that the first PS emission from the 2 → 1 process will not fall within the 2 → 2 configuration phase-space. There are two ways of accomplishing it: either by vetoing all radiated 2 → 1 events that would be within the 2 → 2 configuration phase-space or simply by limiting the phase-space region of the radiated 2 → 1 events to the boundaries defined by the resolution variable. We choose to adopt the later.
In order to follow this approach, one must ensure that the PS mechanism added to the generated events from the Born configuration will populate the region with P T < P match T , which can be assured using a PT-ordered shower [44], available in both current PYTHIA versions 6.4 and 8.1. We therefore assume that the generated events will be showered by a PT-ordered mechanism. Therefore we start by calculating the three cross sections from equation (10), with P cut T = P match T for the 2 → 2 process. For the σ N LO cross section we have used the expressions from [37], where the top quark is on-shell. The tree-level direct top and top+gluon amplitudes were generated with CalcHEP where the top quark and the W decays were included in order to preserve spin correlations. The cross sections are then calculated with the Cuba library. Hence, the K factor is calculated "on the fly" for each sub-process. After extracting the K-factor, the events are generated weighed according to equation 10. The PS starting scale can then be configured to start the branching in P match T for the 2 → 1 events configuration, which in MEtop is done by preparing the .lhe files to be used by PYTHIA. A short remark is in order -in the 2 → 1 configuration, no meaningful FSR from the top quark can be present due to its large mass. Hence, we consider a good approximation to take only ISR into account. In figure 4, the black solid line represents the P T distribution of the top quark in direct top production, after the first branching in ISR, with starting scale of m t . In the same figure, the blue dashed line represents the hard process: top+gluon production. As described previously, P T is the kinematical variable chosen to match the two processes avoiding double counting in the low P T region. In figure 5 we present the P T distribution of the top quark after the first ISR branching with a P match T of 10 GeV (left) and 15 GeV (right). The natural criterion to determine the value of the P T matching parameter in the effective NLO approximation is the smoothness of the transverse momentum distribution. There are no significant differences when the value of P match T is varied in the 5 GeV to 20 GeV range. As can be seen in figure 5 , there is never a completely smooth transition between the two sets of events. This effect should be included as part of the systematic uncertainties. This feature was checked for a large range of P T match. After including the full shower (ISR+FSR) and Multiple Interaction (MI) we have opted for a value of P match T of 10 GeV.  In figure 7 we present the LO and NLO P T (left) and η (right) distributions of the top quark at the partonic level after the full shower and MI. It is clear from figure 7 that the P T and η distributions of LO direct top production are quite different from the corresponding NLO direct top ones. In fact, the distributions show that the use of a constant K factor does not correctly describe the behaviour of direct top at NLO. Hence, a new analysis is needed to improve the accuracy of the bounds on the strong coupling constants κ u and κ c . The direct top NLO p T distribution is shifted to low values of p T as compared to the LO distribution while the η distributions are shifted to higher values of η as compared to the LO one.
The actual experimental analysis is performed by looking at the distributions of the final state particles. Therefore, in figure 8 we present the comparison between LO and NLO P T (left) and η (right) distributions of the lepton from t → bW → blν at the partonic level after the full shower and MI. Again, it is clear that the level of improvement by considering the NLO distributions heavily depends on the particular analysis being performed.  We have just described how we generate a sample of inclusive direct top production at NLO. However, if the goal is to set a limit on the strong FCNC coupling, one needs to add the events generated in the process pp → t + jet composed by the parton level processes gg → t + q and q 1 q 2 → t + q, where q, q 1 , q 2 are quarks other then the top and hermitian conjugate reactions are included. Although the main contribution to top + jet production comes from the reaction where the jet is a gluon, all processes where the strong FCNC operator intervenes should be taken into account in all analysis. As these processes also suffer from infrared and collinear divergences, we have decided to avoid them by using a similar cut to that of P match T , that is, p T > 10 GeV. The complete NLO QCD corrections to the FCNC process of top+jet production were presented in [45]. The corrections can increase the cross section by 10 % to 30 % at the LHC@14TeV.
When generating the top + quark subprocesses we have to decide what is considered as signal in our analysis. pp → tq has three different classes of subprocesses: the ones which are exclusive to the Standard Model, like ud → tb, the ones that are originated exclusively via FCNC interactions, e.g. uu → tu, and the ones where interference occurs, like ub → tb. We define as FCNC signal the contributions from the two later classes of subprocesses. For the pure FCNC processes this poses no problems. However, for the interference terms this procedure leads to the inclusion of a small portion of events that will also be counted as background. However, choosing the effective strong coupling constants as κ u = 0.01 (Λ = 1 TeV) , f u = 1/ √ 2 and h u = 1/ √ 2 and for a CM energy of 7 TeV, the pure FCNC cross section is 8.718 pb, the interference term is 1.205 pb while the SM contribution amounts to only 0.018 pb. Hence, the SM contributions can be safely neglected. In figure 10 we show the P T and η distributions for the direct top at NLO summed with pp → tq for a P match T = 10 GeV and the jet p T > 10 GeV. It is clear that the shape of the distributions do not change much with the inclusion of the pp → tq process but still the pp → tq process gives a contribution of the order of 10 % to the total cross section of the inclusive top production at the LHC at √ s = 7 TeV.

Single top beyond the strong FCNC operators
In the previous sections we have discussed NLO direct top and t + q production when only the strong FCNC operator is considered. We note that the leading order contribution to direct top does not receive contributions from other operators. Therefore, the NLO calculation is again performed with only the strong FCNC operator. As long as no excess is found at the LHC in the single top channel, the procedure described in the previous section gives us the best possible bound on the anomalous strong FCNC coupling when all other operators are discarded.
In the hard P T region, the process pp → t + jet gets contributions from the complete set of independent operators. As these operators are independent from each other (and therefore so are the respective coupling constants) the interference terms between strong and electroweak or 4F could be sizeable. If an excess is found in the single top channel, one has to take into account all possible contributions from the remaining operators. A thorough analysis of the distributions of each individual operator will help us understand which operators could be important for a given experimental analysis. Moreover, even if an excess is not seen in the single top channel, dedicated analysis could most probably help constraining definite sets of operators. We start by considering the strong operator. When all other operators are turned off the P T and η distributions have a very mild dependence on the strong coupling constant κ (κ u to be more precise, and we have set κ c = 0). This is shown in figure 11 where the P T (left) and η (right) distributions of the top quark are shown for three values of κ, 0.001, 0.01 and 0.1 and Λ = 1 TeV. The process is direct top NLO plus pp → tq for √ s = 7 TeV and P cut T = P match T = 10 GeV. It is clear that the shape of the distributions does not vary much making it possible to perform the analysis for one constant and then to extract a bound on the strong operator.
We now move to the study of the electroweak operators. We first consider only one operator O uW φ turned on. In figure 12 we present the P T (left) and η (right) distributions of the top quark for three values of α uW φ , 0.01 , 0.1 and 1 and Λ = 1 TeV. As α uW φ → 0 we recover the pure SM contribution of electroweak origin. The SM cross section for this process and for 7 TeV is σ = 0.019 pb while the total cross section for α uW φ = 0.01, 0.1 and 1 are σ = 0.0020 pb, 0.148 pb and 12.4 pb respectively. Therefore the different shapes of the P T and η distributions are due to the interference with the SM contribution. When α uW φ = 1, the total cross section is almost 100 times larger than its pure SM counterpart. Therefore, this value shows how the distribution behaves when the SM contribution is negligible.
This kind of behaviour can occur for any operator on the list, provided that the coupling constants are such that SM and FCNC cross sections are of the same order of magnitude. Any deviation relative to the SM showing up in the distributions could mean an interference with one or more operators. Understanding the different distribution requires dedicated studies with no assurance however that the responsible operators could be identified. One should emphasise that a thorough study of the P T and η distributions of the top quark could help identifying classes of FCNC operators. A similar discussion applies to the 4F operators case.
One should also note that pp → tq does not include direct top, because the strong operator is turned off. Contrary to strong operator scenario, in this case the distributions change with the value of the electroweak constant. Therefore, any bound based on the production process has to take into account that different coupling constants can lead to different distributions.
Having studied the distributions of a definite operator representative of each class, we will now perform a comparison between classes. In figure 13 we present the P T (left) and η (right) distributions of the top quark when just one operator is taken non-zero at a time. We compare the distributions of the strong FCNC operator with one electroweak operator (with coupling constant α uW φ ), and one 4F operator, (ū γ µ γ L u) (ū γ µ γ L t). It is clear from the figure that the distributions can be quite different and therefore distinguishable to some extent. The ability to distinguish the different operators depends heavily on the relative values of the coupling constants. If an excess in single top production is seen we can try to understand its origin by looking at all possible distributions. However, this will always be a hard task because different operators give similar distributions and therefore only very particular scenarios can be probed.

Conclusions
We have presented a new generator for the study of FCNC top interactions. The generator MEtop comes with different packages, each with a subset of a complete set of dimension six operators. At the moment MEtop can generate events for direct top and for top plus jet production, where the jet can be any quark other then the top or it can be a gluon. The direct top production process is implemented at NLO using an effective NLO approximation. Also, the inclusive contribution to direct top coming from pp → tq can be included in the event generation. We have shown that the top P T and η distributions show clear differences when the events are generated at LO or at NLO. Therefore, the use of a constant K-factor does not provide an accurate description of direct top production at NLO. We conclude that a new experimental analysis is in order to improve the constraints on the strong FCNC coupling constants. The inclusion of the inclusive process pp → tq will further improve this bound. We note that a detailed study of the P T and η distributions of the top quark could help identifying classes of FCNC operators .
At LO, the contributions stemming from the different operators can be compared in the single top production process. In particular, 4F operators can be for the first time constrained at hadron colliders. Constraining the 4F operators can help us understand their role in the asymmetry measured at the Tevatron.
The bounds on BR(t → u(c)Z) and BR(t → u(c)γ) are obtained in the process pp → tt where one of the top-quarks decays as t → bW while the other decays as t → u(c)Z or t → u(c)γ. This means that all electroweak FCNC couplings always appear in the same combination. With MEtop we are able to look for distribution that isolates each electroweak FCNC operator. This way more detailed information can be obtained about each operator.
New final states with FCNC contributions, like for instance pp → tW [47] , are to be included in the next version of MEtop.

Acknowledgments
We thank A. Belyaev, A. Semenov, T. Sjöstrand and Li Lin Yang for discussions. We thank C. Friedrich and D. Hirschbuehl for their help in the debugging and cross-check with existing generators like PROTOS and also for their support. The

A Effective couplings translation
The representation of the effective coupling constants for each operator is arbitrary. In this appendix we will relate the two most common representations of the strong FCNC coupling constants appearing in the literature. The translation is simple and our goal is to clarify the relation between experimental bounds (and theoretical bounds as well) present in the literature. When the effective strong operator comes from a dimension six effective lagrangian it is usually represented in the form where Λ is the scale of new physics. This operator can be written as a dimension five-like operator when the scalar field,φ is replaced by (v/ √ 2 0), resulting in There are in principle four complex constant involved, α it and α ti with i = u, c in a total of four degrees of freedom. On the other hand, several authors adopt to write the same strong operator as a dimension five effective operator. In this case it is usually written as where g S is the strong coupling constant, κ ij is taken as real and positive, f ij and h ij are complex and |f ij | 2 + |h ij | 2 = 1. In both cases T r[λ a λ b ] = δ ab /2, and the vacuum expectation value is v = 246 GeV. Note that in reference [14] (as in other references) the vacuum expectation value is defined as v = 246/ √ 2 GeV. The constant κ ij is real and only two constants are needed, one for each light flavour; the same is true for the complex constants f ij and h ij . Therefore we can use just the light quark index to represent those coupling constants and in this case the index j refers to the top quark. It is now straightforward to find the relation between the two set of operators with i = u, c. In most cases all constants are taken as real. This means that equation 14 can be written as and consequently with i = u, c, all constants are now real and |f i | 2 + |h i | 2 = 1.

B The complete dimension six lagrangian for single top production
As previously discussed, the number of effective dimension six operators is huge. Therefore, no meaningful analysis is possible when all operators are considered simultaneously. Any subset of operators, however small, does not simplify much the task of obtaining information about each individual operator. In order to perform any relevant study involving all different types of 4F operators, we have built a 4F subset to be used as a basis for our study. Since we are dealing with hadron colliders, our first simplification is to consider only processes initiated by u-quarks. This is equivalent to say that the coupling constants in the 4F sector are all of the same order. Using references [14,20,21] we extract nine 4F effective operators that could contribute to single top production at hadron colliders. This set can be written as where i, j, k, l = 1, 2, 3 are flavour indices and the sub-indices a and b indicate the contraction of color indices whenever this pairing is different from the one in spinor contraction.
In [21] this set of operators was simplified to the one presented in table 1. There is a total Table 1: 4F operators contributing to single top production with parton level processes of the type u(ū)u(ū) → t(t)u(ccū).
of 24 different combinations coming from setting two of the indices i, j, k equal to 1 while the remaining one is set to 1 or 2. By forcing the initial state to be composed of u-quarks only, we can further reduce the number of operators to 12 (we only allow for one "FCNCcurrent", (ūΓt) or (cΓt), where Γ stands for a generic Lorentz structure). Therefore the final 12 independent 4F operators are the ones obtained by setting k = j = 1 and i = 1, 2. Because the operators had to be rearranged to allow for the implementation in LanHEP (see discussion below), the final lagrangian is then written as where L 4f u is the 4F lagrangian for anomalous top-up coupling and L 4f c is the 4F lagrangian for anomalous top-charm coupling The operators with the Gell-Mann matrices originate from re-writing the ones where the quark colours indices were explicitly summed. The inclusion of 4F operators in MEtop was done by implementing the 4F effective lagrangian in LanHEP. All 4F operators in table 1 have four coloured particles converging in one point which is a type of interaction LanHEP is not able to handle automatically due to the complex color flow. Therefore, we had to implement these operators using an auxiliary field mechanism [46], where the 4-color vertex is replaced by 3-color vertices that when combined in s,t and u channels, will reconstruct the 4-fermion interaction. These 3-color vertices are implemented by introducing the interaction terms in the initial lagrangian together with a unit mass field with a point-like propagator. An example of how a lagrangian is written is shown in equation 31 with a vectorial auxiliary field where ψ L is a left-handed spinor and X µ is a spin 1 field that does not propagate.
C Using MEtop

C.1 Installation
MEtop is written in C and python and it generates events following the Les Houches Accord format. It can therefore be easily interfaced with PYTHIA or Herwig. In order to compile it, you need a C compiler 2 and python version 2.6 or later. To run the package you must additionally install • Cuba Library version 3.0 • LHAPDF version 5.8.6 • Numpy version 1.3.0 The Cuba and LHAPDF library must be available through the library environment variable (for example).
To install MEtop you just have to execute "make" in the main directory. In MEtop all parameters are set in one file: "param.dat". Table 2 summarizes the definition of each parameter.

C.2.2 Physical processes
In addition to the parameters defined in table 2 there are two more flags in "param.dat" file: "cs" and "Process". The first one dictates whether or not to calculate the cross sections and/or to generate events. The second sets which physical process should be taken into account. If "cs" is set to 0, the cross sections for all sub-processes defined by the "Process" flag will be calculated and no generation will be performed. The result will be stored in the CS folder, in a csX.txt file, where X can be "Dtop","Gtop" and "Lqtop". If "cs" is set to 1, only the event generation will be performed. In this case events are produced according to the calculated cross sections. After generation, the LHE files will be stored in the Events folder together with a file "runinfo.txt" which stores all information related to the event generation.  top quark FCNC interactions were introduced in MEtop through an effective lagrangian. Depending on which operators are "turned on", different physics will be generated. Two different topologies are available: 2 → 1 → 3 and 2 → 2 → 4. 3 The first one concerns "Direct top" production, and the second is related to "top+gluon" and "top+ light quark".  Strong FCNC top interactions are included in MEtop through two equivalent effective operators, one for the top and u-quark interaction, and the other for the interaction of the top with a c-quark. In process 1, only the strong coupling constants are needed.
Process 2 has the same effective operators but due to the infrared divergences appearing in top+gluon production a cut in the top quark transverse momentum has to be set via the variable PTmin. In process 3, top + light quark 4 production, all operators can contribute, strong, electroweak and 4F. It is now possible to choose which operators to include. Again a value for PTmin has to be chosen.
Process 21 is inclusive direct top production at NLO and again only strong operators intervene. The NLO result is obtained by a matching procedure (as described previously) which depends on one variable, PTmatch, to be chosen by the user. The cross section results are written in three files: "csDtopLO.txt", the LO result for direct top, "csDtopNLO.txt", the NLO increment relative to the LO result (σ T otal N LO − σ T otal LO ), and "csGtop.txt", the LO cross section for "top+gluon" process with a top quark transverse momentum above "PTmatch". Therefore the variable "PTmin" is irrelevant for this process. After the generation, the results are stored in one file in the Events folder named "DtopNLO.lhe", containing 2 → 1 → 3 and the 2 → 2 → 4 configurations. These events constitute the inclusive direct top NLO event generation, and must subsequently be showered by PYTHIA using the PT-ordered scheme, in order to complete the matching procedure. Finally, with process 22, MEtop sums process 21 with process 3. The "PTmacth" variable plays the same role as in process 21 and "PTmin" will be the top transverse momentum cut, for the "top + quark" sub-processes.

C.3 Running MEtop
To run the package you just have to execute the command "./run.py" in the main directory. Care should be taken when changing the values of the physical parameters and/or the process you wish to calculate. In such cases you must always recalculate the value of the cross section. In addition, if you change the process used for the generation, you must be sure that all cross sections pertaining the new process are calculated beforehand. This is mandatory because the generation is done using the cs*.txt files saved in the CS folder.

C.4 Available Model files
At the moment there are three different packages available in MEtop, with a different set of effective operators. The main reason to have the different packages is to make the generation of events faster. The lightest version includes only the strong sector. Then there are two other versions one with strong plus electroweak operators and the other one with strong plus 4F operators.
In equation 32 we present the strong FCNC lagrangian as it is written in the package "MEtop S vxx.tar.gz" with O uG = i g s Λū λ a σ µν (f u + hu γ 5 )tG a µν , O cG = i g s Λc λ a σ µν (f c + hc γ 5 )tG a µν (33) 4 Here light quark stands for the set u,ū, d,d, c,c, s,s, b,b. and co1, co2, f u, hu, f c, hc are real constants to be chosen in the file param.dat. The constants f i , h i allow the choice of different chiralities while co i are overall normalization constants. Although it usually considered that f 2 i + h 2 i = 1, this relation has to be implemented by the user by a judicious choice of parameters f i and h i .
The package "MEtop SEW vxx.tar.gz" contains the strong and electroweak sectors. The lagrangian introduced in this package is where the electroweak operators are where the 4F lagrangians were presented in equations 29 and 30. The relation between the parameters in equations 29 and 30 and the corresponding co i parameters in the param.dat file is shown in