|
|
|
Definition at line 50 of file ZMuMuAnalyzer_cynematics.cc.
| ZMuMuAnalyzer_cynematics::ZMuMuAnalyzer_cynematics | ( | const edm::ParameterSet & | pset | ) |
Definition at line 87 of file ZMuMuAnalyzer_cynematics.cc.
References a, b, binEta, binPhi, binPt, etacut_, etamax, etamin, h_OneSta_mass, h_zMuMu_mass_, h_zMuMu_numberOfCand, h_zMuMu_numberOfCand_etapassed, h_zMuMu_numberOfCand_isopassed, h_zMuMu_numberOfCand_masspassed, h_zMuMu_numberOfCand_passed, h_zMuMu_numberOfCand_ptetamasspassed, h_zMuMu_numberOfCand_ptetapassed, h_zMuMu_numberOfCand_ptpassed, h_zMuMuMatched_mass_, h_zMuSingleStandAlone_mass_, h_zMuSingleStandAloneMatched_mass_, h_zMuSingleStandAloneOverlap_mass_, h_zMuSingleStandAloneOverlapMatched_mass_, h_zMuSingleTrack_mass_, h_zMuSingleTrackMatched_mass_, h_zMuSta_MCmatched_numberOfCand_passed, h_zMuSta_numberOfCand, h_zMuSta_numberOfCand_notcontained, h_zMuSta_numberOfCand_passed, h_zMuTrack_MCmatched_numberOfCand_passed, h_zMuTrack_numberOfCand, h_zMuTrack_numberOfCand_notcontained, h_zMuTrack_numberOfCand_passed, hmumu_eta, hmumu_phi, hmumu_pt, hmusta_eta, hmusta_phi, hmusta_pt, hmutrack_eta, hmutrack_phi, hmutrack_pt, i, TFileDirectory::make(), numberOfIntervals, phimax, phimin, ptcut_, ptmax, and ptmin.
00087 : 00088 zMuMu_( pset.getParameter<InputTag>( "zMuMu" ) ), 00089 zMuTrack_( pset.getParameter<InputTag>( "zMuTrack" ) ), 00090 zMuStandAlone_( pset.getParameter<InputTag>( "zMuStandAlone" ) ), 00091 muIso_( pset.getParameter<InputTag>( "muIso" ) ), 00092 trackIso_( pset.getParameter<InputTag>( "trackIso" ) ), 00093 standAloneIso_( pset.getParameter<InputTag>( "standAloneIso" ) ), 00094 zMuMuMap_( pset.getParameter<InputTag>( "zMuMuMap" ) ), 00095 zMuTrackMap_( pset.getParameter<InputTag>( "zMuTrackMap" ) ), 00096 zMuStandAloneMap_( pset.getParameter<InputTag>( "zMuStandAloneMap" ) ), 00097 isocut_( pset.getParameter<double>( "isocut" ) ), 00098 etacut_( pset.getParameter<double>( "etacut" ) ), 00099 ptcut_( pset.getParameter<double>( "ptcut" ) ), 00100 ptSTAcut_( pset.getParameter<double>( "ptSTAcut" ) ), 00101 minZmass_( pset.getParameter<double>( "minZmass" )), 00102 maxZmass_( pset.getParameter<double>( "maxZmass" )) { 00103 00104 Service<TFileService> fs; 00105 h_zMuMu_numberOfCand = fs->make<TH1D>("ZMuMunumberOfCand","number of ZMuMu cand",10, -.5, 9.5); 00106 h_zMuMu_numberOfCand_passed = fs->make<TH1D>("ZMuMunumberOfCandpassed","number of ZMuMu cand selected",10, -.5, 9.5); 00107 h_zMuMu_numberOfCand_ptpassed = fs->make<TH1D>("ZMuMunumberOfCandptpassed","number of ZMuMu cand after pt cut selected",10, -.5, 9.5); 00108 h_zMuMu_numberOfCand_etapassed = fs->make<TH1D>("ZMuMunumberOfCandetapassed","number of ZMuMu cand after eta cut selected",10, -.5, 9.5); 00109 h_zMuMu_numberOfCand_masspassed = fs->make<TH1D>("ZMuMunumberOfCandmasspassed","number of ZMuMu cand after mass cut selected",10, -.5, 9.5); 00110 h_zMuMu_numberOfCand_isopassed = fs->make<TH1D>("ZMuMunumberOfCandisopassed","number of ZMuMu cand after iso cut selected",10, -.5, 9.5); 00111 h_zMuMu_numberOfCand_ptetapassed = fs->make<TH1D>("ZMuMunumberOfCandptetapassed","number of ZMuMu cand after pt & eta cut selected",10, -.5, 9.5); 00112 h_zMuMu_numberOfCand_ptetamasspassed = fs->make<TH1D>("ZMuMunumberOfCandptetamaspassed","number of ZMuMu cand after pt & eta & mass cut selected",10, -.5, 9.5); 00113 00114 00115 h_zMuMu_mass_ = fs->make<TH1D>( "ZMuMumass", "ZMuMu mass(GeV)", 200, 0., 200. ); 00116 h_zMuSingleTrack_mass_ = fs->make<TH1D>( "ZMuSingleTrackmass", "ZMuSingleTrack mass(GeV)", 100, 0., 200. ); 00117 h_zMuSingleStandAlone_mass_ = fs->make<TH1D>( "ZMuSingleStandAlonemass", "ZMuSingleStandAlone mass(GeV)", 50, 0., 200. ); 00118 h_zMuSingleStandAloneOverlap_mass_ = fs->make<TH1D>( "ZMuSingleStandAloneOverlapmass", "ZMuSingleStandAloneOverlap mass(GeV)", 50, 0., 200. ); 00119 00120 00121 h_zMuMuMatched_mass_ = fs->make<TH1D>( "ZMuMuMatchedmass", "ZMuMu Matched mass(GeV)", 200, 0., 200. ); 00122 h_zMuSingleTrackMatched_mass_ = fs->make<TH1D>( "ZMuSingleTrackmassMatched", "ZMuSingleTrackMatched mass(GeV)", 100, 0., 200. ); 00123 h_zMuSingleStandAloneMatched_mass_ = fs->make<TH1D>( "ZMuSingleStandAlonemassMatched", "ZMuSingleStandAloneMatched mass(GeV)", 50, 0., 200. ); 00124 h_zMuSingleStandAloneOverlapMatched_mass_ = fs->make<TH1D>( "ZMuSingleStandAloneOverlapmassMatched", "ZMuSingleStandAloneMatched Overlap mass(GeV)", 50, 0., 200. ); 00125 00126 h_zMuSta_numberOfCand = fs->make<TH1D>("ZMuStanumberOfCand","number of ZMuSta cand (if ZMuMu not selected)",10, -.5, 9.5); 00127 h_OneSta_mass = fs->make<TH2D>("ZOneMuStaMass","inv. mass of ZMuSta1 vs ZMuSta2 when one ZMuSta has been found (if ZMuMu not selected)",100, 0., 400, 100, 0., 400.); 00128 h_zMuSta_numberOfCand_notcontained = fs->make<TH1D>("ZMuStanumberOfCandnotcontained","number of independent ZMuSta cand (if ZMuMu not selected)",10, -.5, 9.5); 00129 h_zMuSta_numberOfCand_passed = fs->make<TH1D>("ZMuStanumberOfCandpassed","number of ZMuSta cand selected (if ZMuMu not selected)",10, -.5, 9.5); 00130 h_zMuSta_MCmatched_numberOfCand_passed = fs->make<TH1D>("ZMuStaMCmatchedNumberOfCandpassed","number of ZMuSta MC matched cand selected (if ZMuMu not selected)",10, -.5, 9.5); 00131 h_zMuTrack_numberOfCand = fs->make<TH1D>("ZMuTranumberOfCand","number of ZMuTrack cand (if ZMuMu and ZMuSTa not selected)",10, -.5, 9.5); 00132 h_zMuTrack_numberOfCand_notcontained = fs->make<TH1D>("ZMuTranumberOfCandnotcontaind","number of indeendent ZMuTrack cand (if ZMuMu and ZMuSTa not selected)",10, -.5, 9.5); 00133 h_zMuTrack_numberOfCand_passed = fs->make<TH1D>("ZMuTranumberOfCandpassed","number of ZMuTrack cand selected (if ZMuMu and ZMuSta not selected)",10, -.5, 9.5); 00134 h_zMuTrack_MCmatched_numberOfCand_passed = fs->make<TH1D>("ZMuTraMCmacthedNumberOfCandpassed","number of ZMuTrack MC matched cand selected (if ZMuMu and ZMuSta not selected)",10, -.5, 9.5); 00135 00136 00137 // creating histograms for each Pt, eta, phi interval 00138 00139 etamin = -etacut_; 00140 etamax = etacut_; 00141 phimin = -3.1415; 00142 phimax = 3.1415; 00143 ptmin = ptcut_; 00144 ptmax = 100; 00145 numberOfIntervals = 8; // number of intervals in which to divide cynematic variables 00146 binEta = (etamax - etamin)/numberOfIntervals; 00147 binPhi = (phimax - phimin)/numberOfIntervals; 00148 binPt = (ptmax - ptmin)/numberOfIntervals; 00149 TFileDirectory etaDirectory = fs->mkdir("etaIntervals"); // in this directory will be saved all the histos of different eta intervals 00150 TFileDirectory phiDirectory = fs->mkdir("phiIntervals"); // in this directory will be saved all the histos of different phi intervals 00151 TFileDirectory ptDirectory = fs->mkdir("ptIntervals"); // in this directory will be saved all the histos of different pt intervals 00152 00153 // eta histograms creation 00154 00155 for (int i=0;i<numberOfIntervals;i++) { 00156 double range0 = etamin + i*binEta; 00157 double range1= range0 + binEta; 00158 char a[30], b[50]; 00159 sprintf(a,"zmumu_etaRange%d",i); 00160 sprintf(b,"zmumu mass eta Range %f to %f",range0,range1); 00161 hmumu_eta.push_back(etaDirectory.make<TH1D>(a,b,200,0.,200.)); 00162 char asta[30], bsta[50]; 00163 sprintf(asta,"zmusta_etaRange%d",i); 00164 sprintf(bsta,"zmusta mass eta Range %f to %f",range0,range1); 00165 hmusta_eta.push_back(etaDirectory.make<TH1D>(asta,bsta,50,0.,200.)); 00166 char atk[30], btk[50]; 00167 sprintf(atk,"zmutrack_etaRange%d",i); 00168 sprintf(btk,"zmutrack mass eta Range %f to %f",range0,range1); 00169 hmutrack_eta.push_back(etaDirectory.make<TH1D>(atk,btk,100,0.,200.)); 00170 } 00171 00172 // phi histograms creation 00173 00174 for (int i=0;i<numberOfIntervals;i++) { 00175 double range0 = phimin + i*binPhi; 00176 double range1= range0 + binPhi; 00177 char a[30], b[50]; 00178 sprintf(a,"zmumu_phiRange%d",i); 00179 sprintf(b,"zmumu mass phi Range %f to %f",range0,range1); 00180 hmumu_phi.push_back(phiDirectory.make<TH1D>(a,b,200,0.,200.)); 00181 char asta[30], bsta[50]; 00182 sprintf(asta,"zmusta_phiRange%d",i); 00183 sprintf(bsta,"zmusta mass phi Range %f to %f",range0,range1); 00184 hmusta_phi.push_back(phiDirectory.make<TH1D>(asta,bsta,50,0.,200.)); 00185 char atk[30], btk[50]; 00186 sprintf(atk,"zmutrack_phiRange%d",i); 00187 sprintf(btk,"zmutrack mass phi Range %f to %f",range0,range1); 00188 hmutrack_phi.push_back(phiDirectory.make<TH1D>(atk,btk,100,0.,200.)); 00189 } 00190 00191 // pt histograms creation 00192 00193 for (int i=0;i<numberOfIntervals;i++) { 00194 double range0 = ptmin + i*binPt; 00195 double range1= range0 + binPt; 00196 char a[30], b[50]; 00197 sprintf(a,"zmumu_ptRange%d",i); 00198 sprintf(b,"zmumu mass pt Range %f to %f",range0,range1); 00199 hmumu_pt.push_back(ptDirectory.make<TH1D>(a,b,200,0.,200.)); 00200 char asta[30], bsta[50]; 00201 sprintf(asta,"zmusta_ptRange%d",i); 00202 sprintf(bsta,"zmusta mass pt Range %f to %f",range0,range1); 00203 hmusta_pt.push_back(ptDirectory.make<TH1D>(asta,bsta,50,0.,200.)); 00204 char atk[30], btk[50]; 00205 sprintf(atk,"zmutrack_ptRange%d",i); 00206 sprintf(btk,"zmutrack mass pt Range %f to %f",range0,range1); 00207 hmutrack_pt.push_back(ptDirectory.make<TH1D>(atk,btk,100,0.,200.)); 00208 } 00209 }
| void ZMuMuAnalyzer_cynematics::analyze | ( | const edm::Event & | event, | |
| const edm::EventSetup & | setup | |||
| ) | [private, virtual] |
Implements edm::EDAnalyzer.
Definition at line 211 of file ZMuMuAnalyzer_cynematics.cc.
References binEta, binPhi, binPt, edm::RefToBase< T >::castTo(), GenMuonPlsPt100GeV_cfg::cout, reco::Candidate::daughter(), lat::endl(), etacut_, etamin, edm::Ref< C, T, F >::get(), h_zMuMu_mass_, h_zMuMu_numberOfCand, h_zMuMu_numberOfCand_etapassed, h_zMuMu_numberOfCand_isopassed, h_zMuMu_numberOfCand_masspassed, h_zMuMu_numberOfCand_passed, h_zMuMu_numberOfCand_ptetamasspassed, h_zMuMu_numberOfCand_ptetapassed, h_zMuMu_numberOfCand_ptpassed, h_zMuMuMatched_mass_, h_zMuSingleStandAlone_mass_, h_zMuSingleStandAloneMatched_mass_, h_zMuSingleTrack_mass_, h_zMuSingleTrackMatched_mass_, h_zMuSta_MCmatched_numberOfCand_passed, h_zMuSta_numberOfCand, h_zMuSta_numberOfCand_notcontained, h_zMuSta_numberOfCand_passed, h_zMuTrack_MCmatched_numberOfCand_passed, h_zMuTrack_numberOfCand, h_zMuTrack_numberOfCand_notcontained, h_zMuTrack_numberOfCand_passed, hmumu_eta, hmumu_phi, hmumu_pt, hmusta_eta, hmusta_phi, hmusta_pt, hmutrack_eta, hmutrack_phi, hmutrack_pt, i, edm::Ref< C, T, F >::id(), edm::Handle< T >::id(), edm::ProductID::id(), isContained(), iso1, isocut_, j, edm::Ref< C, T, F >::key(), m, reco::Particle::mass(), reco::Candidate::masterClone(), maxZmass_, minZmass_, muIso_, numberOfIntervals, phimin, ptcut_, ptmin, ptSTAcut_, standAloneIso_, trackIso_, zMuMu_, zMuMuMap_, zMuStandAlone_, zMuStandAloneMap_, zMuTrack_, and zMuTrackMap_.
00211 { 00212 Handle<CandidateCollection> zMuMu; 00213 event.getByLabel(zMuMu_, zMuMu); 00214 Handle<CandidateCollection> zMuTrack; 00215 event.getByLabel( zMuTrack_, zMuTrack ); 00216 Handle<CandidateCollection> zMuStandAlone; 00217 event.getByLabel( zMuStandAlone_, zMuStandAlone ); 00218 00219 size_t nZMuMu = zMuMu->size(); 00220 size_t nZTrackMu = zMuTrack->size(); 00221 size_t nZStandAloneMu = zMuStandAlone->size(); 00222 // static const double zMass = 91.1876; // PDG Z mass 00223 00224 cout << "++++++++++++++++++++++++++" << endl; 00225 cout << "nZMuMu = " << nZMuMu << endl; 00226 cout << "nZTrackMu = " << nZTrackMu << endl; 00227 cout << "nZStandAloneMu = " << nZStandAloneMu << endl; 00228 cout << "++++++++++++++++++++++++++" << endl; 00229 00230 // ZMuMu counters 00231 00232 int ZMuMu_passed = 0; 00233 int ZMuMu_ptcut_counter = 0; 00234 int ZMuMu_etacut_counter = 0; 00235 int ZMuMu_masscut_counter = 0; 00236 int ZMuMu_isocut_counter = 0; 00237 int ZMuMu_ptetacut_counter = 0; 00238 int ZMuMu_ptetamasscut_counter = 0; 00239 int ZMuMu_allcut_counter = 0; 00240 00241 // ZMuTrack counters 00242 00243 int ZMuTrack_passed = 0; 00244 int ZMuTrack_notcontained = 0; 00245 int ZMuTrack_MCmatched_passed = 0; 00246 00247 // ZMuStandalone counters 00248 int ZMuStandalone_notcontained = 0; 00249 int ZMuStandalone_passed = 0; 00250 int ZMuStandalone_MCmatched_passed = 0; 00251 00252 Handle<CandMatchMap> zMuMuMap; 00253 if( nZMuMu > 0 ) { 00254 event.getByLabel(zMuMuMap_, zMuMuMap); 00255 } 00256 00257 Handle<CandMatchMap> zMuTrackMap; 00258 if( nZTrackMu > 0 ) { 00259 event.getByLabel( zMuTrackMap_, zMuTrackMap ); 00260 } 00261 00262 Handle<CandMatchMap> zMuStandAloneMap; 00263 if( nZStandAloneMu > 0 ) { 00264 event.getByLabel( zMuStandAloneMap_, zMuStandAloneMap ); 00265 } 00266 00267 Handle<IsolationCollection> muIso; 00268 event.getByLabel(muIso_, muIso); 00269 ProductID muIsoId = muIso->keyProduct().id(); 00270 Handle<IsolationCollection> trackIso; 00271 event.getByLabel(trackIso_, trackIso); 00272 ProductID trackIsoId = trackIso->keyProduct().id(); 00273 00274 Handle<IsolationCollection> standAloneIso; 00275 event.getByLabel(standAloneIso_, standAloneIso); 00276 ProductID standAloneIsoId = standAloneIso->keyProduct().id(); 00277 00278 if (nZMuMu > 0) { 00279 // double mass = 1000000.; 00280 for( size_t i = 0; i < nZMuMu; i++ ) { 00281 bool ptcutAccept = false; 00282 bool etacutAccept = false; 00283 bool masscutAccept = false; 00284 bool isocutAccept = false; 00285 const Candidate & zmmCand = (*zMuMu)[ i ]; 00286 CandidateRef CandRef(zMuMu,i); 00287 CandidateRef lep1 = zmmCand.daughter( 0 )->masterClone().castTo<CandidateRef>(); 00288 CandidateRef lep2 = zmmCand.daughter( 1 )->masterClone().castTo<CandidateRef>(); 00289 00290 const double iso1 = muIso->value( lep1.key() ); 00291 const double iso2 = muIso->value( lep2.key() ); 00292 00293 double m = zmmCand.mass(); 00294 // check single cuts 00295 00296 if (lep1->pt()>ptcut_ && lep2->pt()>ptcut_) ptcutAccept = true; 00297 if (fabs(lep1->eta())<etacut_ && fabs(lep2->eta())<etacut_) etacutAccept = true; 00298 if (m>minZmass_ && m<maxZmass_) masscutAccept = true; 00299 if (iso1 < isocut_ && iso2 <isocut_) isocutAccept = true; 00300 00301 00302 if (ptcutAccept) ZMuMu_ptcut_counter++; 00303 if (etacutAccept) ZMuMu_etacut_counter++; 00304 if (masscutAccept) ZMuMu_masscut_counter++; 00305 if (isocutAccept) ZMuMu_isocut_counter++; 00306 00307 // check sequencial cuts 00308 00309 if (ptcutAccept && etacutAccept) { 00310 ZMuMu_ptetacut_counter++; 00311 if (masscutAccept) { 00312 ZMuMu_ptetamasscut_counter++; 00313 if (isocutAccept) { 00314 ZMuMu_passed++;} 00315 } 00316 } 00317 00318 if (ptcutAccept && etacutAccept && masscutAccept && isocutAccept) { 00319 ZMuMu_allcut_counter++; 00320 h_zMuMu_mass_->Fill( m ); 00321 00322 // check the cynematics to fill correct histograms 00323 for (int j=0;j<numberOfIntervals;j++) { 00324 bool statusBinEta = false; 00325 bool statusBinPhi = false; 00326 bool statusBinPt = false; 00327 double range0 = etamin + j*binEta; 00328 double range1= range0 + binEta; 00329 double range0phi = phimin + j*binPhi; 00330 double range1phi= range0phi + binPhi; 00331 double range0pt = ptmin + j*binPt; 00332 double range1pt = range0pt + binPt; 00333 // eta histograms 00334 if (lep1->eta()>=range0 && lep1->eta()<range1) 00335 { 00336 hmumu_eta[j]->Fill(m); 00337 statusBinEta = true; 00338 } 00339 if (lep2->eta()>=range0 && lep2->eta()<range1 && !statusBinEta){ 00340 hmumu_eta[j]->Fill(m); // If eta1 is in the same bin of eta2 fill just once 00341 } 00342 // phi histograms 00343 if (lep1->phi()>=range0phi && lep1->phi()<range1phi) 00344 { 00345 hmumu_phi[j]->Fill(m); 00346 statusBinPhi = true; 00347 } 00348 if (lep2->phi()>=range0phi && lep2->phi()<range1phi && !statusBinPhi){ 00349 hmumu_phi[j]->Fill(m); // If phi1 is in the same bin of phi2 fill just once 00350 } 00351 // pt histograms 00352 if (lep1->pt()>=range0pt && lep1->pt()<range1pt) 00353 { 00354 hmumu_pt[j]->Fill(m); 00355 statusBinPt = true; 00356 } 00357 if (lep2->pt()>=range0pt && lep2->pt()<range1pt && !statusBinPt){ 00358 hmumu_pt[j]->Fill(m); // If pt1 is in the same bin of pt2 fill just once 00359 } 00360 } 00361 00362 CandMatchMap::const_iterator m0 = zMuMuMap->find(CandRef); 00363 if( m0 != zMuMuMap->end()) { // the Z is matched to MC thruth 00364 h_zMuMuMatched_mass_->Fill( m ); 00365 } 00366 } 00367 } 00368 } 00369 00370 h_zMuMu_numberOfCand->Fill(nZMuMu); // number of Z cand found per event 00371 h_zMuMu_numberOfCand_passed->Fill(ZMuMu_allcut_counter); // number of Z cand after all cuts found per event 00372 h_zMuMu_numberOfCand_ptpassed->Fill(ZMuMu_ptcut_counter); // number of Z cand afer pt cut found per event 00373 h_zMuMu_numberOfCand_etapassed->Fill(ZMuMu_etacut_counter); // number of Z cand afer eta cut found per event 00374 h_zMuMu_numberOfCand_masspassed->Fill(ZMuMu_masscut_counter); // number of Z cand afer mass cut found per event 00375 h_zMuMu_numberOfCand_isopassed->Fill(ZMuMu_isocut_counter); // number of Z cand afer iso cut found per event 00376 h_zMuMu_numberOfCand_ptetapassed->Fill(ZMuMu_ptetacut_counter); // number of Z cand afer pt&eta cut found per event 00377 h_zMuMu_numberOfCand_ptetamasspassed->Fill(ZMuMu_ptetamasscut_counter); // number of Z cand afer pt&eta&mass cut found per event 00378 00379 00380 //ZmuSingleStandAlone (check MuStandalone if MuMu has not been selected by cuts) 00381 // cout << "ZMuMuanalyzer : n of zMuMu " << nZMuMu << " passed " << ZMuMu_passed << " n. of zStaMu " << nZStandAloneMu << endl; 00382 00383 if (ZMuMu_passed == 0 && nZStandAloneMu>0 ) { 00384 // size_t index = 1000; 00385 for( size_t j = 0; j < nZStandAloneMu; j++ ) { 00386 const Candidate & zsmCand = (*zMuStandAlone)[ j ]; 00387 bool skipZ = false; 00388 for( size_t i = 0; i < nZMuMu; i++ ) { // chek if the ZMuSTandalone is contained in a ZMuMu 00389 const Candidate & zmmCand = (*zMuMu)[ i ]; // if yes .. the event has to be skipped 00390 if (isContained(zmmCand,zsmCand)) skipZ=true; 00391 } 00392 if (!skipZ) { // ZSMuSTandalone not contained in a ZMuMu 00393 ZMuStandalone_notcontained++; 00394 CandidateRef CandRef(zMuStandAlone,j); 00395 CandidateRef lep1 = zsmCand.daughter( 0 )->masterClone().castTo<CandidateRef>(); 00396 CandidateRef lep2 = zsmCand.daughter( 1 )->masterClone().castTo<CandidateRef>(); 00397 00398 ProductID id1 = lep1.id(); 00399 ProductID id2 = lep2.id(); 00400 double iso1 = -1; 00401 double iso2 = -1; 00402 00403 if( id1 == muIsoId ) 00404 iso1 = muIso->value( lep1.key() ); 00405 else if ( id1 == standAloneIsoId ) 00406 iso1 = standAloneIso->value( lep1.key() ); 00407 00408 if( id2 == muIsoId ) 00409 iso2 = muIso->value( lep2.key() ); 00410 else if ( id2 == standAloneIsoId ) 00411 iso2 = standAloneIso->value( lep2.key() ); 00412 00413 double ms = zsmCand.mass(); 00414 if (lep1->pt()>ptSTAcut_ && lep2->pt()>ptSTAcut_ && 00415 fabs(lep1->eta())<etacut_ && fabs(lep2->eta())<etacut_ && 00416 ms>minZmass_ && ms<maxZmass_ && iso1<isocut_ && iso2 <isocut_) { 00417 h_zMuSingleStandAlone_mass_->Fill( ms ); 00418 ZMuStandalone_passed++; 00419 // check the cynematics to fill correct histograms 00420 for (int j=0;j<numberOfIntervals;j++) { 00421 double range0 = etamin + j*binEta; 00422 double range1= range0 + binEta; 00423 double range0phi = phimin + j*binPhi; 00424 double range1phi= range0phi + binPhi; 00425 double range0pt = ptmin + j*binPt; 00426 double range1pt = range0pt + binPt; 00427 00428 // check which muon is a standalone (standalone means that there is a reference missing.) 00429 if ((lep1->get<TrackRef,reco::StandAloneMuonTag>()).isNull()) 00430 { 00431 if (lep1->eta()>=range0 && lep1->eta()<range1) hmusta_eta[j]->Fill(ms); 00432 if (lep1->phi()>=range0phi && lep1->phi()<range1phi) hmusta_phi[j]->Fill(ms); 00433 if (lep1->pt()>=range0pt && lep1->pt()<range1pt) hmusta_pt[j]->Fill(ms); 00434 } 00435 if ((lep2->get<TrackRef,reco::StandAloneMuonTag>()).isNull()) 00436 { 00437 if (lep2->eta()>=range0 && lep2->eta()<range1) hmusta_eta[j]->Fill(ms); 00438 if (lep2->phi()>=range0phi && lep2->phi()<range1phi) hmusta_phi[j]->Fill(ms); 00439 if (lep2->pt()>=range0pt && lep2->pt()<range1pt) hmusta_pt[j]->Fill(ms); 00440 } 00441 00442 } 00443 CandMatchMap::const_iterator m0 = zMuStandAloneMap->find(CandRef); 00444 if( m0 != zMuStandAloneMap->end()) { 00445 ZMuStandalone_MCmatched_passed++; 00446 h_zMuSingleStandAloneMatched_mass_->Fill( ms ); 00447 } 00448 } 00449 } 00450 } 00451 h_zMuSta_numberOfCand->Fill(nZStandAloneMu); // number of ZMuStandalone cand found per event (no higher priority Z selected) 00452 h_zMuSta_numberOfCand_notcontained->Fill(ZMuStandalone_notcontained); 00453 h_zMuSta_numberOfCand_passed->Fill(ZMuStandalone_passed); // number of ZMuSTa cand after all cuts found per event (no higher prioriy Z selected) 00454 h_zMuSta_MCmatched_numberOfCand_passed->Fill(ZMuStandalone_MCmatched_passed); // number of ZMuSTa MC matched cand after all cuts found per event (no higher prioriy Z selected) 00455 00456 } 00457 00458 //ZmuSingleTRack (check MuTrack if MuMu has not been selected) 00459 if (ZMuMu_passed == 0 && ZMuStandalone_passed == 0 && nZTrackMu>0) { 00460 for( size_t j = 0; j < nZTrackMu; j++ ) { 00461 const Candidate & ztmCand = (*zMuTrack)[ j ]; 00462 bool skipZ = false; 00463 for( size_t i = 0; i < nZMuMu; i++ ) { // chek if the ZMuTrack is contained in a ZMuMu 00464 const Candidate & zmmCand = (*zMuMu)[ i ]; // if yes .. the event has to be skipped 00465 if (isContained(zmmCand,ztmCand)) skipZ=true; 00466 } 00467 if (!skipZ) { 00468 ZMuTrack_notcontained++; 00469 CandidateRef CandRef(zMuTrack,j); 00470 CandidateRef lep1 = ztmCand.daughter( 0 )->masterClone().castTo<CandidateRef>(); 00471 CandidateRef lep2 = ztmCand.daughter( 1 )->masterClone().castTo<CandidateRef>(); 00472 00473 ProductID id1 = lep1.id(); 00474 ProductID id2 = lep2.id(); 00475 double iso1 = -1; 00476 double iso2 = -1; 00477 00478 if( id1 == muIsoId ) 00479 iso1 = muIso->value( lep1.key() ); 00480 else if ( id1 == trackIsoId ) 00481 iso1 = trackIso->value( lep1.key() ); 00482 00483 if( id2 == muIsoId ) 00484 iso2 = muIso->value( lep2.key() ); 00485 else if ( id2 == trackIsoId ) 00486 iso2 = trackIso->value( lep2.key() ); 00487 00488 double mt = ztmCand.mass(); 00489 if (lep1->pt()>ptcut_ && lep2->pt()>ptcut_ && 00490 fabs(lep1->eta())<etacut_ && fabs(lep2->eta())<etacut_ && 00491 mt>minZmass_ && mt<maxZmass_ && iso1<isocut_ && iso2 <isocut_) { 00492 h_zMuSingleTrack_mass_->Fill( mt ); 00493 ZMuTrack_passed++; 00494 00495 // check the cynematics to fill correct histograms 00496 for (int j=0;j<numberOfIntervals;j++) { 00497 double range0 = etamin + j*binEta; 00498 double range1= range0 + binEta; 00499 double range0phi = phimin + j*binPhi; 00500 double range1phi= range0phi + binPhi; 00501 double range0pt = ptmin + j*binPt; 00502 double range1pt = range0pt + binPt; 00503 00504 // check which muon is a track only (track only means that there is a reference missing.) 00505 if ((lep1->get<TrackRef,reco::StandAloneMuonTag>()).isNull()) 00506 { 00507 if (lep1->eta()>=range0 && lep1->eta()<range1) hmutrack_eta[j]->Fill(mt); 00508 if (lep1->phi()>=range0phi && lep1->phi()<range1phi) hmutrack_phi[j]->Fill(mt); 00509 if (lep1->pt()>=range0pt && lep1->pt()<range1pt) hmutrack_pt[j]->Fill(mt); 00510 } 00511 if ((lep2->get<TrackRef,reco::StandAloneMuonTag>()).isNull()) 00512 { 00513 if (lep2->eta()>=range0 && lep2->eta()<range1) hmutrack_eta[j]->Fill(mt); 00514 if (lep2->phi()>=range0phi && lep2->phi()<range1phi) hmutrack_phi[j]->Fill(mt); 00515 if (lep2->pt()>=range0pt && lep2->pt()<range1pt) hmutrack_pt[j]->Fill(mt); 00516 } 00517 } 00518 CandMatchMap::const_iterator m0 = zMuTrackMap->find(CandRef); 00519 if( m0 != zMuTrackMap->end()) { 00520 ZMuTrack_MCmatched_passed++; 00521 h_zMuSingleTrackMatched_mass_->Fill( mt ); 00522 } 00523 } 00524 } 00525 } 00526 h_zMuTrack_numberOfCand->Fill(nZTrackMu); // number of ZMuTrack cand found per event (no higher priority Z selected) 00527 h_zMuTrack_numberOfCand_notcontained->Fill(ZMuTrack_notcontained); // number of ZMuTrack cand not cntained in ZMuMu (no higher priority Z selected) 00528 00529 h_zMuTrack_numberOfCand_passed->Fill(ZMuTrack_passed); // number of ZMuTrack cand after all cuts found per event (no higher priority Z selected) 00530 00531 h_zMuTrack_MCmatched_numberOfCand_passed->Fill(ZMuTrack_MCmatched_passed); 00532 00533 } 00534 }
Reimplemented from edm::EDAnalyzer.
Definition at line 636 of file ZMuMuAnalyzer_cynematics.cc.
References GenMuonPlsPt100GeV_cfg::cout, lat::endl(), h_zMuMu_mass_, h_zMuMu_numberOfCand, h_zMuMu_numberOfCand_etapassed, h_zMuMu_numberOfCand_isopassed, h_zMuMu_numberOfCand_masspassed, h_zMuMu_numberOfCand_passed, h_zMuMu_numberOfCand_ptetamasspassed, h_zMuMu_numberOfCand_ptetapassed, h_zMuMu_numberOfCand_ptpassed, h_zMuMuMatched_mass_, h_zMuSingleStandAlone_mass_, h_zMuSingleStandAloneMatched_mass_, h_zMuSingleStandAloneOverlap_mass_, h_zMuSingleStandAloneOverlapMatched_mass_, h_zMuSingleTrack_mass_, h_zMuSingleTrackMatched_mass_, h_zMuSta_MCmatched_numberOfCand_passed, h_zMuSta_numberOfCand, h_zMuSta_numberOfCand_notcontained, h_zMuSta_numberOfCand_passed, h_zMuTrack_MCmatched_numberOfCand_passed, h_zMuTrack_numberOfCand, h_zMuTrack_numberOfCand_notcontained, and h_zMuTrack_numberOfCand_passed.
00636 { 00637 00638 // candidate analysis 00639 // ZMuMu 00640 double Nzmmc = h_zMuMu_numberOfCand->GetEntries(); 00641 double Nzmmc_0Z = h_zMuMu_numberOfCand->GetBinContent(1); 00642 double Nzmmc_1Z = h_zMuMu_numberOfCand->GetBinContent(2); 00643 double Nzmmc_moreZ = Nzmmc-Nzmmc_0Z-Nzmmc_1Z; 00644 double Nzmmc_passed_0Z = h_zMuMu_numberOfCand_passed->GetBinContent(1); 00645 double Nzmmc_passed_1Z = h_zMuMu_numberOfCand_passed->GetBinContent(2); 00646 double Nzmmc_passed_moreZ = Nzmmc-Nzmmc_passed_0Z-Nzmmc_passed_1Z; 00647 double Nzmmc_ptpassed_0Z = h_zMuMu_numberOfCand_ptpassed->GetBinContent(1); 00648 double Nzmmc_ptpassed_1Z = h_zMuMu_numberOfCand_ptpassed->GetBinContent(2); 00649 double Nzmmc_etapassed_0Z = h_zMuMu_numberOfCand_etapassed->GetBinContent(1); 00650 double Nzmmc_etapassed_1Z = h_zMuMu_numberOfCand_etapassed->GetBinContent(2); 00651 double Nzmmc_masspassed_0Z = h_zMuMu_numberOfCand_masspassed->GetBinContent(1); 00652 double Nzmmc_masspassed_1Z = h_zMuMu_numberOfCand_masspassed->GetBinContent(2); 00653 double Nzmmc_isopassed_0Z = h_zMuMu_numberOfCand_isopassed->GetBinContent(1); 00654 double Nzmmc_isopassed_1Z = h_zMuMu_numberOfCand_isopassed->GetBinContent(2); 00655 double Nzmmc_ptetapassed_0Z = h_zMuMu_numberOfCand_ptetapassed->GetBinContent(1); 00656 double Nzmmc_ptetapassed_1Z = h_zMuMu_numberOfCand_ptetapassed->GetBinContent(2); 00657 double Nzmmc_ptetamasspassed_0Z = h_zMuMu_numberOfCand_ptetamasspassed->GetBinContent(1); 00658 double Nzmmc_ptetamasspassed_1Z = h_zMuMu_numberOfCand_ptetamasspassed->GetBinContent(2); 00659 double Nzmmc_ptpassed_moreZ = Nzmmc-Nzmmc_ptpassed_0Z-Nzmmc_ptpassed_1Z; 00660 double Nzmmc_etapassed_moreZ = Nzmmc-Nzmmc_etapassed_0Z-Nzmmc_etapassed_1Z; 00661 double Nzmmc_masspassed_moreZ = Nzmmc-Nzmmc_masspassed_0Z-Nzmmc_masspassed_1Z; 00662 double Nzmmc_isopassed_moreZ = Nzmmc-Nzmmc_isopassed_0Z-Nzmmc_isopassed_1Z; 00663 double Nzmmc_ptetapassed_moreZ = Nzmmc-Nzmmc_ptetapassed_0Z-Nzmmc_ptetapassed_1Z; 00664 double Nzmmc_ptetamasspassed_moreZ = Nzmmc-Nzmmc_ptetamasspassed_0Z-Nzmmc_ptetamasspassed_1Z; 00665 double Nzmsc = h_zMuSta_numberOfCand->GetEntries(); 00666 double Nzmsc_0Z = h_zMuSta_numberOfCand->GetBinContent(1); 00667 double Nzmsc_1Z = h_zMuSta_numberOfCand->GetBinContent(2); 00668 double Nzmsc_moreZ = Nzmsc - Nzmsc_0Z - Nzmsc_1Z; 00669 double Nzmsc_notcontained_0Z = h_zMuSta_numberOfCand_notcontained->GetBinContent(1); 00670 double Nzmsc_notcontained_1Z = h_zMuSta_numberOfCand_notcontained->GetBinContent(2); 00671 double Nzmsc_notcontained_moreZ = Nzmsc-Nzmsc_notcontained_0Z-Nzmsc_notcontained_1Z; 00672 double Nzmsc_passed_0Z = h_zMuSta_numberOfCand_passed->GetBinContent(1); 00673 double Nzmsc_passed_1Z = h_zMuSta_numberOfCand_passed->GetBinContent(2); 00674 double Nzmsc_passed_moreZ = Nzmsc - Nzmsc_passed_0Z - Nzmsc_passed_1Z; 00675 double Nzmsc_MCmatched_passed_0Z = h_zMuSta_MCmatched_numberOfCand_passed->GetBinContent(1); 00676 double Nzmsc_MCmatched_passed_1Z = h_zMuSta_MCmatched_numberOfCand_passed->GetBinContent(2); 00677 double Nzmsc_MCmatched_passed_moreZ = Nzmsc - Nzmsc_MCmatched_passed_0Z - Nzmsc_MCmatched_passed_1Z; 00678 double Nzmtc = h_zMuTrack_numberOfCand->GetEntries(); 00679 double Nzmtc_0Z = h_zMuTrack_numberOfCand->GetBinContent(1); 00680 double Nzmtc_1Z = h_zMuTrack_numberOfCand->GetBinContent(2); 00681 double Nzmtc_moreZ = Nzmtc - Nzmtc_0Z - Nzmtc_1Z; 00682 double Nzmtc_notcontained_0Z = h_zMuTrack_numberOfCand_notcontained->GetBinContent(1); 00683 double Nzmtc_notcontained_1Z = h_zMuTrack_numberOfCand_notcontained->GetBinContent(2); 00684 double Nzmtc_notcontained_moreZ = Nzmtc-Nzmtc_notcontained_0Z-Nzmtc_notcontained_1Z; 00685 double Nzmtc_passed_0Z = h_zMuTrack_numberOfCand_passed->GetBinContent(1); 00686 double Nzmtc_passed_1Z = h_zMuTrack_numberOfCand_passed->GetBinContent(2); 00687 double Nzmtc_passed_moreZ = Nzmtc - Nzmtc_passed_0Z - Nzmtc_passed_1Z; 00688 double Nzmtc_MCmatched_passed_0Z = h_zMuTrack_MCmatched_numberOfCand_passed->GetBinContent(1); 00689 double Nzmtc_MCmatched_passed_1Z = h_zMuTrack_MCmatched_numberOfCand_passed->GetBinContent(2); 00690 double Nzmtc_MCmatched_passed_moreZ = Nzmtc - Nzmtc_MCmatched_passed_0Z - Nzmtc_MCmatched_passed_1Z; 00691 00692 cout << "--------------- Statistics ----------------------------------------------------------" << endl; 00693 cout << "n of ZMuMu entries ...................................................... " << Nzmmc << endl; 00694 cout << "n of ZMuMu events with 0 cand ............................................. " << Nzmmc_0Z << endl; 00695 cout << "n of ZMuMu events with 1 cand ............................................. " << Nzmmc_1Z << endl; 00696 cout << "n of ZMuMu events with 2 or more cand ..................................... " << Nzmmc_moreZ << endl << endl ; 00697 00698 cout << "n of ZMuMu events not selected by cuts .................................... " << Nzmmc_passed_0Z << endl; 00699 cout << "n of ZMuMu events with 1 cand selected by cuts ............................ " << Nzmmc_passed_1Z << endl; 00700 cout << "n of ZMuMu events with 2 or more cand elected by cuts ..................... " << Nzmmc_passed_moreZ << endl<< endl ; 00701 00702 cout << "n of ZMuMu events not selected by pt cut .................................. " << Nzmmc_ptpassed_0Z << endl; 00703 cout << "n of ZMuMu events with 1 cand selected by pt cut .......................... " << Nzmmc_ptpassed_1Z << endl; 00704 cout << "n of ZMuMu events with 2 or more cand elected by pt cut ................... " << Nzmmc_ptpassed_moreZ << endl<< endl ; 00705 00706 cout << "n of ZMuMu events not selected by eta cut ................................. " << Nzmmc_etapassed_0Z << endl; 00707 cout << "n of ZMuMu events with 1 cand selected by eta cut ......................... " << Nzmmc_etapassed_1Z << endl; 00708 cout << "n of ZMuMu events with 2 or more cand elected by eta cut .................. " << Nzmmc_etapassed_moreZ << endl<< endl ; 00709 00710 cout << "n of ZMuMu events not selected by mass cut ................................ " << Nzmmc_masspassed_0Z << endl; 00711 cout << "n of ZMuMu events with 1 cand selected by mass cut ........................ " << Nzmmc_masspassed_1Z << endl; 00712 cout << "n of ZMuMu events with 2 or more cand elected by mass cut ................. " << Nzmmc_masspassed_moreZ << endl<< endl ; 00713 00714 cout << "n of ZMuMu events not selected by iso cut ................................. " << Nzmmc_isopassed_0Z << endl; 00715 cout << "n of ZMuMu events with 1 cand selected iso cut ............................ " << Nzmmc_isopassed_1Z << endl; 00716 cout << "n of ZMuMu events with 2 or more cand elected iso cut ..................... " << Nzmmc_isopassed_moreZ << endl<< endl ; 00717 00718 cout << "n of ZMuMu events not selected by pt and eta cut .......................... " << Nzmmc_ptetapassed_0Z << endl; 00719 cout << "n of ZMuMu events with 1 cand selected by pt and eta cut .................. " << Nzmmc_ptetapassed_1Z << endl; 00720 cout << "n of ZMuMu events with 2 or more cand elected by pt and eta cut ........... " << Nzmmc_ptetapassed_moreZ << endl<< endl ; 00721 00722 cout << "n of ZMuMu events not selected by pt and eta and mass cut ................. " << Nzmmc_ptetamasspassed_0Z << endl; 00723 cout << "n of ZMuMu events with 1 cand selected by pt and eta and mass cut ......... " << Nzmmc_ptetamasspassed_1Z << endl; 00724 cout << "n of ZMuMu events with 2 or more cand elected by pt and eta and mass cut .. " << Nzmmc_ptetamasspassed_moreZ << endl<< endl ; 00725 00726 cout << "................When No ZMuMu are selected.................................." << endl; 00727 cout << "n of ZMuSta entries ....................................................... " << Nzmsc << endl; 00728 cout << "n of ZMuSta events with 0 cand ............................................ " << Nzmsc_0Z << endl; 00729 cout << "n of ZMuSta events with 1 cand ............................................ " << Nzmsc_1Z << endl; 00730 cout << "n of ZMuSta events with 2 or more cand .................................... " << Nzmsc_moreZ << endl<< endl ; 00731 00732 cout << "n of ZMuSta not contained events with 0 cand .............................. " << Nzmsc_notcontained_0Z << endl; 00733 cout << "n of ZMuSta events not contained with 1 cand .............................. " << Nzmsc_notcontained_1Z << endl; 00734 cout << "n of ZMuSta events no contained with 2 or more cand ....................... " << Nzmsc_notcontained_moreZ << endl<< endl ; 00735 00736 cout << "n of ZMuSta cand not selectd by cuts ...................................... " << Nzmsc_passed_0Z << endl; 00737 cout << "n of ZMuSta events with 1 cand selected by cuts ........................... " << Nzmsc_passed_1Z << endl; 00738 cout << "n of ZMuSta events with 2 or more cand selected by cuts ................... " << Nzmsc_passed_moreZ << endl<< endl ; 00739 00740 cout << "n of ZMuSta MCmatched cand not selectd by cuts ............................ " << Nzmsc_MCmatched_passed_0Z << endl; 00741 cout << "n of ZMuSta MCmatched events with 1 cand selected by cuts ................. " << Nzmsc_MCmatched_passed_1Z << endl; 00742 cout << "n of ZMuSta MCmatched events with 2 or more cand selected by cuts ......... " << Nzmsc_MCmatched_passed_moreZ << endl<< endl ; 00743 00744 cout << "...............When no ZMuMu and ZMuSta are selcted........................." << endl; 00745 cout << "n of ZMuTrack entries ................................................... " << Nzmtc << endl; 00746 cout << "n of ZMuTrack events with 0 cand ..........................................." << Nzmtc_0Z << endl; 00747 cout << "n of ZMuTrack events with 1 cand ..........................................." << Nzmtc_1Z << endl; 00748 cout << "n of ZMuTrack events with 2 or more cand ..................................." << Nzmtc_moreZ << endl<< endl ; 00749 00750 cout << "n of ZMuTrack not contained events with 0 cand ............................ " << Nzmtc_notcontained_0Z << endl; 00751 cout << "n of ZMuTrack events not contained with 1 cand ............................ " << Nzmtc_notcontained_1Z << endl; 00752 cout << "n of ZMuTrack events no contained with 2 or more cand ..................... " << Nzmtc_notcontained_moreZ << endl<< endl ; 00753 00754 cout << "n of ZMuTrack cand not selectd by cuts ....................................." << Nzmtc_passed_0Z << endl; 00755 cout << "n of ZMuTrack events with 1 cand selected by cuts .........................." << Nzmtc_passed_1Z << endl; 00756 cout << "n of ZMuTrack events with 2 or more cand selected by cuts .................." << Nzmtc_passed_moreZ << endl<< endl ; 00757 00758 cout << "n of ZMuTrack MCmatched cand not selectd by cuts .......................... " << Nzmtc_MCmatched_passed_0Z << endl; 00759 cout << "n of ZMuTrcak MCmatched events with 1 cand selected by cuts ............... " << Nzmtc_MCmatched_passed_1Z << endl; 00760 cout << "n of ZMuTrack MCmatched events with 2 or more cand selected by cuts ....... " << Nzmtc_MCmatched_passed_moreZ << endl; 00761 00762 cout << "------------------------------------------------------------------------------------------" << endl; 00763 00764 double Nzmm = h_zMuMu_mass_->GetEntries() ; 00765 double Nzsm = h_zMuSingleStandAlone_mass_->GetEntries() ; 00766 double Nzsnom = h_zMuSingleStandAloneOverlap_mass_->GetEntries() ; 00767 double Nztm = h_zMuSingleTrack_mass_->GetEntries(); 00768 00769 double NzmmMatch = h_zMuMuMatched_mass_->GetEntries() ; 00770 double NzsmMatch = h_zMuSingleStandAloneMatched_mass_->GetEntries() ; 00771 double NzsnomMatch = h_zMuSingleStandAloneOverlapMatched_mass_->GetEntries() ; 00772 double NztmMatch = h_zMuSingleTrackMatched_mass_->GetEntries(); 00773 00774 cout<<"-- N SingleTrackMu = "<<Nztm<<endl; 00775 cout<<"-----N SinglStandAloneMu = "<<Nzsm<<endl; 00776 cout<<"-----N SingleStandAloneOverlapMu = "<<Nzsnom<<endl; 00777 cout<<"------- N MuMu = "<<Nzmm<<endl; 00778 00779 cout<<"-- N SingleTrackMuMatched = "<<NztmMatch<<endl; 00780 cout<<"-----N SinglStandAloneMuMatched = "<<NzsmMatch<<endl; 00781 cout<<"-----N SingleStandAloneOverlapMuMatched = "<<NzsnomMatch<<endl; 00782 cout<<"------- N MuMu Matched = "<<NzmmMatch<<endl; 00783 }
| bool ZMuMuAnalyzer_cynematics::isContained | ( | const Candidate & | obj1, | |
| const Candidate & | obj2 | |||
| ) | [private] |
Definition at line 536 of file ZMuMuAnalyzer_cynematics.cc.
References GenMuonPlsPt100GeV_cfg::cout, reco::Candidate::daughter(), lat::endl(), reco::Candidate::get(), i, edm::Ref< C, T, F >::id(), j, edm::Ref< C, T, F >::key(), and reco::Candidate::numberOfDaughters().
Referenced by analyze().
00537 { 00538 // check if a candidate obj2 is different from obj1 (assume that obj1 is a ZMuMu and obj2 is any other type) 00539 // (for example a Z can be done with two global muons, or with a global muon plus a standalone muon. 00540 // if the standalone muon is part of the second global muon in fact this is the same Z) 00541 00542 const int maxd = 10; 00543 const Candidate * daughters1[maxd]; 00544 const Candidate * daughters2[maxd]; 00545 TrackRef trackerTrack1[maxd]; 00546 TrackRef stAloneTrack1[maxd]; 00547 TrackRef globalTrack1[maxd]; 00548 TrackRef trackerTrack2[maxd]; 00549 TrackRef stAloneTrack2[maxd]; 00550 TrackRef globalTrack2[maxd]; 00551 bool flag; 00552 size_t nd1 = obj1.numberOfDaughters(); 00553 size_t nd2 = obj2.numberOfDaughters(); 00554 size_t matched=0; 00555 00556 for( size_t i = 0; i < nd1; ++ i ) { 00557 daughters1[i] = obj1.daughter( i ); 00558 trackerTrack1[i] = daughters1[i]->get<TrackRef>(); 00559 stAloneTrack1[i] = daughters1[i]->get<TrackRef,reco::StandAloneMuonTag>(); 00560 globalTrack1[i] = daughters1[i]->get<TrackRef,reco::CombinedMuonTag>(); 00561 00562 /*********************************************** just used for debug ******************** 00563 if (trackerTrack1[i].isNull()) 00564 cout << "in ZMuMu daughter " << i << " tracker ref non found " << endl; 00565 else 00566 cout << "in ZMuMu daughter " << i << " tracker ref FOUND" 00567 << " id: " << trackerTrack1[i].id() << ", index: " << trackerTrack1[i].key() 00568 << endl; 00569 if (stAloneTrack1[i].isNull()) 00570 cout << "in ZMuMu daughter " << i << " stalone ref non found " << endl; 00571 else 00572 cout << "in ZMuMu daughter " << i << " stalone ref FOUND" 00573 << " id: " << stAloneTrack1[i].id() << ", index: " << stAloneTrack1[i].key() 00574 << endl; 00575 00576 if (globalTrack1[i].isNull()) 00577 cout << "in ZMuMu daughter " << i << " global ref non found " << endl; 00578 else 00579 cout << "in ZMuMu daughter " << i << " global ref FOUND" 00580 << " id: " << globalTrack1[i].id() << ", index: " << globalTrack1[i].key() 00581 << endl; 00582 */ 00583 } 00584 for( size_t i = 0; i < nd2; ++ i ) { 00585 daughters2[i] = obj2.daughter( i ); 00586 trackerTrack2[i] = daughters2[i]->get<TrackRef>(); 00587 stAloneTrack2[i] = daughters2[i]->get<TrackRef,reco::StandAloneMuonTag>(); 00588 globalTrack2[i] = daughters2[i]->get<TrackRef,reco::CombinedMuonTag>(); 00589 00590 /******************************************** just used for debug ************ 00591 if (trackerTrack2[i].isNull()) 00592 cout << "in ZMuSta daughter " << i << " tracker ref non found " << endl; 00593 else 00594 cout << "in ZMuSta daughter " << i << " tracker ref FOUND" 00595 << " id: " << trackerTrack2[i].id() << ", index: " << trackerTrack2[i].key() 00596 << endl; 00597 if (stAloneTrack2[i].isNull()) 00598 cout << "in ZMuSta daughter " << i << " standalone ref non found " << endl; 00599 else 00600 cout << "in ZMuSta daughter " << i << " standalone ref FOUND" 00601 << " id: " << stAloneTrack2[i].id() << ", index: " << stAloneTrack2[i].key() 00602 << endl; 00603 00604 if (globalTrack2[i].isNull()) 00605 cout << "in ZMuSta daughter " << i << " global ref non found " << endl; 00606 else 00607 cout << "in ZMuSta daughter " << i << " global ref FOUND" 00608 << " id: " << globalTrack2[i].id() << ", index: " << globalTrack2[i].key() 00609 << endl; 00610 00611 */ 00612 } 00613 if (nd1 != nd2) 00614 { 00615 cout << "ZMuMuAnalyzer::isContained WARNING n.of daughters different " << nd1 << " " << nd2 << endl; 00616 } 00617 else 00618 { 00619 for (size_t i = 0; i < nd1; i++) { 00620 flag = false; 00621 for (size_t j = 0; j < nd2; j++) { // if the obj2 is a standalone the trackref is alwais in the trackerTRack position 00622 if (trackerTrack2[i].id()==trackerTrack1[j].id() && trackerTrack2[i].key()==trackerTrack1[j].key() || 00623 trackerTrack2[i].id()==stAloneTrack1[j].id() && trackerTrack2[i].key()==stAloneTrack1[j].key()) { 00624 flag = true; 00625 } 00626 } 00627 if (flag) matched++; 00628 } 00629 } 00630 if (matched==nd1) // return true if all the childrens of the ZMuMu have a children matched in ZMuXX 00631 return true; 00632 else 00633 return false; 00634 }
double ZMuMuAnalyzer_cynematics::binEta [private] |
Definition at line 80 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), and ZMuMuAnalyzer_cynematics().
double ZMuMuAnalyzer_cynematics::binPhi [private] |
Definition at line 80 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), and ZMuMuAnalyzer_cynematics().
double ZMuMuAnalyzer_cynematics::binPt [private] |
Definition at line 80 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), and ZMuMuAnalyzer_cynematics().
double ZMuMuAnalyzer_cynematics::etacut_ [private] |
Definition at line 64 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), and ZMuMuAnalyzer_cynematics().
double ZMuMuAnalyzer_cynematics::etamax [private] |
Definition at line 78 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by ZMuMuAnalyzer_cynematics().
double ZMuMuAnalyzer_cynematics::etamin [private] |
Definition at line 78 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), and ZMuMuAnalyzer_cynematics().
TH2D* ZMuMuAnalyzer_cynematics::h_OneSta_mass [private] |
Definition at line 76 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by ZMuMuAnalyzer_cynematics().
TH1D * ZMuMuAnalyzer_cynematics::h_zMuMu_mass_ [private] |
Definition at line 65 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), endJob(), and ZMuMuAnalyzer_cynematics().
TH1D* ZMuMuAnalyzer_cynematics::h_zMuMu_numberOfCand [private] |
Definition at line 65 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), endJob(), and ZMuMuAnalyzer_cynematics().
TH1D * ZMuMuAnalyzer_cynematics::h_zMuMu_numberOfCand_etapassed [private] |
Definition at line 65 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), endJob(), and ZMuMuAnalyzer_cynematics().
TH1D * ZMuMuAnalyzer_cynematics::h_zMuMu_numberOfCand_isopassed [private] |
Definition at line 65 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), endJob(), and ZMuMuAnalyzer_cynematics().
TH1D * ZMuMuAnalyzer_cynematics::h_zMuMu_numberOfCand_masspassed [private] |
Definition at line 65 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), endJob(), and ZMuMuAnalyzer_cynematics().
TH1D * ZMuMuAnalyzer_cynematics::h_zMuMu_numberOfCand_passed [private] |
Definition at line 65 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), endJob(), and ZMuMuAnalyzer_cynematics().
TH1D * ZMuMuAnalyzer_cynematics::h_zMuMu_numberOfCand_ptetamasspassed [private] |
Definition at line 65 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), endJob(), and ZMuMuAnalyzer_cynematics().
TH1D * ZMuMuAnalyzer_cynematics::h_zMuMu_numberOfCand_ptetapassed [private] |
Definition at line 65 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), endJob(), and ZMuMuAnalyzer_cynematics().
TH1D * ZMuMuAnalyzer_cynematics::h_zMuMu_numberOfCand_ptpassed [private] |
Definition at line 65 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), endJob(), and ZMuMuAnalyzer_cynematics().
TH1D * ZMuMuAnalyzer_cynematics::h_zMuMuMatched_mass_ [private] |
Definition at line 65 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), endJob(), and ZMuMuAnalyzer_cynematics().
TH1D * ZMuMuAnalyzer_cynematics::h_zMuSingleStandAlone_mass_ [private] |
Definition at line 65 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), endJob(), and ZMuMuAnalyzer_cynematics().
TH1D * ZMuMuAnalyzer_cynematics::h_zMuSingleStandAloneMatched_mass_ [private] |
Definition at line 65 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), endJob(), and ZMuMuAnalyzer_cynematics().
TH1D * ZMuMuAnalyzer_cynematics::h_zMuSingleStandAloneOverlap_mass_ [private] |
Definition at line 65 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by endJob(), and ZMuMuAnalyzer_cynematics().
TH1D * ZMuMuAnalyzer_cynematics::h_zMuSingleStandAloneOverlapMatched_mass_ [private] |
Definition at line 65 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by endJob(), and ZMuMuAnalyzer_cynematics().
TH1D * ZMuMuAnalyzer_cynematics::h_zMuSingleTrack_mass_ [private] |
Definition at line 65 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), endJob(), and ZMuMuAnalyzer_cynematics().
TH1D * ZMuMuAnalyzer_cynematics::h_zMuSingleTrackMatched_mass_ [private] |
Definition at line 65 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), endJob(), and ZMuMuAnalyzer_cynematics().
TH1D * ZMuMuAnalyzer_cynematics::h_zMuSta_MCmatched_numberOfCand_passed [private] |
Definition at line 72 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), endJob(), and ZMuMuAnalyzer_cynematics().
TH1D* ZMuMuAnalyzer_cynematics::h_zMuSta_numberOfCand [private] |
Definition at line 72 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), endJob(), and ZMuMuAnalyzer_cynematics().
TH1D * ZMuMuAnalyzer_cynematics::h_zMuSta_numberOfCand_notcontained [private] |
Definition at line 72 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), endJob(), and ZMuMuAnalyzer_cynematics().
TH1D * ZMuMuAnalyzer_cynematics::h_zMuSta_numberOfCand_passed [private] |
Definition at line 72 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), endJob(), and ZMuMuAnalyzer_cynematics().
TH1D * ZMuMuAnalyzer_cynematics::h_zMuTrack_MCmatched_numberOfCand_passed [private] |
Definition at line 72 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), endJob(), and ZMuMuAnalyzer_cynematics().
TH1D * ZMuMuAnalyzer_cynematics::h_zMuTrack_numberOfCand [private] |
Definition at line 72 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), endJob(), and ZMuMuAnalyzer_cynematics().
TH1D * ZMuMuAnalyzer_cynematics::h_zMuTrack_numberOfCand_notcontained [private] |
Definition at line 72 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), endJob(), and ZMuMuAnalyzer_cynematics().
TH1D * ZMuMuAnalyzer_cynematics::h_zMuTrack_numberOfCand_passed [private] |
Definition at line 72 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), endJob(), and ZMuMuAnalyzer_cynematics().
vector<TH1D *> ZMuMuAnalyzer_cynematics::hmumu_eta [private] |
Definition at line 81 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), and ZMuMuAnalyzer_cynematics().
vector<TH1D *> ZMuMuAnalyzer_cynematics::hmumu_phi [private] |
Definition at line 82 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), and ZMuMuAnalyzer_cynematics().
vector<TH1D *> ZMuMuAnalyzer_cynematics::hmumu_pt [private] |
Definition at line 83 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), and ZMuMuAnalyzer_cynematics().
vector<TH1D *> ZMuMuAnalyzer_cynematics::hmusta_eta [private] |
Definition at line 81 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), and ZMuMuAnalyzer_cynematics().
vector<TH1D *> ZMuMuAnalyzer_cynematics::hmusta_phi [private] |
Definition at line 82 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), and ZMuMuAnalyzer_cynematics().
vector<TH1D *> ZMuMuAnalyzer_cynematics::hmusta_pt [private] |
Definition at line 83 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), and ZMuMuAnalyzer_cynematics().
vector<TH1D *> ZMuMuAnalyzer_cynematics::hmutrack_eta [private] |
Definition at line 81 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), and ZMuMuAnalyzer_cynematics().
vector<TH1D *> ZMuMuAnalyzer_cynematics::hmutrack_phi [private] |
Definition at line 82 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), and ZMuMuAnalyzer_cynematics().
vector<TH1D *> ZMuMuAnalyzer_cynematics::hmutrack_pt [private] |
Definition at line 83 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), and ZMuMuAnalyzer_cynematics().
double ZMuMuAnalyzer_cynematics::isocut_ [private] |
double ZMuMuAnalyzer_cynematics::maxZmass_ [private] |
double ZMuMuAnalyzer_cynematics::minZmass_ [private] |
InputTag ZMuMuAnalyzer_cynematics::muIso_ [private] |
Definition at line 79 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), and ZMuMuAnalyzer_cynematics().
Definition at line 58 of file ZMuMuAnalyzer_cynematics.cc.
double ZMuMuAnalyzer_cynematics::phimax [private] |
Definition at line 78 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by ZMuMuAnalyzer_cynematics().
double ZMuMuAnalyzer_cynematics::phimin [private] |
Definition at line 78 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), and ZMuMuAnalyzer_cynematics().
double ZMuMuAnalyzer_cynematics::ptcut_ [private] |
Definition at line 64 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), and ZMuMuAnalyzer_cynematics().
double ZMuMuAnalyzer_cynematics::ptmax [private] |
Definition at line 78 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by ZMuMuAnalyzer_cynematics().
double ZMuMuAnalyzer_cynematics::ptmin [private] |
Definition at line 78 of file ZMuMuAnalyzer_cynematics.cc.
Referenced by analyze(), and ZMuMuAnalyzer_cynematics().
double ZMuMuAnalyzer_cynematics::ptSTAcut_ [private] |
InputTag ZMuMuAnalyzer_cynematics::trackIso_ [private] |
InputTag ZMuMuAnalyzer_cynematics::zMuMu_ [private] |
InputTag ZMuMuAnalyzer_cynematics::zMuMuMap_ [private] |
InputTag ZMuMuAnalyzer_cynematics::zMuTrack_ [private] |
1.5.4