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

Comparing UserCode/MitPhysics/Utils/src/MuonTools.cc (file contents):
Revision 1.9 by loizides, Tue Apr 7 15:37:10 2009 UTC vs.
Revision 1.23 by sixie, Sat Dec 31 23:18:39 2011 UTC

# Line 3 | Line 3
3   #include "MitPhysics/Utils/interface/MuonTools.h"
4   #include <TFile.h>
5  
6 + ClassImp(mithep::MuonTools)
7 +
8   using namespace mithep;
9  
10   //--------------------------------------------------------------------------------------------------
# Line 426 | Line 428 | TH2D *MuonTools::LoadHisto(const char *n
428    ret->SetDirectory(0);
429    return ret;
430   }
431 + //--------------------------------------------------------------------------------------------------
432 + Bool_t MuonTools::PassD0Cut(const Muon *mu, const VertexCol *vertices, Double_t fD0Cut, Int_t nVertex)
433 + {
434 +  Bool_t d0cut = kFALSE;
435 +  const Track *mt = mu->BestTrk();
436 +  if (!mt) return kFALSE;
437 +
438 +  Double_t d0_real = 1e30;
439 +  if(nVertex >= 0) d0_real = TMath::Abs(mt->D0Corrected(*vertices->At(nVertex)));
440 +  else            {
441 +    Double_t distVtx = 999.0;
442 +    Int_t closestVtx = 0;
443 +    for(UInt_t nv=0; nv<vertices->GetEntries(); nv++){
444 +      double dz = TMath::Abs(mt->DzCorrected(*vertices->At(nv)));
445 +      if(dz < distVtx) {
446 +        distVtx    = dz;
447 +        closestVtx = nv;
448 +      }
449 +    }
450 +    d0_real = TMath::Abs(mt->D0Corrected(*vertices->At(closestVtx)));
451 +  }
452 +  if(d0_real < fD0Cut) d0cut = kTRUE;
453 +  
454 +  return d0cut;
455 + }
456 +
457 + //--------------------------------------------------------------------------------------------------
458 + Bool_t MuonTools::PassD0Cut(const Muon *mu, const BeamSpotCol *beamspots, Double_t fD0Cut)
459 + {
460 +  Bool_t d0cut = kFALSE;
461 +  const Track *mt = mu->BestTrk();
462 +  if (!mt) return kFALSE;
463 +
464 +  // d0 cut
465 +  Double_t d0_real = 99999;
466 +  for(UInt_t i0 = 0; i0 < beamspots->GetEntries(); i0++) {
467 +    Double_t pD0 = mt->D0Corrected(*beamspots->At(i0));
468 +    if(TMath::Abs(pD0) < TMath::Abs(d0_real)) d0_real = TMath::Abs(pD0);
469 +  }
470 +  if(d0_real < fD0Cut) d0cut = kTRUE;
471 +  
472 +  return d0cut;
473 + }
474 +
475 + //--------------------------------------------------------------------------------------------------
476 + Bool_t MuonTools::PassDZCut(const Muon *mu, const VertexCol *vertices, Double_t fDZCut, Int_t nVertex)
477 + {
478 +  Bool_t dzcut = kFALSE;
479 +  const Track *mt = mu->BestTrk();
480 +  if (!mt) return kFALSE;
481 +
482 +  Double_t distVtx = 999.0;
483 +  if(nVertex >= 0) distVtx = TMath::Abs(mt->DzCorrected(*vertices->At(nVertex)));
484 +  else {
485 +    for(UInt_t nv=0; nv<vertices->GetEntries(); nv++){
486 +      double dz = TMath::Abs(mt->DzCorrected(*vertices->At(nv)));
487 +      if(dz < distVtx) {
488 +        distVtx = dz;
489 +      }
490 +    }
491 +  }
492 +
493 +  if(distVtx < fDZCut) dzcut = kTRUE;
494 +  
495 +  return dzcut;
496 + }
497 +
498 + //--------------------------------------------------------------------------------------------------
499 + Bool_t MuonTools::PassSoftMuonCut(const Muon *mu, const VertexCol *vertices, const Double_t fDZCut,
500 +                                  const Bool_t applyIso)
501 + {
502 +  if(mu->Pt() <= 3.0) return kFALSE;
503 +
504 +  if(!mu->IsTrackerMuon()) return kFALSE;
505 +  
506 +  if(!mu->Quality().Quality(MuonQuality::TMLastStationAngTight)) return kFALSE;
507 +
508 +  if(mu->BestTrk()->NHits() <= 10) return kFALSE;
509 +
510 +  if(!PassD0Cut(mu, vertices, 0.2, 0)) return kFALSE;
511 +
512 +  if(!PassDZCut(mu, vertices, fDZCut, 0)) return kFALSE;
513 +
514 +  if(applyIso == kTRUE){
515 +    Double_t totalIso = 1.0 * mu->IsoR03SumPt() +
516 +                        1.0 * mu->IsoR03EmEt() +
517 +                        1.0 * mu->IsoR03HadEt();
518 +    if (totalIso < (mu->Pt()*0.10) && mu->Pt() > 20.0) return kFALSE;
519 +  }
520 +
521 +  return kTRUE;
522 + }
523 +
524 + Double_t MuonTools::MuonEffectiveArea(EMuonEffectiveAreaType type, Double_t Eta) {
525 +
526 +  Double_t EffectiveArea = 0;
527 +  if (fabs(Eta) < 1.0) {
528 +    if (type == kMuChargedIso03) EffectiveArea = 0.000;
529 +    if (type == kMuNeutralIso03) EffectiveArea = 0.080;
530 +    if (type == kMuChargedIso04) EffectiveArea = 0.000;
531 +    if (type == kMuNeutralIso04) EffectiveArea = 0.163;
532 +    if (type == kMuHadEnergy)    EffectiveArea = 0.000;
533 +    if (type == kMuHoEnergy)     EffectiveArea = 0.000;
534 +    if (type == kMuEmEnergy)     EffectiveArea = 0.000;
535 +    if (type == kMuHadS9Energy)  EffectiveArea = 0.016;
536 +    if (type == kMuHoS9Energy)   EffectiveArea = 0.000;
537 +    if (type == kMuEmS9Energy)   EffectiveArea = 0.000;
538 +  } else if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) {
539 +    if (type == kMuChargedIso03) EffectiveArea = 0.000;
540 +    if (type == kMuNeutralIso03) EffectiveArea = 0.083;
541 +    if (type == kMuChargedIso04) EffectiveArea = 0.000;
542 +    if (type == kMuNeutralIso04) EffectiveArea = 0.168;
543 +    if (type == kMuHadEnergy)    EffectiveArea = 0.005;
544 +    if (type == kMuHoEnergy)     EffectiveArea = 0.000;
545 +    if (type == kMuEmEnergy)     EffectiveArea = 0.000;
546 +    if (type == kMuHadS9Energy)  EffectiveArea = 0.041;
547 +    if (type == kMuHoS9Energy)   EffectiveArea = 0.000;
548 +    if (type == kMuEmS9Energy)   EffectiveArea = 0.000;
549 +  } else if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) {
550 +    if (type == kMuChargedIso03) EffectiveArea = 0.000;
551 +    if (type == kMuNeutralIso03) EffectiveArea = 0.060;
552 +    if (type == kMuChargedIso04) EffectiveArea = 0.000;
553 +    if (type == kMuNeutralIso04) EffectiveArea = 0.131;
554 +    if (type == kMuHadEnergy)    EffectiveArea = 0.020;
555 +    if (type == kMuHoEnergy)     EffectiveArea = 0.000;
556 +    if (type == kMuEmEnergy)     EffectiveArea = 0.000;
557 +    if (type == kMuHadS9Energy)  EffectiveArea = 0.072;
558 +    if (type == kMuHoS9Energy)   EffectiveArea = 0.000;
559 +    if (type == kMuEmS9Energy)   EffectiveArea = 0.000;
560 +  } else if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.25 ) {
561 +    if (type == kMuChargedIso03) EffectiveArea = 0.000;
562 +    if (type == kMuNeutralIso03) EffectiveArea = 0.066;
563 +    if (type == kMuChargedIso04) EffectiveArea = 0.000;
564 +    if (type == kMuNeutralIso04) EffectiveArea = 0.149;
565 +    if (type == kMuHadEnergy)    EffectiveArea = 0.056;
566 +    if (type == kMuHoEnergy)     EffectiveArea = 0.000;
567 +    if (type == kMuEmEnergy)     EffectiveArea = 0.000;
568 +    if (type == kMuHadS9Energy)  EffectiveArea = 0.148;
569 +    if (type == kMuHoS9Energy)   EffectiveArea = 0.000;
570 +    if (type == kMuEmS9Energy)   EffectiveArea = 0.000;
571 +  } else if (fabs(Eta) >= 2.25 && fabs(Eta) < 2.4 ) {
572 +    if (type == kMuChargedIso03) EffectiveArea = 0.000;
573 +    if (type == kMuNeutralIso03) EffectiveArea = 0.098;
574 +    if (type == kMuChargedIso04) EffectiveArea = 0.000;
575 +    if (type == kMuNeutralIso04) EffectiveArea = 0.200;
576 +    if (type == kMuHadEnergy)    EffectiveArea = 0.093;
577 +    if (type == kMuHoEnergy)     EffectiveArea = 0.000;
578 +    if (type == kMuEmEnergy)     EffectiveArea = 0.000;
579 +    if (type == kMuHadS9Energy)  EffectiveArea = 0.260;
580 +    if (type == kMuHoS9Energy)   EffectiveArea = 0.000;
581 +    if (type == kMuEmS9Energy)   EffectiveArea = 0.000;
582 +  }
583 +  return EffectiveArea;
584 + }
585 +
586  

Diff Legend

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