ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/HbbAnalysis/plugins/HbbTreeMaker.cc
(Generate patch)

Comparing UserCode/HbbAnalysis/plugins/HbbTreeMaker.cc (file contents):
Revision 1.4 by amagnan, Tue Feb 9 14:52:23 2010 UTC vs.
Revision 1.5 by amagnan, Wed Feb 10 13:50:17 2010 UTC

# Line 28 | Line 28
28   #include "UserCode/HbbAnalysis/interface/Jet.hh"
29   #include "UserCode/HbbAnalysis/interface/Met.hh"
30   #include "UserCode/HbbAnalysis/interface/Trigger.hh"
31 + #include "UserCode/HbbAnalysis/interface/Vertex.hh"
32  
33   #include "UserCode/HbbAnalysis/plugins/HbbTreeMaker.hh"
34  
# Line 35 | Line 36 | using namespace HbbAnalysis;
36  
37   HbbTreeMaker::HbbTreeMaker(const edm::ParameterSet & pset):
38    debug_(pset.getParameter<int>("DEBUG")),
39 <  processData_(pset.getParameter<bool>("ProcessData")),
39 >  processData_(false),
40    flavour_(pset.getParameter<unsigned int>("JetFlavour")),
41    doGen_(pset.getParameter<bool>("DOGEN")),
42    genParticleSrc_(pset.getParameter<edm::InputTag>("GenParticles")),
# Line 69 | Line 70 | HbbTreeMaker::~HbbTreeMaker(){//destruct
70  
71  
72  
73 < void HbbTreeMaker::beginJob(const edm::EventSetup&){//beginJob
73 > void HbbTreeMaker::beginJob(){//beginJob
74  
75  
76    edm::Service<TFileService> lFileService;
# Line 83 | Line 84 | void HbbTreeMaker::beginJob(const edm::E
84    if (debug_) std::cout << "Initialising JetFlavour : " << std::endl;
85  
86    lDir = lFileService->mkdir("JetFlavours");
87 <
87 <  if (!processData_) jetFlav_.Initialise(lDir, debug_, flavour_);
87 >  jetFlav_.Initialise(lDir, debug_, flavour_);
88  
89   }//beginJob
90  
# Line 107 | Line 107 | void HbbTreeMaker::analyze(const edm::Ev
107  
108    event_->Clear();
109    event_->event(aEvt.id().event());
110 +  processData_ = aEvt.isRealData();
111 +
112 +  event_->run(aEvt.run());
113 +  event_->lumiBlock(aEvt.luminosityBlock());
114  
115    edm::Handle<reco::GenParticleCollection> lGenParticles;
116    try {
# Line 135 | Line 139 | void HbbTreeMaker::analyze(const edm::Ev
139      std::cout << "AMM: Collection " << vertexSrc_  << " not available! Exception : " << e.what() << ". " << std::endl;
140    }
141  
142 +  HbbVertices(lRecoVertices,event_->vertices());
143 +
144    edm::Handle<std::vector<pat::Electron> > lElectronCollection;
145    
146    try {
# Line 165 | Line 171 | void HbbTreeMaker::analyze(const edm::Ev
171    HbbMuons(lMuonCollection,lRecoVertices,event_->muons());
172  
173  
174 <  edm::Handle<std::vector<pat::Tau> > lTauCollection;
174 >  if (!( caloTauSrc_.label()=="" && caloTauSrc_.instance()=="" )) {
175 >    edm::Handle<std::vector<pat::Tau> > lTauCollection;
176    
177 <  try {
178 <    aEvt.getByLabel(caloTauSrc_,lTauCollection);
179 <    if (!lTauCollection.isValid()){
180 <      edm::LogInfo("ERROR")<< "Error! Can't get caloTau by label. ";
181 <    }
182 <    if (debug_) std::cout << "** caloTaucollection = " << lTauCollection->size() << " elements." << std::endl;
183 <  } catch(cms::Exception& e)  {
184 <    std::cout << "AMM: Collection " << caloTauSrc_  << " not available! Exception : " << e.what() << ". " << std::endl;
185 <  }
177 >    try {
178 >      aEvt.getByLabel(caloTauSrc_,lTauCollection);
179 >      if (!lTauCollection.isValid()){
180 >        edm::LogInfo("ERROR")<< "Error! Can't get caloTau by label. ";
181 >      }
182 >      if (debug_) std::cout << "** caloTaucollection = " << lTauCollection->size() << " elements." << std::endl;
183 >    } catch(cms::Exception& e)  {
184 >      std::cout << "AMM: Collection " << caloTauSrc_  << " not available! Exception : " << e.what() << ". " << std::endl;
185 >    }
186  
187 <  HbbTaus(lTauCollection,lRecoVertices,event_->caloTaus());
187 >    HbbTaus(lTauCollection,lRecoVertices,event_->caloTaus());
188 >  }
189  
190    edm::Handle<std::vector<pat::Tau> > lPFTauCollection;
191    
# Line 687 | Line 695 | void HbbTreeMaker::HbbTaus(const edm::Ha
695            const std::vector<std::pair<std::string,float> > &  lIDs = (*iter).tauIDs();
696            bool lPrint = false;
697            //a few warning if additional discriminants are found:
698 <          if ((isPF && lIDs.size() != 7) ||
698 >          if ((isPF && lIDs.size() != 16) ||
699                (isCalo && lIDs.size() != 3))
700              lPrint = true;
701  
702            if (lPrint) {
703 <            std::cout << "!!!!!!! Discriminants changed, please update histograms !!!!!!!" << std::endl;
703 >            std::cout << "!!!!!!! Discriminants changed, please update tree !!!!!!!" << std::endl;
704              std::cout << "--- isCaloTau = " << isCalo << ", isPFTau = " << isPF << std::endl;
705              std::cout << "------ ID names = " << std::endl;
706            }
# Line 711 | Line 719 | void HbbTreeMaker::HbbTaus(const edm::Ha
719            lpfId.byIsolation = 0;
720            lpfId.againstElectron = 0;
721            lpfId.againstMuon = 0;
722 +          lpfId.byIsolationUsingLeadingPion = 0;
723 +          lpfId.byTaNC = 0;
724 +          lpfId.byTaNCfrHalfPercent = 0;
725 +          lpfId.byTaNCfrOnePercent = 0;
726 +          lpfId.byTaNCfrQuarterPercent = 0;
727 +          lpfId.byTaNCfrTenthPercent = 0;
728 +          lpfId.ecalIsolationUsingLeadingPion = 0;
729 +          lpfId.leadingPionPtCut = 0;
730 +          lpfId.trackIsolationUsingLeadingPion = 0;
731  
732  
733  
# Line 725 | Line 742 | void HbbTreeMaker::HbbTaus(const edm::Ha
742              if (isPF) {
743                if (lName.find("leadingTrackFinding") != lName.npos) lpfId.byLeadingTrackFinding = lDiscri;
744                if (lName.find("leadingTrackPtCut") != lName.npos) lpfId.byLeadingTrackPtCut = lDiscri;
745 <              if (lName.find("trackIsolation") != lName.npos) lpfId.byTrackIsolation = lDiscri;
746 <              if (lName.find("ecalIsolation") != lName.npos) lpfId.byECALIsolation = lDiscri;
747 <              if (lName.find("byIsolation") != lName.npos) lpfId.byIsolation = lDiscri;
745 >              if (lName.find("trackIsolationUsingLeadingPion") != lName.npos) lpfId.trackIsolationUsingLeadingPion = lDiscri;
746 >              if (lName.find("trackIsolation") != lName.npos &&
747 >                  lName.find("trackIsolationUsingLeadingPion") == lName.npos) lpfId.byTrackIsolation = lDiscri;
748 >              if (lName.find("ecalIsolationUsingLeadingPion") != lName.npos) lpfId.ecalIsolationUsingLeadingPion = lDiscri;
749 >              if (lName.find("ecalIsolation") != lName.npos &&
750 >                  lName.find("ecalIsolationUsingLeadingPion") == lName.npos) lpfId.byECALIsolation = lDiscri;
751 >              if (lName.find("byIsolationUsingLeadingPion") != lName.npos) lpfId.byIsolationUsingLeadingPion = lDiscri;
752 >              if (lName.find("byIsolation") != lName.npos &&
753 >                  lName.find("byIsolationUsingLeadingPion") == lName.npos) lpfId.byIsolation = lDiscri;
754                if (lName.find("againstElectron") != lName.npos) lpfId.againstElectron = lDiscri;
755                if (lName.find("againstMuon") != lName.npos) lpfId.againstMuon = lDiscri;
756 +              if (lName.find("byTaNCfrHalfPercent") != lName.npos) lpfId.byTaNCfrHalfPercent = lDiscri;
757 +              if (lName.find("byTaNCfrOnePercent") != lName.npos) lpfId.byTaNCfrOnePercent = lDiscri;
758 +              if (lName.find("byTaNCfrQuarterPercent") != lName.npos) lpfId.byTaNCfrQuarterPercent = lDiscri;
759 +              if (lName.find("byTaNCfrTenthPercent") != lName.npos) lpfId.byTaNCfrTenthPercent = lDiscri;
760 +              if (lName.find("byTaNC") != lName.npos &&
761 +                  lName.find("byTaNCfr") == lName.npos) lpfId.byTaNC = lDiscri;
762 +              if (lName.find("leadingPionPtCut") != lName.npos) lpfId.leadingPionPtCut = lDiscri;
763              }
764              if (isCalo){
765                if (lName.find("byIsolation") != lName.npos) lcaloId.byIsolation = lDiscri;
# Line 925 | Line 955 | void HbbTreeMaker::HbbJets(const edm::Ha
955            lBtag.bProba = (*iter).bDiscriminator("jetBProbabilityBJetTags");
956            lBtag.probability = (*iter).bDiscriminator("jetProbabilityBJetTags");
957            lBtag.sSV = (*iter).bDiscriminator("simpleSecondaryVertexBJetTags");
958 <          lBtag.softElectron = (*iter).bDiscriminator("softElectronBJetTags");
958 >          lBtag.softElectronByPt = (*iter).bDiscriminator("softElectronByPtBJetTags");
959 >          lBtag.softElectronByIP3d = (*iter).bDiscriminator("softElectronByIP3dBJetTags");
960            lBtag.softMuon = (*iter).bDiscriminator("softMuonBJetTags");
961 <          lBtag.softMuonNoIP = (*iter).bDiscriminator("softMuonNoIPBJetTags");
961 >          lBtag.softMuonByPt = (*iter).bDiscriminator("softMuonByPtBJetTags");
962 >          lBtag.softMuonByIP3d = (*iter).bDiscriminator("softMuonByIP3dBJetTags");
963            lBtag.tCHE = (*iter).bDiscriminator("trackCountingHighEffBJetTags");
964            lBtag.tCHP = (*iter).bDiscriminator("trackCountingHighPurBJetTags");
965 <          
965 >
966 >
967            bool isCalo = (*iter).isCaloJet();
968            bool isPF = (*iter).isPFJet();
969  
970            assert (isCalo == !isPF);
971            if (isCalo) {
972 +            HbbAnalysis::JetIDVars lId;
973 +            lId.fHPD = (*iter).jetID().fHPD;
974 +            lId.fRBX = (*iter).jetID().fRBX;
975 +            lId.n90Hits = (*iter).jetID().n90Hits;
976 +            
977              HbbAnalysis::CaloJetVars lCalo;
978              lCalo.maxEInEmTowers = (*iter).maxEInEmTowers();
979              lCalo.maxEInHadTowers = (*iter).maxEInHadTowers();
# Line 952 | Line 990 | void HbbTreeMaker::HbbJets(const edm::Ha
990              lCalo.n90 = (*iter).n90();
991              lCalo.n60 = (*iter).n60();
992  
993 <            HbbAnalysis::Jet lObj(lGen,lGenJet,lReco,lCommon,lCalo,lBtag);
993 >            HbbAnalysis::Jet lObj(lGen,lGenJet,lReco,lCommon,lCalo,lBtag,lId);
994              aVec.push_back(lObj);
995            
996            }
# Line 978 | Line 1016 | void HbbTreeMaker::HbbJets(const edm::Ha
1016            
1017            }
1018  
1019 +
1020 +
1021            iEle++;
1022          }//loop on element
1023      }//non empty
# Line 1115 | Line 1155 | void HbbTreeMaker::HbbParticles(const ed
1155  
1156  
1157   }//genparticles
1158 +
1159 +
1160 +
1161 +
1162 +
1163 + void HbbTreeMaker::HbbVertices(const edm::Handle<std::vector<reco::Vertex> > & aCol,
1164 +                               std::vector<HbbAnalysis::Vertex> & aVec)
1165 + {
1166 +
1167 +  for (std::vector<reco::Vertex>::const_iterator iter = aCol->begin();
1168 +       iter != aCol->end();
1169 +       iter++)
1170 +    {
1171 +
1172 +      if (!((*iter).isValid()) || (*iter).isFake()) continue;
1173 +
1174 +      HbbAnalysis::VertexVars lVtx;
1175 +
1176 +      if ((*iter).tracksSize() > 0) {
1177 +        for (std::vector<reco::TrackBaseRef>::const_iterator lTrk = (*iter).tracks_begin();
1178 +             lTrk != (*iter).tracks_end();
1179 +             lTrk++) {
1180 +          lVtx.trackWeights.push_back((*iter).trackWeight(*lTrk));
1181 +        }
1182 +      }
1183 +
1184 +      if (lVtx.trackWeights.size() != (*iter).tracksSize()) {
1185 +        std::cout<< " -- Problem with tracks, size is not as expected ! "
1186 +                 << std::endl
1187 +                 << " --- Size of recoVertex tracks = " << (*iter).tracksSize()
1188 +                 << std::endl
1189 +                 << " --- Size of trackWeights = " << lVtx.trackWeights.size()
1190 +                 << std::endl;
1191 +      }
1192 +
1193 +      lVtx.chi2 = (*iter).chi2();
1194 +      lVtx.ndof = (*iter).ndof();
1195 +      lVtx.x = (*iter).x();
1196 +      lVtx.y = (*iter).y();
1197 +      lVtx.z = (*iter).z();
1198 +      lVtx.xError = (*iter).xError();
1199 +      lVtx.yError = (*iter).yError();
1200 +      lVtx.zError = (*iter).zError();
1201 +      lVtx.cov01 = (*iter).covariance(0,1);
1202 +      lVtx.cov02 = (*iter).covariance(0,2);
1203 +      lVtx.cov12 = (*iter).covariance(1,2);
1204 +
1205 +      HbbAnalysis::Vertex lObj(lVtx);
1206 +      aVec.push_back(lObj);
1207 +    }
1208 +
1209 + }

Diff Legend

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