\section{Introduction}
The production of jets with large transverse momenta (\pt) in high-energy proton-proton (pp) collisions originates from parton-parton scattering, a process well described by quantum chromodynamics (QCD), the theory of the strong interaction. The cross section is evaluated as the convolution of the partonic cross sections and the parton distribution functions (PDF) in the proton. At the CERN LHC, the inclusive cross section measured for high-\pt jet production~\cite{Aad:2011fc,Chatrchyan:2012bja,Khachatryan:2016wdh} is in good agreement with the predictions of perturbative QCD (pQCD) calculations at next-to-leading order (NLO) accuracy.

Multijet final states allow studies of further features of pQCD. While at leading order (LO) a parton pair (dijet) is produced in a single parton scattering (SPS), additional jets at lower momenta can originate from two other sources. Either they arise from additional gluon radiation from SPS, or they result from double parton scattering (DPS) processes where two different pairs of partons from the two protons collide independently. The SPS processes provide tests of higher-order pQCD calculations as well as of the parton shower evolution. The contributions from DPS processes increase with center-of-mass energies as the gluon density becomes large at low values of longitudinal momentum fraction in the protons. Experimentally, SPS and DPS contributions can be separated by exploiting the different final-state topology of the two processes. Final states arising from SPS exhibit strong azimuthal and \pt correlations among all final jets, while DPS final states predominantly have a back-to-back topology only for each of the independently produced jet pairs.
Measurements of DPS signals have been performed at different collision energies and for different channels~\cite{Abe:1993rv,Abe:1997xk,Abazov:2009gc,Abazov:2014fha,Abazov:2015nnn,Aaij:2012dz,Aaij:2015wpa}. At 7\TeV, exclusive four-jet final states have been measured by CMS~\cite{Chatrchyan:2013qza}, and \PW+dijet production has been studied by ATLAS~\cite{Aad:2013bjm} and CMS~\cite{Chatrchyan:2013xxa}. Various DPS-sensitive final states have also been measured without a direct extraction of the DPS signal by CMS~\cite{Chatrchyan:2013zja,Khachatryan:2015rra} and ATLAS~\cite{Aad:2014rua,Aad:2014kba}. The present study complements the four-jet measurement~\cite{Chatrchyan:2013qza} by selecting events with jets originating from bottom quarks (denoted as ``\PQb jets''). In a four-jet sample, the SPS and DPS contributions can be disentangled by exploiting the differences expected in the angular and momentum correlations of the measured jets, as discussed in Refs.~\cite{Berger:2009cm,Blok:2015rka,Blok:2015afa}. The requirement of \PQb jets allows grouping the four jets into two pairs according to their flavor, and selecting them with lower \pt thresholds than in the untagged case, thereby facilitating the identification of DPS contributions present in the data sample.

This paper presents a measurement of DPS-sensitive observables in heavy-flavor multijet final states. The results are compared to the predictions of various MC event generators using fixed-order NLO matrix elements, and including the contributions of parton showers and multiple parton interactions (MPI). The latter processes are needed, in particular, to describe the softer hadronic production coming from the ``underlying event'' (UE). The MC used implement the DPS component as a high-\pt extension of the modelling of MPI at \pt values of the order of 3--5\GeV~\cite{Khachatryan:2015pea}. The parameters that control the simulation of softer MPI are assumed to be the same for the generation of MPI at higher-\pt scales, \ie, of DPS processes.  This assumption is used for the predictions based on either LO or NLO matrix element calculations. The MC event generators generally simulate MPI starting from the scale corresponding to the hardest parton-parton scattering provided by the matrix element calculation. In LO event generators, such as \PYTHIA and \HERWIGpp, such a scale is the \pt of the partons participating in the hard scattering, while in NLO dijet generators, e.g., \POWHEG, or multijet generators (without NLO virtual corrections), such as \MADGRAPH, the \pt of the additional outgoing partons in the matrix element calculation is also relevant for the definition of the MPI scale.
Comparing the predictions of these generators with DPS-sensitive observables in data is an important step to validate the extrapolation from soft to hard MPI, and thereby the matching of the matrix element calculations to the simulation of the UE.

The paper is organized as follows. In Section 2, a brief detector description is presented along with details of the MC simulations. In Section 3, the event selection and analysis strategy are described, while Section 4 illustrates the corrections applied to the data and the systematic uncertainties that affect the measurement. Section 5 presents the results, which are then summarized in Section 6.

\section{The CMS detector and Monte Carlo simulation}
The central feature of the CMS apparatus is a superconducting solenoid, of 6 m internal diameter and 15 m in length, which provides a magnetic field of 3.8\unit{T}. Charged-particle trajectories are measured using silicon pixel and strip trackers that cover the pseudorapidity region $\abs{\eta}< 2.5$. An electromagnetic crystal calorimeter (ECAL), and a brass/scintillator hadron calorimeter (HCAL) surround the tracking volume and cover the region $\abs{\eta}< 3.0$. A forward quartz-fiber Cherenkov hadron calorimeter extends the coverage to $\abs{\eta} \leq 5.2$. Muons are measured in the range $\abs{\eta} < 2.4$ in gas-ionization detectors embedded in the steel flux-return yoke of the magnet. The CMS experiment uses a two-level trigger system consisting of a level-1 trigger based on custom hardware using signals from the muon detectors and the calorimeters, and a high level trigger (HLT) based on a farm of computers that have access to the full data for each event. A more detailed description of the CMS detector can be found elsewhere~\cite{Chatrchyan:2008aa}.

Samples of multijet events are produced with the following MC event generators:
\begin{itemize}
\item \PYTHIA~6.426~\cite{Sjostrand:2006za}, \PYTHIA~8.185~\cite{Sjostrand:2007gs}, and \HERWIGpp 2.5.0~\cite{Bahr:2008pv}. All of them use LO 2$\to$2 matrix elements. The \PYTHIA~6 and \PYTHIA~8 event generators simulate parton showers ordered in \pt and use the Lund string model~\cite{Andersson:1998tv} for hadronization, while \HERWIGpp assumes parton showers with radiated gluons ordered in emission angle (angular ordering), and uses a cluster fragmentation model~\cite{Webber:1983if} for hadronization. The \PYTHIA and \HERWIGpp samples are generated with transverse momentum of the outgoing partons $\hat{p}_{\mathrm{T}} > 15$\GeV. The contribution of MPI is also simulated in \PYTHIA and \HERWIGpp. The \PYTHIA~6 event generator with tune Z2*~\cite{Chatrchyan:2013gfi} uses a model~\cite{skands2007} where MPI are interleaved with parton showering. Predictions obtained with \PYTHIA~6 and \PYTHIA~8 with the CUETS1 tunes~\cite{Khachatryan:2015pea} are also considered. These use the CTEQ6L1 PDF set~\cite{Pumplin:2002vw} and include an improved set of UE parameters~\cite{Khachatryan:2015pea}. The \HERWIGpp event generator with two tunes to LHC data, UE-EE-3~\cite{Gieseke:2011na} with the MRST LO** PDF set~\cite{Thorne:2009ky,Sherstnev:2007nd} and UE-EE-5-CTEQ6L1~\cite{Seymour:2013qka} with the CTEQ6L1 PDF set, is also used for comparison. The parameters of the hadronization model are determined from LEP data for both \PYTHIA~\cite{Corke:2010yf} and \HERWIGpp~\cite{Gieseke:2011na}.
\item \POWHEG 1.0~\cite{Nason:2004rx,Frixione:2007vw} matched to the \PYTHIA~8 parton showers including a simulation of MPI. The \POWHEG event generator uses NLO dijet matrix elements implemented via 2$\to$2 and 2$\to$3 diagrams. These matrix elements include only LO effects for the four-jet configuration of the present analysis. For the hard-scattering process, the HERAPDF1.5NLO~\cite{Sarkar:2014zua} PDF set is used with a minimum $\hat{p}_\mathrm{T}$ of 5\GeV. The \PQb quarks are treated as massless in the matrix element calculation. The UE provided by \PYTHIA~8 is simulated with the CUETS1 tune, which uses the HERAPDF1.5LO~\cite{Sarkar:2014zua} PDF set and reproduces with very high precision UE and jet observables at various collision energies. Since the \POWHEG predictions contain both real and virtual corrections for the dijet matrix elements, they are used as the reference baseline in the present analysis. Therefore, the full theoretical uncertainty is provided for the \POWHEG simulation, while only the central predictions are provided for the other MC simulations.
\item \MADGRAPH~5.1.5~\cite{Alwall:2011uj} interfaced with \PYTHIA~8. The \MADGRAPH predictions use a LO multijet matrix element with up to four final-state partons, calculated with the CTEQ6L1 PDF, and a simulation of the UE provided by \PYTHIA~8 tune CUETM1~\cite{Khachatryan:2015pea}, which uses the NNPDF2.3LO PDF set~\cite{Ball:2013hta,Ball:2011uy}. The \pt sum of the four partons, $H_\mathrm{T}$, is required to be  $H_\mathrm{T}>50$\GeV, and the \PQb quarks are treated as massless. The matching scale between the matrix element calculations and the parton shower simulation is taken to be 10\GeV, within the \kt-MLM scheme~\cite{Alwall:2007fs}. Underlying event data are well described by this combination of matrix elements plus parton showers with a proper UE tune~\cite{Khachatryan:2015pea}.
\end{itemize}

The detector response is simulated in detail with the \GEANTfour package~\cite{Agostinelli:2002hh}. All simulated samples are processed and reconstructed in the same manner as collision data. The multijet final state can be mimicked by various background sources, such as Drell--Yan and $\PW$ boson production associated to jets, and top-antitop events. The size of these backgrounds is estimated with \PYTHIA~8 and found to be negligible, with a cross section in the measured phase space less than 0.5\% of that for pure QCD multijet events. Therefore, these background sources are neglected in the following.

\section{Event selection}

This analysis uses data from pp collisions at $\sqrt{s} = 7$\TeV recorded with the CMS apparatus in 2010 corresponding to an integrated luminosity of 3\pbinv. The data were collected at low luminosity (${<} 0.2\times10^{33}\percms$), and consequently with low probability of multiple pp interactions in the same bunch crossing (pileup). These running conditions correspond to a fraction of the total integrated luminosity of 36\pbinv collected in 2010.
The mean number of interactions per bunch crossing is around 1.6 for this sample, which results in small pileup effects in the measured distributions. The MC samples are reweighted to the number of interactions in the data in order to match the multiplicity of reconstructed primary vertices.

For the present study, three HLT single-jet trigger sets are analyzed: one with jet \pt threshold of 15\GeV is used for leading jets with  $20 < \pt < 50$\GeV , a second with \pt threshold of 30\GeV for leading jets with  $50 < \pt < 140$\GeV, and a third with \pt threshold of 50\GeV for leading jets with \pt above 140\GeV. In the region $20 < \pt < 80$\GeV, the trigger efficiency is less than 100\%, increasing from 45\% for leading jets with $\pt\approx20\GeV$. A correction is thus applied as a function of the leading jet \pt and $\eta$. For leading jet $\pt > 80$\GeV,  the trigger is fully efficient. The choice of such regions is a compromise between statistics and reliability of the trigger efficiency correction.

The physics objects used in this analysis are particle flow (PF) jets~\cite{CMS:2009nxa}. The PF algorithm~\cite{CMS-PAS-PFT-10-001} combines information from all relevant CMS subdetectors to identify and reconstruct all particle candidates in the event, namely leptons, photons, charged and neutral hadrons. The energy of the muons is obtained from the corresponding track momentum. Charged hadrons are reconstructed from tracks in the tracker. The energy of the electrons is determined from a combination of the track momentum at the main interaction vertex, the corresponding ECAL cluster energy, and the energy sum of all bremsstrahlung photons attached to the track. Photons and neutral hadrons are reconstructed from energy clusters in the ECAL and HCAL, respectively; only clusters far away from the extrapolated position of any track are used.
Jets are reconstructed from the four-momenta of the PF candidates with the anti-\kt algorithm~\cite{Cacciari:2008gp} with a distance parameter of 0.5. A tight quality selection~\cite{CMS-PAS-JME-09-008} is applied to suppress unphysical jets, \ie, jets resulting from noise in the ECAL and/or HCAL. Each jet is required to contain at least two PF candidates, one of which has to be a charged hadron. The jet energy fraction carried by neutral hadrons, photons, muons, and electrons must be less than 90\%. With these criteria, jets are selected with an efficiency greater than 99\% and a misidentification rate (\ie the probability of selecting fake jets, like \eg, originating from leptons or calorimeter noise) smaller than 0.5\% for jet $\pt>20\GeV$. A jet \pt correction is applied to both data and simulation to  account for the nonlinear response of the calorimeters and other instrumental effects. These corrections are based on in situ measurements using dijet, $\gamma$+jet, and \Z{}+jet data samples~\cite{Chatrchyan:2011ds}.

A primary vertex (PV) is identified by a collection of tracks measured in the tracker. If more than one PV is present, the vertex with the highest sum of the squared \pt of the tracks associated to it is selected. The selected vertex is required to be reconstructed from at least five charged-particle tracks and must satisfy a set of quality requirements, including $\abs{z_\mathrm{PV}} < 24\unit{cm}$ and $\rho_{\mathrm{PV}} < 2\unit{cm}$, where $z_{\mathrm{PV}}$ and $\rho_{\mathrm{PV}}$ are the longitudinal and transverse distances of the PV from the nominal interaction point in the CMS detector.

The \PQb jets are identified by using information on the secondary decay vertex of the \PQb hadrons, the impact parameter significance, \ie, the three-dimensional impact parameter divided by its resolution, and the tracks and jet kinematics~\cite{Chatrchyan:2012jua}, through the so-called ``combined secondary vertex'' (CSV) discriminant. A loose selection~\cite{Chatrchyan:2012jua} is used in the \PQb tagging algorithm, which gives a \PQb tagging efficiency on single jets larger than 75\% for jet $\pt>20\GeV$, with a maximum of 85\% at $\pt\approx150\GeV$, as estimated by simulation studies with the \PYTHIA~6 sample. The light-flavor (\PQu, \PQd, \PQs quark or gluon) mistag probability is 20\%, 10\% and 15\% for $\pt \approx20,$ 75 and 300\GeV, respectively, for $\abs{\eta} < 2$, increasing to 35\% for jets in the region $2.0 <\abs{\eta}< 2.4$. This loose selection provides a high-statistics sample, though with relatively few genuine \PQb jets. After requiring the two \PQb tags, the \PQb jet purity, \ie the percentage of selected events where both tagged jets originate from \PQb quarks, is about 12\% for this loose selection. The highest-\pt (leading) \PQb-tagged jet is a genuine \PQb jet in 18\% of the selected events, while the fraction of events where the second-highest-\pt (subleading) \PQb-tagged jet originates from a \PQb quark is about 14\%. There is a high degree of correlation between the purities of the leading and the subleading jets. From simulation studies, about 65\% of the selected events with a true leading \PQb jet also contain a true subleading \PQb jet. The \PQb jet purity of the medium selection for the \PQb tagging algorithm~\cite{Chatrchyan:2012jua} is 58\% for the current analysis. Since the results obtained with the medium selection are consistent with those obtained with the loose selection within the systematic uncertainties, we use the latter results, which have higher statistical accuracy.

The correction for the events with four jets that pass the selection criteria but for which the two \PQb-tagged jets are not genuine \PQb jets is performed through the unfolding procedure employed to obtain stable-particle level distributions (Sec.~4). The amount of this type of background is estimated from the purity of the measured distributions. The measurement of the \PQb jet purity is based on fits of the track counting high efficiency (TCHE) distributions~\cite{Chatrchyan:2012jua} of each \PQb-tagged jet with three different shape templates obtained from MC simulation, corresponding to the TCHE values for light-quark and gluon, charm, and bottom jet flavors. The TCHE discriminant corresponds to the second-highest impact parameter significance among all selected tracks belonging to the considered jet. The \PQb jet purities measured in the data and those in the simulation differ by 2--7\%. Scale factors ($\mathrm{SF}_{\text{\PQb-purity}}$), depending on jet \pt and $\eta$, are applied to the simulation to correct for this difference. By applying $\mathrm{SF}_{\text{\PQb-purity}}$ to the simulated events, the \PQb jet purity of the data sample passing the analysis criteria is consistent with that of the MC. Compatible results are obtained if the CSV discriminant of the \PQb tagged jets is used in the fitting procedure, instead of TCHE distributions. The \PQb jet purity of the selection is estimated in the data separately for leading and subleading \PQb-tagged jets in different bins of \pt and $\eta$.

Additional scale factors ($\mathrm{SF}_{\text{\PQb-tag}}$) are applied to the simulation in order to match the \PQb tagging efficiencies measured in data~\cite{Chatrchyan:2012jua}. They depend on the jet \pt, $\eta$, and flavor, and range between 0.9 and 1.1.

A further reweighting as a function of $\hat{p}_\mathrm{T}$ is applied to the LO generators used for data correction, in order to improve their description of the measured distributions.

Events with at least one PV and at least four jets with $\pt>20\GeV$ are selected for the analysis: two of the four jets are the two \PQb-tagged jets with highest \pt within $\abs{\eta}< 2.4$, while the other two are the remaining highest-\pt jets selected within $\abs{\eta}<4.7$ without any \PQb tagging requirement. If two or more \PQb-tagged jets are present, the two with the highest \pt are taken as the ``b quark jet pair'' (referred to as ``bottom'' hereafter). The ``untagged jet pair'' (referred to as ``light'' hereafter) is taken as the remaining two leading jets. The two different $\eta$ ranges are chosen because the absence of the tracker in the forward region does not allow \PQb jets to be identified for $\abs{\eta}> 2.4$.

About 60\,000 events are left in the data after the offline selection described above. In Fig.~\ref{figControlDistr}, the shapes of the \pt and $\eta$ distributions of the leading \PQb-tagged and the leading untagged jet are compared to predictions of \PYTHIA~6 and \HERWIGpp, before unfolding to the stable-particle level. These shapes are well described by both MCs in the central region and over the whole range of \pt, while there are differences of up to 20--40\% for the most forward pseudorapidities ($\abs{\eta}>$ 3).

\begin{figure*}[htbp]
\centering
    \includegraphics[width=0.49\textwidth]{Figure_001-a.pdf}
    \includegraphics[width=0.49\textwidth]{Figure_001-b.pdf}\\
    \includegraphics[width=0.49\textwidth]{Figure_001-c.pdf}
    \includegraphics[width=0.49\textwidth]{Figure_001-d.pdf}
    \caption{Uncorrected transverse momentum (left) and pseudorapidity (right) distributions of data and simulations (\PYTHIA~6 and \HERWIGpp) for the leading \PQb-tagged (top) and leading untagged (bottom) jets.  Only statistical uncertainties are shown.}
    \label{figControlDistr}
\end{figure*}

Differential cross sections (referred to as ``absolute cross sections'' hereafter) as a function of \pt and $\eta$ of each of the four jets are measured in this analysis. In addition, differential distributions normalized to the total number of selected events (referred to as ``normalized cross sections'') are measured as a function of jet correlation variables very similar to those used in the four-jet analysis of Ref.~\cite{Chatrchyan:2013qza}:
\begin{itemize}
\item the difference in azimuthal angle (in the plane transverse to the beam axis, in radians) between the jets belonging to the light-jet pair:
\begin{equation}
\Delta\phi_\text{light}=\abs{\phi_{\text{light}_1}-\phi_{\text{light}_2}};\\
\end{equation}
\item the balance in \pt of the two light jets:
\begin{equation}
\Delta^\text{rel}_\text{light}\pt = \frac{\abs{\ptvec^{\text{light}_1}+\ptvec^{\text{light}_2}}}{\abs{\ptvec^{\text{light}_1}}+\abs{\ptvec^{\text{light}_2}}};
\end{equation}
\item the azimuthal angle $\Delta$S between the two dijet pairs, defined as:
\begin{equation}
\Delta S=\arccos \left[ \frac{\ptvec(\text{bottom}_1,\text{bottom}_2)\cdot \ptvec(\text{light}_1,\text{light}_2)}{\abs{\ptvec(\text{bottom}_1,\text{bottom}_2)}\cdot \abs{\ptvec(\text{light}_1,\text{light}_2)}} \right],
\end{equation}
\end{itemize}
where $\text{bottom}_1$ ($\text{bottom}_2$) and $\text{light}_1$ ($\text{light}_2$) are the leading (subleading) jets of the bottom and light jet pairs, respectively, and $\ptvec(\text{bottom}_1,\text{bottom}_2)$ and $\ptvec(\text{light}_1,\text{light}_2)$ the momentum vectors of each pair, obtained as the vectorial sum of the momenta of the bottom and light jets, respectively.

Results of the jet correlation observables are presented as distributions normalized to the number of events measured in the selected kinematic region. Such normalized distributions are affected by smaller systematic uncertainties than the absolute cross section measurements.

\section{Corrections and systematic uncertainties}
Particle-level distributions are inferred from the reconstructed data by correcting for selection efficiencies and detector effects. The results are corrected to particle level by applying an iterative unfolding~\cite{D'Agostini:1994zf} as implemented in the \textsc{RooUnfold} package~\cite{Adye:2011gm}. Particles are considered stable if their mean path length $c\tau$ is greater than 10\unit{mm}. MC jets are identified as ``b jets'' at the particle level if a \PQb quark is found within a cone of radius R = $\sqrt{\smash[b]{(\Delta\eta)^2+(\Delta\phi)^2}} = 0.3$ around the jet axis. The background consisting of events with four jets that pass the selection criteria but for which the \PQb-tagged jets are not genuine \PQb jets is corrected for with \PYTHIA~6 tune Z2*, after applying the $\mathrm{SF}_{\text{\PQb-tag}}$ and $\mathrm{SF}_{\text{\PQb-purity}}$ scale factors. The correlation between events selected at the reconstructed and particle levels is then studied by constructing the response matrix. The response matrix quantifies the migration probability between the particle level and reconstructed quantities, as well as the overall reconstruction efficiency. It is obtained for each observable with the \PYTHIA~6 tune Z2* sample. Diagonal terms in the response matrix correspond to particle-level quantities that are reconstructed in the same bin after detector simulation. Off-diagonal terms represent the probability of migration between bins at the particle level and bins at the reconstructed level. As an example, Fig.~\ref{figResponseMatrix} shows the response matrices for the \pt and the $\eta$ of the leading \PQb-tagged and the leading untagged jet. They exhibit a diagonal structure, with off-diagonal terms less than 30--40\%. The bin widths are larger than the detector resolution at each bin.

\begin{figure*}[htbp]
\centering
    \includegraphics[width=0.49\textwidth]{Figure_002-a.pdf}
    \includegraphics[width=0.49\textwidth]{Figure_002-b.pdf}\\
    \includegraphics[width=0.49\textwidth]{Figure_002-c.pdf}
    \includegraphics[width=0.49\textwidth]{Figure_002-d.pdf}
    \caption{Response matrices obtained with the \PYTHIA~6 tune Z2* simulation for the transverse momentum (left) and pseudorapidity (right) of the leading \PQb-tagged (top) and leading untagged (bottom) jets.}
    \label{figResponseMatrix}
\end{figure*}

The response matrix obtained with \PYTHIA~6 is used for the data unfolding. As a cross check, a sample of events generated with \HERWIGpp tune UE-EE-3 is unfolded with the \PYTHIA~6 response matrix. All distributions agree with the generated ones within 9--20\%. The iterative unfolding procedure is regularized by limiting the number of iterations to a certain value for each measured distribution. The optimal number of iterations is determined by minimizing the difference between the distributions measured in the data and the ones obtained by applying backwards the detector effects to the unfolded distributions. The number of iterations ranges between 2 and 4 depending on the observable. As expected, the statistical uncertainties of the unfolded distributions are larger than those of the reconstructed data. The unfolding to particle level includes corrections for jet resolution, flavor misidentification, and pileup effects. The results are presented in the kinematic region defined in Table \ref{tablestable}.

\begin{table}[htbp]
\topcaption{Phase space for the cross section measurement.}
\centering
\begin{scotch}{l l}
At least four jets & $\pt>20$\GeV\\
Two leading \PQb jets & $\abs{\eta} < 2.4$\\
Two leading other jets & $\abs{\eta} < 4.7$\\
\end{scotch}
\label{tablestable}
\end{table}

All significant sources of systematic uncertainties are investigated and the corresponding uncertainty is calculated for each distribution. The total uncertainty is obtained by summing up the individual contributions in quadrature. The following systematic effects are considered:
\begin{description}
\item[Model dependence] the response matrix obtained with \PYTHIA~6 is used for the final correction, and the difference between this and that obtained with \HERWIGpp is taken as a measure of the model dependence of the unfolding, resulting in an uncertainty ranging from 9\% to 20\%.
\item[Jet energy scale (JES)] the momentum of the jets is varied according to the uncertainty associated with the reconstructed \pt~\cite{Chatrchyan:2011ds}. The resulting uncertainty is of the order of 20--25\% (5\%) for the absolute (normalized) cross sections.
\item[Jet energy resolution (JER)] the JER differs for data and simulation by 6--19\%~\cite{Chatrchyan:2011ds} depending on the $\eta$ range, and introduces a systematic uncertainty of 4--8\% in all results.
\item[Pileup reweighting] the effect of the pileup reweighting procedure is evaluated and found to be negligible ($<$ 0.1\%).
\item[$\PB$ tagging scale factor ($\mathrm{SF}_{\text{\PQb-tag}}$)] the values of the scale factors are varied by 10\% for each jet flavor~\cite{Chatrchyan:2012jua}. This variation results in an uncertainty of 15--18\% for absolute cross sections and of 1--2\% for the normalized ones.
\item[$\PB$ jet purity] the \PQb jet purity of the sample is evaluated by fitting separately the TCHE distribution of the leading and of the subleading \PQb-tagged jet in bins of \pt, $\eta$ and $\Delta$S. The difference between the unfolded results when using the $\mathrm{SF}_{\text{\PQb-purity}}$ obtained from the two fits is used as a systematic uncertainty, resulting in values of 10--12\% for the absolute cross sections and 1--2\% for the normalized distributions.
\item[Trigger efficiency] the trigger efficiency correction is varied within its uncertainty and the resulting corrections are applied to the data. These variations result in an uncertainty ranging from 1 to 6\%.
\item[Integrated luminosity] the systematic uncertainty on the luminosity of the 2010 data, affecting the absolute cross sections, is 4\%~\cite{CMS:2011egv}.
\end{description}

The dominant source of uncertainty is the JES, which is considered as correlated among the measured bins. The following aspects of the theoretical uncertainty affecting the \POWHEG predictions are also evaluated:

\begin{description}
\item[PDF uncertainty] the choice of the PDF set influences the theoretical predictions. The uncertainty related to the PDF is determined by generating predictions with various PDF eigenvectors. As central PDF set, the HERAPDF1.5NLO together with the \PYTHIA~6 tune CUETS1 is used.
\item[Scale uncertainty] the default renormalization and the factorization scales ($\mu_\mathrm{R}$ and $\mu_\mathrm{F}$) in the matrix element calculations are chosen to be equal to the leading jet \pt value. The uncertainty related to the $\mu_\mathrm{R}$ and $\mu_\mathrm{F}$ choices is estimated by using \POWHEG interfaced to the UE simulation provided by \PYTHIA~8 tune CUETS1-HERAPDF. Six combinations of the ($\mu_\mathrm{R}/\pt$, $\mu_\mathrm{F}/\pt$) scales: (0.5,0.5), (0.5,1), (1,0.5), (1,2), (2,1), and (2,2), are used. The scale uncertainties are evaluated by taking the envelope of the predictions obtained with the listed scale choices.
\end{description}

A summary of all the systematic effects is given in Table \ref{tableuncert}.

\begin{table*}[htbp]
\topcaption{Systematic and statistical uncertainties affecting the absolute and the normalized cross sections for each measured observable: each source of uncertainty is specified and the value is the average over all the bins of the observable. The 4\% uncertainty from the integrated luminosity is included in the total uncertainty affecting the absolute cross sections. The total uncertainty is obtained by summing the individual experimental uncertainties quadratically. The theoretical uncertainties, listed in the last two columns, affect all the predictions. The systematic uncertainties in the normalized cross sections are smaller than those for the absolute cross sections, since, among others, they are not affected by the migration effects from outside the selected phase space.}
\centering
\cmsTable{
\begin{scotch}{c c c c c c c c c | c c}
 {Measured} & {Model} & {JES} & {JER} & {$\mathrm{SF}_{\text{\PQb-tag}}$} & {$\mathrm{SF}_{\text{\PQb-purity}}$} & {Trigger} & {Stat} & {Total} & {PDF} & {Scale}\\
 {observable} &  &  & & & & {efficiency} & & {incl. int. lumi,} & & \\
\hline\\[0.2ex]
\multicolumn{11}{c}{Absolute cross sections}\\
\hline
{\PQb-tagged jet \pt} & 20\% & 25\% & 4\% & 15\% & 12\% & 6\% & 4\% & 38\% & 10\% & 10\%\\
{Untagged jet \pt} & 10\% & 25\% & 4\% & 15\% & 12\% & 6\% & 4\% & 34\% & 10\% & 10\%\\
{Jet $\abs{\eta}\leq3$} & 10\% & 25\% & 4\% & 15\% & 12\% & 5\% & 4\% & 34\% & 15\% & 10\%\\
{Jet $\abs{\eta}>3$} & 20\% & 35\% & 4\% & 15\% & 12\% & 5\% & 4\% & 45\% & 50\% & 15\% \\
\hline\\[0.2ex]
\multicolumn{11}{c}{Normalized cross sections}\\
\hline
$\Delta\phi^{\text{light}}$ & 13\% & 5\% & 1\% & 2\% & 1\% & 1\% & 4\% & 15\% & 5\% & 2\% \\
$\Delta^{\text{rel}}_{\text{light}}p_\text{T}$ & 13\% & 5\% & 7\% & 2\% & 1\% & 1\% & 4\% & 16\% & 5\% & 2\%\\
$\Delta$S & 20\% & 5\% & 10\% & 2\% & 2\% & 1\% & 4\% & 23\% & 10\% & 2\%\\
\end{scotch}}
\label{tableuncert}
\end{table*}

\section{Results}
The absolute differential cross sections are measured as a function of the jet \pt and $\eta$, along with the normalized cross sections as a function of the jet correlation variables. In Table~\ref{CrossSections}, the cross section is given, and compared to predictions from different event generators at the particle level. The \POWHEG event generator interfaced with \PYTHIA~8 tune CUETS1, referred to in the following as ``\POWHEG'', reproduces the measured cross section best. However, if the MPI simulation is switched off, the same \POWHEG predictions, referred to in the following as ``\POWHEG MPI-off'', underestimate the value of the measured cross section. All predictions are consistent with the data within uncertainties, although \MADGRAPH{}+\PYTHIA~8 tune CUETM1 (``\MADGRAPH'' in the following) tends to underestimate the data, and \PYTHIA~8 to overestimate them.

\begin{table*}[htbp]
\topcaption{Inclusive cross section for $\Pp\Pp \to 2 \PQb + 2 \Pj + \PX$ for jet $\pt>20\GeV$, with \PQb jets within $\abs{\eta}<2.4$, and the other jets within $\abs{\eta}<4.7$. The measurements are compared to the MC predictions.}
\centering
\newcolumntype{x}{D{,}{\,\pm\,}{-1}}
\begin{scotch}{l l x}
\multicolumn{1}{c}{Sample} & \multicolumn{1}{c}{PDF} & \multicolumn{1}{c}{Cross section (nb)}\\
\hline
Data & \NA & \multicolumn{1}{c}{$69\pm3\stat\pm24\syst$}\\
\hline
\POWHEG{}+\PYTHIA~8 tune CUETS1 & HERAPDF1.5 & 65,12\\
\POWHEG{}+\PYTHIA~8 tune CUETS1 MPI off & HERAPDF1.5 & 31,6\\
\hline
\PYTHIA~6 tune Z2* & CTEQ6L1 & 77,15\\
\PYTHIA~6 tune CUETS1 & CTEQ6L1 & 77,15\\
\HERWIGpp tune UE-EE-3 & MRST LO** & 44,8\\
\HERWIGpp tune UE-EE-5C & CTEQ6L1 & 47,9\\
\PYTHIA~8 tune CUETS1 & CTEQ6L1 & 96,18\\
\MADGRAPH+\PYTHIA~8 tune CUETM1 & CTEQ6L1 & 39,7\\
\end{scotch}
\label{CrossSections}
\end{table*}

In Fig.~\ref{figPt}, the absolute differential cross sections as a function of the \pt and $\eta$ of the selected jets are shown compared to predictions from \POWHEG. Figures~\ref{figEta1} and \ref{figEta2} present the same differential cross sections as ratios of theoretical predictions from various MC event generators to the data. The \POWHEG predictions reproduce very well the measurements as a function of \pt and $\eta$ of each jet, in both the central and forward regions. The other MC simulations also describe the data satisfactorily, although \HERWIGpp tune UE-EE-5C and \MADGRAPH are systematically lower than the data.
Similar conclusions about \HERWIGpp and \MADGRAPH have been already drawn for inclusive~\cite{Khachatryan:2015pea} and exclusive four-jet~\cite{Chatrchyan:2013qza} final states.

\begin{figure}[htbp]
\centering
    \includegraphics[width=\cmsSmallerFigWidth]{Figure_003-a.pdf}
    \includegraphics[width=\cmsSmallerFigWidth]{Figure_003-b.pdf}
    \caption{Differential cross sections unfolded to the particle level as a function of the jet transverse momenta \pt (left) and pseudorapidity $\eta$ (right) compared to predictions of \POWHEG{}+\PYTHIA~8 tune CUETS1. Scale factors of $10^8$, $10^6$, and $10^2$ are applied (for clarity) to the measurement of the leading, subleading, and third jet, respectively. The error bars on the data represent the total uncertainties, \ie, statistical and systematic added quadratically. The band represents the theoretical uncertainty due to the choice of the scales and PDFs.}
    \label{figPt}
\end{figure}

\begin{figure}[htbp]
\centering
   \includegraphics[width=\cmsFigWidth]{Figure_004.pdf}
    \caption{Ratios of the absolute cross section predictions of \POWHEG, \MADGRAPH, \PYTHIA~6 (P6), \PYTHIA~8 (P8), and \HERWIGpp over data (unfolded to the particle level) as a function of the jet transverse momenta \pt for each jet. The error bars on the data represent the total uncertainties, \ie, statistical and systematic added quadratically. Data are shown with markers at unity. The band represents the theoretical uncertainty due to the choice of the scales and PDFs (shown only around the \POWHEG ratio for clarity, but affecting all predictions in the same way).}
    \label{figEta1}
\end{figure}

\begin{figure}[htbp]
\centering
    \includegraphics[width=\cmsFigWidth]{Figure_005.pdf}
    \caption{Ratios of the absolute cross section predictions of \POWHEG, \MADGRAPH, \PYTHIA~6 (P6), \PYTHIA~8 (P8), and \HERWIGpp over data (unfolded to the particle level) as a function of the jet pseudorapidity $\eta$ for each jet. The error bars on the data represent the total uncertainties, \ie, statistical and systematic added quadratically. Data are shown with markers at unity. The band represents the theoretical uncertainty due to the choice of the scales and PDFs (shown only around the \POWHEG ratio for clarity, but affecting all predictions in the same way).}
    \label{figEta2}
\end{figure}

Figures \ref{figShapeCorr1}--\ref{figShapeCorr3} show the normalized differential cross sections as a function of the correlation observables, $\Delta\phi^{\text{light}}$, $\Delta^{\text{rel}}_{\text{light}}\pt$, and $\Delta$S. The data are compared to the MC simulations considered previously. In addition, predictions from \POWHEG MPI-off are also shown. All MC simulations that include MPI contributions describe the data well. This is remarkable given that the predictions are based on MPI models tuned to data at softer scales ($\pt \approx3$--5\GeV). Conversely, \POWHEG MPI-off is ruled out by the data, especially at low values of $\Delta^{\text{rel}}_{\text{light}}\pt$ ($<$0.1) and for values of $\Delta$S smaller than 2. This is a clear indication for the need of MPI contributions. The discrepancy between the measurement and the \POWHEG MPI-off predictions goes up to 60\% in the low $\Delta$S region, while for the four-jet events of Ref.~\cite{Chatrchyan:2013qza} the disagreement is of about 40\%. This shows that heavy-flavor multijet production with common jet threshold is more sensitive to a DPS contribution than an untagged four-jet sample with asymmetric \pt thresholds. The fact that the normalized distribution as a function of $\Delta\phi^{\text{light}}$ is also described reasonably well by \POWHEG MPI-off reflects the limited DPS sensitivity of this observable, as already observed for exclusive four-jet final states~\cite{Chatrchyan:2013qza}.

In summary, predictions using LO or NLO dijet matrix elements matched to the simulation of MPI effects reproduce the measured normalized cross sections, whereas those without MPI fail to describe them. This study demonstrates the presence of DPS in the data and confirms the sensitivity to such contributions of the jet correlation variables $\Delta$S and $\Delta^{\text{rel}}_{\text{light}}\pt$.

\begin{figure}[htbp]
\centering
    \includegraphics[width=\cmsSmallerFigWidth]{Figure_006-a.pdf}
    \includegraphics[width=\cmsSmallerFigWidth]{Figure_006-b.pdf}
    \caption{Normalized cross sections unfolded to the particle level as a function of $\Delta\phi^{\text{light}}$, compared to predictions of \POWHEG, \MADGRAPH, \PYTHIA~8 (P8), and \HERWIGpp (\cmsLeft), and of the \POWHEG{}+\PYTHIA~8 tune CUETS1 without MPI (\cmsRight). The lower panels show the ratios of the MC predictions over the data. The error bars on the data represent the total uncertainties, \ie, statistical and systematic added quadratically. Data are shown with markers at unity. The band represents the theoretical uncertainty due to the choice of the scales and PDFs (shown only around the \POWHEG line for clarity, but affecting all predictions in the same way).}
    \label{figShapeCorr1}
\end{figure}

\begin{figure}[htbp]
\centering
    \includegraphics[width=\cmsSmallerFigWidth]{Figure_007-a.pdf}
    \includegraphics[width=\cmsSmallerFigWidth]{Figure_007-b.pdf}
    \caption{Normalized cross sections unfolded to the particle level as a function of $\Delta^{\text{rel}}_{\text{light}}\pt$, compared to predictions of \POWHEG, \MADGRAPH, \PYTHIA~8 (P8), and \HERWIGpp (\cmsLeft), and of the \POWHEG{}+\PYTHIA~8 tune CUETS1 without MPI (\cmsRight). The lower panels show the ratios of the MC predictions over the data. The error bars on the data represent the total uncertainties, \ie, statistical and systematic added quadratically. Data are shown with markers at unity. The band represents the theoretical uncertainty due to the choice of the scales and PDFs (shown only around the \POWHEG line for clarity, but affecting all predictions in the same way).}
    \label{figShapeCorr2}
\end{figure}

\begin{figure}[htbp]
\centering
    \includegraphics[width=\cmsSmallerFigWidth]{Figure_008-a.pdf}
    \includegraphics[width=\cmsSmallerFigWidth]{Figure_008-b.pdf}
    \caption{Normalized cross sections unfolded to the particle level as a function of $\Delta$S, compared to predictions of \POWHEG, \MADGRAPH, \PYTHIA~8 (P8), and \HERWIGpp (\cmsLeft), and of the \POWHEG{}+\PYTHIA~8 tune CUETS1 without MPI (\cmsRight). The lower panels show the ratios of the MC predictions over the data. The error bars on the data represent the total uncertainties, \ie, statistical and systematic added quadratically. Data are shown with markers at unity. The band represents the theoretical uncertainty due to the choice of the scales and PDFs (shown only around the \POWHEG line for clarity, but affecting all predictions in the same way).}
    \label{figShapeCorr3}
\end{figure}

\section{Summary}
A study of events with at least four jets, at least two of which are \PQb jets, in proton-proton collisions at 7\TeV is presented. The data, corresponding to an integrated luminosity of 3\pbinv, were collected with the CMS experiment in 2010. The two \PQb jets must be within pseudorapidity $\abs{\eta}<2.4$, and the two other jets must be within $\abs{\eta}<4.7$. The transverse momenta of all the jets are required to be greater than 20\GeV. The cross section is measured to be $\sigma(\Pp\Pp \to 2 \PQb + 2 \Pj  + \PX) = 69 \pm 3\stat\pm24\syst\unit{nb}$. The differential cross sections as a function of the \pt and $\eta$ of each of the four jets are presented, along with the cross sections as a function of kinematic jet correlation variables. The results are compared to several theoretical predictions with and without contributions from double parton scattering. The models based on leading order or next-to-leading-order dijet matrix element calculations, matched to parton shower and including multiparton interaction (MPI) contributions, describe well the differential cross sections as a function of \pt and $\eta$ in the whole measured region. The differential cross sections as a function of the jet correlation variables are poorly reproduced by models that do not include contributions from MPI. Specifically, the predictions of \POWHEG interfaced with \PYTHIA~8 without the simulation of multiple parton interactions underestimate the cross sections as a function of $\Delta$S and $\Delta^{\text{rel}}_{\text{light}}\pt$ in the regions of the phase space where a double parton scattering (DPS) signal is expected. These results demonstrate, for the first time, the sensitivity of kinematic jet correlation variables, such as $\Delta$S and $\Delta^{\text{rel}}_{\text{light}}\pt$, to DPS processes in multijet final states with heavy-quarks.

\clearpage
\begin{acknowledgments}
\hyphenation{Rachada-pisek}
We congratulate our colleagues in the CERN accelerator departments for the excellent performance of the LHC and thank the technical and administrative staffs at CERN and at other CMS institutes for their contributions to the success of the CMS effort. In addition, we gratefully acknowledge the computing centers and personnel of the Worldwide LHC Computing Grid for delivering so effectively the computing infrastructure essential to our analyses. Finally, we acknowledge the enduring support for the construction and operation of the LHC and the CMS detector provided by the following funding agencies: BMWFW and FWF (Austria); FNRS and FWO (Belgium); CNPq, CAPES, FAPERJ, and FAPESP (Brazil); MES (Bulgaria); CERN; CAS, MoST, and NSFC (China); COLCIENCIAS (Colombia); MSES and CSF (Croatia); RPF (Cyprus); SENESCYT (Ecuador); MoER, ERC IUT and ERDF (Estonia); Academy of Finland, MEC, and HIP (Finland); CEA and CNRS/IN2P3 (France); BMBF, DFG, and HGF (Germany); GSRT (Greece); OTKA and NIH (Hungary); DAE and DST (India); IPM (Iran); SFI (Ireland); INFN (Italy); MSIP and NRF (Republic of Korea); LAS (Lithuania); MOE and UM (Malaysia); BUAP, CINVESTAV, CONACYT, LNS, SEP, and UASLP-FAI (Mexico); MBIE (New Zealand); PAEC (Pakistan); MSHE and NSC (Poland); FCT (Portugal); JINR (Dubna); MON, RosAtom, RAS and RFBR (Russia); MESTD (Serbia); SEIDI and CPAN (Spain); Swiss Funding Agencies (Switzerland); MST (Taipei); ThEPCenter, IPST, STAR and NSTDA (Thailand); TUBITAK and TAEK (Turkey); NASU and SFFR (Ukraine); STFC (United Kingdom); DOE and NSF (USA).

Individuals have received support from the Marie-Curie program and the European Research Council and EPLANET (European Union); the Leventis Foundation; the A. P. Sloan Foundation; the Alexander von Humboldt Foundation; the Belgian Federal Science Policy Office; the Fonds pour la Formation \`a la Recherche dans l'Industrie et dans l'Agriculture (FRIA-Belgium); the Agentschap voor Innovatie door Wetenschap en Technologie (IWT-Belgium); the Ministry of Education, Youth and Sports (MEYS) of the Czech Republic; the Council of Science and Industrial Research, India; the HOMING PLUS program of the Foundation for Polish Science, cofinanced from European Union, Regional Development Fund, the Mobility Plus program of the Ministry of Science and Higher Education, the National Science Center (Poland), contracts Harmonia 2014/14/M/ST2/00428, Opus 2013/11/B/ST2/04202, 2014/13/B/ST2/02543 and 2014/15/B/ST2/03998, Sonata-bis 2012/07/E/ST2/01406; the Thalis and Aristeia programs cofinanced by EU-ESF and the Greek NSRF; the National Priorities Research Program by Qatar National Research Fund; the Programa Clar\'in-COFUND del Principado de Asturias; the Rachadapisek Sompot Fund for Postdoctoral Fellowship, Chulalongkorn University and the Chulalongkorn Academic into Its 2nd Century Project Advancement Project (Thailand); and the Welch Foundation, contract C-1845. \end{acknowledgments}


\section{Introduction}
CERN has around 45~km of underground tunnels containing particle accelerators, beam transfer lines and access tunnels at depths of up to 100~m below the surface \cite{saraiva2013radiation}. During accelerator shutdown periods extensive access is required for human operatives performing maintenance and upgrade activities. A robust emergency lighting system, comprising in excess of 3000 luminaires, exists to ensure safe evacuation from underground areas in the event of an emergency. Currently emergency lighting is provided by legacy lighting systems dating from the original underground accelerator installations (Super Proton Synchrotron and Large Electron Positron Collider) which comprise a mixture of low-pressure sodium and incandescent luminaires. These are now increasingly difficult to maintain due to regulations applicable to lighting systems which limit the supply of replacement lamps, creating the need for the development of a replacement system based on current lighting technology and radiation resistant properties.

\subsection{Emergency lighting systems}
Emergency lighting systems in tunnel accelerators are intended to provide a minimum level of illumination for escape routes, in case of a loss of power to normal luminaires. The purpose of emergency lights is to allow people working in underground areas to escape safely from the tunnels and to facilitate the intervention of emergency services if required.
At CERN, emergency luminaires are supplied from rechargeable Nickle Metal Hydroxide (NiMh) batteries via 48~V DC to 230~V AC inverters. Luminaires are supplied through fire rated cables, with each underground area typically being served by two separate electrical circuits to provide redundancy. In most areas of the Large Hadron Collider (LHC) tunnel, the emergency lighting system makes use of low-pressure sodium discharge lamps retained from the era of the LEP accelerator. Sodium discharge lamps with wire wound inductive ballasts have proved to be highly reliable and radiation resistant throughout the radiation environments created by both the LEP and the LHC \cite{roed2012method}, as they are intrinsically highly tolerant to radiation effects.

In recent years, Light-Emitting Diodes (LED) technology has become an attractive alternative to traditional incandescent, fluorescent and discharge lamps, not only for the general illumination market but also in the field of emergency lighting, thanks to their compact dimensions, instant response and high luminous efficacy. A major role in this transition has been played by national and international regulations, such as the EU directive no. 245/2009, which has led to a ban on certain types of inefficient metal halide, mercury and  high-pressure sodium lamps.
As a combined consequence of legislative pressure and technological progress, the expansion of the LED market appears inevitable while high intensity discharge lamps, including low-pressure sodium lights, are losing market share and ultimately becoming obsolete \cite{de2014solid}.
For this reason, a new generation of LED emergency luminaires is being developed for use in CERN underground facilities.

\subsection{Radiation environment}
\label{subsection:radEnvironment}
The LHC radiation environment is characterised in \cite{roed2012method}, \cite{de2015radiation} and \cite{r2e2016}.
Since the emergency luminaires will eventually be installed on the walls of  all underground areas (tunnels and caverns) of the CERN accelerator complex, they must be radiation hardened to an appropriate level in order to prevent premature failure.
Working within the context of the As Low As Reasonably Achievable (ALARA) principle, the system must be designed to maximise operational lifetime, in order to prevent interventions related to luminaire breakdown and replacement. For this reason, the desired lifetime of emergency lighting luminaires is a minimum of 5 years.

The latest values (2015) of yearly  High Energy Hadrons (HEH) fluence (i.e. fluence of all hadrons with an energy  greater than 20~MeV), 1~MeV equivalent neutron (neq) fluence and  Total Integrated Dose (TID) in the LHC tunnel and its injectors are presented in \cite{de2015radiation}; in addition, this report illustrates that TID values on the tunnel walls are approximately one order of magnitude smaller than close to the corresponding beam lines, as evidenced by from both direct TID measurements and Monte Carlo simulations. Therefore, a suitable target dose for the emergency luminaires under development can be set to 100~kGy, since this value exceeds the expected TID relative to 5 years irradiation in most of the CERN accelerator complex, meeting the requirement of the desired lifetime.

In underground areas, an LED emergency lighting system may also be interlocked to activate only during machine stops (i.e. no beam condition), since no access to these areas is possible for human operators during accelerator runs.%
\footnote{Due to the long warm-up time of low-pressure sodium lamps, the present emergency luminaires in the CERN accelerator complex are always energised in order to avoid conflict with EN 1838 requirements for the rapid activation of emergency lighting systems. Therefore, adoption of LED lighting  provides increased flexibility for energy savings thanks to potential for reduced operational time.}
Hence, LED emergency luminaires in CERN tunnels and caverns will also be studied when not powered, at ambient temperature (for reference, typical LHC tunnel temperature during physics runs is within the range  $17\pm3^{\circ}$C).

\begin{figure}
	\centering
	\includegraphics[scale=0.285]{./figures/MODEL_cropped.pdf}
	\caption{System model of an LED emergency luminaire.}
	\label{fig:model}
\end{figure}

\subsection{LED luminaires}
The block diagram model of a COTS LED emergency luminaire is depicted in Fig. \ref{fig:model}. The system comprises a number of components, including a Power Supply Unit (PSU), which converts the input supply voltage (typically 230~V 50~Hz AC) to the DC value required by the LEDs supplied. Switch-Mode Power Supply (SMPS) topologies, such as buck, boost and flyback converters, are commonly used to supply LEDs \cite{rico2005evaluation} in COTS luminaires. However, these PSUs are not designed to be radiation hard and when tested show very poor performance, causing almost immediate catastrophic failure of LED luminaires when subjected to a typical accelerator radiation environment \cite{adell2002total}. The reliability required for their use in accelerator environments (see, for example, \cite{adell2013radiation} for a review of radiation effects on power converters), can therefore not be guaranteed. Each luminaire may contain one or more high power LED light sources, each with a lumen output of approximately 100~lm or greater. The effect of radiation on LED light output has been studied extensively for optical power levels up to 100~mW \cite{johnston2013radiation}, \cite{jimenez2006proton}, \cite{pearton2016review} however little published research exists for high power LEDs up to 1~W and beyond. The individual LED is generally equipped with a plastic lens to improve the light output distribution (see Fig. \ref{fig:model}); the lens is usually made of PolyMethyl MethAcrylate (PMMA) plastic selected for high light transmission, a property which degrades when exposed to radiation \cite{toh2011pmma}. Depending upon the luminaire design, the LED may also be contained in an IP6X housing protected with glass windows to allow light transmission (see Fig. \ref{fig:model}). Borosilicate retains its optical transmission properties relatively well when irradiated \cite{baydogan2012borosilicate} and has been specified for installations following the tests described in this paper. However, the overall light output of the system may be further improved by careful selection of quartz or high purity fused silica materials \cite{esamats} in the future.

Having developed a simple system level model for the LED luminaires, the failure modes of each component and the corresponding impact on the system performance can be demonstrated (see Table \ref{table:r2e}). Early functional testing on a range of COTS products rapidly indicated that first stage of radiation hardening would be to prevent the catastrophic failure mode of the PSU due to the inclusion of sensitive SMPS. Following this insight, CERN has worked with two luminaire vendors to modify their standard products into radiation hardened versions. The luminaires contain modified PSUs, based on a diode bridge rectifier principle, providing significantly increased radiation hardness, and do not conflict with the requirements of IEC 60598-2-22.

\begin{table}[!h]
	\renewcommand{\arraystretch}{1.3}
	\caption{Radiation effects on luminaire components.}
	\label{table:r2e}
	\centering
	\begin{tabular}{ll}
		\toprule
		Component & Effects of radiation\\ \midrule
		PSU & Catastrophic failure modes for SMPS \\ \midrule
		LED & \begin{tabular}{@{}l@{}}
			Radiation damage (displacement) leading to reduced \\ flux \end{tabular} \\ \midrule
		\begin{tabular}{@{}l@{}} PMMA (plastic \\ lens compound)  \end{tabular}  & \begin{tabular}{@{}l@{}} Free radical formation leading to reduced light \\ transmission \\
			Degraded mechanical properties \end{tabular} \\ \midrule
		Glass & \begin{tabular}{@{}l@{}} Colour centre formation leading to reduced light \\ transmission\end{tabular} \\ \bottomrule
	\end{tabular}
\end{table}

\section{Design response}

\begin{figure}
	\centering
	\subfloat[]{\includegraphics[width=0.88\columnwidth]{./figures/COTSdiode.pdf}\label{fig:COTSdiode}}
	\vspace{1.2cm}
	\subfloat[]{\includegraphics[width=0.97\columnwidth]{./figures/COTSsmps.pdf}\label{fig:SMPS}}
	\caption{(a) Circuit diagram of the power supply using diode bridge rectifier topology. (b) Block diagram of the SMPS used in COTS luminaires under test.}
\end{figure}

\subsection{Preliminary irradiation tests}
\label{sec:preliminary}
Preliminary irradiation tests of two sets of COTS emergency lighting systems have been performed at CNRAD facility \cite{CNRAD2016} at CERN. The goal was to analyse the response to radiation of COTS power supplies in order to identify a suitable basis for future design of a rad-tolerant PSU for LED lighting.  The first set of luminaires under test comprised three lamps with a simple power supply based on the principle of a transformer and a diode bridge rectifier, whose circuit diagram is represented in  Fig. \ref{fig:COTSdiode}. Each luminaire presents a step down transformer (denoted with $T$ on the circuit diagram) supplying two W08M diode bridge rectifiers (manufactured by MIC Electronics, and denoted with $B$), each one feeding two parallel strings of six series connected low power, low current ($<$50~mA) white LEDs, indicated as $L$. Each string is series connected to a 560~$\Omega$ resistance in order to prevent thermal runaway of LEDs. The second batch comprised three luminaires using SMPS converters, whose block diagram is depicted in Fig. \ref{fig:SMPS}. As can be seen, the AC input voltage is first rectified and filtered; an EMI filter is used to mitigate the injection of harmonics into the grid. Then, this stage is followed by a high frequency DC/DC converter composed of an inverter using MOSFET power switches, an high frequency transformer and an AC/DC diode rectifier followed by a filter, to smooth the DC voltage. The power switches of the inverter are controlled via PWM modulation, at a switching frequency of 60~kHz. The control feedback loop is closed by the measurement/detection system, which sends  backs to the control circuitry the DC output voltage, using opto-couplers. In order to prevent overloads on the input and output side, an overcurrent protection system (respectively denoted as O.L.P. and O.I.P. on the input and output side) is embedded into the control system; moreover, and Over Voltage Protection (O.V.P.) on the DC output side is also present.

The irradiation took place into four steps, with cumulative TID of 9, 47.9, 93.7 and 129.3~Gy. The four irradiation steps lasted respectively 230, 900, 1260 and 1095~hours,  separated by intervals with no irradiation of 109, 75 and 105~hours. The irradiation conditions are summarized in Table \ref{table:COTS_PSU}.
	The six luminaires were installed on a test bench (illustrated in Fig. \ref{fig:preliminary_layout}), located in CNRAD test area, with a deported RadMon \cite{masi2012radmon} positioned in the middle of the rack, at about 60~cm from ground, to measure precisely the radiation level of the devices under test. All luminaires were powered while irradiated (active test) and were monitored through a remote rad-hard camera.

\begin{table}[!h]
	\centering
	\caption{Irradiation condition of the COTS luminaires.}
	\label{table:COTS_PSU}
	\renewcommand{\arraystretch}{1.3}
	\begin{tabular}{lccccc}
		\toprule
		& Step 1 & Step 2 & Step 3 & Step 4 & Total \\ \midrule
		TID (Gy) & 9 & 38.9 & 45.8 & 35.6 & 129.3\\
		1 Mev neq fluence $\left(10^{11}\cdot\mbox{cm}^{-2}\right)$ & 0.7 & 3.0 & 3.6 & 2.7 & 10.0  \\
		HEH fluence $\left(10^{11}\cdot\mbox{cm}^{-2}\right)$ & 0.5 & 2.1 & 2.5 & 1.9 & 7.0 \\
		Duration (h)  & 230 & 900 & 1260 & 1095 & 3485 \\
		Time after previous step (h)  & / & 109 & 75 & 105 & /\\ \bottomrule
	\end{tabular}
\end{table}

\begin{figure}
	\centering
	\subfloat[]{\includegraphics[height=7cm]{./figures/layout_cots}%
		\label{fig:preliminary_layout}}
	\hfil
	\subfloat[]{\includegraphics[height=7cm]{./figures/layout_cots_on}%
		\label{fig:preliminary_irrad}}
	\caption{(a) Layout of luminaires during preliminary irradiation tests at CNRAD (CERN); on the left, three luminaires using SMPS, on the right three luminaires with diode bridge rectifiers. (b) Picture from remote camera after 5 hours of irradiation; the three luminaires using diode bridges are on the right, luminaires with SMPS on the left.}
	\label{fig:preliminary}
\end{figure}

The three luminaires using SMPS units failed just after the start of the test, after 1 hour, 4 hours and 18~hours of irradiation, corresponding to TID values of less than 0.1~Gy, 0.1~Gy and 0.4~Gy respectively. These results are summarized in Table \ref{table:SMPS_irrad}.
	Figure \ref{fig:preliminary_irrad} shows the test bench from the remote camera after 5 hours of irradiation, with two of the luminaires using SMPS already failed. The failure of SMPS luminaires can be ascribed to displacement damage effects on the power switches of the high frequency DC/DC converter, on the PWM control circuitry and on opto-couplers.

\begin{table}[!h]
	\centering
	\renewcommand{\arraystretch}{1.5}
	\caption{Irradiation of luminaires with SMPS.}
	\label{table:SMPS_irrad}
	\begin{tabular}{lccc}
		\toprule
		Characteristics & Luminaire 1 & Luminaire 2 & Luminaire 3 \\ \midrule
		\begin{tabular}{@{}l@{}} Time before failure (h) \end{tabular}  & 1 & 4 & 18 \\
		\begin{tabular}{@{}l@{}} TID reached (Gy) \end{tabular}  & $<0.1$ & 0.1 & 0.4 \\
		\begin{tabular}{@{}l@{}} HEH fluence $\left(\mbox{cm}^{-2}\right)$\end{tabular}  & $<5.4\cdot 10^8$ & $5.4\cdot 10^8$ & $2.2\cdot10^9$ \\
		\begin{tabular}{@{}l@{}} 1 MeV neq fluence $\left(\mbox{cm}^{-2}\right)$ \end{tabular}  & $<7.3\cdot10^8$ & $7.3\cdot10^8$ & $2.9\cdot10^9$ \\ \bottomrule
	\end{tabular}
\end{table}

Luminaires using power supplies with diode bridge rectifiers succeeded in completing the irradiation test, working for the full duration without showing any significant degradation in illuminance, as confirmed by visual inspections and luxmeter measurements performed after the end of the irradiation. The total dose reached was of 129~Gy, at a HEH fluence of $7.0\cdot10^{11}\,\mbox{cm}^{-2}$ and 1 MeV neq fluence of $1.0\cdot10^{12}\,\mbox{cm}^{-2}$.

	The irradiation test identified the power supply as the  most critical component in COTS emergency luminaires and confirmed that converters using diode bridge topology offer an higher reliability in radiation environment than SMPS. Following these preliminary results, a prototype of a rad-tolerant PSU using the bridge rectifier topology has been built; this converter includes a GBU8K diode bridge, samples of which have been previously tested for radiation hardness by the LHCb experiment \cite{bager2002lhcb}, up to 1 MeV neq fluence of $4.8\cdot 10^{13}\,\mbox{cm}^{-2}$, corresponding to a TID of 3.3~kGy; GBU8K bridge rectifier makes use of glass passivated junctions.
	The circuit diagram of the prototype power supply is shown in Fig. \ref{fig:proto_diagram}; it includes a 230:6 step-down transformer (denoted with $T_{\scriptsize{\mbox{p}}}$ in Fig. \ref{fig:proto_diagram}) and an electrolytic 220~$\mu$F capacitor (denoted with $C_{\scriptsize{\mbox{p}}}$). The power supply feeds a single Cree XR-E LED, indicated with $L_{\scriptsize{\mbox{p}}}$, producing white light; Cree XR-E LED is in turn composed of a blue InGaN LED and yellow phosphors. A picture of an assembled prototype is shown in Fig. \ref{fig:proto_pic}. 2 PSU prototypes have been tested in CNRAD facility, in a test position which allowed higher TID and fluences than the location of the previous test of COTS luminaires.  Irradiation took place in four steps with cumulated TID of 160.2, 406.1, 760.1 and 1107.6~Gy, with   final cumulated HEH fluence of $1.1\cdot10^{13}\,\mbox{cm}^{-2}$ and cumulated 1 MeV neq fluence of $7.7\cdot10^{12}\,\mbox{cm}^{-2}$. From a comparison with the values presented in \cite{de2015radiation}, the TID reached during this irradiation test corresponds approximately to the dose accumulated on the LHC tunnel walls after 5 years exposure.
	All the PSUs were powered while irradiated (active test).
	Inspections were performed after the last irradiation step, and confirmed that all PSUs were still working. Following this encouraging result, the prototype has been chosen as the basis for a rad-tolerant PSU for LED lighting.

\begin{figure}[!h]
	\centering
	\subfloat[]{\includegraphics[width=0.48\columnwidth]{./figures/PSU_proto.pdf}%
		\label{fig:proto_diagram}}
	\hfil
	\subfloat[]{\includegraphics[width=0.4\columnwidth]{./figures/ProtoPCB}%
		\label{fig:proto_pic}}
	\caption{(a) Circuit diagram of the prototype power supply using GBU8K diode bridge rectifier. (b) Picture of the prototype PSU assembled on a stripboard.}
	\label{fig:proto}
\end{figure}

\begin{table}
	\centering
	\caption{Irradiation condition of the prototype PSU.}
	\label{table:protoPSU}
	\renewcommand{\arraystretch}{1.3}
	\begin{tabular}{lccccc}
		\toprule
		 & Step 1 & Step 2 & Step 3 & Step 4 & Total \\ \midrule
		TID (Gy) & 160.2 & 245.9 & 354 & 347.5 & 1107.6\\
		1 Mev neq fluence $\left(10^{12}\cdot\mbox{cm}^{-2}\right)$ & 1.6 & 2.4 & 3.5 & 3.4 & 11  \\
		HEH fluence $\left(10^{12}\cdot\mbox{cm}^{-2}\right)$ & 1.1 & 1.7 & 2.5 & 2.4 & 7.7 \\ \bottomrule
	\end{tabular}
\end{table}

\subsection{Power Supply Modifications}
Once a suitable power converter topology had been identified, work with luminaire vendors was commenced in order to incorporate the new power converter topology into standard products with optical characteristics suitable for emergency evacuation lighting within underground tunnels. A second vendor with an alternative luminaire design more suited to underground caverns and larger spaces was also approached. The power converter specification was issued to both vendors, which included the same bridge diode rectifier (GBU8K) used in the prototype PSU presented in Section \ref{sec:preliminary}. As a result of these initiatives, two types of radiation hardened emergency luminaire were produced for further radiation testing.
The circuit diagrams of Vendor 1 and Vendor 2 power supplies are shown respectively in Fig. \ref{fig:circuitdiagram1} and Fig. \ref{fig:circuitdiagram2}. The values of the components are listed in Table \ref{table:components}.

\subsection{Components under test}
The luminaires from Vendor 1 utilise three Cree XP-G LEDs as light sources, with a combination of PMMA lenses and glass windows. Vendor 2 uses a single Cree XR-E LED with a PMMA lens to direct the light output. Cree XR-E and XP-G LEDs produce white light; they are composed of blue LEDs (InGaN epitaxial layer on silicon substrate) with phosphors producing yellow light. To the best of our knowledge, no previous irradiation tests of Cree XR-E and XP-G LEDs have been reported. Both luminaires use the GBU8K diode manufactured by Vishay.  Further components include surface mounted multi-layer ceramic capacitors and metal-film resistors.

\begin{figure}[!h]
	\centering
	\subfloat[]{\includegraphics[width=0.57\columnwidth]{./figures/PSU_vendor1_C}%
		\label{fig:circuitdiagram1}}
	\hfil
	\subfloat[]{\includegraphics[width=0.57\columnwidth]{./figures/PSU_vendor2_C2}%
		\label{fig:circuitdiagram2}}
	\caption{(a) Circuit diagram of Vendor 1 power supply. (b) Circuit diagram of Vendor 2 power supply.}
	\label{fig:circuitdiagrams}
\end{figure}

\begin{table}[!h]
	\centering
	\renewcommand{\arraystretch}{1.3}
	\caption{Electrical components of Vendor 1 and Vendor 2 power supplies.}
	\label{table:components}
	\begin{tabular}{lllll}
		\toprule
		Name & Quantity & & Name & Quantity \\\midrule
		$T_1$ &  16~VA transformer, 230$\,:\,$14 &\qquad\qquad$\;\;$ & $C$ &  1~mF ceramic capacitor \\
		$T_2$ &  16~VA transformer, 230$\,:\,$7 & & 	$L_1$ &  Cree XR-E white LED \\
		$R_1$ &  16~$\Omega$ 5\% resistor, 25~W & & $L_2$ &  Cree XP-G white LED \\
		$R_2$ &  5~$\Omega$ 5\% resistor, 25~W  & & $B$ & GBU8K diode bridge rectifier \\\bottomrule
	\end{tabular}
\end{table}

\section{Experimental results}

\subsection{Test methodology}
Gamma ray and neutron irradiation was performed on luminaires from both vendors. Gamma ray irradiation was conducted at the Fraunhofer Institute INS (Germany) to a TID of 100~kGy; two sample luminaires for each vendor were sent to Fraunhofer for test purposes. Neutron irradiation was performed at CEA in France on another set of samples, likewise composed of two luminaires from Vendor 1 and two luminaires from Vendor 2. Analysis of the neutron irradiated samples will be presented in a future work.
In addition to irradiation tests, a trial installation using the proposed LED luminaires has recently been completed in part of the SPS accelerator complex (TDC2 and TCC2 tunnels and target caverns) with a very challenging radiation environment \cite{Strabel2015}, with half of the luminaires permanently on and half interlocked to activate only when there is no beam present in the accelerator.

The gamma ray irradiation has been performed at the BGS low dose test facility at Fraunhofer, using a Cobalt-60 source at a dose rate of 720~Gy(Si)/h, to a final dose value of 100~kGy(Si).
The luminaires were not powered during irradiation (passive test). As detailed in Section \ref{subsection:radEnvironment}, in CERN tunnels LED luminaires may be powered off while irradiated; for this reason, passive irradiation tests are also representative of actual radiation effects in accelerator tunnel environment.
The irradiation was conducted in five steps, with final dose per exposure of 4.1, 15.3, 34.2, 50.0 and 100.0~kGy(Si).
After each step (and within three hours after the completion of the exposure),  preliminary electrical and optical measurements were taken to verify the continued functioning of the luminaire.
For each step, this activity took less than two hours.
The irradiation of the sample luminaires lasted overall 146.7 hours. The time interval between the end of a step and the beginning of the following one was within 5 and 20 hours. The salient data of the gamma ray irradiation test are summarized in Table \ref{table:irradiation}.
\begin{table}
	\renewcommand{\arraystretch}{1.3}
	\caption{Summary of gamma ray irradiation conditions.}
	\label{table:irradiation}
	\centering
	\begin{tabular}{ll}
		\toprule
		Parameter & Value\\ \midrule
		Irradiation facility & BGS low dose \\
		Irradiation source & Co-60 gamma \\
		Dose steps  & 4.1, 15.3, 34.2, 50.0, 100.0 kGy(Si) \\
		Dose rate &  720 Gy(Si)/h \\
		Time of irradiation & 146.7 h\\
		Time between steps & 9, 8, 5, 20 hours \\ \bottomrule
	\end{tabular}
\end{table}

\subsection{Results from goniophotometer analysis}
Following the gamma irradiation exposure, the luminaires were sent for a Type C Moving Mirror goniophotometer analysis at TSI-LUX Ltd. in the UK.
Electrical and photometric measurements have been performed and the results are collected in Table \ref{table:lux}.

\begin{table}[!h]
	\renewcommand{\arraystretch}{1.3}
	\caption{Summary of electrical and photometric measurements.}
	\label{table:lux}
	\centering
	\begin{tabular}{lcc}
		\toprule
		Quantity & Vendor 1 & Vendor 2\\ \midrule
		Frequency & 50 Hz& 50 Hz \\
		Voltage & 230.06 V & 230.04 V \\
		Current & 0.031 A & 0.031 A\\
		Power & 6.50 W & 6.20 W \\
		Power factor & 0.924 & 0.879 \\
		Apparent power & 7.04 VA & 7.09 VA \\
		CRI (Ra) &  71 & 65 \\
		Luminous flux pre-irradiation & 900 lm & 209 lm \\
		Luminous flux post-irradiation & 177 lm & 55 lm \\
		\begin{tabular}{@{}l@{}} Variation in luminous flux \end{tabular}  & -80.33\% & -73.68\% \\
		Luminous efficacy pre-irradiation & 128.57 lm/W & 32.15 lm/W\\
		Luminous efficacy post-irradiation & 27.23 lm/W & 8.87 lm/W\\
		\begin{tabular}{@{}l@{}} Variation in luminous efficacy \end{tabular}  & -78.82\% & -72.41\%\\ \bottomrule
	\end{tabular}
\end{table}

\begin{figure}[!b]
	\centering
	\subfloat[]{\includegraphics[height=5cm]{./figures/Vendor1_detail_4kgy}\label{Fig:vendor1glass4kGy}}
	\hfil
	\subfloat[]{\includegraphics[height=5cm]{./figures/vendor1glasson}\label{Fig:vendor1glasson}}
	\caption{(a) Detail of Vendor 1 luminaire glasses after 4.1~kGy irradiation. (b) Vendor 1 luminaire on, after 50~kGy irradiation.}
\end{figure}

\subsection{Analysis of results}
From a comparison of the luminous flux before and after irradiation, an overall reduction of 80.3\% has been observed for Vendor 1 luminaires. This reduction appears to be predominantly due to the glass, which became dark just after the 4.1~kGy exposure (see Fig. \ref{Fig:vendor1glass4kGy}), and the PMMA lenses of the LEDs, which turned yellow. Figure \ref{Fig:vendor1glasson} shows the luminaire switched on, emitting a yellowish light, after a dose of 50~kGy. Finally, Figures \ref{fig:vend1_step1}-\ref{fig:vend1_step5} show Vendor 1 luminaires after every irradiation step.
As for Vendor 2 product, the drop in the luminous flux is slightly lower than Vendor 1, thanks to the absence of the glass, and can be ascribed to the degradation in the PMMA which started to become yellow after a dose of 15~kGy. Figures \ref{fig:vend2step1}-\ref{fig:vend2step5} illustrate Vendor 2 luminaires after every irradiation step.
The variations in luminous efficacy are close to those in luminous flux for both products, meaning that no significant variation in input power (due to LEDs and diodes degradation) has been observed.

\begin{figure}
	\centering
	\subfloat[]{\includegraphics[width=0.48\columnwidth]{./figures/Vendor1_step1}\label{fig:vend1_step1}}
	\hfil
	\subfloat[]{\includegraphics[width=0.48\columnwidth]{./figures/Vendor1_step2}\label{fig:vend1_step2}}
	\hfil
	\subfloat[]{\includegraphics[width=0.48\columnwidth]{./figures/Vendor1_step3}\label{fig:vend1_step3}}
	\hfil
	\subfloat[]{\includegraphics[width=0.48\columnwidth]{./figures/Vendor1_step4}\label{fig:vend1_step4}}
	\hfil
	\subfloat[]{\includegraphics[width=0.48\columnwidth]{./figures/Vendor1_step5}\label{fig:vend1_step5}}
	\caption{Comparison between irradiated and non-irradiated Vendor 1 luminaires after exposure to: (a) 4.1~kGy, (b) 15.3~kGy, (c) 34.2~kGy, (d) 50~kGy, (e) 100~kGy. In Figure (b) both irradiated luminaires are shown.}
\end{figure}

\begin{figure}
	\centering
	\subfloat[]{\includegraphics[width=0.48\columnwidth, angle=180, origin=c]{./figures/Vendor2_step1}\label{fig:vend2step1}}
	\hfil
	\subfloat[]{\includegraphics[width=0.48\columnwidth]{./figures/Vendor2_step2}\label{fig:vend2step2}}
	\hfil
	\subfloat[]{\includegraphics[width=0.48\columnwidth]{./figures/Vendor2_step3}\label{fig:vend2step3}}
	\hfil
	\subfloat[]{\includegraphics[width=0.48\columnwidth]{./figures/Vendor2_step4}\label{fig:vend2step4}}
	\hfil
	\subfloat[]{\includegraphics[width=0.48\columnwidth, angle=180, origin=c]{./figures/Vendor2_step5}\label{fig:vend2step5}}
	\caption{Comparison between an irradiated (white) and non-irradiated (orange) Vendor 2 luminaire after exposure to: (a) 4.1~kGy, (b) 15.3~kGy, (c) 34.2~kGy, (d) 50~kGy, (e) 100~kGy.}
\end{figure}

The data from the photometric measurement has been integrated into a DIALux model for a lighting simulation of a typical segment of the LHC tunnel. Fixing centres of the luminaires are designed to account for the reduction in luminous flux due to radiation. In the simulation, Vendor 1 luminaires are placed on the tunnel wall at fixing centres of 14~m, while Vendor 2 luminaires are positioned on the ceiling of the tunnel at 10~m fixing centres. The simulation results of Vendor 1 product are summarized in Table \ref{table:dialuxVendor1and2}
and a 3D rendering in false colour is shown in Fig. \ref{Fig:visualafter}.
We recall that EN 1838 standard requires a minimum of 0.5 lx for anti-panic lighting and 1 lx for escape routes (at centre line). For this reason, the plastic and glass types for future luminaires will be carefully selected to improve optical transmission when irradiated. The same analysis has been performed for Vendor 2 luminaire, and results are shown in Table \ref{table:dialuxVendor1and2} and Fig. \ref{Fig:visualafter_vend2}. Figure \ref{Fig:visualnormal} illustrates the natural visual impression of the tunnel segment using luminaires of both vendors after 100~kGy irradiation.

\begin{figure}[!t]
	\centering
	\subfloat[]{\includegraphics[width=0.65\columnwidth]{./figures/beforeafterlegendVendor1_4_right2}\label{Fig:visualafter}}
	\hfil
	\subfloat[]{\includegraphics[width=0.65\columnwidth]{./figures/beforeafterlegendVendor2_3_rightlegend}\label{Fig:visualafter_vend2}}
	\hfil
	\subfloat[]{\includegraphics[width=0.65\columnwidth]{./figures/NormalIrradiatedVision4}\label{Fig:visualnormal}}
	\caption{(a) 3D rendering of the DIALux simulation with Vendor 1 luminaires, before (left) and after (right) 100~kGy exposure. (b) 3D rendering of the DIALux simulation with Vendor 2 luminaires, before (left) and after (right) 100~kGy exposure. (c) Natural visual impression using Vendor 1 (left) and Vendor 2 (right) luminaires after 100~kGy exposure.}
\end{figure}

\begin{table}[!h]
	\renewcommand{\arraystretch}{1.3}
	\caption{DIALux simulations of Vendor 1 and Vendor 2 luminaires in LHC tunnel.}
	\label{table:dialuxVendor1and2}
	\centering
	\begin{tabular}{lcccccc}
		\toprule
		\multirow{2}{*}{Quantity} & \qquad & \multicolumn{2}{c}{Vendor 1}  & \qquad &\multicolumn{2}{c}{Vendor 2} \\
			& \qquad & Before & After  & \qquad &  Before & After \\ \midrule
		Average over the whole surface & \qquad & 9.78~lx & 2.01~lx & \qquad & 3.46~lx & 1.60~lx \\
		Maximum value  over the whole surface & \qquad & 17.0~lx & 9.45~lx & \qquad& 6.01~lx & 2.78~lx\\
		Minimum value along the escape route & \qquad & 5.67~lx & 0.65~lx & \qquad & 1.85~lx & 0.92~lx \\
		Minimum value at center line along escape route & \qquad & 7.48~lx & 1.22~lx & \qquad &  2.02~lx & 1.02~lx  \\ \bottomrule
	\end{tabular}
\end{table}

\section{Future developments}
The need for a new PSU as the first stage in radiation hardening COTS luminaires, together with impending obsolescence of the GBU8K rectifier, has spurred the production of a reference design, released under the CERN Open Hardware License \cite{devine2016psu}. Future tests up to 100~kGy are planned for the new PSU in CHARM facility at CERN \cite{ThorntonCHARM2016} using mixed field radiation produced from 24 GeV protons, investigating also the effect of different load conditions. Dedicated irradiation tests for the power LEDs are scheduled at IRRAD facility at CERN \cite{RavottiIRRAD2014}, so as to identify the impact of displacement damage using 24 GeV protons, under different bias conditions. Likewise, gamma ray irradiation tests for samples of optical components (PMMA lenses and borosilicate windows) are planned. Future studies will address the analysis of photometric data from the neutron irradiated sample luminaires, which will provide a stronger indication on the effect of displacement damage to the LEDs, and the consequent impact on overall system performance. For the purposes of this paper, the effects of annealing within the LED die have been ignored due to a significant delay between irradiation and photometric measurement, however the impact of this phenomena would also be an interesting subject for future study, within the context of the duty cycle for the luminaires and operational practices of particle accelerators.

\section{Conclusions}
In this paper, we presented progress towards radiation hardening of LED emergency luminaires for use in accelerators environments.
After the identification of the PSU as the cause of early catastrophic failures in COTS luminaires, we described a simple AC/DC power supply for LED emergency luminaires, the inclusion of which allows modified COTS luminaires to remain functional up to a dose of 100~kGy. The proposed luminaires are now being installed within limited areas of the CERN accelerator complex. Future TID tests are planned in order to bound the lifetime of the whole LED luminaire; moreover, detailed study of the principal optical components, (power LEDs, plastic lenses and glass) are scheduled. Future studies identified include active and non-active irradiation tests of the reference PSU design up to 100~kGy.

\section*{Acknowledgments}
The authors would like to thank Jean-Marie Foray for commencing the LED luminaire testing programme during his time at CERN from 2009-2014.

\bibliographystyle{ieeetr}