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

Comparing UserCode/MitPhysics/Utils/src/MuonIDMVA.cc (file contents):
Revision 1.9 by sixie, Tue Apr 24 12:48:30 2012 UTC vs.
Revision 1.11 by vdutta, Thu May 3 16:52:55 2012 UTC

# Line 848 | Line 848 | Double_t MuonIDMVA::MVAValue(const Muon
848        //************************************************************
849        // Lepton Footprint Removal
850        //************************************************************            
851 <      for (UInt_t q=0; q < goodElectrons->GetEntries() ; ++q) {
852 <        //if pf candidate matches an electron passing ID cuts, then veto it
853 <        if(pf->GsfTrk() && goodElectrons->At(q)->GsfTrk() &&
854 <           pf->GsfTrk() == goodElectrons->At(q)->GsfTrk()) IsLeptonFootprint = kTRUE;
855 <        if(pf->TrackerTrk() && goodElectrons->At(q)->TrackerTrk() &&
856 <           pf->TrackerTrk() == goodElectrons->At(q)->TrackerTrk()) IsLeptonFootprint = kTRUE;
857 <        //if pf candidate lies in veto regions of electron passing ID cuts, then veto it
858 <        if(pf->BestTrk() && fabs(goodElectrons->At(q)->SCluster()->Eta()) >= 1.479
859 <           && MathUtils::DeltaR(goodElectrons->At(q)->Mom(), pf->Mom()) < 0.015) IsLeptonFootprint = kTRUE;
860 <        if(pf->PFType() == PFCandidate::eGamma && fabs(goodElectrons->At(q)->SCluster()->Eta()) >= 1.479 &&
861 <           MathUtils::DeltaR(goodElectrons->At(q)->Mom(), pf->Mom()) < 0.08) IsLeptonFootprint = kTRUE;
851 >      if(goodElectrons) {
852 >        for (UInt_t q=0; q < goodElectrons->GetEntries() ; ++q) {
853 >          //if pf candidate matches an electron passing ID cuts, then veto it
854 >          if(pf->GsfTrk() && goodElectrons->At(q)->GsfTrk() &&
855 >             pf->GsfTrk() == goodElectrons->At(q)->GsfTrk()) IsLeptonFootprint = kTRUE;
856 >          if(pf->TrackerTrk() && goodElectrons->At(q)->TrackerTrk() &&
857 >             pf->TrackerTrk() == goodElectrons->At(q)->TrackerTrk()) IsLeptonFootprint = kTRUE;
858 >          //if pf candidate lies in veto regions of electron passing ID cuts, then veto it
859 >          if(pf->BestTrk() && fabs(goodElectrons->At(q)->SCluster()->Eta()) >= 1.479
860 >             && MathUtils::DeltaR(goodElectrons->At(q)->Mom(), pf->Mom()) < 0.015) IsLeptonFootprint = kTRUE;
861 >          if(pf->PFType() == PFCandidate::eGamma && fabs(goodElectrons->At(q)->SCluster()->Eta()) >= 1.479 &&
862 >             MathUtils::DeltaR(goodElectrons->At(q)->Mom(), pf->Mom()) < 0.08) IsLeptonFootprint = kTRUE;
863 >        }
864        }
865 <      for (UInt_t q=0; q < goodMuons->GetEntries() ; ++q) {
866 <        //if pf candidate matches an muon passing ID cuts, then veto it
867 <        if(pf->TrackerTrk() && goodMuons->At(q)->TrackerTrk() &&
868 <           pf->TrackerTrk() == goodMuons->At(q)->TrackerTrk()) IsLeptonFootprint = kTRUE;
869 <        //if pf candidate lies in veto regions of muon passing ID cuts, then veto it
870 <        if(pf->BestTrk() && MathUtils::DeltaR(goodMuons->At(q)->Mom(), pf->Mom()) < 0.01) IsLeptonFootprint = kTRUE;
865 >      if(goodMuons) {
866 >        for (UInt_t q=0; q < goodMuons->GetEntries() ; ++q) {
867 >          //if pf candidate matches an muon passing ID cuts, then veto it
868 >          if(pf->TrackerTrk() && goodMuons->At(q)->TrackerTrk() &&
869 >             pf->TrackerTrk() == goodMuons->At(q)->TrackerTrk()) IsLeptonFootprint = kTRUE;
870 >          //if pf candidate lies in veto regions of muon passing ID cuts, then veto it
871 >          if(pf->BestTrk() && MathUtils::DeltaR(goodMuons->At(q)->Mom(), pf->Mom()) < 0.01) IsLeptonFootprint = kTRUE;
872 >        }
873        }
874  
875        if (!IsLeptonFootprint) {
# Line 909 | Line 913 | Double_t MuonIDMVA::MVAValue(const Muon
913      } //in 1.0 dr cone
914    } //loop over PF candidates
915    
916 <  Double_t fMVAVar_ChargedIso_DR0p0To0p1  = 0;
917 <  Double_t fMVAVar_ChargedIso_DR0p1To0p2  = 0;
918 <  Double_t fMVAVar_ChargedIso_DR0p2To0p3  = 0;
919 <  Double_t fMVAVar_ChargedIso_DR0p3To0p4  = 0;
920 <  Double_t fMVAVar_ChargedIso_DR0p4To0p5  = 0;
921 <  Double_t fMVAVar_GammaIso_DR0p0To0p1  = 0;
922 <  Double_t fMVAVar_GammaIso_DR0p1To0p2  = 0;
923 <  Double_t fMVAVar_GammaIso_DR0p2To0p3  = 0;
924 <  Double_t fMVAVar_GammaIso_DR0p3To0p4  = 0;
925 <  Double_t fMVAVar_GammaIso_DR0p4To0p5  = 0;
926 <  Double_t fMVAVar_NeutralHadronIso_DR0p0To0p1  = 0;
927 <  Double_t fMVAVar_NeutralHadronIso_DR0p1To0p2  = 0;
928 <  Double_t fMVAVar_NeutralHadronIso_DR0p2To0p3  = 0;
929 <  Double_t fMVAVar_NeutralHadronIso_DR0p3To0p4  = 0;
930 <  Double_t fMVAVar_NeutralHadronIso_DR0p4To0p5  = 0;
916 > //   Double_t fMVAVar_ChargedIso_DR0p0To0p1  = 0;
917 > //   Double_t fMVAVar_ChargedIso_DR0p1To0p2  = 0;
918 > //   Double_t fMVAVar_ChargedIso_DR0p2To0p3  = 0;
919 > //   Double_t fMVAVar_ChargedIso_DR0p3To0p4  = 0;
920 > //   Double_t fMVAVar_ChargedIso_DR0p4To0p5  = 0;
921 > //   Double_t fMVAVar_GammaIso_DR0p0To0p1  = 0;
922 > //   Double_t fMVAVar_GammaIso_DR0p1To0p2  = 0;
923 > //   Double_t fMVAVar_GammaIso_DR0p2To0p3  = 0;
924 > //   Double_t fMVAVar_GammaIso_DR0p3To0p4  = 0;
925 > //   Double_t fMVAVar_GammaIso_DR0p4To0p5  = 0;
926 > //   Double_t fMVAVar_NeutralHadronIso_DR0p0To0p1  = 0;
927 > //   Double_t fMVAVar_NeutralHadronIso_DR0p1To0p2  = 0;
928 > //   Double_t fMVAVar_NeutralHadronIso_DR0p2To0p3  = 0;
929 > //   Double_t fMVAVar_NeutralHadronIso_DR0p3To0p4  = 0;
930 > //   Double_t fMVAVar_NeutralHadronIso_DR0p4To0p5  = 0;
931  
932    fMVAVar_ChargedIso_DR0p0To0p1   = TMath::Min((tmpChargedIso_DR0p0To0p1)/mu->Pt(), 2.5);
933    fMVAVar_ChargedIso_DR0p1To0p2   = TMath::Min((tmpChargedIso_DR0p1To0p2)/mu->Pt(), 2.5);

Diff Legend

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