
The random grid search (RGS) is a simple, but efficient, stochastic algorithm
to find optimal cuts that was developed in the context of the search for the
top quark at Fermilab in the mid1990s. The algorithm, and associated code,
have been enhanced recently with the introduction of two new cut types, one of
which has been successfully used in searches for supersymmetry at the Large
Hadron Collider. The RGS optimization algorithm is described along with the
recent developments, which are illustrated with two examples from particle
physics. One explores the optimization of the selection of vector boson fusion
events in the fourlepton decay mode of the Higgs boson and the other optimizes
SUSY searches using boosted objects and the razor variables.

Particle physics has an ambitious and broad experimental programme for the
coming decades. This programme requires large investments in detector hardware,
either to build new facilities and experiments, or to upgrade existing ones.
Similarly, it requires commensurate investment in the R&D of software to
acquire, manage, process, and analyse the shear amounts of data to be recorded.
In planning for the HLLHC in particular, it is critical that all of the
collaborating stakeholders agree on the software goals and priorities, and that
the efforts complement each other. In this spirit, this white paper describes
the R&D activities required to prepare for this software upgrade.

This note introduces CutLang, a domain specific language that aims to provide
a clear, human readable way to define analyses in high energy particle physics
(HEP) along with an interpretation framework of that language. A proof of
principle (PoP) implementation of the CutLang interpreter, achieved using C++
as a layer over the CERN data analysis framework ROOT, is presently available.
This PoP implementation permits writing HEP analyses in an unobfuscated manner,
as a set of commands in human readable text files, which are interpreted by the
framework at runtime. We describe the main features of [CutLang] and illustrate
its usage with two analysis examples. Initial experience with CutLang has shown
that a justintime interpretation of a human readable HEP specific language is
a practical alternative to analysis writing using compiled languages such as
C++.

CMS has developed a fast detector simulation package, which serves as a fast
and reliable alternative to the detailed GEANT4based (full) simulation, and
enables efficient simulation of large numbers of standard model and new physics
events. Fast simulation becomes particularly important with the current
increase in the LHC luminosity. Here, I will discuss the basic principles
behind the CMS fast simulation framework, and how they are implemented in the
different detector components in order to simulate and reconstruct sufficiently
accurate physics objects for analysis. I will focus on recent developments in
tracking and geometry interface, which improve the flexibility and emulation
performance of the framework, and allow a better synchronization with the full
simulation. I will then show how these developments have led to an improved
agreement of basic analysis objects and event variables between fast and full
simulation.

We present an extensive study of the MSSM parameter space allowing for
general generation mixing in the squark sector. Employing an MCMC algorithm, we
establish the parameter ranges which are allowed with respect to various
experimental and theoretical constraints. Based on this analysis, we propose
benchmark scenarios for future studies. Moreover, we discuss aspects of
signatures at the LHC.

We present an extensive study of nonminimal flavour violation in the squark
sector in the framework of the Minimal Supersymmetric Standard Model. We
investigate the effects of multiple nonvanishing flavourviolating elements in
the squark mass matrices by means of a Markov Chain Monte Carlo scanning
technique and identify parameter combinations that are favoured by both current
data and theoretical constraints. We then detail the resulting distributions of
the flavourconserving and flavourviolating model parameters. Based on this
analysis, we propose a set of benchmark scenarios relevant for future studies
of nonminimal flavour violation in the Minimal Supersymmetric Standard Model.

We present an extensive study of nonminimally flavour violating (NMFV) terms
in the Lagrangian of the Minimal Supersymmetric Standard Model (MSSM). We
impose a variety of theoretical and experimental constraints and perform a
detailed scan of the parameter space by means of a Markov Chain MonteCarlo
(MCMC) setup. This represents the first study of several nonzero
flavourviolating elements within the MSSM. We present the results of the MCMC
scan with a special focus on the flavourviolating parameters. Based on these
results, we define benchmark scenarios for future studies of NMFV effects at
the LHC.

Natural SUSY suggests the existence of light stop quarks, accessible at the
LHC, which are the focus of a dedicated CMS search program. I present two
inclusive CMS searches that look for TeV scale colored sparticles in final
states with jets, btagged jets and missing transverse energy performed using
up to 19.4fb1 of 8TeV LHC protonproton data. No deviation from the Standard
Model was observed in these searches, and the implications for this was shown
for several simplified model scenarios and phenomenological MSSM.

The LHC has started to constrain supersymmetrybreaking parameters by setting
bounds on possible colored particles at the weak scale. Moreover, constraints
from Higgs physics, flavor physics, the anomalous magnetic moment of the muon,
as well as from searches at LEP and the Tevatron have set additional bounds on
these parameters. Renormalization Group Invariants (RGIs) provide a very useful
way of representing the allowed parameter space by making direct connection
with the values of these parameters at the messenger scale. Using a general
approach, based on the pMSSM parametrization of the soft supersymmetrybreaking
parameters, we analyze the current experimental constraints to determine the
probability distributions for the RGIs. As examples of their application, we
use these distributions to analyze the question of Gaugino Mass Unification and
to probabilistically determine the parameters of General and Minimal Gauge
Mediation with arbitrary Higgs mass parameters at the Messenger Scale.

We present the status of dark matter searches performed by the Compact Muon
Solenoid Experiment using 7 TeV pp data collected by the CERN Large Hadron
Collider in 2010 and 2011. The majority of the results shown here were obtained
using 1.1 fb1 of data. We give highlights from analyses searching for
candidates such as WIMPs, gravitinos, axinos and TeV scale particles. All
observations so far were found to be consistent with the Standard Model
predictions. The search results were used to set exclusion limits on various
new physics scenarios.

The interpretation of data in terms of multiparameter models of new physics,
using the Bayesian approach, requires the construction of multiparameter
priors. We propose a construction that uses elements of Bayesian reference
analysis. Our idea is to initiate the chain of inference with the reference
prior for a likelihood function that depends on a single parameter of interest
that is a function of the parameters of the physics model. The reference
posterior density of the parameter of interest induces on the parameter space
of the physics model a class of posterior densities. We propose to continue the
chain of inference with a particular density from this class, namely, the one
for which indistinguishable models are equiprobable and use it as the prior for
subsequent analysis. We illustrate our method by applying it to the constrained
minimal supersymmetric Standard Model and two nonuniversal variants of it.

We describe the application of model inference based on reference priors to
two concrete examples in high energy physics: the determination of the CKM
matrix parameters rhobar and etabar and the determination of the parameters m_0
and m_1/2 in a simplified version of the CMSSM SUSY model. We show how a
1dimensional reference posterior can be mapped to the ndimensional (nD)
parameter space of the given class of models, under a minimal set of conditions
on the nD function. This referencebased function can be used as a prior for
the next iteration of inference, using Bayes' theorem recursively.

A successful implementation of thermal leptogenesis requires the reheat
temperature after inflation T_R to exceed ~2\times 10^9 GeV. Such a high T_R
value typically leads to an overproduction of gravitinos in the early universe,
which will cause conflicts, mainly with BBN constraints. Asaka and Yanagida
(AY) have proposed that these two issues can be reconciled in the context of
the PecceiQuinn augmented MSSM (PQMSSM) if one adopts a mass hierarchy
m(sparticle)>m(gravitino)>m(axino), with m(axino) keV. We calculate the relic
abundance of mixed axion/axino dark matter in the AY scenario, and investigate
under what conditions a value of T_R sufficient for thermal leptogenesis can be
generated. A high value of PQ breaking scale f_a is needed to suppress
overproduction of axinos, while a small vacuum misalignment angle \theta_i is
needed to suppress overproduction of axions. The large value of f_a results in
late decaying neutralinos. To avoid BBN constraints, the AY scenario requires a
low thermal abundance of neutralinos and high values of neutralino mass. We
include entropy production from late decaying saxions, and find the saxion
needs to be typically at least several times heavier than the gravitino. A
viable AY scenario suggests that LHC should discover a spectrum of SUSY
particles consistent with weak scale supergravity; that the apparent neutralino
abundance is low; that a possible axion detection signal (probably with m_axion
in the submicroeV range) should occur, but no direct or indirect signals for
WIMP dark matter should be observed.

Successful implementation of thermal leptogenesis requires reheat
temperatures T_R\agt 2\times 10^9 GeV, in apparent conflict with SUSY models
with TeVscale gravitinos, which require much lower T_R in order to avoid Big
Bang Nucleosynthesis (BBN) constraints. We show that mixed axion/axino dark
matter can reconcile thermal leptogenesis with the gravitino problem in models
with m_{\tG}\agt 30 TeV, a rather high PecceiQuinn breaking scale and an
initial misalignment angle \theta_i < 1. We calculate axion and axino dark
matter production from four sources, and impose BBN constraints on longlived
gravitinos and neutralinos. Moreover, we discuss several SUSY models which
naturally have gravitino masses of the order of tens of TeV. We find a
reconciliation difficult in Yukawaunified SUSY and in AMSB with a winolike
lightest neutralino. However, T_R\sim 10^{10}10^{12} GeV can easily be
achieved in effective SUSY and in models based on mixed modulianomaly
mediation. Consequences of this scenario include: 1. an LHC SUSY discovery
should be consistent with SUSY models with a large gravitino mass, 2. an
apparent neutralino relic abundance \Omega_{\tz_1}h^2\alt 1, 3. no WIMP direct
or indirect detection signals should be found, and 4. the axion mass should be
less than \sim 10^{6} eV, somewhat below the conventional range which is
explored by microwave cavity axion detection experiments.

We investigate the phenomenology of Effective Supersymmetry (ESUSY) models
wherein electroweak gauginos and third generation scalars have masses up to
about 1~TeV while first and second generation scalars lie in the multiTeV
range. Such models ameliorate the SUSY flavor and CP problems via a decoupling
solution, while at the same time maintaining naturalness. In our analysis, we
assume independent GUT scale mass parameters for third and first/second
generation scalars and for the Higgs scalars, in addition to m_{1/2}, \tan\beta
and A_0, and require radiative electroweak symmetry breaking as usual. We
analyse the parameter space which is consistent with current constraints, by
means of a Markov Chain Monte Carlo scan. The lightest MSSM particle (LMP) is
mostly, but not always the lightest neutralino, and moreover, the thermal relic
density of the neutralino LMP is frequently very large. These models may
phenomenologically be perfectly viable if the LMP before nucleosynthesis decays
into the axino plus SM particles. Dark matter is then an axion/axino mixture.
At the LHC, the most important production mechanisms are gluino production (for
m_{1/2} ~<700~GeV) and third generation squark production, while SUSY events
rich in bjets are the hallmark of the ESUSY scenario. We present a set of
ESUSY benchmark points with characteristic features and discuss their LHC
phenomenology.

Simple SUSY GUT models based on the gauge group SO(10) require tb\tau
Yukawa coupling unification, in addition to gauge coupling and matter
unification. The Yukawa coupling unification places strong constraints on the
expected superparticle mass spectrum, with scalar masses \sim 10 TeV while
gluino masses are much lighter: in the 300500 GeV range. The very heavy
squarks suppress negative interference in the q\bar{q}\to\tg\tg cross section,
leading to a large enhancement in production rates. The gluinos decay almost
always via threebody modes into a pair of bquarks, so we expect at least four
bjets per signal event. We investigate the capability of Fermilab Tevatron
collider experiments to detect gluino pair production in Yukawaunified SUSY.
By requiring events with large missing E_T and \ge 2 or 3 tagged bjets, we
find a 5\sigma reach in excess of m_{\tg}\sim 400 GeV for 5 fb^{1} of data.
This range in m_{\tg} is much further than the conventional Tevatron SUSY
reach, and should cut a significant swath through the most favored region of
parameter space for Yukawaunified SUSY models.

We examine the prospects for testing SO(10) Yukawaunified supersymmetric
models during the first year of LHC running at \sqrt{s}= 7 TeV, assuming
integrated luminosity values of 0.1 to 1 fb^1. We consider two cases: the
Higgs splitting (HS) and the Dterm splitting (DR3) models. Each generically
predicts light gluinos and heavy squarks, with an inverted scalar mass
hierarchy. We hence expect large rates for gluino pair production followed by
decays to final states with large bjet multiplicity. For 0.2 fb^1 of
integrated luminosity, we find a 5 sigma discovery reach of m(gluino) ~ 400 GeV
even if missing transverse energy, E_T^miss, is not a viable cut variable, by
examining the multibjet final state. A corroborating signal should stand out
in the oppositesign (OS) dimuon channel in the case of the HS model; the DR3
model will require higher integrated luminosity to yield a signal in the OS
dimuon channel. This region may also be probed by the Tevatron with 510 fb^1
of data, if a corresponding search in the multib+ E_T^miss channel is
performed. With higher integrated luminosities of ~1 fb^1, using E_T^miss plus
a large multiplicity of bjets, LHC should be able to discover Yukawaunified
SUSY with m(gluino) up to about 630 GeV. Thus, the year 1 LHC reach for
Yukawaunified SUSY should be enough to either claim a discovery of the gluino,
or to very nearly rule out this class of models, since higher values of
m(gluino) lead to rather poor Yukawa unification.

Recent renormalization group calculations of the sparticle mass spectrum in
the Minimal Supersymmetric Standard Model (MSSM) show that tb\tau Yukawa
coupling unification at M_{\rm GUT} is possible when the mass spectra follow
the pattern of a radiatively induced inverted scalar mass hierarchy. The
calculation is entirely consistent with expectations from SO(10) SUSY GUT
theories, with one exception: it seems to require MSSM Higgs soft term mass
splitting at M_{\rm GUT}, dubbed "justso Higgs splitting" (HS) in the
literature, which apparently violates the SO(10) gauge symmetry. Here, we
investigate three alternative effects: {\it i}). SO(10) Dterm splitting, {\it
ii}). inclusion of right hand neutrino in the RG calculation, and {\it iii}).
first/third generation scalar mass splitting. By combining all three effects
(the DR3 model), we find tb\tau Yukawa unification at M_{\rm GUT} can be
achieved at the 2.5% level. In the DR3 case, we expect lighter (and possibly
detectable) third generation and heavy Higgs scalars than in the model with HS.
In addition, the light bottom squark in DR3 should be dominantly a right state,
while in the HS model, it is dominantly a left state.

Supersymmetric models with tb\tau Yukawa unification at M_{GUT}
qualitatively predict a sparticle mass spectrum including first and second
generation scalars at the 315 TeV scale, third generation scalars at the
(few) TeV scale and gluinos in the subTeV range. The neutralino relic density
in these models typically turns out to lie far above the measured dark matter
abundance, prompting the suggestion that instead dark matter is composed of an
axion/axino mixture. We explore the axion and thermal and nonthermal axino
dark matter abundance in Yukawaunified SUSY models. We find in this scenario
that {\it i}). rather large values of PecceiQuinn symmetry breaking scale
f_a\sim 10^{12} GeV are favored and {\it ii}). rather large values of GUT scale
scalar masses \sim 1015 TeV allow for the reheat temperature T_R of the
universe to be T_R\agt 10^6 GeV. This allows in turn a solution to the
gravitino/Big Bang Nucleosynthesis problem while also allowing for baryogenesis
via nonthermal leptogenesis. The large scalar masses for Yukawaunified models
are also favored by data on b\to s\gamma and B_s\to \mu^+\mu^ decay. Testable
consequences from this scenario include a variety of robust LHC signatures, a
possible axion detection at axion search experiments, but null results from
direct and indirect WIMP search experiments.

The requirement of tb\tau Yukawa coupling unification is common in simple
grand unified models based on the gauge group SO(10), and it also places a
severe constraint on the expected spectrum of superpartners. For Yukawaunified
models with \mu >0, the spectrum is characterized by three mass scales: {\it
i}). first and second generation scalars in the multiTeV range, {\it ii}).
third generation scalars, \mu and m_A in the fewTeV range and {\it iii}).
gluinos in the \sim 350500 GeV range with chargino masses around 100160 GeV.
In such a scenario, gluino pair production should occur at large rates at the
CERN LHC, followed by gluino threebody decays into neutralinos or charginos.
Discovery of Yukawaunified SUSY at the LHC should hence be possible with only
1 fb^{1} of integrated luminosity, by tagging multijet events with 23
isolated leptons, without relying on missing E_T. A characteristic dilepton
mass edge should easily be apparent above Standard Model background. Combining
dileptons with bjets, along with the gluino pair production cross section
information, should allow for gluino and neutralino mass reconstruction. A
secondary corroborative signal should be visible at higher integrated
luminosity in the W1Z2> 3\ell channel, and should exhibit the same dilepton
mass edge as in the gluino cascade decay signal.

Unification of GUTscale tb\tau Yukawa couplings is a significant feature
of simple SO(10) SUSY GUTs. Here we present the results of a search that used
the Markov Chain Monte Carlo technique to investigate regions of Yukawa
unification and WMAPcompatible dark matter relic density in SO(10)like MSSM
parameter spaces. We mention the possible LHC signatures of Yukawa unified
scenarios and discuss the consequences for dark matter.

Simple supersymmetric grand unified models based on the gauge group SO(10)
require in addition to gauge and matter unification the unification of
tb\tau Yukawa couplings. Yukawa unification, however, only occurs for very
special values of the soft SUSY breaking parameters. We perform a search using
a Markov Chain Monte Carlo (MCMC) technique to investigate model parameters and
sparticle mass spectra which occur in Yukawaunified SUSY models, where we also
require the relic density of neutralino dark matter to saturate the
WMAPmeasured abundance. We find the spectrum is characterizd by three mass
scales: first/second generation scalars in the multiTeV range, third
generation scalars in the TeV range, and gauginos in the \sim 100 GeV range.
Most solutions give far too high a relic abundance of neutralino dark matter.
The dark matter discrepancy can be rectified by 1. allowing for neutralino
decay to axino plus photon, 2. imposing gaugino mass nonuniversality or 3.
imposing generational nonuniversality. In addition, the MCMC approach finds 4.
a compromise solution where scalar masses are not too heavy, and where
neutralino annihilation occurs via the light Higgs h resonance. By imposing
weak scale Higgs soft term boundary conditions, we are also able to generate 5.
low \mu, m_A solutions with neutralino annihilation via a light A resonance,
though these solutions seem to be excluded by CDF/D0 measurements of the B_s\to
\mu^+\mu^ branching fraction. Based on the dual requirements of Yukawa
coupling unification and dark matter relic density, we predict new physics
signals at the LHC from pair production of 350450 GeV gluinos. The events are
characterized by very high bjet multiplicity and a dilepton mass edge around
mz2mz1 \sim 5075 GeV.