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

#include <Validation/HcalDigis/src/HcalDigisValidation.cc>

Inheritance diagram for HcalDigisValidation:
DQMEDAnalyzer edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > > edm::stream::EDAnalyzerBase edm::EDConsumerBase

Classes

struct  HistLim
 

Public Member Functions

virtual void bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &)
 
virtual void dqmBeginRun (const edm::Run &run, const edm::EventSetup &c)
 
 HcalDigisValidation (const edm::ParameterSet &)
 
 ~HcalDigisValidation ()
 
- Public Member Functions inherited from DQMEDAnalyzer
virtual void beginRun (edm::Run const &, edm::EventSetup const &) final
 
virtual void beginStream (edm::StreamID id) final
 
virtual void dqmBeginRun (edm::Run const &, edm::EventSetup const &)
 
 DQMEDAnalyzer (void)
 
virtual void endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
 
virtual void endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
 
uint32_t streamId () const
 
- Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
 EDAnalyzer ()=default
 
- Public Member Functions inherited from edm::stream::EDAnalyzerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzerBase ()
 
ModuleDescription const & moduleDescription () const
 
virtual ~EDAnalyzerBase ()
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBase & operator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
virtual ~EDConsumerBase () noexcept(false)
 

Private Member Functions

virtual void analyze (const edm::Event &, const edm::EventSetup &)
 
void book1D (DQMStore::IBooker &ib, std::string name, int n, double min, double max)
 
void book1D (DQMStore::IBooker &ib, std::string name, const HistLim &limX)
 
void book2D (DQMStore::IBooker &ib, std::string name, const HistLim &limX, const HistLim &limY)
 
void booking (DQMStore::IBooker &ib, std::string subdetopt, int bnoise, int bmc)
 
void bookPf (DQMStore::IBooker &ib, std::string name, const HistLim &limX, const HistLim &limY)
 
void bookPf (DQMStore::IBooker &ib, std::string name, const HistLim &limX, const HistLim &limY, const char *option)
 
void fill1D (std::string name, double X, double weight=1)
 
void fill2D (std::string name, double X, double Y, double weight=1)
 
void fillPf (std::string name, double X, double Y)
 
MonitorElement * monitor (std::string name)
 
template<class Digi >
void reco (const edm::Event &iEvent, const edm::EventSetup &iSetup, const edm::EDGetTokenT< edm::SortedCollection< Digi > > &tok)
 
template<class dataFrameType >
void reco (const edm::Event &iEvent, const edm::EventSetup &iSetup, const edm::EDGetTokenT< HcalDataFrameContainer< dataFrameType > > &tok)
 
std::string str (int x)
 

Private Attributes

edm::ESHandle< HcalDbService > conditions
 
edm::InputTag dataTPsTag_
 
std::string dirName_
 
edm::InputTag emulTPsTag_
 
edm::ESHandle< CaloGeometry > geometry
 
const HcalDDDRecConstants * hcons
 
bool hep17_
 
edm::ESHandle< HcalTopology > htopo
 
const HcalTopology * htopology
 
edm::InputTag inputTag_
 
int maxDepth_ [5]
 
std::string mc_
 
std::string mode_
 
std::map< std::string, MonitorElement * > * msm_
 
int nChannels_ [5]
 
int nevent1
 
int nevent2
 
int nevent3
 
int nevent4
 
int nevtot
 
int noise_
 
std::string outputFile_
 
edm::InputTag QIE10inputTag_
 
edm::InputTag QIE11inputTag_
 
bool skipDataTPs
 
std::string subdet_
 
bool testNumber_
 
edm::EDGetTokenT< HcalTrigPrimDigiCollection > tok_dataTPs_
 
edm::EDGetTokenT< HcalTrigPrimDigiCollection > tok_emulTPs_
 
edm::EDGetTokenT< HBHEDigiCollection > tok_hbhe_
 
edm::EDGetTokenT< HFDigiCollection > tok_hf_
 
edm::EDGetTokenT< HODigiCollection > tok_ho_
 
edm::EDGetTokenT< edm::PCaloHitContainer > tok_mc_
 
edm::EDGetTokenT< QIE10DigiCollection > tok_qie10_hf_
 
edm::EDGetTokenT< QIE11DigiCollection > tok_qie11_hbhe_
 
std::string zside_
 

Additional Inherited Members

- Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
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
 
- Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr< dqmDetails::NoCache > globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)
 
static std::shared_ptr< dqmDetails::NoCache > globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)
 
static void globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *, dqmDetails::NoCache *)
 
static void globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *, dqmDetails::NoCache *)
 
- Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 59 of file HcalDigisValidation.h.

Constructor & Destructor Documentation

HcalDigisValidation::HcalDigisValidation ( const edm::ParameterSet &  iConfig)
explicit

Definition at line 26 of file HcalDigisValidation.cc.

References dataTPsTag_, dirName_, emulTPsTag_, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hep17_, inputTag_, mc_, mode_, msm_, nevent1, nevent2, nevent3, nevent4, nevtot, outputFile_, QIE10inputTag_, QIE11inputTag_, skipDataTPs, AlCaHLTBitMon_QueryRunRegistry::string, subdet_, testNumber_, tok_dataTPs_, tok_emulTPs_, tok_hbhe_, tok_hf_, tok_ho_, tok_mc_, tok_qie10_hf_, and tok_qie11_hbhe_.

26  {
27 
28  using namespace std;
29 
30  subdet_ = iConfig.getUntrackedParameter<std::string > ("subdetector", "all");
31  outputFile_ = iConfig.getUntrackedParameter<std::string > ("outputFile", "");
32 // inputLabel_ = iConfig.getParameter<std::string > ("digiLabel");
33  inputTag_ = iConfig.getParameter<edm::InputTag > ("digiTag");
34  QIE10inputTag_ = iConfig.getParameter<edm::InputTag > ("QIE10digiTag");
35  QIE11inputTag_ = iConfig.getParameter<edm::InputTag > ("QIE11digiTag");
36  emulTPsTag_ = iConfig.getParameter<edm::InputTag > ("emulTPs");
37  dataTPsTag_ = iConfig.getParameter<edm::InputTag > ("dataTPs");
38  mc_ = iConfig.getUntrackedParameter<std::string > ("mc", "no");
39  mode_ = iConfig.getUntrackedParameter<std::string > ("mode", "multi");
40  dirName_ = iConfig.getUntrackedParameter<std::string > ("dirName", "HcalDigisV/HcalDigiTask");
41  testNumber_= iConfig.getParameter<bool>("TestNumber");
42  hep17_ = iConfig.getParameter<bool>("hep17");
43 
44  // register for data access
45  if (iConfig.exists("simHits"))
46  {
47  tok_mc_ = consumes<edm::PCaloHitContainer>(iConfig.getUntrackedParameter<edm::InputTag>("simHits"));
48  }
49  tok_hbhe_ = consumes< HBHEDigiCollection >(inputTag_);
50  tok_ho_ = consumes< HODigiCollection >(inputTag_);
51  tok_hf_ = consumes< HFDigiCollection >(inputTag_);
52  tok_emulTPs_ = consumes<HcalTrigPrimDigiCollection>(emulTPsTag_);
53  if(dataTPsTag_==edm::InputTag("")) skipDataTPs = true;
54  else {
55  skipDataTPs = false;
56  tok_dataTPs_ = consumes<HcalTrigPrimDigiCollection>(dataTPsTag_);
57  }
58 
59  tok_qie10_hf_ = consumes< QIE10DigiCollection >(QIE10inputTag_);
60  tok_qie11_hbhe_ = consumes< QIE11DigiCollection >(QIE11inputTag_);
61 
62  nevent1 = 0;
63  nevent2 = 0;
64  nevent3 = 0;
65  nevent4 = 0;
66  nevtot = 0;
67 
68  msm_ = new std::map<std::string, MonitorElement*>();
69 
70  if (outputFile_.size() != 0) edm::LogInfo("OutputInfo") << " Hcal Digi Task histograms will be saved to '" << outputFile_.c_str() << "'";
71  else edm::LogInfo("OutputInfo") << " Hcal Digi Task histograms will NOT be saved";
72 
73 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
edm::ParameterSet::exists
bool exists(std::string const &parameterName) const
checks if a parameter exists
Definition: ParameterSet.cc:760
HcalDigisValidation::tok_qie11_hbhe_
edm::EDGetTokenT< QIE11DigiCollection > tok_qie11_hbhe_
Definition: HcalDigisValidation.h:133
HcalDigisValidation::tok_qie10_hf_
edm::EDGetTokenT< QIE10DigiCollection > tok_qie10_hf_
Definition: HcalDigisValidation.h:132
HcalDigisValidation::msm_
std::map< std::string, MonitorElement * > * msm_
Definition: HcalDigisValidation.h:82
HcalDigisValidation::tok_hbhe_
edm::EDGetTokenT< HBHEDigiCollection > tok_hbhe_
Definition: HcalDigisValidation.h:126
HcalDigisValidation::QIE10inputTag_
edm::InputTag QIE10inputTag_
Definition: HcalDigisValidation.h:115
HcalDigisValidation::tok_hf_
edm::EDGetTokenT< HFDigiCollection > tok_hf_
Definition: HcalDigisValidation.h:128
HcalDigisValidation::tok_mc_
edm::EDGetTokenT< edm::PCaloHitContainer > tok_mc_
Definition: HcalDigisValidation.h:125
HcalDigisValidation::tok_ho_
edm::EDGetTokenT< HODigiCollection > tok_ho_
Definition: HcalDigisValidation.h:127
HcalDigisValidation::tok_emulTPs_
edm::EDGetTokenT< HcalTrigPrimDigiCollection > tok_emulTPs_
Definition: HcalDigisValidation.h:129
HcalDigisValidation::QIE11inputTag_
edm::InputTag QIE11inputTag_
Definition: HcalDigisValidation.h:116
HcalDigisValidation::tok_dataTPs_
edm::EDGetTokenT< HcalTrigPrimDigiCollection > tok_dataTPs_
Definition: HcalDigisValidation.h:130
HcalDigisValidation::~HcalDigisValidation ( )

Definition at line 76 of file HcalDigisValidation.cc.

References msm_.

76  {
77  delete msm_;
78 }
HcalDigisValidation::msm_
std::map< std::string, MonitorElement * > * msm_
Definition: HcalDigisValidation.h:82

Member Function Documentation

void HcalDigisValidation::analyze ( const edm::Event &  iEvent,
const edm::EventSetup &  iSetup 
)
privatevirtual

Definition at line 391 of file HcalDigisValidation.cc.

References funct::abs(), edm::SortedCollection< T, SORT >::begin(), EnergyCorrector::c, conditions, HcalDigisParam_cfi::dataTPs, HcalDigisParam_cfi::emulTPs, edm::SortedCollection< T, SORT >::end(), fill1D(), fill2D(), fillPf(), HcalTrigTowerGeometry::firstHFTower(), edm::EventSetup::get(), edm::Event::getByToken(), HcalBarrel, HcalEndcap, CaloTPGTranscoder::hcaletValue(), HcalForward, htopo, iEvent, nevtot, noise_, skipDataTPs, subdet_, tok_dataTPs_, tok_emulTPs_, tok_hbhe_, tok_hf_, tok_ho_, tok_qie10_hf_, and tok_qie11_hbhe_.

391  {
392  using namespace edm;
393  using namespace std;
394 
395  iSetup.get<HcalDbRecord > ().get(conditions);
396 
397  //TP Code
398  ESHandle<CaloTPGTranscoder> decoder;
399  iSetup.get<CaloTPGRecord>().get(decoder);
400 
401  ESHandle<HcalTrigTowerGeometry> tp_geometry;
402  iSetup.get<CaloGeometryRecord>().get(tp_geometry);
403 
404  iSetup.get<HcalRecNumberingRecord>().get(htopo);
405 
406  //Get all handles
407  edm::Handle<HcalTrigPrimDigiCollection> emulTPs;
408  iEvent.getByToken(tok_emulTPs_, emulTPs);
409 
410  edm::Handle<HcalTrigPrimDigiCollection> dataTPs;
411  if(!skipDataTPs) iEvent.getByToken(tok_dataTPs_, dataTPs);
412  //iEvent.getByLabel("hcalDigis", dataTPs);
413 
414  //~TP Code
415 
416  if (subdet_ != "all") {
417  noise_ = 0;
418  if (subdet_ == "HB") reco<HBHEDataFrame > (iEvent, iSetup, tok_hbhe_);
419  if (subdet_ == "HE"){
420  reco<HBHEDataFrame > (iEvent, iSetup, tok_hbhe_);
421  reco<QIE11DataFrame>(iEvent, iSetup, tok_qie11_hbhe_);
422  }
423  if (subdet_ == "HO") reco<HODataFrame > (iEvent, iSetup, tok_ho_);
424  if (subdet_ == "HF"){
425  reco<HFDataFrame > (iEvent, iSetup, tok_hf_);
426  reco<QIE10DataFrame>(iEvent, iSetup, tok_qie10_hf_);
427  }
428 
429  if (subdet_ == "noise") {
430  noise_ = 1;
431  subdet_ = "HB";
432  reco<HBHEDataFrame > (iEvent, iSetup, tok_hbhe_);
433  subdet_ = "HE";
434  reco<HBHEDataFrame > (iEvent, iSetup, tok_hbhe_);
435  reco<QIE11DataFrame>(iEvent, iSetup, tok_qie11_hbhe_);
436  subdet_ = "HO";
437  reco<HODataFrame > (iEvent, iSetup, tok_ho_);
438  subdet_ = "HF";
439  reco<HFDataFrame > (iEvent, iSetup, tok_hf_);
440  reco<QIE10DataFrame>(iEvent, iSetup, tok_qie10_hf_);
441  subdet_ = "noise";
442  }
443  }// all subdetectors
444  else {
445  noise_ = 0;
446 
447  subdet_ = "HB";
448  reco<HBHEDataFrame > (iEvent, iSetup, tok_hbhe_);
449  subdet_ = "HE";
450  reco<HBHEDataFrame > (iEvent, iSetup, tok_hbhe_);
451  reco<QIE11DataFrame>(iEvent, iSetup, tok_qie11_hbhe_);
452  subdet_ = "HO";
453  reco<HODataFrame > (iEvent, iSetup, tok_ho_);
454  subdet_ = "HF";
455  reco<HFDataFrame > (iEvent, iSetup, tok_hf_);
456  reco<QIE10DataFrame>(iEvent, iSetup, tok_qie10_hf_);
457  subdet_ = "all";
458  }
459 
460  fill1D("nevtot", 0);
461  nevtot++;
462 
463  //TP Code
464  //Counters
465  int c = 0, cv0 = 0, cv1 = 0, chb = 0, che = 0, chf = 0, chfv0 = 0, chfv1 = 0;
466 
467  if(skipDataTPs) return;
468 
469  for (HcalTrigPrimDigiCollection::const_iterator itr = dataTPs->begin(); itr != dataTPs->end(); ++itr) {
470  int ieta = itr->id().ieta();
471  int iphi = itr->id().iphi();
472 
473  HcalSubdetector subdet = (HcalSubdetector) 0;
474 
475 
476  if ( abs(ieta) <= 16 )
477  subdet = HcalSubdetector::HcalBarrel ;
478  else if ( abs(ieta) < tp_geometry->firstHFTower(itr->id().version()) )
479  subdet = HcalSubdetector::HcalEndcap ;
480  else if ( abs(ieta) <= 42 )
481  subdet = HcalSubdetector::HcalForward;
482 
483  //Right now, the only case where version matters is in HF
484  //If the subdetector is not HF, set version to -1
485  int tpVersion = (subdet == HcalSubdetector::HcalForward ? itr->id().version() : -1);
486 
487  float en = decoder->hcaletValue(itr->id(), itr->t0());
488 
489  if (en < 0.00001) continue;
490 
491  //Plot the variables
492  //Retain classic behavior (include all tps)
493  //Additional plots that only include HF 3x2 or HF 1x1
494 
495  //Classic
496  fill1D("HcalDigiTask_tp_et",en);
497  fill2D("HcalDigiTask_tp_et_ieta",ieta,en);
498  fill2D("HcalDigiTask_tp_ieta_iphi",ieta,iphi);
499  fillPf("HcalDigiTask_tp_ave_et_ieta",ieta,en);
500  fill1D("HcalDigiTask_tp_ntp_ieta",ieta);
501  fill1D("HcalDigiTask_tp_ntp_iphi",iphi);
502  if ( en > 10. ) fill1D("HcalDigiTask_tp_ntp_10_ieta",ieta);
503 
504  //3x2 Trig Primitives (tpVersion == 0)
505  if( subdet != HcalSubdetector::HcalForward || tpVersion==0){
506  fill1D("HcalDigiTask_tp_et_v0",en);
507  fill2D("HcalDigiTask_tp_et_ieta_v0",ieta,en);
508  fill2D("HcalDigiTask_tp_ieta_iphi_v0",ieta,iphi);
509  fillPf("HcalDigiTask_tp_ave_et_ieta_v0",ieta,en);
510  fill1D("HcalDigiTask_tp_ntp_ieta_v0",ieta);
511  fill1D("HcalDigiTask_tp_ntp_iphi_v0",iphi);
512  if ( en > 10. ) fill1D("HcalDigiTask_tp_ntp_10_ieta_v0",ieta);
513  }
514 
515  //1x1 Trig Primitives (tpVersion == 1)
516  if( subdet != HcalSubdetector::HcalForward || tpVersion==1){
517  fill1D("HcalDigiTask_tp_et_v1",en);
518  fill2D("HcalDigiTask_tp_et_ieta_v1",ieta,en);
519  fill2D("HcalDigiTask_tp_ieta_iphi_v1",ieta,iphi);
520  fillPf("HcalDigiTask_tp_ave_et_ieta_v1",ieta,en);
521  fill1D("HcalDigiTask_tp_ntp_ieta_v1",ieta);
522  fill1D("HcalDigiTask_tp_ntp_iphi_v1",iphi);
523  if ( en > 10. ) fill1D("HcalDigiTask_tp_ntp_10_ieta_v1",ieta);
524  }
525 
526  ++c;
527  if ( subdet == HcalSubdetector::HcalBarrel ) {
528  fill1D("HcalDigiTask_tp_et_HB",en);
529  ++chb;
530  ++cv0;
531  ++cv1;
532  }
533  if ( subdet == HcalSubdetector::HcalEndcap ) {
534  fill1D("HcalDigiTask_tp_et_HE",en);
535  ++che;
536  ++cv0;
537  ++cv1;
538  }
539  if ( subdet == HcalSubdetector::HcalForward ) {
540  fill1D("HcalDigiTask_tp_et_HF",en);
541  ++chf;
542 
543  if(tpVersion == 0){
544  fill1D("HcalDigiTask_tp_et_HF_v0",en);
545  ++chfv0;
546  ++cv0;
547  }
548 
549  if(tpVersion == 1){
550  fill1D("HcalDigiTask_tp_et_HF_v1",en);
551  ++chfv1;
552  ++cv1;
553  }
554 
555  }
556 
557  }//end data TP collection
558 
559  fill1D("HcalDigiTask_tp_ntp",c);
560  fill1D("HcalDigiTask_tp_ntp_v0",cv0);
561  fill1D("HcalDigiTask_tp_ntp_v1",cv1);
562  fill1D("HcalDigiTask_tp_ntp_HB",chb);
563  fill1D("HcalDigiTask_tp_ntp_HE",che);
564  fill1D("HcalDigiTask_tp_ntp_HF",chf);
565  fill1D("HcalDigiTask_tp_ntp_HF_v0",chfv0);
566  fill1D("HcalDigiTask_tp_ntp_HF_v1",chfv1);
567 
568  //~TP Code
569 }
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:460
edm::SortedCollection< HcalTriggerPrimitiveDigi >::const_iterator
std::vector< HcalTriggerPrimitiveDigi >::const_iterator const_iterator
Definition: SortedCollection.h:82
HcalDigisValidation::tok_qie11_hbhe_
edm::EDGetTokenT< QIE11DigiCollection > tok_qie11_hbhe_
Definition: HcalDigisValidation.h:133
HcalDigisValidation::tok_qie10_hf_
edm::EDGetTokenT< QIE10DigiCollection > tok_qie10_hf_
Definition: HcalDigisValidation.h:132
HcalDigisParam_cfi.emulTPs
emulTPs
Definition: HcalDigisParam_cfi.py:13
iEvent
int iEvent
Definition: GenABIO.cc:230
HcalDigisValidation::tok_hbhe_
edm::EDGetTokenT< HBHEDigiCollection > tok_hbhe_
Definition: HcalDigisValidation.h:126
HcalDigisValidation::fillPf
void fillPf(std::string name, double X, double Y)
Definition: HcalDigisValidation.cc:1255
HcalDigisValidation::fill2D
void fill2D(std::string name, double X, double Y, double weight=1)
Definition: HcalDigisValidation.cc:1243
HcalDigisValidation::htopo
edm::ESHandle< HcalTopology > htopo
Definition: HcalDigisValidation.h:140
HcalSubdetector
HcalSubdetector
Definition: HcalAssistant.h:31
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
edm::SortedCollection::end
const_iterator end() const
Definition: SortedCollection.h:310
HcalDigisValidation::tok_hf_
edm::EDGetTokenT< HFDigiCollection > tok_hf_
Definition: HcalDigisValidation.h:128
HcalDigisValidation::fill1D
void fill1D(std::string name, double X, double weight=1)
Definition: HcalDigisValidation.cc:1235
HcalDigisValidation::tok_ho_
edm::EDGetTokenT< HODigiCollection > tok_ho_
Definition: HcalDigisValidation.h:127
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:56
HcalDigisValidation::tok_emulTPs_
edm::EDGetTokenT< HcalTrigPrimDigiCollection > tok_emulTPs_
Definition: HcalDigisValidation.h:129
edm
HLT enums.
Definition: AlignableModifier.h:17
HcalDigisValidation::tok_dataTPs_
edm::EDGetTokenT< HcalTrigPrimDigiCollection > tok_dataTPs_
Definition: HcalDigisValidation.h:130
HcalDigisValidation::conditions
edm::ESHandle< HcalDbService > conditions
Definition: HcalDigisValidation.h:137
HcalDigisParam_cfi.dataTPs
dataTPs
Definition: HcalDigisParam_cfi.py:14
CaloTPGTranscoder::hcaletValue
virtual double hcaletValue(const int &ieta, const int &iphi, const int &version, const int &compressedValue) const =0
edm::SortedCollection::begin
const_iterator begin() const
Definition: SortedCollection.h:303
HcalTrigTowerGeometry::firstHFTower
int firstHFTower(int version) const
Definition: HcalTrigTowerGeometry.h:19
void HcalDigisValidation::book1D ( DQMStore::IBooker &  ib,
std::string  name,
int  n,
double  min,
double  max 
)
private

Definition at line 1227 of file HcalDigisValidation.cc.

References DQMStore::IBooker::book1D(), hpstanc_transforms::max, min(), msm_, gen::n, and dataset::name.

Referenced by bookHistograms(), and booking().

1227  {
1228  if (!msm_->count(name)) (*msm_)[name] = ib.book1D(name.c_str(), name.c_str(), n, min, max);
1229 }
hpstanc_transforms.max
max
Definition: hpstanc_transforms.py:5
HcalDigisValidation::msm_
std::map< std::string, MonitorElement * > * msm_
Definition: HcalDigisValidation.h:82
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
min
T min(T a, T b)
Definition: MathUtil.h:58
void HcalDigisValidation::book1D ( DQMStore::IBooker &  ib,
std::string  name,
const HistLim &  limX 
)
private

Definition at line 1231 of file HcalDigisValidation.cc.

References DQMStore::IBooker::book1D(), HcalDigisValidation::HistLim::max, HcalDigisValidation::HistLim::min, msm_, HcalDigisValidation::HistLim::n, and dataset::name.

1231  {
1232  if (!msm_->count(name)) (*msm_)[name] = ib.book1D(name.c_str(), name.c_str(), limX.n, limX.min, limX.max);
1233 }
HcalDigisValidation::msm_
std::map< std::string, MonitorElement * > * msm_
Definition: HcalDigisValidation.h:82
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
void HcalDigisValidation::book2D ( DQMStore::IBooker &  ib,
std::string  name,
const HistLim &  limX,
const HistLim &  limY 
)
private

Definition at line 1239 of file HcalDigisValidation.cc.

References DQMStore::IBooker::book2D(), HcalDigisValidation::HistLim::max, HcalDigisValidation::HistLim::min, msm_, HcalDigisValidation::HistLim::n, and dataset::name.

Referenced by bookHistograms(), and booking().

1239  {
1240  if (!msm_->count(name)) (*msm_)[name] = ib.book2D(name.c_str(), name.c_str(), limX.n, limX.min, limX.max, limY.n, limY.min, limY.max);
1241 }
HcalDigisValidation::msm_
std::map< std::string, MonitorElement * > * msm_
Definition: HcalDigisValidation.h:82
DQMStore::IBooker::book2D
MonitorElement * book2D(Args &&...args)
Definition: DQMStore.h:133
void HcalDigisValidation::bookHistograms ( DQMStore::IBooker &  ib,
edm::Run const &  run,
edm::EventSetup const &  es 
)
virtual

Implements DQMEDAnalyzer.

Definition at line 113 of file HcalDigisValidation.cc.

References book1D(), book2D(), booking(), bookPf(), dirName_, mc_, DQMStore::IBooker::setCurrentFolder(), skipDataTPs, and subdet_.

114 {
115 
116  ib.setCurrentFolder(dirName_);
117 
118  // book
119  book1D(ib,"nevtot", 1, 0, 1);
120  int bnoise = 0;
121  int bmc = 0;
122  if (subdet_ == "noise") bnoise = 1;
123  if (mc_ == "yes") bmc = 1;
124  if (subdet_ == "noise" || subdet_ == "all") {
125  booking(ib,"HB", bnoise, bmc);
126  booking(ib,"HO", bnoise, bmc);
127  booking(ib,"HF", bnoise, bmc);
128  booking(ib,"HE", bnoise, bmc);
129  } else {
130  booking(ib,subdet_, 0, bmc);
131  }
132 
133  if(skipDataTPs) return;
134 
135  HistLim tp_hl_et(260, -10, 250);
136  HistLim tp_hl_ntp(640, -20, 3180);
137  HistLim tp_hl_ntp_sub(404, -20, 2000);
138  HistLim tp_hl_ieta(85, -42.5, 42.5);
139  HistLim tp_hl_iphi(74,-0.5,73.5);
140 
141 
142  book1D(ib,"HcalDigiTask_tp_et", tp_hl_et);
143  book1D(ib,"HcalDigiTask_tp_et_v0", tp_hl_et);
144  book1D(ib,"HcalDigiTask_tp_et_v1", tp_hl_et);
145  book1D(ib,"HcalDigiTask_tp_et_HB", tp_hl_et);
146  book1D(ib,"HcalDigiTask_tp_et_HE", tp_hl_et);
147  book1D(ib,"HcalDigiTask_tp_et_HF", tp_hl_et);
148  book1D(ib,"HcalDigiTask_tp_et_HF_v0", tp_hl_et);
149  book1D(ib,"HcalDigiTask_tp_et_HF_v1", tp_hl_et);
150  book1D(ib,"HcalDigiTask_tp_ntp", tp_hl_ntp);
151  book1D(ib,"HcalDigiTask_tp_ntp_v0", tp_hl_ntp);
152  book1D(ib,"HcalDigiTask_tp_ntp_v1", tp_hl_ntp);
153  book1D(ib,"HcalDigiTask_tp_ntp_HB", tp_hl_ntp_sub);
154  book1D(ib,"HcalDigiTask_tp_ntp_HE", tp_hl_ntp_sub);
155  book1D(ib,"HcalDigiTask_tp_ntp_HF", tp_hl_ntp_sub);
156  book1D(ib,"HcalDigiTask_tp_ntp_HF_v0", tp_hl_ntp_sub);
157  book1D(ib,"HcalDigiTask_tp_ntp_HF_v1", tp_hl_ntp_sub);
158  book1D(ib,"HcalDigiTask_tp_ntp_ieta", tp_hl_ieta);
159  book1D(ib,"HcalDigiTask_tp_ntp_ieta_v0", tp_hl_ieta);
160  book1D(ib,"HcalDigiTask_tp_ntp_ieta_v1", tp_hl_ieta);
161  book1D(ib,"HcalDigiTask_tp_ntp_iphi", tp_hl_iphi);
162  book1D(ib,"HcalDigiTask_tp_ntp_iphi_v0", tp_hl_iphi);
163  book1D(ib,"HcalDigiTask_tp_ntp_iphi_v1", tp_hl_iphi);
164  book1D(ib,"HcalDigiTask_tp_ntp_10_ieta", tp_hl_ieta);
165  book1D(ib,"HcalDigiTask_tp_ntp_10_ieta_v0", tp_hl_ieta);
166  book1D(ib,"HcalDigiTask_tp_ntp_10_ieta_v1", tp_hl_ieta);
167  book2D(ib,"HcalDigiTask_tp_et_ieta", tp_hl_ieta, tp_hl_et);
168  book2D(ib,"HcalDigiTask_tp_et_ieta_v0", tp_hl_ieta, tp_hl_et);
169  book2D(ib,"HcalDigiTask_tp_et_ieta_v1", tp_hl_ieta, tp_hl_et);
170  book2D(ib,"HcalDigiTask_tp_ieta_iphi", tp_hl_ieta, tp_hl_iphi);
171  book2D(ib,"HcalDigiTask_tp_ieta_iphi_v0", tp_hl_ieta, tp_hl_iphi);
172  book2D(ib,"HcalDigiTask_tp_ieta_iphi_v1", tp_hl_ieta, tp_hl_iphi);
173  bookPf(ib,"HcalDigiTask_tp_ave_et_ieta", tp_hl_ieta, tp_hl_et, " ");
174  bookPf(ib,"HcalDigiTask_tp_ave_et_ieta_v0", tp_hl_ieta, tp_hl_et, " ");
175  bookPf(ib,"HcalDigiTask_tp_ave_et_ieta_v1", tp_hl_ieta, tp_hl_et, " ");
176 
177 }
HcalDigisValidation::bookPf
void bookPf(DQMStore::IBooker &ib, std::string name, const HistLim &limX, const HistLim &limY)
Definition: HcalDigisValidation.cc:1247
HcalDigisValidation::book2D
void book2D(DQMStore::IBooker &ib, std::string name, const HistLim &limX, const HistLim &limY)
Definition: HcalDigisValidation.cc:1239
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:277
HcalDigisValidation::booking
void booking(DQMStore::IBooker &ib, std::string subdetopt, int bnoise, int bmc)
Definition: HcalDigisValidation.cc:179
HcalDigisValidation::book1D
void book1D(DQMStore::IBooker &ib, std::string name, int n, double min, double max)
Definition: HcalDigisValidation.cc:1227
void HcalDigisValidation::booking ( DQMStore::IBooker &  ib,
std::string  subdetopt,
int  bnoise,
int  bmc 
)
private

Definition at line 179 of file HcalDigisValidation.cc.

References book1D(), book2D(), bookPf(), particleFlowClusterECALTimeSelected_cfi::depth, cropTnPTrees::frac, hep17_, trackerHits::histo, maxDepth_, nChannels_, muonCSCDigis_cfi::pedestal, and particleFlowDisplacedVertex_cfi::ratio.

Referenced by bookHistograms().

179  {
180 
181  // Adjust/Optimize binning (JR Dittmann, 16-JUL-2015)
182 
183  HistLim Ndigis(2600, 0., 2600.);
184  HistLim ndigis(520, -20., 1020.);
185  HistLim sime(200, 0., 1.0);
186 
187  HistLim digiAmp(360, -100., 7100.);
188  HistLim digiAmpWide(360, -10000., 710000.);
189  HistLim ratio(2000, -100., 3900.);
190  HistLim sumAmp(100, -500., 1500.);
191 
192  HistLim nbin(11, -0.5, 10.5);
193 
194  HistLim pedestal(80, -1.0, 15.);
195  HistLim pedestalfC(400, -10., 30.);
196 
197  HistLim frac(80, -0.20, 1.40);
198 
199  HistLim pedLim(80, 0., 8.);
200  HistLim pedWidthLim(100, 0., 2.);
201 
202  HistLim gainLim(120, 0., 0.6);
203  HistLim gainWidthLim(160, 0., 0.32);
204 
205  HistLim ietaLim(85, -42.5, 42.5);
206  HistLim iphiLim(74, -0.5, 73.5);
207 
208  HistLim depthLim(15,-0.5,14.5);
209 
210  if (bsubdet == "HB") {
211  Ndigis = HistLim( ((int)(nChannels_[1]/100) + 1)*100, 0., (float)((int)(nChannels_[1]/100) + 1)*100);
212  } else if (bsubdet == "HE") {
213  sime = HistLim(200, 0., 1.0);
214  Ndigis = HistLim( ((int)(nChannels_[2]/100) + 1)*100, 0., (float)((int)(nChannels_[2]/100) + 1)*100);
215  } else if (bsubdet == "HF") {
216  sime = HistLim(100, 0., 100.);
217  pedLim = HistLim(100, 0., 20.);
218  pedWidthLim = HistLim(100, 0., 5.);
219  frac = HistLim(400, -4.00, 4.00);
220  Ndigis = HistLim( ((int)(nChannels_[4]/100) + 1)*100, 0., (float)((int)(nChannels_[4]/100) + 1)*100);
221  } else if (bsubdet == "HO") {
222  sime = HistLim(200, 0., 1.0);
223  gainLim = HistLim(160, 0., 1.6);
224  Ndigis = HistLim( ((int)(nChannels_[3]/100) + 1)*100, 0., (float)((int)(nChannels_[3]/100) + 1)*100);
225  }
226 
227  int isubdet=0;
228  if (bsubdet=="HB") isubdet=1;
229  else if (bsubdet=="HE") isubdet=2;
230  else if (bsubdet=="HO") isubdet=3;
231  else if (bsubdet=="HF") isubdet=4;
232  else edm::LogWarning("HcalDigisValidation") << "HcalDigisValidation Warning: not HB/HE/HF/HO " << bsubdet << std::endl;
233 
234  Char_t histo[100];
235  const char * sub = bsubdet.c_str();
236  if (bnoise == 0) {
237  // number of digis in each subdetector
238  sprintf(histo, "HcalDigiTask_Ndigis_%s", sub);
239  book1D(ib, histo, Ndigis);
240 
241  // maps of occupancies
242  for (int depth = 1; depth <= maxDepth_[isubdet]; depth++) {
243  sprintf(histo, "HcalDigiTask_ieta_iphi_occupancy_map_depth%d_%s", depth, sub);
244  book2D(ib, histo, ietaLim, iphiLim);
245  }
246 
247  //Depths
248  sprintf(histo, "HcalDigiTask_depths_%s",sub);
249  book1D(ib,histo, depthLim);
250 
251  // occupancies vs ieta
252  for (int depth = 1; depth <= maxDepth_[isubdet]; depth++) {
253  sprintf(histo, "HcalDigiTask_occupancy_vs_ieta_depth%d_%s", depth, sub);
254  book1D(ib, histo, ietaLim);
255  }
256 
257 
258  // just 1D of all cells' amplitudes
259  sprintf(histo, "HcalDigiTask_sum_all_amplitudes_%s", sub);
260  book1D(ib, histo, sumAmp);
261 
262  sprintf(histo, "HcalDigiTask_number_of_amplitudes_above_10fC_%s", sub);
263  book1D(ib, histo, ndigis);
264 
265  for (int depth = 1; depth <= maxDepth_[isubdet]; depth++) {
266  sprintf(histo, "HcalDigiTask_ADC0_adc_depth%d_%s", depth, sub);
267  book1D(ib, histo, pedestal);
268  }
269 
270  for (int depth = 1; depth <= maxDepth_[isubdet]; depth++) {
271  sprintf(histo, "HcalDigiTask_ADC0_fC_depth%d_%s", depth, sub);
272  book1D(ib, histo, pedestalfC);
273  }
274 
275  sprintf(histo, "HcalDigiTask_signal_amplitude_%s", sub);
276  book1D(ib, histo, digiAmp);
277 
278  if(hep17_ && bsubdet=="HE"){
279  sprintf(histo, "HcalDigiTask_signal_amplitude_HEP17");
280  book1D(ib, histo, digiAmpWide);
281  }
282  //
283  for (int depth = 1; depth <= maxDepth_[isubdet]; depth++) {
284  sprintf(histo, "HcalDigiTask_signal_amplitude_depth%d_%s", depth, sub);
285  book1D(ib, histo, digiAmp);
286  if(hep17_ && bsubdet=="HE"){
287  sprintf(histo, "HcalDigiTask_signal_amplitude_depth%d_HEP17", depth);
288  book1D(ib, histo, digiAmpWide);
289  }
290  }
291 
292  sprintf(histo, "HcalDigiTask_signal_amplitude_vs_bin_all_depths_%s", sub);
293  book2D(ib, histo, nbin, digiAmp);
294  if(hep17_ && bsubdet=="HE"){
295  sprintf(histo, "HcalDigiTask_signal_amplitude_vs_bin_all_depths_HEP17");
296  book2D(ib, histo, nbin, digiAmpWide);
297  }
298 
299  for (int depth = 1; depth <= maxDepth_[isubdet]; depth++) {
300  sprintf(histo, "HcalDigiTask_all_amplitudes_vs_bin_1D_depth%d_%s", depth, sub);
301  book1D(ib, histo, nbin);
302  if(hep17_ && bsubdet=="HE"){
303  sprintf(histo, "HcalDigiTask_all_amplitudes_vs_bin_1D_depth%d_HEP17", depth);
304  book1D(ib, histo, nbin);
305  }
306  }
307 
308  sprintf(histo, "HcalDigiTask_SOI_frac_%s", sub);
309  book1D(ib, histo, frac);
310  sprintf(histo, "HcalDigiTask_postSOI_frac_%s", sub);
311  book1D(ib, histo, frac);
312 
313  if (bmc == 1) {
314  sprintf(histo, "HcalDigiTask_amplitude_vs_simhits_%s", sub);
315  book2D(ib, histo, sime, digiAmp);
316  for (int depth = 1; depth <= maxDepth_[isubdet]; depth++) {
317  sprintf(histo, "HcalDigiTask_amplitude_vs_simhits_depth%d_%s", depth, sub);
318  book2D(ib, histo, sime, digiAmp);
319  }
320 
321  sprintf(histo, "HcalDigiTask_amplitude_vs_simhits_profile_%s", sub);
322  bookPf(ib, histo, sime, digiAmp);
323  for (int depth = 1; depth <= maxDepth_[isubdet]; depth++) {
324  sprintf(histo, "HcalDigiTask_amplitude_vs_simhits_profile_depth%d_%s", depth, sub);
325  bookPf(ib, histo, sime, digiAmp);
326  }
327 
328  sprintf(histo, "HcalDigiTask_ratio_amplitude_vs_simhits_%s", sub);
329  book1D(ib, histo, ratio);
330  for (int depth = 1; depth <= maxDepth_[isubdet]; depth++) {
331  sprintf(histo, "HcalDigiTask_ratio_amplitude_vs_simhits_depth%d_%s", depth, sub);
332  book1D(ib, histo, ratio);
333  }
334  }//mc only
335 
336  } else { // noise only
337 
338  // EVENT "1" distributions of all cells properties
339 
340  //KH
341  for (int depth = 1; depth <= maxDepth_[isubdet]; depth++) {
342  sprintf(histo, "HcalDigiTask_gain_capId0_Depth%d_%s", depth, sub);
343  book1D(ib, histo, gainLim);
344  sprintf(histo, "HcalDigiTask_gain_capId1_Depth%d_%s", depth, sub);
345  book1D(ib, histo, gainLim);
346  sprintf(histo, "HcalDigiTask_gain_capId2_Depth%d_%s", depth, sub);
347  book1D(ib, histo, gainLim);
348  sprintf(histo, "HcalDigiTask_gain_capId3_Depth%d_%s", depth, sub);
349  book1D(ib, histo, gainLim);
350 
351  sprintf(histo, "HcalDigiTask_gainWidth_capId0_Depth%d_%s", depth, sub);
352  book1D(ib, histo, gainWidthLim);
353  sprintf(histo, "HcalDigiTask_gainWidth_capId1_Depth%d_%s", depth, sub);
354  book1D(ib, histo, gainWidthLim);
355  sprintf(histo, "HcalDigiTask_gainWidth_capId2_Depth%d_%s", depth, sub);
356  book1D(ib, histo, gainWidthLim);
357  sprintf(histo, "HcalDigiTask_gainWidth_capId3_Depth%d_%s", depth, sub);
358  book1D(ib, histo, gainWidthLim);
359 
360  sprintf(histo, "HcalDigiTask_pedestal_capId0_Depth%d_%s", depth, sub);
361  book1D(ib, histo, pedLim);
362  sprintf(histo, "HcalDigiTask_pedestal_capId1_Depth%d_%s", depth, sub);
363  book1D(ib, histo, pedLim);
364  sprintf(histo, "HcalDigiTask_pedestal_capId2_Depth%d_%s", depth, sub);
365  book1D(ib, histo, pedLim);
366  sprintf(histo, "HcalDigiTask_pedestal_capId3_Depth%d_%s", depth, sub);
367  book1D(ib, histo, pedLim);
368 
369  sprintf(histo, "HcalDigiTask_pedestal_width_capId0_Depth%d_%s", depth, sub);
370  book1D(ib, histo, pedWidthLim);
371  sprintf(histo, "HcalDigiTask_pedestal_width_capId1_Depth%d_%s", depth, sub);
372  book1D(ib, histo, pedWidthLim);
373  sprintf(histo, "HcalDigiTask_pedestal_width_capId2_Depth%d_%s", depth, sub);
374  book1D(ib, histo, pedWidthLim);
375  sprintf(histo, "HcalDigiTask_pedestal_width_capId3_Depth%d_%s", depth, sub);
376  book1D(ib, histo, pedWidthLim);
377 
378  }
379 
380  //KH
381  for (int depth = 1; depth <= maxDepth_[isubdet]; depth++) {
382  sprintf(histo, "HcalDigiTask_gainMap_Depth%d_%s", depth, sub);
383  book2D(ib, histo, ietaLim, iphiLim);
384  sprintf(histo, "HcalDigiTask_pwidthMap_Depth%d_%s", depth, sub);
385  book2D(ib, histo, ietaLim, iphiLim);
386  }
387 
388  } //end of noise-only
389 }//book
particleFlowDisplacedVertex_cfi.ratio
ratio
Definition: particleFlowDisplacedVertex_cfi.py:93
particleFlowClusterECALTimeSelected_cfi.depth
depth
Definition: particleFlowClusterECALTimeSelected_cfi.py:8
HcalDigisValidation::bookPf
void bookPf(DQMStore::IBooker &ib, std::string name, const HistLim &limX, const HistLim &limY)
Definition: HcalDigisValidation.cc:1247
HcalDigisValidation::book2D
void book2D(DQMStore::IBooker &ib, std::string name, const HistLim &limX, const HistLim &limY)
Definition: HcalDigisValidation.cc:1239
muonCSCDigis_cfi.pedestal
pedestal
Definition: muonCSCDigis_cfi.py:12
HcalDigisValidation::book1D
void book1D(DQMStore::IBooker &ib, std::string name, int n, double min, double max)
Definition: HcalDigisValidation.cc:1227
void HcalDigisValidation::bookPf ( DQMStore::IBooker &  ib,
std::string  name,
const HistLim &  limX,
const HistLim &  limY 
)
private

Definition at line 1247 of file HcalDigisValidation.cc.

References DQMStore::IBooker::bookProfile(), HcalDigisValidation::HistLim::max, HcalDigisValidation::HistLim::min, msm_, HcalDigisValidation::HistLim::n, and dataset::name.

Referenced by bookHistograms(), and booking().

1247  {
1248  if (!msm_->count(name)) (*msm_)[name] = ib.bookProfile(name.c_str(), name.c_str(), limX.n, limX.min, limX.max, limY.n, limY.min, limY.max);
1249 }
DQMStore::IBooker::bookProfile
MonitorElement * bookProfile(Args &&...args)
Definition: DQMStore.h:157
HcalDigisValidation::msm_
std::map< std::string, MonitorElement * > * msm_
Definition: HcalDigisValidation.h:82
void HcalDigisValidation::bookPf ( DQMStore::IBooker &  ib,
std::string  name,
const HistLim &  limX,
const HistLim &  limY,
const char *  option 
)
private

Definition at line 1251 of file HcalDigisValidation.cc.

References DQMStore::IBooker::bookProfile(), HcalDigisValidation::HistLim::max, HcalDigisValidation::HistLim::min, msm_, HcalDigisValidation::HistLim::n, dataset::name, and TSGForRoadSearch_cfi::option.

1251  {
1252  if (!msm_->count(name)) (*msm_)[name] = ib.bookProfile(name.c_str(), name.c_str(), limX.n, limX.min, limX.max, limY.n, limY.min, limY.max, option);
1253 }
DQMStore::IBooker::bookProfile
MonitorElement * bookProfile(Args &&...args)
Definition: DQMStore.h:157
HcalDigisValidation::msm_
std::map< std::string, MonitorElement * > * msm_
Definition: HcalDigisValidation.h:82
TSGForRoadSearch_cfi.option
option
Definition: TSGForRoadSearch_cfi.py:9
void HcalDigisValidation::dqmBeginRun ( const edm::Run &  run,
const edm::EventSetup &  c 
)
virtual

Definition at line 81 of file HcalDigisValidation.cc.

References edm::EventSetup::get(), HcalGeometry::getHxSize(), HcalDDDRecConstants::getMaxDepth(), CaloGeometry::getSubdetectorGeometry(), DetId::Hcal, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hcons, htopology, maxDepth_, and nChannels_.

81  {
82 
83  edm::ESHandle<HcalDDDRecConstants> pHRNDC;
84  es.get<HcalRecNumberingRecord>().get( pHRNDC );
85  hcons = &(*pHRNDC);
86 
87  htopology = new HcalTopology(hcons);
88 
89  maxDepth_[1] = hcons->getMaxDepth(0); // HB
90  maxDepth_[2] = hcons->getMaxDepth(1); // HE
91  maxDepth_[3] = hcons->getMaxDepth(3); // HO
92  maxDepth_[4] = hcons->getMaxDepth(2); // HF
93  maxDepth_[0] = (maxDepth_[1] > maxDepth_[2] ? maxDepth_[1] : maxDepth_[2]);
94  maxDepth_[0] = (maxDepth_[0] > maxDepth_[3] ? maxDepth_[0] : maxDepth_[3]);
95  maxDepth_[0] = (maxDepth_[0] > maxDepth_[4] ? maxDepth_[0] : maxDepth_[4]); // any of HB/HE/HO/HF
96 
97  es.get<CaloGeometryRecord>().get(geometry);
98  const CaloGeometry* geo = geometry.product();
99  const HcalGeometry* gHB = (HcalGeometry*)(geo->getSubdetectorGeometry(DetId::Hcal,HcalBarrel));
100  const HcalGeometry* gHE = (HcalGeometry*)(geo->getSubdetectorGeometry(DetId::Hcal,HcalEndcap));
101  const HcalGeometry* gHO = (HcalGeometry*)(geo->getSubdetectorGeometry(DetId::Hcal,HcalOuter));
102  const HcalGeometry* gHF = (HcalGeometry*)(geo->getSubdetectorGeometry(DetId::Hcal,HcalForward));
103 
104  nChannels_[1] = gHB->getHxSize(1);
105  nChannels_[2] = gHE->getHxSize(2);
106  nChannels_[3] = gHO->getHxSize(3);
107  nChannels_[4] = gHF->getHxSize(4);
108 
109  nChannels_[0] = nChannels_[1] + nChannels_[2] + nChannels_[3] + nChannels_[4];
110 
111 }
CaloGeometry::getSubdetectorGeometry
const CaloSubdetectorGeometry * getSubdetectorGeometry(const DetId &id) const
access the subdetector geometry for the given subdetector directly
Definition: CaloGeometry.cc:45
HcalGeometry::getHxSize
unsigned int getHxSize(const int type) const
Definition: HcalGeometry.cc:253
HcalDigisValidation::htopology
const HcalTopology * htopology
Definition: HcalDigisValidation.h:150
HcalDDDRecConstants::getMaxDepth
int getMaxDepth(const int type) const
Definition: HcalDDDRecConstants.h:84
HcalDigisValidation::hcons
const HcalDDDRecConstants * hcons
Definition: HcalDigisValidation.h:149
void HcalDigisValidation::fill1D ( std::string  name,
double  X,
double  weight = 1 
)
private

Definition at line 1235 of file HcalDigisValidation.cc.

References msm_.

Referenced by analyze(), and reco().

1235  {
1236  msm_->find(name)->second->Fill(X, weight);
1237 }
X
#define X(str)
Definition: MuonsGrabber.cc:48
weight
Definition: weight.py:1
HcalDigisValidation::msm_
std::map< std::string, MonitorElement * > * msm_
Definition: HcalDigisValidation.h:82
void HcalDigisValidation::fill2D ( std::string  name,
double  X,
double  Y,
double  weight = 1 
)
private

Definition at line 1243 of file HcalDigisValidation.cc.

References msm_.

Referenced by analyze(), and reco().

1243  {
1244  msm_->find(name)->second->Fill(X, Y, weight);
1245 }
X
#define X(str)
Definition: MuonsGrabber.cc:48
weight
Definition: weight.py:1
HcalDigisValidation::msm_
std::map< std::string, MonitorElement * > * msm_
Definition: HcalDigisValidation.h:82
void HcalDigisValidation::fillPf ( std::string  name,
double  X,
double  Y 
)
private

Definition at line 1255 of file HcalDigisValidation.cc.

References msm_.

Referenced by analyze(), and reco().

1255  {
1256  msm_->find(name)->second->Fill(X, Y);
1257 }
X
#define X(str)
Definition: MuonsGrabber.cc:48
HcalDigisValidation::msm_
std::map< std::string, MonitorElement * > * msm_
Definition: HcalDigisValidation.h:82
MonitorElement * HcalDigisValidation::monitor ( std::string  name)
private

Definition at line 1259 of file HcalDigisValidation.cc.

References msm_, and NULL.

1259  {
1260  if (!msm_->count(name)) return NULL;
1261  else return msm_->find(name)->second;
1262 }
NULL
#define NULL
Definition: scimark2.h:8
HcalDigisValidation::msm_
std::map< std::string, MonitorElement * > * msm_
Definition: HcalDigisValidation.h:82
template<class Digi >
void HcalDigisValidation::reco ( const edm::Event &  iEvent,
const edm::EventSetup &  iSetup,
const edm::EDGetTokenT< edm::SortedCollection< Digi > > &  tok 
)
private

Definition at line 571 of file HcalDigisValidation.cc.

References HcalCoderDb::adc2fC(), AlignmentProducer_cff::calibrations, conditions, particleFlowClusterECALTimeSelected_cfi::depth, HcalDetId::depth(), hgcalDigitizer_cfi::digiCollection, fill1D(), fill2D(), fillPf(), muonCSCDigis_cfi::gain, edm::Event::getByToken(), HcalDbService::getGain(), HcalDbService::getGainWidth(), HcalDbService::getHcalCalibrations(), HcalDbService::getHcalCoder(), HcalDbService::getHcalShape(), HcalDbService::getPedestal(), HcalGainWidth::getValue(), HcalGain::getValue(), HcalPedestalWidth::getWidth(), hcons, mps_fire::i, HcalDetId::ieta(), cuy::ii, HcalDetId::iphi(), edm::HandleBase::isValid(), maxDepth_, mc_, mode_, nevent1, nevent2, nevent3, nevent4, noise_, muonCSCDigis_cfi::pedestal, HcalCalibrations::pedestal(), CaloSamples::presamples(), edm::Handle< T >::product(), HcalHitRelabeller::relabel(), CaloSamples::size(), str(), AlCaHLTBitMon_QueryRunRegistry::string, HcalDetId::subdet(), subdet_, testNumber_, tok_mc_, and heppy_batch::val.

571  {
572 
573 
574  // HistLim =============================================================
575 
576  std::string strtmp;
577 
578  // ======================================================================
579  using namespace edm;
580  typename edm::Handle<edm::SortedCollection<Digi> > digiCollection;
581  typename edm::SortedCollection<Digi>::const_iterator digiItr;
582 
583  // ADC2fC
584  HcalCalibrations calibrations;
585  CaloSamples tool;
586  iEvent.getByToken(tok, digiCollection);
587  if (!digiCollection.isValid()) return;
588  int isubdet = 0;
589  if (subdet_ == "HB") isubdet = 1;
590  if (subdet_ == "HE") isubdet = 2;
591  if (subdet_ == "HO") isubdet = 3;
592  if (subdet_ == "HF") isubdet = 4;
593 
594  if (isubdet == 1) nevent1++;
595  if (isubdet == 2) nevent2++;
596  if (isubdet == 3) nevent3++;
597  if (isubdet == 4) nevent4++;
598 
599  int indigis = 0;
600  // amplitude for signal cell at diff. depths
601  std::vector<double> v_ampl_c(maxDepth_[isubdet]+1,0);
602 
603  // is set to 1 if "seed" SimHit is found
604  int seedSimHit = 0;
605 
606  int ieta_Sim = 9999;
607  int iphi_Sim = 9999;
608  double emax_Sim = -9999.;
609 
610 
611  // SimHits MC only
612  if (mc_ == "yes") {
613  edm::Handle<edm::PCaloHitContainer> hcalHits;
614  iEvent.getByToken(tok_mc_, hcalHits);
615  const edm::PCaloHitContainer * simhitResult = hcalHits.product();
616 
617  if (isubdet != 0 && noise_ == 0) { // signal only SimHits
618 
619  for (std::vector<PCaloHit>::const_iterator simhits = simhitResult->begin(); simhits != simhitResult->end(); ++simhits) {
620 
621  unsigned int id_ = simhits->id();
622  int sub, ieta, iphi;
623  HcalDetId hid;
624  if (testNumber_) hid = HcalHitRelabeller::relabel(id_,hcons);
625  else hid = HcalDetId(id_);
626  sub = hid.subdet();
627  ieta = hid.ieta();
628  iphi = hid.iphi();
629 
630  double en = simhits->energy();
631 
632  if (en > emax_Sim && sub == isubdet) {
633  emax_Sim = en;
634  ieta_Sim = ieta;
635  iphi_Sim = iphi;
636  // to limit "seed" SimHit energy in case of "multi" event
637  if (mode_ == "multi" &&
638  ((sub == 4 && en < 100. && en > 1.)
639  || ((sub != 4) && en < 1. && en > 0.02))) {
640  seedSimHit = 1;
641  break;
642  }
643  }
644 
645  } // end of SimHits cycle
646 
647 
648  // found highest-energy SimHit for single-particle
649  if (mode_ != "multi" && emax_Sim > 0.) seedSimHit = 1;
650  } // end of SimHits
651  }// end of mc_ == "yes"
652 
653  // CYCLE OVER CELLS ========================================================
654  int Ndig = 0;
655 
656  for (digiItr = digiCollection->begin(); digiItr != digiCollection->end(); digiItr++) {
657 
658  HcalDetId cell(digiItr->id());
659  int depth = cell.depth();
660  int iphi = cell.iphi();
661  int ieta = cell.ieta();
662  int sub = cell.subdet();
663 
664  if(depth > maxDepth_[isubdet] && sub == isubdet){
665  edm::LogWarning("HcalDetId") << "HcalDetID presents conflicting information. Depth: " << depth << ", iphi: " << iphi << ", ieta: " << ieta << ". Max depth from geometry is: " << maxDepth_[isubdet] << ". TestNumber = " << testNumber_;
666  continue;
667  }
668 
669  // amplitude for signal cell at diff. depths
670  std::vector<double> v_ampl(maxDepth_[isubdet]+1,0);
671 
672  // Gains, pedestals (once !) and only for "noise" case
673  if (((nevent1 == 1 && isubdet == 1) ||
674  (nevent2 == 1 && isubdet == 2) ||
675  (nevent3 == 1 && isubdet == 3) ||
676  (nevent4 == 1 && isubdet == 4)) && noise_ == 1 && sub == isubdet) {
677 
678  HcalGenericDetId hcalGenDetId(digiItr->id());
679  const HcalPedestal* pedestal = conditions->getPedestal(hcalGenDetId);
680  const HcalGain* gain = conditions->getGain(hcalGenDetId);
681  const HcalGainWidth* gainWidth = conditions->getGainWidth(hcalGenDetId);
682  const HcalPedestalWidth* pedWidth = conditions-> getPedestalWidth(hcalGenDetId);
683 
684  for (int i = 0; i < 4; i++) {
685  fill1D("HcalDigiTask_gain_capId" + str(i) + "_Depth" + str(depth) + "_" + subdet_, gain->getValue(i));
686  fill1D("HcalDigiTask_gainWidth_capId" + str(i) + "_Depth" + str(depth) + "_" + subdet_, gainWidth->getValue(i));
687  fill1D("HcalDigiTask_pedestal_capId" + str(i) + "_Depth" + str(depth) + "_" + subdet_, pedestal->getValue(i));
688  fill1D("HcalDigiTask_pedestal_width_capId" + str(i) + "_Depth" + str(depth) + "_" + subdet_, pedWidth->getWidth(i));
689  }
690 
691  fill2D("HcalDigiTask_gainMap_Depth" + str(depth) + "_" + subdet_, double(ieta), double(iphi), gain->getValue(0));
692  fill2D("HcalDigiTask_pwidthMap_Depth" + str(depth) + "_" + subdet_, double(ieta), double(iphi), pedWidth->getWidth(0));
693 
694  }// end of event #1
695 
696  if (sub == isubdet) Ndig++; // subdet number of digi
697 
698  // No-noise case, only single subdet selected ===========================
699 
700  if (sub == isubdet && noise_ == 0) {
701 
702  HcalCalibrations calibrations = conditions->getHcalCalibrations(cell);
703 
704  const HcalQIECoder* channelCoder = conditions->getHcalCoder(cell);
705  const HcalQIEShape* shape = conditions->getHcalShape(channelCoder);
706  HcalCoderDb coder(*channelCoder, *shape);
707  coder.adc2fC(*digiItr, tool);
708 
709  double noiseADC = (*digiItr)[0].adc();
710  double noisefC = tool[0];
711  // noise evaluations from "pre-samples"
712  fill1D("HcalDigiTask_ADC0_adc_depth" + str(depth) + "_" + subdet_, noiseADC);
713  fill1D("HcalDigiTask_ADC0_fC_depth" + str(depth) + "_" + subdet_, noisefC);
714 
715 
716  // OCCUPANCY maps fill
717  fill2D("HcalDigiTask_ieta_iphi_occupancy_map_depth" + str(depth) + "_" + subdet_, double(ieta), double(iphi));
718 
719  fill1D("HcalDigiTask_depths_"+subdet_,double(depth));
720 
721  // Cycle on time slices
722  // - for each Digi
723  // - for one Digi with max SimHits E in subdet
724 
725 
726  int closen = 0; // =1 if 1) seedSimHit = 1 and 2) the cell is the same
727  if (ieta == ieta_Sim && iphi == iphi_Sim) closen = seedSimHit;
728 
729  for (int ii = 0; ii < tool.size(); ii++) {
730  int capid = (*digiItr)[ii].capid();
731  // single ts amplitude
732  double val = (tool[ii] - calibrations.pedestal(capid));
733 
734  if (val > 100.) {
735  fill1D("HcalDigiTask_ADC0_adc_depth" + str(depth) + "_" + subdet_, noiseADC);
736  strtmp = "HcalDigiTask_all_amplitudes_vs_bin_1D_depth" + str(depth) + "_" + subdet_;
737  fill1D(strtmp, double(ii), val);
738  }
739 
740  if (closen == 1) {
741  strtmp = "HcalDigiTask_signal_amplitude_vs_bin_all_depths_" + subdet_;
742  fill2D(strtmp, double(ii), val);
743  }
744 
745 
746  // HB/HE/HO
747  if (isubdet != 4 && ii >= 4 && ii <= 7) {
748  v_ampl[0] += val;
749  v_ampl[depth] += val;
750 
751  if (closen == 1) {
752  v_ampl_c[0] += val;
753  v_ampl_c[depth] += val;
754  }
755  }
756 
757  // HF
758  if (isubdet == 4 && ii >= 2 && ii <= 4) {
759  v_ampl[0] += val;
760  v_ampl[depth] += val;
761  if (closen == 1) {
762  v_ampl_c[0] += val;
763  v_ampl_c[depth] += val;
764  }
765  }
766  }
767  // end of time bucket sample
768 
769  // maps of sum of amplitudes (sum lin.digis(4,5,6,7) - ped) all depths
770  // just 1D of all cells' amplitudes
771  strtmp = "HcalDigiTask_sum_all_amplitudes_" + subdet_;
772  fill1D(strtmp, v_ampl[0]);
773 
774  std::vector<int> v_ampl_sub(v_ampl.begin() + 1, v_ampl.end()); // remove element 0, which is the sum of any depth
775  double ampl_max = *std::max_element(v_ampl_sub.begin(),v_ampl_sub.end());
776  if (ampl_max>10.) indigis++;
777  //KH if (ampl1 > 10. || ampl2 > 10. || ampl3 > 10. || ampl4 > 10.) indigis++;
778 
779  // fraction 5,6 bins if ampl. is big.
780  if (v_ampl[1] > 30. && depth == 1) {
781  int soi = tool.presamples();
782  int lastbin = tool.size() - 1;
783 
784  double fbinSOI = tool[soi] - calibrations.pedestal((*digiItr)[soi].capid());
785  double fbinPS = 0;
786 
787  for(int j = soi+1; j <= lastbin; j++) fbinPS += tool[j] - calibrations.pedestal((*digiItr)[j].capid());
788 
789  fbinSOI /= v_ampl[1];
790  fbinPS /= v_ampl[1];
791  strtmp = "HcalDigiTask_SOI_frac_" + subdet_;
792  fill1D(strtmp, fbinSOI);
793  strtmp = "HcalDigiTask_postSOI_frac_" + subdet_;
794  fill1D(strtmp, fbinPS);
795  }
796 
797  strtmp = "HcalDigiTask_signal_amplitude_" + subdet_;
798  fill1D(strtmp, v_ampl[0]);
799  strtmp = "HcalDigiTask_signal_amplitude_depth" + str(depth) + "_" + subdet_;
800  fill1D(strtmp, v_ampl[depth]);
801  }
802  } // End of CYCLE OVER CELLS =============================================
803 
804  if (isubdet != 0 && noise_ == 0) { // signal only, once per event
805  strtmp = "HcalDigiTask_number_of_amplitudes_above_10fC_" + subdet_;
806  fill1D(strtmp, indigis);
807 
808  // SimHits once again !!!
809  double eps = 1.e-3;
810  std::vector<double> v_ehits(maxDepth_[isubdet]+1,0);
811 
812  if (mc_ == "yes") {
813  edm::Handle<edm::PCaloHitContainer> hcalHits;
814  iEvent.getByToken(tok_mc_, hcalHits);
815  const edm::PCaloHitContainer * simhitResult = hcalHits.product();
816  for (std::vector<PCaloHit>::const_iterator simhits = simhitResult->begin(); simhits != simhitResult->end(); ++simhits) {
817 
818  unsigned int id_ = simhits->id();
819  int sub, depth, ieta, iphi;
820  HcalDetId hid;
821  if (testNumber_) hid = HcalHitRelabeller::relabel(id_,hcons);
822  else hid = HcalDetId(id_);
823  sub = hid.subdet();
824  depth = hid.depth();
825  ieta = hid.ieta();
826  iphi = hid.iphi();
827 
828  if(depth > maxDepth_[isubdet] && sub == isubdet){
829  edm::LogWarning("HcalDetId") << "HcalDetID(SimHit) presents conflicting information. Depth: " << depth << ", iphi: " << iphi << ", ieta: " << ieta << ". Max depth from geometry is: " << maxDepth_[isubdet] << ". TestNumber = " << testNumber_;
830  continue;
831  }
832 
833 
834  // take cell already found to be max energy in a particular subdet
835  if (sub == isubdet && ieta == ieta_Sim && iphi == iphi_Sim) {
836  double en = simhits->energy();
837 
838  v_ehits[0] += en;
839  v_ehits[depth] += en;
840 
841  }
842  } // simhit loop
843 
844  strtmp = "HcalDigiTask_amplitude_vs_simhits_" + subdet_;
845  if (v_ehits[0] > eps) fill2D(strtmp, v_ehits[0], v_ampl_c[0]);
846  for (int depth = 1; depth <= maxDepth_[isubdet]; depth++) {
847  strtmp = "HcalDigiTask_amplitude_vs_simhits_depth" + str(depth) + "_" + subdet_;
848  if (v_ehits[depth] > eps) fill2D(strtmp, v_ehits[depth], v_ampl_c[depth]);
849  }
850 
851  strtmp = "HcalDigiTask_amplitude_vs_simhits_profile_" + subdet_;
852  if (v_ehits[0] > eps) fillPf(strtmp, v_ehits[0], v_ampl_c[0]);
853  for (int depth = 1; depth <= maxDepth_[isubdet]; depth++) {
854  strtmp = "HcalDigiTask_amplitude_vs_simhits_profile_depth" + str(depth) + "_" + subdet_;
855  if (v_ehits[depth] > eps) fillPf(strtmp, v_ehits[depth], v_ampl_c[depth]);
856  }
857 
858  strtmp = "HcalDigiTask_ratio_amplitude_vs_simhits_" + subdet_;
859  if (v_ehits[0] > eps) fill1D(strtmp, v_ampl_c[0] / v_ehits[0]);
860  for (int depth = 1; depth <= maxDepth_[isubdet]; depth++) {
861  strtmp = "HcalDigiTask_amplitude_vs_simhits_profile_depth" + str(depth) + "_" + subdet_;
862  if (v_ehits[depth] > eps) fillPf(strtmp, v_ehits[depth], v_ampl_c[depth]);
863  strtmp = "HcalDigiTask_ratio_amplitude_vs_simhits_depth" + str(depth) + "_" + subdet_;
864  if (v_ehits[depth] > eps) fill1D(strtmp, v_ampl_c[depth] / v_ehits[depth]);
865  }
866 
867 
868  } // end of if(mc_ == "yes")
869 
870  strtmp = "HcalDigiTask_Ndigis_" + subdet_;
871  fill1D(strtmp, double(Ndig));
872 
873  } // end of if( subdet != 0 && noise_ == 0) { // signal only
874 
875 }
HcalDbService::getGainWidth
const HcalGainWidth * getGainWidth(const HcalGenericDetId &fId) const
Definition: HcalDbService.cc:237
edm::PCaloHitContainer
std::vector< PCaloHit > PCaloHitContainer
Definition: PCaloHitContainer.h:8
mps_fire.i
i
Definition: mps_fire.py:156
HcalDetId::subdet
HcalSubdetector subdet() const
get the subdetector
Definition: HcalDetId.h:49
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:460
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
edm::SortedCollection::const_iterator
std::vector< T >::const_iterator const_iterator
Definition: SortedCollection.h:82
CaloSamples::presamples
int presamples() const
access presample information
Definition: CaloSamples.h:36
HcalCalibrations::pedestal
double pedestal(int fCapId) const
get pedestal for capid=0..3
Definition: HcalCalibrations.h:20
HcalGain::getValue
float getValue(int fCapId) const
get value for capId = 0..3
Definition: HcalGain.h:22
HcalGainWidth::getValue
float getValue(int fCapId) const
get value for capId = 0..3
Definition: HcalGainWidth.h:21
hgcalDigitizer_cfi.digiCollection
digiCollection
Definition: hgcalDigitizer_cfi.py:12
HcalDetId::depth
int depth() const
get the tower depth
Definition: HcalDetId.cc:108
particleFlowClusterECALTimeSelected_cfi.depth
depth
Definition: particleFlowClusterECALTimeSelected_cfi.py:8
HcalDigisValidation::fillPf
void fillPf(std::string name, double X, double Y)
Definition: HcalDigisValidation.cc:1255
HcalDigisValidation::fill2D
void fill2D(std::string name, double X, double Y, double weight=1)
Definition: HcalDigisValidation.cc:1243
HcalDetId::ieta
int ieta() const
get the cell ieta
Definition: HcalDetId.h:56
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:74
HcalDigisValidation::hcons
const HcalDDDRecConstants * hcons
Definition: HcalDigisValidation.h:149
cuy.ii
ii
Definition: cuy.py:588
muonCSCDigis_cfi.gain
gain
Definition: muonCSCDigis_cfi.py:36
HcalDetId::iphi
int iphi() const
get the cell iphi
Definition: HcalDetId.cc:103
AlignmentProducer_cff.calibrations
calibrations
Definition: AlignmentProducer_cff.py:59
edm::Handle::product
T const * product() const
Definition: Handle.h:81
HcalDigisValidation::fill1D
void fill1D(std::string name, double X, double weight=1)
Definition: HcalDigisValidation.cc:1235
CaloSamples::size
int size() const
get the size
Definition: CaloSamples.h:24
HcalDigisValidation::tok_mc_
edm::EDGetTokenT< edm::PCaloHitContainer > tok_mc_
Definition: HcalDigisValidation.h:125
HcalDbService::getGain
const HcalGain * getGain(const HcalGenericDetId &fId) const
Definition: HcalDbService.cc:230
HcalPedestalWidth::getWidth
float getWidth(int fCapId) const
get width (sqrt(sigma_i_i)) for capId = 0..3
Definition: HcalPedestalWidth.cc:30
HcalDbService::getHcalCoder
const HcalQIECoder * getHcalCoder(const HcalGenericDetId &fId) const
Definition: HcalDbService.cc:244
HcalDbService::getHcalShape
const HcalQIEShape * getHcalShape(const HcalGenericDetId &fId) const
Definition: HcalDbService.cc:251
HcalDigisValidation::str
std::string str(int x)
Definition: HcalDigisValidation.cc:1264
muonCSCDigis_cfi.pedestal
pedestal
Definition: muonCSCDigis_cfi.py:12
edm
HLT enums.
Definition: AlignableModifier.h:17
HcalDigisValidation::conditions
edm::ESHandle< HcalDbService > conditions
Definition: HcalDigisValidation.h:137
HcalHitRelabeller::relabel
DetId relabel(const uint32_t testId) const
Definition: HcalHitRelabeller.cc:48
HcalDbService::getHcalCalibrations
const HcalCalibrations & getHcalCalibrations(const HcalGenericDetId &fId) const
Definition: HcalDbService.cc:50
HcalDbService::getPedestal
const HcalPedestal * getPedestal(const HcalGenericDetId &fId) const
Definition: HcalDbService.cc:216
heppy_batch.val
val
Definition: heppy_batch.py:349
template<class dataFrameType >
void HcalDigisValidation::reco ( const edm::Event &  iEvent,
const edm::EventSetup &  iSetup,
const edm::EDGetTokenT< HcalDataFrameContainer< dataFrameType > > &  tok 
)
private

Definition at line 876 of file HcalDigisValidation.cc.

References ecalMGPA::adc(), HcalCoderDb::adc2fC(), edm::DataFrameContainer::begin(), AlignmentProducer_cff::calibrations, conditions, particleFlowClusterECALTimeSelected_cfi::depth, HcalDetId::depth(), hgcalDigitizer_cfi::digiCollection, edm::DataFrameContainer::end(), fill1D(), fill2D(), fillPf(), muonCSCDigis_cfi::gain, edm::Event::getByToken(), HcalDbService::getGain(), HcalDbService::getGainWidth(), HcalDbService::getHcalCalibrations(), HcalDbService::getHcalCoder(), HcalDbService::getHcalShape(), HcalDbService::getPedestal(), HcalGainWidth::getValue(), HcalGain::getValue(), HcalPedestalWidth::getWidth(), hcons, hep17_, mps_fire::i, HcalDetId::ieta(), cuy::ii, HcalDetId::iphi(), edm::HandleBase::isValid(), maxDepth_, mc_, mode_, nevent1, nevent2, nevent3, nevent4, noise_, muonCSCDigis_cfi::pedestal, HcalCalibrations::pedestal(), CaloSamples::presamples(), edm::Handle< T >::product(), HcalHitRelabeller::relabel(), CaloSamples::size(), str(), AlCaHLTBitMon_QueryRunRegistry::string, HcalDetId::subdet(), subdet_, testNumber_, tok_mc_, and heppy_batch::val.

876  {
877 
878 
879  // HistLim =============================================================
880 
881  std::string strtmp;
882 
883  // ======================================================================
884  using namespace edm;
885  typename edm::Handle<HcalDataFrameContainer<dataFrameType> > digiHandle;
886  //typename HcalDataFrameContainer<dataFrameType>::const_iterator digiItr;
887 
888 
889  // ADC2fC
890  HcalCalibrations calibrations;
891  CaloSamples tool;
892  iEvent.getByToken(tok, digiHandle);
893  if (!digiHandle.isValid()) return;
894  const HcalDataFrameContainer<dataFrameType> *digiCollection = digiHandle.product();
895  int isubdet = 0;
896  if (subdet_ == "HB") isubdet = 1;
897  if (subdet_ == "HE") isubdet = 2;
898  if (subdet_ == "HO") isubdet = 3;
899  if (subdet_ == "HF") isubdet = 4;
900 
901  if (isubdet == 1) nevent1++;
902  if (isubdet == 2) nevent2++;
903  if (isubdet == 3) nevent3++;
904  if (isubdet == 4) nevent4++;
905 
906  int indigis = 0;
907  // amplitude for signal cell at diff. depths
908  std::vector<double> v_ampl_c(maxDepth_[isubdet]+1,0);
909 
910  // is set to 1 if "seed" SimHit is found
911  int seedSimHit = 0;
912 
913  int ieta_Sim = 9999;
914  int iphi_Sim = 9999;
915  double emax_Sim = -9999.;
916 
917 
918  // SimHits MC only
919  if (mc_ == "yes") {
920  edm::Handle<edm::PCaloHitContainer> hcalHits;
921  iEvent.getByToken(tok_mc_, hcalHits);
922  const edm::PCaloHitContainer * simhitResult = hcalHits.product();
923 
924  if (isubdet != 0 && noise_ == 0) { // signal only SimHits
925 
926  for (std::vector<PCaloHit>::const_iterator simhits = simhitResult->begin(); simhits != simhitResult->end(); ++simhits) {
927 
928  unsigned int id_ = simhits->id();
929  int sub, ieta, iphi;
930  HcalDetId hid;
931  if (testNumber_) hid = HcalHitRelabeller::relabel(id_,hcons);
932  else hid = HcalDetId(id_);
933  sub = hid.subdet();
934  ieta = hid.ieta();
935  iphi = hid.iphi();
936 
937  double en = simhits->energy();
938 
939  if (en > emax_Sim && sub == isubdet) {
940  emax_Sim = en;
941  ieta_Sim = ieta;
942  iphi_Sim = iphi;
943  // to limit "seed" SimHit energy in case of "multi" event
944  if (mode_ == "multi" &&
945  ((sub == 4 && en < 100. && en > 1.)
946  || ((sub != 4) && en < 1. && en > 0.02))) {
947  seedSimHit = 1;
948  break;
949  }
950  }
951 
952  } // end of SimHits cycle
953 
954 
955  // found highest-energy SimHit for single-particle
956  if (mode_ != "multi" && emax_Sim > 0.) seedSimHit = 1;
957  } // end of SimHits
958  }// end of mc_ == "yes"
959 
960  // CYCLE OVER CELLS ========================================================
961  int Ndig = 0;
962 
963  for (typename HcalDataFrameContainer<dataFrameType>::const_iterator digiItr = digiCollection->begin(); digiItr != digiCollection->end(); digiItr++) {
964 
965  dataFrameType dataFrame = *digiItr;
966 
967  HcalDetId cell(digiItr->id());
968  int depth = cell.depth();
969  int iphi = cell.iphi();
970  int ieta = cell.ieta();
971  int sub = cell.subdet();
972 
973  //Is this in HEP17
974  bool isHEP17 = (iphi>=63)&&(iphi<=66)&&(ieta>0)&&(sub == 2);
975 
976 
977  if(depth > maxDepth_[isubdet] && sub == isubdet){
978  edm::LogWarning("HcalDetId") << "HcalDetID presents conflicting information. Depth: " << depth << ", iphi: " << iphi << ", ieta: " << ieta << ". Max depth from geometry is: " << maxDepth_[isubdet] << ". TestNumber = " << testNumber_;
979  continue;
980  }
981 
982  // amplitude for signal cell at diff. depths
983  std::vector<double> v_ampl(maxDepth_[isubdet]+1,0);
984 
985  // Gains, pedestals (once !) and only for "noise" case
986  if (((nevent1 == 1 && isubdet == 1) ||
987  (nevent2 == 1 && isubdet == 2) ||
988  (nevent3 == 1 && isubdet == 3) ||
989  (nevent4 == 1 && isubdet == 4)) && noise_ == 1 && sub == isubdet) {
990 
991  HcalGenericDetId hcalGenDetId(digiItr->id());
992  const HcalPedestal* pedestal = conditions->getPedestal(hcalGenDetId);
993  const HcalGain* gain = conditions->getGain(hcalGenDetId);
994  const HcalGainWidth* gainWidth = conditions->getGainWidth(hcalGenDetId);
995  const HcalPedestalWidth* pedWidth = conditions-> getPedestalWidth(hcalGenDetId);
996 
997  for (int i = 0; i < 4; i++) {
998  fill1D("HcalDigiTask_gain_capId" + str(i) + "_Depth" + str(depth) + "_" + subdet_, gain->getValue(i));
999  fill1D("HcalDigiTask_gainWidth_capId" + str(i) + "_Depth" + str(depth) + "_" + subdet_, gainWidth->getValue(i));
1000  fill1D("HcalDigiTask_pedestal_capId" + str(i) + "_Depth" + str(depth) + "_" + subdet_, pedestal->getValue(i));
1001  fill1D("HcalDigiTask_pedestal_width_capId" + str(i) + "_Depth" + str(depth) + "_" + subdet_, pedWidth->getWidth(i));
1002  }
1003 
1004  fill2D("HcalDigiTask_gainMap_Depth" + str(depth) + "_" + subdet_, double(ieta), double(iphi), gain->getValue(0));
1005  fill2D("HcalDigiTask_pwidthMap_Depth" + str(depth) + "_" + subdet_, double(ieta), double(iphi), pedWidth->getWidth(0));
1006 
1007  }// end of event #1
1008  //std::cout << "==== End of event noise block in cell cycle" << std::endl;
1009 
1010  if (sub == isubdet) Ndig++; // subdet number of digi
1011 
1012  // No-noise case, only single subdet selected ===========================
1013 
1014  if (sub == isubdet && noise_ == 0) {
1015 
1016  HcalCalibrations calibrations = conditions->getHcalCalibrations(cell);
1017 
1018  const HcalQIECoder* channelCoder = conditions->getHcalCoder(cell);
1019  const HcalQIEShape* shape = conditions->getHcalShape(channelCoder);
1020  HcalCoderDb coder(*channelCoder, *shape);
1021  coder.adc2fC(dataFrame, tool);
1022 
1023  double noiseADC = (dataFrame)[0].adc();
1024  double noisefC = tool[0];
1025  // noise evaluations from "pre-samples"
1026  fill1D("HcalDigiTask_ADC0_adc_depth" + str(depth) + "_" + subdet_, noiseADC);
1027  fill1D("HcalDigiTask_ADC0_fC_depth" + str(depth) + "_" + subdet_, noisefC);
1028 
1029 
1030  // OCCUPANCY maps fill
1031  fill2D("HcalDigiTask_ieta_iphi_occupancy_map_depth" + str(depth) + "_" + subdet_, double(ieta), double(iphi));
1032 
1033  fill1D("HcalDigiTask_depths_"+subdet_,double(depth));
1034 
1035  // Cycle on time slices
1036  // - for each Digi
1037  // - for one Digi with max SimHits E in subdet
1038 
1039 
1040  int closen = 0; // =1 if 1) seedSimHit = 1 and 2) the cell is the same
1041  if (ieta == ieta_Sim && iphi == iphi_Sim) closen = seedSimHit;
1042 
1043  for (int ii = 0; ii < tool.size(); ii++) {
1044  int capid = (dataFrame)[ii].capid();
1045  // single ts amplitude
1046  double val = (tool[ii] - calibrations.pedestal(capid));
1047 
1048  if (val > 100.) {
1049  fill1D("HcalDigiTask_ADC0_adc_depth" + str(depth) + "_" + subdet_, noiseADC);
1050  if(hep17_){
1051  if(!isHEP17){
1052  strtmp = "HcalDigiTask_all_amplitudes_vs_bin_1D_depth" + str(depth) + "_" + subdet_;
1053  fill1D(strtmp, double(ii), val);
1054  } else {
1055  strtmp = "HcalDigiTask_all_amplitudes_vs_bin_1D_depth" + str(depth) + "_HEP17";
1056  fill1D(strtmp, double(ii), val);
1057  }
1058  } else {
1059  strtmp = "HcalDigiTask_all_amplitudes_vs_bin_1D_depth" + str(depth) + "_" + subdet_;
1060  fill1D(strtmp, double(ii), val);
1061  }
1062  }
1063 
1064  if (closen == 1) {
1065  if(hep17_){
1066  if(!isHEP17){
1067  strtmp = "HcalDigiTask_signal_amplitude_vs_bin_all_depths_" + subdet_;
1068  fill2D(strtmp, double(ii), val);
1069  }else{
1070  strtmp = "HcalDigiTask_signal_amplitude_vs_bin_all_depths_HEP17";
1071  fill2D(strtmp, double(ii), val);
1072  }
1073  } else {
1074  strtmp = "HcalDigiTask_signal_amplitude_vs_bin_all_depths_" + subdet_;
1075  fill2D(strtmp, double(ii), val);
1076  }
1077  }
1078 
1079 
1080  // HB/HE/HO
1081  if (isubdet != 4 && ii >= 4 && ii <= 7) {
1082  v_ampl[0] += val;
1083  v_ampl[depth] += val;
1084 
1085  if (closen == 1) {
1086  v_ampl_c[0] += val;
1087  v_ampl_c[depth] += val;
1088  }
1089  }
1090 
1091  // HF
1092  if (isubdet == 4 && ii >= 2 && ii <= 4) {
1093  v_ampl[0] += val;
1094  v_ampl[depth] += val;
1095  if (closen == 1) {
1096  v_ampl_c[0] += val;
1097  v_ampl_c[depth] += val;
1098  }
1099  }
1100  }
1101  // end of time bucket sample
1102 
1103  // just 1D of all cells' amplitudes
1104  strtmp = "HcalDigiTask_sum_all_amplitudes_" + subdet_;
1105  fill1D(strtmp, v_ampl[0]);
1106 
1107  std::vector<int> v_ampl_sub(v_ampl.begin() + 1, v_ampl.end()); // remove element 0, which is the sum of any depth
1108  double ampl_max = *std::max_element(v_ampl_sub.begin(),v_ampl_sub.end());
1109  if (ampl_max>10.) indigis++;
1110  //KH if (ampl1 > 10. || ampl2 > 10. || ampl3 > 10. || ampl4 > 10.) indigis++;
1111 
1112  // fraction 5,6 bins if ampl. is big.
1113  //histogram names have not been changed, but it should be understood that bin_5 is soi, and bin_6_7 is latter TS'
1114  if (v_ampl[1] > 30. && depth == 1) {
1115  int soi = tool.presamples();
1116  int lastbin = tool.size() - 1;
1117 
1118  double fbinSOI = tool[soi] - calibrations.pedestal((dataFrame)[soi].capid());
1119  double fbinPS = 0;
1120 
1121  for(int j = soi+1; j <= lastbin; j++) fbinPS += tool[j] - calibrations.pedestal((dataFrame)[j].capid());
1122 
1123 
1124  fbinSOI /= v_ampl[1];
1125  fbinPS /= v_ampl[1];
1126  strtmp = "HcalDigiTask_SOI_frac_" + subdet_;
1127  fill1D(strtmp, fbinSOI);
1128  strtmp = "HcalDigiTask_postSOI_frac_" + subdet_;
1129  fill1D(strtmp, fbinPS);
1130  }
1131 
1132 
1133  if(hep17_){
1134  if(!isHEP17){
1135  strtmp = "HcalDigiTask_signal_amplitude_" + subdet_;
1136  fill1D(strtmp, v_ampl[0]);
1137  strtmp = "HcalDigiTask_signal_amplitude_depth" + str(depth) + "_" + subdet_;
1138  fill1D(strtmp, v_ampl[depth]);
1139  } else {
1140  strtmp = "HcalDigiTask_signal_amplitude_HEP17";
1141  fill1D(strtmp, v_ampl[0]);
1142  strtmp = "HcalDigiTask_signal_amplitude_depth" + str(depth) + "_HEP17";
1143  fill1D(strtmp, v_ampl[depth]);
1144  }
1145  }else{
1146  strtmp = "HcalDigiTask_signal_amplitude_" + subdet_;
1147  fill1D(strtmp, v_ampl[0]);
1148  strtmp = "HcalDigiTask_signal_amplitude_depth" + str(depth) + "_" + subdet_;
1149  fill1D(strtmp, v_ampl[depth]);
1150  }
1151  }
1152  } // End of CYCLE OVER CELLS =============================================
1153 
1154  if (isubdet != 0 && noise_ == 0) { // signal only, once per event
1155  strtmp = "HcalDigiTask_number_of_amplitudes_above_10fC_" + subdet_;
1156  fill1D(strtmp, indigis);
1157 
1158  // SimHits once again !!!
1159  double eps = 1.e-3;
1160  std::vector<double> v_ehits(maxDepth_[isubdet]+1,0);
1161 
1162  if (mc_ == "yes") {
1163  edm::Handle<edm::PCaloHitContainer> hcalHits;
1164  iEvent.getByToken(tok_mc_, hcalHits);
1165  const edm::PCaloHitContainer * simhitResult = hcalHits.product();
1166  for (std::vector<PCaloHit>::const_iterator simhits = simhitResult->begin(); simhits != simhitResult->end(); ++simhits) {
1167 
1168  unsigned int id_ = simhits->id();
1169  int sub, depth, ieta, iphi;
1170  HcalDetId hid;
1171  if (testNumber_) hid = HcalHitRelabeller::relabel(id_,hcons);
1172  else hid = HcalDetId(id_);
1173  sub = hid.subdet();
1174  depth = hid.depth();
1175  ieta = hid.ieta();
1176  iphi = hid.iphi();
1177 
1178  if(depth > maxDepth_[isubdet] && sub == isubdet){
1179  edm::LogWarning("HcalDetId") << "HcalDetID(SimHit) presents conflicting information. Depth: " << depth << ", iphi: " << iphi << ", ieta: " << ieta << ". Max depth from geometry is: " << maxDepth_[isubdet] << ". TestNumber = " << testNumber_;
1180  continue;
1181  }
1182 
1183 
1184  // take cell already found to be max energy in a particular subdet
1185  if (sub == isubdet && ieta == ieta_Sim && iphi == iphi_Sim) {
1186  double en = simhits->energy();
1187 
1188  v_ehits[0] += en;
1189  v_ehits[depth] += en;
1190 
1191  }
1192  } // simhit loop
1193 
1194  strtmp = "HcalDigiTask_amplitude_vs_simhits_" + subdet_;
1195  if (v_ehits[0] > eps) fill2D(strtmp, v_ehits[0], v_ampl_c[0]);
1196  for (int depth = 1; depth <= maxDepth_[isubdet]; depth++) {
1197  strtmp = "HcalDigiTask_amplitude_vs_simhits_depth" + str(depth) + "_" + subdet_;
1198  if (v_ehits[depth] > eps) fill2D(strtmp, v_ehits[depth], v_ampl_c[depth]);
1199  }
1200 
1201  strtmp = "HcalDigiTask_amplitude_vs_simhits_profile_" + subdet_;
1202  if (v_ehits[0] > eps) fillPf(strtmp, v_ehits[0], v_ampl_c[0]);
1203  for (int depth = 1; depth <= maxDepth_[isubdet]; depth++) {
1204  strtmp = "HcalDigiTask_amplitude_vs_simhits_profile_depth" + str(depth) + "_" + subdet_;
1205  if (v_ehits[depth] > eps) fillPf(strtmp, v_ehits[depth], v_ampl_c[depth]);
1206  }
1207 
1208  strtmp = "HcalDigiTask_ratio_amplitude_vs_simhits_" + subdet_;
1209  if (v_ehits[0] > eps) fill1D(strtmp, v_ampl_c[0] / v_ehits[0]);
1210  for (int depth = 1; depth <= maxDepth_[isubdet]; depth++) {
1211  strtmp = "HcalDigiTask_amplitude_vs_simhits_profile_depth" + str(depth) + "_" + subdet_;
1212  if (v_ehits[depth] > eps) fillPf(strtmp, v_ehits[depth], v_ampl_c[depth]);
1213  strtmp = "HcalDigiTask_ratio_amplitude_vs_simhits_depth" + str(depth) + "_" + subdet_;
1214  if (v_ehits[depth] > eps) fill1D(strtmp, v_ampl_c[depth] / v_ehits[depth]);
1215  }
1216 
1217 
1218  } // end of if(mc_ == "yes")
1219 
1220  strtmp = "HcalDigiTask_Ndigis_" + subdet_;
1221  fill1D(strtmp, double(Ndig));
1222 
1223  } // end of if( subdet != 0 && noise_ == 0) { // signal only
1224 
1225 } //HcalDataFrameContainer
ecalMGPA::adc
int adc(sample_type sample)
get the ADC sample (12 bits)
Definition: EcalMGPASample.h:11
HcalDbService::getGainWidth
const HcalGainWidth * getGainWidth(const HcalGenericDetId &fId) const
Definition: HcalDbService.cc:237
edm::PCaloHitContainer
std::vector< PCaloHit > PCaloHitContainer
Definition: PCaloHitContainer.h:8
mps_fire.i
i
Definition: mps_fire.py:156
edm::DataFrameContainer::const_iterator
boost::transform_iterator< IterHelp, boost::counting_iterator< int > > const_iterator
Definition: DataFrameContainer.h:58
HcalDetId::subdet
HcalSubdetector subdet() const
get the subdetector
Definition: HcalDetId.h:49
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:460
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
CaloSamples::presamples
int presamples() const
access presample information
Definition: CaloSamples.h:36
HcalCalibrations::pedestal
double pedestal(int fCapId) const
get pedestal for capid=0..3
Definition: HcalCalibrations.h:20
edm::DataFrameContainer::begin
const_iterator begin() const
Definition: DataFrameContainer.h:159
HcalGain::getValue
float getValue(int fCapId) const
get value for capId = 0..3
Definition: HcalGain.h:22
HcalGainWidth::getValue
float getValue(int fCapId) const
get value for capId = 0..3
Definition: HcalGainWidth.h:21
hgcalDigitizer_cfi.digiCollection
digiCollection
Definition: hgcalDigitizer_cfi.py:12
HcalDetId::depth
int depth() const
get the tower depth
Definition: HcalDetId.cc:108
particleFlowClusterECALTimeSelected_cfi.depth
depth
Definition: particleFlowClusterECALTimeSelected_cfi.py:8
HcalDigisValidation::fillPf
void fillPf(std::string name, double X, double Y)
Definition: HcalDigisValidation.cc:1255
HcalDigisValidation::fill2D
void fill2D(std::string name, double X, double Y, double weight=1)
Definition: HcalDigisValidation.cc:1243
HcalDetId::ieta
int ieta() const
get the cell ieta
Definition: HcalDetId.h:56
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:74
HcalDigisValidation::hcons
const HcalDDDRecConstants * hcons
Definition: HcalDigisValidation.h:149
cuy.ii
ii
Definition: cuy.py:588
muonCSCDigis_cfi.gain
gain
Definition: muonCSCDigis_cfi.py:36
HcalDetId::iphi
int iphi() const
get the cell iphi
Definition: HcalDetId.cc:103
AlignmentProducer_cff.calibrations
calibrations
Definition: AlignmentProducer_cff.py:59
edm::Handle::product
T const * product() const
Definition: Handle.h:81
HcalDigisValidation::fill1D
void fill1D(std::string name, double X, double weight=1)
Definition: HcalDigisValidation.cc:1235
CaloSamples::size
int size() const
get the size
Definition: CaloSamples.h:24
HcalDigisValidation::tok_mc_
edm::EDGetTokenT< edm::PCaloHitContainer > tok_mc_
Definition: HcalDigisValidation.h:125
HcalDbService::getGain
const HcalGain * getGain(const HcalGenericDetId &fId) const
Definition: HcalDbService.cc:230
HcalPedestalWidth::getWidth
float getWidth(int fCapId) const
get width (sqrt(sigma_i_i)) for capId = 0..3
Definition: HcalPedestalWidth.cc:30
HcalDbService::getHcalCoder
const HcalQIECoder * getHcalCoder(const HcalGenericDetId &fId) const
Definition: HcalDbService.cc:244
HcalDbService::getHcalShape
const HcalQIEShape * getHcalShape(const HcalGenericDetId &fId) const
Definition: HcalDbService.cc:251
HcalDigisValidation::str
std::string str(int x)
Definition: HcalDigisValidation.cc:1264
muonCSCDigis_cfi.pedestal
pedestal
Definition: muonCSCDigis_cfi.py:12
edm::DataFrameContainer::end
const_iterator end() const
Definition: DataFrameContainer.h:163
edm
HLT enums.
Definition: AlignableModifier.h:17
HcalDigisValidation::conditions
edm::ESHandle< HcalDbService > conditions
Definition: HcalDigisValidation.h:137
HcalHitRelabeller::relabel
DetId relabel(const uint32_t testId) const
Definition: HcalHitRelabeller.cc:48
HcalDbService::getHcalCalibrations
const HcalCalibrations & getHcalCalibrations(const HcalGenericDetId &fId) const
Definition: HcalDbService.cc:50
HcalDbService::getPedestal
const HcalPedestal * getPedestal(const HcalGenericDetId &fId) const
Definition: HcalDbService.cc:216
heppy_batch.val
val
Definition: heppy_batch.py:349
std::string HcalDigisValidation::str ( int  x)
private

Definition at line 1264 of file HcalDigisValidation.cc.

References DEFINE_FWK_MODULE, MillePedeFileConverter_cfg::out, and x.

Referenced by reco().

1264  {
1265  std::stringstream out;
1266  out << x;
1267  return out.str();
1268 }
MillePedeFileConverter_cfg.out
out
Definition: MillePedeFileConverter_cfg.py:31

Member Data Documentation

edm::ESHandle<HcalDbService> HcalDigisValidation::conditions
private

Definition at line 137 of file HcalDigisValidation.h.

Referenced by analyze(), and reco().

edm::InputTag HcalDigisValidation::dataTPsTag_
private

Definition at line 118 of file HcalDigisValidation.h.

Referenced by HcalDigisValidation().

std::string HcalDigisValidation::dirName_
private

Definition at line 112 of file HcalDigisValidation.h.

Referenced by bookHistograms(), and HcalDigisValidation().

edm::InputTag HcalDigisValidation::emulTPsTag_
private

Definition at line 117 of file HcalDigisValidation.h.

Referenced by HcalDigisValidation().

edm::ESHandle<CaloGeometry> HcalDigisValidation::geometry
private

Definition at line 135 of file HcalDigisValidation.h.

Referenced by Vispa.Gui.ConnectableWidget.ConnectableWidget::leaveEvent().

const HcalDDDRecConstants* HcalDigisValidation::hcons
private

Definition at line 149 of file HcalDigisValidation.h.

Referenced by dqmBeginRun(), and reco().

bool HcalDigisValidation::hep17_
private

Definition at line 123 of file HcalDigisValidation.h.

Referenced by booking(), HcalDigisValidation(), and reco().

edm::ESHandle<HcalTopology> HcalDigisValidation::htopo
private

Definition at line 140 of file HcalDigisValidation.h.

Referenced by analyze().

const HcalTopology* HcalDigisValidation::htopology
private

Definition at line 150 of file HcalDigisValidation.h.

Referenced by dqmBeginRun().

edm::InputTag HcalDigisValidation::inputTag_
private

Definition at line 114 of file HcalDigisValidation.h.

Referenced by HcalDigisValidation().

int HcalDigisValidation::maxDepth_[5]
private

Definition at line 152 of file HcalDigisValidation.h.

Referenced by booking(), dqmBeginRun(), and reco().

std::string HcalDigisValidation::mc_
private

Definition at line 120 of file HcalDigisValidation.h.

Referenced by bookHistograms(), HcalDigisValidation(), and reco().

std::string HcalDigisValidation::mode_
private

Definition at line 119 of file HcalDigisValidation.h.

Referenced by HcalDigisValidation(), and reco().

std::map<std::string, MonitorElement*>* HcalDigisValidation::msm_
private

Definition at line 82 of file HcalDigisValidation.h.

Referenced by book1D(), book2D(), bookPf(), fill1D(), fill2D(), fillPf(), HcalDigisValidation(), monitor(), and ~HcalDigisValidation().

int HcalDigisValidation::nChannels_[5]
private

Definition at line 153 of file HcalDigisValidation.h.

Referenced by booking(), and dqmBeginRun().

int HcalDigisValidation::nevent1
private

Definition at line 143 of file HcalDigisValidation.h.

Referenced by HcalDigisValidation(), and reco().

int HcalDigisValidation::nevent2
private

Definition at line 144 of file HcalDigisValidation.h.

Referenced by HcalDigisValidation(), and reco().

int HcalDigisValidation::nevent3
private

Definition at line 145 of file HcalDigisValidation.h.

Referenced by HcalDigisValidation(), and reco().

int HcalDigisValidation::nevent4
private

Definition at line 146 of file HcalDigisValidation.h.

Referenced by HcalDigisValidation(), and reco().

int HcalDigisValidation::nevtot
private

Definition at line 147 of file HcalDigisValidation.h.

Referenced by analyze(), and HcalDigisValidation().

int HcalDigisValidation::noise_
private

Definition at line 121 of file HcalDigisValidation.h.

Referenced by analyze(), and reco().

std::string HcalDigisValidation::outputFile_
private

Definition at line 109 of file HcalDigisValidation.h.

Referenced by HcalDigisValidation().

edm::InputTag HcalDigisValidation::QIE10inputTag_
private

Definition at line 115 of file HcalDigisValidation.h.

Referenced by HcalDigisValidation().

edm::InputTag HcalDigisValidation::QIE11inputTag_
private

Definition at line 116 of file HcalDigisValidation.h.

Referenced by HcalDigisValidation().

bool HcalDigisValidation::skipDataTPs
private

Definition at line 155 of file HcalDigisValidation.h.

Referenced by analyze(), bookHistograms(), and HcalDigisValidation().

std::string HcalDigisValidation::subdet_
private

Definition at line 110 of file HcalDigisValidation.h.

Referenced by analyze(), bookHistograms(), HcalDigisValidation(), and reco().

bool HcalDigisValidation::testNumber_
private

Definition at line 122 of file HcalDigisValidation.h.

Referenced by HcalDigisValidation(), and reco().

edm::EDGetTokenT<HcalTrigPrimDigiCollection> HcalDigisValidation::tok_dataTPs_
private

Definition at line 130 of file HcalDigisValidation.h.

Referenced by analyze(), and HcalDigisValidation().

edm::EDGetTokenT<HcalTrigPrimDigiCollection> HcalDigisValidation::tok_emulTPs_
private

Definition at line 129 of file HcalDigisValidation.h.

Referenced by analyze(), and HcalDigisValidation().

edm::EDGetTokenT< HBHEDigiCollection > HcalDigisValidation::tok_hbhe_
private

Definition at line 126 of file HcalDigisValidation.h.

Referenced by analyze(), and HcalDigisValidation().

edm::EDGetTokenT< HFDigiCollection > HcalDigisValidation::tok_hf_
private

Definition at line 128 of file HcalDigisValidation.h.

Referenced by analyze(), and HcalDigisValidation().

edm::EDGetTokenT< HODigiCollection > HcalDigisValidation::tok_ho_
private

Definition at line 127 of file HcalDigisValidation.h.

Referenced by analyze(), and HcalDigisValidation().

edm::EDGetTokenT<edm::PCaloHitContainer> HcalDigisValidation::tok_mc_
private

Definition at line 125 of file HcalDigisValidation.h.

Referenced by HcalDigisValidation(), and reco().

edm::EDGetTokenT< QIE10DigiCollection > HcalDigisValidation::tok_qie10_hf_
private

Definition at line 132 of file HcalDigisValidation.h.

Referenced by analyze(), and HcalDigisValidation().

edm::EDGetTokenT< QIE11DigiCollection > HcalDigisValidation::tok_qie11_hbhe_
private

Definition at line 133 of file HcalDigisValidation.h.

Referenced by analyze(), and HcalDigisValidation().

std::string HcalDigisValidation::zside_
private

Definition at line 111 of file HcalDigisValidation.h.

Generated for CMSSW Reference Manual by doxygen 1.8.11