Measurement of Beauty and Charm Photoproduction using Semi-muonic Decays in Dijet Events at HERA

Measurements of cross sections for beauty and charm events with dijets and a muon in the photoproduction regime at HERA are presented. The data were collected with the H1 detector and correspond to an integrated luminosity of 179 pb^{-1}. Events with dijets of transverse momentum P_T^{jet1}>7 GeV and P_T^{jet2}>6 GeV in the pseudorapidity range -1.5<\eta^{jet}<2.5 in the laboratory frame are selected in the kinematic region of photon virtuality Q^2<2.5 GeV^2 and inelasticity 0.22.5 GeV in the pseudorapidity range -1.3<\eta^{\mu}<1.5. The fractions of beauty and charm events are determined using the impact parameters of the muon tracks with respect to the primary vertex and their transverse momentum relative to the axis of the associated jet. Both variables are reconstructed using the H1 vertex detector. The measurements are in agreement with QCD predictions at leading and next-to-leading order.


Introduction
The photoproduction of heavy quarks at the HERA electron-proton collider is of particular interest for testing calculations in the framework of perturbative quantum chromodynamics (pQCD). The masses m b and m c of the beauty and charm quarks, as well as the transverse momentum of a jet, provide a hard scale, which is necessary for the calculation of pQCD predictions [1]. The dominant process for beauty and charm production at HERA is boson-gluon fusion (BGF): γg → QQX, with Q = b, c. In the kinematic regime of photoproduction, where the photon virtuality Q 2 is small, two classes of processes contribute to BGF. In direct photon processes, the photon emitted from the electron 1 enters the hard BGF process directly. In resolved photon processes, the photon fluctuates into a hadronic state before the hard interaction and acts as a source of partons, one of which takes part in the hard interaction. At leading order (LO) pQCD resolved photon processes are expected to contribute significantly in the photoproduction region.
Heavy quark photoproduction has been investigated at HERA using different analysis techniques in different regions of phase space. The lifetime or mass of heavy flavoured hadrons [2][3][4][5][6][7], semileptonic decays [5,6,8] or the full reconstruction of a D or D * meson [9,10] are exploited to perform the measurements. In general, the measured cross sections agree reasonably well with the theoretical predictions. In the H1 measurement performed with HERA I data [5], where beauty photoproduction was investigated using two jets and a muon in the final state, the data were found to be described reasonably well by the pQCD calculations at next-to-leading order (NLO), except in the region of low transverse momentum of the muon 2.5 < P µ T < 3.3 GeV and of low leading jet transverse momentum 7.0 < P jet1 T < 10.0 GeV, where the NLO prediction was lower than the data. Comparable measurements were performed by the ZEUS collaboration [6,7], covering a similar phase space. Here, a good description by the NLO QCD predictions is found, also at low transverse momenta of the leading jet and the muon. This paper reports on beauty and charm measurements of cross sections for photoproduction of events with two jets and a muon, where the muon is associated with one of the jets. The analysis exploits the lifetime and the mass of heavy flavoured hadrons as in the former HERA I beauty production analysis [5]. The measurement is additionally performed for charm photoproduction. Increased statistics due to increased luminosity and extended phase space, a better understanding and description of the H1 vertex detector [11], which is crucial for this analysis, as well as reduced systematic uncertainties compared to the previous H1 analysis make a more detailed test of pQCD predictions possible. Going beyond the earlier HERA I analysis, the cross sections as a function of the azimuthal angular difference ∆φ jj between the two leading jets, which are sensitive to higher order corrections, are also measured. Furthermore, cross sections are measured in two different regions of phase space, which are either enriched by resolved or direct photon processes.

QCD Calculations
The data presented here are compared with LO calculations supplemented by parton showers as well as with NLO calculations. The calculations are performed using either collinear 1 In this paper the term 'electron' is used generically to refer to both electrons and positrons. 4 factorisation, which is based on the DGLAP evolution equations [12], or the k t -factorisation approach, which employs the CCFM [13] evolution equations. In the collinear approach transverse momenta obtained through initial state QCD evolution are neglected and all the transverse momenta are generated in the hard scattering process, i.e. the partons entering the hard interaction are collinear with the proton. Effects from non-zero transverse momenta of the gluons enter only at NLO. In the k t -factorisation ansatz the transverse momenta of gluons k t entering the hard interaction are already included at leading order both in the off-shell matrix element and the k t -dependent unintegrated gluon density [14]. Therefore, corrections appearing only at higher order in collinear factorisation are partially included at LO in the k t -factorisation approach.
For beauty and charm photoproduction two classes of processes occur, the direct photon processes and the resolved photon processes. The distinction between these two classes depends on the factorisation scheme and the order in which the calculation is performed.
The production of heavy quarks is calculated either in the massive scheme, where heavy quarks are produced only perturbatively via boson gluon fusion, or in the massless scheme, where heavy quarks are treated as massless partons. These two schemes are expected to be appropriate in different regions of phase space [15]: the massive scheme is expected to be reliable when the transverse momentum P T of the heavy quarks is of similar size compared to the heavy quark mass m Q , whereas the massless scheme is expected to be valid for P T ≫ m Q .
An overview of the parameters used in the Monte Carlo (MC) generators is given in table 1. The following MC generators are used: PYTHIA: The MC program PYTHIA 6.4 [16] is based on LO QCD matrix elements with leading-log parton showers in the collinear factorisation approach. PYTHIA includes both direct photon gluon fusion and resolved photon processes. In the resolved photon processes either a beauty or a charm quark or a gluon from the photon enters the hard scattering. In the inclusive mode of PYTHIA used here beauty and charm quarks are treated as massless partons in all steps of the calculation in both types of processes. Three samples are generated containing photoproduction events for the processes ep → bµX, ep → cµX and ep → qjX where q is a light quark of flavour u, d or s and j denotes a jet. The latter sample is generated without specifically requiring a muon in order to use it for studying the background arising from muon candidates which originate from sources other than beauty or charm processes. The hadronisation process is simulated according to the Lund string model [17]. For the systematic uncertainty arising from the fragmentation model, additional samples are generated using the Peterson fragmentation function [18] for heavy quarks.
CASCADE: The MC program CASCADE 2.0 [19] is used for simulating events based on LO QCD calculations in the k t -factorisation approach. Only the direct boson gluon fusion process is implemented using off-shell matrix elements. Higher order QCD corrections are simulated with initial state parton showers applying the CCFM evolution [13]. Here, two samples containing the processes ep → bµX and ep → cµX are generated. The unintegrated PDFs of the proton from set A0 [20] Table 1: Parameters used in the QCD calculations PYTHIA, CASCADE, HERWIG, and MC@NLO. Here, P T stands for the transverse momentum, m 1 , m 2 are the masses of the two outgoing partons from the hard process, and m Q denotes the heavy quark mass. The variableŝ represents the invariant mass of the QQ system and Q ⊥ stands for its transverse momentum.
HERWIG: The MC program HERWIG 6.510 [21] is used to simulate events based on collinear factorisation and massless LO QCD calculations. HERWIG includes both direct photon gluon fusion and resolved photon processes. The hadronisation of partons is performed using the cluster fragmentation [22].
MC@NLO: The MC@NLO program for HERA [1] combines a collinear NLO calculation of heavy flavour production in the massive approach [23] with parton showers and hadronisation as described in [24]. The direct and resolved parts of the cross section are calculated separately. MC@NLO uses parton showers applying the DGLAP evolution to simulate higher order contributions and the cluster fragmentation [22] as implemented in HERWIG [21]. The quark masses are set to m c = 1.5 GeV and m b = 4.75 GeV. The central value of the renormalisation scale µ R is set to m 2 Q + P 2 T , while the factorisation scale µ F is 2µ R . As an estimate of the theoretical uncertainties on the NLO QCD predictions the scales µ R and µ F are varied separately by factors of 0.5 and 2, while m b is changed by ±0.25 GeV and m c is changed by ±0.2 GeV. The resulting variations are added in quadrature to obtain the resulting theoretical uncertainty.
PYTHIA and CASCADE are used to simulate detector effects in order to determine the acceptance and the efficiency 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 simulation program [25] and are processed using the same reconstruction and analysis chain as is used for the data.

H1 Detector
Only a short description of the H1 detector is given here including the most relevant detector components for this analysis. A more complete description may be found elsewhere [30,31]. A right-handed coordinate system is employed at H1, with its origin at the nominal interaction vertex, its z-axis pointing in the proton beam direction and its x(y) axis pointing in the horizontal (vertical) direction. Polar (θ) and azimuthal (φ) angles are measured with respect to this reference system. The pseudorapidity η is related to the polar angle θ by η = − ln tan(θ/2).
Charged particles are measured in the central tracking detector (CTD) with a transverse momentum resolution of σ(P T )/P T ≈ 0.5%P T /GeV ⊕ 1.5% [32]. This device consists of two cylindrical drift chambers (CJC) interspersed with a drift chamber designed to improve the z-coordinate reconstruction. A multiwire proportional chamber mainly used for triggering is located in front of the inner CJC. The CTD is operated in a uniform solenoidal 1.16 T magnetic field, enabling the momentum measurement of charged particles over the polar angular range 20 • < θ < 160 • . The efficiency for finding tracks in the CTD is greater than 99%.
The CTD tracks are linked to hits in the vertex detector, the central silicon tracker (CST) [33], to provide precise spatial track reconstruction. The CST consists of two layers of doublesided silicon strip detectors surrounding the beam pipe, covering an angular range of 30 • < θ < 150 • for tracks passing through both layers. The information of the z-coordinate of the CST hits is not used in the analysis presented in this paper. For CTD tracks with CST hits in both layers the transverse distance of closest approach (DCA) to the nominal vertex in x − y, averaged over the azimuthal angle, is measured to have a resolution of 43 µm ⊕ 51 µm/(P T /GeV), where the first term represents the intrinsic resolution (including alignment uncertainty) and the second term is the contribution from multiple scattering in the beam pipe and the CST. The efficiency for linking hits in both layers of the CST to a CTD track is around 84%.
The track detectors are surrounded in the forward and central directions (4 • < θ < 154 • ) by a finely grained liquid argon calorimeter (LAr) and in the backward region (153 • < θ < 178 • ) by a lead-scintillating fibre calorimeter (SpaCal) both with electromagnetic and hadronic sections. These calorimeters provide energy and angular reconstruction for final state particles from the hadronic system. In the LAr electromagnetic shower energies are measured with a precision of σ(E)/E = 11%/ E/GeV ⊕ 1% and hadronic energies with σ(E)/E = 50%/ E/GeV ⊕ 1%, as determined in test beam measurements. The energy resolution for electromagnetic showers in the SpaCal is σ(E)/E = 7.1%/ E/GeV ⊕ 1%, as determined in test beam measurements [34].

Experimental Method
The data were collected with the H1 detector at the HERA collider during the years 2006 and 2007 and correspond to an integrated luminosity of L = 179 pb −1 . The beam energies were E e = 27.6 GeV and E p = 920 GeV for electrons and protons, respectively, resulting in a centreof-mass energy of √ s ≈ 320 GeV. The trigger requires a track segment in the muon system and track activity in the central jet chamber. A detailed account of the present analysis can be found in [35]. A summary of the kinematic range and the definition of the measurement is given in table 2.

Photoproduction Event Selection
Events in the photoproduction regime are selected by requiring that no isolated high energy electromagnetic cluster, consistent with a signal from a scattered electron, is detected in the LAr and SpaCal calorimeters. This limits the photon virtuality to values of Q 2 < 2.5 GeV 2 . The inelasticity y is reconstructed using the relation y = h (E − P z )/2E e [36]. Here, the sum includes all particles of the hadronic final state (HFS), while E denotes their energies and P z stands for the z-components of their momenta. The HFS particles are reconstructed using a combination of tracks and calorimeter deposits in an energy flow algorithm that avoids double counting [37]. The inelasticity in this analysis is restricted to 0.2 < y < 0.8.

Muon Reconstruction and Selection
Muon candidates are identified as track segments in the barrel and endcap parts of the instrumented iron. The iron track segments must be well matched to a track reconstructed in the CTD. At least two CST hits in the r − φ plane have to be associated with the muon track. The combined CTD-CST track in r − φ is required to have a fit probability of at least 10%. The muon momentum is reconstructed using the CTD-CST track information. The CST hit requirements for the muon track restrict the allowed range of ep interaction vertices along the z-axis to |z vtx | 20 cm. Events are selected with at least one muon candidate reconstructed in the instrumented iron having a pseudorapidity within −1.3 < η µ < 1.5 and a transverse momentum of P µ T > 2.5 GeV. If more than one muon candidate is found, the one with the highest transverse momentum is selected and other candidates are ignored. In 5.4% of the events after the full selection more than one muon is found.

Jet Reconstruction and Selection
Jets are reconstructed using the inclusive longitudinally invariant k T algorithm in the massless P T recombination scheme and with the distance parameter R 0 = 1 in the η − φ plane [38]. The algorithm is applied in the laboratory frame using all reconstructed HFS particles including the muon candidate. A jet is defined as a µ-jet if the selected muon candidate lies within a cone of radius 1 about the jet axis in the η − φ plane. The efficiency for this matching amounts to about 90% and is consistent for data and MCs. The jet with the highest P T is referred to as jet1, while the second highest is called jet2. Events with at least two jets are selected, where the leading two jets are required to be in the angular range −1.5 < η jet < 2.5 and to have a transverse momentum of P jet1 T > 7 GeV and P jet2 T > 6 GeV. One of the two selected jets must be classified as a µ-jet. The Monte Carlo simulation is used to define hadron level jets, which consist of stable particles including neutrinos, but excluding the scattered electron, before they are passed through the simulation of the detector response.

Separation of Direct and Resolved Processes
The fraction of the photon energy entering the hard interaction is estimated using the observable x obs γ : where the sums in the numerator run over the particles associated with the two jets and the one in the denominator over all detected hadronic final state particles. For direct processes x obs γ approaches unity, as the hadronic final state consists only of the two hard jets selected in the present analysis and the proton remnant in the forward region only has a minor contribution to HF S (E − P z ). In resolved processes x obs γ can have smaller values.

Flavour Separation
The flavour of an event is defined as the hadron flavour of the µ-jet. The measured cross sections are proportional to the rate of events with a muon and a dijet system rather than the rate of muons or jets.
The separation of b, c and light quark (uds) events is only briefly described here. The procedure closely follows that described in [5]. The separation is performed using the properties of the muon track associated to the µ-jet. The impact parameter δ of a muon track is the transverse distance of closest approach (DCA) of the muon track to the beam spot point, which is the position of the beam interaction region in x and y. The beam spot is derived from tracks with CST hits averaging over many events and is updated regularly to account for drifts during beam storage. Muon tracks with δ > 0.1 cm are rejected to suppress contributions from the decays of long-lived strange particles. If the angle α between the azimuthal angle of the µ-jet and the line joining the primary vertex to the point of muon closest approach is less than 90 • , the impact parameter is defined as positive. It is defined as negative otherwise. The transverse momentum P rel T of the muon relative to the µ-jet axis is also sensitive to the quark content of the event sample and used together with the impact parameter for the flavour separation.
The fractions of events with beauty, charm and light quarks, f b , f c and f l , are obtained by a binned likelihood fit [39] in the δ − P rel T plane. Following [39], a likelihood ratio is calculated based on Poisson statistics. The fit is performed separately for each individual measurement bin, while the total cross sections are determined using the fractions obtained from a fit to the complete event sample. The uds (light), c and b PYTHIA Monte Carlo simulation samples are used as templates. Only the statistical errors of the data and the Monte Carlo simulations are considered in the fit. As a cross check, all fits are also performed using one-dimensional distributions of P rel T and δ separately. These two one-dimensional fits give a compatible b fraction in all measurement bins within the statistical error. The one-dimensional P rel T fit does not allow a determination of the charm fraction.
The two distributions that are used in the flavour separation are shown in figure 1. The distribution of the impact parameter δ shown in figure 1 (a) is symmetric for uds events while b and c events contribute more at large positive values of δ. Therefore, the fit of this variable can distinguish the three different quark contributions. As shown in figure 1 (a), the sum of all three fitted contributions in the Monte Carlo simulation is able to describe the data quite well. This description is achieved by a better understanding of the detector and an improved detector simulation with regard to signal heights, noise levels and dead strips in the CST [11], the inclusion of effects from alignment imperfections, and the description of the dead material in front of the CST and CJC. Therefore, a further smearing of measured track parameters in the simulation is not necessary, as it was done in the former H1 analysis using HERA I data [5]. The P rel T distribution is shown in figure 1 (b). The uds and c distributions are very similar and peak at low values of P rel T , while the b events contribute more at higher values of P rel T . The sum of all quark contributions in the Monte Carlo simulation is able to describe the data reasonably well.
The fitted parameters f b , f c and f l for the whole kinematic range are: The χ 2 /ndf is found to be 0.76 for the total sample.
Control distributions for the data sample in comparison to the Monte Carlo simulations are shown in figure 2. All selection cuts are applied. The data are compared with the MC contributions from beauty, charm and light quark events and their sum with the relative fractions taken from the fit as discussed above. The number of events in the simulation is normalised to the one of the data. It is observed that the shapes of the MC contributions are rather similar for beauty, charm and light quark events for the distributions shown here. For the determination of the detector corrections a reweighting of ∆φ jj and the transverse momentum of the leading selected track in the event is performed on hadron level in the Monte Carlo simulation to provide a better description of the data. Only small deviations between data and MC are observed, such as in the forward η µ region in figure 2 (b) and in the very small ∆φ jj region in 2 (d).

Cross Section Determination
Total and differential visible beauty and charm ep cross sections are measured in the photoproduction regime. The fitted fractions f b and f c are converted to cross sections σ b(c) in each bin using stands for the events selected in the Monte Carlo simulation on hadron level and L denotes the luminosity of the data. The differential cross sections are obtained by dividing by the bin width.

Fake Muon Rate
All backgrounds to semi-muonic b and c decays are called fake muons here. These contributions are modelled using MC simulation and originate mainly from uds events, with a small fraction from b and c events. Three sources of fake muons are considered: • Hadrons which reach the muon detector and are misidentified as muons. According to the fully inclusive PYTHIA MC, 1.6% of the selected b events and 2.3% of the selected c events originate from hadron misidentification.
• Muons which do not originate from b or c hadron decays, but from other hadrons such as kaons and pions. This background source is denoted as inflight decay in the following. According to the full inclusive PYTHIA MC, 0.9% of the b and 0.7% of the selected c events originate from inflight decays.
• Cosmic ray muons which coincide in time with real ep events. About 1% of the selected muon candidates are rejected as cosmic ray muons based on timing information from the CTD [35]. The remaining background from cosmic ray muons is negligible.

11
The probability to fake a muon depends on the particle species. The kinematic distributions of different particle species differ, but the uds MC is used as one single template. Therefore the fake probabilities for the most important particle species are studied in data and MC. The fake muon probabilities from misidentification and inflight decays P µ h for hadron h are defined as: The fake muon probability is investigated for pions originating from K 0 S → π + π − decays and amounts to P µ π ≈ 0.05 in the data. For protons the decay channels Λ → pπ − and Λ →pπ + are used and the fake muon probability is found to be P µ p ≈ 0.04 in the data. For K ± mesons from the decay D * ± → D 0 π ± slow → (K ∓ π ± )π ± slow a fake muon probability of P µ K± ≈ 0.01 is measured in the data. It is observed that the pion and proton fake muon probabilities in the data are not described by the Monte Carlo simulation. They are reweighted in the Monte Carlo simulation by factors of 2.0 and 1.9, respectively, to match the data. The K ± fake muon probability in the data and the simulation are in agreement. The misidentified muon events remain in the event sample.

Systematic Uncertainties
The following uncertainties are taken into account in order to evaluate the systematic errors.
• The trigger efficiencies are determined using independent trigger channels in the DIS regime since no independent triggers exist in photoproduction. The uncertainty is estimated by the difference between the efficiency found in the data and the simulation. The Monte Carlo simulation is reweighted to match the data trigger efficiency.
• The efficiency for the identification of the muons is determined using a high statistics sample of events of elastically produced J/ψ mesons [40]. The efficiency is known to a precision of 4%.
• The track efficiency of the CTD is known to ±1% and that of the CST to ±2%. The uncertainty due to the track efficiencies is estimated by varying the efficiencies of the CTD and CST correspondingly.
• The integrated luminosity is known to a precision of 4%.
• The uncertainty due to the resolution of the impact parameter δ of the muon tracks is estimated by varying the resolution by an amount that encompasses any difference between the data and the simulation. This is achieved by applying an additional Gaussian smearing in the Monte Carlo simulation of 200 µm to 5% of randomly selected tracks and 12 µm to the rest.
• The uncertainty on the cross section arising from the uncertainty on the reconstruction of φ jet is estimated by shifting its value by ±2 • .  Table 3: Summary of the systematic uncertainties of the beauty and charm cross sections.
• The uncertainty arising from the hadronic energy scale is estimated by changing it by ±1% for the complete hadronic final state.
• The dependence of the measurement on the physics model used for the templates representing different QCD evolution schemes is estimated by replacing the PYTHIA b and c Monte Carlo templates with CASCADE.
• The uncertainty on the cross section arising from the uncertainty of the parton fragmentation model is estimated by replacing the PYTHIA Monte Carlo samples using the Lund fragmentation function with samples based on the Peterson fragmentation function.
• The uncertainty arising from fake muon background is estimated by not applying the weights that have been found in the fake rate probabilities for K and Λ decays.
• The impact of the reweighting on the cross sections are investigated and found to be negligible.
The effect of the listed experimental uncertainties are estimated by varying the relevant variables in the Monte Carlo simulation or by modifying the corresponding efficiencies in the cross section calculation. The difference between the obtained cross sections with and without the change result in the measurement systematic uncertainties, which are summarised in table 3. The individual effects of the above experimental uncertainties are combined in quadrature, yielding a total uncertainty of 10.5% and 10.4% on the measured b and c cross sections, respectively. The systematic uncertainties as derived from the integrated sample are applied to each analysis bin in order to avoid statistical fluctuations. The largest contribution to this uncertainty for the c measurement arises from systematics attributed to the hadronic energy scale (5%). The systematic errors of the b analysis are not dominated by a single source.

Results
The cross sections for b and c in photoproduction using semi-muonic decays in dijet events are measured. The cross sections are determined for the phase space defined by the kinematic range and the event selection cuts presented in table 2. The measured cross sections are compared to the expectations of the MC programs PYTHIA, CASCADE, HERWIG, and MC@NLO. The total measured and predicted cross sections are listed in table 4. PYTHIA shows the highest normalisation of the three LO MCs, while the normalisation of CASCADE is below the one of PYTHIA and HERWIG has the lowest normalisation. For the beauty measurement, the PYTHIA prediction is closest to the data and gives the best description of the three LO MC predictions. The beauty and charm data cross sections tend to be underestimated by the MC@NLO predictions but are in agreement within the errors. The precision of the measured cross sections are much higher than the ones of the theory predictions shown here.
The beauty and charm cross sections are measured differentially as a function of the transverse momentum of the leading jet P jet1 T and of the muon P µ T , the pseudorapidity of the muon η µ , the momentum fraction x obs γ carried by the photon entering the hard interaction and the azimuthal angular difference ∆φ jj between the two leading jets. The measurements are performed for the full sample, as well as for direct and resolved enriched processes separately. The distinction is performed by the variable x obs γ , which leads to enriched resolved processes in the region x obs γ ≤ 0.75 and direct photon enriched processes for x obs γ > 0.75. The cross sections for the beauty measurements are given in tables 5-7 and shown in figures 3-5.
In the case of beauty production the models provide a good description of the measured cross sections in terms of shape in all distributions. For the LO MCs this is true for the full sample, as well as for the direct and resolved enriched regions. The cross sections for direct enriched processes are well described in shape, but tend to be underestimated by MC@NLO, while for resolved enriched processes a reasonable agreement is observed both in shape and normalisation.
In the analysis of semi-muonic b decays in dijet events with HERA I data [5] an excess of data compared to the NLO predictions of the FMNR program was observed in the first bin of P µ T and P jet1 T . In this analysis the NLO predictions are provided by MC@NLO which is based on the FMNR parton level calculations. Also in this analysis, the NLO predictions lie below the data in the first bin of P µ T and P jet1 T , but they are consistent with the data within 2σ of the experimental and theoretical uncertainty. Whereas in enriched direct processes, the data tend to be underestimated by the MC@NLO prediction in the first bin of P µ T and the first bin of P jet1 T , in enriched resolved processes, no such effect is visible.
The cross sections as a function ∆φ jj show a significant contribution away from the backto-back configuration at ∆φ jj ≃ 180 • . Such a configuration can be described by models which include significant contributions from higher order QCD radiation or a transverse momentum of the gluon in the initial state. This distribution is reasonably well described by all models. In direct and resolved enriched processes, this observation also holds.
The measured charm cross sections are presented in tables 8-10 and figures 6-8. The distributions are reasonably well described by all models. Similar to recent observations in H1 measurements of the photoproduction of D * mesons [41], the central value of the MC@NLO calculations tend to be lower than the measured charm cross sections.

Conclusions
Beauty and charm photoproduction cross sections for events with dijets and a muon are measured using the data collected by the H1 detector at HERA. Compared to the previous H1 beauty measurement [5], the analysis profits from a three times larger luminosity of the data sample, an extended phase space as well as improved understanding of the H1 vertex detector. The flavour composition of the event sample is determined by the transverse momentum of the muon relative to the jet axis of its associated jet P rel T and by its impact parameter δ. Total visible and differential cross sections are measured, and the results are compared to leading order QCD models provided by PYTHIA, CASCADE, and HERWIG, as well as to the next-to-leading order calculations provided by MC@NLO. At low values of P µ T and P jet1 T , the present beauty measurement does not show a significant excess as observed by the previous H1 measurement with respect to the NLO calculation. In general the predictions are in reasonable agreement with the beauty and charm measurements.

H1 Beauty Dijet Muon Cross Sections
σ vis (ep → ebbX → ejjµX) [GeV] [pb/GeV] 2. 5 3  Table 5: Bin averaged differential cross sections for beauty photoproduction of dijet events using semi-muonic decays in bins of P µ T , η µ , P jet1 T , ∆φ jj , and x obs γ with their statistical and systematic uncertainties. The fit parameter f b is given including its statistical error.

H1 Beauty Dijet Muon Cross Sections
[GeV] [pb/GeV] 2. 5 3  Table 6: Bin averaged differential cross sections for beauty photoproduction of dijet events using semi-muonic decays for x obs γ > 0.75 in bins of P µ T , η µ , P jet1 T , and ∆φ jj with their statistical and systematic uncertainties. The fit parameter f b is given including its statistical error.

H1 Beauty Dijet Muon Cross Sections
[GeV] [pb/GeV] 2. 5 3  Table 7: Bin averaged differential cross sections for beauty photoproduction of dijet events using semi-muonic decays for x obs γ ≤ 0.75 in bins of P µ T , η µ , P jet1 T , and ∆φ jj with their statistical and systematic uncertainties. The fit parameter f b is given including its statistical error.

H1 Charm Dijet Muon Cross Sections
σ vis (ep → eccX → ejjµX) [GeV] [pb/GeV] 2. 5 3.  Table 8: Bin averaged differential cross sections for charm photoproduction of dijet events using semi-muonic decays in bins of P µ T , η µ , P jet1 T , ∆φ jj , and x obs γ with their statistical and systematic uncertainties. The fit parameter f c is given including its statistical error.  Table 9: Bin averaged differential cross sections for charm photoproduction of dijet events using semi-muonic decays for x obs γ > 0.75 in bins of P µ T , η µ , P jet1 T , and ∆φ jj with their statistical and systematic uncertainties. The fit parameter f c is given including its statistical error.  Table 10: Bin averaged differential cross sections for charm photoproduction of dijet events using semi-muonic decays for x obs γ ≤ 0.75 in bins of P µ T , η µ , P jet1 T , and ∆φ jj with their statistical and systematic uncertainties. The fit parameter f c is given including its statistical error.   (c) (d) Figure 5: The differential cross sections for beauty photoproduction of dijet events using semimuonic decays for x obs γ ≤ 0.75 as a function of P µ T , η µ , P jet1 T , and ∆φ jj . For details see caption of figure 3.  Figure 6: The differential cross sections for charm photoproduction of dijet events using semimuonic decays as a function of P µ T , η µ , P jet1 T , ∆φ jj , and x obs γ . For details see caption of figure 3.    Figure 8: The differential cross sections for charm photoproduction of dijet events using semimuonic decays for x obs γ ≤ 0.75 as a function of P µ T , η µ , P jet1 T , and ∆φ jj . For details see caption of figure 3.