ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MPIAnalyzer/src/MPIntuple.cc
(Generate patch)

Comparing UserCode/MPIAnalyzer/src/MPIntuple.cc (file contents):
Revision 1.18 by naodell, Wed Oct 27 14:21:20 2010 UTC vs.
Revision 1.22 by naodell, Fri Nov 5 13:47:05 2010 UTC

# Line 126 | Line 126 | class MPIntuple : public edm::EDAnalyzer
126    TH1F * h1_trackDxy;
127    TH1F * h1_allTrackDeltaZ_PV;
128    TH1F * h1_allTrackDeltaZ_LJ;
129 –  TH1F * h1_ptHat;
129    TH2F * h2_nAssTracksVsJetPt;
130    TProfile * p1_nVtcs;
131  
# Line 143 | Line 142 | MPIntuple::MPIntuple(const edm::Paramete
142    triggerHLT_ = iConfig.getUntrackedParameter<bool>("triggerHLT");
143    rootfilename = iConfig.getUntrackedParameter<string>("rootfilename");
144    hlTriggerResults_ = iConfig.getUntrackedParameter<string>("HLTriggerResults","TriggerResults");
145 <  hltName_ = iConfig.getUntrackedParameter<string>("hltName","RECO");
145 >  hltName_ = iConfig.getUntrackedParameter<string>("hltName","HLT");
146   }
147  
148   MPIntuple::~MPIntuple()
# Line 229 | Line 228 | void MPIntuple::analyze(const edm::Event
228      for(PFJetCollection::const_iterator jet_iter = PFJets->begin(); jet_iter!= PFJets->end(); ++jet_iter){
229        
230        reco::PFJet myJet = reco::PFJet(*jet_iter);
231 <      if(myJet.et() > 5){
231 >      if(myJet.pt() > 5){
232  
233          for(jet_it_1 = bTags1.begin(); jet_it_1 != bTags1.end(); jet_it_1++) {
234             if( (fabs(jet_it_1->first->eta() - myJet.eta()) < .005) && (deltaPhi(jet_it_1->first->phi(),myJet.phi()) < .005)) break;
# Line 286 | Line 285 | void MPIntuple::analyze(const edm::Event
285  
286            for(TrackRefVector::const_iterator iTrack = tracks.begin(); iTrack != tracks.end(); ++iTrack){
287  
288 <            const reco::Track& myTrack = **iTrack;
289 <
290 <            sumTrackZ += myTrack.vz();
291 <            sumTrackIP += myTrack.dxy(vertexBeamSpot.position());
288 >            const reco::Track    &myJetTrack = **iTrack;
289 >            //const reco::TrackRef &myTrackRef = *iTrack;
290 >
291 >            sumTrackZ += myJetTrack.vz();
292 >            sumTrackIP += myJetTrack.dxy(vertexBeamSpot.position());
293              ++nAssociatedTracks;
294 <            
294 >
295 >            for(VertexCollection::const_iterator vtx_iter = primaryVtcs->begin(); vtx_iter!= primaryVtcs->end(); ++vtx_iter){
296 >              
297 >              reco::Vertex myVtx = reco::Vertex(*vtx_iter);
298 >              if(!myVtx.isValid() || myVtx.isFake()) continue;
299 >              
300 >              for(Vertex::trackRef_iterator iTrackRef = myVtx.tracks_begin(); iTrackRef != myVtx.tracks_end(); ++iTrackRef){
301 >
302 >                const edm::RefToBase<reco::Track> &myTrackRef = *iTrackRef;
303 >
304 >                if(myTrackRef.isAvailable()){
305 >
306 >                  const reco::Track &myVertexTrack = *myTrackRef.get();
307 >                                
308 >                
309 >                }
310 >              }      
311 >            }              
312            }
313            
314            meanTrackZ = sumTrackZ/nAssociatedTracks;
# Line 300 | Line 317 | void MPIntuple::analyze(const edm::Event
317            
318            for(TrackRefVector::const_iterator iTrack = tracks.begin(); iTrack != tracks.end(); ++iTrack){
319              
320 <            const reco::Track& myTrack = **iTrack;
320 >            const reco::Track &myJetTrack = **iTrack;
321              
322 <            h1_trackZ->Fill(meanTrackZ - myTrack.vz());
323 <            h1_trackDxy->Fill(myTrack.dxy(vertexBeamSpot.position()));
322 >            if(nAssociatedTracks == 3) h1_trackZ->Fill(meanTrackZ - myJetTrack.vz());
323 >            h1_trackDxy->Fill(myJetTrack.dxy(vertexBeamSpot.position()));
324              
325            }    
326            
# Line 343 | Line 360 | void MPIntuple::analyze(const edm::Event
360        qScale = GenInfoHandle->qScale();
361        ptHat  = (GenInfoHandle->hasBinningValues() ? GenInfoHandle->binningValues()[0] : 0.0);
362  
346 –      h1_ptHat->Fill(ptHat);
347 –
363      }
364          
365      for(GenJetCollection::const_iterator jet_iter = GenJets->begin(); jet_iter!= GenJets->end(); ++jet_iter){
# Line 422 | Line 437 | void  MPIntuple::beginJob()
437    h1_trackDxy          = new TH1F("h1_trackDxy", "dxy for all associated tracks", 50, -0.4, 0.4);
438    h1_allTrackDeltaZ_PV = new TH1F("h1_allTrackDeltaZ_PV", "#Delta z between jet and primary vertex", 60, -0.3, 0.3);
439    h1_allTrackDeltaZ_LJ = new TH1F("h1_allTrackDeltaZ_LJ", "#Delta z between leading jet and other jet vertices", 60, -0.3, 0.3);
425 –  h1_ptHat             = new TH1F("h1_ptHat", "ptHat", 50, 10, 100);
440    h2_nAssTracksVsJetPt = new TH2F("h2_nAssTracksVsJetPt", "nTracks vs Jet p_{T}", 50, 0, 100, 20, -0.5, 19.5);
441    p1_nVtcs             = new TProfile("p1_nVtcs", "<N> per run", 5200.0, 132300.0, 137500.0, 0.0, 6.0);
442  
# Line 456 | Line 470 | void MPIntuple::endJob()
470    h1_trackDxy->Write();
471    h1_allTrackDeltaZ_PV->Write();
472    h1_allTrackDeltaZ_LJ->Write();
459 –  h1_ptHat->Write();
473    h2_nAssTracksVsJetPt->Write();
474    p1_nVtcs->Write();
475  

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines