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.41 by tboccali, Mon Sep 19 10:34:54 2011 UTC vs.
Revision 1.45 by tboccali, Mon Sep 19 14:46:12 2011 UTC

# Line 420 | Line 420 | BTagSFContainer btagSFs;
420      //     if(jet_iter->pt()>50)
421      //       njetscounter++;
422      VHbbEvent::SimpleJet sj;
423 +    //    std::cout <<" sj1"<<std::endl;
424      fillSimpleJet(sj,jet_iter);
425      //    if(!runOnMC_)
426  
# Line 435 | Line 436 | BTagSFContainer btagSFs;
436        //genJet
437        const reco::GenJet *gJ = jet_iter->genJet();
438        //physical parton for mother info ONLY
439 <      if( (jet_iter->genParton())
440 <          and (jet_iter->genParton()->mother()) )
441 <        sj.bestMCmomid=jet_iter->genParton()->mother()->pdgId();
439 >      if( (jet_iter->genParton()) ){
440 >        sj.bestMCid = jet_iter->genParton()->pdgId();
441 >        if( (jet_iter->genParton()->mother()) )
442 >          sj.bestMCmomid=jet_iter->genParton()->mother()->pdgId();
443 >      }
444        TLorentzVector gJp4;
445        if(gJ){
446          gJp4.SetPtEtaPhiE(gJ->pt(),gJ->eta(),gJ->phi(),gJ->energy());
447 <        sj.bestMCp4mom = gJp4;
445 <        std::cout <<" KKKKK 1 "<< sj.bestMCp4mom.Pt() <<std::endl;
447 >        sj.bestMCp4 = gJp4;
448          if(verbose_){
449            std::clog << "genJet matched Pt = " << gJp4.Pt() << std::endl;
450            std::clog << "genJet matched eta = " << gJp4.Eta() << std::endl;
# Line 461 | Line 463 | BTagSFContainer btagSFs;
463      //     if(jet_iter->pt()>50)
464      //       njetscounter++;
465      VHbbEvent::SimpleJet sj;
466 +    //    std::cout <<" sj3"<<std::endl;
467      fillSimpleJet(sj,jet_iter);
468      //    if(!runOnMC_)  
469    setJecUnc(sj,jecUnc);
# Line 476 | Line 479 | BTagSFContainer btagSFs;
479        //genJet
480        const reco::GenJet *gJ = jet_iter->genJet();
481        //physical parton for mother info ONLY
482 <      if( (jet_iter->genParton())
483 <          and (jet_iter->genParton()->mother()) )
484 <        sj.bestMCmomid=jet_iter->genParton()->mother()->pdgId();
482 >      if( (jet_iter->genParton()) ){
483 >        sj.bestMCid = jet_iter->genParton()->pdgId();
484 >        if( (jet_iter->genParton()->mother()) )
485 >          sj.bestMCmomid=jet_iter->genParton()->mother()->pdgId();
486 >      }
487        TLorentzVector gJp4;
488        if(gJ){
489          gJp4.SetPtEtaPhiE(gJ->pt(),gJ->eta(),gJ->phi(),gJ->energy());
490 <        sj.bestMCp4mom = gJp4;
486 <        std::cout <<" KKKKK 2 "<< (sj.bestMCp4mom.Pt()) <<std::endl;
490 >        sj.bestMCp4 = gJp4;
491          if(verbose_){
492            std::clog << "genJet matched Pt = " << gJp4.Pt() << std::endl;
493            std::clog << "genJet matched eta = " << gJp4.Eta() << std::endl;
# Line 502 | Line 506 | BTagSFContainer btagSFs;
506      //     if(jet_iter->pt()>50)
507      //       njetscounter++;
508      VHbbEvent::SimpleJet sj;
509 +    //    std::cout <<" sj4"<<std::endl;
510      fillSimpleJet(sj,jet_iter);
511      //    if(!runOnMC_)  
512      setJecUnc(sj,jecUnc);
# Line 517 | Line 522 | BTagSFContainer btagSFs;
522        //genJet
523        const reco::GenJet *gJ = jet_iter->genJet();
524        //physical parton for mother info ONLY
525 <      if( (jet_iter->genParton())
526 <          and (jet_iter->genParton()->mother()) )
527 <        sj.bestMCmomid=jet_iter->genParton()->mother()->pdgId();
525 >      if( (jet_iter->genParton()) ){
526 >        sj.bestMCid = jet_iter->genParton()->pdgId();
527 >        if( (jet_iter->genParton()->mother()) )
528 >          sj.bestMCmomid=jet_iter->genParton()->mother()->pdgId();
529 >      }
530        TLorentzVector gJp4;
531        if(gJ){
532          gJp4.SetPtEtaPhiE(gJ->pt(),gJ->eta(),gJ->phi(),gJ->energy());
533 <        sj.bestMCp4mom = gJp4;
527 <        std::cout <<" KKKKK 3 "<< sj.bestMCp4mom.Pt() <<std::endl;
533 >        sj.bestMCp4 = gJp4;
534          if(verbose_){
535            std::clog << "genJet matched Pt = " << gJp4.Pt() << std::endl;
536            std::clog << "genJet matched eta = " << gJp4.Eta() << std::endl;
# Line 543 | Line 549 | BTagSFContainer btagSFs;
549    for(edm::View<pat::Jet>::const_iterator jet_iter = simplejets2.begin(); jet_iter!=simplejets2.end(); ++jet_iter){
550      
551      VHbbEvent::SimpleJet sj;
552 +    //    std::cout <<" sj2"<<std::endl;
553      fillSimpleJet(sj,jet_iter);    
554      //  if(!runOnMC_)  
555   setJecUnc(sj,jecUnc);
# Line 598 | Line 605 | BTagSFContainer btagSFs;
605        //genJet
606        const reco::GenJet *gJ = jet_iter->genJet();
607        //physical parton for mother info ONLY
608 <      if( (jet_iter->genParton())
609 <          and (jet_iter->genParton()->mother()) )
610 <        sj.bestMCmomid=jet_iter->genParton()->mother()->pdgId();
608 >      if( (jet_iter->genParton()) ){
609 >        sj.bestMCid = jet_iter->genParton()->pdgId();
610 >        if( (jet_iter->genParton()->mother()) )
611 >          sj.bestMCmomid=jet_iter->genParton()->mother()->pdgId();
612 >      }
613        TLorentzVector gJp4;
614        if(gJ){
615          gJp4.SetPtEtaPhiE(gJ->pt(),gJ->eta(),gJ->phi(),gJ->energy());
616 <        sj.bestMCp4mom = gJp4;
608 <        std::cout <<" KKKKK 4 "<< sj.bestMCp4mom.Pt() <<std::endl;
616 >        sj.bestMCp4 = gJp4;
617          if(verbose_){
618            std::clog << "genJet matched Pt = " << gJp4.Pt() << std::endl;
619            std::clog << "genJet matched eta = " << gJp4.Eta() << std::endl;
# Line 685 | Line 693 | BTagSFContainer btagSFs;
693                              << "," << subjet_iter->bDiscriminator("combinedSecondaryVertexBJetTags") << "\n";}
694  
695      VHbbEvent::SimpleJet sj;
696 +    //    std::cout <<" sub jet "<<std::endl;
697      fillSimpleJet(sj,subjet_iter);
698      //  if(!runOnMC_)  
699      setJecUnc(sj,jecUnc);
# Line 805 | Line 814 | BTagSFContainer btagSFs;
814      double acop = deltaphi.value();
815      mf.acop=acop;
816  
817 +    mf.nMatches = mu->numberOfMatches();
818 +
819 +
820      mf.ipDb=mu->dB();
821      mf.ipErrDb=mu->edB();
822      mf.cat=0;
# Line 1199 | Line 1211 | void HbbAnalyzerNew ::fillSimpleJet (VHb
1211      sj.jpb=jet_iter->bDiscriminator("jetBProbabilityBJetTags");
1212      sj.ssvhe=jet_iter->bDiscriminator("simpleSecondaryVertexHighEffBJetTags");
1213      sj.csv=jet_iter->bDiscriminator("combinedSecondaryVertexBJetTags");
1214 <    sj.csvmva=jet_iter->bDiscriminator("combinedSecondaryVertexMVABJetTags");
1214 >   sj.csvmva=jet_iter->bDiscriminator("combinedSecondaryVertexMVABJetTags");
1215      sj.charge=jet_iter->jetCharge();
1216      sj.ntracks=jet_iter->associatedTracks().size();
1217      sj.p4=GENPTOLORP(jet_iter);
1218 +    //    std::cout << " ECCO "<<sj.csv<< " "<< sj.p4.Pt()<<std::endl;
1219      sj.chargedTracksFourMomentum=(getChargedTracksMomentum(&*(jet_iter)));
1220      sj.SF_CSVL=1;
1221      sj.SF_CSVM=1;

Diff Legend

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