Sensitivity of the SHiP experiment to dark photons decaying to a pair of charged particles

Dark photons are hypothetical massive vector particles that could mix with ordinary photons. The simplest theoretical model is fully characterised by only two parameters: the mass of the dark photon m$_{\gamma^{\mathrm{D}}}$ and its mixing parameter with the photon, $\varepsilon$. The sensitivity of the SHiP detector is reviewed for dark photons in the mass range between 0.002 and 10 GeV. Different production mechanisms are simulated, with the dark photons decaying to pairs of visible fermions, including both leptons and quarks. Exclusion contours are presented and compared with those of past experiments. The SHiP detector is expected to have a unique sensitivity for m$_{\gamma^{\mathrm{D}}}$ ranging between 0.8 and 3.3$^{+0.2}_{-0.5}$ GeV, and $\varepsilon^2$ ranging between $10^{-11}$ and $10^{-17}$.

in the minimal dark photon model with decays to charged particles is given in Section 4 for the three production modes studied. Finally Section 5 provides a conclusion.

The SHiP detector and simulation
SHiP [15] is a new general purpose fixed-target experiment intended to exploit the proposed BDF to search for particles present in hidden portal models. The 400 GeV proton beam extracted from the SPS will be dumped on a high density target with the aim of accumulating 2 × 10 20 protons on target during 5 years of operation. A dedicated detector, based on a long vacuum tank followed by a spectrometer and particle identification detectors, will allow probing a variety of models with light long-lived exotic particles and masses below O(10) GeV. A critical component of SHiP is the muon shield, which deflects the high flux of muons produced in the target [18,19], that would represent a serious background in the search for hidden-sector particles. The detector is designed to fully reconstruct the exclusive decays of hidden particles and to reject the background down to below 0.1 events in the sample of 2 × 10 20 protons on target [20].
The detector consists of a large magnetic spectrometer located downstream of a 50 m-long and up to 5×11 m-wide cone-shaped decay volume [21]. To suppress the background from neutrinos interacting in the fiducial volume, the decay volume is maintained under a vacuum. The spectrometer tracker is designed to accurately reconstruct the decay vertex, mass and impact parameter of the decaying particle. A set of calorimeters followed by muon chambers provide identification of electrons, photons, muons and charged hadrons. A dedicated timing detector measures the coincidence of the decay products, which allows the rejection of combinatorial backgrounds. The decay volume is surrounded by background taggers to tag neutrino and muon inelastic scattering in the surrounding structures, which may produce long-lived SM V 0 particles, such as K L , that have topologies similar to the expected signals.
The spectrometer tracker is a crucial component in the reconstruction of the charged particles produced by the decay of dark photons. The baseline layout consists of four tracking stations (T1 to T4) symmetrically arranged around a dipole magnet as shown in Fig. 1. The transverse size of the tracker stations matches the size of the magnet. Each station consists of 9072 straw tubes which are arranged in four views (Y-U-V-Y), as shown in Fig. 2. The Y view has straws horizontally aligned. The U and V views are rotated by an angle of θ stereo = ±5 o . The x coordinate is hence measured with an accuracy of 1/sin(θ stereo ), directly impacting the measurement of the decay vertex, of the opening angle of the daughter particles (which enters the invariant mass), and of the impact parameter at the production target. In order to provide good spatial resolution and minimise the contribution from multiple scattering, the straw tubes are made of thin polyethylene terephthalate (PET). More detail about the initial design of the straw detector can be found in Refs. [20,22]. The pattern recognition algorithms applied to the hits on the straw spectrometer are described in Ref. [23], and the algorithms for particle identification are presented in Ref. [24].
In the simulation, proton fixed-target collisions are generated by PYTHIA 8.2 [25], inelastic neutrino interactions by GENIE [26] and inelastic muon interactions by PYTHIA 6 [27]. The heavy-flavour cascade production is also taken into account [28]. The SHiP detector response is simulated in the GEANT4 [29] framework. The simulation is done within FairShip, which is based on the FairRoot framework [30].

Dark photon production and decay
The minimal dark photon model contains an additional U(1) gauge group to the SM, A µ , whose vector gauge boson is called the dark photon γ D . A kinetic mixing term between the dark photon with field strength F µν and the SM U(1) gauge bosons with field strength F µν Y is allowed [8], with a reduced strength parameterised by a coupling ε, also called the kinetic mixing parameter. The corresponding   terms in the Lagrangian can hence be written as: In its simplest form, the knowledge of the mass of the dark photon m γ D and the kinetic mixing parameter ε is enough to characterise the model and calculate production cross section and decay properties.
Three different mechanisms are possible for the production of such new particles at a fixed-target experiment. All of them are studied in this paper.
The initial 400 GeV proton beam interacts with the nucleons from the target material, producing mesons. For meson decay channels involving photons, the photon can mix with the dark photon, as described in Section 3.1. This mode is opened only for dark photon masses below 0.9 GeV, as for mesons with masses above this threshold the main decay channels do not involve photons anymore.
The proton-nucleon interaction could also lead to the radiation of a dark photon via a bremsstrahlung process, as described in Section 3.2. This mode is heavily suppressed when the dark photon mass exceeds that of the proton, and remains relevant only for masses below 2 GeV.
The third production mode is via a Drell-Yan like process in Quantum Chromodynamic (QCD), i.e. quark-antiquark annihilation into the dark photon, as described in Section 3.3. This process is relevant for dark photon masses in the range O(1-10) GeV. Using the parton model with a factorisation scale below the GeV scale cannot give sensible results, as expected from the range of validity of parton distribution functions, and hence this region of the parameter space has not been considered for this production mechanism.
In this paper, the assumption is made that only the initial proton interacts. In reality, similar interactions could also happen with protons or mesons coming from cascade decays happening in the target material. For electromagnetic processes (electron bremsstrahlung of photons mixing with the dark photon), it has been shown in Ref. [31] that their contribution is negligible compared to the main production mechanisms described above. The study however remains to be done for hadronic interactions in the cascade decays, and will be the subject of future work. Hence the expectations presented here are conservative and the sensitivity could be improved in the future when this contribution is added.
As a final state, the dark photon decay to pairs of leptons or quarks as described in Section 3.4 is considered.

Production in meson decay
The PYTHIA 8.2 [25] Monte Carlo (MC) generator is used to produce inclusive QCD events in protonproton (p-p) collisions, through all available non-diffractive processes. Diffractive processes are less important in meson production, expected to increase the number of mesons produced by about 15% according to PYTHIA simulation. Because the diffractive processes also suffer from larger theoretical uncertainties, they have not been considered. The leading-order (LO) NNPDF2.3 PDF set [32] has been used with the default Monash 2013 tune [33], and the strong coupling constant set to α s = 0.13. One proton beam momentum is set to 400 GeV and the other to 0 (protons or neutrons from the fixed-target material). The mesons that are produced are then used as sources of dark photons, if they have decay channels to photons and their decay to a dark photon of mass m γ D is kinematically allowed. Four processes are found dominant (with other contributions neglected) and shown in Table 1. The decay tables of these four mesons are reset to having only one decay channel allowed with 100% branching ratio (π 0 → γγ, η → γγ, ω → π 0 γ, η → γγ). All relevant processes are then added together.
The branching ratios of the mesons to these new decay channels are functions of the m γ D , the kinetic mixing parameter ε, the meson type, pseudo-scalar or vector, and the meson mass [31,34,35]. For pseudo-scalar mesons P (π 0 , η 0 and η ), the branching ratio to γ D γ is given by: For vector mesons V (ω), the branching ratio to a γ D and a pseudo-scalar meson P is given by: For the branching ratios of the mesons to γγ or γπ 0 , the same values as implemented in PYTHIA 8.2 are used. The average number of mesons produced per pp interaction, n meson / pp, is shown for each meson type in the last column of Table 1, from non-diffractive pp collisions simulated with PYTHIA 8.2 [25], with its associated statistical uncertainty.
The cross section for the production of dark photons via meson decays produced in non-diffractive primary interactions of the proton beam is then computed as: using Eqs. (2) and (3) and values reported in Table 1. The Θ(m meson − m γ D ) factor is a step function ensuring that only the mesons in the accessible mass range are considered. To take into account the fact that the nucleon is bound in the target, and not free as assumed by PYTHIA in our simulation of this process, the total normalisation is taken using the inelastic proton-nucleon cross section corresponding to the SHiP target, σ inel SHiP (see Section 4). The cross section is proportional to ε 2 , from the dependency of Br(meson → γ D + X) in Eqs. (2) and (3).

Production in proton bremsstrahlung
In analogy with ordinary photon bremsstrahlung of scattering protons, the same process is used for dark photon production by scattering of the incoming 400 GeV proton beam on the target protons. Following Refs. [31,36], the differential γ D production rate can be expressed as: where σ pp (s/s ) are the total proton-proton cross sections evaluated for the incoming/outgoing proton energy scales, m p is the proton mass (set to m p = 0.938272 GeV [37]), P and E p are the proton beam initial momentum and energy respectively, p and E γ D are the momentum and energy of the generated dark photon respectively, p ⊥ and p are the components of the γ D momentum orthogonal and parallel to the direction of the incoming proton respectively, z is the fraction of the proton momentum carried away by the dark photon in the beam direction, α QED is the fine structure constant of Quantum Electro Dynamic (QED), set to 1/137, In this formulation, the nuclear effects from having bound rather than free protons in the target material cancel in the ratio σpp(s ) σpp(s) . However, the above formula does not take into account possible QCD contributions when the mass of the emitted γ D exceeds that of the proton, and the bremsstrahlung process starts to depend on the internal partons. It does not take into account the possibility of enhancement in the cross section due to nuclear resonances in the so-called vector meson dominance (VMD) model either. In consequence, two independent approaches are followed, leading to two different estimates of the final cross section.
In the first approach, when the mass of the dark photon is larger than 1 GeV, the standard dipole form factor [38] is included in the proton-γ D vertex, leading to a penalty factor that models the strong suppression of the bremsstrahlung production: According to Ref. [31], this form factor is conservative and probably underestimates the rates. The direct parton-parton QCD production will dominate above 1.5 GeV and is described in Section 3.3.
In the second approach, the VMD form factor from Refs. [39,40] is used, leading to an enhancement of the cross section by a factor 10 4 around the ρ and ω meson mass of 0.8 GeV, and still up to a factor 10 in the tail due to also considering resonances of masses 1.25 and 1.45 GeV following the description in Ref. [40].
Reformulating Eq. (5) as a function of the γ D angle θ to the beam line and its total momentum p, a two-dimensional normalised probability density function (PDF) f (p, θ) is extracted, and shown in Fig. 3 for two representative choices of m γ D . Note that due to the simple dependency of the production rate scaling as ε 2 , the normalised PDF is independent of ε. The dark photons are generated with maximum probability on each side of the beam axis (θ close to 0) with a factor of 5 more chance to have p < 100 GeV compared to p > 200 GeV, for the low masses, and increased probability to have high momentum as the mass increases. Events are generated using a PYTHIA 8 particle gun with the γ D as particle, randomly choosing the γ D (p, θ) values according to the normalised 2D PDF f (p, θ), extracted for each m γ D point studied.
The integral of d 2 N dpdθ × FF, with FF the penalty dipole form factor or the VMD form factor, in the range of momenta and solid angle kinematically allowed, provides an estimate of the dark photon production rate per pp interaction through proton bremsstrahlung, scaling as ε 2 . The production cross sections using the dipole form factor and VMD form factor methods are expressed by: and shown in Fig. 4. The conditions of validity of the approximation used to derive Eq. (5) [42,43] require a lower momentum bound for the γ D at p min = 0.1P p [36], and an upper bound at p max = 0.9P p , as well as an upper bound on p ⊥ < 4 GeV, giving θ max 0.1 rad.

Drell-Yan production
For production of the dark photon in parton-parton scattering, the generic implementation of a resonance that couples both to SM fermion pairs and hidden particles is used, as implemented in PYTHIA 8.2 under the "HiddenValley" Z model [44]. A cross-check has been done that similar kinematic distributions for the dark photons are found using another Z implementation in PYTHIA from the "New Gauge Bosons" class of processes [45].
The dark photons are generated in the mass range 1.4 < m γ D < 10 GeV. Below 1.4 GeV one leaves the domain of perturbative QCD and the parton model cannot be used anymore.
The LO cross section given by PYTHIA when the new particle has the properties of the dark photon is shown in Fig. 5. The nuclear effects are neglected and the parton-parton cross section from PYTHIA is used without modification. Like for the meson and proton bremsstrahlung processes, it is found to scale as ε 2 . The LO NNPDF2.3 PDF set [32] has been used with the default Monash 2013 tune [33], and the strong coupling constant set to α s = 0.13. An empirical function is extracted to parameterise the cross section (in mb) as a function of the γ D mass (in GeV) in a continuous way, described in Eq. (9). The impact of several sources of theoretical uncertainties (PDF choice [46], QCD scales, α s ) are studied and shown in Fig. 5 (see also Section 4.4). The impact of nuclear effects is checked using the nuclear modification factors available in PYTHIA, with the most recent nuclear PDF set EPPS16 [47], using the two atomic masses (A = 84 and 117) available around the SHiP target material one (A = 96). Both give very similar results, with a cross section varying within ±6% from the NNPDF2.3 proton PDF one, depending on the γ D mass. The alternative generator Madgraph5_aMC@NLO v2.7.2 [48] is also used to cross check the cross section calculation. The parameterisation of the width of the resonance, dependent on the branching ratios to fermion pairs, is a little different and explains the difference seen.
Higher-order contributions to the process could lead to a sizable increase of the cross section at such low masses. Using the MATRIX v1.0.5 program [49][50][51][52][53][54], the ratio of next-to-next-to-leading-order (NNLO) over LO differential cross sections for standard Drell-Yan production at √ s = 27.43 GeV as a function of the dilepton invariant mass M is found to be rather flat at 1.7 ± 0.17 (stat) for M between 1.4 and 5 GeV, increasing up to 2.2 at M = 10 GeV. The QCD scale uncertainties on the ratio are found to be +20% −11% in the range 2 <M < 5 GeV increasing above ±30% for M above 10 GeV or below 2 GeV. The MSTW2008 NNLO PDF set [55] has been used for all calculations. These calculations are also found to be in fair agreement (within 10%) with a study performed more specifically in protonantiproton collisions, with a special interest for low dilepton masses and in particular for √ s = 30 GeV, from Ref. [56], and considering soft-gluon resummations at all orders. Given the lack of experimental data at these low masses and low √ s to confirm the size of the expected correction and the impact from PDF and non-perturbative effects on the actual dark photon production, a final k-factor of 1.7 ± 0.7 is applied to the LO PYTHIA cross section from Eq. (9).
The relative contribution from each process is shown in Fig. 6, as a function of m γ D , for the three production modes, in the two scenarios considered for the proton bremsstrahlung mode.

Dark photon decays
Except for the meson production mode, in which the new particle couples to the parent meson via mixing with the photon and hence cannot be a resonance from PYTHIA's point-of-view, in QCD and proton bremsstrahlung the γ D is implemented as a resonance. In all cases, the decay channels are implemented as follows.
The partial decay width of the dark photon into a lepton pair is given by [36]:   where m is the lepton mass, for electron, muon or tau leptons, if kinematically allowed. Following the approach used by the authors of Ref. [61], the partial decay width into quark pairs is computed as: where is the energy-dependent R-ratio quantifying the hadronic annihilation in e + e − collisions [62], tabulated from 0.3 to 10.29 GeV.
The lifetime of the γ D is then naturally set to the inverse of its total width, summing all the kinematically-allowed channels for calculating the total width. It is proportional to 1/ε 2 . The branching ratios to individual channels are set to the ratio of the partial over total width, and are hence independent of ε. For separating the hadronic channels into the different quark-flavoured pairs allowed kinematically, the coupling is assumed to be proportional to the quark charge q as n C × q 2 [63], with n C = 3 the number of coloured charges. When the γ D is implemented as a resonance in PYTHIA, the decay goes explicitly through the pair of quarks, before hadronisation. Otherwise the hadrons are found as direct decay products of the γ D .
The branching ratio of the γ D into pairs of leptons or quarks is shown in Fig. 7 as a function of m γ D . The hadronic decays become available above the pion mass threshold. The expected lifetime of the γ D as a function of its mass and ε mixing parameter is shown in Fig. 8.

SHiP sensitivity
In order to maximise the statistical power of the limited number of events produced with PYTHIA in the different production modes, the γ D decay vertex position is randomly assigned to be inside the decay vessel of length L Vessel = 50.760 m, and the associated probability of this happening is given as a function of the γ D four-momentum (p, E γ D ) and lifetime cτ : with γ = E γ D / E 2 γ D − p 2 , β = p/E γ D , L 0 is the distance from the target to the entrance of the decay vessel, and is randomly distributed between 0 and L Vessel with a flat prior.
The total event rate expected is then extracted from the cross sections σ prod defined in Section 3, i.e. Eqs. (4), (8) and (9) for the meson, proton bremsstrahlung and QCD productions, respectively, normalising to the luminosity expected from the N = 2 × 10 20 proton-on-target events that will be collected by the end of the SHiP physics program. The expected rate is taking into account the detector acceptance and the efficiency to reconstruct the decay products in the SHiP detector, P vessel and P reco described in detail in Sections 4.2 and 4.3, respectively, and following Eq. (14): The SHiP luminosity is defined as: , using an inelastic proton-nucleon cross section of σ inel SHiP = 10.7 mb [15], which directly corresponds to the SHiP target material (Molybdenum) nuclear interaction length and density.
The strategy of the analysis relies on identifying the decays of the γ D into at least two charged particles, γ D → ch + ch. The reconstructed charged tracks must originate from a common vertex. These requirements are enough to ensure that almost no background event will survive the selection, as demonstrated in Refs. [5,15]. The 90% confidence level (CL) limits on the existence of a γ D with given (m γ D , ε) are hence set by excluding regions where more than N γ D = 2.3 events are expected.

Decay channels
The following final states are considered, whenever available for a given m γ D : e + e − , µ + µ − , τ + τ − , and any hadronic decay channels leading to charged particles (e.g. π + π − + X, K + K − + X). The branching ratio to the different final states is shown in Fig. 9 for all the simulated (m γ D , ε) points in the three different production modes, as a function of m γ D , calculating the mean value over the different ε samples. All events classified under "e + e − ", "µ + µ − ", "τ + τ − " and "charged hadrons" have at least two charged particles, their sum is represented as "ch + ch". Only the events classified under "neutral hadrons" are lost due to the analysis selection described in Section 4.3. Compared to Fig. 7, Fig. 9 highlights the mass scan actually simulated, and the separation of the hadronic final states into the charged and neutral ones.

Vessel acceptance
For events which have two charged particles, the γ D decay vertex is further required to be within the vessel volume. The efficiency of this selection, P vessel is defined as the ratio of the sum of the weights w vtx ( ) of events passing the vertex selection described in Table 2 over the total number of events with a dark photon decaying to at least two charged particles. This efficiency is shown in Fig. 10 as a function of (m γ D , ε), for the three production modes. It is mostly driven by the lifetime of the γ D , and the kinematics of the γ D produced in the target. Its maximum is around 5% for the production via meson decay, 10% for the proton bremsstrahlung production, and for higher masses in QCD production.

Reconstruction of the decay products
The strategy employed in this analysis relies uniquely on the reconstruction of charged particles by the SHiP straw tracker. Future extensions of this work could consider also calorimeter deposits (with the possibility to fully reconstruct π 0 decays to two photons) and muon detectors. Events are retained if two tracks are found passing the criteria summarised in Table 2, namely that the two tracks are within the fiducial area of the detector up to the fourth layer after the magnet, the fit converged with good quality requirements (χ 2 /NDF < 5 with NDF the number of degrees of freedom of the fit). The tracks are required to have an impact parameter (IP) less than 0.1 m in the (x,y) plane, a momentum p above 1 GeV, and a distance of closest approach (DOCA) below 1 cm. Criteria on the number of hits (NDF > 25) or presence of hits before/after the magnet are meant to reduce backgrounds which could come from   Fig. 10: Efficiency of requiring the γ D decay vertex to be inside the decay vessel volume, for the three production modes studied. An interpolation between the simulated (m γ D , ε) values is performed using a linear interpolation via Delaunay triangulation.
particles re-entering the detector volume due to the magnetic field. At the moment, the resolution of the timing detector is neglected, and MC truth information is used instead.
The efficiency of having two good tracks passing the selection for events which had two charged particles and γ D vertex in the decay volume, P reco is shown in Fig. 11. Once the γ D decays in the vessel volume, the reconstruction efficiency is above 80% in most of the parameter space. For production via meson decay, a dependency on ε is observed, with the efficiency dropping to below 50% as ε decreases. This is found to be related to the wider angular distribution of dark photons produced in meson decays, introducing a dependency on the position of the decay vertex. x-y within vessel volume and at least 5 cm away from its inner walls. Straw tracker hits in each layer -before and after magnet -up to tracking station 4 Tracks ≥ 2 tracks NDF > 25, χ 2 /NDF < 5, DOCA < 1 cm, p> 1 GeV, IP< 0.1 m

Systematic uncertainties
The following sources of systematic uncertainties from theory are investigated, for the three production modes. The missing contributions from cascade decays will be the subject of future work and is not considered.
For the meson production, the overall rate is affected by the following uncertainties: -Branching ratios of the mesons to decay channels with photons from Table 1: from Ref. [37], the uncertainties on the measurement of these branching ratios are 0.03%, 0.5%, 3.4% and 3.6% for π 0 → γγ, η 0 → γγ, ω → π 0 γ and η → γγ respectively, translating directly to the final rate. -Uncertainty on the meson multiplicities and shape of their kinematics properties: PYTHIA 8.2 has been compared with data in several existing publications. In Ref. [64], a comparison to NA27 and NA56 data is made for the inclusive production of π 0 mesons, and reasonable agreement is found, within 30% in the kinematic regions targeted by our measurement. PHENIX and ALICE also measured inclusive π 0 , η and ω production and ratios [65][66][67], and showed global agreement within about 20% with the PYTHIA 8 (Monash 2013 Tune) simulation.
Adding the different sources in quadrature, this results in a total systematic uncertainty of ±30%.
For the proton bremsstrahlung, the theory systematic uncertainties concern: -uncertainties on the inelastic p-p cross section σ pp (s), which will mostly cancel in the ratio σpp(s ) σpp(s) , are neglected. -Dipole form factor versus VMD form factor: the two scenarios are presented separately in the final exclusion limits. -Contribution from protons undergoing elastic scattering before radiating the γ D : an upper bound is derived using a factor 1 1−P el = 1.34 [68], with P el the probability for an incoming proton to generate an elastic scattering, P el = ±2 GeV, the total rate is changed by +15% −30% . Varying the lower and upper bounds p min (p max ) by ±0.04 (∓0.04), the total rate is changed by +40% −25% .
A total systematic uncertainty of +50% −40% is assumed to cover these sources. For the QCD production, the theory systematic uncertainties concern the parameterisation of the LO cross section, the choice of NNLO k-factor and the impact from QCD scales and PDFs. Figure 5 shows the relative contributions from QCD scales and PDF on PYTHIA's LO cross section. The choice of PDF set is giving large variations in normalisation, but not affecting the overall shape of the cross section versus mass. The PDF set chosen is conservatively the one giving the lowest cross section. As discussed in Section 3.3, in the end the uncertainty is dominated by the NNLO k-factor of 1.7 ± 0.7. The total systematic uncertainty is hence taken as ±40%.
Experimental systematic uncertainties concern the measurement of the luminosity, the modeling of the tracking efficiency and the assumptions entering the 0-background estimate. They have been neglected in this study, as they are expected to be small compared to the theoretical uncertainties.

Extraction of the limit
Events are generated following a discrete grid in (m γ D ,ε) values, and passed through the full simulation of the SHiP detector and reconstruction algorithms. The γ D mass is varied between the electron-pair production threshold, and 10 GeV, in 0.001 to 1 GeV steps. The kinetic mixing parameter ε is varied between 10 −4 and 10 −9 in varying-size steps in log(ε).
To find the ε values that allow to reach 2.3 expected events, the expected rate is studied as a function of ε for the discrete mass points, with a linear interpolation between fully-simulated values. Between mass points, a linear interpolation is also performed. The rate of events is driven by two aspects. For large ε values, larger cross sections are expected but the detection efficiency decreases rapidly due to small lifetimes and decays happening before the decay vessel. As ε decreases, the cross section decreases as ε 2 but the events have more and more probability to reach the vessel and the rate increases, up to a turning point where the decay vertex happens after the decay vessel and/or the cross section becomes too small. Hence the 90% CL exclusion region is contained inside a lower and upper limits on ε 2 for each mass point. The dependency of the excluded region on the mass is driven by the kinematic properties of the γ D and its decay products, affecting the detector acceptance and selection efficiency.
As shown in Fig. 12 for representative mass points, for all processes, the upper bounds have little dependency on the absolute normalisation of the rate (so in particular systematic uncertainties on the cross sections and other quantities affecting the overall rate), due to the very steep dependency of the rate as a function of ε. The lower bounds are however more sensitive.
The SHiP experiment is expected to have a unique sensitivity in the mass region m γ D ranging between 0.8 and 3.3 +0.2 −0.5 GeV, and ε 2 ranging between 10 −11 and 10 −17 .

Conclusion
The sensitivity of the SHiP detector has been investigated for the simplest vector portal model, in which the only hidden-sector particle connecting to SM particles is a dark photon. The model is fully parameterised by only two parameters, the mass of the dark photon m γ D and the kinetic mixing parameter ε. Three different production mechanisms have been investigated, namely the production via meson decays from non-diffractive proton-nucleon interactions, by proton bremsstrahlung and by QCD parton-parton interaction. Different sources of systematic uncertainties have been considered, dominated by theory predictions on the cross section times branching ratios (meson decays), two scenarios of nuclear form factor (proton bremsstrahlung) and higher-order corrections (QCD scattering). Only the primary protonnucleon interaction is taken into account, secondaries from hadronic interactions in cascade decays could lead to an improvement in the sensitivity and will be the object of future work. The dark photon is assumed to decay to pairs of leptons or quarks, and only decay channels producing at least two charged particles coming from a common vertex are used. With the selection applied, backgrounds are neglected and 90% CL exclusion contours are extracted and compared with those from past experiments. The SHiP detector is expected to have a unique sensitivity for m γ D ranging between 0.8 and 3.3 +0.2 −0.5 GeV, and ε 2 ranging between 10 −11 and 10 −17 .