Search for supersymmetry with photons in pp collisions at s p 8 TeV

Two searches for physics beyond the standard model in events containing photons are presented. The data sample used corresponds to an integrated luminosity of 19 . 7 fb − 1 of proton-proton collisions at ﬃﬃﬃ s p ¼ 8 TeV, collected with the CMS experiment at the CERN LHC. The analyses pursue different inclusive search strategies. One analysis requires at least one photon, at least two jets, and a large amount of transverse momentum imbalance, while the other selects events with at least two photons and at least one jet, and uses the razor variables to search for signal events. The background expected from standard model processes is evaluated mainly from data. The results are interpreted in the context of general gauge-mediated supersymmetry, with the next-to-lightest supersymmetric particle either a bino- or wino-like neutralino, and within simplified model scenarios. Upper limits at the 95% confidence level are obtained for cross sections as functions of the masses of the intermediate supersymmetric particles.


I. INTRODUCTION
Supersymmetry (SUSY) [1-7] is a popular extension of the standard model, which offers a solution to the hierarchy problem [8] by introducing a supersymmetric partner for each standard model particle. In models with conserved R-parity [9,10], as are considered here, SUSY particles are produced in pairs and the lightest supersymmetric particle (LSP) is stable. If the LSP is weakly interacting, it escapes without detection, resulting in events with an imbalancẽ p miss T in transverse momentum. Models of SUSY with gauge-mediated symmetry breaking [11][12][13][14][15][16][17] predict that the gravitino (G) is the LSP. If the next-to-lightest SUSY particle is a neutralino (χ 0 1 ) with a bino or wino component, photons with large transverse momenta (p T ) may be produced inχ 0 1 → γG decays. The event contains jets if theχ 0 1 originates from the cascade decay of a strongly coupled SUSY particle (a squark or a gluino).
In this paper, we present two searches for gaugemediated SUSY particles in proton-proton (pp) collisions: a search for events with at least one isolated high-p T photon and at least two jets, and a search for events with at least two isolated high-p T photons and at least one jet. The discriminating variables are E miss T for the single-photon analysis, and the razor variables M R and R 2 [18,19] for the double-photon analysis, where E miss T is the magnitude of p miss T . These studies are based on a sample of pp collision events collected with the CMS experiment at the CERN LHC at a center-of-mass energy of 8 TeV. The integrated luminosity of the data sample is 19.7 fb −1 .
Searches for new physics with similar signatures were previously reported by the ATLAS and CMS collaborations at ffiffi ffi s p ¼ 7 TeV, using samples of data no larger than around 5 fb −1 [20][21][22][23]. No evidence for a signal was found, and models with production cross sections larger than ≈10 fb −1 were excluded in the context of general gauge-mediation (GGM) SUSY scenarios [24][25][26][27][28][29]. This paper is organized as follows. In Sec. II we describe the CMS detector, in Sec. III the benchmark signal models, and in Sec. IV the part of the event reconstruction strategy that is common to the two analyses. The specific aspects of the single-and double-photon searches are discussed in Secs. V and VI, respectively. The results of the analyses are presented in Sec. VIII. A summary is given in Sec. IX.

II. CMS DETECTOR
The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8 T. Within the superconducting solenoid volume are a silicon pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter (ECAL), and a brass and scintillator hadron calorimeter (HCAL), each composed of a barrel and two end cap sections. Muons are measured in gas-ionization detectors embedded in the steel flux-return yoke outside the solenoid. Extensive forward calorimetry complements the coverage provided by the barrel and end cap detectors.
Events are recorded using a trigger that requires the presence of at least one high-energy photon. This trigger is utilized both for the selection of signal events, and for the selection of control samples used for the background determination. The specific trigger requirements for the two analyses are described below. Corrections are applied to account for trigger inefficiencies, which are evaluated using samples of data collected with orthogonal trigger conditions. A more detailed description of the CMS detector, together with a definition of the coordinate system used and the relevant kinematic variables, can be found in Ref. [30].

III. SUSY BENCHMARK MODELS
The two searches are interpreted in the context of GGM SUSY scenarios [24][25][26][27][28][29], and in terms of simplified model spectra (SMS) scenarios [31][32][33][34] inspired by GGM models. In these scenarios, R-parity is conserved and the LSP is a gravitino with negligible mass. Four models are considered: GGMbino model.-In this model, squarks (q) and gluinos (g) are produced and decay to a final state with jets and a bino-likeχ 0 1 . This production process dominates over electroweak production in the squark-and gluino-mass region accessible to the analyses. Theχ 0 1 mass is set to 375 GeV, leading to aχ 0 1 →Gγ branching fraction of about 80% [26]. The events are examined as a function of the squark and gluino masses. All other SUSY particles have masses set to 5 TeV, which renders them too heavy to participate in the interactions. In most cases, the final state contains two photons, jets, and E miss T . GGMwino model.-This model is similar to the GGMbino model, except that it contains mass-degenerate wino-likeχ 0 1 andχ AE 1 particles instead of a bino-likeχ 0 1 . The common mass of theχ 0 1 andχ AE 1 is set to 375 GeV. The final state contains a γγ, γV, or VV combination in addition to jets and E miss T , where V is a Z or W boson. With aχ 0 1 →Gγ branching fraction of about 28%, approximately 48% of all events contain at least one photon. T5gg model.-This SMS model is based on gluino pair production, with the gluinos undergoing a three-body decay to qqχ 0 1 , followed byχ 0 1 →Gγ. All decays occur with a branching fraction of 100%. The final state contains at least two photons, jets, and E miss T . T5wg model.-This SMS model is also based on gluino pair production, with one gluino undergoing a three-body decay to qqχ 0 1 , followed byχ 0 1 →Gγ, and the other gluino undergoing a three-body decay to qqχ AE 1 , followed bỹ χ AE 1 →GW AE . All decays occur with a branching fraction of 100%. The final state contains at least one photon, jets, and E miss T . Typical Feynman diagrams corresponding to these processes are shown in Fig. 1. Note that for the two GGM models, the events can proceed through the production of gluino-gluino, gluino-squark, or squark-squark pairs.
Signal events for the GGM models are simulated using the PYTHIA 6 [35] event generator. The squark and gluino masses are varied between 400 and 2000 GeV. Eight massdegenerate squarks of different flavor (u, d, s, and c) and chirality (left and right) are considered. The production cross sections are normalized to next-to-leading order (NLO) in quantum chromodynamics, determined using the PROSPINO [36] program, and is dominated by gluinogluino, gluino-squark, and squark-squark production.
The SMS signal events are simulated with the MadGraph 5 [37] Monte Carlo (MC) event generator in association with FIG. 1. Typical Feynman diagrams for the general gauge-mediation model with bino-(top left) and wino-like (top right) neutralino mixing scenarios. Here, theχ 0 1 can decay toGγ orGZ, with the branching fraction dependent on theχ 0 1 mass. The diagrams for the T5gg (bottom left) and T5wg (bottom right) simplified model spectra are also shown. up to two additional partons. The decays of SUSY particles, the parton showers, and the hadronization of partons, are described using the PYTHIA6 program. Matching of the parton shower with the MADGRAPH5 matrix element calculation is performed using the MLM [38] procedure. The gluino pair-production cross section is described to NLO þ NLL accuracy [36,[39][40][41][42], where NLL refers to next-toleading-logarithm calculations. All SUSY particles except the gluino, squark, LSP, andχ 1 states are assumed to be too heavy to participate in the interactions. The NLO þ NLL cross section and the associated theoretical uncertainty [43] are taken as a reference to derive exclusion limits on SUSY particle masses. Gluino masses of 400 (800) to 1600 GeV, andχ 1 masses up to 1575 GeV, are probed in the T5wg (T5gg) model.
For all the signal models, detector effects are simulated through a fast simulation of the CMS experiment [44].

IV. EVENT RECONSTRUCTION
The events selected in this study are required to have at least one high quality reconstructed interaction vertex. The primary vertex is defined as the one with the highest sum of the p 2 T values of the associated tracks. A set of detector-and beam-related noise cleaning algorithms is applied to remove events with spurious signals, which can mimic signal events with high energetic particles or large E miss T [45,46]. Events are reconstructed using the particle-flow algorithm [47,48], which combines information from various detector components to identify all particles in the event. Individual particles are reconstructed and classified in five categories: muons, electrons, photons, charged hadrons, and neutral hadrons. All neutral particles, and charged particles with a track pointing to the primary vertex, are clustered into jets using the anti-k T clustering algorithm [49], as implemented in the Fast Jet package [50], with a distance parameter of 0.5. The momenta of the jets are corrected for the response of the detector and for the effects of multiple interactions in the same bunch crossing (pileup) [51]. Jets are required to satisfy loose quality criteria that remove candidates caused by detector noise.
Photons are reconstructed from clusters of energy in the ECAL [52]. The lateral distribution of the cluster energy is required to be consistent with that expected from a photon, and the energy detected in the HCAL behind the photon shower cannot exceed 5% of the ECAL cluster energy. A veto is applied to photon candidates that match hit patterns consistent with a track in the pixel detector (pixel seeds), to reduce spurious photon candidates originating from electrons. Spurious photon candidates originating from quark and gluon jets are suppressed by requiring each photon candidate to be isolated from other reconstructed particles. In a cone of radius ΔR ≡ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi ðΔηÞ 2 þ ðΔϕÞ 2 p ¼ 0.3 around the candidate's direction, the scalar p T sums of charged hadrons (I π ), neutral hadrons (I n ), and other electromagnetic objects (I γ ) are separately formed, excluding the contribution from the candidate itself. Each momentum sum is corrected for the pileup contribution, computed for each event from the estimated energy density in the ðη; ϕÞ plane. Selected photons are required to be isolated according to criteria imposed on I π , I n , and I γ as defined in Ref. [52].

V. SINGLE-PHOTON SEARCH
The single-photon analysis is based on a trigger requiring the presence of at least one photon candidate with p T ≥ 70 GeV. The trigger also requires H T > 400 GeV, where H T is the scalar sum of jet p T values for jets with p T ≥ 40 GeV and jηj ≤ 3, including photons that are misreconstructed as jets.
In the subsequent analysis, we make use of the variable p Ã T , which is defined by considering the photon candidate and nearby reconstructed particles, clustered as a jet as described in Sec. IV. If a jet (possibly including the photon) is reconstructed within ΔR < 0.2 of the photon candidate and the p T value of the jet is less than 3 times that of the photon candidate itself, it is referred to as the "photon jet." If such a jet is found, p Ã T is defined as the p T value of the photon jet. Otherwise, p Ã T is the p T value of the photon candidate. We require photon candidates to satisfy p Ã T > 110 GeV and jηj < 1. 44. Also, in the subsequent analysis, we make use of the variable H Ã T , defined as for H T in the previous paragraph but including the p Ã T values of all selected photon candidates. The variables p Ã T and H Ã T reduce differences between photon candidates selected with different isolation requirements compared to the unmodified variables p T and H T . We require events to satisfy H Ã T ≥ 500 GeV. The sample of events with isolated photons so selected is referred to as the γ tight sample. The trigger efficiency for the selected events to enter the sample is determined to be 97%, independent of p Ã T and H Ã T . We require at least two jets with p T ≥ 30 GeV and jηj ≤ 2.5. The jets must be separated by ΔR ≥ 0.3 from all photon candidates, to prevent double counting. In addition, the requirement E miss T ≥ 100 GeV is imposed and events with isolated electrons or isolated muons are vetoed. The selection is summarized in Table I. Note that 0.16% of the selected events contain more than one photon candidate.
The relevant sources of background to the single-photon search are: (i) multijet events with large E miss T , originating from the mismeasured momenta of some of the reconstructed jets. This class of events contains both genuine photons and spurious photon candidates from jets. This is by far the largest contribution to the background.
(ii) events with genuine E miss T originating from the leptonic decay of W bosons, both directly produced and originating from top quark decays, which we refer to as electroweak (EW) background. (iii) rare processes with initial-or final-state photon radiation (ISR/FSR), such as γW, γZ (especially γZ → γνν), and γtt production. The kinematic properties of the multijet background are estimated from a control sample of photon candidates with isolation-variable values (I π , I n , I γ ) too large to satisfy the signal photon selection. We refer to these events as the γ loose sample. Photon candidates of this kind typically originate from jets with anomalous fractions of energy deposited in the ECAL. Other than the orthogonal requirement of a γ loose rather than a γ tight candidate, events in this control sample are selected with the same requirements as the γ tight sample, as summarized in Table I. Despite the different isolation requirement, this sample has properties similar to those of the γ tight sample, due to the use of p Ã T rather than photon p T in the definition of the event kinematic variables. Moreover, events in the γ loose control sample are corrected for a residual difference with respect to the γ tight sample in the distributions of p Ã T and hadronic recoil p T , estimated from events with E miss T < 100 GeV. The corrected distribution of a given kinematic property (e.g., E miss T ) for γ loose events provides an estimate of the corresponding distribution for γ tight events. The uncertainty in the correction factors, propagated to the prediction, is fully correlated among bins in the signal region and is treated as a systematic uncertainty in the background yield. The limited statistical precision of the control sample dominates the total systematic uncertainty. Figure 2 (left) shows the E miss T distribution from the γ tight sample and the corresponding prediction from the γ loose sample, for simulated multijet and γ þ jet events. No discrepancy is observed within the quoted uncertainties.
The EW background is characterized by the presence of an electron misidentified as a photon. The kinematic properties of this background are evaluated from a second control sample, denoted the γ pixel sample, defined by requiring at least one pixel seed matching the photon candidate but otherwise using the γ tight selection criteria, as summarized in Table I. Events in the γ pixel sample are weighted by the probability f e→γ for an electron to be Tests of the background estimation method for the single-photon analysis using simulated events in bins of E miss T . The direct simulation of γ tight events is compared to the prediction of the multijet background from simulated γ loose events (left). Simulated events with γ tight originating from generated electrons are compared to the simulated prediction using the EW background method (right). The blue hatched area represents the total uncertainty and is the quadratic sum of the statistical (red vertical bars) and systematic (red hatched area) uncertainties. In the bottom panels, the ratio of the direct simulation to the prediction is shown.
misidentified as a photon, which is measured as a function of the γ candidate p T , the number of tracks associated with the primary vertex, and the number of reconstructed vertices in an event by determining the rate of events with reconstructed eγ pixel and eγ tight combinations in a sample of Z → e þ e − events. The event-by-event misidentification rate is about 1.5%, with a weak dependence on the number of vertices. A systematic uncertainty of 11% is assigned to f e→γ to account for the uncertainty in the shape of the function and for differences between the control sample in which the misidentification rate is calculated and the control sample to which it is applied. The predicted E miss T distribution for the EW background, obtained from a simulated sample of W boson and tt events, is shown in Fig. 2 (right) in comparison with the results from the direct simulation of events with γ tight originating from electrons. The distributions agree within the quoted uncertainties.
The contribution of ISR/FSR background events is estimated from simulation using leading-order results from the MadGraph 5 MC event generator with up to two additional partons, scaled by a factor of 1.50 AE 0.75 including NLO corrections determined with the MCFM [53,54] program.
The measured E miss T spectrum in the γ tight sample is shown in Fig. 3 in comparison with the predicted standard model background. A SUSY signal would appear as an excess at large E miss T above the standard model expectation. Figure 3 includes, as an example, the simulated distribution for a benchmark GGMwino model with a squark mass of 1700 GeV, a gluino mass of 720 GeV, and a total NLO cross section of 0.32 pb.
For purposes of interpretation, we divide the data into six bins of E miss T , indicated in Table II. For each bin, Table II lists the number of observed events, the number of predicted standard model events, the acceptance for the benchmark signal model, and the number of background events introduced by the predicted signal contributions to the control regions, where this latter quantity is normalized to the corresponding signal yield.
No significant excess of events is observed. An exclusion limit on the signal yield is derived at 95% confidence level (CL), using the CL s method [55][56][57]. For a given signal hypothesis, the six E miss T signal regions are combined in a multichannel counting experiment to derive an upper limit on the production cross section. The results, presented in Sec. VIII, account for the possible contribution of signal events to the two control samples, which lowers the effective acceptance by 10%-20% depending on the assumed SUSY mass values.

VI. DOUBLE-PHOTON SEARCH
Events considered for the double-photon search are collected using triggers developed for the discovery of the Higgs boson in diphoton events [58][59][60]. These triggers use complementary kinematic selections: (i) two photons with p T > 18 GeV, where the highest p T photon is required to have p T > 26 GeV, while the diphoton invariant mass is required to be larger than 70 GeV. (ii) two photons with p T > 22 GeV, where the highest p T photon is required to have p T > 36 GeV. In addition, each photon must satisfy at least one of two requirements: a high value of the shower shape variable R 9 [52] or loose calorimetric identification. For the targeted signals, the combination of the two triggers is found to be 99% efficient.
In the subsequent analysis, at least two photon candidates with p T > 22 GeV and jηj < 2.5 are required. Events are selected if the highest p T photon has p T > 30 GeV. Jets must have p T > 40 GeV and jηj < 2.5, with each jet required to lie a distance ΔR > 0.5 from an identified photon. Only events with at least one selected jet are considered.
The background is dominated by multijet events, which mostly consist of events with at least one genuine photon. Due to the requirement of two photons in the event, the EW and ISR/FSR backgrounds are negligible.
The razor variables M R and R 2 [18,19] are used to distinguish a potential signal from background. To evaluate these variables, the selected jets and photons are grouped into two exclusive groups, referred to as "megajets" [19]. The four-momentum of a megajet is computed as the vector sum of the four-momenta of its constituents. Among all possible megajet pairs in an event, we select the pair with the smallest sum of squared invariant masses of the megajets. Although not explicitly required, the two photons are associated with different megajets in more than 80% of the selected signal events.
The variable M R is defined as wherep j i and p j i z are, respectively, the momentum of the ith megajet and the magnitude of its component along the beam axis. The p T imbalance in the event is quantified by the variable M R T , defined as wherep j i T is the transverse component ofp j i . The razor ratio R is defined as For squark pair production in R-parity conserving models in which both squarks decay to a quark and LSP, the M R distribution peaks at M Δ ¼ ðm 2 q − m 2 LSP Þ=mq, where mq ðm LSP Þ is the squark (LSP) mass. Figure 4  demonstrates that M R also peaks for gluino pair production (left) and in the GGMbino model (right).
The (M R , R 2 ) plane is divided into two regions: (i) a signal region with M R > 600 GeV and R 2 > 0.02, and (ii) a control region with M R > 600 GeV and 0.01 < R 2 ≤ 0.02. The control region is defined such that any potential signal contribution to the control region is less than 10% of the expected number of signal events, producing a negligible bias on the background shape determination, corresponding to less than a 2% shift in the predicted number of background events for 20 expected signal events.
The background shape is determined through a maximum likelihood fit of the M R distribution in the data control region, using the empirical template function with fitted parameters k, M 0 R , and n. The best-fit shape is used to describe the M R background distribution in the signal region, fixing the overall normalization to the observed yield in the signal region. This implicitly assumes a negligible contribution of signal events to the overall normalization. We have studied the impact of the resulting bias and found it to be negligible for the expected signal distributions and magnitudes. The covariance matrix derived from the fit in the control region is used to sample an ensemble of alternative M R background shapes. For each bin of the M R distribution, a   6 (color online). Distribution of M R in the double-photon search for a control sample of jets misreconstructed as photons to which a simulated sample of GGMbino events has been added. The squark and gluino masses are respectively set to mq ¼ 1400 GeV and m~g ¼ 1820 GeV, and the production cross section is fixed to σ ¼ 2.7 fb. The signal contribution is shown by the red histogram. The representations of the uncertainty bands, data points, and the information shown in the bottom panel are the same as in Fig. 5. probability distribution for the yield is derived using pseudoexperiments. The uncertainty in each bin is defined by requiring 68% of the pseudoexperiments to be contained within the uncertainty band.
This background prediction method is tested by applying it to a control sample of events in which jets are misidentified as photons, obtained by selecting photon candidates that fail the requirement on the cluster shape or the photon isolation. The remainder of the photon-selection criteria are the same as for the signal sample. In Fig. 5 we show the fit result in the control region (left) and the extrapolation to the signal region (right).
The contribution of the EW and ISR/FSR backgrounds, characterized by genuine E miss T , is evaluated from simulated events and is found to be negligible compared to the systematic uncertainty associated with the multijet background method, and is accordingly ignored.
A signal originating from heavy squarks or gluinos would result in a wide peak in the M R distribution. This is shown in Fig. 6, where a GGMbino signal sample is added to the control sample of jets misreconstructed as photons, and the background prediction method is applied. The contribution of signal events to the control region is negligible and does not alter the background shape of Fig. 5 (left). The signal is visible as a peak at around 2 TeV. Figure 7 (left) shows the result of the fit and the associated uncertainty band, compared to the data in the control region. The fit result is then used to derive the background prediction in the signal region. The comparison of the prediction to the observed data distribution is shown in Fig. 7 (right). No evidence for a signal is found. The largest positive and negative deviations from the predictions are observed for M R ≳ 2.3 TeV and 1.1≲ M R ≲ 1.9 TeV, respectively, each corresponding to a local significance of ≈1.5 standard deviations.

VII. SIGNAL MODEL SYSTEMATIC UNCERTAINTIES
Systematic uncertainties in the description of the signals are listed in Table III. Differences between the simulation and data for the photon reconstruction, identification, and isolation efficiencies are listed as Data/MC photon scale factors. The uncertainty associated with the parton distribution functions (PDF) is estimated using the difference in the acceptance when different sets of PDFs are used [61][62][63][64][65]. Similarly, different sets of PDFs and different choices for the renormalization scales yield different predictions for the expected production cross section.   [66] 2.6 2.6 Jet energy scale corrections [67] 2-3 shape (bin by bin) 2-5 Initial-state radiation 3-5 < 1 Acceptance due to PDF 1-3 1 -3 Signal yield due to PDF and scales 5-20 1-50

VIII. INTERPRETATION OF THE RESULTS
The result of the single-photon analysis is used to extract a limit on the production cross sections of the GGM and SMS models. Comparing the excluded cross section to the corresponding predicted value, a mass limit is derived in the squark versus gluino mass plane. This procedure allows comparisons with previous results [23]. In the SMS, the limits are derived in the gluino versus gaugino mass plane. The resulting cross section upper limits and the corresponding exclusion contours are shown in Fig. 8. Figure 9 shows the excluded mass regions and the cross section upper limits for the GGMbino and T5gg models obtained from the double-photon analysis.
The single-and double-photon analyses are complementary with respect to the event selection and the search strategy. While the former is a multichannel counting experiment based on the absolute prediction of the standard model backgrounds, the latter uses kinematic information about the razor variables to perform a shape analysis. The best individual sensitivity is in the wino-and the bino-like   neutralino mixing scenario, respectively. The doublephoton analysis performs slightly better compared to the single-photon search in the bino scenario, because of the high-H T trigger requirement in the single-photon selection.

IX. SUMMARY
Two searches for gauge-mediated supersymmetry are presented: a search based on events with at least one photon and at least two jets, and a search based on events with at least two photons and at least one jet. The single-photon search characterizes a potential signal as an excess in the tail of the E miss T spectrum beyond 100 GeV, while the double-photon search exploits the razor variables M R and R 2 . These searches are based on pp collision data collected with the CMS experiment at a center-of-mass energy of ffiffi ffi s p ¼ 8 TeV, corresponding to an integrated luminosity of 19.7 fb −1 . No evidence for supersymmetry production is found, and 95% CL upper limits are set on the production cross sections, in the context of simplified models of gauge-mediated supersymmetry breaking and general gauge-mediation (GGM) models. Lower limits from the double-photon razor analysis range beyond 1.3 TeV for the gluino mass and beyond 1.5 TeV for the squark mass for bino-like neutralino mixings in the studied GGM phase space, extending previous limits [23] by up to 300 and 500 GeV, respectively. The limits from the singlephoton analysis for wino-like neutralino mixings range beyond 0.8 TeV for the gluino mass and 1 TeV for the squark mass in the same GGM phase space, extending previous limits by about 100 and 200 GeV. Within the discussed supersymmetry scenarios, these results represent the current most stringent limits.

ACKNOWLEDGMENTS
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