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.21 by peiffer, Tue Jun 26 08:13:28 2012 UTC vs.
Revision 1.22 by peiffer, Wed Jul 25 09:56:57 2012 UTC

# Line 143 | Line 143 | NtupleWriter::NtupleWriter(const edm::Pa
143      }
144    }
145    if(doGenInfo){
146 +    genparticle_source= iConfig.getParameter<edm::InputTag>("genparticle_source");
147      tr->Branch("genInfo","GenInfo",&genInfo);
148      tr->Branch("GenParticles","std::vector<GenParticle>", &genps);
149    }
# Line 292 | Line 293 | NtupleWriter::analyze(const edm::Event&
293       }
294  
295       edm::Handle<reco::GenParticleCollection> genPartColl;
296 <     iEvent.getByLabel(edm::InputTag("genParticles"), genPartColl);
296 >     iEvent.getByLabel(genparticle_source, genPartColl);
297       int index=-1;
298       for(reco::GenParticleCollection::const_iterator iter = genPartColl->begin(); iter != genPartColl->end(); ++ iter){
299         index++;
# Line 383 | Line 384 | NtupleWriter::analyze(const edm::Event&
384           ele.set_EcalEnergy(pat_ele.ecalEnergy());
385           ele.set_mvaTrigV0(pat_ele.electronID("mvaTrigV0"));
386           ele.set_mvaNonTrigV0(pat_ele.electronID("mvaNonTrigV0"));
387 +         float AEff03 = 0.00;
388 +         if(isRealData){
389 +           AEff03 = ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleGammaAndNeutralHadronIso03, pat_ele.superCluster()->eta(), ElectronEffectiveArea::kEleEAData2011);
390 +         }else{
391 +           AEff03 = ElectronEffectiveArea::GetElectronEffectiveArea(ElectronEffectiveArea::kEleGammaAndNeutralHadronIso03, pat_ele.superCluster()->eta(), ElectronEffectiveArea::kEleEAFall11MC);
392 +         }
393 +         ele.set_AEff(AEff03);
394  
395           eles[j].push_back(ele);
396         }
# Line 592 | Line 600 | NtupleWriter::analyze(const edm::Event&
600           jet.set_btag_jetBProbability(pat_jet.bDiscriminator("jetBProbabilityBJetTags"));
601           jet.set_btag_jetProbability(pat_jet.bDiscriminator("jetProbabilityBJetTags"));
602  
603 +        
604           const reco::GenJet *genj = pat_jet.genJet();
605           if(genj){
606             jet.set_genjet_pt(genj->pt());
# Line 611 | Line 620 | NtupleWriter::analyze(const edm::Event&
620               if(jet.genparticles_indices().size()!= jetgenps.size())
621                 std::cout << "WARNING: Found only " << jet.genparticles_indices().size() << " from " << jetgenps.size() << " gen particles of this jet"<<std::endl;
622             }
623 +          
624           }
625 +        
626           jets[j].push_back(jet);
627         }
628       }

Diff Legend

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