ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/VHbbAnalysis/HbbAnalyzer/plugins/HbbAnalyzerNew.cc
(Generate patch)

Comparing UserCode/VHbbAnalysis/HbbAnalyzer/plugins/HbbAnalyzerNew.cc (file contents):
Revision 1.37 by degrutto, Thu Sep 15 14:00:01 2011 UTC vs.
Revision 1.40 by tboccali, Mon Sep 19 10:14:16 2011 UTC

# Line 441 | Line 441 | BTagSFContainer btagSFs;
441        TLorentzVector gJp4;
442        if(gJ){
443          gJp4.SetPtEtaPhiE(gJ->pt(),gJ->eta(),gJ->phi(),gJ->energy());
444 <        sj. bestMCp4mom = gJp4;
444 >        sj.bestMCp4mom = gJp4;
445          if(verbose_){
446            std::clog << "genJet matched Pt = " << gJp4.Pt() << std::endl;
447            std::clog << "genJet matched eta = " << gJp4.Eta() << std::endl;
# Line 481 | Line 481 | BTagSFContainer btagSFs;
481        TLorentzVector gJp4;
482        if(gJ){
483          gJp4.SetPtEtaPhiE(gJ->pt(),gJ->eta(),gJ->phi(),gJ->energy());
484 <        sj. bestMCp4mom = gJp4;
484 >        sj.bestMCp4mom = gJp4;
485          if(verbose_){
486            std::clog << "genJet matched Pt = " << gJp4.Pt() << std::endl;
487            std::clog << "genJet matched eta = " << gJp4.Eta() << std::endl;
# Line 521 | Line 521 | BTagSFContainer btagSFs;
521        TLorentzVector gJp4;
522        if(gJ){
523          gJp4.SetPtEtaPhiE(gJ->pt(),gJ->eta(),gJ->phi(),gJ->energy());
524 <        sj. bestMCp4mom = gJp4;
524 >        sj.bestMCp4mom = gJp4;
525          if(verbose_){
526            std::clog << "genJet matched Pt = " << gJp4.Pt() << std::endl;
527            std::clog << "genJet matched eta = " << gJp4.Eta() << std::endl;
# Line 601 | Line 601 | BTagSFContainer btagSFs;
601        TLorentzVector gJp4;
602        if(gJ){
603          gJp4.SetPtEtaPhiE(gJ->pt(),gJ->eta(),gJ->phi(),gJ->energy());
604 <        sj. bestMCp4mom = gJp4;
604 >        sj.bestMCp4mom = gJp4;
605          if(verbose_){
606            std::clog << "genJet matched Pt = " << gJp4.Pt() << std::endl;
607            std::clog << "genJet matched eta = " << gJp4.Eta() << std::endl;
# Line 732 | Line 732 | BTagSFContainer btagSFs;
732      if (verbose_)     std::cout <<" METTC "<<     hbbInfo->tcmet.metSig <<" " <<     hbbInfo->tcmet.sumEt<<std::endl;
733    }
734    
735 +  edm::Handle<edm::View<reco::MET> > pfMETNoPUHandle;
736 +  iEvent.getByLabel("pfMETNoPU",pfMETNoPUHandle);
737 +  edm::View<reco::MET> metspfMETNoPU = *pfMETNoPUHandle;
738 +  if(metspfMETNoPU.size()){
739 +    hbbInfo->metNoPU.sumEt=(metspfMETNoPU[0]).sumEt();
740 +    hbbInfo->metNoPU.metSig=(metspfMETNoPU[0]).significance();
741 +    hbbInfo->metNoPU.eLong=(metspfMETNoPU[0]).e_longitudinal();
742 +    hbbInfo->metNoPU.p4=GENPTOLOR((metspfMETNoPU[0]));
743 +    if (verbose_)     std::cout <<" pfMETNoPU "<<     hbbInfo->metNoPU.metSig <<" " <<     hbbInfo->metNoPU.sumEt<<std::endl;
744 +  }
745  
746    edm::Handle<edm::View<reco::MET> > mHTHandle;
747    iEvent.getByLabel("patMETsHT",mHTHandle);
# Line 744 | Line 754 | BTagSFContainer btagSFs;
754      if (verbose_)     std::cout <<" METHT "<<     hbbInfo->mht.metSig <<" " <<     hbbInfo->mht.sumEt<<std::endl;
755    }
756    
757 <  edm::Handle<edm::View<pat::MET> > metHandle;
757 >  edm::Handle<edm::View<reco::MET> > metHandle;
758    iEvent.getByLabel(metLabel_,metHandle);
759 <  edm::View<pat::MET> mets = *metHandle;
759 >  edm::View<reco::MET> mets = *metHandle;
760    
761    if(mets.size()){
762      hbbInfo->calomet.sumEt=(mets[0]).sumEt();
# Line 813 | Line 823 | BTagSFContainer btagSFs;
823        const reco::HitPattern& q = gTrack->hitPattern();
824        mf.globChi2=gTrack.get()->normalizedChi2();
825        mf.globNHits=q.numberOfValidMuonHits();
826 <      mf. validMuStations = q. muonStationsWithValidHits();
826 >      mf.validMuStations = q.muonStationsWithValidHits();
827      }else{
828        mf.globChi2=-99;
829        mf.globNHits=-99;
# Line 876 | Line 886 | BTagSFContainer btagSFs;
886      ef.pfPhoIso=elec->photonIso();
887      ef.pfNeuIso=elec->neutralHadronIso();
888  
889 <    Geom::Phi<double> deltaphi(elec->superCluster()->phi()-atan2(hbbInfo->calomet.p4.Py(),hbbInfo->calomet.p4.Px()));
889 >    Geom::Phi<double> deltaphi(elec->superCluster()->phi()-atan2(hbbInfo->pfmet.p4.Py(),hbbInfo->pfmet.p4.Px()));
890      ef.acop = deltaphi.value();
891      //
892      // fill eleids
# Line 935 | Line 945 | BTagSFContainer btagSFs;
945      tf.tIso=tau->trackIso();
946      tf.eIso=tau->ecalIso();
947      tf.hIso=tau->hcalIso();
948 <    Geom::Phi<double> deltaphi(tau->phi()-atan2(hbbInfo->calomet.p4.Py(),hbbInfo->calomet.p4.Px()));
948 >    Geom::Phi<double> deltaphi(tau->phi()-atan2(hbbInfo->pfmet.p4.Py(),hbbInfo->pfmet.p4.Px()));
949      double acop = deltaphi.value();
950      tf.acop=acop;
951      tf.idbyIso=tau->tauID("byIsolation");
# Line 1204 | Line 1214 | void HbbAnalyzerNew ::fillSimpleJet (VHb
1214        sj.neutralHadronEFraction = jet_iter-> neutralHadronEnergyFraction ();
1215        sj.chargedEmEFraction = jet_iter-> chargedEmEnergyFraction ();
1216        sj.neutralEmEFraction = jet_iter-> neutralEmEnergyFraction ();
1217 <      sj. nConstituents = jet_iter->getPFConstituents().size();
1217 >      sj.nConstituents = jet_iter->getPFConstituents().size();
1218        
1219      }
1220      //

Diff Legend

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