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.36 by arizzi, Thu Sep 15 10:50:35 2011 UTC vs.
Revision 1.41 by tboccali, Mon Sep 19 10:34:54 2011 UTC

# Line 176 | Line 176 | HbbAnalyzerNew::produce(edm::Event& iEve
176    edm::Handle<double> rhoHandle;
177    iEvent.getByLabel(edm::InputTag("kt6PFJets", "rho"),rhoHandle);  
178    auxInfo->puInfo.rho = *rhoHandle;
179 <  
179 >  
180 >  edm::Handle<double> rho25Handle;
181 >  iEvent.getByLabel(edm::InputTag("kt6PFJets25", "rho"),rho25Handle);  
182 >  auxInfo->puInfo.rho25 = *rho25Handle;
183 >
184    edm::Handle<std::vector< PileupSummaryInfo> > puHandle;
185  
186    if (runOnMC_){
# Line 437 | 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 >        std::cout <<" KKKKK 1 "<< sj.bestMCp4mom.Pt() <<std::endl;
446          if(verbose_){
447            std::clog << "genJet matched Pt = " << gJp4.Pt() << std::endl;
448            std::clog << "genJet matched eta = " << gJp4.Eta() << std::endl;
# Line 477 | Line 482 | BTagSFContainer btagSFs;
482        TLorentzVector gJp4;
483        if(gJ){
484          gJp4.SetPtEtaPhiE(gJ->pt(),gJ->eta(),gJ->phi(),gJ->energy());
485 <        sj. bestMCp4mom = gJp4;
485 >        sj.bestMCp4mom = gJp4;
486 >        std::cout <<" KKKKK 2 "<< (sj.bestMCp4mom.Pt()) <<std::endl;
487          if(verbose_){
488            std::clog << "genJet matched Pt = " << gJp4.Pt() << std::endl;
489            std::clog << "genJet matched eta = " << gJp4.Eta() << std::endl;
# Line 517 | Line 523 | BTagSFContainer btagSFs;
523        TLorentzVector gJp4;
524        if(gJ){
525          gJp4.SetPtEtaPhiE(gJ->pt(),gJ->eta(),gJ->phi(),gJ->energy());
526 <        sj. bestMCp4mom = gJp4;
526 >        sj.bestMCp4mom = gJp4;
527 >        std::cout <<" KKKKK 3 "<< sj.bestMCp4mom.Pt() <<std::endl;
528          if(verbose_){
529            std::clog << "genJet matched Pt = " << gJp4.Pt() << std::endl;
530            std::clog << "genJet matched eta = " << gJp4.Eta() << std::endl;
# Line 597 | Line 604 | BTagSFContainer btagSFs;
604        TLorentzVector gJp4;
605        if(gJ){
606          gJp4.SetPtEtaPhiE(gJ->pt(),gJ->eta(),gJ->phi(),gJ->energy());
607 <        sj. bestMCp4mom = gJp4;
607 >        sj.bestMCp4mom = gJp4;
608 >        std::cout <<" KKKKK 4 "<< sj.bestMCp4mom.Pt() <<std::endl;
609          if(verbose_){
610            std::clog << "genJet matched Pt = " << gJp4.Pt() << std::endl;
611            std::clog << "genJet matched eta = " << gJp4.Eta() << std::endl;
# Line 728 | Line 736 | BTagSFContainer btagSFs;
736      if (verbose_)     std::cout <<" METTC "<<     hbbInfo->tcmet.metSig <<" " <<     hbbInfo->tcmet.sumEt<<std::endl;
737    }
738    
739 +  edm::Handle<edm::View<reco::MET> > pfMETNoPUHandle;
740 +  iEvent.getByLabel("pfMETNoPU",pfMETNoPUHandle);
741 +  edm::View<reco::MET> metspfMETNoPU = *pfMETNoPUHandle;
742 +  if(metspfMETNoPU.size()){
743 +    hbbInfo->metNoPU.sumEt=(metspfMETNoPU[0]).sumEt();
744 +    hbbInfo->metNoPU.metSig=(metspfMETNoPU[0]).significance();
745 +    hbbInfo->metNoPU.eLong=(metspfMETNoPU[0]).e_longitudinal();
746 +    hbbInfo->metNoPU.p4=GENPTOLOR((metspfMETNoPU[0]));
747 +    if (verbose_)     std::cout <<" pfMETNoPU "<<     hbbInfo->metNoPU.metSig <<" " <<     hbbInfo->metNoPU.sumEt<<std::endl;
748 +  }
749  
750    edm::Handle<edm::View<reco::MET> > mHTHandle;
751    iEvent.getByLabel("patMETsHT",mHTHandle);
# Line 740 | Line 758 | BTagSFContainer btagSFs;
758      if (verbose_)     std::cout <<" METHT "<<     hbbInfo->mht.metSig <<" " <<     hbbInfo->mht.sumEt<<std::endl;
759    }
760    
761 <  edm::Handle<edm::View<pat::MET> > metHandle;
761 >  edm::Handle<edm::View<reco::MET> > metHandle;
762    iEvent.getByLabel(metLabel_,metHandle);
763 <  edm::View<pat::MET> mets = *metHandle;
763 >  edm::View<reco::MET> mets = *metHandle;
764    
765    if(mets.size()){
766      hbbInfo->calomet.sumEt=(mets[0]).sumEt();
# Line 780 | Line 798 | BTagSFContainer btagSFs;
798      mf.eIso=mu->ecalIso();
799      mf.hIso=mu->hcalIso();
800      mf.pfChaIso=mu->chargedHadronIso();
801 +    mf.pfChaPUIso=mu->userIso(5);
802      mf.pfPhoIso=mu->photonIso();
803      mf.pfNeuIso=mu->neutralHadronIso();
804      Geom::Phi<double> deltaphi(mu->phi()-atan2(mf.p4.Px(), mf.p4.Py()));
# Line 808 | Line 827 | BTagSFContainer btagSFs;
827        const reco::HitPattern& q = gTrack->hitPattern();
828        mf.globChi2=gTrack.get()->normalizedChi2();
829        mf.globNHits=q.numberOfValidMuonHits();
830 <      mf. validMuStations = q. muonStationsWithValidHits();
830 >      mf.validMuStations = q.muonStationsWithValidHits();
831      }else{
832        mf.globChi2=-99;
833        mf.globNHits=-99;
# Line 867 | Line 886 | BTagSFContainer btagSFs;
886      ef.eIso=elec->ecalIso();
887      ef.hIso=elec->hcalIso();
888      ef.pfChaIso=elec->chargedHadronIso();
889 +    ef.pfChaPUIso=elec->userIso(5);
890      ef.pfPhoIso=elec->photonIso();
891      ef.pfNeuIso=elec->neutralHadronIso();
892  
893 <    Geom::Phi<double> deltaphi(elec->superCluster()->phi()-atan2(hbbInfo->calomet.p4.Py(),hbbInfo->calomet.p4.Px()));
893 >    Geom::Phi<double> deltaphi(elec->superCluster()->phi()-atan2(hbbInfo->pfmet.p4.Py(),hbbInfo->pfmet.p4.Px()));
894      ef.acop = deltaphi.value();
895      //
896      // fill eleids
# Line 929 | Line 949 | BTagSFContainer btagSFs;
949      tf.tIso=tau->trackIso();
950      tf.eIso=tau->ecalIso();
951      tf.hIso=tau->hcalIso();
952 <    Geom::Phi<double> deltaphi(tau->phi()-atan2(hbbInfo->calomet.p4.Py(),hbbInfo->calomet.p4.Px()));
952 >    Geom::Phi<double> deltaphi(tau->phi()-atan2(hbbInfo->pfmet.p4.Py(),hbbInfo->pfmet.p4.Px()));
953      double acop = deltaphi.value();
954      tf.acop=acop;
955      tf.idbyIso=tau->tauID("byIsolation");
# Line 1198 | Line 1218 | void HbbAnalyzerNew ::fillSimpleJet (VHb
1218        sj.neutralHadronEFraction = jet_iter-> neutralHadronEnergyFraction ();
1219        sj.chargedEmEFraction = jet_iter-> chargedEmEnergyFraction ();
1220        sj.neutralEmEFraction = jet_iter-> neutralEmEnergyFraction ();
1221 <      sj. nConstituents = jet_iter->getPFConstituents().size();
1221 >      sj.nConstituents = jet_iter->getPFConstituents().size();
1222        
1223      }
1224      //

Diff Legend

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