-
Notifications
You must be signed in to change notification settings - Fork 1
Expand file tree
/
Copy pathRtuple.cc
More file actions
executable file
·133 lines (111 loc) · 5.47 KB
/
Copy pathRtuple.cc
File metadata and controls
executable file
·133 lines (111 loc) · 5.47 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
#include "Rtuple.h"
void setRtupleAddresses(TTree* tr, Rtuple& vr ) {
tr->SetBranchAddress("runId", &vr.runId);
tr->SetBranchAddress("eventId", &vr.eventId);
tr->SetBranchAddress("lumiSection", &vr.lumi);
tr->SetBranchAddress("triggered", &vr.triggered);
tr->SetBranchAddress("L1a", &vr.L1a);
tr->SetBranchAddress("nPhotons", &vr.nPhotons);
tr->SetBranchAddress("nJets", &vr.nJets);
tr->SetBranchAddress("nMuons", &vr.nMuons);
tr->SetBranchAddress("nElectrons", &vr.nElectrons);
tr->SetBranchAddress("nVertices", &vr.nVertices );
tr->SetBranchAddress("totalNVtx", &vr.totalNVtx);
tr->SetBranchAddress("t_metPx", &vr.t_metPx );
tr->SetBranchAddress("t_metPy", &vr.t_metPy );
tr->SetBranchAddress("t_met", &vr.t_met );
tr->SetBranchAddress("metPx", &vr.metPx );
tr->SetBranchAddress("metPy", &vr.metPy );
tr->SetBranchAddress("met", &vr.metE );
tr->SetBranchAddress("met_dx1", &vr.met_dx1 );
tr->SetBranchAddress("met_dy1", &vr.met_dy1 );
tr->SetBranchAddress("met_dx2", &vr.met_dx2 );
tr->SetBranchAddress("met_dy2", &vr.met_dy2 );
tr->SetBranchAddress("met_dx3", &vr.met_dx3 );
tr->SetBranchAddress("met_dy3", &vr.met_dy3 );
tr->SetBranchAddress("jetPx", &vr.jetPx );
tr->SetBranchAddress("jetPy", &vr.jetPy );
tr->SetBranchAddress("jetPz", &vr.jetPz );
tr->SetBranchAddress("jetE", &vr.jetE );
tr->SetBranchAddress("jetNDau", &vr.jetNDau );
tr->SetBranchAddress("jetCM", &vr.jetCM );
tr->SetBranchAddress("jetCEF", &vr.jetCEF );
tr->SetBranchAddress("jetCHF", &vr.jetCHF );
tr->SetBranchAddress("jetNHF", &vr.jetNHF );
tr->SetBranchAddress("jetNEF", &vr.jetNEF );
tr->SetBranchAddress("jecUnc", &vr.jecUnc );
tr->SetBranchAddress("jerUnc", &vr.jerUnc );
tr->SetBranchAddress("jetTime", &vr.jetTime );
tr->SetBranchAddress("jetTimeErr", &vr.jetTimeErr );
tr->SetBranchAddress("phoPx", &vr.phoPx );
tr->SetBranchAddress("phoPy", &vr.phoPy );
tr->SetBranchAddress("phoPz", &vr.phoPz );
tr->SetBranchAddress("phoE", &vr.phoE );
tr->SetBranchAddress("phoEcalIso", &vr.phoEcalIso );
tr->SetBranchAddress("phoHcalIso", &vr.phoHcalIso );
tr->SetBranchAddress("phoTrkIso", &vr.phoTrkIso );
tr->SetBranchAddress("cHadIso", &vr.cHadIso );
tr->SetBranchAddress("nHadIso", &vr.nHadIso );
tr->SetBranchAddress("photIso", &vr.photIso );
tr->SetBranchAddress("sigmaIeta", &vr.sigmaIeta );
tr->SetBranchAddress("dR_TrkPho", &vr.dR_TrkPho );
tr->SetBranchAddress("phoHoverE", &vr.phoHovE );
tr->SetBranchAddress("sMinPho", &vr.sMinPho );
tr->SetBranchAddress("sMajPho", &vr.sMajPho );
tr->SetBranchAddress("seedTime", &vr.seedTime );
tr->SetBranchAddress("aveTime", &vr.aveTime );
tr->SetBranchAddress("cscdPhi", &vr.cscdPhi );
tr->SetBranchAddress("cscRho", &vr.cscRho );
tr->SetBranchAddress("cscTime", &vr.cscTime );
tr->SetBranchAddress("dtdPhi", &vr.dtdPhi );
tr->SetBranchAddress("dtdEta", &vr.dtdEta );
tr->SetBranchAddress("nXtals", &vr.nXtals );
tr->SetBranchAddress("timeChi2", &vr.timeChi2 );
tr->SetBranchAddress("seedSwissX", &vr.seedSwissX );
tr->SetBranchAddress("seedE", &vr.seedE );
tr->SetBranchAddress("elePx", &vr.elePx );
tr->SetBranchAddress("elePy", &vr.elePy );
tr->SetBranchAddress("elePz", &vr.elePz );
tr->SetBranchAddress("eleE", &vr.eleE );
tr->SetBranchAddress("eleEcalIso", &vr.eleEcalIso );
tr->SetBranchAddress("eleHcalIso", &vr.eleHcalIso );
tr->SetBranchAddress("eleTrkIso", &vr.eleTrkIso );
tr->SetBranchAddress("eleNLostHits", &vr.eleNLostHits );
tr->SetBranchAddress("muPx", &vr.muPx );
tr->SetBranchAddress("muPy", &vr.muPy );
tr->SetBranchAddress("muPz", &vr.muPz );
tr->SetBranchAddress("muE", &vr.muE );
//tr->SetBranchAddress("muEcalIso", muEcalIso );
//tr->SetBranchAddress("muHcalIso", muHcalIso );
//tr->SetBranchAddress("muTrkIso", muTrkIso );
//tr->SetBranchAddress("vtxX", vtxX );
//tr->SetBranchAddress("vtxY", vtxY );
tr->SetBranchAddress("vtxZ", &vr.vtxZ );
tr->SetBranchAddress("vtxChi2", &vr.vtxChi2 );
tr->SetBranchAddress("vtxNdof", &vr.vtxNdof );
tr->SetBranchAddress("nGen", &vr.nGen);
tr->SetBranchAddress("pdgId", &vr.pdgId );
tr->SetBranchAddress("momId", &vr.momId );
}
/*
void setRtupleBranches(TTree* chain, Rtuple& treeVars) {
chain -> Branch("nJets", &treeVars.nJets, "nJets/I");
chain -> Branch("nPhotons", &treeVars.nPhotons, "nPhotons/I");
chain -> Branch("met", &treeVars.met, "met/F");
chain -> Branch("g1Pt", &treeVars.g1Pt, "g1Pt/F");
chain -> Branch("g2Pt", &treeVars.g2Pt, "g2Pt/F");
chain -> Branch("g1Time", &treeVars.g1Time, "g1Time/F");
chain -> Branch("g2Time", &treeVars.g2Time, "g2Time/F");
chain -> Branch("bgTime", &treeVars.bgTime, "bgTime/F");
}
void initRtuple(TTree* chain, Rtuple& treeVars) {
treeVars.nJets = 0 ;
treeVars.nPhotons = 0 ;
treeVars.met = 0 ;
treeVars.g1Pt = 0 ;
treeVars.g2Pt = 0 ;
treeVars.g1Time = -25 ;
treeVars.g2Time = -25 ;
treeVars.bgTime = -25 ;
}
*/