ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/UHHAnalysis/NtupleWriter/src/NtupleWriter.cc
(Generate patch)

Comparing UserCode/UHHAnalysis/NtupleWriter/src/NtupleWriter.cc (file contents):
Revision 1.25 by peiffer, Thu Jun 6 07:38:31 2013 UTC vs.
Revision 1.26 by peiffer, Mon Jun 10 13:33:05 2013 UTC

# Line 311 | Line 311 | NtupleWriter::analyze(const edm::Event&
311       for(reco::GenParticleCollection::const_iterator iter = genPartColl->begin(); iter != genPartColl->end(); ++ iter){
312         index++;
313        
314 <       //write out only top quarks and status 3 particles (works fine only for MadGraph)
315 <       if(abs(iter->pdgId())==6 || iter->status()==3 || doAllGenParticles){
314 >       //write out only top quarks,final state leptons and status 3 particles (works fine only for MadGraph)
315 >       bool islepton = iter->status()==1 && abs(iter->pdgId())>=11 && abs(iter->pdgId())<=16 ;
316 >       if(abs(iter->pdgId())==6 || iter->status()==3 || islepton ||  doAllGenParticles){
317           GenParticle genp;
318           genp.set_charge(iter->charge());
319           genp.set_pt(iter->p4().pt());
# Line 527 | Line 528 | NtupleWriter::analyze(const edm::Event&
528           tau.set_phi( pat_tau.phi());
529           tau.set_energy( pat_tau.energy());
530           tau.set_decayModeFinding ( pat_tau.tauID("decayModeFinding")>0.5);
531 <         tau.set_byVLooseCombinedIsolationDeltaBetaCorr  ( pat_tau.tauID("byVLooseCombinedIsolationDeltaBetaCorr")>0.5);
531 >         //tau.set_byVLooseCombinedIsolationDeltaBetaCorr  ( pat_tau.tauID("byVLooseCombinedIsolationDeltaBetaCorr")>0.5);
532           tau.set_byLooseCombinedIsolationDeltaBetaCorr ( pat_tau.tauID("byLooseCombinedIsolationDeltaBetaCorr")>0.5);
533           tau.set_byMediumCombinedIsolationDeltaBetaCorr ( pat_tau.tauID("byMediumCombinedIsolationDeltaBetaCorr")>0.5);
534           tau.set_byTightCombinedIsolationDeltaBetaCorr ( pat_tau.tauID("byTightCombinedIsolationDeltaBetaCorr")>0.5);
535 <         tau.set_againstElectronLoose  ( pat_tau.tauID("againstElectronLoose")>0.5);
536 <         tau.set_againstElectronMedium ( pat_tau.tauID("againstElectronMedium")>0.5);
537 <         tau.set_againstElectronTight ( pat_tau.tauID("againstElectronTight")>0.5);
538 <         tau.set_againstElectronMVA  ( pat_tau.tauID("againstElectronMVA")>0.5);
539 <         tau.set_againstMuonLoose ( pat_tau.tauID("againstMuonLoose")>0.5);
540 <         tau.set_againstMuonMedium ( pat_tau.tauID("againstMuonMedium")>0.5);
541 <         tau.set_againstMuonTight ( pat_tau.tauID("againstMuonTight")>0.5);
535 >         tau.set_byLooseIsolationMVA( pat_tau.tauID("byLooseIsolationMVA")>0.5);
536 >         tau.set_byMediumIsolationMVA( pat_tau.tauID("byMediumIsolationMVA")>0.5);
537 >         tau.set_byTightIsolationMVA( pat_tau.tauID("byTightIsolationMVA")>0.5);
538 >         tau.set_byLooseIsolationMVA2( pat_tau.tauID("byLooseIsolationMVA2")>0.5);
539 >         tau.set_byMediumIsolationMVA2( pat_tau.tauID("byMediumIsolationMVA2")>0.5);
540 >         tau.set_byTightIsolationMVA2( pat_tau.tauID("byTightIsolationMVA2")>0.5);
541 >         tau.set_byLooseCombinedIsolationDeltaBetaCorr3Hits(  pat_tau.tauID("byLooseCombinedIsolationDeltaBetaCorr3Hits")>0.5);
542 >         tau.set_byMediumCombinedIsolationDeltaBetaCorr3Hits ( pat_tau.tauID("byMediumCombinedIsolationDeltaBetaCorr3Hits")>0.5);
543 >         tau.set_byTightCombinedIsolationDeltaBetaCorr3Hits ( pat_tau.tauID("byTightCombinedIsolationDeltaBetaCorr3Hits")>0.5);
544 >         tau.set_againstElectronLooseMVA3  ( pat_tau.tauID("againstElectronLooseMVA3")>0.5);
545 >         tau.set_againstElectronMediumMVA3 ( pat_tau.tauID("againstElectronMediumMVA3")>0.5);
546 >         tau.set_againstElectronTightMVA3 ( pat_tau.tauID("againstElectronTightMVA3")>0.5);
547 >         tau.set_againstElectronVTightMVA3 ( pat_tau.tauID("againstElectronVTightMVA3")>0.5);
548 >         tau.set_againstMuonLoose2 ( pat_tau.tauID("againstMuonLoose2")>0.5);
549 >         tau.set_againstMuonMedium2 ( pat_tau.tauID("againstMuonMedium2")>0.5);
550 >         tau.set_againstMuonTight2 ( pat_tau.tauID("againstMuonTight2")>0.5);
551 >         tau.set_byIsolationMVAraw(  pat_tau.tauID("byIsolationMVAraw"));
552 >         tau.set_byIsolationMVA2raw(  pat_tau.tauID("byIsolationMVA2raw"));
553 >         tau.set_decayMode( pat_tau.decayMode() );
554 >         tau.set_byCombinedIsolationDeltaBetaCorrRaw( pat_tau.tauID("byCombinedIsolationDeltaBetaCorrRaw"));
555 >         tau.set_byCombinedIsolationDeltaBetaCorrRaw3Hits( pat_tau.tauID("byCombinedIsolationDeltaBetaCorrRaw3Hits"));
556  
557 + //       std::cout << pat_tau.tauID("byCombinedIsolationDeltaBetaCorrRaw3Hits") << std::endl;
558 +        
559   //       reco::PFCandidateRef leadPFCand = pat_tau.leadPFCand();
560   //       if(!leadPFCand.isNull()){
561   //         tau.set_leadPFCand_px ( leadPFCand->px());

Diff Legend

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