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

#include <RPCEfficiencySecond.h>

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

int rollY (std::string shortname, const std::vector< std::string > &rollNames)
 
 RPCEfficiencySecond (const edm::ParameterSet &)
 
 ~RPCEfficiencySecond () override
 
- 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 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< 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 beginJob () override
 
void dqmEndJob (DQMStore::IBooker &, DQMStore::IGetter &) override
 
void dqmEndLuminosityBlock (DQMStore::IBooker &, DQMStore::IGetter &, edm::LuminosityBlock const &, edm::EventSetup const &) override
 
- 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)
 

Private Member Functions

void myBooker (DQMStore::IBooker &)
 

Private Attributes

MonitorElement * DiskSummary [10]
 
MonitorElement * EffDistroD [10]
 
MonitorElement * EffDistroW [5]
 
MonitorElement * EffGlobD [10]
 
MonitorElement * EffGlobW [5]
 
MonitorElement * ExGregR2D [10]
 
MonitorElement * ExGregR3D [10]
 
MonitorElement * ExpLayerW [5]
 
MonitorElement * ExsectorEffW [5]
 
std::string folderPath
 
MonitorElement * GregR2D [10]
 
MonitorElement * GregR3D [10]
 
MonitorElement * histoCSC
 
MonitorElement * histoCSC2
 
MonitorElement * histoDT
 
MonitorElement * histoDT2
 
MonitorElement * histoeffIdRPC_CSC
 
MonitorElement * histoeffIdRPC_DT
 
MonitorElement * histoPRO
 
MonitorElement * histoRealRPC
 
MonitorElement * histoRealRPC2
 
MonitorElement * histoRPC
 
MonitorElement * histoRPC2
 
bool init_
 
int innermostRings_
 
std::map< int, std::map< std::string, MonitorElement * > > meCollection
 
int numberOfDisks_
 
MonitorElement * ObsLayerW [5]
 
MonitorElement * OcGregR2D [10]
 
MonitorElement * OcGregR3D [10]
 
MonitorElement * OcsectorEffW [5]
 
edm::ESHandle< RPCGeometry > rpcGeo_
 
MonitorElement * sectorEffW [5]
 
MonitorElement * WheelSummary [5]
 

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)
 
- Protected Attributes inherited from DQMEDHarvester
edm::EDPutTokenT< DQMToken > lumiToken_
 
edm::EDPutTokenT< DQMToken > runToken_
 

Detailed Description

author Camilo Carrillo (Uniandes)

Definition at line 19 of file RPCEfficiencySecond.h.

Constructor & Destructor Documentation

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

Definition at line 19 of file RPCEfficiencySecond.cc.

References folderPath, edm::ParameterSet::getUntrackedParameter(), init_, innermostRings_, numberOfDisks_, and AlCaHLTBitMon_QueryRunRegistry::string.

19  {
20 
21  folderPath = iConfig.getUntrackedParameter<std::string>("folderPath","RPC/RPCEfficiency/");
22  numberOfDisks_ = iConfig.getUntrackedParameter<int>("NumberOfEndcapDisks", 4);
23  innermostRings_ = iConfig.getUntrackedParameter<int>("NumberOfInnermostEndcapRings", 2);
24 
25  init_ = false;
26 
27 }
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
RPCEfficiencySecond::~RPCEfficiencySecond ( )
override

Definition at line 29 of file RPCEfficiencySecond.cc.

29 {}

Member Function Documentation

void RPCEfficiencySecond::beginJob ( void  )
overrideprotectedvirtual

Reimplemented from edm::one::EDProducerBase.

Definition at line 31 of file RPCEfficiencySecond.cc.

31 {}
void RPCEfficiencySecond::dqmEndJob ( DQMStore::IBooker &  ibooker,
DQMStore::IGetter &  igetter 
)
overrideprotectedvirtual

Implements DQMEDHarvester.

Definition at line 41 of file RPCEfficiencySecond.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), edmScanValgrind::buffer, edmIntegrityCheck::d, rpcdqm::utils::detId2RollNr(), RPCGeometry::dets(), DiskSummary, EffDistroD, EffDistroW, EffGlobD, EffGlobW, ExGregR2D, ExGregR3D, ExpLayerW, ExsectorEffW, MonitorElement::Fill(), printsummarytable::folder, folderPath, RPCBookFolderStructure::folderStructure(), DQMStore::IGetter::get(), MonitorElement::getBinContent(), MonitorElement::getTH1F(), GregR2D, GregR3D, histoCSC, histoDT, histoRPC, mps_fire::i, innermostRings_, createfilelist::int, gen::k, rpcdqm::utils::labelXAxisSector(), rpcdqm::utils::labelXAxisSegment(), rpcdqm::utils::labelYAxisRing(), rpcdqm::utils::labelYAxisRoll(), RPCDetId::layer(), LogDebug, N, RPCGeomServ::name(), nullptr, numberOfDisks_, connectstrParser::o, ObsLayerW, OcGregR2D, OcGregR3D, OcsectorEffW, AlCaHLTBitMon_ParallelJobs::p, alignCSCRings::r, DetId::rawId(), RPCDetId::region(), relativeConstraints::ring, RPCDetId::ring(), RPCChamber::rolls(), rpcGeo_, RPCDetId::sector(), sectorEffW, RPCGeomServ::segment(), MonitorElement::setAxisRange(), MonitorElement::setAxisTitle(), MonitorElement::setBinContent(), MonitorElement::setBinError(), MonitorElement::setBinLabel(), DQMStore::IBooker::setCurrentFolder(), RPCGeomServ::shortname(), mathSSE::sqrt(), relativeConstraints::station, RPCDetId::station(), AlCaHLTBitMon_QueryRunRegistry::string, w, makeMuonMisalignmentScenario::wheel, WheelSummary, x, y, and cuy::yBin.

41  {
42 
43 
44  rpcdqm::utils rpcUtils;
45  //Barrel
46  std::stringstream folderName;
47  std::stringstream MeName; std::stringstream MeTitle;
48  for (int w =-2; w<=2; w++){
49  folderName.str("");
50  folderName<<folderPath<<"Wheel_"<<w;
51  ibooker.setCurrentFolder(folderName.str());
52  MeName.str("");
53  MeName<<"EffDistroWheel_"<<w;
54  MeTitle.str("");
55  MeTitle<<"Efficiency Distribution for Wheel "<<w;
56  EffDistroW[w+2]=ibooker.book1D(MeName.str(),MeTitle.str(),20,0.5,100.5);
57  MeName.str("");
58  MeName<<"GlobEfficiencyWheel_"<<w;
59  MeTitle.str("");
60  MeTitle<<"Efficiency Wheel "<<w;
61  EffGlobW[w+2]=ibooker.book1D(MeName.str(),MeTitle.str(),206,0.5,206.5);
62  ibooker.setCurrentFolder(folderPath);
63  MeName.str("");
64  MeName<<"Efficiency_Roll_vs_Sector_Wheel_"<<w;
65  WheelSummary[w+2] = ibooker.book2D(MeName.str(),MeName.str(),12, 0.5,12.5, 21, 0.5, 21.5);
66  rpcUtils.labelXAxisSector( WheelSummary[w+2] );
67  rpcUtils.labelYAxisRoll( WheelSummary[w+2], 0, w, true);
68  ibooker.setCurrentFolder(folderPath+"Azimutal/");
69  MeName.str("");
70  MeName<<"AzimutalDistroW"<<w;
71  MeTitle.str("");
72  MeTitle<<"Efficiency per Sector Wheel "<<w;
73  sectorEffW[w+2]= ibooker.book1D(MeName.str(),MeTitle.str(),12,0.5,12.5);
74  MeName.str("");
75  MeName<<"AzimutalDistroW"<<w<<"Ex";
76  MeTitle.str("");
77  MeTitle<<"Expected per Sector Wheel "<<w;
78  ExsectorEffW[w+2]= ibooker.book1D(MeName.str(),MeTitle.str(),12,0.5,12.5);
79  MeName.str("");
80  MeName<<"AzimutalDistroW"<<w<<"Oc";
81  MeTitle.str("");
82  MeTitle<<"Occupancy per Sector Wheel "<<w;
83  OcsectorEffW[w+2]= ibooker.book1D(MeName.str(),MeTitle.str(),12,0.5,12.5);
84  ibooker.setCurrentFolder(folderPath+"BarrelPerLayer/");
85  MeName.str("");
86  MeName<<"ExpLayerW"<<w;
87  MeTitle.str("");
88  MeTitle<<"Expected Wheel "<<w;
89  ExpLayerW[w+2]= ibooker.book1D(MeName.str(),MeTitle.str(),6,0.5,6.5);
90  MeName.str("");
91  MeName<<"ObsLayerW"<<w;
92  MeTitle.str("");
93  MeTitle<<"Observed Wheel "<<w;
94  ObsLayerW[w+2]= ibooker.book1D(MeName.str(),MeTitle.str(),6,0.5,6.5);
95  }
96  //EndCap
97  int index = 0;
98  for (int d = (-1 *numberOfDisks_); d<=numberOfDisks_; d++){
99  if (d==0) {continue;}
100  folderName.str("");
101  folderName<<folderPath<<"Disk_"<<d;
102  ibooker.setCurrentFolder(folderName.str());
103  MeName.str("");
104  MeName<<"EffDistroDisk_"<<d;
105  MeTitle.str("");
106  MeTitle<<"Efficiency Distribution Disk "<<d;
107  EffDistroD[index]=ibooker.book1D(MeName.str(),MeTitle.str(),20,0.5,100.5);
108  MeName.str("");
109  MeName<<"GlobEfficiencyDisk_"<<d;
110  MeTitle.str("");
111  MeTitle<<"Efficiency Disk "<<d;
112  EffGlobD[index] = ibooker.book1D(MeName.str(),MeTitle.str(),218,0.5,218.5);
113  ibooker.setCurrentFolder(folderPath);
114  MeName.str("");
115  MeName<<"Efficiency_Roll_vs_Segment_Disk_"<<d;
116  DiskSummary[index] = ibooker.book2D(MeName.str(),MeName.str(),36,0.5,36.5,6,0.5,6.5);
117  rpcUtils.labelXAxisSegment(DiskSummary[index]);
118  rpcUtils.labelYAxisRing(DiskSummary[index], innermostRings_ , true);
119  ibooker.setCurrentFolder(folderPath+"Azimutal/");
120  MeName.str("");
121  MeName<<"GregDistroR2D"<<d;
122  MeTitle.str("");
123  MeTitle<<"Efficiency for Station "<<d<<" Ring 2";
124  GregR2D[index]= ibooker.book1D(MeName.str(),MeTitle.str(),36,0.5,36.5);
125  MeName.str("");
126  MeName<<"GregDistroR3D"<<d;
127  MeTitle.str("");
128  MeTitle<<"Efficiency for Station "<<d<<" Ring 3";
129  GregR3D[index]= ibooker.book1D(MeName.str(),MeTitle.str(),36,0.5,36.5);
130  MeName.str("");
131  MeName<<"OcGregDistroR2D"<<d;
132  MeTitle.str("");
133  MeTitle<<"Occupancy Distribution for Station "<<d<<" Ring 2";
134  OcGregR2D[index]= ibooker.book1D(MeName.str(),MeTitle.str(),36,0.5,36.5);
135  MeName.str("");
136  MeName<<"OcGregDistroR3D"<<d;
137  MeTitle.str("");
138  MeTitle<<"Occupancy Distribution for Station "<<d<<" Ring 3";
139  OcGregR3D[index]= ibooker.book1D(MeName.str(),MeTitle.str(),36,0.5,36.5);
140  MeName.str("");
141  MeName<<"ExGregDistroR2D"<<d;
142  MeTitle.str("");
143  MeTitle<<"Expected Distribution for Station "<<d<<" Ring 2";
144  ExGregR2D[index]= ibooker.book1D(MeName.str(),MeTitle.str(),36,0.5,36.5);
145  MeName.str("");
146  MeName<<"ExGregDistroR3D"<<d;
147  MeTitle.str("");
148  MeTitle<<"Expected Distribution for Station "<<d<<" Ring 3";
149  ExGregR3D[index]= ibooker.book1D(MeName.str(),MeTitle.str(),36,0.5,36.5);
150  index++;
151  }
152 
153 
154 
155  LogDebug("rpcefficiencysecond")<<"Getting the RPC Geometry";
156 
157 
158  //Setting Labels in Summary Label.
159  std::stringstream binLabel;
160  for (int w = -2; w<=2 ;w++){
161  for(int x = 1;x<=12;x++){
162  for(int y = 1;y<=21;y++){
163  WheelSummary[w+2]->setBinContent(x,y,-1);
164  }
165  }
166  }
167  for (int d = 0 ; d<(numberOfDisks_*2); d++){
168  for(int x = 1;x<=36;x++){
169  for(int y = 1;y<=6;y++){
170  DiskSummary[d]->setBinContent(x,y,-1);
171  }
172  }
173  }
174  binLabel.str("");
175  int indexWheel[5];
176  for(int j=0;j<5;j++){
177  indexWheel[j]=0;
178  }
179  int indexDisk[10];
180  for(int j=0;j<10;j++){
181  indexDisk[j]=0;
182  }
183  for(TrackingGeometry::DetContainer::const_iterator it=rpcGeo_->dets().begin(); it!=rpcGeo_->dets().end();it++){
184  if(dynamic_cast< const RPCChamber* >( *it ) != nullptr ){
185  const RPCChamber* ch = dynamic_cast< const RPCChamber* >( *it );
186  std::vector< const RPCRoll*> roles = (ch->rolls());
187  for(std::vector<const RPCRoll*>::const_iterator r = roles.begin();r != roles.end(); ++r){
188  RPCDetId rpcId = (*r)->id();
189  RPCGeomServ rpcsrv(rpcId);
190  // std::string nameRoll = rpcsrv.name();
191  std::string camera = rpcsrv.name();
192  //Breaking down the geometry
193  int region = rpcId.region();
194  int wheel = rpcId.ring(); int ring = rpcId.ring();
195  int sector = rpcId.sector();
196  int station = rpcId.station();
197  int geolayer = rpcId.layer();
198 
199  RPCBookFolderStructure * folderStr = new RPCBookFolderStructure();
200  std::string folder = folderPath + "MuonSegEff/" + folderStr->folderStructure(rpcId);
201 
202 
203  if(region==0){//Barrel
204  std::stringstream meIdRPC, meIdDT; //, bxDistroId;
205 
206 
207  meIdRPC<<folder<<"/RPCDataOccupancyFromDT_"<<rpcId.rawId();
208  meIdDT<<folder<<"/ExpectedOccupancyFromDT_"<<rpcId.rawId();
209  histoRPC = igetter.get(meIdRPC.str());
210  histoDT = igetter.get(meIdDT.str());
211  int NumberWithOutPrediction=0;
212  double p = 0.;
213  double o = 0.;
214  double ef = 0.;
215  double er = 0.;
216  double buffef = 0.;
217  double buffer = 0.;
218  double sumbuffef = 0.;
219  double sumbuffer = 0.;
220  double averageeff = 0.;
221  int NumberStripsPointed = 0;
222  if(histoRPC && histoDT){
223  for(int i=1;i<=int((*r)->nstrips());++i){
224  if(histoDT->getBinContent(i)!=0){
225  LogDebug("rpcefficiencysecond")<<"Inside the If";
226  buffef = double(histoRPC->getBinContent(i))/double(histoDT->getBinContent(i));
227  buffer = sqrt(buffef*(1.-buffef)/double(histoDT->getBinContent(i)));
228  sumbuffef=sumbuffef+buffef;
229  sumbuffer = sumbuffer + buffer*buffer;
230  NumberStripsPointed++;
231  }else{
232  NumberWithOutPrediction++;
233  }
234  LogDebug("rpcefficiencysecond")<<"Strip="<<i<<" RPC="<<histoRPC->getBinContent(i)<<" DT="<<histoDT->getBinContent(i)<<" buffef="<<buffef<<" sumbuffef="<<sumbuffef<<" NumberStripsPointed="<<NumberStripsPointed<<" NumberWithOutPrediction"<<NumberWithOutPrediction;
235  }
236  p=histoDT->getTH1F()->Integral();
237  o=histoRPC->getTH1F()->Integral();
238  if(NumberStripsPointed!=0){
239  averageeff = (sumbuffef/double(NumberStripsPointed))*100.;
240  //averageerr = sqrt(sumbuffer/double(NumberStripsPointed))*100.;
241  }
242  }
243  // int Ring = rpcId.ring();
244  if(p!=0){
245  ef = double(o)/double(p);
246  er = sqrt(ef*(1.-ef)/double(p));
247  }
248  ef=ef*100;
249  er=er*100;
250  //Filling azimutal Wheel Histograms
251  int layer = 0;
252  if(station==1&&geolayer==1) layer = 1;
253  else if(station==1&&geolayer==2) layer = 2;
254  else if(station==2&&geolayer==1) layer = 3;
255  else if(station==2&&geolayer==2) layer = 4;
256  else if(station==3) layer = 5;
257  else if(station==4) layer = 6;
258  ExsectorEffW[wheel+2]->Fill(sector,p); OcsectorEffW[wheel+2]->Fill(sector,o);
259  ExpLayerW[wheel+2]->Fill(layer, p); ObsLayerW[wheel+2]->Fill(layer, o);
260  LogDebug("rpcefficiencysecond")<<"Pigi "<<camera<<" "<<rpcsrv.shortname()<<" "<<(*r)->id();
261  if(p > 100){//We need at least 100 predictions to fill the summary plot
262  int xBin,yBin;
263  xBin= (*r)->id().sector();
264  rpcdqm::utils rollNumber;
265  yBin = rollNumber.detId2RollNr((*r)->id());
266  WheelSummary[wheel+2]->setBinContent(xBin,yBin,averageeff);
267  }
268  EffDistroW[wheel+2]->Fill(averageeff);
269  indexWheel[wheel+2]++;
270  EffGlobW[wheel+2]->setBinContent(indexWheel[wheel+2],ef);
271  EffGlobW[wheel+2]->setBinError(indexWheel[wheel+2],er);
272  EffGlobW[wheel+2]->setBinLabel(indexWheel[wheel+2],camera,1);
273  }else{//EndCap
274  std::stringstream meIdRPC,meIdCSC; //, bxDistroId;
275  std::string meIdPRO;
276  meIdRPC<<folder<<"/RPCDataOccupancyFromCSC_"<<rpcId.rawId();
277  meIdCSC<<folder<<"/ExpectedOccupancyFromCSC_"<<rpcId.rawId();
278  histoRPC= igetter.get(meIdRPC.str());
279  histoCSC= igetter.get(meIdCSC.str());
280 
281  int NumberWithOutPrediction=0;
282  double p = 0;
283  double o = 0;
284  double ef =0;
285  double er =0;
286  double buffef = 0;
287  double buffer = 0;
288  double sumbuffef = 0;
289  double sumbuffer = 0;
290  double averageeff = 0;
291  int NumberStripsPointed = 0;
292  if(histoRPC && histoCSC) {
293  LogDebug("rpcefficiencysecond")<<rpcsrv.name();
294  for(int i=1;i<=int((*r)->nstrips());++i){
295  if(histoCSC->getBinContent(i)!=0){
296  LogDebug("rpcefficiencysecond")<<"Inside the If";
297  buffef = double(histoRPC->getBinContent(i))/double(histoCSC->getBinContent(i));
298  buffer = sqrt(buffef*(1.-buffef)/double(histoCSC->getBinContent(i)));
299  sumbuffef=sumbuffef+buffef;
300  sumbuffer = sumbuffer + buffer*buffer;
301  NumberStripsPointed++;
302  }else{
303  NumberWithOutPrediction++;
304  }
305  LogDebug("rpcefficiencysecond")<<"Strip="<<i<<" RPC="<<histoRPC->getBinContent(i)<<" CSC="<<histoCSC->getBinContent(i)<<" buffef="<<buffef<<" buffer="<<buffer<<" sumbuffef="<<sumbuffef<<" sumbuffer="<<sumbuffer<<" NumberStripsPointed="<<NumberStripsPointed<<" NumberWithOutPrediction"<<NumberWithOutPrediction;
306  }
307  p=histoCSC->getTH1F()->Integral();
308  o=histoRPC->getTH1F()->Integral();
309  if(NumberStripsPointed!=0){
310  averageeff = (sumbuffef/double(NumberStripsPointed))*100.;
311  }
312  }
313  int Disk = station*region;
314  if(Disk > numberOfDisks_ || Disk<-numberOfDisks_){continue;} //remove strange disk numbers!!
315  int dIndex = -1;
316  if(region == -1){
317  dIndex = Disk+numberOfDisks_;
318  }else if(region == 1){
319  dIndex = Disk+numberOfDisks_-1;
320  }
321  if( dIndex<0 || dIndex>= numberOfDisks_*2){continue;} //extra check on disk numering
322  if(p!=0){
323  ef = double(o)/double(p);
324  er = sqrt(ef*(1.-ef)/double(p));
325  }
326  ef=ef*100;
327  er=er*100;
328  if(ring==2){
329  if (ExGregR2D[dIndex] && OcGregR2D[dIndex] ){
330  ExGregR2D[dIndex]->Fill(rpcsrv.segment(),p);
331  OcGregR2D[dIndex]->Fill(rpcsrv.segment(),o);
332  }
333  } else if(ring==3){
334  ExGregR3D[dIndex]->Fill(rpcsrv.segment(),p);
335  OcGregR3D[dIndex]->Fill(rpcsrv.segment(),o);
336  }
337  if(p > 100){ //We need at least 100 predictions to fill the summary plot
338  RPCGeomServ RPCServ((*r)->id());
339  int xBin = RPCServ.segment();
340  int yBin= (ring-1)*3-(*r)->id().roll()+1;
341  DiskSummary[dIndex]->setBinContent(xBin, yBin, averageeff);
342  }
343  if ( EffDistroD[dIndex]){EffDistroD[dIndex]->Fill(averageeff);}
344  indexDisk[dIndex]++;
345  if ( EffGlobD[dIndex]){
346  EffGlobD[dIndex]->setBinContent(indexDisk[dIndex],ef);
347  EffGlobD[dIndex]->setBinError(indexDisk[dIndex],er);
348  EffGlobD[dIndex]->setBinLabel(indexDisk[dIndex],camera,1);
349  }
350  }
351  delete folderStr;
352  }
353  }
354 
355 
356  }
357  double eff,N,err;
358  int k;
359  for (int d = 0; d<(numberOfDisks_*2); d++){
360  for(k=1;k<=36;k++){
361  err=0; eff=0; N=ExGregR2D[d]->getBinContent(k);
362  if(N!=0.){ eff = OcGregR2D[d]->getBinContent(k)/N; err=sqrt(eff*(1-eff)/N);}
363  GregR2D[d]->setBinContent(k,eff); GregR2D[d]->setBinError(k,err);
364  err=0; eff=0; N=ExGregR3D[d]->getBinContent(k);
365  if(N!=0.){ eff = OcGregR3D[d]->getBinContent(k)/N; err=sqrt(eff*(1-eff)/N);}
366  GregR3D[d]->setBinContent(k,eff); GregR3D[d]->setBinError(k,err);
367  }
368  }
369  for (int w =-2; w<=2; w++){
370  for(k=1;k<=12;k++){
371  err=0; eff=0; N=ExsectorEffW[w+2]->getBinContent(k);
372  if(N!=0.){ eff = OcsectorEffW[w+2]->getBinContent(k)/N;err=sqrt(eff*(1-eff)/N);}
373  sectorEffW[w+2]->setBinContent(k,eff); sectorEffW[w+2]->setBinError(k,err);
374  }
375  }
376  //Ranges for Both
377  //Barrel
378  for (int w=-2; w<=2; w++){
379  EffGlobW[w+2]->setAxisRange(-4.,100.,2);
380  EffGlobW[w+2]->setAxisTitle("%",2);
381  }
382  for (int d=0; d<(numberOfDisks_*2); d++){
383  EffGlobD[d]->setAxisRange(-4.,100.,2);
384  EffGlobD[d]->setAxisTitle("%",2);
385  }
386 
387 
388 }
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:670
RPCEfficiencySecond::EffGlobW
MonitorElement * EffGlobW[5]
Definition: RPCEfficiencySecond.h:51
MonitorElement::setBinContent
void setBinContent(int binx, double content)
set content of bin (1-D)
Definition: MonitorElement.cc:866
mps_fire.i
i
Definition: mps_fire.py:338
AlCaHLTBitMon_ParallelJobs.p
p
Definition: AlCaHLTBitMon_ParallelJobs.py:153
w
const double w
Definition: UKUtility.cc:23
RPCEfficiencySecond::histoCSC
MonitorElement * histoCSC
Definition: RPCEfficiencySecond.h:36
MonitorElement::setAxisRange
void setAxisRange(double xmin, double xmax, int axis=1)
set x-, y- or z-axis range (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:964
MonitorElement::getTH1F
TH1F * getTH1F() const
Definition: MonitorElement.cc:1483
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
nullptr
#define nullptr
Definition: GCC11Compatibility.h:37
DetId::rawId
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:50
RPCEfficiencySecond::OcGregR3D
MonitorElement * OcGregR3D[10]
Definition: RPCEfficiencySecond.h:75
RPCEfficiencySecond::histoDT
MonitorElement * histoDT
Definition: RPCEfficiencySecond.h:34
MonitorElement::setBinLabel
void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:942
RPCEfficiencySecond::OcsectorEffW
MonitorElement * OcsectorEffW[5]
Definition: RPCEfficiencySecond.h:68
rpcdqm::utils::labelXAxisSegment
void labelXAxisSegment(MonitorElement *myMe)
Definition: utils.h:264
rpcdqm::utils::labelYAxisRoll
void labelYAxisRoll(MonitorElement *myMe, int region, int ring, bool useRollInfo)
Definition: utils.h:283
RPCEfficiencySecond::sectorEffW
MonitorElement * sectorEffW[5]
Definition: RPCEfficiencySecond.h:67
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:166
RPCEfficiencySecond::rpcGeo_
edm::ESHandle< RPCGeometry > rpcGeo_
Definition: RPCEfficiencySecond.h:83
Disk
Definition: BoundDisk.h:19
rpcdqm::utils::labelYAxisRing
void labelYAxisRing(MonitorElement *myMe, int numberOfRings, bool useRollInfo)
Definition: utils.h:311
RPCEfficiencySecond::ObsLayerW
MonitorElement * ObsLayerW[5]
Definition: RPCEfficiencySecond.h:81
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:268
RPCDetId::ring
int ring() const
Definition: RPCDetId.h:72
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:106
RPCEfficiencySecond::EffGlobD
MonitorElement * EffGlobD[10]
Definition: RPCEfficiencySecond.h:57
RPCChamber::rolls
const std::vector< const RPCRoll * > & rolls() const
Return the Rolls.
Definition: RPCChamber.cc:68
RPCEfficiencySecond::EffDistroD
MonitorElement * EffDistroD[10]
Definition: RPCEfficiencySecond.h:60
RPCEfficiencySecond::ExpLayerW
MonitorElement * ExpLayerW[5]
Definition: RPCEfficiencySecond.h:80
MonitorElement::setBinError
void setBinError(int binx, double error)
set uncertainty on content of bin (1-D)
Definition: MonitorElement.cc:891
RPCEfficiencySecond::OcGregR2D
MonitorElement * OcGregR2D[10]
Definition: RPCEfficiencySecond.h:74
rpcdqm::utils::labelXAxisSector
void labelXAxisSector(MonitorElement *myMe)
Definition: utils.h:249
DQMStore::IGetter::get
MonitorElement * get(std::string const &path)
Definition: DQMStore.cc:303
gen::k
int k[5][pyjets_maxn]
Definition: Cascade2Hadronizer.cc:79
RPCDetId::layer
int layer() const
Definition: RPCDetId.h:108
RPCGeomServ::segment
virtual int segment()
Definition: RPCGeomServ.cc:469
rpcdqm::utils::detId2RollNr
int detId2RollNr(const RPCDetId &_id)
Definition: utils.h:31
N
#define N
Definition: blowfish.cc:9
RPCBookFolderStructure::folderStructure
std::string folderStructure(RPCDetId detId)
Definition: RPCBookFolderStructure.h:18
DQMStore::IBooker::book2D
MonitorElement * book2D(Args &&...args)
Definition: DQMStore.h:109
RPCGeometry::dets
const DetContainer & dets() const override
Returm a vector of all GeomDet (including all GeomDetUnits)
Definition: RPCGeometry.cc:33
RPCEfficiencySecond::WheelSummary
MonitorElement * WheelSummary[5]
Definition: RPCEfficiencySecond.h:63
RPCEfficiencySecond::ExGregR2D
MonitorElement * ExGregR2D[10]
Definition: RPCEfficiencySecond.h:77
RPCDetId::sector
int sector() const
Sector id: the group of chambers at same phi (and increasing r)
Definition: RPCDetId.h:102
RPCEfficiencySecond::DiskSummary
MonitorElement * DiskSummary[10]
Definition: RPCEfficiencySecond.h:64
MonitorElement::getBinContent
double getBinContent(int binx) const
get content of bin (1-D)
Definition: MonitorElement.cc:765
RPCEfficiencySecond::ExsectorEffW
MonitorElement * ExsectorEffW[5]
Definition: RPCEfficiencySecond.h:69
RPCEfficiencySecond::GregR2D
MonitorElement * GregR2D[10]
Definition: RPCEfficiencySecond.h:71
RPCEfficiencySecond::ExGregR3D
MonitorElement * ExGregR3D[10]
Definition: RPCEfficiencySecond.h:78
MonitorElement::setAxisTitle
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:973
RPCEfficiencySecond::histoRPC
MonitorElement * histoRPC
Definition: RPCEfficiencySecond.h:33
RPCEfficiencySecond::GregR3D
MonitorElement * GregR3D[10]
Definition: RPCEfficiencySecond.h:72
cuy.yBin
yBin
Definition: cuy.py:893
RPCEfficiencySecond::EffDistroW
MonitorElement * EffDistroW[5]
Definition: RPCEfficiencySecond.h:54
RPCDetId::region
int region() const
Region id: 0 for Barrel, +/-1 For +/- Endcap.
Definition: RPCDetId.h:63
RPCDetId::station
int station() const
Definition: RPCDetId.h:96
void RPCEfficiencySecond::dqmEndLuminosityBlock ( DQMStore::IBooker &  ibooker,
DQMStore::IGetter &  igetter,
edm::LuminosityBlock const &  lb,
edm::EventSetup const &  iSetup 
)
overrideprotectedvirtual

Reimplemented from DQMEDHarvester.

Definition at line 33 of file RPCEfficiencySecond.cc.

References edm::EventSetup::get(), init_, LogDebug, and rpcGeo_.

33  {
34 
35  LogDebug("rpcefficiencysecond")<<"Getting the RPC Geometry";
36  iSetup.get<MuonGeometryRecord>().get(rpcGeo_);
37  init_= true;
38 }
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:670
RPCEfficiencySecond::rpcGeo_
edm::ESHandle< RPCGeometry > rpcGeo_
Definition: RPCEfficiencySecond.h:83
void RPCEfficiencySecond::myBooker ( DQMStore::IBooker &  )
private
int RPCEfficiencySecond::rollY ( std::string  shortname,
const std::vector< std::string > &  rollNames 
)

Member Data Documentation

MonitorElement* RPCEfficiencySecond::DiskSummary[10]
private

Definition at line 64 of file RPCEfficiencySecond.h.

Referenced by dqmEndJob().

MonitorElement* RPCEfficiencySecond::EffDistroD[10]
private

Definition at line 60 of file RPCEfficiencySecond.h.

Referenced by dqmEndJob().

MonitorElement* RPCEfficiencySecond::EffDistroW[5]
private

Definition at line 54 of file RPCEfficiencySecond.h.

Referenced by dqmEndJob().

MonitorElement* RPCEfficiencySecond::EffGlobD[10]
private

Definition at line 57 of file RPCEfficiencySecond.h.

Referenced by dqmEndJob().

MonitorElement* RPCEfficiencySecond::EffGlobW[5]
private

Definition at line 51 of file RPCEfficiencySecond.h.

Referenced by dqmEndJob().

MonitorElement* RPCEfficiencySecond::ExGregR2D[10]
private

Definition at line 77 of file RPCEfficiencySecond.h.

Referenced by dqmEndJob().

MonitorElement* RPCEfficiencySecond::ExGregR3D[10]
private

Definition at line 78 of file RPCEfficiencySecond.h.

Referenced by dqmEndJob().

MonitorElement* RPCEfficiencySecond::ExpLayerW[5]
private

Definition at line 80 of file RPCEfficiencySecond.h.

Referenced by dqmEndJob().

MonitorElement* RPCEfficiencySecond::ExsectorEffW[5]
private

Definition at line 69 of file RPCEfficiencySecond.h.

Referenced by dqmEndJob().

std::string RPCEfficiencySecond::folderPath
private

Definition at line 91 of file RPCEfficiencySecond.h.

Referenced by dqmEndJob(), and RPCEfficiencySecond().

MonitorElement* RPCEfficiencySecond::GregR2D[10]
private

Definition at line 71 of file RPCEfficiencySecond.h.

Referenced by dqmEndJob().

MonitorElement* RPCEfficiencySecond::GregR3D[10]
private

Definition at line 72 of file RPCEfficiencySecond.h.

Referenced by dqmEndJob().

MonitorElement* RPCEfficiencySecond::histoCSC
private

Definition at line 36 of file RPCEfficiencySecond.h.

Referenced by dqmEndJob().

MonitorElement* RPCEfficiencySecond::histoCSC2
private

Definition at line 46 of file RPCEfficiencySecond.h.

MonitorElement* RPCEfficiencySecond::histoDT
private

Definition at line 34 of file RPCEfficiencySecond.h.

Referenced by dqmEndJob().

MonitorElement* RPCEfficiencySecond::histoDT2
private

Definition at line 44 of file RPCEfficiencySecond.h.

MonitorElement* RPCEfficiencySecond::histoeffIdRPC_CSC
private

Definition at line 40 of file RPCEfficiencySecond.h.

MonitorElement* RPCEfficiencySecond::histoeffIdRPC_DT
private

Definition at line 39 of file RPCEfficiencySecond.h.

MonitorElement* RPCEfficiencySecond::histoPRO
private

Definition at line 38 of file RPCEfficiencySecond.h.

MonitorElement* RPCEfficiencySecond::histoRealRPC
private

Definition at line 35 of file RPCEfficiencySecond.h.

MonitorElement* RPCEfficiencySecond::histoRealRPC2
private

Definition at line 45 of file RPCEfficiencySecond.h.

MonitorElement* RPCEfficiencySecond::histoRPC
private

Definition at line 33 of file RPCEfficiencySecond.h.

Referenced by dqmEndJob().

MonitorElement* RPCEfficiencySecond::histoRPC2
private

Definition at line 43 of file RPCEfficiencySecond.h.

bool RPCEfficiencySecond::init_
private

Definition at line 89 of file RPCEfficiencySecond.h.

Referenced by dqmEndLuminosityBlock(), and RPCEfficiencySecond().

int RPCEfficiencySecond::innermostRings_
private

Definition at line 93 of file RPCEfficiencySecond.h.

Referenced by dqmEndJob(), and RPCEfficiencySecond().

std::map<int, std::map<std::string, MonitorElement*> > RPCEfficiencySecond::meCollection
private

Definition at line 87 of file RPCEfficiencySecond.h.

int RPCEfficiencySecond::numberOfDisks_
private

Definition at line 92 of file RPCEfficiencySecond.h.

Referenced by dqmEndJob(), and RPCEfficiencySecond().

MonitorElement* RPCEfficiencySecond::ObsLayerW[5]
private

Definition at line 81 of file RPCEfficiencySecond.h.

Referenced by dqmEndJob().

MonitorElement* RPCEfficiencySecond::OcGregR2D[10]
private

Definition at line 74 of file RPCEfficiencySecond.h.

Referenced by dqmEndJob().

MonitorElement* RPCEfficiencySecond::OcGregR3D[10]
private

Definition at line 75 of file RPCEfficiencySecond.h.

Referenced by dqmEndJob().

MonitorElement* RPCEfficiencySecond::OcsectorEffW[5]
private

Definition at line 68 of file RPCEfficiencySecond.h.

Referenced by dqmEndJob().

edm::ESHandle<RPCGeometry> RPCEfficiencySecond::rpcGeo_
private

Definition at line 83 of file RPCEfficiencySecond.h.

Referenced by dqmEndJob(), and dqmEndLuminosityBlock().

MonitorElement* RPCEfficiencySecond::sectorEffW[5]
private

Definition at line 67 of file RPCEfficiencySecond.h.

Referenced by dqmEndJob().

MonitorElement* RPCEfficiencySecond::WheelSummary[5]
private

Definition at line 63 of file RPCEfficiencySecond.h.

Referenced by dqmEndJob().

Generated for CMSSW Reference Manual by doxygen 1.8.11