257 bool passed_electron_HLT =
true;
258 bool passed_muon_HLT =
true;
267 int vertex_number = vertexCollection.size();
268 VertexCollection::const_iterator
v = vertexCollection.begin();
269 double vertex_chi2 = v->normalizedChi2();
270 double vertex_d0 =
sqrt(v->x() * v->x() + v->y() * v->y());
272 double vertex_numTrks = v->tracksSize();
273 double vertex_sumTrks = 0.0;
275 vertex_curTrack != v->tracks_end(); vertex_curTrack++) {
276 vertex_sumTrks += (*vertex_curTrack)->pt();
291 if (electronCollection.
isValid()) {
292 int posEle = 0, negEle = 0;
295 if (passed_electron_HLT) {
296 for (reco::GsfElectronCollection::const_iterator recoElectron =
297 electronCollection->begin();
298 recoElectron != electronCollection->end(); recoElectron++) {
304 if (recoElectron->charge() == 1) {
306 }
else if (recoElectron->charge() == -1) {
318 nEle = posEle + negEle;
319 if (nEle > 9.) nEle = 9.;
323 unsigned int eleCollectionSize = electronCollection->size();
324 for (
unsigned int i = 0;
i < eleCollectionSize;
i++) {
326 double pt = ele.
pt();
328 for (
unsigned int j =
i + 1; j < eleCollectionSize; j++) {
329 const GsfElectron& ele2 = electronCollection->at(j);
330 double pt2 = ele2.
pt();
333 ele.
px() + ele2.
px(), ele.
py() + ele2.
py(),
334 ele.
pz() + ele2.
pz(), ele.
p() + ele2.
p());
348 if (muonCollection.
isValid()) {
350 int posMu = 0, negMu = 0;
351 TLorentzVector m1, m2;
352 if (passed_muon_HLT) {
353 for (reco::MuonCollection::const_iterator
recoMuon =
354 muonCollection->begin();
362 }
else if (
recoMuon->isGlobalMuon() &&
368 }
else if (
recoMuon->isGlobalMuon()) {
372 }
else if (
recoMuon->isTrackerMuon() &&
378 }
else if (
recoMuon->isStandAloneMuon()) {
385 }
else if (
recoMuon->charge() == -1) {
390 if (nMu > 9.) nMu = 9.;
395 unsigned int muonCollectionSize = muonCollection->size();
396 for (
unsigned int i = 0;
i < muonCollectionSize;
i++) {
397 const Muon&
mu = muonCollection->at(
i);
401 for (
unsigned int j =
i + 1; j < muonCollectionSize; j++) {
402 const Muon& mu2 = muonCollection->at(j);
403 double pt2 = mu2.
pt();
408 mu.
px() + mu2.
px(), mu.
py() + mu2.
py(), mu.
pz() + mu2.
pz(),
415 mu.
px() + mu2.
px(), mu.
py() + mu2.
py(), mu.
pz() + mu2.
pz(),
422 mu.
px() + mu2.
px(), mu.
py() + mu2.
py(), mu.
pz() + mu2.
pz(),
437 if (caloJetCollection.
isValid()) {
442 float jet2_et = -9.0;
445 for (CaloJetCollection::const_iterator i_calojet =
446 caloJetCollection->begin();
447 i_calojet != caloJetCollection->end(); i_calojet++) {
448 float jet_current_et = i_calojet->et();
455 if (jet_current_et < 15)
continue;
457 if (jet_current_et > jet_et) {
461 jet_et = i_calojet->et();
465 }
else if (jet_current_et > jet2_et) {
466 jet2_et = i_calojet->et();
482 if (caloMETCollection.
isValid()) {
483 float caloMet = caloMETCollection->begin()->et();
484 float caloMet_phi = caloMETCollection->begin()->phi();
490 if (pfMETCollection.
isValid()) {
491 float pfMet = pfMETCollection->begin()->et();
492 float pfMet_phi = pfMETCollection->begin()->phi();
500 if (nMu + nEle > 2 && nMu + nEle < 10) {
551 if ((nMu + nEle) >= 10)
552 LogDebug(
"HiggsDQM") <<
"WARNING: " << nMu + nEle
553 <<
" leptons in this event: run=" <<
e.id().run()
554 <<
", event=" <<
e.id().event() <<
"\n";
edm::EDGetTokenT< reco::CaloMETCollection > theCaloMETCollectionToken_
edm::EDGetTokenT< reco::GsfElectronCollection > theElectronCollectionToken_
MonitorElement * h_mPhi_TM
MonitorElement * h_mPt_TM
MonitorElement * h_vertex_numTrks
MonitorElement * h_mEta_GMTM
double px() const final
x coordinate of momentum vector
MonitorElement * h_mPhi_GM
MonitorElement * h_mEta_TM
edm::EDGetTokenT< reco::VertexCollection > vertexToken_
double pt() const final
transverse momentum
std::vector< Vertex > VertexCollection
collection of Vertex objects
MonitorElement * h_mPt_STAM
MonitorElement * h_vertex_d0
bool isTrackerMuon() const override
MonitorElement * h_dimumass_GMGM
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
MonitorElement * h_jet_count
MonitorElement * h_mPt_GMPT
bool isGlobalMuon() const override
double pz() const final
z coordinate of momentum vector
MonitorElement * h_pfMet_phi
pfMet
____________________________________________________________________________||
float segmentCompatibility(const reco::Muon &muon, reco::Muon::ArbitrationType arbitrationType=reco::Muon::SegmentAndTrackArbitration)
MonitorElement * h_mPhi_GMPT
edm::EDGetTokenT< reco::CaloJetCollection > theCaloJetCollectionToken_
caloMet
____________________________________________________________________________||
MonitorElement * h_jet_et
bool isGoodMuon(const reco::Muon &muon, SelectionType type, reco::Muon::ArbitrationType arbitrationType=reco::Muon::SegmentAndTrackArbitration)
main GoodMuon wrapper call
double p() const final
magnitude of momentum vector
T const * product() const
MonitorElement * h_eMultiplicity
MonitorElement * h_vertex_chi2
MonitorElement * h_vertex_number
MonitorElement * h_dimumass_GMTM
double py() const final
y coordinate of momentum vector
MonitorElement * h_mEta_GMPT
MonitorElement * h_vertex_sumTrks
MonitorElement * h_caloMet
MonitorElement * h_mPt_GMTM
MonitorElement * h_caloMet_phi
MonitorElement * h_mPhi_STAM
edm::EDGetTokenT< reco::PFMETCollection > thePfMETCollectionToken_
std::vector< TrackBaseRef >::const_iterator trackRef_iterator
The iteratator for the vector<TrackRef>
MonitorElement * h_mEta_STAM
MonitorElement * h_mEta_GM
MonitorElement * h_dielemass
MonitorElement * h_mPt_GM
MonitorElement * h_mPhi_GMTM
MonitorElement * h_dimumass_TMTM
edm::EDGetTokenT< reco::MuonCollection > theMuonCollectionToken_
MonitorElement * h_lepcounts
MonitorElement * h_mMultiplicity
Protected Member Functions inherited from edm::EDConsumerBase
1.8.11
