CMS 3D CMS Logo
List of all members | Public Member Functions | Private Member Functions | Private Attributes
EgammaSuperClusters Class Reference

#include <EgammaSuperClusters.h>

Inheritance diagram for EgammaSuperClusters:
DQMEDAnalyzer edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >

Public Member Functions

void analyze (const edm::Event &, const edm::EventSetup &) override
 
void bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
 
 EgammaSuperClusters (const edm::ParameterSet &)
 
 ~EgammaSuperClusters () override
 
- Public Member Functions inherited from DQMEDAnalyzer
void accumulate (edm::Event const &event, edm::EventSetup const &setup) final
 
void beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
 
void beginRun (edm::Run const &run, edm::EventSetup const &setup) final
 
void beginStream (edm::StreamID id) final
 
virtual void dqmBeginRun (edm::Run const &, edm::EventSetup const &)
 
 DQMEDAnalyzer ()
 
void endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
 
void endRun (edm::Run const &run, edm::EventSetup const &setup) final
 
virtual bool getCanSaveByLumi ()
 
- Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
 EDProducer ()=default
 
bool hasAbilityToProduceInBeginLumis () const final
 
bool hasAbilityToProduceInBeginProcessBlocks () const final
 
bool hasAbilityToProduceInBeginRuns () const final
 
bool hasAbilityToProduceInEndLumis () const final
 
bool hasAbilityToProduceInEndProcessBlocks () const final
 
bool hasAbilityToProduceInEndRuns () const final
 

Private Member Functions

void closestMCParticle (HepMC::GenEvent const *, reco::SuperCluster const &, double &, double &) const
 
float ecalEta (float, float, float) const
 

Private Attributes

edm::EDGetTokenT< reco::SuperClusterCollection > barrelCorSuperClusterCollectionToken_
 
edm::EDGetTokenT< reco::SuperClusterCollection > barrelRawSuperClusterCollectionToken_
 
edm::EDGetTokenT< EcalRecHitCollection > barrelRecHitCollectionToken_
 
EcalClusterLazyTools::ESGetTokens ecalClusterToolsESGetTokens_
 
edm::EDGetTokenT< reco::SuperClusterCollection > endcapCorSuperClusterCollectionToken_
 
edm::EDGetTokenT< reco::SuperClusterCollection > endcapPreSuperClusterCollectionToken_
 
edm::EDGetTokenT< reco::SuperClusterCollection > endcapRawSuperClusterCollectionToken_
 
edm::EDGetTokenT< EcalRecHitCollection > endcapRecHitCollectionToken_
 
MonitorElement * hist_EB_CorSC_deltaR_
 
MonitorElement * hist_EB_CorSC_EoverTruth_
 
MonitorElement * hist_EB_CorSC_ET_
 
MonitorElement * hist_EB_CorSC_ET_vs_Eta_
 
MonitorElement * hist_EB_CorSC_ET_vs_Phi_
 
MonitorElement * hist_EB_CorSC_Eta_
 
MonitorElement * hist_EB_CorSC_etaWidth_
 
MonitorElement * hist_EB_CorSC_NumBC_
 
MonitorElement * hist_EB_CorSC_Phi_
 
MonitorElement * hist_EB_CorSC_phiWidth_
 
MonitorElement * hist_EB_CorSC_S1toS9_
 
MonitorElement * hist_EB_CorSC_S25toE_
 
MonitorElement * hist_EB_CorSC_Size_
 
MonitorElement * hist_EB_RawSC_deltaR_
 
MonitorElement * hist_EB_RawSC_EoverTruth_
 
MonitorElement * hist_EB_RawSC_ET_
 
MonitorElement * hist_EB_RawSC_Eta_
 
MonitorElement * hist_EB_RawSC_NumBC_
 
MonitorElement * hist_EB_RawSC_Phi_
 
MonitorElement * hist_EB_RawSC_S1toS9_
 
MonitorElement * hist_EB_RawSC_S25toE_
 
MonitorElement * hist_EB_RawSC_Size_
 
MonitorElement * hist_EE_CorSC_deltaR_
 
MonitorElement * hist_EE_CorSC_EoverTruth_
 
MonitorElement * hist_EE_CorSC_ET_
 
MonitorElement * hist_EE_CorSC_ET_vs_Eta_
 
MonitorElement * hist_EE_CorSC_ET_vs_Phi_
 
MonitorElement * hist_EE_CorSC_ET_vs_R_
 
MonitorElement * hist_EE_CorSC_Eta_
 
MonitorElement * hist_EE_CorSC_etaWidth_
 
MonitorElement * hist_EE_CorSC_NumBC_
 
MonitorElement * hist_EE_CorSC_Phi_
 
MonitorElement * hist_EE_CorSC_phiWidth_
 
MonitorElement * hist_EE_CorSC_preshowerE_
 
MonitorElement * hist_EE_CorSC_S1toS9_
 
MonitorElement * hist_EE_CorSC_S25toE_
 
MonitorElement * hist_EE_CorSC_Size_
 
MonitorElement * hist_EE_PreSC_deltaR_
 
MonitorElement * hist_EE_PreSC_EoverTruth_
 
MonitorElement * hist_EE_PreSC_ET_
 
MonitorElement * hist_EE_PreSC_Eta_
 
MonitorElement * hist_EE_PreSC_NumBC_
 
MonitorElement * hist_EE_PreSC_Phi_
 
MonitorElement * hist_EE_PreSC_preshowerE_
 
MonitorElement * hist_EE_PreSC_S1toS9_
 
MonitorElement * hist_EE_PreSC_S25toE_
 
MonitorElement * hist_EE_PreSC_Size_
 
MonitorElement * hist_EE_RawSC_deltaR_
 
MonitorElement * hist_EE_RawSC_EoverTruth_
 
MonitorElement * hist_EE_RawSC_ET_
 
MonitorElement * hist_EE_RawSC_Eta_
 
MonitorElement * hist_EE_RawSC_NumBC_
 
MonitorElement * hist_EE_RawSC_Phi_
 
MonitorElement * hist_EE_RawSC_S1toS9_
 
MonitorElement * hist_EE_RawSC_S25toE_
 
MonitorElement * hist_EE_RawSC_Size_
 
HistSpec hsdeltaR_
 
HistSpec hsEoverTruth_
 
HistSpec hsET_
 
HistSpec hsEta_
 
HistSpec hsetaWidth_
 
HistSpec hsNumBC_
 
HistSpec hsPhi_
 
HistSpec hsphiWidth_
 
HistSpec hspreshowerE_
 
HistSpec hsR_
 
HistSpec hsS1toS9_
 
HistSpec hsS25toE_
 
HistSpec hsSize_
 
edm::EDGetTokenT< edm::HepMCProduct > MCTruthCollectionToken_
 

Additional Inherited Members

- Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore DQMStore
 
typedef dqm::reco::MonitorElement MonitorElement
 
- Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
typedef CacheContexts< T... > CacheTypes
 
typedef CacheTypes::GlobalCache GlobalCache
 
typedef AbilityChecker< T... > HasAbility
 
typedef CacheTypes::LuminosityBlockCache LuminosityBlockCache
 
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache > LuminosityBlockContext
 
typedef CacheTypes::LuminosityBlockSummaryCache LuminosityBlockSummaryCache
 
typedef CacheTypes::RunCache RunCache
 
typedef RunContextT< RunCache, GlobalCache > RunContext
 
typedef CacheTypes::RunSummaryCache RunSummaryCache
 
- Static Public Member Functions inherited from DQMEDAnalyzer
static void globalEndJob (DQMEDAnalyzerGlobalCache const *)
 
static void globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
 
static void globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
 
static std::unique_ptr< DQMEDAnalyzerGlobalCache > initializeGlobalCache (edm::ParameterSet const &)
 
- Protected Member Functions inherited from DQMEDAnalyzer
uint64_t meId () const
 
- Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken > lumiToken_
 
edm::EDPutTokenT< DQMToken > runToken_
 
unsigned int streamId_
 

Detailed Description

Description: SVSuite Super Cluster Validation

Implementation: \

Author
: Michael A. Balazs, Nov 2006

Definition at line 27 of file EgammaSuperClusters.h.

Constructor & Destructor Documentation

◆ EgammaSuperClusters()

EgammaSuperClusters::EgammaSuperClusters ( const edm::ParameterSet &  ps)
explicit

Definition at line 16 of file EgammaSuperClusters.cc.

17  : MCTruthCollectionToken_(consumes<edm::HepMCProduct>(ps.getParameter<edm::InputTag>("MCTruthCollection"))),
18  barrelRawSuperClusterCollectionToken_(
19  consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("barrelRawSuperClusterCollection"))),
20  barrelCorSuperClusterCollectionToken_(
21  consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("barrelCorSuperClusterCollection"))),
22  endcapRawSuperClusterCollectionToken_(
23  consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("endcapRawSuperClusterCollection"))),
24  endcapPreSuperClusterCollectionToken_(
25  consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("endcapPreSuperClusterCollection"))),
26  endcapCorSuperClusterCollectionToken_(
27  consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("endcapCorSuperClusterCollection"))),
28  barrelRecHitCollectionToken_(
29  consumes<EcalRecHitCollection>(ps.getParameter<edm::InputTag>("barrelRecHitCollection"))),
30  endcapRecHitCollectionToken_(
31  consumes<EcalRecHitCollection>(ps.getParameter<edm::InputTag>("endcapRecHitCollection"))),
32  ecalClusterToolsESGetTokens_{consumesCollector()},
33  hsSize_(ps, "Size"),
34  hsNumBC_(ps, "NumBC"),
35  hsET_(ps, "ET"),
36  hsEta_(ps, "Eta"),
37  hsPhi_(ps, "Phi"),
38  hsS1toS9_(ps, "S1toS9"),
39  hsS25toE_(ps, "S25toE"),
40  hsEoverTruth_(ps, "EoverTruth"),
41  hsdeltaR_(ps, "deltaR"),
42  hsphiWidth_(ps, "phiWidth"),
43  hsetaWidth_(ps, "etaWidth"),
44  hspreshowerE_(ps, "preshowerE"),
45  hsR_(ps, "R"),
46  hist_EB_RawSC_Size_(nullptr),
47  hist_EE_RawSC_Size_(nullptr),
48  hist_EB_CorSC_Size_(nullptr),
49  hist_EE_CorSC_Size_(nullptr),
50  hist_EE_PreSC_Size_(nullptr),
51  hist_EB_RawSC_NumBC_(nullptr),
52  hist_EE_RawSC_NumBC_(nullptr),
53  hist_EB_CorSC_NumBC_(nullptr),
54  hist_EE_CorSC_NumBC_(nullptr),
55  hist_EE_PreSC_NumBC_(nullptr),
56  hist_EB_RawSC_ET_(nullptr),
57  hist_EE_RawSC_ET_(nullptr),
58  hist_EB_CorSC_ET_(nullptr),
59  hist_EE_CorSC_ET_(nullptr),
60  hist_EE_PreSC_ET_(nullptr),
61  hist_EB_RawSC_Eta_(nullptr),
62  hist_EE_RawSC_Eta_(nullptr),
63  hist_EB_CorSC_Eta_(nullptr),
64  hist_EE_CorSC_Eta_(nullptr),
65  hist_EE_PreSC_Eta_(nullptr),
66  hist_EB_RawSC_Phi_(nullptr),
67  hist_EE_RawSC_Phi_(nullptr),
68  hist_EB_CorSC_Phi_(nullptr),
69  hist_EE_CorSC_Phi_(nullptr),
70  hist_EE_PreSC_Phi_(nullptr),
71  hist_EB_RawSC_S1toS9_(nullptr),
72  hist_EE_RawSC_S1toS9_(nullptr),
73  hist_EB_CorSC_S1toS9_(nullptr),
74  hist_EE_CorSC_S1toS9_(nullptr),
75  hist_EE_PreSC_S1toS9_(nullptr),
76  hist_EB_RawSC_S25toE_(nullptr),
77  hist_EE_RawSC_S25toE_(nullptr),
78  hist_EB_CorSC_S25toE_(nullptr),
79  hist_EE_CorSC_S25toE_(nullptr),
80  hist_EE_PreSC_S25toE_(nullptr),
81  hist_EB_RawSC_EoverTruth_(nullptr),
82  hist_EE_RawSC_EoverTruth_(nullptr),
83  hist_EB_CorSC_EoverTruth_(nullptr),
84  hist_EE_CorSC_EoverTruth_(nullptr),
85  hist_EE_PreSC_EoverTruth_(nullptr),
86  hist_EB_RawSC_deltaR_(nullptr),
87  hist_EE_RawSC_deltaR_(nullptr),
88  hist_EB_CorSC_deltaR_(nullptr),
89  hist_EE_CorSC_deltaR_(nullptr),
90  hist_EE_PreSC_deltaR_(nullptr),
91  hist_EE_PreSC_preshowerE_(nullptr),
92  hist_EE_CorSC_preshowerE_(nullptr),
93  hist_EE_CorSC_phiWidth_(nullptr),
94  hist_EB_CorSC_phiWidth_(nullptr),
95  hist_EE_CorSC_etaWidth_(nullptr),
96  hist_EB_CorSC_etaWidth_(nullptr),
97  hist_EB_CorSC_ET_vs_Eta_(nullptr),
98  hist_EB_CorSC_ET_vs_Phi_(nullptr),
99  hist_EE_CorSC_ET_vs_Eta_(nullptr),
100  hist_EE_CorSC_ET_vs_Phi_(nullptr),
101  hist_EE_CorSC_ET_vs_R_(nullptr) {}

◆ ~EgammaSuperClusters()

EgammaSuperClusters::~EgammaSuperClusters ( )
override

Definition at line 103 of file EgammaSuperClusters.cc.

103 {}

Member Function Documentation

◆ analyze()

void EgammaSuperClusters::analyze ( const edm::Event &  evt,
const edm::EventSetup &  es 
)
overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 371 of file EgammaSuperClusters.cc.

371  {
372  bool skipMC = false;
373  bool skipBarrel = false;
374  bool skipEndcap = false;
375 
376  //
377  // Get MCTRUTH
378  //
379  edm::Handle<edm::HepMCProduct> pMCTruth;
380  evt.getByToken(MCTruthCollectionToken_, pMCTruth);
381  if (!pMCTruth.isValid()) {
382  edm::LogError("EgammaSuperClusters") << "Error! can't get MC collection ";
383  skipMC = true;
384  }
385  const HepMC::GenEvent *genEvent = pMCTruth->GetEvent();
386 
387  if (skipMC)
388  return;
389 
390  //
391  // Get the BARREL products
392  //
393  edm::Handle<reco::SuperClusterCollection> pBarrelRawSuperClusters;
394  evt.getByToken(barrelRawSuperClusterCollectionToken_, pBarrelRawSuperClusters);
395  if (!pBarrelRawSuperClusters.isValid()) {
396  edm::LogError("EgammaSuperClusters") << "Error! can't get collection Raw SC";
397  skipBarrel = true;
398  }
399 
400  edm::Handle<reco::SuperClusterCollection> pBarrelCorSuperClusters;
401  evt.getByToken(barrelCorSuperClusterCollectionToken_, pBarrelCorSuperClusters);
402  if (!pBarrelCorSuperClusters.isValid()) {
403  edm::LogError("EgammaSuperClusters") << "Error! can't get collection Cor SC";
404  skipBarrel = true;
405  }
406 
407  edm::Handle<EBRecHitCollection> pBarrelRecHitCollection;
408  evt.getByToken(barrelRecHitCollectionToken_, pBarrelRecHitCollection);
409  if (!pBarrelRecHitCollection.isValid()) {
410  skipBarrel = true;
411  }
412  edm::Handle<EERecHitCollection> pEndcapRecHitCollection;
413  evt.getByToken(endcapRecHitCollectionToken_, pEndcapRecHitCollection);
414  if (!pEndcapRecHitCollection.isValid()) {
415  skipEndcap = true;
416  }
417 
418  if (skipBarrel || skipEndcap)
419  return;
420 
421  EcalClusterLazyTools lazyTool(
422  evt, ecalClusterToolsESGetTokens_.get(es), barrelRecHitCollectionToken_, endcapRecHitCollectionToken_);
423 
424  // Get the BARREL collections
425  const reco::SuperClusterCollection *barrelRawSuperClusters = pBarrelRawSuperClusters.product();
426  const reco::SuperClusterCollection *barrelCorSuperClusters = pBarrelCorSuperClusters.product();
427 
428  // Number of entries in collections
429  hist_EB_RawSC_Size_->Fill(barrelRawSuperClusters->size());
430  hist_EB_CorSC_Size_->Fill(barrelCorSuperClusters->size());
431 
432  // Do RAW BARREL SuperClusters
433  for (reco::SuperClusterCollection::const_iterator aClus = barrelRawSuperClusters->begin();
434  aClus != barrelRawSuperClusters->end();
435  aClus++) {
436  // kinematics
437  hist_EB_RawSC_NumBC_->Fill(aClus->clustersSize());
438  hist_EB_RawSC_ET_->Fill(aClus->energy() / std::cosh(aClus->position().eta()));
439  hist_EB_RawSC_Eta_->Fill(aClus->position().eta());
440  hist_EB_RawSC_Phi_->Fill(aClus->position().phi());
441 
442  // cluster shape
443  const reco::CaloClusterPtr seed = aClus->seed();
444  hist_EB_RawSC_S1toS9_->Fill(lazyTool.eMax(*seed) / lazyTool.e3x3(*seed));
445  hist_EB_RawSC_S25toE_->Fill(lazyTool.e5x5(*seed) / aClus->energy());
446 
447  // truth
448  double dRClosest = 999.9;
449  double energyClosest = 0;
450  closestMCParticle(genEvent, *aClus, dRClosest, energyClosest);
451 
452  if (dRClosest < 0.1) {
453  hist_EB_RawSC_EoverTruth_->Fill(aClus->energy() / energyClosest);
454  hist_EB_RawSC_deltaR_->Fill(dRClosest);
455  }
456  }
457 
458  // Do CORRECTED BARREL SuperClusters
459  for (reco::SuperClusterCollection::const_iterator aClus = barrelCorSuperClusters->begin();
460  aClus != barrelCorSuperClusters->end();
461  aClus++) {
462  // kinematics
463  hist_EB_CorSC_NumBC_->Fill(aClus->clustersSize());
464  hist_EB_CorSC_ET_->Fill(aClus->energy() / std::cosh(aClus->position().eta()));
465  hist_EB_CorSC_Eta_->Fill(aClus->position().eta());
466  hist_EB_CorSC_Phi_->Fill(aClus->position().phi());
467 
468  hist_EB_CorSC_ET_vs_Eta_->Fill(aClus->energy() / std::cosh(aClus->position().eta()), aClus->eta());
469  hist_EB_CorSC_ET_vs_Phi_->Fill(aClus->energy() / std::cosh(aClus->position().eta()), aClus->phi());
470 
471  // cluster shape
472  const reco::CaloClusterPtr seed = aClus->seed();
473  hist_EB_CorSC_S1toS9_->Fill(lazyTool.eMax(*seed) / lazyTool.e3x3(*seed));
474  hist_EB_CorSC_S25toE_->Fill(lazyTool.e5x5(*seed) / aClus->energy());
475 
476  // correction variables
477  hist_EB_CorSC_phiWidth_->Fill(aClus->phiWidth());
478  hist_EB_CorSC_etaWidth_->Fill(aClus->etaWidth());
479 
480  // truth
481  double dRClosest = 999.9;
482  double energyClosest = 0;
483  closestMCParticle(genEvent, *aClus, dRClosest, energyClosest);
484 
485  if (dRClosest < 0.1) {
486  hist_EB_CorSC_EoverTruth_->Fill(aClus->energy() / energyClosest);
487  hist_EB_CorSC_deltaR_->Fill(dRClosest);
488  }
489  }
490 
491  //
492  // Get the ENDCAP products
493  //
494  edm::Handle<reco::SuperClusterCollection> pEndcapRawSuperClusters;
495  evt.getByToken(endcapRawSuperClusterCollectionToken_, pEndcapRawSuperClusters);
496  if (!pEndcapRawSuperClusters.isValid()) {
497  edm::LogError("EgammaSuperClusters") << "Error! can't get collection Raw EE SC";
498  }
499 
500  edm::Handle<reco::SuperClusterCollection> pEndcapPreSuperClusters;
501  evt.getByToken(endcapPreSuperClusterCollectionToken_, pEndcapPreSuperClusters);
502  if (!pEndcapPreSuperClusters.isValid()) {
503  edm::LogError("EgammaSuperClusters") << "Error! can't get collection Pre EE SC";
504  }
505 
506  edm::Handle<reco::SuperClusterCollection> pEndcapCorSuperClusters;
507  evt.getByToken(endcapCorSuperClusterCollectionToken_, pEndcapCorSuperClusters);
508  if (!pEndcapCorSuperClusters.isValid()) {
509  edm::LogError("EgammaSuperClusters") << "Error! can't get collection Cor EE SC";
510  }
511 
512  // Get the ENDCAP collections
513  const reco::SuperClusterCollection *endcapRawSuperClusters = pEndcapRawSuperClusters.product();
514  const reco::SuperClusterCollection *endcapPreSuperClusters = pEndcapPreSuperClusters.product();
515  const reco::SuperClusterCollection *endcapCorSuperClusters = pEndcapCorSuperClusters.product();
516 
517  // Number of entries in collections
518  hist_EE_RawSC_Size_->Fill(endcapRawSuperClusters->size());
519  hist_EE_PreSC_Size_->Fill(endcapPreSuperClusters->size());
520  hist_EE_CorSC_Size_->Fill(endcapCorSuperClusters->size());
521 
522  // Do RAW ENDCAP SuperClusters
523  for (reco::SuperClusterCollection::const_iterator aClus = endcapRawSuperClusters->begin();
524  aClus != endcapRawSuperClusters->end();
525  aClus++) {
526  hist_EE_RawSC_NumBC_->Fill(aClus->clustersSize());
527  hist_EE_RawSC_ET_->Fill(aClus->energy() / std::cosh(aClus->position().eta()));
528  hist_EE_RawSC_Eta_->Fill(aClus->position().eta());
529  hist_EE_RawSC_Phi_->Fill(aClus->position().phi());
530 
531  const reco::CaloClusterPtr seed = aClus->seed();
532  hist_EE_RawSC_S1toS9_->Fill(lazyTool.eMax(*seed) / lazyTool.e3x3(*seed));
533  hist_EE_RawSC_S25toE_->Fill(lazyTool.e5x5(*seed) / aClus->energy());
534 
535  // truth
536  double dRClosest = 999.9;
537  double energyClosest = 0;
538  closestMCParticle(genEvent, *aClus, dRClosest, energyClosest);
539 
540  if (dRClosest < 0.1) {
541  hist_EE_RawSC_EoverTruth_->Fill(aClus->energy() / energyClosest);
542  hist_EE_RawSC_deltaR_->Fill(dRClosest);
543  }
544  }
545 
546  // Do ENDCAP SuperClusters with PRESHOWER
547  for (reco::SuperClusterCollection::const_iterator aClus = endcapPreSuperClusters->begin();
548  aClus != endcapPreSuperClusters->end();
549  aClus++) {
550  hist_EE_PreSC_NumBC_->Fill(aClus->clustersSize());
551  hist_EE_PreSC_ET_->Fill(aClus->energy() / std::cosh(aClus->position().eta()));
552  hist_EE_PreSC_Eta_->Fill(aClus->position().eta());
553  hist_EE_PreSC_Phi_->Fill(aClus->position().phi());
554  hist_EE_PreSC_preshowerE_->Fill(aClus->preshowerEnergy());
555 
556  const reco::CaloClusterPtr seed = aClus->seed();
557  hist_EE_PreSC_S1toS9_->Fill(lazyTool.eMax(*seed) / lazyTool.e3x3(*seed));
558  hist_EE_PreSC_S25toE_->Fill(lazyTool.e5x5(*seed) / aClus->energy());
559 
560  // truth
561  double dRClosest = 999.9;
562  double energyClosest = 0;
563  closestMCParticle(genEvent, *aClus, dRClosest, energyClosest);
564 
565  if (dRClosest < 0.1) {
566  hist_EE_PreSC_EoverTruth_->Fill(aClus->energy() / energyClosest);
567  hist_EE_PreSC_deltaR_->Fill(dRClosest);
568  }
569  }
570 
571  // Do CORRECTED ENDCAP SuperClusters
572  for (reco::SuperClusterCollection::const_iterator aClus = endcapCorSuperClusters->begin();
573  aClus != endcapCorSuperClusters->end();
574  aClus++) {
575  hist_EE_CorSC_NumBC_->Fill(aClus->clustersSize());
576  hist_EE_CorSC_ET_->Fill(aClus->energy() / std::cosh(aClus->position().eta()));
577  hist_EE_CorSC_Eta_->Fill(aClus->position().eta());
578  hist_EE_CorSC_Phi_->Fill(aClus->position().phi());
579  hist_EE_CorSC_preshowerE_->Fill(aClus->preshowerEnergy());
580 
581  hist_EE_CorSC_ET_vs_Eta_->Fill(aClus->energy() / std::cosh(aClus->position().eta()), aClus->eta());
582  hist_EE_CorSC_ET_vs_Phi_->Fill(aClus->energy() / std::cosh(aClus->position().eta()), aClus->phi());
583  hist_EE_CorSC_ET_vs_R_->Fill(aClus->energy() / std::cosh(aClus->position().eta()),
584  std::sqrt(std::pow(aClus->x(), 2) + std::pow(aClus->y(), 2)));
585 
586  // correction variables
587  hist_EE_CorSC_phiWidth_->Fill(aClus->phiWidth());
588  hist_EE_CorSC_etaWidth_->Fill(aClus->etaWidth());
589 
590  const reco::CaloClusterPtr seed = aClus->seed();
591  hist_EE_CorSC_S1toS9_->Fill(lazyTool.eMax(*seed) / lazyTool.e3x3(*seed));
592  hist_EE_CorSC_S25toE_->Fill(lazyTool.e5x5(*seed) / aClus->energy());
593 
594  // truth
595  double dRClosest = 999.9;
596  double energyClosest = 0;
597  closestMCParticle(genEvent, *aClus, dRClosest, energyClosest);
598 
599  if (dRClosest < 0.1) {
600  hist_EE_CorSC_EoverTruth_->Fill(aClus->energy() / energyClosest);
601  hist_EE_CorSC_deltaR_->Fill(dRClosest);
602  }
603  }
604 }

References barrelCorSuperClusterCollectionToken_, barrelRawSuperClusterCollectionToken_, barrelRecHitCollectionToken_, closestMCParticle(), ecalClusterToolsESGetTokens_, endcapCorSuperClusterCollectionToken_, endcapPreSuperClusterCollectionToken_, endcapRawSuperClusterCollectionToken_, endcapRecHitCollectionToken_, dqm::impl::MonitorElement::Fill(), genWeightsTable_cfi::genEvent, EcalClusterLazyToolsBase::ESGetTokens::get(), edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), hist_EB_CorSC_deltaR_, hist_EB_CorSC_EoverTruth_, hist_EB_CorSC_ET_, hist_EB_CorSC_ET_vs_Eta_, hist_EB_CorSC_ET_vs_Phi_, hist_EB_CorSC_Eta_, hist_EB_CorSC_etaWidth_, hist_EB_CorSC_NumBC_, hist_EB_CorSC_Phi_, hist_EB_CorSC_phiWidth_, hist_EB_CorSC_S1toS9_, hist_EB_CorSC_S25toE_, hist_EB_CorSC_Size_, hist_EB_RawSC_deltaR_, hist_EB_RawSC_EoverTruth_, hist_EB_RawSC_ET_, hist_EB_RawSC_Eta_, hist_EB_RawSC_NumBC_, hist_EB_RawSC_Phi_, hist_EB_RawSC_S1toS9_, hist_EB_RawSC_S25toE_, hist_EB_RawSC_Size_, hist_EE_CorSC_deltaR_, hist_EE_CorSC_EoverTruth_, hist_EE_CorSC_ET_, hist_EE_CorSC_ET_vs_Eta_, hist_EE_CorSC_ET_vs_Phi_, hist_EE_CorSC_ET_vs_R_, hist_EE_CorSC_Eta_, hist_EE_CorSC_etaWidth_, hist_EE_CorSC_NumBC_, hist_EE_CorSC_Phi_, hist_EE_CorSC_phiWidth_, hist_EE_CorSC_preshowerE_, hist_EE_CorSC_S1toS9_, hist_EE_CorSC_S25toE_, hist_EE_CorSC_Size_, hist_EE_PreSC_deltaR_, hist_EE_PreSC_EoverTruth_, hist_EE_PreSC_ET_, hist_EE_PreSC_Eta_, hist_EE_PreSC_NumBC_, hist_EE_PreSC_Phi_, hist_EE_PreSC_preshowerE_, hist_EE_PreSC_S1toS9_, hist_EE_PreSC_S25toE_, hist_EE_PreSC_Size_, hist_EE_RawSC_deltaR_, hist_EE_RawSC_EoverTruth_, hist_EE_RawSC_ET_, hist_EE_RawSC_Eta_, hist_EE_RawSC_NumBC_, hist_EE_RawSC_Phi_, hist_EE_RawSC_S1toS9_, hist_EE_RawSC_S25toE_, hist_EE_RawSC_Size_, edm::HandleBase::isValid(), MCTruthCollectionToken_, funct::pow(), edm::Handle< T >::product(), fileCollector::seed, and mathSSE::sqrt().

◆ bookHistograms()

void EgammaSuperClusters::bookHistograms ( DQMStore::IBooker &  _ibooker,
edm::Run const &  ,
edm::EventSetup const &   
)
overridevirtual

Implements DQMEDAnalyzer.

Definition at line 105 of file EgammaSuperClusters.cc.

105  {
106  _ibooker.setCurrentFolder("EcalClusterV/EcalSuperClusters/");
107 
108  // Number of SuperClusters
109  //
110  hist_EB_RawSC_Size_ =
111  _ibooker.book1D("hist_EB_RawSC_Size_", "# Raw SuperClusters in Barrel", hsSize_.bins, hsSize_.min, hsSize_.max);
112  hist_EE_RawSC_Size_ =
113  _ibooker.book1D("hist_EE_RawSC_Size_", "# Raw SuperClusters in Endcap", hsSize_.bins, hsSize_.min, hsSize_.max);
114  hist_EB_CorSC_Size_ = _ibooker.book1D(
115  "hist_EB_CorSC_Size_", "# Corrected SuperClusters in Barrel", hsSize_.bins, hsSize_.min, hsSize_.max);
116  hist_EE_CorSC_Size_ = _ibooker.book1D(
117  "hist_EE_CorSC_Size_", "# Corrected SuperClusters in Endcap", hsSize_.bins, hsSize_.min, hsSize_.max);
118  hist_EE_PreSC_Size_ = _ibooker.book1D(
119  "hist_EE_PreSC_Size_", "# SuperClusters with Preshower in Endcap", hsSize_.bins, hsSize_.min, hsSize_.max);
120 
121  // Number of BasicClusters in SuperCluster
122  //
123  hist_EB_RawSC_NumBC_ = _ibooker.book1D("hist_EB_RawSC_NumBC_",
124  "# of Basic Clusters in Raw Super Clusters in Barrel",
125  hsNumBC_.bins,
126  hsNumBC_.min,
127  hsNumBC_.max);
128  hist_EE_RawSC_NumBC_ = _ibooker.book1D("hist_EE_RawSC_NumBC_",
129  "# of Basic Clusters in Raw Super Clusters in Endcap",
130  hsNumBC_.bins,
131  hsNumBC_.min,
132  hsNumBC_.max);
133  hist_EB_CorSC_NumBC_ = _ibooker.book1D("hist_EB_CorSC_NumBC_",
134  "# of Basic Clusters in Corrected SuperClusters in Barrel",
135  hsNumBC_.bins,
136  hsNumBC_.min,
137  hsNumBC_.max);
138  hist_EE_CorSC_NumBC_ = _ibooker.book1D("hist_EE_CorSC_NumBC_",
139  "# of Basic Clusters in Corrected SuperClusters in Endcap",
140  hsNumBC_.bins,
141  hsNumBC_.min,
142  hsNumBC_.max);
143  hist_EE_PreSC_NumBC_ = _ibooker.book1D("hist_EE_PreSC_NumBC_",
144  "# of Basic Clusters in SuperClusters with Preshower in Endcap",
145  hsNumBC_.bins,
146  hsNumBC_.min,
147  hsNumBC_.max);
148 
149  // ET distribution of SuperClusters
150  //
151  hist_EB_RawSC_ET_ =
152  _ibooker.book1D("hist_EB_RawSC_ET_", "ET of Raw SuperClusters in Barrel", hsET_.bins, hsET_.min, hsET_.max);
153  hist_EE_RawSC_ET_ =
154  _ibooker.book1D("hist_EE_RawSC_ET_", "ET of Raw SuperClusters in Endcap", hsET_.bins, hsET_.min, hsET_.max);
155  hist_EB_CorSC_ET_ =
156  _ibooker.book1D("hist_EB_CorSC_ET_", "ET of Corrected SuperClusters in Barrel", hsET_.bins, hsET_.min, hsET_.max);
157  hist_EE_CorSC_ET_ =
158  _ibooker.book1D("hist_EE_CorSC_ET_", "ET of Corrected SuperClusters in Endcap", hsET_.bins, hsET_.min, hsET_.max);
159  hist_EE_PreSC_ET_ = _ibooker.book1D(
160  "hist_EE_PreSC_ET_", "ET of SuperClusters with Preshower in Endcap", hsET_.bins, hsET_.min, hsET_.max);
161 
162  // Eta distribution of SuperClusters
163  //
164  hist_EB_RawSC_Eta_ =
165  _ibooker.book1D("hist_EB_RawSC_Eta_", "Eta of Raw SuperClusters in Barrel", hsEta_.bins, hsEta_.min, hsEta_.max);
166  hist_EE_RawSC_Eta_ =
167  _ibooker.book1D("hist_EE_RawSC_Eta_", "Eta of Raw SuperClusters in Endcap", hsEta_.bins, hsEta_.min, hsEta_.max);
168  hist_EB_CorSC_Eta_ = _ibooker.book1D(
169  "hist_EB_CorSC_Eta_", "Eta of Corrected SuperClusters in Barrel", hsEta_.bins, hsEta_.min, hsEta_.max);
170  hist_EE_CorSC_Eta_ = _ibooker.book1D(
171  "hist_EE_CorSC_Eta_", "Eta of Corrected SuperClusters in Endcap", hsEta_.bins, hsEta_.min, hsEta_.max);
172  hist_EE_PreSC_Eta_ = _ibooker.book1D(
173  "hist_EE_PreSC_Eta_", "Eta of SuperClusters with Preshower in Endcap", hsEta_.bins, hsEta_.min, hsEta_.max);
174 
175  // Phi distribution of SuperClusters
176  //
177  hist_EB_RawSC_Phi_ =
178  _ibooker.book1D("hist_EB_RawSC_Phi_", "Phi of Raw SuperClusters in Barrel", hsPhi_.bins, hsPhi_.min, hsPhi_.max);
179  hist_EE_RawSC_Phi_ =
180  _ibooker.book1D("hist_EE_RawSC_Phi_", "Phi of Raw SuperClusters in Endcap", hsPhi_.bins, hsPhi_.min, hsPhi_.max);
181  hist_EB_CorSC_Phi_ = _ibooker.book1D(
182  "hist_EB_CorSC_Phi_", "Phi of Corrected SuperClusters in Barrel", hsPhi_.bins, hsPhi_.min, hsPhi_.max);
183  hist_EE_CorSC_Phi_ = _ibooker.book1D(
184  "hist_EE_CorSC_Phi_", "Phi of Corrected SuperClusters in Endcap", hsPhi_.bins, hsPhi_.min, hsPhi_.max);
185  hist_EE_PreSC_Phi_ = _ibooker.book1D(
186  "hist_EE_PreSC_Phi_", "Phi of SuperClusters with Preshower in Endcap", hsPhi_.bins, hsPhi_.min, hsPhi_.max);
187 
188  // S1/S9 distribution of SuperClusters
189  //
190  hist_EB_RawSC_S1toS9_ = _ibooker.book1D(
191  "hist_EB_RawSC_S1toS9_", "S1/S9 of Raw Super Clusters in Barrel", hsS1toS9_.bins, hsS1toS9_.min, hsS1toS9_.max);
192  hist_EE_RawSC_S1toS9_ = _ibooker.book1D(
193  "hist_EE_RawSC_S1toS9_", "S1/S9 of Raw Super Clusters in Endcap", hsS1toS9_.bins, hsS1toS9_.min, hsS1toS9_.max);
194  hist_EB_CorSC_S1toS9_ = _ibooker.book1D("hist_EB_CorSC_S1toS9_",
195  "S1/S9 of Corrected SuperClusters in Barrel",
196  hsS1toS9_.bins,
197  hsS1toS9_.min,
198  hsS1toS9_.max);
199  hist_EE_CorSC_S1toS9_ = _ibooker.book1D("hist_EE_CorSC_S1toS9_",
200  "S1/S9 of Corrected SuperClusters in Endcap",
201  hsS1toS9_.bins,
202  hsS1toS9_.min,
203  hsS1toS9_.max);
204  hist_EE_PreSC_S1toS9_ = _ibooker.book1D("hist_EE_PreSC_S1toS9_",
205  "S1/S9 of SuperClusters with Preshower in Endcap",
206  hsS1toS9_.bins,
207  hsS1toS9_.min,
208  hsS1toS9_.max);
209 
210  // S25/E distribution of SuperClusters
211  //
212  hist_EB_RawSC_S25toE_ = _ibooker.book1D(
213  "hist_EB_RawSC_S25toE_", "S25/E of Raw Super Clusters in Barrel", hsS25toE_.bins, hsS25toE_.min, hsS25toE_.max);
214  hist_EE_RawSC_S25toE_ = _ibooker.book1D(
215  "hist_EE_RawSC_S25toE_", "S25/E of Raw Super Clusters in Endcap", hsS25toE_.bins, hsS25toE_.min, hsS25toE_.max);
216  hist_EB_CorSC_S25toE_ = _ibooker.book1D("hist_EB_CorSC_S25toE_",
217  "S25/E of Corrected SuperClusters in Barrel",
218  hsS25toE_.bins,
219  hsS25toE_.min,
220  hsS25toE_.max);
221  hist_EE_CorSC_S25toE_ = _ibooker.book1D("hist_EE_CorSC_S25toE_",
222  "S25/E of Corrected SuperClusters in Endcap",
223  hsS25toE_.bins,
224  hsS25toE_.min,
225  hsS25toE_.max);
226  hist_EE_PreSC_S25toE_ = _ibooker.book1D("hist_EE_PreSC_S25toE_",
227  "S25/E of SuperClusters with Preshower in Endcap",
228  hsS25toE_.bins,
229  hsS25toE_.min,
230  hsS25toE_.max);
231 
232  // E/E(true) distribution of SuperClusters
233  //
234  hist_EB_RawSC_EoverTruth_ = _ibooker.book1D("hist_EB_RawSC_EoverTruth_",
235  "E/True E of Raw SuperClusters in Barrel",
236  hsEoverTruth_.bins,
237  hsEoverTruth_.min,
238  hsEoverTruth_.max);
239  hist_EE_RawSC_EoverTruth_ = _ibooker.book1D("hist_EE_RawSC_EoverTruth_",
240  "E/True E of Raw SuperClusters in Endcap",
241  hsEoverTruth_.bins,
242  hsEoverTruth_.min,
243  hsEoverTruth_.max);
244  hist_EB_CorSC_EoverTruth_ = _ibooker.book1D("hist_EB_CorSC_EoverTruth_",
245  "E/True E of Corrected SuperClusters in Barrel",
246  hsEoverTruth_.bins,
247  hsEoverTruth_.min,
248  hsEoverTruth_.max);
249  hist_EE_CorSC_EoverTruth_ = _ibooker.book1D("hist_EE_CorSC_EoverTruth_",
250  "E/True E of Corrected SuperClusters in Endcap",
251  hsEoverTruth_.bins,
252  hsEoverTruth_.min,
253  hsEoverTruth_.max);
254  hist_EE_PreSC_EoverTruth_ = _ibooker.book1D("hist_EE_PreSC_EoverTruth_",
255  "E/True E of SuperClusters with Preshower in Endcap",
256  hsEoverTruth_.bins,
257  hsEoverTruth_.min,
258  hsEoverTruth_.max);
259 
260  // dR distribution of SuperClusters from truth
261  //
262  hist_EB_RawSC_deltaR_ = _ibooker.book1D("hist_EB_RawSC_deltaR_",
263  "dR to MC truth of Raw Super Clusters in Barrel",
264  hsdeltaR_.bins,
265  hsdeltaR_.min,
266  hsdeltaR_.max);
267  hist_EE_RawSC_deltaR_ = _ibooker.book1D("hist_EE_RawSC_deltaR_",
268  "dR to MC truth of Raw Super Clusters in Endcap",
269  hsdeltaR_.bins,
270  hsdeltaR_.min,
271  hsdeltaR_.max);
272  hist_EB_CorSC_deltaR_ = _ibooker.book1D("hist_EB_CorSC_deltaR_",
273  "dR to MC truth of Corrected SuperClusters in Barrel",
274  hsdeltaR_.bins,
275  hsdeltaR_.min,
276  hsdeltaR_.max);
277  hist_EE_CorSC_deltaR_ = _ibooker.book1D("hist_EE_CorSC_deltaR_",
278  "dR to MC truth of Corrected SuperClusters in Endcap",
279  hsdeltaR_.bins,
280  hsdeltaR_.min,
281  hsdeltaR_.max);
282  hist_EE_PreSC_deltaR_ = _ibooker.book1D("hist_EE_PreSC_deltaR_",
283  "dR to MC truth of SuperClusters with Preshower in Endcap",
284  hsdeltaR_.bins,
285  hsdeltaR_.min,
286  hsdeltaR_.max);
287 
288  // phi width stored in corrected SuperClusters
289  hist_EB_CorSC_phiWidth_ = _ibooker.book1D("hist_EB_CorSC_phiWidth_",
290  "phiWidth of Corrected Super Clusters in Barrel",
291  hsphiWidth_.bins,
292  hsphiWidth_.min,
293  hsphiWidth_.max);
294  hist_EE_CorSC_phiWidth_ = _ibooker.book1D("hist_EE_CorSC_phiWidth_",
295  "phiWidth of Corrected Super Clusters in Endcap",
296  hsphiWidth_.bins,
297  hsphiWidth_.min,
298  hsphiWidth_.max);
299 
300  // eta width stored in corrected SuperClusters
301  hist_EB_CorSC_etaWidth_ = _ibooker.book1D("hist_EB_CorSC_etaWidth_",
302  "etaWidth of Corrected Super Clusters in Barrel",
303  hsetaWidth_.bins,
304  hsetaWidth_.min,
305  hsetaWidth_.max);
306  hist_EE_CorSC_etaWidth_ = _ibooker.book1D("hist_EE_CorSC_etaWidth_",
307  "etaWidth of Corrected Super Clusters in Endcap",
308  hsetaWidth_.bins,
309  hsetaWidth_.min,
310  hsetaWidth_.max);
311 
312  // preshower energy
313  hist_EE_PreSC_preshowerE_ = _ibooker.book1D("hist_EE_PreSC_preshowerE_",
314  "preshower energy in Super Clusters with Preshower in Endcap",
315  hspreshowerE_.bins,
316  hspreshowerE_.min,
317  hspreshowerE_.max);
318  hist_EE_CorSC_preshowerE_ = _ibooker.book1D("hist_EE_CorSC_preshowerE_",
319  "preshower energy in Corrected Super Clusters with Preshower in Endcap",
320  hspreshowerE_.bins,
321  hspreshowerE_.min,
322  hspreshowerE_.max);
323 
324  //
325  hist_EB_CorSC_ET_vs_Eta_ = _ibooker.book2D("hist_EB_CorSC_ET_vs_Eta_",
326  "Corr Super Cluster ET versus Eta in Barrel",
327  hsET_.bins,
328  hsET_.min,
329  hsET_.max,
330  hsEta_.bins,
331  hsEta_.min,
332  hsEta_.max);
333 
334  hist_EB_CorSC_ET_vs_Phi_ = _ibooker.book2D("hist_EB_CorSC_ET_vs_Phi_",
335  "Corr Super Cluster ET versus Phi in Barrel",
336  hsET_.bins,
337  hsET_.min,
338  hsET_.max,
339  hsPhi_.bins,
340  hsPhi_.min,
341  hsPhi_.max);
342 
343  hist_EE_CorSC_ET_vs_Eta_ = _ibooker.book2D("hist_EE_CorSC_ET_vs_Eta_",
344  "Corr Super Cluster ET versus Eta in Endcap",
345  hsET_.bins,
346  hsET_.min,
347  hsET_.max,
348  hsEta_.bins,
349  hsEta_.min,
350  hsEta_.max);
351 
352  hist_EE_CorSC_ET_vs_Phi_ = _ibooker.book2D("hist_EE_CorSC_ET_vs_Phi_",
353  "Corr Super Cluster ET versus Phi in Endcap",
354  hsET_.bins,
355  hsET_.min,
356  hsET_.max,
357  hsPhi_.bins,
358  hsPhi_.min,
359  hsPhi_.max);
360 
361  hist_EE_CorSC_ET_vs_R_ = _ibooker.book2D("hist_EE_CorSC_ET_vs_R_",
362  "Corr Super Cluster ET versus Radius in Endcap",
363  hsET_.bins,
364  hsET_.min,
365  hsET_.max,
366  hsR_.bins,
367  hsR_.min,
368  hsR_.max);
369 }

References HistSpec::bins, dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), hist_EB_CorSC_deltaR_, hist_EB_CorSC_EoverTruth_, hist_EB_CorSC_ET_, hist_EB_CorSC_ET_vs_Eta_, hist_EB_CorSC_ET_vs_Phi_, hist_EB_CorSC_Eta_, hist_EB_CorSC_etaWidth_, hist_EB_CorSC_NumBC_, hist_EB_CorSC_Phi_, hist_EB_CorSC_phiWidth_, hist_EB_CorSC_S1toS9_, hist_EB_CorSC_S25toE_, hist_EB_CorSC_Size_, hist_EB_RawSC_deltaR_, hist_EB_RawSC_EoverTruth_, hist_EB_RawSC_ET_, hist_EB_RawSC_Eta_, hist_EB_RawSC_NumBC_, hist_EB_RawSC_Phi_, hist_EB_RawSC_S1toS9_, hist_EB_RawSC_S25toE_, hist_EB_RawSC_Size_, hist_EE_CorSC_deltaR_, hist_EE_CorSC_EoverTruth_, hist_EE_CorSC_ET_, hist_EE_CorSC_ET_vs_Eta_, hist_EE_CorSC_ET_vs_Phi_, hist_EE_CorSC_ET_vs_R_, hist_EE_CorSC_Eta_, hist_EE_CorSC_etaWidth_, hist_EE_CorSC_NumBC_, hist_EE_CorSC_Phi_, hist_EE_CorSC_phiWidth_, hist_EE_CorSC_preshowerE_, hist_EE_CorSC_S1toS9_, hist_EE_CorSC_S25toE_, hist_EE_CorSC_Size_, hist_EE_PreSC_deltaR_, hist_EE_PreSC_EoverTruth_, hist_EE_PreSC_ET_, hist_EE_PreSC_Eta_, hist_EE_PreSC_NumBC_, hist_EE_PreSC_Phi_, hist_EE_PreSC_preshowerE_, hist_EE_PreSC_S1toS9_, hist_EE_PreSC_S25toE_, hist_EE_PreSC_Size_, hist_EE_RawSC_deltaR_, hist_EE_RawSC_EoverTruth_, hist_EE_RawSC_ET_, hist_EE_RawSC_Eta_, hist_EE_RawSC_NumBC_, hist_EE_RawSC_Phi_, hist_EE_RawSC_S1toS9_, hist_EE_RawSC_S25toE_, hist_EE_RawSC_Size_, hsdeltaR_, hsEoverTruth_, hsET_, hsEta_, hsetaWidth_, hsNumBC_, hsPhi_, hsphiWidth_, hspreshowerE_, hsR_, hsS1toS9_, hsS25toE_, hsSize_, HistSpec::max, HistSpec::min, and dqm::implementation::NavigatorBase::setCurrentFolder().

◆ closestMCParticle()

void EgammaSuperClusters::closestMCParticle ( HepMC::GenEvent const *  genEvent,
reco::SuperCluster const &  sc,
double &  dRClosest,
double &  energyClosest 
) const
private

Definition at line 609 of file EgammaSuperClusters.cc.

612  {
613  // SuperCluster eta, phi
614  double scEta = sc.eta();
615  double scPhi = sc.phi();
616 
617  // initialize dRClosest to a large number
618  dRClosest = 999.9;
619 
620  // loop over the MC truth particles to find the
621  // closest to the superCluster in dR space
622  for (HepMC::GenEvent::particle_const_iterator currentParticle = genEvent->particles_begin();
623  currentParticle != genEvent->particles_end();
624  currentParticle++) {
625  if ((*currentParticle)->status() == 1) {
626  // need GenParticle in ECAL co-ordinates
627  HepMC::FourVector vtx = (*currentParticle)->production_vertex()->position();
628  double phiTrue = (*currentParticle)->momentum().phi();
629  double etaTrue = ecalEta((*currentParticle)->momentum().eta(), vtx.z() / 10., vtx.perp() / 10.);
630 
631  double dPhi = reco::deltaPhi(phiTrue, scPhi);
632  double dEta = scEta - etaTrue;
633  double deltaR = std::sqrt(dPhi * dPhi + dEta * dEta);
634 
635  if (deltaR < dRClosest) {
636  dRClosest = deltaR;
637  energyClosest = (*currentParticle)->momentum().e();
638  }
639 
640  } // end if stable particle
641 
642  } // end loop on get particles
643 }

References reco::deltaPhi(), PbPb_ZMuSkimMuonDPG_cff::deltaR, HLT_FULL_cff::dEta, HLT_FULL_cff::dPhi, ecalEta(), reco::CaloCluster::eta(), genWeightsTable_cfi::genEvent, reco::CaloCluster::phi(), mathSSE::sqrt(), and extraflags_cff::vtx.

Referenced by analyze().

◆ ecalEta()

float EgammaSuperClusters::ecalEta ( float  EtaParticle,
float  Zvertex,
float  plane_Radius 
) const
private

Definition at line 648 of file EgammaSuperClusters.cc.

648  {
649  const float R_ECAL = 136.5;
650  const float Z_Endcap = 328.0;
651  const float etaBarrelEndcap = 1.479;
652 
653  if (EtaParticle != 0.) {
654  float Theta = 0.0;
655  float ZEcal = (R_ECAL - plane_Radius) * sinh(EtaParticle) + Zvertex;
656 
657  if (ZEcal != 0.0)
658  Theta = atan(R_ECAL / ZEcal);
659  if (Theta < 0.0)
660  Theta = Theta + Geom::pi();
661 
662  float ETA = -log(tan(0.5 * Theta));
663 
664  if (fabs(ETA) > etaBarrelEndcap) {
665  float Zend = Z_Endcap;
666  if (EtaParticle < 0.0)
667  Zend = -Zend;
668  float Zlen = Zend - Zvertex;
669  float RR = Zlen / sinh(EtaParticle);
670  Theta = atan((RR + plane_Radius) / Zend);
671  if (Theta < 0.0)
672  Theta = Theta + Geom::pi();
673  ETA = -log(tan(0.5 * Theta));
674  }
675 
676  return ETA;
677  } else {
678  edm::LogWarning("") << "[EgammaSuperClusters::ecalEta] Warning: Eta equals "
679  "to zero, not correcting";
680  return EtaParticle;
681  }
682 }

References ETA, etaBarrelEndcap, dqm-mbProfile::log, Geom::pi(), R_ECAL, funct::tan(), Z_Endcap, and ZEcal.

Referenced by closestMCParticle().

Member Data Documentation

◆ barrelCorSuperClusterCollectionToken_

edm::EDGetTokenT<reco::SuperClusterCollection> EgammaSuperClusters::barrelCorSuperClusterCollectionToken_
private

Definition at line 41 of file EgammaSuperClusters.h.

Referenced by analyze().

◆ barrelRawSuperClusterCollectionToken_

edm::EDGetTokenT<reco::SuperClusterCollection> EgammaSuperClusters::barrelRawSuperClusterCollectionToken_
private

Definition at line 40 of file EgammaSuperClusters.h.

Referenced by analyze().

◆ barrelRecHitCollectionToken_

edm::EDGetTokenT<EcalRecHitCollection> EgammaSuperClusters::barrelRecHitCollectionToken_
private

Definition at line 49 of file EgammaSuperClusters.h.

Referenced by analyze().

◆ ecalClusterToolsESGetTokens_

EcalClusterLazyTools::ESGetTokens EgammaSuperClusters::ecalClusterToolsESGetTokens_
private

Definition at line 52 of file EgammaSuperClusters.h.

Referenced by analyze().

◆ endcapCorSuperClusterCollectionToken_

edm::EDGetTokenT<reco::SuperClusterCollection> EgammaSuperClusters::endcapCorSuperClusterCollectionToken_
private

Definition at line 46 of file EgammaSuperClusters.h.

Referenced by analyze().

◆ endcapPreSuperClusterCollectionToken_

edm::EDGetTokenT<reco::SuperClusterCollection> EgammaSuperClusters::endcapPreSuperClusterCollectionToken_
private

Definition at line 45 of file EgammaSuperClusters.h.

Referenced by analyze().

◆ endcapRawSuperClusterCollectionToken_

edm::EDGetTokenT<reco::SuperClusterCollection> EgammaSuperClusters::endcapRawSuperClusterCollectionToken_
private

Definition at line 44 of file EgammaSuperClusters.h.

Referenced by analyze().

◆ endcapRecHitCollectionToken_

edm::EDGetTokenT<EcalRecHitCollection> EgammaSuperClusters::endcapRecHitCollectionToken_
private

Definition at line 50 of file EgammaSuperClusters.h.

Referenced by analyze().

◆ hist_EB_CorSC_deltaR_

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_deltaR_
private

Definition at line 110 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_CorSC_EoverTruth_

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_EoverTruth_
private

Definition at line 105 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_CorSC_ET_

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_ET_
private

Definition at line 80 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_CorSC_ET_vs_Eta_

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_ET_vs_Eta_
private

Definition at line 119 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_CorSC_ET_vs_Phi_

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_ET_vs_Phi_
private

Definition at line 120 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_CorSC_Eta_

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_Eta_
private

Definition at line 85 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_CorSC_etaWidth_

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_etaWidth_
private

Definition at line 118 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_CorSC_NumBC_

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_NumBC_
private

Definition at line 75 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_CorSC_Phi_

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_Phi_
private

Definition at line 90 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_CorSC_phiWidth_

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_phiWidth_
private

Definition at line 116 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_CorSC_S1toS9_

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_S1toS9_
private

Definition at line 95 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_CorSC_S25toE_

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_S25toE_
private

Definition at line 100 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_CorSC_Size_

MonitorElement* EgammaSuperClusters::hist_EB_CorSC_Size_
private

Definition at line 70 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_RawSC_deltaR_

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_deltaR_
private

Definition at line 108 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_RawSC_EoverTruth_

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_EoverTruth_
private

Definition at line 103 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_RawSC_ET_

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_ET_
private

Definition at line 78 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_RawSC_Eta_

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_Eta_
private

Definition at line 83 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_RawSC_NumBC_

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_NumBC_
private

Definition at line 73 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_RawSC_Phi_

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_Phi_
private

Definition at line 88 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_RawSC_S1toS9_

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_S1toS9_
private

Definition at line 93 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_RawSC_S25toE_

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_S25toE_
private

Definition at line 98 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EB_RawSC_Size_

MonitorElement* EgammaSuperClusters::hist_EB_RawSC_Size_
private

Definition at line 68 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_deltaR_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_deltaR_
private

Definition at line 111 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_EoverTruth_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_EoverTruth_
private

Definition at line 106 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_ET_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_ET_
private

Definition at line 81 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_ET_vs_Eta_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_ET_vs_Eta_
private

Definition at line 121 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_ET_vs_Phi_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_ET_vs_Phi_
private

Definition at line 122 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_ET_vs_R_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_ET_vs_R_
private

Definition at line 123 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_Eta_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_Eta_
private

Definition at line 86 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_etaWidth_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_etaWidth_
private

Definition at line 117 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_NumBC_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_NumBC_
private

Definition at line 76 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_Phi_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_Phi_
private

Definition at line 91 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_phiWidth_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_phiWidth_
private

Definition at line 115 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_preshowerE_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_preshowerE_
private

Definition at line 114 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_S1toS9_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_S1toS9_
private

Definition at line 96 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_S25toE_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_S25toE_
private

Definition at line 101 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_CorSC_Size_

MonitorElement* EgammaSuperClusters::hist_EE_CorSC_Size_
private

Definition at line 71 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_PreSC_deltaR_

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_deltaR_
private

Definition at line 112 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_PreSC_EoverTruth_

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_EoverTruth_
private

Definition at line 107 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_PreSC_ET_

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_ET_
private

Definition at line 82 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_PreSC_Eta_

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_Eta_
private

Definition at line 87 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_PreSC_NumBC_

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_NumBC_
private

Definition at line 77 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_PreSC_Phi_

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_Phi_
private

Definition at line 92 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_PreSC_preshowerE_

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_preshowerE_
private

Definition at line 113 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_PreSC_S1toS9_

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_S1toS9_
private

Definition at line 97 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_PreSC_S25toE_

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_S25toE_
private

Definition at line 102 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_PreSC_Size_

MonitorElement* EgammaSuperClusters::hist_EE_PreSC_Size_
private

Definition at line 72 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_RawSC_deltaR_

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_deltaR_
private

Definition at line 109 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_RawSC_EoverTruth_

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_EoverTruth_
private

Definition at line 104 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_RawSC_ET_

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_ET_
private

Definition at line 79 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_RawSC_Eta_

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_Eta_
private

Definition at line 84 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_RawSC_NumBC_

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_NumBC_
private

Definition at line 74 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_RawSC_Phi_

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_Phi_
private

Definition at line 89 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_RawSC_S1toS9_

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_S1toS9_
private

Definition at line 94 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_RawSC_S25toE_

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_S25toE_
private

Definition at line 99 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hist_EE_RawSC_Size_

MonitorElement* EgammaSuperClusters::hist_EE_RawSC_Size_
private

Definition at line 69 of file EgammaSuperClusters.h.

Referenced by analyze(), and bookHistograms().

◆ hsdeltaR_

HistSpec EgammaSuperClusters::hsdeltaR_
private

Definition at line 62 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

◆ hsEoverTruth_

HistSpec EgammaSuperClusters::hsEoverTruth_
private

Definition at line 61 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

◆ hsET_

HistSpec EgammaSuperClusters::hsET_
private

Definition at line 56 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

◆ hsEta_

HistSpec EgammaSuperClusters::hsEta_
private

Definition at line 57 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

◆ hsetaWidth_

HistSpec EgammaSuperClusters::hsetaWidth_
private

Definition at line 64 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

◆ hsNumBC_

HistSpec EgammaSuperClusters::hsNumBC_
private

Definition at line 55 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

◆ hsPhi_

HistSpec EgammaSuperClusters::hsPhi_
private

Definition at line 58 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

◆ hsphiWidth_

HistSpec EgammaSuperClusters::hsphiWidth_
private

Definition at line 63 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

◆ hspreshowerE_

HistSpec EgammaSuperClusters::hspreshowerE_
private

Definition at line 65 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

◆ hsR_

HistSpec EgammaSuperClusters::hsR_
private

Definition at line 66 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

◆ hsS1toS9_

HistSpec EgammaSuperClusters::hsS1toS9_
private

Definition at line 59 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

◆ hsS25toE_

HistSpec EgammaSuperClusters::hsS25toE_
private

Definition at line 60 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

◆ hsSize_

HistSpec EgammaSuperClusters::hsSize_
private

Definition at line 54 of file EgammaSuperClusters.h.

Referenced by bookHistograms().

◆ MCTruthCollectionToken_

edm::EDGetTokenT<edm::HepMCProduct> EgammaSuperClusters::MCTruthCollectionToken_
private

Definition at line 37 of file EgammaSuperClusters.h.

Referenced by analyze().

EgammaSuperClusters::hist_EE_RawSC_ET_
MonitorElement * hist_EE_RawSC_ET_
Definition: EgammaSuperClusters.h:79
EgammaSuperClusters::hist_EB_RawSC_deltaR_
MonitorElement * hist_EB_RawSC_deltaR_
Definition: EgammaSuperClusters.h:108
EgammaSuperClusters::barrelRawSuperClusterCollectionToken_
edm::EDGetTokenT< reco::SuperClusterCollection > barrelRawSuperClusterCollectionToken_
Definition: EgammaSuperClusters.h:40
EgammaSuperClusters::hist_EB_CorSC_etaWidth_
MonitorElement * hist_EB_CorSC_etaWidth_
Definition: EgammaSuperClusters.h:118
EgammaSuperClusters::hist_EB_CorSC_phiWidth_
MonitorElement * hist_EB_CorSC_phiWidth_
Definition: EgammaSuperClusters.h:116
EgammaSuperClusters::hist_EB_RawSC_NumBC_
MonitorElement * hist_EB_RawSC_NumBC_
Definition: EgammaSuperClusters.h:73
edm::Handle::product
T const * product() const
Definition: Handle.h:70
EgammaSuperClusters::endcapPreSuperClusterCollectionToken_
edm::EDGetTokenT< reco::SuperClusterCollection > endcapPreSuperClusterCollectionToken_
Definition: EgammaSuperClusters.h:45
EgammaSuperClusters::hist_EE_CorSC_phiWidth_
MonitorElement * hist_EE_CorSC_phiWidth_
Definition: EgammaSuperClusters.h:115
EgammaSuperClusters::hsEoverTruth_
HistSpec hsEoverTruth_
Definition: EgammaSuperClusters.h:61
EgammaSuperClusters::hist_EE_CorSC_NumBC_
MonitorElement * hist_EE_CorSC_NumBC_
Definition: EgammaSuperClusters.h:76
reco::deltaPhi
constexpr double deltaPhi(double phi1, double phi2)
Definition: deltaPhi.h:26
EgammaSuperClusters::hist_EB_CorSC_NumBC_
MonitorElement * hist_EB_CorSC_NumBC_
Definition: EgammaSuperClusters.h:75
EgammaSuperClusters::hist_EE_CorSC_deltaR_
MonitorElement * hist_EE_CorSC_deltaR_
Definition: EgammaSuperClusters.h:111
EgammaSuperClusters::hsSize_
HistSpec hsSize_
Definition: EgammaSuperClusters.h:54
EgammaSuperClusters::hist_EB_RawSC_Phi_
MonitorElement * hist_EB_RawSC_Phi_
Definition: EgammaSuperClusters.h:88
EgammaSuperClusters::MCTruthCollectionToken_
edm::EDGetTokenT< edm::HepMCProduct > MCTruthCollectionToken_
Definition: EgammaSuperClusters.h:37
EgammaSuperClusters::endcapRawSuperClusterCollectionToken_
edm::EDGetTokenT< reco::SuperClusterCollection > endcapRawSuperClusterCollectionToken_
Definition: EgammaSuperClusters.h:44
dqm::implementation::NavigatorBase::setCurrentFolder
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:32
EgammaSuperClusters::barrelCorSuperClusterCollectionToken_
edm::EDGetTokenT< reco::SuperClusterCollection > barrelCorSuperClusterCollectionToken_
Definition: EgammaSuperClusters.h:41
EgammaSuperClusters::hist_EE_PreSC_Phi_
MonitorElement * hist_EE_PreSC_Phi_
Definition: EgammaSuperClusters.h:92
HistSpec::bins
int bins
Definition: HistSpec.h:9
EgammaSuperClusters::hsS1toS9_
HistSpec hsS1toS9_
Definition: EgammaSuperClusters.h:59
EgammaSuperClusters::hsetaWidth_
HistSpec hsetaWidth_
Definition: EgammaSuperClusters.h:64
EgammaSuperClusters::hist_EE_PreSC_S1toS9_
MonitorElement * hist_EE_PreSC_S1toS9_
Definition: EgammaSuperClusters.h:97
EgammaSuperClusters::hist_EE_PreSC_Size_
MonitorElement * hist_EE_PreSC_Size_
Definition: EgammaSuperClusters.h:72
HLT_FULL_cff.dPhi
dPhi
Definition: HLT_FULL_cff.py:13702
edm::Handle< edm::HepMCProduct >
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:122
EgammaSuperClusters::hist_EE_CorSC_ET_vs_Phi_
MonitorElement * hist_EE_CorSC_ET_vs_Phi_
Definition: EgammaSuperClusters.h:122
EgammaSuperClusters::hist_EB_RawSC_Eta_
MonitorElement * hist_EB_RawSC_Eta_
Definition: EgammaSuperClusters.h:83
HepMC::GenEvent
Definition: hepmc_rootio.cc:9
EgammaSuperClusters::hist_EE_RawSC_S1toS9_
MonitorElement * hist_EE_RawSC_S1toS9_
Definition: EgammaSuperClusters.h:94
EgammaSuperClusters::hist_EE_CorSC_EoverTruth_
MonitorElement * hist_EE_CorSC_EoverTruth_
Definition: EgammaSuperClusters.h:106
fileCollector.seed
seed
Definition: fileCollector.py:127
reco::SuperClusterCollection
std::vector< SuperCluster > SuperClusterCollection
collection of SuperCluser objectr
Definition: SuperClusterFwd.h:9
EgammaSuperClusters::hist_EE_PreSC_ET_
MonitorElement * hist_EE_PreSC_ET_
Definition: EgammaSuperClusters.h:82
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
dqm::impl::MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:290
EgammaSuperClusters::hist_EB_CorSC_ET_vs_Eta_
MonitorElement * hist_EB_CorSC_ET_vs_Eta_
Definition: EgammaSuperClusters.h:119
EgammaSuperClusters::hist_EE_CorSC_Phi_
MonitorElement * hist_EE_CorSC_Phi_
Definition: EgammaSuperClusters.h:91
HistSpec::max
double max
Definition: HistSpec.h:8
EgammaSuperClusters::hist_EB_CorSC_EoverTruth_
MonitorElement * hist_EB_CorSC_EoverTruth_
Definition: EgammaSuperClusters.h:105
Geom::pi
constexpr double pi()
Definition: Pi.h:31
EgammaSuperClusters::hist_EB_RawSC_S1toS9_
MonitorElement * hist_EB_RawSC_S1toS9_
Definition: EgammaSuperClusters.h:93
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:531
EgammaSuperClusters::hist_EE_RawSC_EoverTruth_
MonitorElement * hist_EE_RawSC_EoverTruth_
Definition: EgammaSuperClusters.h:104
EgammaSuperClusters::ecalEta
float ecalEta(float, float, float) const
Definition: EgammaSuperClusters.cc:648
EgammaSuperClusters::hist_EE_PreSC_Eta_
MonitorElement * hist_EE_PreSC_Eta_
Definition: EgammaSuperClusters.h:87
EgammaSuperClusters::hist_EB_RawSC_ET_
MonitorElement * hist_EB_RawSC_ET_
Definition: EgammaSuperClusters.h:78
EgammaSuperClusters::hist_EE_CorSC_ET_
MonitorElement * hist_EE_CorSC_ET_
Definition: EgammaSuperClusters.h:81
EcalClusterLazyTools
PbPb_ZMuSkimMuonDPG_cff.deltaR
deltaR
Definition: PbPb_ZMuSkimMuonDPG_cff.py:63
EgammaSuperClusters::hist_EE_CorSC_preshowerE_
MonitorElement * hist_EE_CorSC_preshowerE_
Definition: EgammaSuperClusters.h:114
EgammaSuperClusters::hsphiWidth_
HistSpec hsphiWidth_
Definition: EgammaSuperClusters.h:63
EgammaSuperClusters::hist_EE_CorSC_Size_
MonitorElement * hist_EE_CorSC_Size_
Definition: EgammaSuperClusters.h:71
EgammaSuperClusters::endcapRecHitCollectionToken_
edm::EDGetTokenT< EcalRecHitCollection > endcapRecHitCollectionToken_
Definition: EgammaSuperClusters.h:50
EgammaSuperClusters::closestMCParticle
void closestMCParticle(HepMC::GenEvent const *, reco::SuperCluster const &, double &, double &) const
Definition: EgammaSuperClusters.cc:609
EgammaSuperClusters::hist_EE_RawSC_deltaR_
MonitorElement * hist_EE_RawSC_deltaR_
Definition: EgammaSuperClusters.h:109
EgammaSuperClusters::hsNumBC_
HistSpec hsNumBC_
Definition: EgammaSuperClusters.h:55
EgammaSuperClusters::hist_EB_CorSC_Phi_
MonitorElement * hist_EB_CorSC_Phi_
Definition: EgammaSuperClusters.h:90
funct::tan
Tan< T >::type tan(const T &t)
Definition: Tan.h:22
EgammaSuperClusters::hsdeltaR_
HistSpec hsdeltaR_
Definition: EgammaSuperClusters.h:62
EgammaSuperClusters::hist_EE_PreSC_EoverTruth_
MonitorElement * hist_EE_PreSC_EoverTruth_
Definition: EgammaSuperClusters.h:107
EgammaSuperClusters::hist_EE_RawSC_S25toE_
MonitorElement * hist_EE_RawSC_S25toE_
Definition: EgammaSuperClusters.h:99
edm::HepMCProduct::GetEvent
const HepMC::GenEvent * GetEvent() const
Definition: HepMCProduct.h:34
EgammaSuperClusters::hist_EB_CorSC_deltaR_
MonitorElement * hist_EB_CorSC_deltaR_
Definition: EgammaSuperClusters.h:110
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:123
EgammaSuperClusters::hist_EE_PreSC_S25toE_
MonitorElement * hist_EE_PreSC_S25toE_
Definition: EgammaSuperClusters.h:102
EgammaSuperClusters::hist_EE_RawSC_Eta_
MonitorElement * hist_EE_RawSC_Eta_
Definition: EgammaSuperClusters.h:84
EgammaSuperClusters::ecalClusterToolsESGetTokens_
EcalClusterLazyTools::ESGetTokens ecalClusterToolsESGetTokens_
Definition: EgammaSuperClusters.h:52
EgammaSuperClusters::hspreshowerE_
HistSpec hspreshowerE_
Definition: EgammaSuperClusters.h:65
EcalClusterLazyToolsBase::ESGetTokens::get
ESData get(edm::EventSetup const &eventSetup) const
Definition: EcalClusterLazyTools.h:62
EgammaSuperClusters::hist_EB_CorSC_Eta_
MonitorElement * hist_EB_CorSC_Eta_
Definition: EgammaSuperClusters.h:85
edm::Ptr< CaloCluster >
HistSpec::min
double min
Definition: HistSpec.h:7
EgammaSuperClusters::hsR_
HistSpec hsR_
Definition: EgammaSuperClusters.h:66
EgammaSuperClusters::hist_EE_CorSC_S1toS9_
MonitorElement * hist_EE_CorSC_S1toS9_
Definition: EgammaSuperClusters.h:96
EgammaSuperClusters::hist_EE_PreSC_deltaR_
MonitorElement * hist_EE_PreSC_deltaR_
Definition: EgammaSuperClusters.h:112
EgammaSuperClusters::hist_EB_CorSC_S25toE_
MonitorElement * hist_EB_CorSC_S25toE_
Definition: EgammaSuperClusters.h:100
EgammaSuperClusters::hist_EE_PreSC_preshowerE_
MonitorElement * hist_EE_PreSC_preshowerE_
Definition: EgammaSuperClusters.h:113
EgammaSuperClusters::hist_EE_CorSC_S25toE_
MonitorElement * hist_EE_CorSC_S25toE_
Definition: EgammaSuperClusters.h:101
extraflags_cff.vtx
vtx
Definition: extraflags_cff.py:18
genWeightsTable_cfi.genEvent
genEvent
Definition: genWeightsTable_cfi.py:4
HLT_FULL_cff.dEta
dEta
Definition: HLT_FULL_cff.py:13701
dqm::implementation::IBooker::book2D
MonitorElement * book2D(TString const &name, TString const &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, FUNC onbooking=NOOP())
Definition: DQMStore.h:177
EgammaSuperClusters::barrelRecHitCollectionToken_
edm::EDGetTokenT< EcalRecHitCollection > barrelRecHitCollectionToken_
Definition: EgammaSuperClusters.h:49
EgammaSuperClusters::hist_EE_CorSC_etaWidth_
MonitorElement * hist_EE_CorSC_etaWidth_
Definition: EgammaSuperClusters.h:117
ETA
#define ETA
Definition: GenericBenchmark.cc:28
EgammaSuperClusters::hsET_
HistSpec hsET_
Definition: EgammaSuperClusters.h:56
EgammaSuperClusters::hist_EB_CorSC_ET_
MonitorElement * hist_EB_CorSC_ET_
Definition: EgammaSuperClusters.h:80
EgammaSuperClusters::hist_EB_RawSC_EoverTruth_
MonitorElement * hist_EB_RawSC_EoverTruth_
Definition: EgammaSuperClusters.h:103
EgammaSuperClusters::hist_EE_CorSC_ET_vs_Eta_
MonitorElement * hist_EE_CorSC_ET_vs_Eta_
Definition: EgammaSuperClusters.h:121
EgammaSuperClusters::hist_EE_CorSC_Eta_
MonitorElement * hist_EE_CorSC_Eta_
Definition: EgammaSuperClusters.h:86
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
EgammaSuperClusters::hsPhi_
HistSpec hsPhi_
Definition: EgammaSuperClusters.h:58
dqm-mbProfile.log
log
Definition: dqm-mbProfile.py:17
EgammaSuperClusters::hist_EE_RawSC_Phi_
MonitorElement * hist_EE_RawSC_Phi_
Definition: EgammaSuperClusters.h:89
EgammaSuperClusters::hist_EB_RawSC_S25toE_
MonitorElement * hist_EB_RawSC_S25toE_
Definition: EgammaSuperClusters.h:98
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:29
EgammaSuperClusters::endcapCorSuperClusterCollectionToken_
edm::EDGetTokenT< reco::SuperClusterCollection > endcapCorSuperClusterCollectionToken_
Definition: EgammaSuperClusters.h:46
EgammaSuperClusters::hist_EE_RawSC_NumBC_
MonitorElement * hist_EE_RawSC_NumBC_
Definition: EgammaSuperClusters.h:74
EgammaSuperClusters::hist_EB_RawSC_Size_
MonitorElement * hist_EB_RawSC_Size_
Definition: EgammaSuperClusters.h:68
EgammaSuperClusters::hist_EB_CorSC_ET_vs_Phi_
MonitorElement * hist_EB_CorSC_ET_vs_Phi_
Definition: EgammaSuperClusters.h:120
ZEcal
static constexpr float ZEcal
Definition: L1TkEmParticleProducer.cc:40
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
edm::InputTag
Definition: InputTag.h:15
EgammaSuperClusters::hist_EE_CorSC_ET_vs_R_
MonitorElement * hist_EE_CorSC_ET_vs_R_
Definition: EgammaSuperClusters.h:123
Z_Endcap
static constexpr float Z_Endcap
Definition: ECALPositionCalculator.cc:11
EgammaSuperClusters::hsS25toE_
HistSpec hsS25toE_
Definition: EgammaSuperClusters.h:60
EgammaSuperClusters::hist_EE_RawSC_Size_
MonitorElement * hist_EE_RawSC_Size_
Definition: EgammaSuperClusters.h:69
R_ECAL
static constexpr float R_ECAL
Definition: ECALPositionCalculator.cc:10
EgammaSuperClusters::hsEta_
HistSpec hsEta_
Definition: EgammaSuperClusters.h:57
dqm::implementation::IBooker::book1D
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX, FUNC onbooking=NOOP())
Definition: DQMStore.h:98
EgammaSuperClusters::hist_EB_CorSC_S1toS9_
MonitorElement * hist_EB_CorSC_S1toS9_
Definition: EgammaSuperClusters.h:95
etaBarrelEndcap
static constexpr float etaBarrelEndcap
Definition: ECALPositionCalculator.cc:12
EgammaSuperClusters::hist_EB_CorSC_Size_
MonitorElement * hist_EB_CorSC_Size_
Definition: EgammaSuperClusters.h:70
EgammaSuperClusters::hist_EE_PreSC_NumBC_
MonitorElement * hist_EE_PreSC_NumBC_
Definition: EgammaSuperClusters.h:77
Generated for CMSSW Reference Manual by doxygen 1.8.16