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

Comparing UserCode/Vuko/WZAnalysis/src/MuonAnalyzer.cc (file contents):
Revision 1.4 by senka, Fri Apr 4 09:08:11 2008 UTC vs.
Revision 1.6 by senka, Fri Apr 4 10:07:27 2008 UTC

# Line 300 | Line 300 | MuonAnalyzer::analyze(const edm::Event&
300      muon_mother_id=0;
301      W_mother_id=0;
302  
303 +    muon_global=0;
304  
305      muon_global_chi2=-10.;
306      muon_global_p=-10.;
# Line 317 | Line 318 | MuonAnalyzer::analyze(const edm::Event&
318      muon_global_outerPosition_x=-10.;
319      muon_global_outerPosition_y=-10.;
320      muon_global_outerPosition_z=-10.;
321 +
322 +    muon_STA=0;
323      
324      muon_STA_chi2=-10.;
325      muon_STA_p=-10.;
# Line 335 | Line 338 | MuonAnalyzer::analyze(const edm::Event&
338      muon_STA_outerPosition_y=-10.;
339      muon_STA_outerPosition_z=-10.;
340      
341 +    muon_track=0;
342 +
343      muon_track_chi2=-10.;
344      muon_track_p=-10.;
345      muon_track_pt=-10.;
# Line 373 | Line 378 | MuonAnalyzer::analyze(const edm::Event&
378        muon_Pt = matchedmuon->pt();
379        
380        const reco::Muon* matchedmuon_recoMuon = dynamic_cast<const reco::Muon *> (matchedmuon); // conversion of matchedmuon from reco::Candidate to reco::Muon
381 + //       reco::TrackRef combined_muon = new reco::Track::Track();
382 + //       combined_muon = matchedmuon_recoMuon->combinedMuon();
383        
384 <      muon_global_chi2 = matchedmuon_recoMuon->combinedMuon()->chi2();
385 <      muon_global_p = matchedmuon_recoMuon->combinedMuon()->p();
386 <      muon_global_pt = matchedmuon_recoMuon->combinedMuon()->pt();
387 <      muon_global_outerP = matchedmuon_recoMuon->combinedMuon()->outerP();
388 <      muon_global_outerPt = matchedmuon_recoMuon->combinedMuon()->outerPt();
389 <      muon_global_ndof = matchedmuon_recoMuon->combinedMuon()->ndof();
390 <      muon_global_normalizedChi2 = matchedmuon_recoMuon->combinedMuon()->normalizedChi2();
391 <      muon_global_recHitsSize = matchedmuon_recoMuon->combinedMuon()->recHitsSize();
392 <      muon_global_numberOfLostHits = matchedmuon_recoMuon->combinedMuon()->numberOfLostHits();
393 <      muon_global_numberOfValidHits = matchedmuon_recoMuon->combinedMuon()->numberOfValidHits();
394 <      muon_global_innerPosition_x = matchedmuon_recoMuon->combinedMuon()->innerPosition().x();
395 <      muon_global_innerPosition_y = matchedmuon_recoMuon->combinedMuon()->innerPosition().y();
396 <      muon_global_innerPosition_z = matchedmuon_recoMuon->combinedMuon()->innerPosition().z();
397 <      muon_global_outerPosition_x = matchedmuon_recoMuon->combinedMuon()->outerPosition().x();
398 <      muon_global_outerPosition_y = matchedmuon_recoMuon->combinedMuon()->outerPosition().y();
399 <      muon_global_outerPosition_z = matchedmuon_recoMuon->combinedMuon()->outerPosition().z();
400 <      
401 <      muon_STA_chi2 = matchedmuon_recoMuon->standAloneMuon()->chi2();
402 <      muon_STA_p = matchedmuon_recoMuon->standAloneMuon()->p();
403 <      muon_STA_pt = matchedmuon_recoMuon->standAloneMuon()->pt();
404 <      muon_STA_outerP = matchedmuon_recoMuon->standAloneMuon()->outerP();
405 <      muon_STA_outerPt = matchedmuon_recoMuon->standAloneMuon()->outerPt();
406 <      muon_STA_ndof = matchedmuon_recoMuon->standAloneMuon()->ndof();
407 <      muon_STA_normalizedChi2 = matchedmuon_recoMuon->standAloneMuon()->normalizedChi2();
408 <      muon_STA_recHitsSize = matchedmuon_recoMuon->standAloneMuon()->recHitsSize();
409 <      muon_STA_numberOfLostHits = matchedmuon_recoMuon->standAloneMuon()->numberOfLostHits();
410 <      muon_STA_numberOfValidHits = matchedmuon_recoMuon->standAloneMuon()->numberOfValidHits();
411 <      muon_STA_innerPosition_x = matchedmuon_recoMuon->standAloneMuon()->innerPosition().x();
412 <      muon_STA_innerPosition_y = matchedmuon_recoMuon->standAloneMuon()->innerPosition().y();
413 <      muon_STA_innerPosition_z = matchedmuon_recoMuon->standAloneMuon()->innerPosition().z();
414 <      muon_STA_outerPosition_x = matchedmuon_recoMuon->standAloneMuon()->outerPosition().x();
415 <      muon_STA_outerPosition_y = matchedmuon_recoMuon->standAloneMuon()->outerPosition().y();
416 <      muon_STA_outerPosition_z = matchedmuon_recoMuon->standAloneMuon()->outerPosition().z();
417 <      
418 <      muon_track_chi2 = matchedmuon_recoMuon->track()->chi2();
419 <      muon_track_p = matchedmuon_recoMuon->track()->p();
420 <      muon_track_pt = matchedmuon_recoMuon->track()->pt();
421 <      muon_track_outerP = matchedmuon_recoMuon->track()->outerP();
422 <      muon_track_outerPt = matchedmuon_recoMuon->track()->outerPt();
423 <      muon_track_ndof = matchedmuon_recoMuon->track()->ndof();
424 <      muon_track_normalizedChi2 = matchedmuon_recoMuon->track()->normalizedChi2();
425 <      muon_track_recHitsSize = matchedmuon_recoMuon->track()->recHitsSize();
426 <      muon_track_numberOfLostHits = matchedmuon_recoMuon->track()->numberOfLostHits();
427 <      muon_track_numberOfValidHits = matchedmuon_recoMuon->track()->numberOfValidHits();
428 <      muon_track_innerPosition_x = matchedmuon_recoMuon->track()->innerPosition().x();
429 <      muon_track_innerPosition_y = matchedmuon_recoMuon->track()->innerPosition().y();
430 <      muon_track_innerPosition_z = matchedmuon_recoMuon->track()->innerPosition().z();
431 <      muon_track_outerPosition_x = matchedmuon_recoMuon->track()->outerPosition().x();
432 <      muon_track_outerPosition_y = matchedmuon_recoMuon->track()->outerPosition().y();
433 <      muon_track_outerPosition_z = matchedmuon_recoMuon->track()->outerPosition().z();
434 <      
384 >      if (&(matchedmuon_recoMuon->combinedMuon())!=0){
385 >      //      if (combined_muon==0){
386 >        muon_global=1;
387 >        muon_global_chi2 = matchedmuon_recoMuon->combinedMuon()->chi2();
388 >        muon_global_p = matchedmuon_recoMuon->combinedMuon()->p();
389 >        muon_global_pt = matchedmuon_recoMuon->combinedMuon()->pt();
390 >        muon_global_outerP = matchedmuon_recoMuon->combinedMuon()->outerP();
391 >        muon_global_outerPt = matchedmuon_recoMuon->combinedMuon()->outerPt();
392 >        muon_global_ndof = matchedmuon_recoMuon->combinedMuon()->ndof();
393 >        muon_global_normalizedChi2 = matchedmuon_recoMuon->combinedMuon()->normalizedChi2();
394 >        muon_global_recHitsSize = matchedmuon_recoMuon->combinedMuon()->recHitsSize();
395 >        muon_global_numberOfLostHits = matchedmuon_recoMuon->combinedMuon()->numberOfLostHits();
396 >        muon_global_numberOfValidHits = matchedmuon_recoMuon->combinedMuon()->numberOfValidHits();
397 >        muon_global_innerPosition_x = matchedmuon_recoMuon->combinedMuon()->innerPosition().x();
398 >        muon_global_innerPosition_y = matchedmuon_recoMuon->combinedMuon()->innerPosition().y();
399 >        muon_global_innerPosition_z = matchedmuon_recoMuon->combinedMuon()->innerPosition().z();
400 >        muon_global_outerPosition_x = matchedmuon_recoMuon->combinedMuon()->outerPosition().x();
401 >        muon_global_outerPosition_y = matchedmuon_recoMuon->combinedMuon()->outerPosition().y();
402 >        muon_global_outerPosition_z = matchedmuon_recoMuon->combinedMuon()->outerPosition().z();
403 >      }
404 >
405 >      if (&(matchedmuon_recoMuon->standAloneMuon())!=0){
406 >        muon_STA=1;
407 >        muon_STA_chi2 = matchedmuon_recoMuon->standAloneMuon()->chi2();
408 >        muon_STA_p = matchedmuon_recoMuon->standAloneMuon()->p();
409 >        muon_STA_pt = matchedmuon_recoMuon->standAloneMuon()->pt();
410 >        muon_STA_outerP = matchedmuon_recoMuon->standAloneMuon()->outerP();
411 >        muon_STA_outerPt = matchedmuon_recoMuon->standAloneMuon()->outerPt();
412 >        muon_STA_ndof = matchedmuon_recoMuon->standAloneMuon()->ndof();
413 >        muon_STA_normalizedChi2 = matchedmuon_recoMuon->standAloneMuon()->normalizedChi2();
414 >        muon_STA_recHitsSize = matchedmuon_recoMuon->standAloneMuon()->recHitsSize();
415 >        muon_STA_numberOfLostHits = matchedmuon_recoMuon->standAloneMuon()->numberOfLostHits();
416 >        muon_STA_numberOfValidHits = matchedmuon_recoMuon->standAloneMuon()->numberOfValidHits();
417 >        muon_STA_innerPosition_x = matchedmuon_recoMuon->standAloneMuon()->innerPosition().x();
418 >        muon_STA_innerPosition_y = matchedmuon_recoMuon->standAloneMuon()->innerPosition().y();
419 >        muon_STA_innerPosition_z = matchedmuon_recoMuon->standAloneMuon()->innerPosition().z();
420 >        muon_STA_outerPosition_x = matchedmuon_recoMuon->standAloneMuon()->outerPosition().x();
421 >        muon_STA_outerPosition_y = matchedmuon_recoMuon->standAloneMuon()->outerPosition().y();
422 >        muon_STA_outerPosition_z = matchedmuon_recoMuon->standAloneMuon()->outerPosition().z();
423 >      }
424 >
425 >      if (&(matchedmuon_recoMuon->track())!=0){
426 >        muon_track=1;
427 >        muon_track_chi2 = matchedmuon_recoMuon->track()->chi2();
428 >        muon_track_p = matchedmuon_recoMuon->track()->p();
429 >        muon_track_pt = matchedmuon_recoMuon->track()->pt();
430 >        muon_track_outerP = matchedmuon_recoMuon->track()->outerP();
431 >        muon_track_outerPt = matchedmuon_recoMuon->track()->outerPt();
432 >        muon_track_ndof = matchedmuon_recoMuon->track()->ndof();
433 >        muon_track_normalizedChi2 = matchedmuon_recoMuon->track()->normalizedChi2();
434 >        muon_track_recHitsSize = matchedmuon_recoMuon->track()->recHitsSize();
435 >        muon_track_numberOfLostHits = matchedmuon_recoMuon->track()->numberOfLostHits();
436 >        muon_track_numberOfValidHits = matchedmuon_recoMuon->track()->numberOfValidHits();
437 >        muon_track_innerPosition_x = matchedmuon_recoMuon->track()->innerPosition().x();
438 >        muon_track_innerPosition_y = matchedmuon_recoMuon->track()->innerPosition().y();
439 >        muon_track_innerPosition_z = matchedmuon_recoMuon->track()->innerPosition().z();
440 >        muon_track_outerPosition_x = matchedmuon_recoMuon->track()->outerPosition().x();
441 >        muon_track_outerPosition_y = matchedmuon_recoMuon->track()->outerPosition().y();
442 >        muon_track_outerPosition_z = matchedmuon_recoMuon->track()->outerPosition().z();
443 >      }
444        //        if (fabs(MatchedMuonWithJet(&(*jetCands)[i])->eta())<2.5) {
445        //          cout <<"   muon prolazi eta cut" <<endl;
446        
# Line 641 | Line 657 | void MuonAnalyzer::initialiseTree() {
657  
658    muonTree->Branch("W_mother_id",&W_mother_id,"W_mother_id/I");
659  
660 +  muonTree->Branch("muon_global", &muon_global , "muon_global/I");
661    muonTree->Branch("muon_global_chi2", &muon_global_chi2 , "muon_global_chi2/F");
662    muonTree->Branch("muon_global_p",&muon_global_p," muon_global_p/F");
663    muonTree->Branch("muon_global_pt",&muon_global_pt,  "muon_global_pt/F");
# Line 658 | Line 675 | void MuonAnalyzer::initialiseTree() {
675    muonTree->Branch("muon_global_outerPosition_y",& muon_global_outerPosition_y,"muon_global_outerPosition_y/F");
676    muonTree->Branch("muon_global_outerPosition_z",& muon_global_outerPosition_z,"muon_global_outerPosition_z/F");
677  
678 +  muonTree->Branch("muon_STA", &muon_STA , "muon_STA/I");
679    muonTree->Branch("muon_STA_chi2", &muon_STA_chi2 , "muon_STA_chi2/F");
680    muonTree->Branch("muon_STA_p",&muon_STA_p," muon_STA_p/F");
681    muonTree->Branch("muon_STA_pt",&muon_STA_pt,  "muon_STA_pt/F");
# Line 675 | Line 693 | void MuonAnalyzer::initialiseTree() {
693    muonTree->Branch("muon_STA_outerPosition_y",& muon_STA_outerPosition_y,"muon_STA_outerPosition_y/F");
694    muonTree->Branch("muon_STA_outerPosition_z",& muon_STA_outerPosition_z,"muon_STA_outerPosition_z/F");
695  
696 +  muonTree->Branch("muon_track", &muon_track , "muon_track/I");
697    muonTree->Branch("muon_track_chi2", &muon_track_chi2 , "muon_track_chi2/F");
698    muonTree->Branch("muon_track_p",&muon_track_p," muon_track_p/F");
699    muonTree->Branch("muon_track_pt",&muon_track_pt,  "muon_track_pt/F");

Diff Legend

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