Jpp
Main Page
Namespaces
Classes
Files
File List
File Members
All
Classes
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Pages
examples
JPhysics
JGeane.cc
Go to the documentation of this file.
1
#include <string>
2
#include <iostream>
3
4
#include "TROOT.h"
5
#include "TFile.h"
6
#include "TH1D.h"
7
#include "TH2D.h"
8
9
#include "
JPhysics/JGeane.hh
"
10
11
#include "
Jeep/JParser.hh
"
12
#include "
Jeep/JMessage.hh
"
13
14
15
/**
16
* \file
17
*
18
* Example program to histogram muon energy loss.
19
* \author mdejong
20
*/
21
int
main
(
int
argc,
char
* argv[])
22
{
23
using namespace
std;
24
using namespace
JPP;
25
26
string
outputFile
;
27
bool
debug
;
28
29
try
{
30
31
JParser<>
zap(
"Example program to histogram muon energy loss."
);
32
33
zap[
'o'
] =
make_field
(
outputFile
) =
"geane.root"
;
34
zap[
'd'
] =
make_field
(
debug
);
35
36
zap(argc, argv);
37
}
38
catch
(
const
exception &error) {
39
FATAL
(error.what() << endl);
40
}
41
42
43
44
TFile out(
outputFile
.c_str(),
"recreate"
);
45
46
TH1D
h1
(
"h1"
, NULL, 900, -1.0, +8.0);
47
TH2D h2(
"h2"
, NULL, 90, -1.0, +8.0, 100, 0.0, 5.0);
48
49
for
(
int
i = 1; i <= h1.GetNbinsX(); ++i) {
50
51
const
double
x = h1.GetBinCenter(i);
52
const
double
E
= pow(10.0, x);
53
54
h1.Fill(x,
gWater
(E));
55
}
56
57
for
(
int
ix = 1; ix <= h2.GetXaxis()->GetNbins(); ++ix) {
58
for
(
int
iy = 1; iy <= h2.GetYaxis()->GetNbins(); ++iy) {
59
60
const
double
x = h2.GetXaxis()->GetBinCenter(ix);
61
const
double
y = h2.GetYaxis()->GetBinCenter(iy);
62
63
const
double
E
= pow(10.0, x);
64
const
double
dx = pow(10.0, y);
65
66
h2.Fill(x, y,
gWater
(E, dx));
67
}
68
}
69
70
out.Write();
71
out.Close();
72
}
JPARSER::JParser
Utility class to parse command line options.
Definition:
JParser.hh:1493
JGeane.hh
Energy loss of muon.
h1
then for HISTOGRAM in h0 h1
Definition:
JMatrixNZ.sh:69
JPHYSICS::gWater
static const JGeaneWater gWater
Function object for energy loss of muon in sea water.
Definition:
JGeane.hh:328
outputFile
string outputFile
Definition:
JDAQTimesliceSelector.cc:37
make_field
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition:
JParser.hh:1954
debug
int debug
debug level
Definition:
JSirene.cc:61
JMessage.hh
General purpose messaging.
FATAL
#define FATAL(A)
Definition:
JMessage.hh:67
JParser.hh
Utility class to parse command line options.
E
then usage $script[input file[working directory[option]]] nWhere option can be E
Definition:
JMuonPostfit.sh:37
main
int main(int argc, char *argv[])
Definition:
Main.cpp:15
Generated by
1.8.5