Panel on future challenges in modeling methodology

This panel paper presents the views of six researchers and practitioners of simulation modeling. Collectively we attempt to address a range of key future challenges to modeling methodology. It is hoped that the views of this paper, and the presentations made by the panelists at the 2004 Winter Simulation Conference will raise awareness and stimulate further discussion on the future of modeling methodology in areas such as modeling problems in business applications, human factors and geographically dispersed networks; rapid model development and maintenance; legacy modeling approaches; markup languages; virtual interactive process design and simulation; standards; and Grid computing

A comparison of CMB- and HLA-based approaches to type I interoperability reference model problems for COTS-based distributed simulation

Commercial-off-the-shelf (COTS) simulation packages (CSPs) are software used by many simulation modellers to build and experiment with models of various systems in domains such as manufacturing, health, logistics and commerce. COTS distributed simulation deals with the interoperation of CSPs and their models. Such interoperability has been classified into six interoperability reference models. As part of an on-going standardisation effort, this paper introduces the COTS Simulation Package Emulator, a proposed benchmark that can be used to investigate Type I interoperability problems in COTS distributed simulation. To demonstrate its use, two approaches to this form of interoperability are discussed, an implementation of the CMB conservative algorithm, an example of a so-called “light” approach, and an implementation of the HLA TAR algorithm, an example of a so-called “heavy” approach. Results from experimentation over four federation topologies are presented and it is shown the HLA approach out performs the CMB approach in almost all cases. The paper concludes that the CSPE benchmark is a valid basis from which the most efficient approach to Type I interoperability problems for COTS distributed simulation can be discovered

Facilitating the analysis of a UK national blood service supply chain using distributed simulation

In an attempt to investigate blood unit ordering policies, researchers have created a discrete-event model of the UK National Blood Service (NBS) supply chain in the Southampton area of the UK. The model has been created using Simul8, a commercial-off-the-shelf discrete-event simulation package (CSP). However, as more hospitals were added to the model, it was discovered that the length of time needed to perform a single simulation severely increased. It has been claimed that distributed simulation, a technique that uses the resources of many computers to execute a simulation model, can reduce simulation runtime. Further, an emerging standardized approach exists that supports distributed simulation with CSPs. These CSP Interoperability (CSPI) standards are compatible with the IEEE 1516 standard The High Level Architecture, the defacto interoperability standard for distributed simulation. To investigate if distributed simulation can reduce the execution time of NBS supply chain simulation, this paper presents experiences of creating a distributed version of the CSP Simul8 according to the CSPI/HLA standards. It shows that the distributed version of the simulation does indeed run faster when the model reaches a certain size. Further, we argue that understanding the relationship of model features is key to performance. This is illustrated by experimentation with two different protocols implementations (using Time Advance Request (TAR) and Next Event Request (NER)). Our contribution is therefore the demonstration that distributed simulation is a useful technique in the timely execution of supply chains of this type and that careful analysis of model features can further increase performance

Measurements of Transverse Energy Flow in Deep-Inelastic Scattering at HERA

Measurements of transverse energy flow are presented for neutral current deep-inelastic scattering events produced in positron-proton collisions at HERA. The kinematic range covers squared momentum transfers Q^2 from 3.2 to 2,200 GeV^2, the Bjorken scaling variable x from 8.10^{-5} to 0.11 and the hadronic mass W from 66 to 233 GeV. The transverse energy flow is measured in the hadronic centre of mass frame and is studied as a function of Q^2, x, W and pseudorapidity. A comparison is made with QCD based models. The behaviour of the mean transverse energy in the central pseudorapidity region and an interval corresponding to the photon fragmentation region are analysed as a function of Q^2 and W.Comment: 26 pages, 8 figures, submitted to Eur. Phys.

Searches at HERA for Squarks in R-Parity Violating Supersymmetry

A search for squarks in R-parity violating supersymmetry is performed in e^+p collisions at HERA at a centre of mass energy of 300 GeV, using H1 data corresponding to an integrated luminosity of 37 pb^(-1). The direct production of single squarks of any generation in positron-quark fusion via a Yukawa coupling lambda' is considered, taking into account R-parity violating and conserving decays of the squarks. No significant deviation from the Standard Model expectation is found. The results are interpreted in terms of constraints within the Minimal Supersymmetric Standard Model (MSSM), the constrained MSSM and the minimal Supergravity model, and their sensitivity to the model parameters is studied in detail. For a Yukawa coupling of electromagnetic strength, squark masses below 260 GeV are excluded at 95% confidence level in a large part of the parameter space. For a 100 times smaller coupling strength masses up to 182 GeV are excluded.Comment: 32 pages, 14 figures, 3 table

Measurement of D* Meson Cross Sections at HERA and Determination of the Gluon Density in the Proton using NLO QCD

With the H1 detector at the ep collider HERA, D* meson production cross sections have been measured in deep inelastic scattering with four-momentum transfers Q^2>2 GeV2 and in photoproduction at energies around W(gamma p)~ 88 GeV and 194 GeV. Next-to-Leading Order QCD calculations are found to describe the differential cross sections within theoretical and experimental uncertainties. Using these calculations, the NLO gluon momentum distribution in the proton, x_g g(x_g), has been extracted in the momentum fraction range 7.5x10^{-4}< x_g <4x10^{-2} at average scales mu^2 =25 to 50 GeV2. The gluon momentum fraction x_g has been obtained from the measured kinematics of the scattered electron and the D* meson in the final state. The results compare well with the gluon distribution obtained from the analysis of scaling violations of the proton structure function F_2.Comment: 27 pages, 9 figures, 2 tables, submitted to Nucl. Phys.

Deep-Inelastic Inclusive ep Scattering at Low x and a Determination of alpha_s

A precise measurement of the inclusive deep-inelastic e^+p scattering cross section is reported in the kinematic range 1.5<= Q^2 <=150 GeV^2 and 3*10^(-5)<= x <=0.2. The data were recorded with the H1 detector at HERA in 1996 and 1997, and correspond to an integrated luminosity of 20 pb^(-1). The double differential cross section, from which the proton structure function F_2(x,Q^2) and the longitudinal structure function F_L(x,Q^2) are extracted, is measured with typically 1% statistical and 3% systematic uncertainties. The measured partial derivative (dF_2(x,Q^2)/dln Q^2)_x is observed to rise continuously towards small x for fixed Q^2. The cross section data are combined with published H1 measurements at high Q^2 for a next-to-leading order DGLAP QCD analysis.The H1 data determine the gluon momentum distribution in the range 3*10^(-4)<= x <=0.1 to within an experimental accuracy of about 3% for Q^2 =20 GeV^2. A fit of the H1 measurements and the mu p data of the BCDMS collaboration allows the strong coupling constant alpha_s and the gluon distribution to be simultaneously determined. A value of alpha _s(M_Z^2)=0.1150+-0.0017 (exp) +0.0009-0.0005 (model) is obtained in NLO, with an additional theoretical uncertainty of about +-0.005, mainly due to the uncertainty of the renormalisation scale.Comment: 68 pages, 24 figures and 18 table

Forward pi^0 Production and Associated Transverse Energy Flow in Deep-Inelastic Scattering at HERA

Deep-inelastic positron-proton interactions at low values of Bjorken-x down to x \approx 4.10^-5 which give rise to high transverse momentum pi^0 mesons are studied with the H1 experiment at HERA. The inclusive cross section for pi^0 mesons produced at small angles with respect to the proton remnant (the forward region) is presented as a function of the transverse momentum and energy of the pi^0 and of the four-momentum transfer Q^2 and Bjorken-x. Measurements are also presented of the transverse energy flow in events containing a forward pi^0 meson. Hadronic final state calculations based on QCD models implementing different parton evolution schemes are confronted with the data.Comment: 27 pages, 8 figures and 3 table

Measurement of Leading Proton and Neutron Production in Deep Inelastic Scattering at HERA

Deep--inelastic scattering events with a leading baryon have been detected by the H1 experiment at HERA using a forward proton spectrometer and a forward neutron calorimeter. Semi--inclusive cross sections have been measured in the kinematic region 2 <= Q^2 <= 50 GeV^2, 6.10^-5 <= x <= 6.10^-3 and baryon p_T <= MeV, for events with a final state proton with energy 580 <= E' <= 740 GeV, or a neutron with energy E' >= 160 GeV. The measurements are used to test production models and factorization hypotheses. A Regge model of leading baryon production which consists of pion, pomeron and secondary reggeon exchanges gives an acceptable description of both semi-inclusive cross sections in the region 0.7 <= E'/E_p <= 0.9, where E_p is the proton beam energy. The leading neutron data are used to estimate for the first time the structure function of the pion at small Bjorken--x.Comment: 30 pages, 9 figures, 2 tables, submitted to Eur. Phys.