arxiv: v3 [hep-ph] 17 Oct 2017

Transcription

1 Photon and Pomeron induced production of Dijets in pp, p and collisions E. Basso, V. P. Goncalves,. K. Kohara, and M. S. Rangel Faculdade de Ciências Exatas e ecnologias, Universidade Federal da Grande Dourados (UFGD), Caixa Postal, Dourados, CEP , MS, Brazil. Instituto de Física e Matemática, Universidade Federal de Pelotas (UFPel), Caixa Postal, CEP , Pelotas, RS, Brazil Instituto de Física, Universidade Federal do Rio de Janeiro (UFRJ), Caixa Postal 88, CEP 9-97, Rio de Janeiro, RJ, Brazil and Centro Brasileiro de Pesquisas Físicas (CBPF), CEP 90-80, Rio de Janeiro, RJ, Brazil arxiv:70.088v [hep-ph] 7 Oct 07 In this paper we present a detailed comparison of the dijet production by photon photon, photon pomeron and pomeron pomeron interactions in pp, p and collisions at the LHC energy. he transverse momentum, pseudo rapidity and angular dependencies of the cross sections are calculated at LHC energy using the Forward Physics Monte Carlo (FPMC), which allows to obtain realistic predictions for the dijet production with two leading intact hadrons. We obtain that γip channel is dominant at forward rapidities in pp collisions and in the full kinematical range in the nuclear collisions of heavy nuclei. Our results indicate that the analysis of dijet production at the LHC can be useful to test the Resolved Pomeron model as well as to constrain the magnitude of the absorption effects. PCS numbers:.0.nn,.8.ni,.8.qk,.87.ce I. INRODUCION he experimental results from evatron, RHIC and LHC for exclusive processes, characterized by a low hadronic multiplicity, intact hadrons and rapidity gaps in final state, has demonstrated that the study of these processes is feasible and that the data can be used to improve our understanding of the strong interactions theory as well constrain possible scenarios for the beyond Standard Model physics (For a recent review see, e.g. Ref. []). In particular, it is expected that the forthcoming data can be used to discriminate between different approaches for the pomeron, which is a long-standing puzzle in the Particle Physics []. his object, with the vacuum quantum numbers, was introduced phenomenologically in the Regge theory as a simple moving pole in the complex angular momentum plane, to describe the high-energy behaviour of the total and elastic cross-sections of the hadronic reactions []. Due to its zero color charge, the pomeron is associated with diffractive events, characterized by the presence of large rapidity gaps in the hadronic final state. One good testing ground for diffractive physics and for the nature of the pomeron (IP ), is the dijet production in hadronic collisions. his process provides important tests of perturbative QCD and is one of the most important backgrounds to new physics processes. hese aspects have motivated the development of an extensive phenomenology for this process in the last years [ ]. In particular, dijet production by photon - pomeron interactions in ultraperipheral pp/p/ collisions, characterized by two intact hadrons and two rapidity gaps in the final state, has been recently investigated in Ref. [], considering the Resolved Pomeron model, in which the pomeron is assumed to have a partonic structure, as proposed by Ingelman and Schlein [] many years ago. hey have obtained large values for the cross sections as a function of various variables. he promising results presented in Ref. [] motivate a more detailed analysis of the dijet production, taking into account the contribution of other processes that are characterized by the same topology. In what follows, we will estimate the dijet production in photon photon and pomeron pomeron interactions present in pp/p/ collisions and compare the predictions with those for the dijet production in photon pomeron interactions. hese different processes are represented in Fig., and as emphasized before, they are characterized by two hadron intacts in the final state as well as two rapidity gaps. One importance difference between the dijet production by γγ interactions and the other processes, is that in pomeron induced processes, the Resolved Pomeron model predicts the existence of particles accompanying the dijet, with the associated rapidity gaps becoming, in general, smaller than in the γγ case. dditionally, the photon and pomeron induced processes are expected to generate emerging hadrons with different transverse momentum distributions, with those associated to pomeron induced having larger transverse momentum. Consequently, in principle, it is possible to introduce a selection criteria to separate these different contributions for the dijet production. lthough these distinct processes have been studied separately by several groups in the last years, the calculations have been performed considering

2 γ P X γ j P Y j j j j j j γ B j X X γ B P B P B B B B B B B B (a) (b) (c) (d) B FIG. : Dijet production by (a) photon photon, (b) pomeron photon, (c) photon pomeron and (d) pomeron pomeron interactions in hadronic collisions. different approximations and assumptions, which makes difficult the direct comparison between its predictions. Our goal in this paper is to estimate these processes considering the same set of assumptions for the pomeron and for the photon flux and obtain realistic predictions for the dijet production in photon and pomeron induced interactions including experimental cuts in the calculations. In order to do that, we will use the Forward Physics Monte Carlo (FPMC), proposed some years ago [7] to treat pomeron pomeron and photon photon interactions in hadronic collisions and recently improved to also include photon pomeron interactions in pp collisions [8]. Here we generalize this Monte Carlo to treat γγ, γip and IP IP interactions in p and collisions. s a consequence, it is possible to estimate the contribution of the different processes presented in Fig. in a common framework. In this paper we will perform a comprehensive analysis of the transverse momentum and pseudo rapidity distributions for the different processes. he content of this paper is organized as follows. In the next section we present a brief review of the formalism for the dijet production in photon and pomeron induced interactions in hadronic collisions. In Section III we present our predictions for the pseudo rapidity and transverse momentum distributions for the dijet production in pp/p/ collisions at LHC energies, considering the contributions associated to γγ, γip and IP IP interactions. Finally, in Section IV we summarize our main conclusions. II. DIJE PRODUCION IN PHOON ND POMERON INDUCED INERCIONS t high energies, a ultra relativistic charged hadron (proton or nuclei) give rise to strong electromagnetic fields, such that the photon stemming from the electromagnetic field of one of the two colliding hadrons can interact with one photon of the other hadron (photon - photon process) or can interact directly with the other hadron (photon - hadron process) [9, 0]. In these processes the total cross section can be factorized in terms of the equivalent flux of photons into the hadron projectiles and the photon-photon or photon-target production cross section. In particular, the dijet production by γγ interactions at high energies in hadronic collisions, represented in Fig. (a), can be described at leading order by the following expression σ(h h B h j j h B ) = dx dx B γ (x, µ ) γ B (x B, µ ) ˆσ(γγ j j ), () where γ i (x i, µ ) is the equivalent photon distribution of the hadron i, with x i being the fraction of the hadron energy carried by the photon and µ has to be identified with a hard scale of the process. Moreover, represents the presence of a rapidity gap in the final state and ˆσ is the partonic cross section for the γ γ B j j subprocess. On the other hand, the cross section for the dijet production in photon pomeron interactions, represented in Figs. (b) and (c), is given by σ(h h B h j j X h B ) = dx dx B [g D (x, µ ) γ B (x B, µ ) + γ (x, µ ) gb D (x B, µ )] ˆσ(γg j j ), () where g D i (x i, µ ) is the diffractive gluon distribution of the hadron i with a momentum fraction x i and we take into account that both incident hadrons can be a source of photons and pomerons. Finally, the cross section for the dijet

3 production in double diffractive processes, represented in Fig. (d), can be expressed by σ(h h B h Xj j Y h B ) = dx dx B g D (x, µ ) gb D (x B, µ ) ˆσ(gg j j ), () where, for simplicity, we assumed that the dominant subprocess is the gg j j interaction, which is a good approximation at high energies. However, in the numerical calculations, the contribution associated to the q q j j subprocess also have been included. he basic ingredients in the analysis of these photon and pomeron induced processes are the equivalent photon distribution of the incident hadrons γ(x, µ ) and its diffractive gluon distributions g D (x, µ ). s our goal is to calculate the cross sections for the processes presented in Fig. considering pp, p and collisions, we should to specify the associated models used in the proton and nuclear cases. Initially, lets present the models used for the photon distribution. he equivalent photon approximation of a charged point-like fermion was formulated many years ago by Fermi [] and developed by Williams [] and Weizsacker []. In contrast, the calculation of the photon distribution of the hadrons still is a subject of debate, due to the fact that they are not point-like particles. In this case it is necessary to distinguish between the elastic and inelastic components. he elastic component, γ el, can be estimated analysing the transition h γh taking into account the effects of the hadronic form factors, with the hadron remaining intact in the final state [0, ]. In contrast, the inelastic contribution, γ inel, is associated to the transition h γx, with X h, and can be estimated taking into account the partonic structure of the hadrons, which can be a source of photons. In what follows we will consider the contribution associated to elastic processes, where the incident hadron remains intact after the photon emission (For a related discussion about this subject see Refs. [, ]). For the proton case, a detailed derivation of the elastic photon distribution was presented in Ref. []. lthough an analytical expression for the elastic component is presented in Ref. [], it is common to found in the literature the study of photon - induced processes considering an approximated expression for the photon distribution of the proton proposed in Ref. [7], which can be obtained from the full expression by disregarding the contribution of the magnetic dipole moment and the corresponding magnetic form factor. s demonstrated in Ref. [8] the difference between the full and the approximated expressions is smaller than % at low-x. Consequently, in what follows we will use the expression proposed in Ref. [7], where the elastic photon distribution of the proton is given by γp el (x) = α ( ) [ em x + 0.x ln(ω) π x + Ω Ω + ] Ω, () where Ω = + (0.7GeV )/Q min and Q min (xm) /( x). On the other hand, the equivalent photon flux of a nuclei is assumed to be given by [9] γ el (x) = α em Z π [ η K0 ( η) K ( η) η U( η) ], () x where η = x s (R h + R hb )/γ L, R hi is the hadron radius and U( η) = K ( η) K0 ( η). One have that γ el (x) is enhanced by a factor Z in comparison to the proton one. Lets now discuss the modelling of the diffractive gluon distributions for the proton and nucleus. In order to describe the diffractive processes we will consider in what follows the Resolved Pomeron model [], which assumes that the diffractive parton distributions can be expressed in terms of parton distributions in the pomeron and a Regge parametrization of the flux factor describing the pomeron emission by the hadron. he parton distributions have evolution given by the DGLP evolution equations and should be determined from events with a rapidity gap or a intact hadron. In the proton case, the diffractive gluon distribution, gp D (x, µ ), is defined as a convolution of the pomeron flux emitted by the proton, f p IP (x IP ), and the gluon distribution in the pomeron, g IP (β, µ ), where β is the momentum fraction carried by the partons inside the pomeron. he pomeron flux is given by f p IP (x IP ) = tmax t min dt f IP /p (x IP, t) = tmax t min dt IP e BIP t x α IP (t) IP, () where t min, t max are kinematic boundaries. he pomeron flux factor is motivated by Regge theory, where the pomeron trajectory is assumed to be linear, α IP (t) = α IP (0) + α IP t, and the parameters B IP, α IP and their uncertainties are obtained from fits to H data [9]. he slope of the pomeron flux is B IP = GeV, the Regge trajectory of the pomeron is α P (t) = α P (0) + α P t with α P(0) =. ± and α P = GeV. he t integration boundaries are t max = m px IP /( x IP ) (m p denotes the proton mass) and t min = GeV. Finally, the normalization factor P =.7 is chosen such that x IP t max t min dt f IP /p (x IP, t) = at x IP = he diffractive gluon distribution of

4 the proton is then given by gp D (x, µ ) = dx IP dβ δ(x x IP β) f p IP (x IP ) g IP (β, µ ) = x dx IP x IP f p IP (x IP ) g IP ( ) x, µ. (7) x IP Similar definition can be established for the diffractive quark distributions. In our analysis we use the diffractive gluon distribution obtained by the H Collaboration at DESY-HER, denoted fit B in Ref. [9]. However, we checked that similar results are obtained using the fit. In order to specify the diffractive gluon distribution for a nucleus g D(x, µ ), we will follow the approach proposed in Ref. [] (See also Ref. [0]), which estimate g D taking into account the nuclear effects associated to the nuclear coherence and the leading twist nuclear shadowing. he basic assumption is that the pomeron - nucleus coupling is proportional to the mass number []. s the associated pomeron flux depends on the square of this coupling, this model predicts that when the pomerons are coherently emitted by the nucleus, f IP / is proportional to. Consequently, the nuclear diffractive gluon distribution can be expressed as follows (For details see Ref. []) g D (x, µ ) = R g x dx IP x IP [ ] ( ) x dt f IP /p (x IP, t) F(t) g IP, µ, (8) x IP where R g is the suppression factor associated to the nuclear shadowing and F (t) is the nuclear form factor. In what follows we will assume that R g = 0. as in Ref. [] and that F (t) e R t/, with R being the nuclear radius. One important open question in the treatment of photon and pomeron induced is if the cross sections for the associated processes are not somewhat modified by soft interactions which lead to an extra production of particles that destroy the rapidity gaps in the final state []. s these effects have nonperturbative nature, they are difficult to treat and its magnitude is strongly model dependent (For recent reviews see Refs. [, ]). In the case of IP IP interactions in pp/p p collisions, the experimental results obtained at EVRON [] and LHC [, 7] have demonstrated that one should take into account of these additional absorption effects that imply the violation of the QCD hard scattering factorization theorem for diffraction [8]. In general, these effects are parametrized in terms of a rapidity gap survival probability, S, which corresponds to the probability of the scattered proton not to dissociate due to the secondary interactions. Different approaches have been proposed to calculate these effects giving distinct predictions (See, e.g. Ref. [9]). n usual approach in the literature is the calculation of an average probability S and after to multiply the cross section by this value. s previous studies for the double diffractive production [9,, 8, 0 ] we also follow this simplified approach assuming S = 0.0 for the dijet production by IP IP interactions in pp collisions. It is important to emphasize that this choice is somewhat arbitrary, and mainly motivated by the possibility to compare our predictions with those obtained in other analysis. Recent studies from the CMS Collaboration [7] indicate that this factor can be larger than this value by a factor. he magnitude of S for IP IP interactions in p and collisions is still more uncertain [, ]. In what follows we will consider the approach proposed in Ref. [] for coherent double exchange processes in nuclear collisions. he basic idea in this approach is to express the IP IP cross section in the impact parameter space, which implies that the double pomeron exchange process becomes dependent on the magnitude of the geometrical overlap of the two nuclei during the collision. s a consequence, it is possible to take into account the centrality of the incident particles and estimate S by requiring that the colliding nuclei remain intact, which is equivalent to suppress the interactions at small impact parameters (b < R + R B ). In order to obtain predictions for S in p and collisions at LHC energies, we have updated and improved the model proposed in Ref. [] and obtained the values presented in able I. In ppendix we give a brief explanation of the model for S calculation. detailed discussion of the model will be presented in a separated publication. One have that the predicted values for S are larger than those obtained in Ref. [] using a Glauber approach and in Ref. [] assuming that the nuclear suppression factor is given by S = S pp /(. ). Consequently, our predictions for the dijet production by IP IP interactions in p/ collisions may be considered an upper bound for the magnitude of the cross sections. In the case of γγ and γip interactions, we will assume S =, motivated by results obtained e.g. in Refs. [7, 8], which verified that the recent LHC data for the exclusive vector meson production in photon induced interactions can be described without the inclusion of a normalization factor associated to absorption effects. However, it is important to emphasize that the magnitude of the rapidity gap survival probability in γip still is an open question. For example, in Ref. [9] the authors have estimated S for the exclusive photoproduction of J/Ψ in pp/p p collisions, obtaining that it is and depends on the rapidity of the vector meson (See also Refs. [, 0, ]). herefore, similarly to our IP IP predictions, the results for the dijet production by γγ and γip interactions also may be considered an upper bound.

5 O 0 Ca 08 P b p BLE I: Gap survival probability S for IP IP interactions in p and collisions at s =.0 ev. dσ/dη [pb] 0 η(jet) dσ/dϕ [pb] ϕ /π jj [pb] dσ/dp p (jet)[gev] γ IP p p IP p IP p γ p γ p FIG. : Differential cross sections as function of η(jet) (left), ϕ(jet) (center) and p (jet) (right) for the dijet production by γγ, γip and IP IP interactions in pp collisions. III. RESULS In what follows we present our results for the dijet production by photon photon, photon pomeron and pomeron pomeron interactions in pp, p and collisions at the LHC energy (For a similar analysis for the heavy quark production Refs. [8, ]). s discussed in the Introduction, these processes are characterized by two rapidity gaps and intact hadrons in the final state. he experimental separation of these events using the two rapidity gaps to tag the event is not an easy task at the LHC due to the non - negligible pile-up present in the normal runs. n alternative is the detection of the outgoing intact hadrons. Recently, the LS, CMS and OEM Collaborations have proposed the setup of forward detectors [ ], which will enhance the kinematic coverage for such investigations. Moreover, the LHCb experiment can study diffractive events by requiring forward regions void of particle production. < η < 8.0 []. In our analysis we will assume pp, ppb and PbPb collisions at a common center of mass energy ( s =.0 ev) in order to estimate their relative contributions as well as how the different channels of production are modified by increasing the atomic number. Moreover, we have reconstructed the jet using the anti k algorithm [] with distance parameter R = 0. as implemented in the Fastjet software package [7] and selected jets with p > 0 GeV and η <.0. In the case of γγ interactions we include the requirement that the impact parameter of the colliding hadrons should be larger than the sum of its radii. he cross sections for the partonic subprocesses are calculated at leading order in FPMC using HERWIG.. Finally, the subleading contribution for the diffractive interactions associated to the Reggeon exchange will be disregarded in our calculations. s demonstrated e.g. in Refs. [, ], such contribution can be similar to the Pomeron one in some regions of the phase space. However, as the description of the Reggeon exchange still is an open question, we postpone the analysis of its impact for the dijet production for a future publication. In Fig. we show the distributions of the pseudo rapidity η(jet) (left panel) and transverse momentum p (jet) (right panel) of the highest-p jet, and the azimuthal angular distance between the two highest-p jets ϕ jj. he predictions are presented separately for the dijet production by γγ, γip and IP IP interactions in pp collisions. Initially, lets analyze the η(jet) distribution. We have that the contribution of the IP IP process dominates at central pseudo rapidities, being a factor ( ) larger than the γip (γγ) one. However, the dijet production by γip interactions implies a broader pseudo rapidity distribution. s a consequence, this process becomes dominant for η(jet).8. In particular, in the kinematical region probed by the LHCb detector, the dijet production will be dominated by γip interactions. Consequently, the analysis of this process by the LHCb Collaboration can be an important test for the QCD treatment of the photoproduction of dijets in terms of the Resolved Pomeron model. On the other hand, the main contribution for the p (jet) distribution comes from the IP IP interactions, which is directly associated to the

6 dσ/dη [pb] 0 η(jet) dσ/dϕ [pb] ϕ /π jj [pb / GeV] dσ/dp p [GeV] IP Pb IP p γ IP Pb p γ IP p Pb γ γ Pb p FIG. : Differential cross sections as function of η(jet) (left), ϕ(jet) (center) and p (jet) (right) for the dijet production by γγ, γip and IP IP interactions in ppb collisions. he contributions for the γip interactions associated to a photon emitted by the nucleus ( γ P b IP p) and by the proton ( γ pip P b ) are presented separately. dσ/dη [pb] 0 η(jet) dσ/dϕ [pb] ϕ /π jj [pb / GeV] dσ/dp p [GeV] γ IP Pb Pb IP Pb IP Pb γ γ Pb Pb FIG. : Differential cross sections as function of η(jet) (left), ϕ(jet) (center) and p (jet) (right) for the dijet production by γγ, γip and IP IP interactions in PbPb collisions. dominance of central rapidities in the calculation of this distribution. Similarly, the IP IP interaction is dominant in the angular distribution of the dijets, with the main contribution being associated to back - to - back configurations. he results for the dijet production in ppb collisions are presented in Fig.. In this case we obtain asymmetric η(jet) distributions, which is expected since the photon and pomeron fluxes are different for a proton and a nucleus. In order to demonstrate it, we show separately in Fig. the γ P b IP p and γ p IP P b contributions, which are associated to a photon emitted by a proton and a nucleus, respectively. s in pp collisions, the IP IP contribution is dominant at central pseudo - rapidities and in the p range considered. Moreover, events are characterized by back - to - back configurations for the dijets. However, differently form the pp case, the IP IP contribution is larger than the γip one in all range of η(jet) considered. In particular, for η(jet), it dominates by a factor, which implies that analysis of the dijet production in this kinematical region can be useful to test the description of IP IP interactions in nuclear reactions. It is important to emphasize that this conclusion is not modified even if our prediction for S is reduced by two orders of magnitude, as predicted in alternative models for the calculation of the gap survival probability in nuclear collisions [, ]. In Fig. we present our predictions for the dijet production in PbPb collisions. In this case we have that the γip contribution dominates in all pseudo - rapidity and transverse momentum ranges considered. In particular, at central pseudo rapidities, we predict that the difference between the predictions is of the order of.his result is directly associated to the large suppression of the diffractive interactions in nuclear collisions due to the soft re-scattering processes that imply the dissociation of the incident nuclei and generate new particles that populate the rapidity gaps in the final state. s a consequence, we have a very small value for the the gap survival probability in PbPb collisions (See able I). lthough the γγ interactions in nuclear collisions are enhanced by a factor Z in comparison to pp case, our results indicate that this channel is only competitive for the dijet production at η(jet) 0. In order to estimate the atomic number dependence of the relative contribution between the IP IP and γip channels for the dijet production, in Fig. we present our predictions for the ratio between IP IP and γip distributions considering collisions and different values of. For comparison, the prediction for pp collisions also is presented. Our results indicate that the IP IP contribution increases at ligher nuclei and become dominant in the dijet production at central rapidities

7 7 Ratio Ratio Ratio proton O Ca Pb 0 η(jet) ϕ /π jj p [GeV] FIG. : Nuclear contributions on the cross section ratio of the IP IP to the γip channels as function of η(jet) and ϕ(jet). hree nuclear species are used: 08 Pb, 0 Ca, and O. in pp collisions. On the other hand, the γip channel is dominant in CaCa and PbPb collisions. In principle, this conclusion should not modified by more elaborated models for the calculation of S. s discussed before, S in γip interactions is expected to be of the order of the unity, while IP IP one the alternative models for S in nuclear reactions predict smaller values than that used in our analysis. herefore, we believe, that the analysis of the dijet production in nuclear collisions with heavy nuclei can be useful to study the photon - pomeron mechanism at high energies. Before to summarize our results, some comments are in order. In our calculations we have estimate the dijet production in γip interactions considering the leading order subprocesses present in the HERWIG. Monte Carlo. he contribution of the next - to - leading order (NLO) corrections for this process is large [8], being approximately a factor. he comparison of these predictions with the recent H and ZEUS data indicates that the NLO QCD calculations overestimate the data by approximately 0 0% (For a recent review see, e.g. Ref. [9]), with the origin of this suppression being a theme of intense debate (See Ref. [0] for a recent discussion). s a consequence, we believe that the leading order predictions are a reasonable first approximation for the dijet photoproduction. However, the inclusion of the NLO corrections and a suppression model for γip interactions is an important aspect that deserves a more detailed analysis in the future. nother important shortcoming in our study is associated to the fact that we only have considered the direct component of the photon for the dijet photoproduction, where a point like photon interacts with a parton from the Pomeron. In other words, we have disregarded the resolved component, where the photon behaves as a source of partons, which subsequently interacts with partons from the Pomeron. In principle, these two components can be separated by measuring the photon momentum participating in the production of the dijet system, denoted by x γ. heoretically, one expect the dominance of the direct (resolved) processes at high (low) values of x γ. Experimentally, this separation is not so simple due to hadronization and detector resolution and acceptances, but still feasible. herefore, our calculations for the dijet production in γip interactions are realistic for events with large values of x γ. However, as the resolved processes are predicted to be important at small p and large η(jet) [], it is possible to analyse the expected impact of the resolved contribution in our main conclusions. In the case of pp collisions (See Fig. ), the resolved processes should to increase the γip prediction for the pseudo - rapidity distribution in the region of large values of η(jet), where the γip interaction is dominant. Consequently, our conclusion that the production of dijets by γip interactions can be studied in pp collisions by the analysis of the large - η(jet) region is not expected to be modified by the inclusion of the resolved processes. Similarly, by the analysis from Fig., we have that the dominance in P bp b collisions of the γip interactions in the full η(jet) range should not be modified. Finally, in the case of pp b collisions, as the dijet production by IP IP interactions is a factor than the direct γip prediction (See Fig. ), we also do not expected that this dominance to be modified by the inclusion of the resolved γip contribution. herefore, we believe that our main conclusions must not be strongly modified if this process is included in the analysis. However, we also believe that the resolved contribution for the dijet production is an important aspect that deserves to be considered and we plan to include this contribution in the FPMC generator. IV. SUMMRY s a summary, in this paper we have presented a detailed analysis for the dijet production in pp/p/ collisions at the LHC. In particular, the comparison between the predictions for the dijet production by photon photon, photon pomeron and pomeron pomeron interactions was presented considering a common framework implemented in the Forward Physics Monte Carlo. We have generalized this Monte Carlo for nuclear reactions and performed a detailed

8 comparison between the IP IP, γip and γγ predictions for the dijet production in pp/ppb/pbpb collisions at s =.0 ev. For the pomeron - induced processes in pp collisions, we have considered the framework of the Resolved Pomeron model corrected for absorption effects, as used in the estimation of several other diffractive processes. In the case of nuclear collisions, we have generalized this model, following Refs. [, ]. Moreover, the absorption effects also have been included in our estimates for the dijet production by IP IP interactions in nuclear collisions. Our results indicate that in pp collisions the IP IP channel is dominant at central rapidities, being suppressed at forward rapidities. In particular, in the kinematical range probed by the LHCb detector, we predict that the main contribution for the dijet production comes from γip interactions. In the case of ppb collisions, the IP IP interactions are dominant. In contrast, our results indicated that in collisions with heavy nuclei, the dijet production by γip interactions is dominant, which indicates that this process can be used to test the Pomeron Resolved Model and its generalization for nuclei. Finally, our results indicate that the experimental analysis of the dijet production would help to constrain the underlying model for the pomeron and the absorption corrections, which are important open questions in Particle Physics. 8 cknowledgments Useful discussions with Marek asevsky are gratefully acknowledged. CPES, FPERJ and FPERGS, Brazil. his research was supported by CNPq, ppendix : Modelling S in the impact parameter space In order to estimate the survival gap probability in the impact parameter space, we will consider an approach similar to that proposed in Refs. [, ]. In this ppendix we present the basic aspects of this approach and postpone for a future publication a detailed discussion of the assumptions and uncertainties present in our calculations. Initially, lets to express Eq.() in terms of the Pomeron-Pomeron cross section σ(h h B h Xj j Y h B ) = dx dx B σ(ip IP B Xj j Y ), () where σ(ip IP B Xj j Y ) = x dx IP x IP x B dx IP B x IP B g IP / ( ) ( ) x, µ xb g IP /B, µ f IP / (x IP x IP x ) f IP /B (x IP B ) ˆσ(gg j j )(). IP B aking into account the transferred momentum q dependence of the partons emitted by the Pomeron, the above equation can be written in terms of the Pomeron flux in the momentum space f IP /i (x IP, q ) as follows σ(ip IP B Xj j Y ) = d dx ( ) ( ) IP b dx IP B x g IP /, µ xb g IP /B, µ x x IP x B x IP B x IP x IP B d q (π) eiq.b fip / (x IP, q ) f IP /B (x IP B, q ) ˆσ(gg j j ). () with f IP / (x IP, q ) = R g d k (π) F (k )F ((k q) ) IP (x IP, k ) IP (x IP, (k q) ), () and f IP /B (x IP, q ) = R g B d k (π) F B(k )F B ((k + q) ) IP (x IP, k ) IP (x IP, (k + q) ). () he functions F i and IP characterize the hadronic form factors and Pomeron propagators, respectively. In what follows we will assume that F i (k ) = / IP er i k /. ()

9 9 and IP (x IP, k ) = x ɛ IP e α log x IP k, (7) with the parameters being those obtained by the HER H experiment [9]. s both the form-factor and the propagator are Gaussians in k, the integrals over the Pomeron momentum k can be performed. Using the above forms in Eqs.(,) we write f IP / (x IP, q ) = R g IP x ɛ d k IP (π) e R k e R (k q), (8) and f IP /B (x IP, q ) = R g B IP x ɛ IP d k (π) e R B k e R B (k+q), (9) with R = R α log x IP, and R B = R B α log x IP. Using Eqs. (8), (9) and () in Eq. () and performing the integrations over k,k B and q variables we obtain where dσ d b (IP IP B Xj j Y, b) = π Q e Q b / σ(ip IP B Xj j Y ), () Q R +. () R B In order to calculate the integrated cross section taking into account the absorptive effects we multiply the above by the probability of not having strong interactions e Ω(s,b), where Ω(s, b) is the nuclear/proton opacity. In the pp case, we assume that the proton elastic profile can be describe by a Gaussian form, which implies that the opacity in proton-proton collisions is given by Ω pp (s, b) = σ tot(s) π B soft (s) e b /B soft (s), () where σ tot is the total pp cross section and B soft is the elastic scattering effective slope. We take these parametrizations from []. In order to derive a similar expression for the nuclear case, which is simple and can be used in analytical calculations, we have adjusted the Wood - Saxon distribution for the nuclei by a Gaussian one e Q 0 b /, with Q 0 being an effective parameter fitted to each nuclei. similar procedure was proposed in Ref. []. s a consequence, we obtain for collisions that Ω (s, b) = σ tot(s) π Q 0e Q 0 b /, () where is the atomic number and Q 0 is obtained from the nuclear form factor [7]. On the other hand, for p collisions we consider that the opacity can be expressed by Ω p (s, b) = σ tot(s) e b /B eff, () πb eff where B eff = + B Q soft 0. Using the above opacities, it is possible to calculate the integrated cross section for a general h a h b collision (a = p or ), which will be given by σ(ip a IP b Xj j Y ) = 0 d b dσ d b (IP aip b Xj j Y, b)e Ω(s,b) = S IP aip b σ(ip a IP b Xj j Y ) Ω=0. () For collisions, the suppression factor can be expressed by S IP IP (s) = Q 0 ( db e Q b / exp σ tot(s) π Q 0e Q 0 b / ), ()

10 O (s) > < S Ca (s) > p < S O Ca Pb Pb 0 8 s (ev) 0 8 s (ev) FIG. : Energy dependence of the survival gap factors S P P (s) and S PpP (s) for (left panel) and p (right panel) collisions considering different nuclei. with Q Q 0. s the above integral is of the type I = 0 dx e αx e λe βx = α 0 dy exp( λ y β/α ) = λ α/β β ( α ) γ β, λ, (7) where γ is the incomplete gamma function, α = Q /, λ = σtot π Q 0 and β = Q 0/, S IP IP (s) can be written as follows S IP IP (s) = Q Q 0 ( σ ) tot Q /Q 0 π Q 0 γ ( Q Similarly, we can obtain the suppression factor for p collisions, which is given by Q 0, σ ) tot π Q 0. (8) S IP pip (s) = Q ( eff B eff σ ) Q tot B eff ( γ πb Q B eff, σ ) tot, (9) eff πb eff where Q eff = /(R + B IP ). Finally, the expression for pp collisions is given by S IP pip p (s) = B ( soft σtot ) Bsoft/BIP ( Bsoft σ ) tot γ,, (0) B IP πb soft B IP πb soft which is similar to the expression derived in Ref. [] using a distinct approach. Using the above expressions we have derived the values for the survival gap probabilities for Pomeron Pomeron interactions in P bp b and pp b collisions presented in able I. In Fig. we present our predictions for its energy dependence considering and p collisions and different nuclei. It is important to emphasize that our prediction for pp collisions at the LHC energy is 0.0, as used in several phenomenological analysis in the literature [9,, 8, 0 ]. [] K. kiba et al. [LHC Forward Physics Working Group], J. Phys. G, 0 (0) []. Hebecker, Phys. Rept., (000); L. Schoeffel, Prog. Part. Nucl. Phys., 9 (0); M. G. lbrow,. D. Coughlin and J. R. Forshaw, Prog. Part. Nucl. Phys., 9 (0) [] P. D. B. Collins, n Introduction to Regge theory and high energy physics (Cambridge University Press, Cambridge, England, 977). [] V.. Khoze,. D. Martin and M. G. Ryskin, Eur. Phys. J. C, (00) [] R. J. M. Covolan and M. S. Soares, Phys. Rev. D 7, 070 (00) [] J. R. Cudell,. Dechambre, O. F. Hernandez and I. P. Ivanov, Eur. Phys. J. C, 9 (009) [7] C. Marquet, C. Royon, M. rzebiski and R. lebk, Phys. Rev. D 87, 00 (0)

11 [8] R. Maciula, R. Pasechnik and. Szczurek, Phys. Rev. D 8, (0); Phys. Rev. D 8, (0). [9] C. Brenner Mariotto and V. P. Goncalves, Phys. Rev. D 88, no. 7, 070 (0) [] C. Marquet, C. Royon, M. Saimpert and D. Werder, Phys. Rev. D 88, no. 7, 0709 (0) []. K. Kohara and C. Marquet, Phys. Lett. B 77, 9 (0) [] C. O. Rasmussen and. Sjstrand, JHEP 0, (0) [] C. Marquet, D. E. Martins,. V. Pereira, M. Rangel and C. Royon, Phys. Lett. B 7, (07) [] M. Luszczak, R. Maciula,. Szczurek and I. Babiarz, JHEP 70, 089 (07). [] V. Guzey and M. Klasen, JHEP 0, 8 (0) [] G. Ingelman and P.E. Schlein, Phys. Lett. B, (98). [7] M. Boonekamp,. Dechambre, V. Juranek, O. Kepka, M. Rangel, C. Royon and R. Staszewski, arxiv:. [hep-ph]. [8] V. P. Goncalves, C. Potterat and M. S. Rangel, Phys. Rev. D 9, no., 008 (0) [9] C.. Bertulani and G. Baur, Phys. Rep., 99 (988); G. Baur, K. Hencken, D. rautmann, S. Sadovsky, Y. Kharlov, Phys. Rep., 9 (00); V. P. Goncalves and M. V.. Machado, Mod. Phys. Lett. 9, (00); C.. Bertulani, S. R. Klein and J. Nystrand, nn. Rev. Nucl. Part. Sci., 7 (00); K. Hencken et al., Phys. Rept. 8, (008). [0] V. M. Budnev, I. F. Ginzburg, G. V. Meledin and V. G. Serbo, Phys. Rept., 8 (97). [] E. Fermi, Z. Phys. 9, (9). [] E. J. Williams, Phys. Rev., 79 (9). [] C. F. von Weizsacker, Z. Phys. 88, (9). [] B.. Kniehl, Phys. Lett. B, 7 (99). [] V. P. Goncalves and G. G. da Silveira, Phys. Rev. D 9, no., 00 (0) [] G. G. da Silveira and V. P. Goncalves, Phys. Rev. D 9, no., 0 (0) [7] M. Drees and D. Zeppenfeld, Phys. Rev. D 9, (989). [8] V. P. Goncalves, D.. da Silva and W. K. Sauter, Phys. Rev. C 87, 080 (0) [9] H Collab.,. ktas et al., Eur. Phys. J. C8, 7 (00). [0] L. Frankfurt, V. Guzey and M. Strikman, Phys. Rept., (0) [] B. Muller and. J. Schramm, Nucl. Phys., 77 (99). [] J. D. Bjorken, Phys. Rev. D 7, (99). [] V.. Khoze,. D. Martin and M. G. Ryskin, Int. J. Mod. Phys. 0 (0) no.08, 0 [] E. Gotsman, E. Levin and U. Maor, Int. J. Mod. Phys. 0, no. 08, 0 (0) []. ffolder et al. [CDF Collaboration], Phys. Rev. Lett. 8, (000):. altonen et al. [CDF Collaboration], Phys. Rev. D 8, 0009 (0) [] G. ad et al. [LS Collaboration], Phys. Lett. B 7, (0) [7] S. Chatrchyan et al. [CMS Collaboration], Phys. Rev. D 87, no., 00 (0) [8] J. C. Collins, Phys. Rev. D 7, 0 (998) Erratum: [Phys. Rev. D, 0990 (000)] [9] M. G. Ryskin,. D. Martin, V.. Khoze and. G. Shuvaev, J. Phys. G, 0900 (009) [0] M.V.. Machado, Phys. Rev. D 7, 000 (007); M. B. Gay Ducati, M. M. Machado, M. V.. Machado, Phys. Rev. D8, 00 (0); M. B. Gay Ducati, M. M. Machado and M. V.. Machado, Phys. Rev. C 8, 090 (0) [] M. Luszczak, R. Maciula and. Szczurek, Phys. Rev. D 8, 8 (0). [] M. Luszczak, R. Maciula and. Szczurek, Phys. Rev. D 9, no., 00 (0). [] C. Brenner Mariotto and V. P. Goncalves, Phys. Rev. D 9, no., 0 (0). []. B. Kaidalov, V.. Khoze,. D. Martin and M. G. Ryskin, cta Phys. Polon. B, (00) [] E. Levin and J. Miller, arxiv:080.9 [hep-ph]. [] V. P. Goncalves and W. K. Sauter, Phys. Rev. D 8, 0009 (0) [7] V. P. Goncalves, B. D. Moreira and F. S. Navarra, Phys. Rev. C 90, no., 0 (0); Phys. Lett. B 7, 7 (0). [8] V. P. Goncalves, B. D. Moreira and F. S. Navarra, Phys. Rev. D 9, no., 00 (07) [9] W. Schafer and. Szczurek, Phys. Rev. D 7, 090 (007) [0] V. Guzey and M. Zhalov, JHEP, 07 (0); JHEP, 0 (0). [] S. P. Jones,. D. Martin, M. G. Ryskin and. eubner, JHEP, 08 (0) [] M. G. lbrow et al. [FP0 R and D Collaborations], JINS, 00 (009) [] he CMS and OEM Collaborations, CMS-OEM Precision Proton Spectrometer echnical Design Report, http: //cds.cern.ch/record/779. [] M. asevsky [LS Collaboration], IP Conf. Proc., (0). [] LHCb Collaboration, CERN-LHCb-CONF [] M. Cacciari, G. P. Salam and G. Soyez, JHEP 080, 0 (008) [7] M. Cacciari, G. P. Salam and G. Soyez, Eur. Phys. J. C 7, 89 (0) [8] M. Klasen,. Kleinwort and G. Kramer, Eur. Phys. J. direct, no., (998) [9]. Valkarova, Int. J. Mod. Phys. 0, no. 08, 0 (0). [0] V. Guzey and M. Klasen, Eur. Phys. J. C 7, no. 8, 7 (0) [] J. Rau, B. Muller, W. Greiner, G. Soff, J. Phys. G, (990) []. K. Kohara, E. Ferreira,. Kodama, Eur. Phys. J. C 7, 7 (0). [] E. Gotsman, H. Kowalski, E. Levin et al. Eur. Phys. J. C 7, (00). [] M. Klasen, Rev. Mod. Phys. 7, (00) [] M. Drees, J. Ellis and D. Zeppenfeld, Phys. Lett. B, (989).