CMS 3D CMS Logo

List of all members | Public Member Functions | Protected Member Functions | Protected Attributes | Private Member Functions | Private Attributes
DTLocalTriggerLutTest Class Reference

#include <DTLocalTriggerLutTest.h>

Inheritance diagram for DTLocalTriggerLutTest:
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

 DTLocalTriggerLutTest (const edm::ParameterSet &ps)
 Constructor. More...
 
 ~DTLocalTriggerLutTest () 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
 
SerialTaskQueueglobalLuminosityBlocksQueue () final
 
SerialTaskQueueglobalRunsQueue () final
 
bool hasAbilityToProduceInLumis () const final
 
bool hasAbilityToProduceInRuns () 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< 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
 
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
 
EDConsumerBaseoperator= (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 runClientDiagnostic (DQMStore::IBooker &ibooker, DQMStore::IGetter &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 >
TgetHisto (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 (bool isTM)
 Get top folder name. More...
 
- 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 baseFolderDDU
 
std::string baseFolderTM
 
std::map< std::string, MonitorElement * > cmsME
 
std::string hwSource
 
std::vector< std::string > hwSources
 
edm::ESHandle< DTGeometrymuonGeom
 
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< DQMTokenlumiToken_
 
edm::EDPutTokenT< DQMTokenrunToken_
 

Private Member Functions

void fillWhPlot (MonitorElement *plot, int sect, int stat, float value, bool lessIsBest=true)
 Fill summary plots managing double MB4 chambers. More...
 
int performLutTest (double mean, double RMS, double thresholdMean, double thresholdRMS)
 Perform Lut Test logical operations. More...
 

Private Attributes

bool bookingdone
 
bool doCorrStudy
 
double thresholdPhibMean
 
double thresholdPhibRMS
 
double thresholdPhiMean
 
double thresholdPhiRMS
 

Additional Inherited Members

- 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 21 of file DTLocalTriggerLutTest.h.

Constructor & Destructor Documentation

DTLocalTriggerLutTest::DTLocalTriggerLutTest ( const edm::ParameterSet ps)

Constructor.

Definition at line 37 of file DTLocalTriggerLutTest.cc.

References edm::ParameterSet::getUntrackedParameter().

37  {
38 
39  setConfig(ps,"DTLocalTriggerLut");
40  baseFolderTM = "DT/03-LocalTrigger-TM/";
41  baseFolderDDU = "DT/04-LocalTrigger-DDU/";
42  thresholdPhiMean = ps.getUntrackedParameter<double>("thresholdPhiMean",1.5);
43  thresholdPhiRMS = ps.getUntrackedParameter<double>("thresholdPhiRMS",.5);
44  thresholdPhibMean = ps.getUntrackedParameter<double>("thresholdPhibMean",1.5);
45  thresholdPhibRMS = ps.getUntrackedParameter<double>("thresholdPhibRMS",.8);
46  doCorrStudy = ps.getUntrackedParameter<bool>("doCorrelationStudy",false);
47 
48  bookingdone = false;
49 
50 }
T getUntrackedParameter(std::string const &, T const &) const
void setConfig(const edm::ParameterSet &ps, std::string name)
Set configuration variables.
DTLocalTriggerLutTest::~DTLocalTriggerLutTest ( )
override

Destructor.

Definition at line 53 of file DTLocalTriggerLutTest.cc.

53  {
54 
55 }

Member Function Documentation

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

BeginRun.

Definition at line 112 of file DTLocalTriggerLutTest.cc.

References DTLocalTriggerBaseTest::beginRun().

112  {
113 
115 }
void beginRun(edm::Run const &run, edm::EventSetup const &context) override
BeginRun.
void DTLocalTriggerLutTest::Bookings ( DQMStore::IBooker ibooker,
DQMStore::IGetter igetter 
)
protected

Definition at line 58 of file DTLocalTriggerLutTest.cc.

58  {
59 
60  vector<string>::const_iterator iTr = trigSources.begin();
61  vector<string>::const_iterator trEnd = trigSources.end();
62  vector<string>::const_iterator iHw = hwSources.begin();
63  vector<string>::const_iterator hwEnd = hwSources.end();
64 
65  //Booking
66  if(parameters.getUntrackedParameter<bool>("staticBooking", true)){
67  for (; iTr != trEnd; ++iTr){
68  trigSource = (*iTr);
69  for (; iHw != hwEnd; ++iHw){
70  hwSource = (*iHw);
71  // Loop over the TriggerUnits
72  for (int wh=-2; wh<=2; ++wh){
73  bookWheelHistos(ibooker,wh,"PhiResidualMean");
74  bookWheelHistos(ibooker,wh,"PhiResidualRMS");
75  bookWheelHistos(ibooker,wh,"PhibResidualMean");
76  bookWheelHistos(ibooker,wh,"PhibResidualRMS");
77  if (doCorrStudy) {
78 
79  bookWheelHistos(ibooker,wh,"PhiTkvsTrigSlope");
80  bookWheelHistos(ibooker,wh,"PhiTkvsTrigIntercept");
81  bookWheelHistos(ibooker,wh,"PhiTkvsTrigCorr");
82  bookWheelHistos(ibooker,wh,"PhibTkvsTrigSlope");
83  bookWheelHistos(ibooker,wh,"PhibTkvsTrigIntercept");
84  bookWheelHistos(ibooker,wh,"PhibTkvsTrigCorr");
85  }
86  }
87  }
88  }
89  }
90 
91  // Summary test histo booking (only static)
92  for (iTr = trigSources.begin(); iTr != trEnd; ++iTr){
93  trigSource = (*iTr);
94  for (iHw = hwSources.begin(); iHw != hwSources.end(); ++iHw){
95  hwSource = (*iHw);
96  // Loop over the TriggerUnits
97  for (int wh=-2; wh<=2; ++wh){
98 
99  bookWheelHistos(ibooker,wh,"PhiLutSummary");
100  bookWheelHistos(ibooker,wh,"PhibLutSummary");
101  }
102 
103  bookCmsHistos(ibooker,"PhiLutSummary");
104  bookCmsHistos(ibooker,"PhibLutSummary");
105  }
106  }
107 
108  bookingdone = true;
109 
110 }
std::vector< std::string > trigSources
std::vector< std::string > hwSources
void bookWheelHistos(DQMStore::IBooker &, int wheel, std::string hTag, std::string folder="")
Book the new MEs (for each wheel)
void bookCmsHistos(DQMStore::IBooker &, std::string hTag, std::string folder="", bool isGlb=false)
Book the new MEs (CMS summary)
void DTLocalTriggerLutTest::fillWhPlot ( MonitorElement plot,
int  sect,
int  stat,
float  value,
bool  lessIsBest = true 
)
private

Fill summary plots managing double MB4 chambers.

Definition at line 332 of file DTLocalTriggerLutTest.cc.

References MonitorElement::getBinContent(), and MonitorElement::setBinContent().

332  {
333 
334  if (sect>12) {
335  int scsect = sect==13 ? 4 : 10;
336  if ( (fabs(value)>fabs(plot->getBinContent(scsect,stat)))==lessIsBest) {
337  plot->setBinContent(scsect,stat,value);
338  }
339  }
340  else {
341  plot->setBinContent(sect,stat,value);
342  }
343 
344  return;
345 
346 }
void setBinContent(int binx, double content)
set content of bin (1-D)
Definition: value.py:1
double getBinContent(int binx) const
get content of bin (1-D)
int DTLocalTriggerLutTest::performLutTest ( double  mean,
double  RMS,
double  thresholdMean,
double  thresholdRMS 
)
private

Perform Lut Test logical operations.

Definition at line 323 of file DTLocalTriggerLutTest.cc.

323  {
324 
325  bool meanErr = fabs(mean)>thresholdMean;
326  bool rmsErr = RMS>thresholdRMS;
327 
328  return (meanErr || rmsErr) ? 2+(meanErr!=rmsErr) : 0 ;
329 
330 }
void DTLocalTriggerLutTest::runClientDiagnostic ( DQMStore::IBooker ibooker,
DQMStore::IGetter igetter 
)
overrideprotectedvirtual

Run client analysis.

Implements DTLocalTriggerBaseTest.

Definition at line 117 of file DTLocalTriggerLutTest.cc.

References taus_updatedMVAIds_cff::category, spr::find(), DQMStore::IGetter::get(), edm::second(), trackingPlots::stat, and DTChamberId::wheel().

117  {
118 
119  if (!bookingdone) Bookings(ibooker,igetter);
120 
121  // Loop over Trig & Hw sources
122  for (vector<string>::const_iterator iTr = trigSources.begin(); iTr != trigSources.end(); ++iTr){
123  trigSource = (*iTr);
124  for (vector<string>::const_iterator iHw = hwSources.begin(); iHw != hwSources.end(); ++iHw){
125  hwSource = (*iHw);
126  vector<const DTChamber*>::const_iterator chIt = muonGeom->chambers().begin();
127  vector<const DTChamber*>::const_iterator chEnd = muonGeom->chambers().end();
128  for (; chIt != chEnd; ++chIt) {
129  DTChamberId chId((*chIt)->id());
130  int wh = chId.wheel();
131  int sect = chId.sector();
132  int stat = chId.station();
133 
134 
135  if (doCorrStudy) {
136  // Perform Correlation Plots analysis (TM + segment Phi)
137 
138  TH2F * TrackPhitkvsPhitrig = getHisto<TH2F>(igetter.get(getMEName("PhitkvsPhitrig","Segment", chId)));
139 
140  if (TrackPhitkvsPhitrig && TrackPhitkvsPhitrig->GetEntries()>10) {
141 
142  // Fill client histos
143  if( whME[wh].find(fullName("PhiTkvsTrigCorr")) == whME[wh].end() ){
144 
145  bookWheelHistos(ibooker,wh,"PhiTkvsTrigSlope");
146  bookWheelHistos(ibooker,wh,"PhiTkvsTrigIntercept");
147  bookWheelHistos(ibooker,wh,"PhiTkvsTrigCorr");
148  }
149 
150  TProfile* PhitkvsPhitrigProf = TrackPhitkvsPhitrig->ProfileX();
151  double phiInt = 0;
152  double phiSlope = 0;
153  double phiCorr = 0;
154  try {
155  TF1 ffPhi("mypol1","pol1");
156  PhitkvsPhitrigProf->Fit(&ffPhi,"CQO");
157  phiInt = ffPhi.GetParameter(0);
158  phiSlope = ffPhi.GetParameter(1);
159  phiCorr = TrackPhitkvsPhitrig->GetCorrelationFactor();
160  } catch (cms::Exception& iException) {
161  edm::LogError(category()) << "[" << testName << "Test]: Error fitting PhitkvsPhitrig for Wheel " << wh
162  <<" Sector " << sect << " Station " << stat;
163  }
164 
165  std::map<std::string,MonitorElement*> &innerME = whME[wh];
166  fillWhPlot(innerME.find(fullName("PhiTkvsTrigSlope"))->second,sect,stat,phiSlope-1);
167  fillWhPlot(innerME.find(fullName("PhiTkvsTrigIntercept"))->second,sect,stat,phiInt);
168  fillWhPlot(innerME.find(fullName("PhiTkvsTrigCorr"))->second,sect,stat,phiCorr,false);
169 
170  }
171 
172  // Perform Correlation Plots analysis (TM + segment Phib)
173  TH2F * TrackPhibtkvsPhibtrig = getHisto<TH2F>(igetter.get(getMEName("PhibtkvsPhibtrig","Segment", chId)));
174 
175  if (stat != 3 && TrackPhibtkvsPhibtrig && TrackPhibtkvsPhibtrig->GetEntries()>10) {// station 3 has no meaningful MB3 phi bending information
176 
177  // Fill client histos
178  if( whME[wh].find(fullName("PhibTkvsTrigCorr")) == whME[wh].end() ){
179  bookWheelHistos(ibooker,wh,"PhibTkvsTrigSlope");
180  bookWheelHistos(ibooker,wh,"PhibTkvsTrigIntercept");
181  bookWheelHistos(ibooker,wh,"PhibTkvsTrigCorr");
182  }
183 
184  TProfile* PhibtkvsPhibtrigProf = TrackPhibtkvsPhibtrig->ProfileX();
185  double phibInt = 0;
186  double phibSlope = 0;
187  double phibCorr = 0;
188  try {
189  TF1 ffPhib("ffPhib","pol1");
190  PhibtkvsPhibtrigProf->Fit(&ffPhib,"CQO");
191  phibInt = ffPhib.GetParameter(0);
192  phibSlope = ffPhib.GetParameter(1);
193  phibCorr = TrackPhibtkvsPhibtrig->GetCorrelationFactor();
194  } catch (cms::Exception& iException) {
195  edm::LogError(category()) << "[" << testName << "Test]: Error fitting PhibtkvsPhibtrig for Wheel " << wh
196  <<" Sector " << sect << " Station " << stat;
197  }
198 
199  std::map<std::string,MonitorElement*> &innerME = whME[wh];
200  fillWhPlot(innerME.find(fullName("PhibTkvsTrigSlope"))->second,sect,stat,phibSlope-1);
201  fillWhPlot(innerME.find(fullName("PhibTkvsTrigIntercept"))->second,sect,stat,phibInt);
202  fillWhPlot(innerME.find(fullName("PhibTkvsTrigCorr"))->second,sect,stat,phibCorr,false);
203 
204  }
205 
206  }
207 
208  // Make Phi Residual Summary
209 
210  TH1F * PhiResidual = getHisto<TH1F>(igetter.get(getMEName("PhiResidual","Segment", chId)));
211  int phiSummary = 1;
212 
213  if (PhiResidual && PhiResidual->GetEffectiveEntries()>10) {
214 
215  // Fill client histos
216  if( whME[wh].find(fullName("PhiResidualMean")) == whME[wh].end() ){
217  bookWheelHistos(ibooker,wh,"PhiResidualMean");
218  bookWheelHistos(ibooker,wh,"PhiResidualRMS");
219  }
220 
221  double peak = PhiResidual->GetBinCenter(PhiResidual->GetMaximumBin());
222  double phiMean = 0;
223  double phiRMS = 0;
224  try {
225  TF1 ffPhi("ffPhi","gaus");
226  PhiResidual->Fit(&ffPhi,"CQO","",peak-5,peak+5);
227  phiMean = ffPhi.GetParameter(1);
228  phiRMS = ffPhi.GetParameter(2);
229  } catch (cms::Exception& iException) {
230  edm::LogError(category()) << "[" << testName << "Test]: Error fitting PhiResidual for Wheel " << wh
231  <<" Sector " << sect << " Station " << stat;
232  }
233 
234  std::map<std::string,MonitorElement*> &innerME = whME[wh];
235  fillWhPlot(innerME.find(fullName("PhiResidualMean"))->second,sect,stat,phiMean);
236  fillWhPlot(innerME.find(fullName("PhiResidualRMS"))->second,sect,stat,phiRMS);
237 
238  phiSummary = performLutTest(phiMean,phiRMS,thresholdPhiMean,thresholdPhiRMS);
239 
240  }
241  fillWhPlot(whME[wh].find(fullName("PhiLutSummary"))->second,sect,stat,phiSummary);
242 
243  // Make Phib Residual Summary
244  TH1F * PhibResidual = getHisto<TH1F>(igetter.get(getMEName("PhibResidual","Segment", chId)));
245  int phibSummary = stat==3 ? 0 : 1; // station 3 has no meaningful MB3 phi bending information
246 
247  if (stat != 3 && PhibResidual && PhibResidual->GetEffectiveEntries()>10) {// station 3 has no meaningful MB3 phi bending information
248 
249  // Fill client histos
250  if( whME[wh].find(fullName("PhibResidualMean")) == whME[wh].end() ){
251  bookWheelHistos(ibooker,wh,"PhibResidualMean");
252  bookWheelHistos(ibooker,wh,"PhibResidualRMS");
253  }
254 
255  double peak = PhibResidual->GetBinCenter(PhibResidual->GetMaximumBin());
256  double phibMean = 0;
257  double phibRMS = 0;
258  try {
259  TF1 ffPhib("ffPhib","gaus");
260  PhibResidual->Fit(&ffPhib,"CQO","",peak-5,peak+5);
261  phibMean = ffPhib.GetParameter(1);
262  phibRMS = ffPhib.GetParameter(2);
263  } catch (cms::Exception& iException) {
264  edm::LogError(category()) << "[" << testName << "Test]: Error fitting PhibResidual for Wheel " << wh
265  <<" Sector " << sect << " Station " << stat;
266  }
267 
268  std::map<std::string,MonitorElement*> &innerME = whME[wh];
269  fillWhPlot(innerME.find(fullName("PhibResidualMean"))->second,sect,stat,phibMean);
270  fillWhPlot(innerME.find(fullName("PhibResidualRMS"))->second,sect,stat,phibRMS);
271 
272  phibSummary = performLutTest(phibMean,phibRMS,thresholdPhibMean,thresholdPhibRMS);
273 
274  }
275  fillWhPlot(whME[wh].find(fullName("PhibLutSummary"))->second,sect,stat,phibSummary);
276 
277  }
278  }
279  }
280 
281  // Barrel Summary Plots
282  for (vector<string>::const_iterator iTr = trigSources.begin(); iTr != trigSources.end(); ++iTr){
283  trigSource = (*iTr);
284  for (vector<string>::const_iterator iHw = hwSources.begin(); iHw != hwSources.end(); ++iHw){
285  hwSource = (*iHw);
286  for (int wh=-2; wh<=2; ++wh){
287  std::map<std::string,MonitorElement*> *innerME = &(whME[wh]);
288 
289  TH2F* phiWhSummary = getHisto<TH2F>(innerME->find(fullName("PhiLutSummary"))->second);
290  TH2F* phibWhSummary = getHisto<TH2F>(innerME->find(fullName("PhibLutSummary"))->second);
291  for (int sect=1; sect<=12; ++sect){
292  int phiErr = 0;
293  int phibErr = 0;
294  int phiNoData = 0;
295  int phibNoData = 0;
296  for (int stat=1; stat<=4; ++stat){
297  switch (static_cast<int>(phiWhSummary->GetBinContent(sect,stat))) {
298  case 1:
299  phiNoData++;
300  case 2:
301  case 3:
302  phiErr++;
303  }
304  switch (static_cast<int>(phibWhSummary->GetBinContent(sect,stat))) {
305  case 1:
306  phibNoData++;
307  case 2:
308  case 3:
309  phibErr++;
310  }
311  }
312  if (phiNoData == 4) phiErr = 5;
313  if (phibNoData == 3) phibErr = 5; // MB3 has no phib information
314  cmsME.find(fullName("PhiLutSummary"))->second->setBinContent(sect,wh+wheelArrayShift,phiErr);
315  cmsME.find(fullName("PhibLutSummary"))->second->setBinContent(sect,wh+wheelArrayShift,phibErr);
316  }
317  }
318  }
319  }
320 
321 }
std::vector< std::string > trigSources
const std::vector< const DTChamber * > & chambers() const
Return a vector of all Chamber.
Definition: DTGeometry.cc:102
std::map< int, std::map< std::string, MonitorElement * > > whME
int performLutTest(double mean, double RMS, double thresholdMean, double thresholdRMS)
Perform Lut Test logical operations.
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:20
edm::ESHandle< DTGeometry > muonGeom
void fillWhPlot(MonitorElement *plot, int sect, int stat, float value, bool lessIsBest=true)
Fill summary plots managing double MB4 chambers.
U second(std::pair< T, U > const &p)
MonitorElement * get(std::string const &path)
Definition: DQMStore.cc:303
std::string category()
Get message logger name.
void Bookings(DQMStore::IBooker &, DQMStore::IGetter &)
std::string getMEName(std::string histoTag, std::string subfolder, const DTChamberId &chambid)
Get the ME name (by chamber)
std::map< std::string, MonitorElement * > cmsME
std::vector< std::string > hwSources
int wheel() const
Return the wheel number.
Definition: DTChamberId.h:45
void bookWheelHistos(DQMStore::IBooker &, int wheel, std::string hTag, std::string folder="")
Book the new MEs (for each wheel)
std::string fullName(std::string htype)
Create fullname from histo partial name.

Member Data Documentation

bool DTLocalTriggerLutTest::bookingdone
private

Definition at line 55 of file DTLocalTriggerLutTest.h.

bool DTLocalTriggerLutTest::doCorrStudy
private

Definition at line 53 of file DTLocalTriggerLutTest.h.

double DTLocalTriggerLutTest::thresholdPhibMean
private

Definition at line 51 of file DTLocalTriggerLutTest.h.

double DTLocalTriggerLutTest::thresholdPhibRMS
private

Definition at line 52 of file DTLocalTriggerLutTest.h.

double DTLocalTriggerLutTest::thresholdPhiMean
private

Definition at line 51 of file DTLocalTriggerLutTest.h.

double DTLocalTriggerLutTest::thresholdPhiRMS
private

Definition at line 52 of file DTLocalTriggerLutTest.h.

const int DTLocalTriggerLutTest::wheelArrayShift = 3
protected

Definition at line 41 of file DTLocalTriggerLutTest.h.