CMS 3D CMS Logo

List of all members | Public Member Functions | Private Attributes
cms::Analyzer_minbias Class Reference
Inheritance diagram for cms::Analyzer_minbias:
edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::SharedResources > edm::one::EDAnalyzerBase edm::EDConsumerBase

Public Member Functions

void analyze (edm::Event const &, edm::EventSetup const &) override
 
 Analyzer_minbias (const edm::ParameterSet &)
 
void beginJob () override
 
void beginRun (edm::Run const &, edm::EventSetup const &) override
 
void endJob () override
 
void endRun (edm::Run const &, edm::EventSetup const &) override
 
 ~Analyzer_minbias () override
 
- Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::SharedResources >
 EDAnalyzer ()=default
 
 EDAnalyzer (const EDAnalyzer &)=delete
 
SerialTaskQueueglobalLuminosityBlocksQueue () final
 
SerialTaskQueueglobalRunsQueue () final
 
const EDAnalyzeroperator= (const EDAnalyzer &)=delete
 
bool wantsGlobalLuminosityBlocks () const final
 
bool wantsGlobalRuns () const final
 
bool wantsInputProcessBlocks () const final
 
bool wantsProcessBlocks () const final
 
- Public Member Functions inherited from edm::one::EDAnalyzerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzerBase ()
 
ModuleDescription const & moduleDescription () const
 
bool wantsStreamLuminosityBlocks () const
 
bool wantsStreamRuns () const
 
 ~EDAnalyzerBase () override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfoconsumesInfo () 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
 
ESResolverIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
std::vector< ESResolverIndex > const & esGetTokenIndicesVector (edm::Transition iTrans) const
 
std::vector< ESRecordIndex > const & esGetTokenRecordIndicesVector (edm::Transition iTrans) const
 
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::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBaseoperator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
void selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Private Attributes

int depth
 
float eta
 
std::string fOutputFileName_
 
TH1F * hbheNoiseE
 
TH1F * hbheSignalE
 
std::string hcalfile_
 
TH1F * hCalo1 [73][43]
 
TH1F * hCalo1mom2 [73][43]
 
TH1F * hCalo2 [73][43]
 
TH1F * hCalo2mom2 [73][43]
 
TH1F * hfNoiseE
 
TH1F * hfSignalE
 
TH2F * hHBHEsize_vs_run
 
TH2F * hHFsize_vs_run
 
TFile * hOutputFile
 
int ieta
 
int iphi
 
double meannoise_min [73][43]
 
double meannoise_pl [73][43]
 
float mom0_Diff
 
float mom0_MB
 
float mom0_Noise
 
float mom1_Diff
 
float mom1_MB
 
float mom1_Noise
 
float mom2_Diff
 
float mom2_MB
 
float mom2_Noise
 
float mom3_Diff
 
float mom3_MB
 
float mom3_Noise
 
float mom4_Diff
 
float mom4_MB
 
float mom4_Noise
 
int mydet
 
std::ofstream * myout_hcal
 
int mysubd
 
TTree * myTree
 
double nevent
 
int nevent_run
 
double noise_min [73][43]
 
double noise_pl [73][43]
 
float occup
 
float phi
 
double theDFFillDetMapMin0 [5][5][73][43]
 
double theDFFillDetMapMin1 [5][5][73][43]
 
double theDFFillDetMapMin2 [5][5][73][43]
 
double theDFFillDetMapPl0 [5][5][73][43]
 
double theDFFillDetMapPl1 [5][5][73][43]
 
double theDFFillDetMapPl2 [5][5][73][43]
 
double theMBFillDetMapMin0 [5][5][73][43]
 
double theMBFillDetMapMin1 [5][5][73][43]
 
double theMBFillDetMapMin2 [5][5][73][43]
 
double theMBFillDetMapMin4 [5][5][73][43]
 
double theMBFillDetMapPl0 [5][5][73][43]
 
double theMBFillDetMapPl1 [5][5][73][43]
 
double theMBFillDetMapPl2 [5][5][73][43]
 
double theMBFillDetMapPl4 [5][5][73][43]
 
double theNSFillDetMapMin0 [5][5][73][43]
 
double theNSFillDetMapMin1 [5][5][73][43]
 
double theNSFillDetMapMin2 [5][5][73][43]
 
double theNSFillDetMapMin4 [5][5][73][43]
 
double theNSFillDetMapPl0 [5][5][73][43]
 
double theNSFillDetMapPl1 [5][5][73][43]
 
double theNSFillDetMapPl2 [5][5][73][43]
 
double theNSFillDetMapPl4 [5][5][73][43]
 
bool theRecalib_
 
const edm::EDGetTokenT< FEDRawDataCollectiontok_data_
 
const edm::EDGetTokenT< L1GlobalTriggerReadoutRecordtok_gtRec_
 
const edm::EDGetTokenT< HBHERecHitCollectiontok_hbhe_
 
const edm::EDGetTokenT< HBHERecHitCollectiontok_hbheNoise_
 
const edm::EDGetTokenT< HBHERecHitCollectiontok_hbheNorm_
 
const edm::EDGetTokenT< HFRecHitCollectiontok_hf_
 
const edm::EDGetTokenT< HFRecHitCollectiontok_hfNoise_
 
const edm::EDGetTokenT< HORecHitCollectiontok_ho_
 
const edm::EDGetTokenT< HORecHitCollectiontok_hoNoise_
 
const edm::ESGetToken< L1GtTriggerMenu, L1GtTriggerMenuRcdtok_l1gt_
 
const edm::ESGetToken< HcalRespCorrs, HcalRespCorrsRcdtok_respCorr_
 

Additional Inherited Members

- Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::one::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)
 
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< Bconsumes (edm::InputTag tag) noexcept
 
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<Transition Tr = Transition::Event>
constexpr auto esConsumes ()
 
template<Transition Tr = Transition::Event>
auto esConsumes (ESInputTag tag)
 
template<Transition Tr = Transition::Event>
ESGetTokenGeneric esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
 Used with EventSetupRecord::doGet. More...
 
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)
 
void resetItemsToGetFrom (BranchType iType)
 

Detailed Description

Definition at line 69 of file Analyzer_minbias.cc.

Constructor & Destructor Documentation

◆ Analyzer_minbias()

cms::Analyzer_minbias::Analyzer_minbias ( const edm::ParameterSet iConfig)
explicit

Definition at line 167 of file Analyzer_minbias.cc.

References mps_fire::i, dqmiolumiharvest::j, TFileService::kSharedResource, noise_min, and noise_pl.

168  : fOutputFileName_(iConfig.getUntrackedParameter<std::string>("HistOutFile")),
169  theRecalib_(iConfig.getParameter<bool>("Recalib")),
170  tok_hbhe_(consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("hbheInputMB"))),
171  tok_ho_(consumes<HORecHitCollection>(iConfig.getParameter<edm::InputTag>("hoInputMB"))),
172  tok_hf_(consumes<HFRecHitCollection>(iConfig.getParameter<edm::InputTag>("hfInputMB"))),
173  tok_data_(consumes<FEDRawDataCollection>(edm::InputTag(iConfig.getParameter<std::string>("InputLabel")))),
174  tok_hbheNoise_(consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("hbheInputNoise"))),
175  tok_hoNoise_(consumes<HORecHitCollection>(iConfig.getParameter<edm::InputTag>("hoInputNoise"))),
176  tok_hfNoise_(consumes<HFRecHitCollection>(iConfig.getParameter<edm::InputTag>("hfInputNoise"))),
177  tok_gtRec_(consumes<L1GlobalTriggerReadoutRecord>(edm::InputTag("gtDigisAlCaMB"))),
178  tok_hbheNorm_(consumes<HBHERecHitCollection>(edm::InputTag("hbhereco"))),
179  tok_respCorr_(esConsumes<HcalRespCorrs, HcalRespCorrsRcd>()),
180  tok_l1gt_(esConsumes<L1GtTriggerMenu, L1GtTriggerMenuRcd>()) {
181  usesResource(TFileService::kSharedResource);
182  // get name of output file with histogramms
183  // get names of modules, producing object collections
184  // some of the label names are hardcodded..
185  //
186  for (int i = 0; i < 73; i++) {
187  for (int j = 0; j < 43; j++) {
188  noise_min[i][j] = 0.;
189  noise_pl[i][j] = 0.;
190  }
191  }
192  }
static const std::string kSharedResource
Definition: TFileService.h:76
const edm::EDGetTokenT< HFRecHitCollection > tok_hfNoise_
T getParameter(std::string const &) const
Definition: ParameterSet.h:307
const edm::ESGetToken< L1GtTriggerMenu, L1GtTriggerMenuRcd > tok_l1gt_
const edm::EDGetTokenT< HBHERecHitCollection > tok_hbheNorm_
const edm::EDGetTokenT< HBHERecHitCollection > tok_hbheNoise_
const edm::EDGetTokenT< L1GlobalTriggerReadoutRecord > tok_gtRec_
const edm::EDGetTokenT< HFRecHitCollection > tok_hf_
T getUntrackedParameter(std::string const &, T const &) const
const edm::EDGetTokenT< HORecHitCollection > tok_ho_
const edm::EDGetTokenT< HBHERecHitCollection > tok_hbhe_
const edm::EDGetTokenT< FEDRawDataCollection > tok_data_
const edm::ESGetToken< HcalRespCorrs, HcalRespCorrsRcd > tok_respCorr_
const edm::EDGetTokenT< HORecHitCollection > tok_hoNoise_

◆ ~Analyzer_minbias()

cms::Analyzer_minbias::~Analyzer_minbias ( )
override

Definition at line 194 of file Analyzer_minbias.cc.

194  {
195  // do anything here that needs to be done at desctruction time
196  // (e.g. close files, deallocate resources etc.)
197  }

Member Function Documentation

◆ analyze()

void cms::Analyzer_minbias::analyze ( edm::Event const &  iEvent,
edm::EventSetup const &  iSetup 
)
overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 421 of file Analyzer_minbias.cc.

References funct::abs(), edm::SortedCollection< T, SORT >::begin(), HcalDetId::depth(), edm::SortedCollection< T, SORT >::end(), CaloRecHit::energy(), HcalObjRepresent::Fill(), nano_mu_digi_cff::float, edm::EventSetup::getData(), HcalRespCorr::getValue(), HcalCondObjectContainer< Item >::getValues(), hbheNoiseE, hbheSignalE, hCalo1, hCalo1mom2, hCalo2, hCalo2mom2, hfNoiseE, hfSignalE, hHBHEsize_vs_run, hHFsize_vs_run, mps_fire::i, HcalDetId::ieta(), iEvent, HcalDetId::iphi(), edm::HandleBase::isValid(), dqmiolumiharvest::j, dqmdumpme::k, MainPageGenerator::l, nevent, nevent_run, funct::pow(), edm::Handle< T >::product(), crab_script::provenance, DetId::rawId(), nano_mu_digi_cff::rawId, edm::SortedCollection< T, SORT >::size(), DetId::subdetId(), theDFFillDetMapMin1, theDFFillDetMapMin2, theDFFillDetMapPl1, theDFFillDetMapPl2, theMBFillDetMapMin0, theMBFillDetMapMin1, theMBFillDetMapMin2, theMBFillDetMapMin4, theMBFillDetMapPl0, theMBFillDetMapPl1, theMBFillDetMapPl2, theMBFillDetMapPl4, theNSFillDetMapMin0, theNSFillDetMapMin1, theNSFillDetMapMin2, theNSFillDetMapMin4, theNSFillDetMapPl0, theNSFillDetMapPl1, theNSFillDetMapPl2, theNSFillDetMapPl4, theRecalib_, tok_hbhe_, tok_hbheNoise_, tok_hbheNorm_, tok_hf_, tok_hfNoise_, and tok_respCorr_.

421  {
422  edm::LogVerbatim("AnalyzerMB") << " Start Analyzer_minbias::analyze " << nevent;
423  nevent++;
424  nevent_run++;
425 
426  float rnnum = (float)iEvent.run();
427 
428  std::vector<edm::StableProvenance const*> theProvenance;
429  iEvent.getAllStableProvenance(theProvenance);
430 
431  for (auto const& provenance : theProvenance) {
432  edm::LogVerbatim("AnalyzerMB") << " Print all process/modulelabel/product names " << provenance->processName()
433  << " , " << provenance->moduleLabel() << " , "
434  << provenance->productInstanceName();
435  }
436  const HcalRespCorrs* myRecalib = nullptr;
437  if (theRecalib_) {
438  myRecalib = &iSetup.getData(tok_respCorr_);
439  } // theRecalib_
440 
441  // Noise part for HB HE
442 
443  double tmpNSFillDetMapPl1[5][5][73][43];
444  double tmpNSFillDetMapMin1[5][5][73][43];
445 
446  for (int i = 0; i < 5; i++) {
447  for (int j = 0; j < 5; j++) {
448  for (int k = 0; k < 73; k++) {
449  for (int l = 0; l < 43; l++) {
450  tmpNSFillDetMapPl1[i][j][k][l] = 0.;
451  tmpNSFillDetMapMin1[i][j][k][l] = 0.;
452  }
453  }
454  }
455  }
456 
457  const edm::Handle<HBHERecHitCollection> hbheNormal = iEvent.getHandle(tok_hbheNorm_);
458  if (!hbheNormal.isValid()) {
459  edm::LogWarning("AnalyzerMB") << " hbheNormal failed ";
460  } else {
461  edm::LogVerbatim("AnalyzerMB") << " The size of the normal collection " << hbheNormal->size();
462  }
463 
464  const edm::Handle<HBHERecHitCollection> hbheNS = iEvent.getHandle(tok_hbheNoise_);
465 
466  if (!hbheNS.isValid()) {
467  edm::LogWarning("AnalyzerMB") << "HcalCalibAlgos: Error! can't get hbhe"
468  << " product! No HBHE MS ";
469  return;
470  }
471 
472  const HBHERecHitCollection HithbheNS = *(hbheNS.product());
473  edm::LogVerbatim("AnalyzerMB") << " HBHE NS size of collection " << HithbheNS.size();
474  hHBHEsize_vs_run->Fill(rnnum, (float)HithbheNS.size());
475 
476  if (HithbheNS.size() != 5184) {
477  edm::LogWarning("AnalyzerMB") << " HBHE problem " << rnnum << " " << HithbheNS.size();
478  }
479  const edm::Handle<HBHERecHitCollection> hbheMB = iEvent.getHandle(tok_hbhe_);
480 
481  if (!hbheMB.isValid()) {
482  edm::LogWarning("AnalyzerMB") << "HcalCalibAlgos: Error! can't get hbhe"
483  << " product! No HBHE MB";
484  }
485 
486  const HBHERecHitCollection HithbheMB = *(hbheMB.product());
487  edm::LogVerbatim("AnalyzerMB") << " HBHE MB size of collection " << HithbheMB.size();
488  if (HithbheMB.size() != 5184) {
489  edm::LogWarning("AnalyzerMB") << " HBHE problem " << rnnum << " " << HithbheMB.size();
490  }
491 
492  const edm::Handle<HFRecHitCollection> hfNS = iEvent.getHandle(tok_hfNoise_);
493 
494  if (!hfNS.isValid()) {
495  edm::LogWarning("AnalyzerMB") << "HcalCalibAlgos: Error! can't get hf"
496  << " product! No HF NS ";
497  }
498 
499  const HFRecHitCollection HithfNS = *(hfNS.product());
500  edm::LogVerbatim("AnalyzerMB") << " HFE NS size of collection " << HithfNS.size();
501  hHFsize_vs_run->Fill(rnnum, (float)HithfNS.size());
502  if (HithfNS.size() != 1728) {
503  edm::LogWarning("AnalyzerMB") << " HF problem " << rnnum << " " << HithfNS.size();
504  }
505 
506  const edm::Handle<HFRecHitCollection> hfMB = iEvent.getHandle(tok_hf_);
507 
508  if (!hfMB.isValid()) {
509  edm::LogWarning("AnalyzerMB") << "HcalCalibAlgos: Error! can't get hf"
510  << " product! No HF MB";
511  }
512 
513  const HFRecHitCollection HithfMB = *(hfMB.product());
514  edm::LogVerbatim("AnalyzerMB") << " HF MB size of collection " << HithfMB.size();
515  if (HithfMB.size() != 1728) {
516  edm::LogWarning("AnalyzerMB") << " HF problem " << rnnum << " " << HithfMB.size();
517  }
518 
519  for (HBHERecHitCollection::const_iterator hbheItr = HithbheNS.begin(); hbheItr != HithbheNS.end(); hbheItr++) {
520  // Recalibration of energy
521  float icalconst = 1.;
522  DetId mydetid = hbheItr->id().rawId();
523  if (theRecalib_)
524  icalconst = myRecalib->getValues(mydetid)->getValue();
525 
526  HBHERecHit aHit(hbheItr->id(), hbheItr->energy() * icalconst, hbheItr->time());
527 
528  double energyhit = aHit.energy();
529 
530  DetId id = (*hbheItr).detid();
531  HcalDetId hid = HcalDetId(id);
532 
533  int mysu = ((hid).rawId() >> 25) & 0x7;
534  if (hid.ieta() > 0) {
535  theNSFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
536  theNSFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()] + 1.;
537  theNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
538  theNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] + energyhit;
539  theNSFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
540  theNSFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] + pow(energyhit, 2);
541  theNSFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
542  theNSFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()] + pow(energyhit, 4);
543 
544  tmpNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] = energyhit;
545 
546  } else {
547  theNSFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
548  theNSFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + 1.;
549  theNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
550  theNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + energyhit;
551  theNSFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
552  theNSFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + pow(energyhit, 2);
553  theNSFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
554  theNSFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + pow(energyhit, 4);
555 
556  tmpNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()] = energyhit;
557  }
558 
559  if (hid.depth() == 1) {
560  hbheNoiseE->Fill(energyhit);
561 
562  if (energyhit < -2.)
563  edm::LogVerbatim("AnalyzerMB") << " Run " << rnnum << " ieta,iphi " << hid.ieta() << " " << hid.iphi()
564  << energyhit;
565 
566  } // depth=1
567 
568  } // HBHE_NS
569 
570  // Signal part for HB HE
571 
572  for (HBHERecHitCollection::const_iterator hbheItr = HithbheMB.begin(); hbheItr != HithbheMB.end(); hbheItr++) {
573  // Recalibration of energy
574  float icalconst = 1.;
575  DetId mydetid = hbheItr->id().rawId();
576  if (theRecalib_)
577  icalconst = myRecalib->getValues(mydetid)->getValue();
578 
579  HBHERecHit aHit(hbheItr->id(), hbheItr->energy() * icalconst, hbheItr->time());
580 
581  double energyhit = aHit.energy();
582 
583  DetId id = (*hbheItr).detid();
584  HcalDetId hid = HcalDetId(id);
585 
586  int mysu = ((hid).rawId() >> 25) & 0x7;
587  if (hid.ieta() > 0) {
588  theMBFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()] += 1.;
589  theMBFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] += energyhit;
590  theMBFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] += pow(energyhit, 2);
591  theMBFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()] += pow(energyhit, 4);
592  float mydiff = energyhit - tmpNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()];
593 
594  theDFFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
595  theDFFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] + mydiff;
596  theDFFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
597  theDFFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] + pow(mydiff, 2);
598  } else {
599  theMBFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
600  theMBFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + 1.;
601  theMBFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
602  theMBFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + energyhit;
603  theMBFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
604  theMBFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + pow(energyhit, 2);
605  theMBFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
606  theMBFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + pow(energyhit, 4);
607 
608  float mydiff = energyhit - tmpNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()];
609  theDFFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
610  theDFFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()] + mydiff;
611  theDFFillDetMapMin2[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
612  theDFFillDetMapMin2[mysu][hid.depth()][hid.iphi()][hid.ieta()] + pow(mydiff, 2);
613  }
614 
615  if (hid.depth() == 1) {
616  hbheSignalE->Fill(energyhit);
617 
618  if (hid.ieta() > 0) {
619  hCalo1[hid.iphi()][hid.ieta()]->Fill(energyhit);
620  hCalo1mom2[hid.iphi()][hid.ieta()]->Fill(pow(energyhit, 2));
621  } else {
622  hCalo2[hid.iphi()][abs(hid.ieta())]->Fill(energyhit);
623  hCalo2mom2[hid.iphi()][abs(hid.ieta())]->Fill(pow(energyhit, 2));
624  } // eta><0
625 
626  } // depth=1
627 
628  } // HBHE_MB
629 
630  // HF
631 
632  for (HFRecHitCollection::const_iterator hbheItr = HithfNS.begin(); hbheItr != HithfNS.end(); hbheItr++) {
633  // Recalibration of energy
634  float icalconst = 1.;
635  DetId mydetid = hbheItr->id().rawId();
636  if (theRecalib_)
637  icalconst = myRecalib->getValues(mydetid)->getValue();
638 
639  HFRecHit aHit(hbheItr->id(), hbheItr->energy() * icalconst, hbheItr->time());
640 
641  double energyhit = aHit.energy();
642  //
643  // Remove PMT hits
644  //
645  DetId id = (*hbheItr).detid();
646  HcalDetId hid = HcalDetId(id);
647 
648  if (fabs(energyhit) > 40.)
649  continue;
650 
651  int mysu = hid.subdetId();
652  if (hid.ieta() > 0) {
653  theNSFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
654  theNSFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()] + 1.;
655  theNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
656  theNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] + energyhit;
657  theNSFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
658  theNSFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] + pow(energyhit, 2);
659  theNSFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
660  theNSFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()] + pow(energyhit, 4);
661 
662  tmpNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] = energyhit;
663 
664  } else {
665  theNSFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
666  theNSFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + 1.;
667  theNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
668  theNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + energyhit;
669  theNSFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
670  theNSFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + pow(energyhit, 2);
671  theNSFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
672  theNSFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + pow(energyhit, 4);
673 
674  tmpNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()] = energyhit;
675  }
676 
677  if (hid.depth() == 1) {
678  hfNoiseE->Fill(energyhit);
679 
680  } // depth=1
681 
682  } // HBHE_NS
683 
684  // Signal part for HB HE
685 
686  for (HFRecHitCollection::const_iterator hbheItr = HithfMB.begin(); hbheItr != HithfMB.end(); hbheItr++) {
687  // Recalibration of energy
688  float icalconst = 1.;
689  DetId mydetid = hbheItr->id().rawId();
690  if (theRecalib_)
691  icalconst = myRecalib->getValues(mydetid)->getValue();
692 
693  HFRecHit aHit(hbheItr->id(), hbheItr->energy() * icalconst, hbheItr->time());
694 
695  double energyhit = aHit.energy();
696  //
697  // Remove PMT hits
698  //
699  if (fabs(energyhit) > 40.)
700  continue;
701 
702  DetId id = (*hbheItr).detid();
703  HcalDetId hid = HcalDetId(id);
704 
705  int mysu = ((hid).rawId() >> 25) & 0x7;
706  if (hid.ieta() > 0) {
707  theMBFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
708  theMBFillDetMapPl0[mysu][hid.depth()][hid.iphi()][hid.ieta()] + 1.;
709  theMBFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
710  theMBFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] + energyhit;
711  theMBFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
712  theMBFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] + pow(energyhit, 2);
713  theMBFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
714  theMBFillDetMapPl4[mysu][hid.depth()][hid.iphi()][hid.ieta()] + pow(energyhit, 4);
715 
716  theDFFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
717  theDFFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()] + energyhit -
718  tmpNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()];
719  theDFFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
720  theDFFillDetMapPl2[mysu][hid.depth()][hid.iphi()][hid.ieta()] +
721  pow((energyhit - tmpNSFillDetMapPl1[mysu][hid.depth()][hid.iphi()][hid.ieta()]), 2);
722  } else {
723  theMBFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
724  theMBFillDetMapMin0[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + 1.;
725  theMBFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
726  theMBFillDetMapMin1[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + energyhit;
727  theMBFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
728  theMBFillDetMapMin2[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + pow(energyhit, 2);
729  theMBFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] =
730  theMBFillDetMapMin4[mysu][hid.depth()][hid.iphi()][abs(hid.ieta())] + pow(energyhit, 4);
731 
732  theDFFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
733  theDFFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()] + energyhit -
734  tmpNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()];
735  theDFFillDetMapMin2[mysu][hid.depth()][hid.iphi()][hid.ieta()] =
736  theDFFillDetMapMin2[mysu][hid.depth()][hid.iphi()][hid.ieta()] +
737  pow((energyhit - tmpNSFillDetMapMin1[mysu][hid.depth()][hid.iphi()][hid.ieta()]), 2);
738  }
739 
740  if (hid.depth() == 1) {
741  hfSignalE->Fill(energyhit);
742 
743  if (hid.ieta() > 0) {
744  hCalo1[hid.iphi()][hid.ieta()]->Fill(energyhit);
745  hCalo1mom2[hid.iphi()][hid.ieta()]->Fill(pow(energyhit, 2));
746  } else {
747  hCalo2[hid.iphi()][abs(hid.ieta())]->Fill(energyhit);
748  hCalo2mom2[hid.iphi()][abs(hid.ieta())]->Fill(pow(energyhit, 2));
749  } // eta><0
750 
751  } // depth=1
752 
753  } // HF_MB
754 
755  edm::LogVerbatim("AnalyzerMB") << " Event is finished ";
756  }
double theMBFillDetMapMin0[5][5][73][43]
double theMBFillDetMapMin2[5][5][73][43]
Log< level::Info, true > LogVerbatim
double theDFFillDetMapPl2[5][5][73][43]
const edm::EDGetTokenT< HFRecHitCollection > tok_hfNoise_
double theMBFillDetMapMin1[5][5][73][43]
const edm::EDGetTokenT< HBHERecHitCollection > tok_hbheNorm_
double theNSFillDetMapPl4[5][5][73][43]
double theNSFillDetMapPl0[5][5][73][43]
double theNSFillDetMapMin0[5][5][73][43]
double theMBFillDetMapPl0[5][5][73][43]
size_type size() const
TH1F * hCalo2mom2[73][43]
const edm::EDGetTokenT< HBHERecHitCollection > tok_hbheNoise_
T const * product() const
Definition: Handle.h:70
std::vector< T >::const_iterator const_iterator
double theDFFillDetMapPl1[5][5][73][43]
double theNSFillDetMapPl2[5][5][73][43]
double theDFFillDetMapMin2[5][5][73][43]
TH1F * hCalo1mom2[73][43]
const edm::EDGetTokenT< HFRecHitCollection > tok_hf_
const Item * getValues(DetId fId, bool throwOnFail=true) const
constexpr float energy() const
Definition: CaloRecHit.h:29
double theNSFillDetMapPl1[5][5][73][43]
double theMBFillDetMapPl1[5][5][73][43]
int iEvent
Definition: GenABIO.cc:224
const edm::EDGetTokenT< HBHERecHitCollection > tok_hbhe_
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
constexpr int ieta() const
get the cell ieta
Definition: HcalDetId.h:155
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
const edm::ESGetToken< HcalRespCorrs, HcalRespCorrsRcd > tok_respCorr_
double theMBFillDetMapPl4[5][5][73][43]
constexpr int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:48
double theNSFillDetMapMin1[5][5][73][43]
double theNSFillDetMapMin4[5][5][73][43]
const_iterator begin() const
double theMBFillDetMapPl2[5][5][73][43]
const_iterator end() const
Definition: DetId.h:17
float getValue() const
Definition: HcalRespCorr.h:19
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:57
double theMBFillDetMapMin4[5][5][73][43]
bool isValid() const
Definition: HandleBase.h:70
Log< level::Warning, false > LogWarning
double theNSFillDetMapMin2[5][5][73][43]
double theDFFillDetMapMin1[5][5][73][43]
constexpr int iphi() const
get the cell iphi
Definition: HcalDetId.h:157
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:29
constexpr int depth() const
get the tower depth
Definition: HcalDetId.h:164

◆ beginJob()

void cms::Analyzer_minbias::beginJob ( )
overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 204 of file Analyzer_minbias.cc.

References depth, eta, compareTotals::fs, hbheNoiseE, hbheSignalE, hCalo1, hCalo1mom2, hCalo2, hCalo2mom2, hfNoiseE, hfSignalE, hHBHEsize_vs_run, hHFsize_vs_run, mps_fire::i, ieta, iphi, dqmiolumiharvest::j, dqmdumpme::k, MainPageGenerator::l, meannoise_min, meannoise_pl, mom0_Diff, mom0_MB, mom0_Noise, mom1_Diff, mom1_MB, mom1_Noise, mom2_Diff, mom2_MB, mom2_Noise, mom4_MB, mom4_Noise, mydet, myout_hcal, mysubd, myTree, nevent, occup, phi, AlCaHLTBitMon_QueryRunRegistry::string, theDFFillDetMapMin0, theDFFillDetMapMin1, theDFFillDetMapMin2, theDFFillDetMapPl0, theDFFillDetMapPl1, theDFFillDetMapPl2, theMBFillDetMapMin0, theMBFillDetMapMin1, theMBFillDetMapMin2, theMBFillDetMapMin4, theMBFillDetMapPl0, theMBFillDetMapPl1, theMBFillDetMapPl2, theMBFillDetMapPl4, theNSFillDetMapMin0, theNSFillDetMapMin1, theNSFillDetMapMin2, theNSFillDetMapMin4, theNSFillDetMapPl0, theNSFillDetMapPl1, theNSFillDetMapPl2, and theNSFillDetMapPl4.

204  {
206  myTree = fs->make<TTree>("RecJet", "RecJet Tree");
207  myTree->Branch("mydet", &mydet, "mydet/I");
208  myTree->Branch("mysubd", &mysubd, "mysubd/I");
209  myTree->Branch("depth", &depth, "depth/I");
210  myTree->Branch("ieta", &ieta, "ieta/I");
211  myTree->Branch("iphi", &iphi, "iphi/I");
212  myTree->Branch("eta", &eta, "eta/F");
213  myTree->Branch("phi", &phi, "phi/F");
214 
215  myTree->Branch("mom0_MB", &mom0_MB, "mom0_MB/F");
216  myTree->Branch("mom1_MB", &mom1_MB, "mom1_MB/F");
217  myTree->Branch("mom2_MB", &mom2_MB, "mom2_MB/F");
218  myTree->Branch("mom4_MB", &mom4_MB, "mom4_MB/F");
219 
220  myTree->Branch("mom0_Noise", &mom0_Noise, "mom0_Noise/F");
221  myTree->Branch("mom1_Noise", &mom1_Noise, "mom1_Noise/F");
222  myTree->Branch("mom2_Noise", &mom2_Noise, "mom2_Noise/F");
223  myTree->Branch("mom4_Noise", &mom4_Noise, "mom4_Noise/F");
224 
225  myTree->Branch("mom0_Diff", &mom0_Diff, "mom0_Diff/F");
226  myTree->Branch("mom1_Diff", &mom1_Diff, "mom1_Diff/F");
227  myTree->Branch("mom2_Diff", &mom2_Diff, "mom2_Diff/F");
228 
229  myTree->Branch("occup", &occup, "occup/F");
230 
231  edm::LogVerbatim("AnalyzerMB") << " Before ordering Histos ";
232 
233  char str0[32];
234  char str1[32];
235 
236  char str10[32];
237  char str11[32];
238 
239  int k = 0;
240  nevent = 0;
241  // Size of collections
242 
244  fs->make<TH2F>("hHBHEsize_vs_run", "hHBHEsize_vs_run", 500, 111500., 112000., 6101, -100.5, 6000.5);
245  hHFsize_vs_run = fs->make<TH2F>("hHFsize_vs_run", "hHFsize_vs_run", 500, 111500., 112000., 6101, -100.5, 6000.5);
246 
247  for (int i = 1; i < 73; i++) {
248  for (int j = 1; j < 43; j++) {
249  meannoise_pl[i][j] = 0.;
250  meannoise_min[i][j] = 0.;
251 
252  k = i * 1000 + j;
253  sprintf(str0, "mpl%d", k);
254  sprintf(str1, "mmin%d", k);
255 
256  sprintf(str10, "vpl%d", k);
257  sprintf(str11, "vmin%d", k);
258  if (j < 30) {
259  // first order moment
260  hCalo1[i][j] = fs->make<TH1F>(str0, "h0", 320, -10., 10.);
261  hCalo2[i][j] = fs->make<TH1F>(str1, "h1", 320, -10., 10.);
262 
263  // second order moment
264  hCalo1mom2[i][j] = fs->make<TH1F>(str10, "h10", 320, 0., 20.);
265  hCalo2mom2[i][j] = fs->make<TH1F>(str11, "h11", 320, 0., 20.);
266  } else {
267  // HF
268  // first order moment
269  if (j < 40) {
270  hCalo1[i][j] = fs->make<TH1F>(str0, "h0", 320, -10., 10.);
271  hCalo2[i][j] = fs->make<TH1F>(str1, "h1", 320, -10., 10.);
272  //
273  // second order moment
274  hCalo1mom2[i][j] = fs->make<TH1F>(str10, "h10", 320, 0., 40.);
275  hCalo2mom2[i][j] = fs->make<TH1F>(str11, "h11", 320, 0., 40.);
276  } else {
277  hCalo1[i][j] = fs->make<TH1F>(str0, "h0", 320, -10., 10.);
278  hCalo2[i][j] = fs->make<TH1F>(str1, "h1", 320, -10., 10.);
279 
280  // second order moment
281  hCalo1mom2[i][j] = fs->make<TH1F>(str10, "h10", 320, 0., 120.);
282  hCalo2mom2[i][j] = fs->make<TH1F>(str11, "h11", 320, 0., 120.);
283  }
284  } // HE/HF boundary
285  } // j
286  } // i
287 
288  hbheNoiseE = fs->make<TH1F>("hbheNoiseE", "hbheNoiseE", 320, -10., 10.);
289  hfNoiseE = fs->make<TH1F>("hfNoiseE", "hfNoiseE", 320, -10., 10.);
290  hbheSignalE = fs->make<TH1F>("hbheSignalE", "hbheSignalE", 320, -10., 10.);
291  hfSignalE = fs->make<TH1F>("hfSignalE", "hfSignalE", 320, -10., 10.);
292 
293  edm::LogVerbatim("AnalyzerMB") << " After ordering Histos ";
294 
295  std::string ccc = "noise_0.dat";
296 
297  myout_hcal = new std::ofstream(ccc.c_str());
298  if (!myout_hcal)
299  edm::LogVerbatim("AnalyzerMB") << " Output file not open!!! ";
300 
301  //
302  for (int i = 0; i < 5; i++) {
303  for (int j = 0; j < 5; j++) {
304  for (int k = 0; k < 73; k++) {
305  for (int l = 0; l < 43; l++) {
306  theMBFillDetMapPl0[i][j][k][l] = 0.;
307  theMBFillDetMapPl1[i][j][k][l] = 0.;
308  theMBFillDetMapPl2[i][j][k][l] = 0.;
309  theMBFillDetMapPl4[i][j][k][l] = 0.;
310 
311  theMBFillDetMapMin0[i][j][k][l] = 0.;
312  theMBFillDetMapMin1[i][j][k][l] = 0.;
313  theMBFillDetMapMin2[i][j][k][l] = 0.;
314  theMBFillDetMapMin4[i][j][k][l] = 0.;
315 
316  theNSFillDetMapPl0[i][j][k][l] = 0.;
317  theNSFillDetMapPl1[i][j][k][l] = 0.;
318  theNSFillDetMapPl2[i][j][k][l] = 0.;
319  theNSFillDetMapPl4[i][j][k][l] = 0.;
320 
321  theNSFillDetMapMin0[i][j][k][l] = 0.;
322  theNSFillDetMapMin1[i][j][k][l] = 0.;
323  theNSFillDetMapMin2[i][j][k][l] = 0.;
324  theNSFillDetMapMin4[i][j][k][l] = 0.;
325 
326  theDFFillDetMapPl0[i][j][k][l] = 0.;
327  theDFFillDetMapPl1[i][j][k][l] = 0.;
328  theDFFillDetMapPl2[i][j][k][l] = 0.;
329  theDFFillDetMapMin0[i][j][k][l] = 0.;
330  theDFFillDetMapMin1[i][j][k][l] = 0.;
331  theDFFillDetMapMin2[i][j][k][l] = 0.;
332  }
333  }
334  }
335  }
336 
337  return;
338  }
double theMBFillDetMapMin0[5][5][73][43]
double theMBFillDetMapMin2[5][5][73][43]
Log< level::Info, true > LogVerbatim
double theDFFillDetMapPl2[5][5][73][43]
double theMBFillDetMapMin1[5][5][73][43]
double theNSFillDetMapPl4[5][5][73][43]
double theNSFillDetMapPl0[5][5][73][43]
double theNSFillDetMapMin0[5][5][73][43]
double theMBFillDetMapPl0[5][5][73][43]
TH1F * hCalo2mom2[73][43]
double theDFFillDetMapPl1[5][5][73][43]
std::ofstream * myout_hcal
double meannoise_pl[73][43]
double theNSFillDetMapPl2[5][5][73][43]
double theDFFillDetMapMin2[5][5][73][43]
TH1F * hCalo1mom2[73][43]
double theNSFillDetMapPl1[5][5][73][43]
double theMBFillDetMapPl1[5][5][73][43]
double theDFFillDetMapPl0[5][5][73][43]
double theMBFillDetMapPl4[5][5][73][43]
double theNSFillDetMapMin1[5][5][73][43]
double theNSFillDetMapMin4[5][5][73][43]
double theDFFillDetMapMin0[5][5][73][43]
double theMBFillDetMapPl2[5][5][73][43]
double theMBFillDetMapMin4[5][5][73][43]
double meannoise_min[73][43]
double theNSFillDetMapMin2[5][5][73][43]
double theDFFillDetMapMin1[5][5][73][43]

◆ beginRun()

void cms::Analyzer_minbias::beginRun ( edm::Run const &  ,
edm::EventSetup const &   
)
override

Definition at line 199 of file Analyzer_minbias.cc.

References nevent_run.

◆ endJob()

void cms::Analyzer_minbias::endJob ( void  )
overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 342 of file Analyzer_minbias.cc.

References depth, hCalo1, hCalo1mom2, hCalo2, hCalo2mom2, mps_fire::i, ieta, cuy::ii, iphi, dqmiolumiharvest::j, dqmdumpme::k, MainPageGenerator::l, mom0_Diff, mom0_MB, mom0_Noise, mom1_Diff, mom1_MB, mom1_Noise, mom2_Diff, mom2_MB, mom2_Noise, mom4_MB, mom4_Noise, mysubd, myTree, theDFFillDetMapMin0, theDFFillDetMapMin1, theDFFillDetMapMin2, theDFFillDetMapPl0, theDFFillDetMapPl1, theDFFillDetMapPl2, theMBFillDetMapMin0, theMBFillDetMapMin1, theMBFillDetMapMin2, theMBFillDetMapMin4, theMBFillDetMapPl0, theMBFillDetMapPl1, theMBFillDetMapPl2, theMBFillDetMapPl4, theNSFillDetMapMin0, theNSFillDetMapMin1, theNSFillDetMapMin2, theNSFillDetMapMin4, theNSFillDetMapPl0, theNSFillDetMapPl1, theNSFillDetMapPl2, and theNSFillDetMapPl4.

Referenced by o2olib.O2ORunMgr::executeJob().

342  {
343  int ii = 0;
344 
345  for (int i = 1; i < 5; i++) {
346  for (int j = 1; j < 5; j++) {
347  for (int k = 1; k < 73; k++) {
348  for (int l = 1; l < 43; l++) {
349  if (theMBFillDetMapPl0[i][j][k][l] > 0) {
361 
362  mysubd = i;
363  depth = j;
364  ieta = l;
365  iphi = k;
366  edm::LogVerbatim("AnalyzerMB") << " Result Plus= " << mysubd << " " << ieta << " " << iphi << " mom0 "
367  << mom0_MB << " mom1 " << mom1_MB << " mom2 " << mom2_MB;
368  myTree->Fill();
369  ii++;
370  } // Pl > 0
371 
372  if (theMBFillDetMapMin0[i][j][k][l] > 0) {
384 
385  mysubd = i;
386  depth = j;
387  ieta = -1 * l;
388  iphi = k;
389  edm::LogVerbatim("AnalyzerMB") << " Result Minus= " << mysubd << " " << ieta << " " << iphi << " mom0 "
390  << mom0_MB << " mom1 " << mom1_MB << " mom2 " << mom2_MB;
391  myTree->Fill();
392  ii++;
393 
394  } // Min>0
395  } // ieta
396  } // iphi
397  } // depth
398  } //subd
399 
400  edm::LogVerbatim("AnalyzerMB") << " Number of cells " << ii;
401 
402  for (int i = 1; i < 73; i++) {
403  for (int j = 1; j < 43; j++) {
404  hCalo1[i][j]->Write();
405  hCalo2[i][j]->Write();
406  hCalo1mom2[i][j]->Write();
407  hCalo2mom2[i][j]->Write();
408  }
409  }
410 
411  edm::LogVerbatim("AnalyzerMB") << " File is closed ";
412 
413  return;
414  }
double theMBFillDetMapMin0[5][5][73][43]
double theMBFillDetMapMin2[5][5][73][43]
Log< level::Info, true > LogVerbatim
double theDFFillDetMapPl2[5][5][73][43]
double theMBFillDetMapMin1[5][5][73][43]
double theNSFillDetMapPl4[5][5][73][43]
double theNSFillDetMapPl0[5][5][73][43]
double theNSFillDetMapMin0[5][5][73][43]
double theMBFillDetMapPl0[5][5][73][43]
TH1F * hCalo2mom2[73][43]
double theDFFillDetMapPl1[5][5][73][43]
double theNSFillDetMapPl2[5][5][73][43]
double theDFFillDetMapMin2[5][5][73][43]
TH1F * hCalo1mom2[73][43]
double theNSFillDetMapPl1[5][5][73][43]
double theMBFillDetMapPl1[5][5][73][43]
double theDFFillDetMapPl0[5][5][73][43]
double theMBFillDetMapPl4[5][5][73][43]
double theNSFillDetMapMin1[5][5][73][43]
double theNSFillDetMapMin4[5][5][73][43]
double theDFFillDetMapMin0[5][5][73][43]
double theMBFillDetMapPl2[5][5][73][43]
ii
Definition: cuy.py:589
double theMBFillDetMapMin4[5][5][73][43]
double theNSFillDetMapMin2[5][5][73][43]
double theDFFillDetMapMin1[5][5][73][43]

◆ endRun()

void cms::Analyzer_minbias::endRun ( edm::Run const &  r,
edm::EventSetup const &   
)
override

Definition at line 200 of file Analyzer_minbias.cc.

References nevent_run.

200  {
201  edm::LogVerbatim("AnalyzerMB") << " Runnumber " << r.run() << " Nevents " << nevent_run;
202  }
Log< level::Info, true > LogVerbatim

Member Data Documentation

◆ depth

int cms::Analyzer_minbias::depth
private

◆ eta

float cms::Analyzer_minbias::eta
private

◆ fOutputFileName_

std::string cms::Analyzer_minbias::fOutputFileName_
private

Definition at line 82 of file Analyzer_minbias.cc.

◆ hbheNoiseE

TH1F* cms::Analyzer_minbias::hbheNoiseE
private

Definition at line 94 of file Analyzer_minbias.cc.

Referenced by analyze(), and beginJob().

◆ hbheSignalE

TH1F* cms::Analyzer_minbias::hbheSignalE
private

Definition at line 95 of file Analyzer_minbias.cc.

Referenced by analyze(), and beginJob().

◆ hcalfile_

std::string cms::Analyzer_minbias::hcalfile_
private

Definition at line 84 of file Analyzer_minbias.cc.

◆ hCalo1

TH1F* cms::Analyzer_minbias::hCalo1[73][43]
private

Definition at line 90 of file Analyzer_minbias.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ hCalo1mom2

TH1F* cms::Analyzer_minbias::hCalo1mom2[73][43]
private

Definition at line 92 of file Analyzer_minbias.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ hCalo2

TH1F* cms::Analyzer_minbias::hCalo2[73][43]
private

Definition at line 91 of file Analyzer_minbias.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ hCalo2mom2

TH1F* cms::Analyzer_minbias::hCalo2mom2[73][43]
private

Definition at line 93 of file Analyzer_minbias.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ hfNoiseE

TH1F* cms::Analyzer_minbias::hfNoiseE
private

Definition at line 96 of file Analyzer_minbias.cc.

Referenced by analyze(), and beginJob().

◆ hfSignalE

TH1F* cms::Analyzer_minbias::hfSignalE
private

Definition at line 97 of file Analyzer_minbias.cc.

Referenced by analyze(), and beginJob().

◆ hHBHEsize_vs_run

TH2F* cms::Analyzer_minbias::hHBHEsize_vs_run
private

Definition at line 99 of file Analyzer_minbias.cc.

Referenced by analyze(), and beginJob().

◆ hHFsize_vs_run

TH2F* cms::Analyzer_minbias::hHFsize_vs_run
private

Definition at line 100 of file Analyzer_minbias.cc.

Referenced by analyze(), and beginJob().

◆ hOutputFile

TFile* cms::Analyzer_minbias::hOutputFile
private

Definition at line 88 of file Analyzer_minbias.cc.

◆ ieta

int cms::Analyzer_minbias::ieta
private

Definition at line 103 of file Analyzer_minbias.cc.

Referenced by beginJob(), and endJob().

◆ iphi

int cms::Analyzer_minbias::iphi
private

Definition at line 103 of file Analyzer_minbias.cc.

Referenced by beginJob(), and endJob().

◆ meannoise_min

double cms::Analyzer_minbias::meannoise_min[73][43]
private

Definition at line 111 of file Analyzer_minbias.cc.

Referenced by beginJob().

◆ meannoise_pl

double cms::Analyzer_minbias::meannoise_pl[73][43]
private

Definition at line 111 of file Analyzer_minbias.cc.

Referenced by beginJob().

◆ mom0_Diff

float cms::Analyzer_minbias::mom0_Diff
private

Definition at line 107 of file Analyzer_minbias.cc.

Referenced by beginJob(), and endJob().

◆ mom0_MB

float cms::Analyzer_minbias::mom0_MB
private

Definition at line 105 of file Analyzer_minbias.cc.

Referenced by beginJob(), and endJob().

◆ mom0_Noise

float cms::Analyzer_minbias::mom0_Noise
private

Definition at line 106 of file Analyzer_minbias.cc.

Referenced by beginJob(), and endJob().

◆ mom1_Diff

float cms::Analyzer_minbias::mom1_Diff
private

Definition at line 107 of file Analyzer_minbias.cc.

Referenced by beginJob(), and endJob().

◆ mom1_MB

float cms::Analyzer_minbias::mom1_MB
private

Definition at line 105 of file Analyzer_minbias.cc.

Referenced by beginJob(), and endJob().

◆ mom1_Noise

float cms::Analyzer_minbias::mom1_Noise
private

Definition at line 106 of file Analyzer_minbias.cc.

Referenced by beginJob(), and endJob().

◆ mom2_Diff

float cms::Analyzer_minbias::mom2_Diff
private

Definition at line 107 of file Analyzer_minbias.cc.

Referenced by beginJob(), and endJob().

◆ mom2_MB

float cms::Analyzer_minbias::mom2_MB
private

Definition at line 105 of file Analyzer_minbias.cc.

Referenced by beginJob(), and endJob().

◆ mom2_Noise

float cms::Analyzer_minbias::mom2_Noise
private

Definition at line 106 of file Analyzer_minbias.cc.

Referenced by beginJob(), and endJob().

◆ mom3_Diff

float cms::Analyzer_minbias::mom3_Diff
private

Definition at line 107 of file Analyzer_minbias.cc.

◆ mom3_MB

float cms::Analyzer_minbias::mom3_MB
private

Definition at line 105 of file Analyzer_minbias.cc.

◆ mom3_Noise

float cms::Analyzer_minbias::mom3_Noise
private

Definition at line 106 of file Analyzer_minbias.cc.

◆ mom4_Diff

float cms::Analyzer_minbias::mom4_Diff
private

Definition at line 107 of file Analyzer_minbias.cc.

◆ mom4_MB

float cms::Analyzer_minbias::mom4_MB
private

Definition at line 105 of file Analyzer_minbias.cc.

Referenced by beginJob(), and endJob().

◆ mom4_Noise

float cms::Analyzer_minbias::mom4_Noise
private

Definition at line 106 of file Analyzer_minbias.cc.

Referenced by beginJob(), and endJob().

◆ mydet

int cms::Analyzer_minbias::mydet
private

Definition at line 103 of file Analyzer_minbias.cc.

Referenced by beginJob().

◆ myout_hcal

std::ofstream* cms::Analyzer_minbias::myout_hcal
private

Definition at line 85 of file Analyzer_minbias.cc.

Referenced by beginJob().

◆ mysubd

int cms::Analyzer_minbias::mysubd
private

Definition at line 103 of file Analyzer_minbias.cc.

Referenced by beginJob(), and endJob().

◆ myTree

TTree* cms::Analyzer_minbias::myTree
private

Definition at line 89 of file Analyzer_minbias.cc.

Referenced by beginJob(), and endJob().

◆ nevent

double cms::Analyzer_minbias::nevent
private

Definition at line 116 of file Analyzer_minbias.cc.

Referenced by analyze(), and beginJob().

◆ nevent_run

int cms::Analyzer_minbias::nevent_run
private

Definition at line 102 of file Analyzer_minbias.cc.

Referenced by analyze(), beginRun(), and endRun().

◆ noise_min

double cms::Analyzer_minbias::noise_min[73][43]
private

Definition at line 112 of file Analyzer_minbias.cc.

Referenced by Analyzer_minbias().

◆ noise_pl

double cms::Analyzer_minbias::noise_pl[73][43]
private

Definition at line 112 of file Analyzer_minbias.cc.

Referenced by Analyzer_minbias().

◆ occup

float cms::Analyzer_minbias::occup
private

Definition at line 105 of file Analyzer_minbias.cc.

Referenced by beginJob().

◆ phi

float cms::Analyzer_minbias::phi
private

◆ theDFFillDetMapMin0

double cms::Analyzer_minbias::theDFFillDetMapMin0[5][5][73][43]
private

Definition at line 140 of file Analyzer_minbias.cc.

Referenced by beginJob(), and endJob().

◆ theDFFillDetMapMin1

double cms::Analyzer_minbias::theDFFillDetMapMin1[5][5][73][43]
private

Definition at line 141 of file Analyzer_minbias.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ theDFFillDetMapMin2

double cms::Analyzer_minbias::theDFFillDetMapMin2[5][5][73][43]
private

Definition at line 142 of file Analyzer_minbias.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ theDFFillDetMapPl0

double cms::Analyzer_minbias::theDFFillDetMapPl0[5][5][73][43]
private

Definition at line 137 of file Analyzer_minbias.cc.

Referenced by beginJob(), and endJob().

◆ theDFFillDetMapPl1

double cms::Analyzer_minbias::theDFFillDetMapPl1[5][5][73][43]
private

Definition at line 138 of file Analyzer_minbias.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ theDFFillDetMapPl2

double cms::Analyzer_minbias::theDFFillDetMapPl2[5][5][73][43]
private

Definition at line 139 of file Analyzer_minbias.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ theMBFillDetMapMin0

double cms::Analyzer_minbias::theMBFillDetMapMin0[5][5][73][43]
private

Definition at line 122 of file Analyzer_minbias.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ theMBFillDetMapMin1

double cms::Analyzer_minbias::theMBFillDetMapMin1[5][5][73][43]
private

Definition at line 123 of file Analyzer_minbias.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ theMBFillDetMapMin2

double cms::Analyzer_minbias::theMBFillDetMapMin2[5][5][73][43]
private

Definition at line 124 of file Analyzer_minbias.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ theMBFillDetMapMin4

double cms::Analyzer_minbias::theMBFillDetMapMin4[5][5][73][43]
private

Definition at line 125 of file Analyzer_minbias.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ theMBFillDetMapPl0

double cms::Analyzer_minbias::theMBFillDetMapPl0[5][5][73][43]
private

Definition at line 117 of file Analyzer_minbias.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ theMBFillDetMapPl1

double cms::Analyzer_minbias::theMBFillDetMapPl1[5][5][73][43]
private

Definition at line 118 of file Analyzer_minbias.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ theMBFillDetMapPl2

double cms::Analyzer_minbias::theMBFillDetMapPl2[5][5][73][43]
private

Definition at line 119 of file Analyzer_minbias.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ theMBFillDetMapPl4

double cms::Analyzer_minbias::theMBFillDetMapPl4[5][5][73][43]
private

Definition at line 120 of file Analyzer_minbias.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ theNSFillDetMapMin0

double cms::Analyzer_minbias::theNSFillDetMapMin0[5][5][73][43]
private

Definition at line 132 of file Analyzer_minbias.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ theNSFillDetMapMin1

double cms::Analyzer_minbias::theNSFillDetMapMin1[5][5][73][43]
private

Definition at line 133 of file Analyzer_minbias.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ theNSFillDetMapMin2

double cms::Analyzer_minbias::theNSFillDetMapMin2[5][5][73][43]
private

Definition at line 134 of file Analyzer_minbias.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ theNSFillDetMapMin4

double cms::Analyzer_minbias::theNSFillDetMapMin4[5][5][73][43]
private

Definition at line 135 of file Analyzer_minbias.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ theNSFillDetMapPl0

double cms::Analyzer_minbias::theNSFillDetMapPl0[5][5][73][43]
private

Definition at line 127 of file Analyzer_minbias.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ theNSFillDetMapPl1

double cms::Analyzer_minbias::theNSFillDetMapPl1[5][5][73][43]
private

Definition at line 128 of file Analyzer_minbias.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ theNSFillDetMapPl2

double cms::Analyzer_minbias::theNSFillDetMapPl2[5][5][73][43]
private

Definition at line 129 of file Analyzer_minbias.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ theNSFillDetMapPl4

double cms::Analyzer_minbias::theNSFillDetMapPl4[5][5][73][43]
private

Definition at line 130 of file Analyzer_minbias.cc.

Referenced by analyze(), beginJob(), and endJob().

◆ theRecalib_

bool cms::Analyzer_minbias::theRecalib_
private

Definition at line 83 of file Analyzer_minbias.cc.

Referenced by analyze().

◆ tok_data_

const edm::EDGetTokenT<FEDRawDataCollection> cms::Analyzer_minbias::tok_data_
private

Definition at line 148 of file Analyzer_minbias.cc.

◆ tok_gtRec_

const edm::EDGetTokenT<L1GlobalTriggerReadoutRecord> cms::Analyzer_minbias::tok_gtRec_
private

Definition at line 155 of file Analyzer_minbias.cc.

◆ tok_hbhe_

const edm::EDGetTokenT<HBHERecHitCollection> cms::Analyzer_minbias::tok_hbhe_
private

Definition at line 144 of file Analyzer_minbias.cc.

Referenced by analyze().

◆ tok_hbheNoise_

const edm::EDGetTokenT<HBHERecHitCollection> cms::Analyzer_minbias::tok_hbheNoise_
private

Definition at line 150 of file Analyzer_minbias.cc.

Referenced by analyze().

◆ tok_hbheNorm_

const edm::EDGetTokenT<HBHERecHitCollection> cms::Analyzer_minbias::tok_hbheNorm_
private

Definition at line 156 of file Analyzer_minbias.cc.

Referenced by analyze().

◆ tok_hf_

const edm::EDGetTokenT<HFRecHitCollection> cms::Analyzer_minbias::tok_hf_
private

Definition at line 146 of file Analyzer_minbias.cc.

Referenced by analyze().

◆ tok_hfNoise_

const edm::EDGetTokenT<HFRecHitCollection> cms::Analyzer_minbias::tok_hfNoise_
private

Definition at line 152 of file Analyzer_minbias.cc.

Referenced by analyze().

◆ tok_ho_

const edm::EDGetTokenT<HORecHitCollection> cms::Analyzer_minbias::tok_ho_
private

Definition at line 145 of file Analyzer_minbias.cc.

◆ tok_hoNoise_

const edm::EDGetTokenT<HORecHitCollection> cms::Analyzer_minbias::tok_hoNoise_
private

Definition at line 151 of file Analyzer_minbias.cc.

◆ tok_l1gt_

const edm::ESGetToken<L1GtTriggerMenu, L1GtTriggerMenuRcd> cms::Analyzer_minbias::tok_l1gt_
private

Definition at line 159 of file Analyzer_minbias.cc.

◆ tok_respCorr_

const edm::ESGetToken<HcalRespCorrs, HcalRespCorrsRcd> cms::Analyzer_minbias::tok_respCorr_
private

Definition at line 158 of file Analyzer_minbias.cc.

Referenced by analyze().