Page Menu
Home
HEPForge
Search
Configure Global Search
Log In
Files
F8308787
CDF_2001_S4517016.cc
No One
Temporary
Actions
View File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Flag For Later
Size
2 KB
Subscribers
None
CDF_2001_S4517016.cc
View Options
// -*- C++ -*-
#include
"Rivet/Analysis.hh"
#include
"Rivet/RivetAIDA.hh"
#include
"Rivet/Tools/Logging.hh"
#include
"Rivet/Tools/BinnedHistogram.hh"
#include
"Rivet/Projections/FinalState.hh"
#include
"Rivet/Projections/FastJets.hh"
namespace
Rivet
{
/// @brief CDF two jet tripply-differential cross-section
class
CDF_2001_S4517016
:
public
Analysis
{
public
:
/// @name Constructors etc.
//@{
/// Constructor
CDF_2001_S4517016
()
:
Analysis
(
"CDF_2001_S4517016"
)
{
setBeams
(
PROTON
,
ANTIPROTON
);
setNeedsCrossSection
(
true
);
}
//@}
public
:
/// @name Analysis methods
//@{
/// Book histograms and initialise projections before the run
void
init
()
{
FinalState
fs
(
-
4.2
,
4.2
);
addProjection
(
FastJets
(
fs
,
FastJets
::
CDFJETCLU
,
0.7
),
"Jets"
);
_h_ET
.
addHistogram
(
0.1
,
0.7
,
bookHistogram1D
(
1
,
1
,
1
));
_h_ET
.
addHistogram
(
0.7
,
1.4
,
bookHistogram1D
(
2
,
1
,
1
));
_h_ET
.
addHistogram
(
1.4
,
2.1
,
bookHistogram1D
(
3
,
1
,
1
));
_h_ET
.
addHistogram
(
2.1
,
3.0
,
bookHistogram1D
(
4
,
1
,
1
));
}
/// Perform the per-event analysis
void
analyze
(
const
Event
&
event
)
{
const
double
weight
=
event
.
weight
();
Jets
jets
=
applyProjection
<
FastJets
>
(
event
,
"Jets"
).
jetsByEt
(
10.0
*
GeV
);
if
(
jets
.
size
()
<
2
)
{
vetoEvent
;
}
FourMomentum
jet1
=
jets
[
0
].
momentum
();
FourMomentum
jet2
=
jets
[
1
].
momentum
();
double
eta1
=
fabs
(
jet1
.
eta
());
double
eta2
=
fabs
(
jet2
.
eta
());
double
ET1
=
jet1
.
Et
();
double
ET2
=
jet2
.
Et
();
if
(
eta1
<
0.1
||
eta1
>
0.7
||
ET1
<
40.0
*
GeV
)
{
vetoEvent
;
}
if
(
eta2
<
0.1
||
eta2
>
3.0
)
{
vetoEvent
;
}
_h_ET
.
fill
(
eta2
,
ET1
,
weight
);
if
(
eta2
<
0.7
&&
ET2
>
40.0
*
GeV
)
_h_ET
.
fill
(
eta1
,
ET2
,
weight
);
}
/// Normalise histograms etc., after the run
void
finalize
()
{
double
deta1
=
1.2
;
_h_ET
.
scale
(
crossSection
()
/
nanobarn
/
sumOfWeights
()
/
deta1
/
2.0
,
this
);
}
//@}
private
:
/// @name Histograms
//@{
BinnedHistogram
<
double
>
_h_ET
;
//@}
};
// This global object acts as a hook for the plugin system
AnalysisBuilder
<
CDF_2001_S4517016
>
plugin_CDF_2001_S4517016
;
}
File Metadata
Details
Attached
Mime Type
text/x-c
Expires
Sat, Dec 21, 1:09 PM (1 d, 41 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4022897
Default Alt Text
CDF_2001_S4517016.cc (2 KB)
Attached To
rRIVETSVN rivetsvn
Event Timeline
Log In to Comment