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.52 by tboccali, Thu Oct 6 09:19:18 2011 UTC

# Line 171 | Line 171 | HbbAnalyzerNew::produce(edm::Event& iEve
171    
172    auxInfo->pvInfo.firstPVInPT2 = TVector3(RecVtxFirst.x(), RecVtxFirst.y(), RecVtxFirst.z());
173    auxInfo->pvInfo.firstPVInProb = TVector3(RecVtx.x(), RecVtx.y(), RecVtx.z());
174 <
174 >  
175 >  (auxInfo->pvInfo).efirstPVInPT2 = (RecVtxFirst.error());
176 >  (auxInfo->pvInfo).efirstPVInProb = RecVtx.error();
177      
178    edm::Handle<double> rhoHandle;
179    iEvent.getByLabel(edm::InputTag("kt6PFJets", "rho"),rhoHandle);  
# Line 420 | Line 422 | BTagSFContainer btagSFs;
422      //     if(jet_iter->pt()>50)
423      //       njetscounter++;
424      VHbbEvent::SimpleJet sj;
425 +    //    std::cout <<" sj1"<<std::endl;
426      fillSimpleJet(sj,jet_iter);
427      //    if(!runOnMC_)
428  
# Line 435 | Line 438 | BTagSFContainer btagSFs;
438        //genJet
439        const reco::GenJet *gJ = jet_iter->genJet();
440        //physical parton for mother info ONLY
441 <      if( (jet_iter->genParton())
442 <          and (jet_iter->genParton()->mother()) )
443 <        sj.bestMCmomid=jet_iter->genParton()->mother()->pdgId();
441 >      if( (jet_iter->genParton()) ){
442 >        sj.bestMCid = jet_iter->genParton()->pdgId();
443 >        if( (jet_iter->genParton()->mother()) )
444 >          sj.bestMCmomid=jet_iter->genParton()->mother()->pdgId();
445 >      }
446        TLorentzVector gJp4;
447        if(gJ){
448          gJp4.SetPtEtaPhiE(gJ->pt(),gJ->eta(),gJ->phi(),gJ->energy());
449 <        sj.bestMCp4mom = gJp4;
445 <        std::cout <<" KKKKK 1 "<< sj.bestMCp4mom.Pt() <<std::endl;
449 >        sj.bestMCp4 = gJp4;
450          if(verbose_){
451            std::clog << "genJet matched Pt = " << gJp4.Pt() << std::endl;
452            std::clog << "genJet matched eta = " << gJp4.Eta() << std::endl;
# Line 461 | Line 465 | BTagSFContainer btagSFs;
465      //     if(jet_iter->pt()>50)
466      //       njetscounter++;
467      VHbbEvent::SimpleJet sj;
468 +    //    std::cout <<" sj3"<<std::endl;
469      fillSimpleJet(sj,jet_iter);
470      //    if(!runOnMC_)  
471    setJecUnc(sj,jecUnc);
# Line 476 | Line 481 | BTagSFContainer btagSFs;
481        //genJet
482        const reco::GenJet *gJ = jet_iter->genJet();
483        //physical parton for mother info ONLY
484 <      if( (jet_iter->genParton())
485 <          and (jet_iter->genParton()->mother()) )
486 <        sj.bestMCmomid=jet_iter->genParton()->mother()->pdgId();
484 >      if( (jet_iter->genParton()) ){
485 >        sj.bestMCid = jet_iter->genParton()->pdgId();
486 >        if( (jet_iter->genParton()->mother()) )
487 >          sj.bestMCmomid=jet_iter->genParton()->mother()->pdgId();
488 >      }
489        TLorentzVector gJp4;
490        if(gJ){
491          gJp4.SetPtEtaPhiE(gJ->pt(),gJ->eta(),gJ->phi(),gJ->energy());
492 <        sj.bestMCp4mom = gJp4;
486 <        std::cout <<" KKKKK 2 "<< (sj.bestMCp4mom.Pt()) <<std::endl;
492 >        sj.bestMCp4 = gJp4;
493          if(verbose_){
494            std::clog << "genJet matched Pt = " << gJp4.Pt() << std::endl;
495            std::clog << "genJet matched eta = " << gJp4.Eta() << std::endl;
# Line 502 | Line 508 | BTagSFContainer btagSFs;
508      //     if(jet_iter->pt()>50)
509      //       njetscounter++;
510      VHbbEvent::SimpleJet sj;
511 +    //    std::cout <<" sj4"<<std::endl;
512      fillSimpleJet(sj,jet_iter);
513      //    if(!runOnMC_)  
514      setJecUnc(sj,jecUnc);
# Line 517 | Line 524 | BTagSFContainer btagSFs;
524        //genJet
525        const reco::GenJet *gJ = jet_iter->genJet();
526        //physical parton for mother info ONLY
527 <      if( (jet_iter->genParton())
528 <          and (jet_iter->genParton()->mother()) )
529 <        sj.bestMCmomid=jet_iter->genParton()->mother()->pdgId();
527 >      if( (jet_iter->genParton()) ){
528 >        sj.bestMCid = jet_iter->genParton()->pdgId();
529 >        if( (jet_iter->genParton()->mother()) )
530 >          sj.bestMCmomid=jet_iter->genParton()->mother()->pdgId();
531 >      }
532        TLorentzVector gJp4;
533        if(gJ){
534          gJp4.SetPtEtaPhiE(gJ->pt(),gJ->eta(),gJ->phi(),gJ->energy());
535 <        sj.bestMCp4mom = gJp4;
527 <        std::cout <<" KKKKK 3 "<< sj.bestMCp4mom.Pt() <<std::endl;
535 >        sj.bestMCp4 = gJp4;
536          if(verbose_){
537            std::clog << "genJet matched Pt = " << gJp4.Pt() << std::endl;
538            std::clog << "genJet matched eta = " << gJp4.Eta() << std::endl;
# Line 543 | Line 551 | BTagSFContainer btagSFs;
551    for(edm::View<pat::Jet>::const_iterator jet_iter = simplejets2.begin(); jet_iter!=simplejets2.end(); ++jet_iter){
552      
553      VHbbEvent::SimpleJet sj;
554 +    //    std::cout <<" sj2"<<std::endl;
555      fillSimpleJet(sj,jet_iter);    
556      //  if(!runOnMC_)  
557   setJecUnc(sj,jecUnc);
# Line 598 | Line 607 | BTagSFContainer btagSFs;
607        //genJet
608        const reco::GenJet *gJ = jet_iter->genJet();
609        //physical parton for mother info ONLY
610 <      if( (jet_iter->genParton())
611 <          and (jet_iter->genParton()->mother()) )
612 <        sj.bestMCmomid=jet_iter->genParton()->mother()->pdgId();
610 >      if( (jet_iter->genParton()) ){
611 >        sj.bestMCid = jet_iter->genParton()->pdgId();
612 >        if( (jet_iter->genParton()->mother()) )
613 >          sj.bestMCmomid=jet_iter->genParton()->mother()->pdgId();
614 >      }
615        TLorentzVector gJp4;
616        if(gJ){
617          gJp4.SetPtEtaPhiE(gJ->pt(),gJ->eta(),gJ->phi(),gJ->energy());
618 <        sj.bestMCp4mom = gJp4;
608 <        std::cout <<" KKKKK 4 "<< sj.bestMCp4mom.Pt() <<std::endl;
618 >        sj.bestMCp4 = gJp4;
619          if(verbose_){
620            std::clog << "genJet matched Pt = " << gJp4.Pt() << std::endl;
621            std::clog << "genJet matched eta = " << gJp4.Eta() << std::endl;
# Line 685 | Line 695 | BTagSFContainer btagSFs;
695                              << "," << subjet_iter->bDiscriminator("combinedSecondaryVertexBJetTags") << "\n";}
696  
697      VHbbEvent::SimpleJet sj;
698 +    //    std::cout <<" sub jet "<<std::endl;
699      fillSimpleJet(sj,subjet_iter);
700      //  if(!runOnMC_)  
701      setJecUnc(sj,jecUnc);
# Line 805 | Line 816 | BTagSFContainer btagSFs;
816      double acop = deltaphi.value();
817      mf.acop=acop;
818  
819 +    mf.emEnergy = mu->calEnergy().em;
820 +    mf.hadEnergy = mu->calEnergy().had;
821 +
822 +    mf.nMatches = mu->numberOfMatches();
823 +
824      mf.ipDb=mu->dB();
825      mf.ipErrDb=mu->edB();
826      mf.cat=0;
# Line 821 | Line 837 | BTagSFContainer btagSFs;
837        TrackRef iTrack1 = mu->innerTrack();
838        const reco::HitPattern& p1 = iTrack1->hitPattern();
839        mf.nPixelHits=p1.pixelLayersWithMeasurement();
840 <    }
840 >
841 >      mf.nValidTracker = p1.numberOfValidTrackerHits();
842 >      mf.nValidPixel = p1.numberOfValidPixelHits();
843 >
844 >
845 >
846 >   }
847      if(mu->isGlobalMuon()){
848        TrackRef gTrack = mu->globalTrack();
849        const reco::HitPattern& q = gTrack->hitPattern();
# Line 890 | Line 912 | BTagSFContainer btagSFs;
912      ef.pfPhoIso=elec->photonIso();
913      ef.pfNeuIso=elec->neutralHadronIso();
914  
915 +    //
916 +    // ip info
917 +    //
918 +
919 +    ef.ipDb=elec->dB();
920 +    ef.ipErrDb=elec->edB();
921 +    
922 +
923 +
924      Geom::Phi<double> deltaphi(elec->superCluster()->phi()-atan2(hbbInfo->pfmet.p4.Py(),hbbInfo->pfmet.p4.Px()));
925      ef.acop = deltaphi.value();
926      //
927 +    ef.sihih = elec->sigmaIetaIeta();
928 +    ef.Dphi = elec->deltaPhiSuperClusterTrackAtVtx();
929 +    ef.Deta = elec->deltaEtaSuperClusterTrackAtVtx();
930 +    ef.HoE = elec->hadronicOverEm();
931 +    ef.convDist = elec->convDist();
932 +    ef.convDcot = elec->convDcot();
933 +    if(elec->gsfTrack().isNonnull()) ef.innerHits = elec->gsfTrack()->trackerExpectedHitsInner().numberOfHits();  
934 +    ef.isEB = elec->isEB();
935 +    ef.isEE = elec->isEE();
936 +    //
937      // fill eleids
938      //    
939   /*    ef.id95 = elec->electronID("simpleEleId95cIso");
# Line 1179 | Line 1220 | void HbbAnalyzerNew::setJecUnc(VHbbEvent
1220    // test
1221    //
1222  
1223 <  return;
1223 >  //  return;
1224    double eta = sj.p4.Eta();
1225    double pt = sj.p4.Pt();
1226    
# Line 1199 | Line 1240 | void HbbAnalyzerNew ::fillSimpleJet (VHb
1240      sj.jpb=jet_iter->bDiscriminator("jetBProbabilityBJetTags");
1241      sj.ssvhe=jet_iter->bDiscriminator("simpleSecondaryVertexHighEffBJetTags");
1242      sj.csv=jet_iter->bDiscriminator("combinedSecondaryVertexBJetTags");
1243 <    sj.csvmva=jet_iter->bDiscriminator("combinedSecondaryVertexMVABJetTags");
1243 >   sj.csvmva=jet_iter->bDiscriminator("combinedSecondaryVertexMVABJetTags");
1244      sj.charge=jet_iter->jetCharge();
1245      sj.ntracks=jet_iter->associatedTracks().size();
1246      sj.p4=GENPTOLORP(jet_iter);
1247 +    //    std::cout << " ECCO "<<sj.csv<< " "<< sj.p4.Pt()<<std::endl;
1248      sj.chargedTracksFourMomentum=(getChargedTracksMomentum(&*(jet_iter)));
1249      sj.SF_CSVL=1;
1250      sj.SF_CSVM=1;

Diff Legend

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