17 |
|
vector<TParticle> Photons_ = ufs->photonList(); |
18 |
|
vector<TParticle> Electrons_ = ufs->electronList(); |
19 |
|
vector<TParticle> Muons_ = ufs->muonList(); |
20 |
+ |
vector<TParticle> MuonsStdGeom_ = ufs->muonListStdGeom(); |
21 |
|
vector<TParticle> Taus_ = ufs->tauList(); |
22 |
|
vector<TParticle> ChargedHadrons_ = ufs->chargedHadronList(); |
23 |
|
vector<TParticle> NeutralHadrons_ = ufs->neutralHadronList(); |
24 |
< |
vector<fastjet::PseudoJet> Jets_ = ufs->jetList(); |
25 |
< |
vector<fastjet::PseudoJet> BJets_ = ufs->bJetList(); |
26 |
< |
vector<fastjet::PseudoJet> BJetsStdGeom_ = ufs->bJetListStdGeom(); |
24 |
> |
vector<TLorentzVector> Jets_ = ufs->jetList(); |
25 |
> |
vector<TLorentzVector> BJets_ = ufs->bJetList(); |
26 |
> |
vector<TLorentzVector> BJetsStdGeom_ = ufs->bJetListStdGeom(); |
27 |
|
|
28 |
|
myevent_->clear(); |
29 |
|
|
93 |
|
obj.phi=it->Phi(); |
94 |
|
(myevent_->Muons).push_back(obj); |
95 |
|
} |
96 |
+ |
for(vector<TParticle>::const_iterator it=MuonsStdGeom_.begin();it!=MuonsStdGeom_.end();it++){ |
97 |
+ |
myobject obj; |
98 |
+ |
obj.px=it->Px(); |
99 |
+ |
obj.py=it->Py(); |
100 |
+ |
obj.pz=it->Pz(); |
101 |
+ |
obj.E=it->Energy(); |
102 |
+ |
obj.pt=it->Pt(); |
103 |
+ |
obj.eta=it->Eta(); |
104 |
+ |
obj.phi=it->Phi(); |
105 |
+ |
(myevent_->MuonsStdGeom).push_back(obj); |
106 |
+ |
} |
107 |
|
for(vector<TParticle>::const_iterator it=Taus_.begin();it!=Taus_.end();it++){ |
108 |
|
myobject obj; |
109 |
|
obj.px=it->Px(); |
137 |
|
obj.phi=it->Phi(); |
138 |
|
(myevent_->NeutralHadrons).push_back(obj); |
139 |
|
} |
140 |
< |
for(vector<fastjet::PseudoJet>::const_iterator it=Jets_.begin();it!=Jets_.end();it++){ |
140 |
> |
for(vector<TLorentzVector>::const_iterator it=Jets_.begin();it!=Jets_.end();it++){ |
141 |
|
myobject obj; |
142 |
< |
obj.px=it->px(); |
143 |
< |
obj.py=it->py(); |
144 |
< |
obj.pz=it->pz(); |
142 |
> |
obj.px=it->Px(); |
143 |
> |
obj.py=it->Py(); |
144 |
> |
obj.pz=it->Pz(); |
145 |
|
obj.E=it->E(); |
146 |
< |
obj.pt=it->perp(); |
147 |
< |
obj.eta=it->eta(); |
148 |
< |
obj.phi=it->phi(); |
146 |
> |
obj.pt=it->Pt(); |
147 |
> |
obj.eta=it->Eta(); |
148 |
> |
obj.phi=it->Phi(); |
149 |
|
(myevent_->Jets).push_back(obj); |
150 |
|
} |
151 |
< |
for(vector<fastjet::PseudoJet>::const_iterator it=BJets_.begin();it!=BJets_.end();it++){ |
151 |
> |
for(vector<TLorentzVector>::const_iterator it=BJets_.begin();it!=BJets_.end();it++){ |
152 |
|
myobject obj; |
153 |
< |
obj.px=it->px(); |
154 |
< |
obj.py=it->py(); |
155 |
< |
obj.pz=it->pz(); |
153 |
> |
obj.px=it->Px(); |
154 |
> |
obj.py=it->Py(); |
155 |
> |
obj.pz=it->Pz(); |
156 |
|
obj.E=it->E(); |
157 |
< |
obj.pt=it->perp(); |
158 |
< |
obj.eta=it->eta(); |
159 |
< |
obj.phi=it->phi(); |
157 |
> |
obj.pt=it->Pt(); |
158 |
> |
obj.eta=it->Eta(); |
159 |
> |
obj.phi=it->Phi(); |
160 |
|
(myevent_->BJets).push_back(obj); |
161 |
|
} |
162 |
< |
for(vector<fastjet::PseudoJet>::const_iterator it=BJetsStdGeom_.begin();it!=BJetsStdGeom_.end();it++){ |
162 |
> |
for(vector<TLorentzVector>::const_iterator it=BJetsStdGeom_.begin();it!=BJetsStdGeom_.end();it++){ |
163 |
|
myobject obj; |
164 |
< |
obj.px=it->px(); |
165 |
< |
obj.py=it->py(); |
166 |
< |
obj.pz=it->pz(); |
164 |
> |
obj.px=it->Px(); |
165 |
> |
obj.py=it->Py(); |
166 |
> |
obj.pz=it->Pz(); |
167 |
|
obj.E=it->E(); |
168 |
< |
obj.pt=it->perp(); |
169 |
< |
obj.eta=it->eta(); |
170 |
< |
obj.phi=it->phi(); |
168 |
> |
obj.pt=it->Pt(); |
169 |
> |
obj.eta=it->Eta(); |
170 |
> |
obj.phi=it->Phi(); |
171 |
|
(myevent_->BJetsStdGeom).push_back(obj); |
172 |
|
} |
173 |
|
|