Measurement of beauty photoproduction near threshold using di-electron events with the H1 detector at HERA

The cross section for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$ep \rightarrow e\, b\bar{b} X$\end{document} in photoproduction is measured with the H1 detector at the ep-collider HERA. The decay channel \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$b\bar{b} \rightarrow ee X'$\end{document} is selected by identifying the semi-electronic decays of the b-quarks. The total production cross section is measured in the kinematic range given by the photon virtuality Q2≤1 GeV2, the inelasticity 0.05≤y≤0.65 and the pseudorapidity of the b-quarks \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$|\eta(b)|, |\eta(\bar{b})|\leq2$\end{document}. The differential production cross section is measured as a function of the average transverse momentum of the beauty quarks 〈PT(b)〉 down to the threshold. The results are compared to next-to-leading-order QCD predictions.

Abstract The cross section for ep → e bbX in photoproduction is measured with the H1 detector at the ep-collider HERA. The decay channel bb → eeX is selected by identifying the semi-electronic decays of the b-quarks. The total production cross section is measured in the kinematic range given by the photon virtuality Q 2 ≤ 1 GeV 2 , the inelasticity 0.05 ≤ y ≤ 0.65 and the pseudorapidity of the b-quarks |η(b)|, |η(b)| ≤ 2. The differential production cross section is measured as a function of the average transverse momentum of the beauty quarks P T (b)  results are compared to next-to-leading-order QCD predictions.

Introduction
In ep collisions at HERA beauty quarks are mainly produced as bb pairs via the fusion of a quasi-real photon emitted by the incoming electron (or positron) and a gluon of the proton as depicted in Fig. 1a. This process is referred to as direct or pointlike and can be calculated using perturbative  Generic leading order diagrams for bb production in ep collisions. The diagram (a) is referred to as direct or pointlike, the diagram (b) is referred to as resolved or hadronlike quantum chromodynamics (QCD) due to the large scale provided by the mass of the heavy b-quark and the correspondingly small coupling α s . Resolved processes where the photon fluctuates into a hadronic state before undergoing a hard collision, as indicated in Fig. 1b, are expected to be largely suppressed compared to the direct production process, because of the large b-quark mass. Due to the dominance of the direct process over the resolved process, the production of b-quarks in ep collisions at HERA is an excellent testing ground for QCD predictions.
Theory uncertainties in the prediction of the cross section, which are mainly related to the renormalisation and factorisation scales, are expected to be smaller for beauty production than for charm production. The study of beauty photoproduction near threshold is of particular theoretical interest as the only hard scale in this process is provided by the b-quark mass, and other scales like the photon virtuality (Q 2 ≈ 0 GeV 2 in photoproduction) or the transverse momentum of the b-quark can be neglected.
At HERA the beauty cross section in photoproduction ep → e bbX has been measured by the H1 [1][2][3][4][5] and ZEUS [6][7][8][9][10][11][12] collaborations and compared to calculations [13][14][15] at next-to-leading order (NLO) QCD, performed in the fixed flavour number scheme in which the beauty quark is treated as massive. In general the predictions using the factorisation and renormalisation scale 2 do not agree well with the data. In particular at low values of the transverse momentum of the beauty quarks P T (b) ≈ 0 GeV, i.e. in the phase space region where the only hard scale involved is the b-quark mass m b , the measurements show a tendency to lie above the prediction. The choice of a lower scale, 2 , leads to a better agreement of the prediction with the data [16].
In the present analysis a measurement of the differential beauty cross section at HERA in photoproduction as function of the quadratically averaged transverse momentum of the produced beauty quarks, dσ /d P T (b) , is made down to the bb-production threshold, using a novel technique based on low momentum electron identification.
Most of the previous beauty measurements at HERA in photoproduction and deep-inelastic scattering (DIS, Q 2 1 GeV 2 ) identified jets of b-quarks using single leptons tags [1,4,5,7,[9][10][11][17][18][19][20] or displaced vertices [2,6,[21][22][23][24]. Jet-based b-tagging algorithms are in general very efficient at high transverse momenta of b-quarks, P T (b) > 6 − 7 GeV, but degrade significantly for lower values due to the absence of the boost and the short decay length. In addition b-quarks almost at rest lead to isotropic decay topologies of the final state where jet finders usually fail. A second class of analyses used double tags to select bb pairs either by reconstructing two muons [8] or a muon and a D * meson [3,10], utilising the semileptonic decay channel b → μX and the decay channel b → D * X , respectively. Lower values of the transverse momentum of the b-quarks become accessible by the use of lepton tags without requiring jets, where the minimum P T (b) value is determined by the minimum transverse momentum cut on the lepton. For muons this cut is typically at P T (μ) ≈ 2 GeV, and therefore too high to measure efficiently the production cross sections of b-quarks near threshold.
In the present analysis the differential beauty cross section is measured using electron pairs, exploiting the doublesemileptonic decay bb → eeX , with online and offline P T (e) thresholds for the electron identification of about 1 GeV. The events were recorded by identifying low momentum electrons already online using a dedicated trigger, which recorded data in the year 2007 with a corresponding integrated luminosity of 48.1 pb −1 . This low cut on the transverse electron momentum, P T (e), improves not only the total acceptance but also makes the low P T (b) phase space experimentally accessible.

Monte Carlo simulations and QCD calculations
The Monte Carlo generators PYTHIA [25,26] and CAS-CADE [27][28][29] are used to determine the signal efficiency and the detector acceptance for the process ep → e bbX → e eeX , and to simulate the production of charm quarks. Differences in the predictions are taken into account as systematic uncertainty, see Sect. 5.4. For the production of J /ψ mesons only CASCADE and for the production of light quarks in photoproduction only PYTHIA is used. Deepinelastic scattering is simulated using the Monte Carlo generator RAPGAP [30].
In PYTHIA leading order matrix elements are implemented taking into account the mass of the heavy quarks. The CTEQ6L [31] set of proton parton density functions is used. The parton shower evolution in PYTHIA is based on the DGLAP equations [32][33][34]. In addition to the direct process, the resolved photon component is calculated by using the photon parton density function SAS 2D [35].
For the CASCADE simulation the direct γ * p → bb and γ * p → cc processes are implemented using off-shell matrix elements, which are convoluted with k T -unintegrated 1 proton parton density functions. The A0 [36] set of parton density functions is used. The parton evolution in CASCADE is based on the CCFM evolution equation [37][38][39][40] for the initial state parton shower.
In PYTHIA, CASCADE and RAPGAP higher order QCD corrections are included by simulating parton showers in the initial and final state. These Monte Carlo generators use the Lund String Model [41,42] for simulating the hadronisation of light quarks. For the hadronisation of heavy quarks the Bowler fragmentation model [43] is employed with parameters as used in a previous analysis [21].
In order to correct for detector effects and to estimate the systematic uncertainties associated with the measurement, the generated events are passed through a detailed simulation of the detector response based on the GEANT program [44] and through the same reconstruction and analysis software as is used for the data.
Theory cross sections are calculated in NLO QCD in the fixed flavour number scheme using the program FMNR [13][14][15] in order to compare with the data. These calculations are expected to give reliable results in the kinematic region considered here, where the transverse momentum of the heavy quark is of the same order of magnitude as its mass. The calculations are performed as a function of the quadratically averaged transverse momentum of the produced beauty pair The prediction of FMNR is evaluated for the direct and resolved photon processes. For the proton the CTEQ6M [31] set and for the photon the GRV-HO [45] set of parton density functions are used. In this analysis, the renormalisation and factorisation scales are chosen to be equal, μ R = μ F = μ 0 , with μ 0 = 1/2 m 2 b + P T (b) 2 and m b = 4.75 GeV. The value used for the QCD scale Λ QCD corresponds to the value of the strong coupling constant α s (M Z ) = 0.118. The theoretical uncertainty of the prediction is evaluated by varying the scales μ R and μ F simultaneously in the window μ 0 /2 < μ R,F < 2μ 0 and the beauty mass in the range 4.5 < m b < 5.0 GeV. By recalculating the cross section with different parton density functions the theoretical uncertainty due to the choice of the photon and proton parton density functions is found to be much smaller than the theoretical uncertainties and thus is neglected. 1 k T denotes the transverse momentum of the parton.

H1 detector
A detailed description of the H1 detector can be found elsewhere [46][47][48]. In the following, only detector components relevant to this analysis are briefly discussed. The origin of the H1 coordinate system is the nominal ep interaction point, with the direction of the proton beam defining the positive z axis (forward direction). Transverse momenta are measured in the x-y plane. Polar (ϑ ) and azimuthal (φ) angles are measured with respect to this reference system. The pseudorapidity is defined to be η = − ln tan(ϑ/2).
In the central region (15 • <ϑ <165 • ) the interaction point is surrounded by the central tracking detector (CTD). The CTD comprises two large cylindrical jet chambers (CJC1 and CJC2) and the silicon vertex detector [49]. The CJCs are separated by a drift chamber which improves the z coordinate reconstruction. The CTD detectors are arranged concentrically around the interaction region in a solenoidal magnetic field of 1.16 T. The trajectories of the charged particles are measured with a transverse momentum resolution of σ (p T )/p T ≈ 0.2 % p T /GeV ⊕ 1.5 %. In addition the CJCs provide a measurement of the specific ionisation energy loss dE/dx of charged particles with a relative resolution of 6.5 % for long tracks. A set of five cylindrical multiwire proportional chambers [50] mainly used for first level triggering are situated inside the inner CJC1 covering the polar angular region 11 • < ϑ < 169 • . The forward tracking detector and the backward proportional chamber measure tracks of charged particles at smaller (7 • <ϑ <25 • ) and larger (155 • <ϑ <175 • ) polar angles than the central tracker, respectively.
The liquid argon (LAr) sampling calorimeter [51] surrounds the tracking chambers and has a polar angle coverage of 4 • <ϑ <154 • . It consists of an inner electromagnetic section with lead absorbers and an outer hadronic section with steel absorbers. The LAr calorimeter is divided into eight wheels along the beam axis. The electromagnetic and the hadronic sections are highly segmented in the transverse and the longitudinal directions. Energies of electromagnetic showers are measured with a precision of σ (E)/E = 12 %/ √ E/ GeV ⊕ 1 % and energies of hadronic showers with σ (E)/E = 50 %/ √ E/ GeV ⊕ 2 %, as determined in test beam experiments [52,53]. In the backward region (153 • <ϑ <178 • ), particle energies are measured by a lead-scintillating fibre calorimeter (SpaCal) [48].
The luminosity is determined from the rate of the elastic QED Compton process ep → e γp, with the electron and the photon detected in the SpaCal calorimeter, and the rate of DIS events measured in the SpaCal calorimeter [54].
For data collection a four level trigger system is employed, of which the first two levels are implemented in hardware. The first level trigger (L1) is based on various sub-detector components, which are combined and refined at the second level (L2). The third level (L3) is a software based trigger using combined L1 and L2 trigger information from various subdetector components. Fully reconstructed events are subject to an additional selection at the software filter farm (L4).
The data used for this measurement were recorded by the Fast Track Trigger (FTT) [55][56][57][58][59][60] which, based on hit information provided by the CJCs, reconstructs tracks with subsequently refined granularity at the first two trigger levels, first in the x-y plane at L1 and then in three dimensions at L2. Of special importance is the third trigger level integrated in the FTT [60], which identifies low energy electrons (E > 1 GeV) [61,62] by combining FTT tracks with energy depositions reconstructed in the LAr calorimeter by the Jet Trigger (JT) [63].

Experimental method
The data sample used for this analysis was recorded in the year 2007, when positrons at an energy of 27.6 GeV collided with protons at 920 GeV, and when all trigger levels of the FTT and the JT were in operation. The recorded data corresponds to a total integrated luminosity of 48.1 pb −1 .
In this analysis the measurement of bb photoproduction is based on the identification and selection of two electrons in the LAr calorimeter at low transverse momentum P T (e) > 1 GeV to tag the semi-electronic decays of the bquarks. In about 2 % of all bb-decays two electrons originate from the same b-quark from the decay chain b → c e −ν e → s e − e +ν e ν e . In about 4 % of all bb-decays the two electrons originate from decays of different b-quarks, where they are either produced directly in the semi-electronic bdecays or in the subsequent semi-electronic c-decays. The electrons can be either of opposite charge (combinations bb → e + e − X and cc → e + e − X) or of same charge (combinations bc → e − e − X andbc → e + e + X). These charge relations hold only in the case of no BB mixing. In the following all possible combinations including BB mixing are considered in order to discriminate bb decays against semielectronic decays of cc events. Electron pairs from J /ψ decays are distinguished from those from b-decays by reconstructing their invariant mass. Misidentified electrons originating mainly from the light quark background are con-strained by varying the cuts on the electron identification described in Sect. 5.

Online electron identification
Events containing several tracks and one or two electron candidates compatible with the signature of semi-electronic b-decays are triggered, using the FTT on the trigger levels L1 to L3. On the first trigger level more than five tracks with transverse momentum thresholds in the range 0.1-1.8 GeV are required. These high multiplicities are verified at the second trigger level, exploiting the higher track resolution available at this level. On the third trigger level the track information as determined by FTT-L2 is combined with the energy depositions as measured in the LAr calorimeter by the Jet Trigger [60][61][62] to identify electrons. Electron candidates are required to fulfil a geometrical track-cluster matching condition using the distance variables ϑ = |ϑ FTT − ϑ JT | and ϕ = |ϕ FTT − ϕ JT |. In addition the transverse momentum P T ,FTT as measured with the FTT-L2 has to be compatible with the associated transverse energy E T ,JT measured in the LAr calorimeter by the JT. A lower cut on the quantity E T ,JT /P T ,FTT is used to discriminate electrons against hadrons, which deposit significantly less energy in the noncompensating LAr calorimeter.
For this analysis three subtriggers are used, which have identical L1 and similar L2 trigger conditions, but different conditions on L3 as summarised in Table 1. The subtrigger with the lowest transverse momentum threshold of P T ,FTT > 1.2 GeV requires events with at least two electron candidates. The other two subtriggers select events with a minimum of one electron candidate with P T ,FTT thresholds of 1.5 and 2.0 GeV. The three data sets recorded by these FTT-JT based subtriggers cover an overlapping kinematic phase space, but correspond to different integrated luminosities due to different trigger prescale factors. The three data sets are combined using a weighting method [61] to account for correlated triggers with prescales. The individual prescale corrected luminosities are also given in Table 1.

Offline electron identification
Electrons in the polar angle range of 20 • < ϑ(e) < 140 • and with a transverse momentum of P T (e) > 1 GeV, with Table 1 L3-online cuts used to trigger electron candidates. Explanations to the cuts are given in the text. The last column contains the prescale corrected integrated luminosity of each subtrigger. The medium-P T subtrigger was commissioned at a later stage 20 33.5 P T (e) and ϑ(e) measured from the electron track, are identified using energy depositions in the LAr calorimeter and specific ionisation loss measured in the CJCs. Two estimators, D calo and D dE/dx , are defined to discriminate electrons from background. The background, which is mainly due to pions misidentified as electrons and to a lesser extent due to kaons and anti-protons, is largely suppressed by combining the two independent estimators into a combined estimator D ele , as explained in Appendix A. The three estimators are defined such that D = 1 for genuine electrons and D = 0 for pion background.
The calorimeter based electron identification [61] is track seeded, which means the cluster shape estimators are calculated from energy deposits in LAr calorimeter cells lying within a cylinder of 30 cm around the extrapolated track trajectory. The cluster energies are corrected for energy losses in the dead material in front of the LAr calorimeter. Electron candidates with energy depositions close to inactive regions between LAr calorimeter modules are rejected. Five estimators are defined: four cluster shape variables and the ratio of the energy deposited in the electromagnetic part of the LAr calorimeter to the momentum of the corresponding track. These estimators together with the logarithm of the total energy and the z position of the cluster, are mapped onto one single estimator D calo using the artificial neural network Multilayer Perceptron [64].
The measured specific ionisation loss of the track, dE/dx, is translated into χ 2 -probabilities of corresponding particle hypotheses P (dE/dx, e) for electrons and P (dE/dx, π) for pions, which constitute the main background. From both probabilities the estimator is constructed. The simulation of the specific ionisation was studied in detail in order to describe precisely the measured energy losses [65]. The performance of both discriminator variables is validated using Monte Carlo and data samples of identified electrons and pions in the transverse momentum range of interest, 1 < P T (e) < 5 GeV, selected in decays J /ψ → e + e − and K 0 s → π + π − , by means of the "tag and probe method" [61].
The simulation describes well the distribution of the discriminators D calo and D dE/dx as measured in data, as can be seen in Fig. 2. The deviations of the simulation from the data at small D values in the electron sample are due to a small remaining pion contamination in the data. Also the combined estimator D ele is found to be well described by the simulation and shows an excellent separation of the electron signal from the pion background. Isolated electrons are selected for D ele > 0.825 with an efficiency of more than 90 % for a pion background rejection of about 99 %.

Event selection
A di-electron sample is obtained by selecting events with two or more offline reconstructed electron candidates, requiring D ele > 0.825. To account for the P T resolution of the third trigger level, the P T cut on electron tracks reconstructed offline is raised by 100 MeV above the trigger threshold of the respective subtrigger (see Table 1), which recorded the event. This implies two electrons with P T > 1.3 GeV for the low-P T subtrigger and one electron with P T > 1.6 GeV (P T > 2.1 GeV) for the medium-P T (high-P T ) subtrigger. If the event is triggered by the medium-P T or high-P T subtrigger, the transverse momentum of the second electron has to fulfil P T > 1.0 GeV.
In order to remove background from non-ep sources, the measurement of a good event vertex is required. The event vertex is reconstructed from all charged tracks of an event and its position along the beamline has to be within 30 cm around the nominal interaction point. In addition, timing vetoes are applied to further reduce non-ep interaction induced backgrounds.
In order to reject background from DIS, events with a positron in the LAr calorimeter identified by the standard electron identification [67] and with E(e + ) > 8 GeV are rejected. As the P T (e)-distribution of semileptonically decaying b-quarks falls steeply, almost all b-decay positrons are at low energies and thus not affected by this cut. DIS events with an electromagnetic cluster in the SpaCal calorimeter with an energy above 8 GeV consistent with originating from the scattered beam positron [66] are also rejected. Events with Q 2 2 GeV 2 are not rejected by these cuts, since the beam positrons leave the detector undetected along the beam pipe.
Only events with measured inelasticities in the phase space region of this measurement, 0.05 < y h < 0.65, are accepted. The inelasticity variable is reconstructed from the sum over all final state particles y h = i (E i − P z,i )/ (2E e + beam ), where E e + beam denotes the energy of the beam positron. Particles belonging to the hadronic final state (HFS) are reconstructed using a combination of tracks and calorimeter deposits in an energy flow algorithm that avoids double counting [68][69][70]. E i and P z,i denote the energies and longitudinal momenta of all final state particles, which correspond to the visible hadronic final state in case of photoproduction, and in case of DIS background also includes the scattered positron. The upper cut on the inelasticity suppresses effectively remaining DIS events.
The beauty signal is further enriched by rejecting electron candidates, which are in a dense hadronic environment. For this purpose the variable R E,cone is defined as the ratio of the summed energy of all HFS particles in a cone of 18 • around the electron track direction, E cone , to the electron energy E e , which must not exceed an upper threshold: The effect of this cut is twofold: First, it reduces misidentified electron candidates resulting from overlapping showers in the LAr calorimeter. Second, it enriches electrons from semileptonic beauty decays, which are in general isolated from hadrons due to the large b-mass. Finally, electrons from photon conversions are rejected by the three following cuts. First, the distance of closest approach in the transverse plane, dca e , of the electron tracks to the beam line is restricted to be smaller than 0.2 cm. Second, a photon conversion finder searching for displaced vertices is used to identify electrons originating from the photon conversion process γ → e + e − . Third, the invariant mass of the selected electron pairs is required to be m e1,e2 > 1.2 GeV. This cut rejects e + e − pairs from Dalitz decays and most of the remaining background from photon conversions.
The selection cuts are summarised in Table 2. After applying all cuts about 1500 electron pairs are selected. In the rare cases with more than two selected electrons per event all pair combinations are considered in the analysis.

Data analysis
The selected di-electron sample is dominated by events from inelastic J /ψ-meson production. While decays of J /ψmesons can be easily identified by kinematic reconstruction of the J /ψ mass peak, the separation of the bb signal events from the other backgrounds originating from the production of light quarks and open charm production is more difficult.
In the following, the reconstruction of the transverse momentum of the produced b-quarks and the flavour separation of the different processes are described.

Reconstruction of b-quarks
The transverse momentum of b-quarks is reconstructed for the measurement of the differential cross section dσ/d P T (b) , where P T (b) is the quadratically averaged transverse momentum of the b andb quark as defined in (1). The bb cross section is largest at small transverse momentum at P T (b) ≈ m b , a kinematic region where standard jet finders cannot be used due to isotropic decay topologies. Therefore an alternative, referenced as the hemisphere method, is exploited. This method was applied in a previous analysis [72] to reconstruct the directions and momenta of charm quarks in the production of cc-pairs in DIS, and is also well suited to reconstruct the transverse momenta of b-quarks in bb production [61]. Table 2 Overview of the online and offline selection cuts. More details on the selection procedure can be found in [61] Overview of the selection cuts  Table 1 and text)

Background rejection and further cuts
Rejection of non ep-background: • good vertex, timing vetoes Rejection of DIS events: • no identified scattered beam positron • 0.05 < y h < 0.65

Rejection of photon conversions and Dalitz decays:
• m e1,e2 > 1.2 GeV • no converted photon • dca e < 0.2 cm As illustrated in Fig. 3, an event is divided into hemispheres, using the thrust-axis which is calculated in the laboratory frame in the plane transverse to the beam directions (x-y plane). Using the transverse momenta from all particles of the HFS, the thrust-axis in the transverse plane is given by the vector a maximising the sum of the projected transverse momenta onto it, A plane perpendicular to the thrust-axis defines two hemispheres, one of them containing the fragmentation products of the b-quark, and the other one containing the fragmentation products of theb-quark. Two observables P T ,hem.I and P T ,hem.II are used to reconstruct the mean transverse momentum of the b (b) quark produced in the hard interaction. These observables, which are derived from the HFS particles assigned to the corresponding hemispheres, show a good correlation to the transverse momentum of the b (b) quarks in the hard process. However, the hadronic final state also contains particles from the so called proton remnant, leaving the interaction in the positive z-direction of the detector and thus deteriorating the above correlation. Simulation studies show that the correlation with the b-quark transverse momentum is improved by excluding particles in the forward direction at polar angles below 15 degrees. The transverse momenta of the b (b)- Fig. 3 Schematic illustration of the determination of the thrust axis in the plane transverse to the ep beams. The transverse thrust axis, indicated by the dashed arrow, maximises the sum of momenta projected onto it in this plane. The thrust axis allows the event to be divided into two hemispheres, each containing the decay products of a beauty quark, used to reconstruct the average transverse beauty mass m T ,rec (b) as defined in (7) quarks are therefore approximated by: This reconstruction method is very reliable at large P T ,b , where two hard jets are measured in the final state. At small P T (b) the transverse momenta of HFS particles in the hemispheres are mainly generated by the b andb-hadron decays themselves and are related to the mass of the b-quark: | P T ,hem.I | ≈ | P T ,hem.II | ≈ m b . In order to allow for a good reconstruction of P T (b) down to the bb production threshold, i.e. P T (b) ≈ 0 GeV, the average transverse beauty mass is used: Detailed studies [61] demonstrated that the average transverse beauty mass can be well reconstructed from the experimental observables P T ,hem.I and P T ,hem.II using the relation: with α being a constant parameter set to α = 1.09, such that the correlation between generated and reconstructed m T (b) is maximised. This correlation as obtained by simulation is shown in Fig. 4. For values of m b it in the range 4.5 < m b < 5.0 GeV the dependence of this correlation on m b is negligible.

Quark flavour separation
For the discrimination of the bb signal against remaining background from misidentified electrons and for the separation of the different quark-flavour components contributing to the di-electron signature, a template method is used. Several independent phase space regions are defined such that individual background sources are enhanced in certain regions of the phase space and can be tested while other contributions are suppressed. Finally the b-signal ("beauty") and the background contributions are obtained by an unfolding procedure. Background sources determined by this method are the production of light quarks ("uds"), open charm production ("charm") and the production of J /ψmesons ("J /ψ"). The uds background contains also a small fraction of charm and beauty events, where at least one electron candidate does not originate from a semi-electronic heavy quark decay.

Fraction of light quarks
In order to determine the background contributions due to misidentified electrons the data sample is grouped in four regions B1, B2, B3 and S using different electron quality criteria on D min(e1,e2) ele and R max(e1,e2) E,cone , see Table 3.   these three background bins mainly constrain the uds background fraction.

Heavy quark fractions
In the signal enhanced region S, the individual contributions from beauty, charm and J /ψ can be disentangled by investigating the charge product, q e1 · q e2 , of the e ± -candidates, their azimuthal separation φ e1,e2 = |φ e1 − φ e2 |, and their invariant mass m e1,e2 . Templates of the different background sources and of the beauty signal, which are all restricted to the signal enhanced region S, are shown in Fig. 6 as function of the invariant mass m e1,e2 and the signed azimuthal separation φ e1,e2 · q e1 · q e2 . The different templates show specific characteristics: J /ψ events have oppositely charged electrons and cluster at m e1,e2 = m J/ψ , whereas background from open charm production covers a large mass range. Electrons from open charm decays are found mostly back-to-back and with opposite charge sign, whereas electron pairs from beauty decays populate all φ e1,e2 values with both charge sign combinations. Both charge products are also found in the uds background, which however populate on average regions with smaller m e1,e2 values. Large values of m e1,e2 are solely populated by beauty decays.  Table 3 and text These distinct signatures of the individual background sources, i.e. uds, J /ψ → ee and cc → ee, are exploited by dividing the signal enhanced region S into 12 subregions (S1 to S12) as shown in Fig. 6. In the following the three background enhanced bins B1-B3 and the 12 signal enhanced bins S1-S12 are referred to as "Flavour Separator", for which templates are derived.

Unfolding
Using an unfolding procedure the number of background events N uds , N J/ψ , N charm and the number of beauty events N beauty,i in four bins of P T (b) are derived. A regularised unfolding procedure is used with a smoothness condition. The procedure is explained in Appendix B. All efficiency corrections and migration effects are described by the response matrix A, which correlates the number of reconstructed events in the Flavour Separator distribution in bins of m T ,rec , represented by the vector y, with the distribution x on parton level via the matrix equation The vector x, defined as x T = (x T beauty , x charm , x J/ψ , x uds ), contains contributions from beauty binned in P T (b) , charm, J /ψ and uds. The contribution from beauty (x beauty ) is defined according to the phase space given in Table 4. The vector b contains the background contribution from DIS events, which is taken from simulation. All other background contributions are incorporated in the response matrix and are determined by unfolding.
Signal and background templates as function of nine m T ,rec bins are generated by Monte Carlo simulations and fitted to the data. The unfolding procedure uses in total N m T ,rec × N Flavour Separator = 9 × 15 input bins and determines the three background fractions and the number of beauty events N beauty,i in four P T (b) bins. A schematic representation of the procedure is shown in Fig. 7. In this procedure the m T ,rec dependences of the different background contributions from uds, J /ψ → ee and cc → ee are fixed by the Monte Carlo predictions. The latter is motivated by recent measurements of the differential cross sections of charm production at HERA, which were found to be consistent with theoretical models and Monte Carlo programs used in this analysis [73,74].
The fitted beauty signal and background contributions are shown in Fig. 8 in the three background and in the signal enhanced regions. The event numbers resulting from the fitted fractions show very good agreement with the data considering statistical errors only. A clear enhancement of the genuine electron signal due to the tightening of the electron identification cuts is seen when going from the first Fig. 6 Templates in the plane φ e1,e2 · q e1 · q e2 , m e1,e2 and restricted to the signal enhanced region S used to separate the heavy quark flavours as obtained by the Monte Carlo simulation. Also shown is the bin numbering S1-S12 of the 12 subregions of S. For the definition of the signal enhanced region S see Table 3 and text. The two vertical lines indicate the peak invariant mass region of the J/ψ → e + e − decays background enhanced region (B1), which contains more than 80 % uds background to the signal enhanced region (S1-S12) with less than 20 % of uds background.
The correlations between the beauty signal and the background sources, which are largest between beauty and charm, are given in Table 5 together with the determined fractions of the selected data sample.
The distribution of the data as a function of the Flavour Separator is shown in Fig. 9 together with the result from the fit of the beauty and the various background contributions. Good agreement is found considering statistical errors only.
Control distributions of electron variables are presented in Fig. 10 for the electron enriched signal region (S1-S12). The data are compared to the simulated beauty signal and Fig. 7 Structure of the response matrix A used to solve the matrix equation y = A · x + b by unfolding. The indicated numbers specify the number of used bins. The sub-binning in y given by the Flavour Separator allows the discrimination of the beauty signal from the uds, charm and J/ψ backgrounds. See text for details Table 5 Correlations between the signal (beauty) and the different background contributions, and the determined relative fractions with their errors for the data sample. The fraction of DIS events (not given in the background distributions using the quark flavour decomposition determined by the unfolding procedure. The main characteristics of the signed variables φ e1,e2 · q e1 · q e2 and m e1,e2 · q e1 · q e2 , and P T (e) and ϑ(e) are well described by the Monte Carlo simulation. In Fig. 11 additional control distributions are presented showing the P T -spectra of the three highest P T tracks. These distributions are strongly dependent on the track trigger conditions used, and imperfections of the trigger simulation would be visible here.
Reasonable agreement between the data and the Monte Carlo simulation is obtained in all distributions which gives confidence that the Monte Carlo simulation is able to correctly model the detector response used for the unfolding procedure.

Cross section determination and systematic uncertainties
The visible cross section is measured for the phase space as defined in Table 4. The bin-averaged differential cross section is obtained as where L is the luminosity, P T ,i (b) the bin width, N beauty,i the number of unfolded signal events in the corresponding bin and BR = 6.17 % the effective branching fraction computed from [71] for a bb pair decaying into at least two electrons. For the calculation of cross section uncertainties correlations between bins are taken into account.
The systematic uncertainties related to the measurement of the number of bb signal events are listed in the following. The effect on N beauty,i is calculated by varying the sources of uncertainties in the simulation and by propagating these variations to the measurement through the response matrix A and the background term b in (8).
• The uncertainty on the electron identification is determined using J /ψ → e + e − events (see Fig. 2 Table 3. Data are represented as points with the statistical error indicated by the error bars. Also shown in colour is the decompositions of the event yields as determined by the unfolding procedure trigger uncertainty is due to the uncertainty of the calibration constants of the JT used at L3. To quantify this uncertainty, the JT calibration constants used in the simulation are varied by scaling the default calibration constants by 15 % [61]. The systematic error on the total beauty cross section due to the uncertainty on the trigger efficiency is determined to be ±8.6 %. • The uncertainty of the charm contribution is evaluated from the relative difference between the Monte Carlo generators CASCADE and PYTHIA in a similar way as for the beauty signal. The systematic error on the extracted total beauty cross section due to the charm model is determined to be ±3.6 %. • The uncertainty due to the fragmentation function of the heavy quarks is estimated by reweighting the events according to the longitudinal string momentum fraction z carried by the heavy hadron in the Lund model using weights of (1 ∓ 0.7) · (1 − z) + z · (1 ± 0.7) for charm quarks and by (1 ∓ 0.5) · (1 − z) + z · (1 ± 0.5) for beauty quarks [22]. The corresponding systematic error on the  Fig. 6 and Table 3 is used total beauty cross section is determined to be ±3.4 % resulting from the charm and ±2.2 % from the beauty fragmentation uncertainty. • The uncertainty on the contribution from the remaining uds background due to misidentified and real electrons, was determined by varying their relative contributions by a factor two up and down. The corresponding systematic error on the total beauty cross section is determined to be ±3.4 %. • CASCADE does not fully simulate the radiative tail of J /ψ → ee events. To estimate the uncertainty on the modelling of it, weights are applied, which are obtained from an elastic J /ψ → ee simulation with radiative QED corrections [75]. The systematic uncertainty on the total beauty cross section is estimated to be ±3.5 %. • The uncertainty of the DIS-background, represented in (8) by the vector b, is taken to be 100 % and results in an error on the total beauty cross section of ±4.5 %.
In addition, a global normalisation uncertainty of 4.1 % is included with contributions from the integrated luminosity Fig. 10 Control distributions of the electron candidates compared to Monte Carlo simulations using the quark flavour decomposition determined by the unfolding procedure: (a) signed azimuthal separation φ e1,e2 · q e1 · q e2 defined by the charges multiplied with the azimuthal angle difference of the two electron candidates, (b) signed invariant mass m e1,e2 · q e1 · q e2 defined by the charges multiplied with the invariant mass of the two electron candidates, (c) polar angle of the electron candidates and (d) transverse momentum of the electron candidates. Data are represented as points with the statistical uncertainties indicated by the error bars. The distributions are restricted to the electron enriched region (S) Fig. 11 Control distributions for the three highest P T tracks of the hadronic final state as function of the track P T . Data are compared to the Monte Carlo simulations using the quark flavour decomposition determined by the unfolding procedure. Data are represented as points with the statistical uncertainties indicated by the error bars uncertainty of 2.7 % and from the uncertainty on the semielectronic branching fractions of 3.0 %.
Adding all above contributions in quadrature gives a total systematic error of 15.4 % on the total beauty cross section.

Results
The differential cross section dσ (ep → e bbX)/d P T (b) is measured in the phase space defined in Table 4 using the unfolding procedure as described in Sect. 5.3.
The result is shown in Table 6 together with statistical and total errors and the coefficients describing the statistical correlations between bins. In order to cross check the unfolding procedure the cross section extraction is repeated without regularisation condition. The results obtained with and without regularisation are found to be consistent within the uncertainties.
The measured differential beauty cross section is compared in Fig. 12 with an NLO QCD prediction in the fixed flavour number scheme as calculated by the program FMNR. The figure also shows the ratio of the measured cross section and the NLO QCD cross section. The uncertainties of the measurement are smallest at low P T (b) , where the cross section is largest. The theoretical prediction of the differential cross section agrees with the measurement within the large experimental and theoretical uncertainties. The prediction has a tendency to be below the data, a trend also observed in previous beauty cross section measurements at large transverse beauty momenta. By integrating the differential cross section the total inclusive beauty photoproduction cross section is measured as: σ (ep → e bbX) = 3.79 ± 0.53 (stat.) ± 0.58 (sys.) nb, (10) to be compared with the NLO prediction obtained from FMNR of σ (ep → e bbX) = 2.40 +0. 55 −0.49 nb. The measured cross section is higher, but within the large experimental and theoretical uncertainty consistent with the NLO expectation. Table 6 Differential cross sections for the phase space defined in Table 4 obtained from unfolding with regularisation condition. Also given are the statistical and total errors, the coefficients of the statistical correlations, the used bin boundaries in P T (b) and the corresponding bin centres [76] P T ,bc (b) . The remaining columns list the bin-to-bin correlated systematic uncertainties in the cross section measurement due to uncertainties of the beauty (δ b sys ) and charm (δ c sys ) modelling, the uds background (δ uds sys ), the electron identification (δ e−Id sys ), the beauty (δ fr.b sys ) and charm (δ fr.c sys ) fragmentation, the modelling of the radiative tail of J/ψ

Conclusions
The inclusive and differential cross section of beauty photoproduction was measured in the di-electron final state, using the H1 detector at the HERA collider. The cross section is measured as function of the quadratically averaged transverse momentum of the produced beauty quarks P T (b) , with a special focus on the low P T (b) regime. Background from uds, charm and J /ψ production is determined exploiting angular, charge and mass correlations of electron pairs in an unfolding procedure. The measured cross section is compared to a QCD prediction at NLO performed in the fixed flavour number scheme and evaluated with μ R = μ F = 1/2 m 2 b + P T (b) 2 as choice for the renormalisation and factorisation scale. The NLO prediction lies below the data but within the large experimental and theoretical uncertainty they agree.
This measurement is in good agreement with previous measurements of beauty photoproduction at HERA and it extends the previously experimentally accessible phase space towards the beauty production threshold.
This equation describes the minimisation of the unfolding problem χ 2 A (x) with the two side conditions given by χ 2 L (x) and χ 2 N (x). The actual minimisation problem is defined by the standard χ 2 function: with V = cov(y i , y j ) being the covariance matrix of the data. This function minimises the deviation of the estimator Ax from the measured, and background subtracted vector y − b.
The additional constraints are given by: with L being the regularisation matrix, m the number of reconstructed bins and n obs the total number of observed events after background subtraction, which ensures that the total number of events is conserved. Both functions enter equation 13 with the parameters τ and μ, where τ is often denoted as regularisation parameter and μ as Lagrange Multiplier. The χ 2 L (x) function is a measure for the smoothness of the result. The matrix L is chosen such that the second derivative ofx between bins describing beauty production is minimised. The regularisation parameter τ determines the strength of the smoothness constraint. For the regularised unfolding τ is chosen such that the correlations of the covariance matrix of the unfolded distributionx are minimised [78].