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

#include <DTLocalTriggerTest.h>

Inheritance diagram for DTLocalTriggerTest:
DTLocalTriggerBaseTest DQMEDHarvester edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns, edm::one::SharedResources > edm::one::EDProducerBase edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

 DTLocalTriggerTest (const edm::ParameterSet &ps)
 Constructor. More...
 
 ~DTLocalTriggerTest () override
 Destructor. More...
 
- Public Member Functions inherited from DTLocalTriggerBaseTest
 DTLocalTriggerBaseTest ()
 Constructor. More...
 
 ~DTLocalTriggerBaseTest () override
 Destructor. More...
 
- Public Member Functions inherited from DQMEDHarvester
void accumulate (edm::Event const &ev, edm::EventSetup const &es) final
 
void beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final
 
void beginRun (edm::Run const &, edm::EventSetup const &) override
 
 DQMEDHarvester ()
 
void endJob () final
 
void endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final
 
void endLuminosityBlockProduce (edm::LuminosityBlock &, edm::EventSetup const &) final
 
void endRun (edm::Run const &, edm::EventSetup const &) override
 
void endRunProduce (edm::Run &run, edm::EventSetup const &setup) override
 
 ~DQMEDHarvester () override=default
 
- Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns, edm::one::SharedResources >
 EDProducer ()=default
 
SerialTaskQueue * globalLuminosityBlocksQueue () final
 
SerialTaskQueue * globalRunsQueue () final
 
bool hasAbilityToProduceInBeginLumis () const final
 
bool hasAbilityToProduceInBeginRuns () const final
 
bool hasAbilityToProduceInEndLumis () const final
 
bool hasAbilityToProduceInEndRuns () const final
 
bool wantsGlobalLuminosityBlocks () const final
 
bool wantsGlobalRuns () const final
 
- Public Member Functions inherited from edm::one::EDProducerBase
 EDProducerBase ()
 
ModuleDescription const & moduleDescription () const
 
bool wantsStreamLuminosityBlocks () const
 
bool wantsStreamRuns () const
 
 ~EDProducerBase () override
 
- Public Member Functions inherited from edm::ProducerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
std::vector< edm::ProductResolverIndex > const & indiciesForPutProducts (BranchType iBranchType) const
 
 ProducerBase ()
 
std::vector< edm::ProductResolverIndex > const & putTokenIndexToProductResolverIndex () const
 
void registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
 
std::function< void(BranchDescription const &)> registrationCallback () const
 used by the fwk to register list of products More...
 
void resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
 
 ~ProducerBase () noexcept(false) override
 
- 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
 
ESProxyIndex const * esGetTokenIndices (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::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
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Protected Member Functions

void beginRun (const edm::Run &r, const edm::EventSetup &c) override
 BeginRun. More...
 
void Bookings (DQMStore::IBooker &, DQMStore::IGetter &)
 
void fillGlobalSummary (DQMStore::IGetter &)
 
void runClientDiagnostic (DQMStore::IBooker &, DQMStore::IGetter &) override
 Run client analysis. More...
 
- Protected Member Functions inherited from DTLocalTriggerBaseTest
void beginRun (edm::Run const &run, edm::EventSetup const &context) override
 BeginRun. More...
 
void bookCmsHistos (DQMStore::IBooker &, std::string hTag, std::string folder="", bool isGlb=false)
 Book the new MEs (CMS summary) More...
 
void bookSectorHistos (DQMStore::IBooker &, int wheel, int sector, std::string hTag, std::string folder="")
 Book the new MEs (for each sector) More...
 
void bookWheelHistos (DQMStore::IBooker &, int wheel, std::string hTag, std::string folder="")
 Book the new MEs (for each wheel) More...
 
std::string category ()
 Get message logger name. More...
 
void dqmEndJob (DQMStore::IBooker &, DQMStore::IGetter &) override
 
void dqmEndLuminosityBlock (DQMStore::IBooker &, DQMStore::IGetter &, edm::LuminosityBlock const &, edm::EventSetup const &) override
 Perform client diagnostic in online. More...
 
void endRun (edm::Run const &run, edm::EventSetup const &context) override
 Perform client diagnostic in offline. More...
 
std::string fullName (std::string htype)
 Create fullname from histo partial name. More...
 
template<class T >
T * getHisto (MonitorElement *me)
 Convert ME to Histogram fo type T. More...
 
std::string getMEName (std::string histoTag, std::string subfolder, const DTChamberId &chambid)
 Get the ME name (by chamber) More...
 
std::string getMEName (std::string histoTag, std::string subfolder, int wh)
 Get the ME name (by wheel) More...
 
std::pair< float, float > phiRange (const DTChamberId &id)
 Calculate phi range for histograms. More...
 
void setConfig (const edm::ParameterSet &ps, std::string name)
 Set configuration variables. More...
 
std::string & topFolder ()
 Get top folder name. More...
 
- Protected Member Functions inherited from edm::ProducerBase
ProducesCollector producesCollector ()
 
- 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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
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)
 

Protected Attributes

const int wheelArrayShift = 3
 
- Protected Attributes inherited from DTLocalTriggerBaseTest
std::string baseFolderTM
 
std::map< std::string, MonitorElement * > cmsME
 
std::string hwSource
 
std::vector< std::string > hwSources
 
edm::ESHandle< DTGeometry > muonGeom
 
int nevents
 
unsigned int nLumiSegs
 
edm::ParameterSet parameters
 
int prescaleFactor
 
int run
 
bool runOnline
 
std::map< int, std::map< std::string, MonitorElement * > > secME
 
std::string sourceFolder
 
std::string testName
 
std::string trigSource
 
std::vector< std::string > trigSources
 
std::map< int, std::map< std::string, MonitorElement * > > whME
 
- Protected Attributes inherited from DQMEDHarvester
edm::EDPutTokenT< DQMToken > lumiToken_
 
edm::EDPutTokenT< DQMToken > runToken_
 

Private Attributes

bool bookingdone
 
int nMinEvts
 

Additional Inherited Members

- Public Types inherited from DQMEDHarvester
typedef dqm::harvesting::DQMStore DQMStore
 
typedef dqm::harvesting::MonitorElement MonitorElement
 
- Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase ModuleType
 
- Public Types inherited from edm::ProducerBase
using ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
 
typedef ProductRegistryHelper::TypeLabelList TypeLabelList
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 

Detailed Description

Definition at line 14 of file DTLocalTriggerTest.h.

Constructor & Destructor Documentation

DTLocalTriggerTest::DTLocalTriggerTest ( const edm::ParameterSet &  ps)

Constructor.

Definition at line 33 of file DTLocalTriggerTest.cc.

References edm::ParameterSet::getUntrackedParameter().

33  {
34  setConfig(ps, "DTLocalTrigger");
35  baseFolderTM = "DT/03-LocalTrigger-TM/";
36  nMinEvts = ps.getUntrackedParameter<int>("nEventsCert", 5000);
37 
38  bookingdone = false;
39 }
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
DTLocalTriggerBaseTest::setConfig
void setConfig(const edm::ParameterSet &ps, std::string name)
Set configuration variables.
Definition: DTLocalTriggerBaseTest.cc:78
DTLocalTriggerTest::~DTLocalTriggerTest ( )
override

Destructor.

Definition at line 41 of file DTLocalTriggerTest.cc.

41 {}

Member Function Documentation

void DTLocalTriggerTest::beginRun ( const edm::Run &  r,
const edm::EventSetup &  c 
)
overrideprotected

BeginRun.

Definition at line 105 of file DTLocalTriggerTest.cc.

References DTLocalTriggerBaseTest::beginRun().

105  {
106  DTLocalTriggerBaseTest::beginRun(r, c);
107 }
DTLocalTriggerBaseTest::beginRun
void beginRun(edm::Run const &run, edm::EventSetup const &context) override
BeginRun.
Definition: DTLocalTriggerBaseTest.cc:37
void DTLocalTriggerTest::Bookings ( DQMStore::IBooker &  ibooker,
DQMStore::IGetter &  igetter 
)
protected

Definition at line 43 of file DTLocalTriggerTest.cc.

References dtFineDelayCorr_cfi::hwSource, and dtFineDelayCorr_cfi::hwSources.

43  {
44  vector<string>::const_iterator iTr = trigSources.begin();
45  vector<string>::const_iterator trEnd = trigSources.end();
46  vector<string>::const_iterator iHw = hwSources.begin();
47  vector<string>::const_iterator hwEnd = hwSources.end();
48 
49  //Booking
50  if (parameters.getUntrackedParameter<bool>("staticBooking", true)) {
51  for (; iTr != trEnd; ++iTr) {
52  trigSource = (*iTr);
53  for (; iHw != hwEnd; ++iHw) {
54  hwSource = (*iHw);
55  // Loop over the TriggerUnits
56  for (int wh = -2; wh <= 2; ++wh) {
57  for (int sect = 1; sect <= 12; ++sect) {
58  bookSectorHistos(ibooker, wh, sect, "BXDistribPhiIn");
59  bookSectorHistos(ibooker, wh, sect, "QualDistribPhiIn");
60  bookSectorHistos(ibooker, wh, sect, "BXDistribPhiOut");
61  bookSectorHistos(ibooker, wh, sect, "QualDistribPhiOut");
62  }
63 
64  bookWheelHistos(ibooker, wh, "CorrectBXPhiIn");
65  bookWheelHistos(ibooker, wh, "ResidualBXPhiIn");
66  bookWheelHistos(ibooker, wh, "CorrFractionPhiIn");
67  bookWheelHistos(ibooker, wh, "2ndFractionPhiIn");
68  bookWheelHistos(ibooker, wh, "TriggerInclusivePhiIn");
69 
70  bookWheelHistos(ibooker, wh, "CorrectBXPhiOut");
71  bookWheelHistos(ibooker, wh, "ResidualBXPhiOut");
72  bookWheelHistos(ibooker, wh, "CorrFractionPhiOut");
73  bookWheelHistos(ibooker, wh, "2ndFractionPhiOut");
74  bookWheelHistos(ibooker, wh, "TriggerInclusivePhiOut");
75  }
76  }
77  }
78  }
79  // Summary test histo booking (only static)
80  for (iTr = trigSources.begin(); iTr != trEnd; ++iTr) {
81  trigSource = (*iTr);
82  for (iHw = hwSources.begin(); iHw != hwSources.end(); ++iHw) {
83  hwSource = (*iHw);
84  // Loop over the TriggerUnits
85  for (int wh = -2; wh <= 2; ++wh) {
86  bookWheelHistos(ibooker, wh, "CorrFractionSummaryIn", "Summaries");
87  bookWheelHistos(ibooker, wh, "2ndFractionSummaryIn", "Summaries");
88  bookWheelHistos(ibooker, wh, "CorrFractionSummaryOut", "Summaries");
89  bookWheelHistos(ibooker, wh, "2ndFractionSummaryOut", "Summaries");
90  }
91  bookCmsHistos(ibooker, "CorrFractionSummaryIn");
92  bookCmsHistos(ibooker, "2ndFractionSummaryIn");
93  bookCmsHistos(ibooker, "CorrFractionSummaryOut");
94  bookCmsHistos(ibooker, "2ndFractionSummaryOut");
95 
96  if (hwSource == "TM") {
97  bookCmsHistos(ibooker, "TrigGlbSummary", "", true);
98  }
99  }
100  }
101 
102  bookingdone = true;
103 }
DTLocalTriggerBaseTest::trigSources
std::vector< std::string > trigSources
Definition: DTLocalTriggerBaseTest.h:103
DTLocalTriggerBaseTest::bookSectorHistos
void bookSectorHistos(DQMStore::IBooker &, int wheel, int sector, std::string hTag, std::string folder="")
Book the new MEs (for each sector)
Definition: DTLocalTriggerBaseTest.cc:134
parameters
parameters
Definition: BeamSpot_PayloadInspector.cc:14
DTLocalTriggerBaseTest::hwSources
std::vector< std::string > hwSources
Definition: DTLocalTriggerBaseTest.h:104
DTLocalTriggerBaseTest::bookWheelHistos
void bookWheelHistos(DQMStore::IBooker &, int wheel, std::string hTag, std::string folder="")
Book the new MEs (for each wheel)
Definition: DTLocalTriggerBaseTest.cc:211
DTLocalTriggerBaseTest::bookCmsHistos
void bookCmsHistos(DQMStore::IBooker &, std::string hTag, std::string folder="", bool isGlb=false)
Book the new MEs (CMS summary)
Definition: DTLocalTriggerBaseTest.cc:194
void DTLocalTriggerTest::fillGlobalSummary ( DQMStore::IGetter &  igetter)
protected

Definition at line 430 of file DTLocalTriggerTest.cc.

References taus_updatedMVAIds_cff::category, corr, newFWLiteAna::fullName, dqm::dqmstoreimpl::DQMStore::IGetter::get(), dqm::impl::MonitorElement::getBinContent(), dqm::impl::MonitorElement::getFloatValue(), dtFineDelayCorr_cfi::hwSource, LogTrace, edm::min(), edm::second(), and L1TBPTX_cfi::testName.

430  {
431  float glbPerc[5] = {1., 0.9, 0.6, 0.3, 0.01};
432  trigSource = "";
433  hwSource = "TM";
434 
435  int nSecReadout = 0;
436 
437  for (int wh = -2; wh <= 2; ++wh) {
438  for (int sect = 1; sect <= 12; ++sect) {
439  float maxErr = 8.;
440  int corr = cmsME.find(fullName("CorrFractionSummaryIn"))->second->getBinContent(sect, wh + wheelArrayShift);
441  int second = cmsME.find(fullName("2ndFractionSummaryIn"))->second->getBinContent(sect, wh + wheelArrayShift);
442  int lut = 0;
443  MonitorElement* lutsME = igetter.get(topFolder() + "Summaries/TrigLutSummary");
444  if (lutsME) {
445  lut = lutsME->getBinContent(sect, wh + wheelArrayShift);
446  maxErr += 4;
447  } else {
448  LogTrace(category()) << "[" << testName << "Test]: TM Lut test Summary histo not found." << endl;
449  }
450  (corr < 5 || second < 5) && nSecReadout++;
451  int errcode = ((corr < 5 ? corr : 4) + (second < 5 ? second : 4) + (lut < 5 ? lut : 4));
452  errcode = min(int((errcode / maxErr + 0.01) * 5), 5);
453  cmsME.find("TrigGlbSummary")->second->setBinContent(sect, wh + wheelArrayShift, glbPerc[errcode]);
454  }
455  }
456 
457  if (!nSecReadout)
458  cmsME.find("TrigGlbSummary")->second->Reset(); // white histo id TM is not RO
459 
460  string nEvtsName = "DT/EventInfo/Counters/nProcessedEventsTrigger";
461  MonitorElement* meProcEvts = igetter.get(nEvtsName);
462 
463  if (meProcEvts) {
464  int nProcEvts = meProcEvts->getFloatValue();
465  cmsME.find("TrigGlbSummary")->second->setEntries(nProcEvts < nMinEvts ? 10. : nProcEvts);
466  } else {
467  cmsME.find("TrigGlbSummary")->second->setEntries(nMinEvts + 1);
468  LogVerbatim(category()) << "[" << testName << "Test]: ME: " << nEvtsName << " not found!" << endl;
469  }
470 }
DTLocalTriggerBaseTest::topFolder
std::string & topFolder()
Get top folder name.
Definition: DTLocalTriggerBaseTest.h:93
edm::second
U second(std::pair< T, U > const &p)
Definition: ParameterSet.cc:215
min
T min(T a, T b)
Definition: MathUtil.h:58
dqm::impl::MonitorElement::getBinContent
virtual double getBinContent(int binx) const
get content of bin (1-D)
Definition: MonitorElement.cc:668
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:671
corr
JetCorrectorParameters corr
Definition: classes.h:5
DTLocalTriggerBaseTest::category
std::string category()
Get message logger name.
Definition: DTLocalTriggerBaseTest.h:96
dqm::impl::MonitorElement::getFloatValue
virtual double getFloatValue() const
Definition: MonitorElement.cc:938
DTLocalTriggerBaseTest::cmsME
std::map< std::string, MonitorElement * > cmsME
Definition: DTLocalTriggerBaseTest.h:115
DTLocalTriggerBaseTest::fullName
std::string fullName(std::string htype)
Create fullname from histo partial name.
Definition: DTLocalTriggerBaseTest.cc:99
void DTLocalTriggerTest::runClientDiagnostic ( DQMStore::IBooker &  ibooker,
DQMStore::IGetter &  igetter 
)
overrideprotectedvirtual

Run client analysis.

Implements DTLocalTriggerBaseTest.

Definition at line 109 of file DTLocalTriggerTest.cc.

References L1TStage2uGTEmulatorClient_cff::BX, end, spr::find(), newFWLiteAna::fullName, dqm::dqmstoreimpl::DQMStore::IGetter::get(), dtFineDelayCorr_cfi::hwSource, dtFineDelayCorr_cfi::hwSources, and hgcalPlots::stat.

109  {
110  if (!bookingdone)
111  Bookings(ibooker, igetter);
112 
113  // Loop over Trig & Hw sources
114  for (vector<string>::const_iterator iTr = trigSources.begin(); iTr != trigSources.end(); ++iTr) {
115  trigSource = (*iTr);
116 
117  for (vector<string>::const_iterator iHw = hwSources.begin(); iHw != hwSources.end(); ++iHw) {
118  hwSource = (*iHw);
119  // Loop over the TriggerUnits
120  for (int stat = 1; stat <= 4; ++stat) {
121  for (int wh = -2; wh <= 2; ++wh) {
122  for (int sect = 1; sect <= 12; ++sect) {
123  DTChamberId chId(wh, stat, sect);
124  int sector_id = (wh + wheelArrayShift) + (sect - 1) * 5;
125 
126  // IN part
127  TH2F* BXvsQual = getHisto<TH2F>(igetter.get(getMEName("BXvsQual_In", "LocalTriggerPhiIn", chId)));
128  TH1F* BestQual = getHisto<TH1F>(igetter.get(getMEName("BestQual_In", "LocalTriggerPhiIn", chId)));
129  TH2F* Flag1stvsQual = getHisto<TH2F>(igetter.get(getMEName("Flag1stvsQual_In", "LocalTriggerPhiIn", chId)));
130  if (BXvsQual && Flag1stvsQual && BestQual) {
131  int corrSummary = 1;
132  int secondSummary = 1;
133  //default values for histograms
134  double BX_OK = 51.;
135  double BXMean = 51.;
136  double corrFrac = 0.;
137  double secondFrac = 0.;
138  double besttrigs = 0.;
139  if (BestQual->GetEntries() > 1) {
140  TH1D* BXHH = BXvsQual->ProjectionY("", 6, 7, "");
141  TH1D* Flag1st = Flag1stvsQual->ProjectionY();
142  int BXOK_bin = BXHH->GetEntries() >= 1 ? BXHH->GetMaximumBin() : 51;
143  BXMean = BXHH->GetEntries() >= 1 ? BXHH->GetMean() : 51;
144  BX_OK = BXvsQual->GetYaxis()->GetBinCenter(BXOK_bin);
145  double trigsFlag2nd = Flag1st->GetBinContent(2);
146  double trigs = Flag1st->GetEntries();
147  besttrigs = BestQual->GetEntries();
148  double besttrigsCorr = BestQual->Integral(5, 7, "");
149  delete BXHH;
150  delete Flag1st;
151 
152  corrFrac = besttrigsCorr / besttrigs;
153  secondFrac = trigsFlag2nd / trigs;
154  if (corrFrac < parameters.getUntrackedParameter<double>("corrFracError", .5)) {
155  corrSummary = 2;
156  } else if (corrFrac < parameters.getUntrackedParameter<double>("corrFracWarning", .6)) {
157  corrSummary = 3;
158  } else {
159  corrSummary = 0;
160  }
161  if (secondFrac > parameters.getUntrackedParameter<double>("secondFracError", .2)) {
162  secondSummary = 2;
163  } else if (secondFrac > parameters.getUntrackedParameter<double>("secondFracWarning", .1)) {
164  secondSummary = 3;
165  } else {
166  secondSummary = 0;
167  }
168 
169  if (secME[sector_id].find(fullName("BXDistribPhiIn")) == secME[sector_id].end()) {
170  bookSectorHistos(ibooker, wh, sect, "QualDistribPhiIn");
171  bookSectorHistos(ibooker, wh, sect, "BXDistribPhiIn");
172  }
173 
174  TH1D* BXDistr = BXvsQual->ProjectionY();
175  TH1D* QualDistr = BXvsQual->ProjectionX();
176  std::map<std::string, MonitorElement*>* innerME = &(secME[sector_id]);
177 
178  int nbinsBX = BXDistr->GetNbinsX();
179  int firstBinCenter = static_cast<int>(BXDistr->GetBinCenter(1));
180  int lastBinCenter = static_cast<int>(BXDistr->GetBinCenter(nbinsBX));
181  int iMin = firstBinCenter > -4 ? firstBinCenter : -4;
182  int iMax = lastBinCenter < 20 ? lastBinCenter : 20;
183  for (int ibin = iMin + 5; ibin <= iMax + 5; ++ibin) {
184  innerME->find(fullName("BXDistribPhiIn"))
185  ->second->setBinContent(ibin, stat, BXDistr->GetBinContent(ibin - 5 - firstBinCenter + 1));
186  }
187  for (int ibin = 1; ibin <= 7; ++ibin) {
188  innerME->find(fullName("QualDistribPhiIn"))
189  ->second->setBinContent(ibin, stat, QualDistr->GetBinContent(ibin));
190  }
191 
192  delete BXDistr;
193  delete QualDistr;
194  }
195 
196  std::map<std::string, MonitorElement*>* innerME = &(secME[sector_id]);
197 
198  if (whME[wh].find(fullName("CorrectBXPhiIn")) == whME[wh].end()) {
199  bookWheelHistos(ibooker, wh, "ResidualBXPhiIn");
200  bookWheelHistos(ibooker, wh, "CorrectBXPhiIn");
201  bookWheelHistos(ibooker, wh, "CorrFractionPhiIn");
202  bookWheelHistos(ibooker, wh, "2ndFractionPhiIn");
203  bookWheelHistos(ibooker, wh, "TriggerInclusivePhiIn");
204  }
205 
206  innerME = &(whME[wh]);
207  innerME->find(fullName("CorrectBXPhiIn"))->second->setBinContent(sect, stat, BX_OK + 0.00001);
208  innerME->find(fullName("ResidualBXPhiIn"))
209  ->second->setBinContent(sect, stat, round(25. * (BXMean - BX_OK)) + 0.00001);
210  innerME->find(fullName("CorrFractionPhiIn"))->second->setBinContent(sect, stat, corrFrac);
211  innerME->find(fullName("TriggerInclusivePhiIn"))->second->setBinContent(sect, stat, besttrigs);
212  innerME->find(fullName("2ndFractionPhiIn"))->second->setBinContent(sect, stat, secondFrac);
213 
214  whME[wh].find(fullName("CorrFractionSummaryIn"))->second->setBinContent(sect, stat, corrSummary);
215  whME[wh].find(fullName("2ndFractionSummaryIn"))->second->setBinContent(sect, stat, secondSummary);
216 
217  } // closes BXvsQual && Flag1stvsQual && BestQual
218 
219  if (hwSource == "TM") {
220  //Out part
221 
222  TH2F* BXvsQual = getHisto<TH2F>(igetter.get(getMEName("BXvsQual_Out", "LocalTriggerPhiOut", chId)));
223  TH1F* BestQual = getHisto<TH1F>(igetter.get(getMEName("BestQual_Out", "LocalTriggerPhiOut", chId)));
224  TH2F* Flag1stvsQual =
225  getHisto<TH2F>(igetter.get(getMEName("Flag1stvsQual_Out", "LocalTriggerPhiOut", chId)));
226  if (BXvsQual && Flag1stvsQual && BestQual) {
227  int corrSummary = 1;
228  int secondSummary = 1;
229  //default values for histograms
230  double BX_OK = 51.;
231  double BXMean = 51.;
232  double corrFrac = 0.;
233  double secondFrac = 0.;
234  double besttrigs = 0.;
235 
236  if (BestQual->GetEntries() > 1) {
237  TH1D* BXHH = BXvsQual->ProjectionY("", 6, 7, "");
238  TH1D* Flag1st = Flag1stvsQual->ProjectionY();
239  int BXOK_bin = BXHH->GetEntries() >= 1 ? BXHH->GetMaximumBin() : 51;
240  BXMean = BXHH->GetEntries() >= 1 ? BXHH->GetMean() : 51;
241  BX_OK = BXvsQual->GetYaxis()->GetBinCenter(BXOK_bin);
242  double trigsFlag2nd = Flag1st->GetBinContent(2);
243  double trigs = Flag1st->GetEntries();
244  besttrigs = BestQual->GetEntries();
245  double besttrigsCorr = BestQual->Integral(5, 7, "");
246  delete BXHH;
247  delete Flag1st;
248 
249  corrFrac = besttrigsCorr / besttrigs;
250  secondFrac = trigsFlag2nd / trigs;
251  if (corrFrac < parameters.getUntrackedParameter<double>("corrFracError", .5)) {
252  corrSummary = 2;
253  } else if (corrFrac < parameters.getUntrackedParameter<double>("corrFracWarning", .6)) {
254  corrSummary = 3;
255  } else {
256  corrSummary = 0;
257  }
258  if (secondFrac > parameters.getUntrackedParameter<double>("secondFracError", .2)) {
259  secondSummary = 2;
260  } else if (secondFrac > parameters.getUntrackedParameter<double>("secondFracWarning", .1)) {
261  secondSummary = 3;
262  } else {
263  secondSummary = 0;
264  }
265 
266  if (secME[sector_id].find(fullName("BXDistribPhiOut")) == secME[sector_id].end()) {
267  bookSectorHistos(ibooker, wh, sect, "QualDistribPhiOut");
268  bookSectorHistos(ibooker, wh, sect, "BXDistribPhiOut");
269  }
270 
271  TH1D* BXDistr = BXvsQual->ProjectionY();
272  TH1D* QualDistr = BXvsQual->ProjectionX();
273  std::map<std::string, MonitorElement*>* innerME = &(secME[sector_id]);
274 
275  int nbinsBX = BXDistr->GetNbinsX();
276  int firstBinCenter = static_cast<int>(BXDistr->GetBinCenter(1));
277  int lastBinCenter = static_cast<int>(BXDistr->GetBinCenter(nbinsBX));
278  int iMin = firstBinCenter > -4 ? firstBinCenter : -4;
279  int iMax = lastBinCenter < 20 ? lastBinCenter : 20;
280  for (int ibin = iMin + 5; ibin <= iMax + 5; ++ibin) {
281  innerME->find(fullName("BXDistribPhiOut"))
282  ->second->setBinContent(ibin, stat, BXDistr->GetBinContent(ibin - 5 - firstBinCenter + 1));
283  }
284  for (int ibin = 1; ibin <= 7; ++ibin) {
285  innerME->find(fullName("QualDistribPhiOut"))
286  ->second->setBinContent(ibin, stat, QualDistr->GetBinContent(ibin));
287  }
288 
289  delete BXDistr;
290  delete QualDistr;
291  }
292 
293  std::map<std::string, MonitorElement*>* innerME = &(secME[sector_id]);
294 
295  if (whME[wh].find(fullName("CorrectBXPhiOut")) == whME[wh].end()) {
296  bookWheelHistos(ibooker, wh, "ResidualBXPhiOut");
297  bookWheelHistos(ibooker, wh, "CorrectBXPhiOut");
298  bookWheelHistos(ibooker, wh, "CorrFractionPhiOut");
299  bookWheelHistos(ibooker, wh, "2ndFractionPhiOut");
300  bookWheelHistos(ibooker, wh, "TriggerInclusivePhiOut");
301  }
302 
303  innerME = &(whME[wh]);
304  innerME->find(fullName("CorrectBXPhiOut"))->second->setBinContent(sect, stat, BX_OK + 0.00001);
305  innerME->find(fullName("ResidualBXPhiOut"))
306  ->second->setBinContent(sect, stat, round(25. * (BXMean - BX_OK)) + 0.00001);
307  innerME->find(fullName("CorrFractionPhiOut"))->second->setBinContent(sect, stat, corrFrac);
308  innerME->find(fullName("TriggerInclusivePhiOut"))->second->setBinContent(sect, stat, besttrigs);
309  innerME->find(fullName("2ndFractionPhiOut"))->second->setBinContent(sect, stat, secondFrac);
310 
311  whME[wh].find(fullName("CorrFractionSummaryOut"))->second->setBinContent(sect, stat, corrSummary);
312  whME[wh].find(fullName("2ndFractionSummaryOut"))->second->setBinContent(sect, stat, secondSummary);
313 
314  } // closes BXvsQual && Flag1stvsQual && BestQual
315 
316  } // Check on TM source
317  //Theta part
318  if (hwSource == "TM") {
319  // Perform TM plot analysis (Theta ones)
320  TH2F* ThetaPosvsBX = getHisto<TH2F>(igetter.get(getMEName("PositionvsBX", "LocalTriggerTheta", chId)));
321  double BX_OK = 48;
322  // no theta triggers in stat 4!
323  if (ThetaPosvsBX && stat < 4 && ThetaPosvsBX->GetEntries() > 1) {
324  TH1D* BX = ThetaPosvsBX->ProjectionX();
325  int BXOK_bin = BX->GetEffectiveEntries() >= 1 ? BX->GetMaximumBin() : 10;
326  BX_OK = ThetaPosvsBX->GetXaxis()->GetBinCenter(BXOK_bin);
327  delete BX;
328 
329  if (whME[wh].find(fullName("CorrectBXTheta")) == whME[wh].end()) {
330  bookWheelHistos(ibooker, wh, "CorrectBXTheta");
331  }
332  std::map<std::string, MonitorElement*>* innerME = &(whME.find(wh)->second);
333  innerME->find(fullName("CorrectBXTheta"))->second->setBinContent(sect, stat, BX_OK + 0.00001);
334  }
335  // Adding trigger info to compute H fraction (11/10/2016) M.C.Fouz
336  TH2F* ThetaBXvsQual = getHisto<TH2F>(igetter.get(getMEName("ThetaBXvsQual", "LocalTriggerTheta", chId)));
337  TH1F* ThetaBestQual = getHisto<TH1F>(igetter.get(getMEName("ThetaBestQual", "LocalTriggerTheta", chId)));
338  if (ThetaBXvsQual && ThetaBestQual && stat < 4 && ThetaBestQual->GetEntries() > 1) {
339  double trigs = ThetaBestQual->GetEntries();
340  double trigsH = ThetaBestQual->GetBinContent(
341  2); // Note that for the new plots H is at bin=2 and not 4 as in DDU!!!!
342  if (whME[wh].find(fullName("HFractionTheta")) == whME[wh].end()) {
343  bookWheelHistos(ibooker, wh, "HFractionTheta");
344  }
345  std::map<std::string, MonitorElement*>* innerME = &(whME.find(wh)->second);
346  innerME->find(fullName("HFractionTheta"))->second->setBinContent(sect, stat, trigsH / trigs);
347  }
348  // END ADDING H Fraction info
349  }
350  }
351  }
352  }
353  }
354  }
355 
356  for (vector<string>::const_iterator iTr = trigSources.begin(); iTr != trigSources.end(); ++iTr) {
357  trigSource = (*iTr);
358  for (vector<string>::const_iterator iHw = hwSources.begin(); iHw != hwSources.end(); ++iHw) {
359  hwSource = (*iHw);
360  for (int wh = -2; wh <= 2; ++wh) {
361  std::map<std::string, MonitorElement*>* innerME = &(whME[wh]);
362  // In part
363  TH2F* corrWhSummaryIn = getHisto<TH2F>(innerME->find(fullName("CorrFractionSummaryIn"))->second);
364  TH2F* secondWhSummaryIn = getHisto<TH2F>(innerME->find(fullName("2ndFractionSummaryIn"))->second);
365  for (int sect = 1; sect <= 12; ++sect) {
366  int corrErr = 0;
367  int secondErr = 0;
368  int corrNoData = 0;
369  int secondNoData = 0;
370  for (int stat = 1; stat <= 4; ++stat) {
371  switch (static_cast<int>(corrWhSummaryIn->GetBinContent(sect, stat))) {
372  case 1:
373  corrNoData++;
374  [[fallthrough]];
375  case 2:
376  corrErr++;
377  }
378  switch (static_cast<int>(secondWhSummaryIn->GetBinContent(sect, stat))) {
379  case 1:
380  secondNoData++;
381  [[fallthrough]];
382  case 2:
383  secondErr++;
384  }
385  }
386  if (corrNoData == 4)
387  corrErr = 5;
388  if (secondNoData == 4)
389  secondErr = 5;
390  cmsME.find(fullName("CorrFractionSummaryIn"))->second->setBinContent(sect, wh + wheelArrayShift, corrErr);
391  cmsME.find(fullName("2ndFractionSummaryIn"))->second->setBinContent(sect, wh + wheelArrayShift, secondErr);
392  }
393  // Out part
394  TH2F* corrWhSummaryOut = getHisto<TH2F>(innerME->find(fullName("CorrFractionSummaryOut"))->second);
395  TH2F* secondWhSummaryOut = getHisto<TH2F>(innerME->find(fullName("2ndFractionSummaryOut"))->second);
396  for (int sect = 1; sect <= 12; ++sect) {
397  int corrErr = 0;
398  int secondErr = 0;
399  int corrNoData = 0;
400  int secondNoData = 0;
401  for (int stat = 1; stat <= 4; ++stat) {
402  switch (static_cast<int>(corrWhSummaryOut->GetBinContent(sect, stat))) {
403  case 1:
404  corrNoData++;
405  [[fallthrough]];
406  case 2:
407  corrErr++;
408  }
409  switch (static_cast<int>(secondWhSummaryOut->GetBinContent(sect, stat))) {
410  case 1:
411  secondNoData++;
412  [[fallthrough]];
413  case 2:
414  secondErr++;
415  }
416  }
417  if (corrNoData == 4)
418  corrErr = 5;
419  if (secondNoData == 4)
420  secondErr = 5;
421  cmsME.find(fullName("CorrFractionSummaryOut"))->second->setBinContent(sect, wh + wheelArrayShift, corrErr);
422  cmsME.find(fullName("2ndFractionSummaryOut"))->second->setBinContent(sect, wh + wheelArrayShift, secondErr);
423  }
424  }
425  }
426  }
427  fillGlobalSummary(igetter);
428 }
DTLocalTriggerBaseTest::trigSources
std::vector< std::string > trigSources
Definition: DTLocalTriggerBaseTest.h:103
L1TStage2uGTEmulatorClient_cff.BX
BX
Definition: L1TStage2uGTEmulatorClient_cff.py:9
DTLocalTriggerTest::Bookings
void Bookings(DQMStore::IBooker &, DQMStore::IGetter &)
Definition: DTLocalTriggerTest.cc:43
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:19
DTLocalTriggerTest::fillGlobalSummary
void fillGlobalSummary(DQMStore::IGetter &)
Definition: DTLocalTriggerTest.cc:430
DTLocalTriggerBaseTest::whME
std::map< int, std::map< std::string, MonitorElement * > > whME
Definition: DTLocalTriggerBaseTest.h:114
end
#define end
Definition: vmac.h:39
DTLocalTriggerBaseTest::getMEName
std::string getMEName(std::string histoTag, std::string subfolder, const DTChamberId &chambid)
Get the ME name (by chamber)
DTLocalTriggerBaseTest::bookSectorHistos
void bookSectorHistos(DQMStore::IBooker &, int wheel, int sector, std::string hTag, std::string folder="")
Book the new MEs (for each sector)
Definition: DTLocalTriggerBaseTest.cc:134
DTLocalTriggerBaseTest::cmsME
std::map< std::string, MonitorElement * > cmsME
Definition: DTLocalTriggerBaseTest.h:115
parameters
parameters
Definition: BeamSpot_PayloadInspector.cc:14
DTLocalTriggerBaseTest::hwSources
std::vector< std::string > hwSources
Definition: DTLocalTriggerBaseTest.h:104
DTLocalTriggerBaseTest::secME
std::map< int, std::map< std::string, MonitorElement * > > secME
Definition: DTLocalTriggerBaseTest.h:113
DTLocalTriggerBaseTest::bookWheelHistos
void bookWheelHistos(DQMStore::IBooker &, int wheel, std::string hTag, std::string folder="")
Book the new MEs (for each wheel)
Definition: DTLocalTriggerBaseTest.cc:211
DTLocalTriggerBaseTest::fullName
std::string fullName(std::string htype)
Create fullname from histo partial name.
Definition: DTLocalTriggerBaseTest.cc:99

Member Data Documentation

bool DTLocalTriggerTest::bookingdone
private

Definition at line 39 of file DTLocalTriggerTest.h.

int DTLocalTriggerTest::nMinEvts
private

Definition at line 37 of file DTLocalTriggerTest.h.

const int DTLocalTriggerTest::wheelArrayShift = 3
protected

Definition at line 34 of file DTLocalTriggerTest.h.

Generated for CMSSW Reference Manual by doxygen 1.8.11