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.8 by tboccali, Tue Jun 21 15:49:01 2011 UTC vs.
Revision 1.9 by tboccali, Thu Jun 30 08:56:27 2011 UTC

# Line 41 | Line 41 | HbbAnalyzerNew::HbbAnalyzerNew(const edm
41    dimuLabel_(iConfig.getParameter<edm::InputTag>("dimuTag")),
42    dielecLabel_(iConfig.getParameter<edm::InputTag>("dielecTag")),
43    hltResults_(iConfig.getParameter<edm::InputTag>("hltResultsTag")),
44 <  runOnMC_(iConfig.getParameter<bool>("runOnMC")) {
44 >  runOnMC_(iConfig.getParameter<bool>("runOnMC")), verbose_(iConfig.getUntrackedParameter<bool>("verbose")) {
45  
46    //
47    // put the setwhatproduced etc etc
# Line 679 | Line 679 | HbbAnalyzerNew::produce(edm::Event& iEve
679    // met is calomet
680    //
681  
682 –
683 –  for(edm::View<pat::MET>::const_iterator met = mets.begin(); met!=mets.end(); ++met){
684 –    hbbInfo->calomet.sumEt=met->sumEt();
685 –    hbbInfo->calomet.metSig=met->mEtSig();
686 –    hbbInfo->calomet.eLong=met->e_longitudinal();
687 –    hbbInfo->calomet.fourMomentum=GENPTOLORP(met);
688 –  }
689 –
682    edm::Handle<edm::View<pat::MET> > metTCHandle;
683    iEvent.getByLabel("patMETsTC",metTCHandle);
684    edm::View<pat::MET> metsTC = *metTCHandle;
685 <  for(edm::View<pat::MET>::const_iterator metTC = metsTC.begin(); metTC!=metsTC.end(); ++metTC){
686 <    hbbInfo->calomet.sumEt=metTC->sumEt();
687 <    hbbInfo->calomet.metSig=metTC->mEtSig();
688 <    hbbInfo->calomet.eLong=metTC->e_longitudinal();
689 <    hbbInfo->calomet.fourMomentum=GENPTOLORP(metTC);
685 >  if(metsTC.size()){
686 >    hbbInfo->tcmet.sumEt=(metsTC[0]).sumEt();
687 >    hbbInfo->tcmet.metSig=(metsTC[0]).mEtSig();
688 >    hbbInfo->tcmet.eLong=(metsTC[0]).e_longitudinal();
689 >    hbbInfo->tcmet.fourMomentum=GENPTOLOR((metsTC[0]));
690 >    if (verbose_)     std::cout <<" METTC "<<     hbbInfo->tcmet.metSig <<" " <<     hbbInfo->tcmet.sumEt<<std::endl;
691 >  }
692 >  
693 >  if(mets.size()){
694 >    hbbInfo->calomet.sumEt=(mets[0]).sumEt();
695 >    hbbInfo->calomet.metSig=(mets[0]).mEtSig();
696 >    hbbInfo->calomet.eLong=(mets[0]).e_longitudinal();
697 >    hbbInfo->calomet.fourMomentum=GENPTOLOR((mets[0]));
698 >    if (verbose_)     std::cout <<" METTC "<<     hbbInfo->calomet.metSig <<" " <<     hbbInfo->calomet.sumEt<<std::endl;
699    }
700 +
701 +  
702    edm::Handle<edm::View<pat::MET> > metPFHandle;
703    iEvent.getByLabel("patMETsPF",metPFHandle);
704    edm::View<pat::MET> metsPF = *metPFHandle;
705 <  for(edm::View<pat::MET>::const_iterator metPF = metsPF.begin(); metPF!=metsPF.end(); ++metPF){
706 <    hbbInfo->calomet.sumEt=metPF->sumEt();
707 <    hbbInfo->calomet.metSig=metPF->mEtSig();
708 <    hbbInfo->calomet.eLong=metPF->e_longitudinal();
709 <    hbbInfo->calomet.fourMomentum=GENPTOLORP(metPF);
705 >
706 >  if(metsPF.size()){
707 >    hbbInfo->pfmet.sumEt=(metsPF[0]).sumEt();
708 >    hbbInfo->pfmet.metSig=(metsPF[0]).mEtSig();
709 >    hbbInfo->pfmet.eLong=(metsPF[0]).e_longitudinal();
710 >    hbbInfo->pfmet.fourMomentum=GENPTOLOR((metsPF[0]));
711 >    if (verbose_)     std::cout <<" METTC "<<     hbbInfo->pfmet.metSig <<" " <<     hbbInfo->pfmet.sumEt<<std::endl;
712    }
713  
714  
715 +  if(verbose_){
716 +    std::cout << "METs: calomet "<<mets.size()<<" tcmet "<<metsTC.size()<<" pfmet "<<metsPF.size()<<std::endl;  
717 +  }
718 +
719 +  std::cout << " INPUT MUONS "<<muons.size()<<std::endl;
720  
721    for(edm::View<pat::Muon>::const_iterator mu = muons.begin(); mu!=muons.end(); ++mu){
722      VHbbEvent::MuonInfo mf;
# Line 767 | Line 777 | HbbAnalyzerNew::produce(edm::Event& iEve
777      hbbInfo->muInfo.push_back(mf);
778    }
779  
780 <
780 >  std::cout << " INPUT electrons "<<electrons.size()<<std::endl;
781    for(edm::View<pat::Electron>::const_iterator elec = electrons.begin(); elec!=electrons.end(); ++elec){
782      VHbbEvent::ElectronInfo ef;
783      ef.fourMomentum=GENPTOLORP(elec);
# Line 801 | Line 811 | HbbAnalyzerNew::produce(edm::Event& iEve
811      hbbInfo->eleInfo.push_back(ef);
812    }
813  
814 +
815 +  std::cout << " INPUT taus "<<taus.size()<<std::endl;
816    for(edm::View<pat::Tau>::const_iterator tau = taus.begin(); tau!=taus.end(); ++tau){
817      VHbbEvent::TauInfo tf;
818      tf.fourMomentum=GENPTOLORP(tau);
# Line 949 | Line 961 | HbbAnalyzerNew::produce(edm::Event& iEve
961      hbbInfo->diElectronInfo.push_back(df);
962      
963    }
964 +   if (verbose_){
965 +     std::cout <<" Pushing hbbInfo "<<std::endl;
966 +     std::cout <<" SimpleJets1 = "<<hbbInfo->simpleJets.size()<<std::endl<<
967 +       " SimpleJets2 = "<<hbbInfo->simpleJets2.size()<<std::endl<<
968 +       " SubJets = "<<hbbInfo->subJets.size()<<std::endl<<
969 +       " HardJets = "<<hbbInfo->hardJets.size()<<std::endl<<
970 +       " Muons = "<<hbbInfo->muInfo.size()<<std::endl<<
971 +       " Electrons = "<<hbbInfo->eleInfo.size()<<std::endl<<
972 +       " Taus = "<<hbbInfo->tauInfo.size()<<std::endl<<
973 +       " Electrons = "<<hbbInfo->eleInfo.size()<<std::endl<<
974 +       "--------------------- "<<std::endl;
975 +  }
976 +
977 +
978    iEvent.put(hbbInfo);
979  
980 +
981   }
982    
983   void

Diff Legend

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