ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/dasu/UltraFastSim/ZHAnalysis.cc
(Generate patch)

Comparing UserCode/dasu/UltraFastSim/ZHAnalysis.cc (file contents):
Revision 1.2 by dasu, Thu Feb 24 07:34:29 2011 UTC vs.
Revision 1.4 by brownson, Mon Feb 28 20:54:48 2011 UTC

# Line 17 | Line 17 | bool ZHAnalysis::Run(UltraFastSim * ufs)
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  
# Line 92 | Line 93 | bool ZHAnalysis::Run(UltraFastSim * ufs)
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();
# Line 125 | Line 137 | bool ZHAnalysis::Run(UltraFastSim * ufs)
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  

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines