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.6 by tboccali, Tue Jun 14 12:49:07 2011 UTC vs.
Revision 1.8 by tboccali, Tue Jun 21 15:49:01 2011 UTC

# Line 27 | Line 27 | Implementation:
27   #define GENPTOLORP(a) TLorentzVector((a)->px(), (a)->py(), (a)->pz(), (a)->energy())
28  
29   HbbAnalyzerNew::HbbAnalyzerNew(const edm::ParameterSet& iConfig):
30 <  eleLabel_(iConfig.getUntrackedParameter<edm::InputTag>("electronTag")),
31 <  muoLabel_(iConfig.getUntrackedParameter<edm::InputTag>("muonTag")),
32 <  jetLabel_(iConfig.getUntrackedParameter<edm::InputTag>("jetTag")),
33 <  subjetLabel_(iConfig.getUntrackedParameter<edm::InputTag>("subjetTag")),
34 <  simplejet1Label_(iConfig.getUntrackedParameter<edm::InputTag>("simplejet1Tag")),
35 <  simplejet2Label_(iConfig.getUntrackedParameter<edm::InputTag>("simplejet2Tag")),
36 <  simplejet3Label_(iConfig.getUntrackedParameter<edm::InputTag>("simplejet3Tag")),
37 <  simplejet4Label_(iConfig.getUntrackedParameter<edm::InputTag>("simplejet4Tag")),
38 <  tauLabel_(iConfig.getUntrackedParameter<edm::InputTag>("tauTag")),
39 <  metLabel_(iConfig.getUntrackedParameter<edm::InputTag>("metTag")),
40 <  phoLabel_(iConfig.getUntrackedParameter<edm::InputTag>("photonTag")),
41 <  dimuLabel_(iConfig.getUntrackedParameter<edm::InputTag>("dimuTag")),
42 <  dielecLabel_(iConfig.getUntrackedParameter<edm::InputTag>("dielecTag")),
43 <  hltResults_(iConfig.getUntrackedParameter<edm::InputTag>("hltResultsTag")),
30 >  eleLabel_(iConfig.getParameter<edm::InputTag>("electronTag")),
31 >  muoLabel_(iConfig.getParameter<edm::InputTag>("muonTag")),
32 >  jetLabel_(iConfig.getParameter<edm::InputTag>("jetTag")),
33 >  subjetLabel_(iConfig.getParameter<edm::InputTag>("subjetTag")),
34 >  simplejet1Label_(iConfig.getParameter<edm::InputTag>("simplejet1Tag")),
35 >  simplejet2Label_(iConfig.getParameter<edm::InputTag>("simplejet2Tag")),
36 >  simplejet3Label_(iConfig.getParameter<edm::InputTag>("simplejet3Tag")),
37 >  simplejet4Label_(iConfig.getParameter<edm::InputTag>("simplejet4Tag")),
38 >  tauLabel_(iConfig.getParameter<edm::InputTag>("tauTag")),
39 >  metLabel_(iConfig.getParameter<edm::InputTag>("metTag")),
40 >  phoLabel_(iConfig.getParameter<edm::InputTag>("photonTag")),
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")) {
45  
46    //
# Line 740 | Line 740 | HbbAnalyzerNew::produce(edm::Event& iEve
740        const reco::HitPattern& q = gTrack->hitPattern();
741        mf.globChi2=gTrack.get()->normalizedChi2();
742        mf.globNHits=q.numberOfValidMuonHits();
743 +      mf. validMuStations = q. muonStationsWithValidHits();
744      }else{
745        mf.globChi2=-99;
746        mf.globNHits=-99;
747      }
748 +    //
749 +    // add stamuon
750 +
751 +    //    if (mu->isStandAloneMuon()) {
752 +    //      reco::TrackRef sta = mu->standAloneMuon();
753 +    //      
754 +    //    }
755 +
756  
757      //     int muInfo[12];
758      //     fillMuBlock(mu,  muInfo);
# Line 771 | Line 780 | HbbAnalyzerNew::produce(edm::Event& iEve
780      ef.hIso=elec->hcalIso();
781      Geom::Phi<double> deltaphi(elec->superCluster()->phi()-atan2(hbbInfo->calomet.fourMomentum.Py(),hbbInfo->calomet.fourMomentum.Px()));
782      ef.acop = deltaphi.value();
783 <    
783 >    //
784 >    // fill eleids
785 >    //    
786 >    ef.id95 = elec->electronID("simpleEleId95cIso");
787 >    ef.id85 = elec->electronID("simpleEleId85cIso");
788 >    ef.id70 = elec->electronID("simpleEleId70cIso");
789 >    ef.id95r = elec->electronID("simpleEleId95relIso");
790 >    ef.id70r = elec->electronID("simpleEleId70relIso");
791 >    ef.id85r = elec->electronID("simpleEleId85relIso");
792 >
793      if(runOnMC_){
794        const GenParticle* elecMc = elec->genLepton();
795        if(elecMc!=0){
# Line 872 | Line 890 | HbbAnalyzerNew::produce(edm::Event& iEve
890          const reco::HitPattern& q = gTrack->hitPattern();
891          df.daughter1.globNHits=q.numberOfValidMuonHits();
892          df.daughter1.globChi2=gTrack.get()->normalizedChi2();
893 +        df.daughter1.validMuStations = q. muonStationsWithValidHits();
894        }
895        if(muonDau1.isGlobalMuon()){
896          TrackRef gTrack = muonDau1.globalTrack();
897          const reco::HitPattern& q = gTrack->hitPattern();
898          df.daughter2.globNHits=q.numberOfValidMuonHits();
899          df.daughter2.globChi2=gTrack.get()->normalizedChi2();
900 +        df.daughter2.validMuStations = q. muonStationsWithValidHits();
901        }
902    
903      }
# Line 910 | Line 930 | HbbAnalyzerNew::produce(edm::Event& iEve
930      df.daughter1.hIso = elecDau0.hcalIso();
931      df.daughter2.hIso = elecDau1.hcalIso();
932      
933 <    
933 >    // ids
934 >    df.daughter1.id95 = elecDau0.electronID("simpleEleId95cIso");
935 >    df.daughter1.id85 = elecDau0.electronID  ("simpleEleId85cIso");
936 >    df.daughter1.id70 = elecDau0.electronID  ("simpleEleId70cIso");
937 >    df.daughter1.id95r = elecDau0.electronID ("simpleEleId95relIso");
938 >    df.daughter1.id85r = elecDau0.electronID ("simpleEleId85relIso");
939 >    df.daughter1.id70r = elecDau0.electronID ("simpleEleId70relIso");
940 >
941 >
942 >    df.daughter2.id95 = elecDau1.electronID("simpleEleId95cIso");
943 >    df.daughter2.id85 = elecDau1.electronID  ("simpleEleId85cIso");
944 >    df.daughter2.id70 = elecDau1.electronID  ("simpleEleId70cIso");
945 >    df.daughter2.id95r = elecDau1.electronID ("simpleEleId95relIso");
946 >    df.daughter2.id85r = elecDau1.electronID ("simpleEleId85relIso");
947 >    df.daughter2.id70r = elecDau1.electronID ("simpleEleId70relIso");
948 >
949      hbbInfo->diElectronInfo.push_back(df);
950      
951    }
# Line 961 | Line 996 | TVector2 HbbAnalyzerNew::getTvect( const
996    double r_mag = 1e10;
997    unsigned int nOfconst = 0;
998  
964 –  //  std::cout <<" ECCCCCCOOOOO "<<patJet->isPFJet()<<std::endl;
999  
1000    if (patJet->isPFJet() == false) {
1001      return t_Vect;

Diff Legend

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