test
CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Private Attributes
RPCRecHitValidClient Class Reference

#include <RPCRecHitValidClient.h>

Inheritance diagram for RPCRecHitValidClient:
edm::EDAnalyzer edm::EDConsumerBase

Public Member Functions

void analyze (const edm::Event &event, const edm::EventSetup &eventSetup)
 
void endRun (const edm::Run &run, const edm::EventSetup &eventSetup)
 
 RPCRecHitValidClient (const edm::ParameterSet &pset)
 
 ~RPCRecHitValidClient ()
 
- Public Member Functions inherited from edm::EDAnalyzer
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzer ()
 
ModuleDescription const & moduleDescription () const
 
std::string workerType () const
 
virtual ~EDAnalyzer ()
 
- Public Member Functions inherited from edm::EDConsumerBase
 EDConsumerBase ()
 
ProductHolderIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
std::vector
< ProductHolderIndexAndSkipBit >
const & 
itemsToGetFromEvent () const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
 
bool registeredToConsume (ProductHolderIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
 
virtual ~EDConsumerBase ()
 

Private Attributes

std::string subDir_
 

Additional Inherited Members

- Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer ModuleType
 
- Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &)
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Detailed Description

Definition at line 18 of file RPCRecHitValidClient.h.

Constructor & Destructor Documentation

RPCRecHitValidClient::RPCRecHitValidClient ( const edm::ParameterSet pset)
RPCRecHitValidClient::~RPCRecHitValidClient ( )
inline

Definition at line 22 of file RPCRecHitValidClient.h.

22 {};

Member Function Documentation

void RPCRecHitValidClient::analyze ( const edm::Event event,
const edm::EventSetup eventSetup 
)
inlinevirtual

Implements edm::EDAnalyzer.

Definition at line 24 of file RPCRecHitValidClient.h.

24 {};
void RPCRecHitValidClient::endRun ( const edm::Run run,
const edm::EventSetup eventSetup 
)
virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 17 of file RPCRecHitValidClient.cc.

References newFWLiteAna::bin, DQMStore::book1D(), alignCSCRings::e, MonitorElement::Fill(), DQMStore::get(), MonitorElement::getTH1(), min(), nEvent, cppFunctionSkipper::operator, and DQMStore::setCurrentFolder().

18 {
20  if ( !dbe ) return;
21 
22  dbe->setCurrentFolder(subDir_);
23  MEP me_rollEfficiencyBarrel_eff = dbe->book1D("RollEfficiencyBarrel_eff", "Roll efficiency in Barrel;Efficiency [%]", 50+2, -2, 100+2);
24  MEP me_rollEfficiencyEndcap_eff = dbe->book1D("RollEfficiencyEndcap_eff", "Roll efficiency in Endcap;Efficiency [%]", 50+2, -2, 100+2);
25  MEP me_rollEfficiencyStatCutOffBarrel_eff = dbe->book1D("RollEfficiencyCutOffBarrel_eff", "Roll efficiency in Barrel without low stat chamber;Efficiency [%]", 50+2, -2, 100+2);
26  MEP me_rollEfficiencyStatCutOffEndcap_eff = dbe->book1D("RollEfficiencyCutOffEndcap_eff", "Roll efficiency in Endcap without low stat chamber;Efficiency [%]", 50+2, -2, 100+2);
27 
28  const double maxNoise = 1e-7;
29  MEP me_rollNoiseBarrel_noise = dbe->book1D("RollNoiseBarrel_noise", "Roll noise in Barrel;Noise level [Event^{-1}cm^{-2}]", 25+2, -maxNoise/25, maxNoise+maxNoise/25);
30  MEP me_rollNoiseEndcap_noise = dbe->book1D("RollNoiseEndcap_noise", "Roll noise in Endcap;Noise level [Event^{-1}cm^{-2}]", 25+2, -maxNoise/25, maxNoise+maxNoise/25);
31 
32  MEP me_matchOccupancyBarrel_detId = dbe->get(subDir_+"/Occupancy/MatchOccupancyBarrel_detId");
33  MEP me_matchOccupancyEndcap_detId = dbe->get(subDir_+"/Occupancy/MatchOccupancyEndcap_detId");
34  MEP me_refOccupancyBarrel_detId = dbe->get(subDir_+"/Occupancy/RefOccupancyBarrel_detId");
35  MEP me_refOccupancyEndcap_detId = dbe->get(subDir_+"/Occupancy/RefOccupancyEndcap_detId");
36 
37  if ( me_matchOccupancyBarrel_detId and me_refOccupancyBarrel_detId )
38  {
39  TH1* h_matchOccupancyBarrel_detId = me_matchOccupancyBarrel_detId->getTH1();
40  TH1* h_refOccupancyBarrel_detId = me_refOccupancyBarrel_detId->getTH1();
41 
42  for ( int bin = 1, nBin = h_matchOccupancyBarrel_detId->GetNbinsX(); bin <= nBin; ++bin )
43  {
44  const double nRec = h_matchOccupancyBarrel_detId->GetBinContent(bin);
45  const double nRef = h_refOccupancyBarrel_detId->GetBinContent(bin);
46 
47  const double eff = nRef ? nRec/nRef*100 : -1;
48 
49  me_rollEfficiencyBarrel_eff->Fill(eff);
50  if ( nRef >= 20 ) me_rollEfficiencyStatCutOffBarrel_eff->Fill(eff);
51  }
52  }
53 
54  if ( me_matchOccupancyEndcap_detId and me_refOccupancyEndcap_detId )
55  {
56  TH1* h_matchOccupancyEndcap_detId = me_matchOccupancyEndcap_detId->getTH1();
57  TH1* h_refOccupancyEndcap_detId = me_refOccupancyEndcap_detId->getTH1();
58 
59  for ( int bin = 1, nBin = h_matchOccupancyEndcap_detId->GetNbinsX(); bin <= nBin; ++bin )
60  {
61  const double nRec = h_matchOccupancyEndcap_detId->GetBinContent(bin);
62  const double nRef = h_refOccupancyEndcap_detId->GetBinContent(bin);
63 
64  const double eff = nRef ? nRec/nRef*100 : -1;
65 
66  me_rollEfficiencyEndcap_eff->Fill(eff);
67  if ( nRef >= 20 ) me_rollEfficiencyStatCutOffEndcap_eff->Fill(eff);
68  }
69  }
70 
71  MEP me_eventCount = dbe->get(subDir_+"/Occupancy/EventCount");
72  const double nEvent = me_eventCount ? me_eventCount->getTH1()->GetBinContent(1) : 1;
73  MEP me_noiseOccupancyBarrel_detId = dbe->get(subDir_+"/Occupancy/NoiseOccupancyBarrel_detId");
74  MEP me_rollAreaBarrel_detId = dbe->get(subDir_+"/Occupancy/RollAreaBarrel_detId");
75  if ( me_noiseOccupancyBarrel_detId and me_rollAreaBarrel_detId )
76  {
77  TH1* h_noiseOccupancyBarrel_detId = me_noiseOccupancyBarrel_detId->getTH1();
78  TH1* h_rollAreaBarrel_detId = me_rollAreaBarrel_detId->getTH1();
79 
80  for ( int bin = 1, nBin = h_noiseOccupancyBarrel_detId->GetNbinsX(); bin <= nBin; ++bin )
81  {
82  const double noiseCount = h_noiseOccupancyBarrel_detId->GetBinContent(bin);
83  const double area = h_rollAreaBarrel_detId->GetBinContent(bin);
84  const double noiseLevel = area > 0 ? noiseCount/area/nEvent : 0;
85  if ( noiseLevel == 0. ) me_rollNoiseBarrel_noise->Fill(-maxNoise/50); // Fill underflow bin if noise is exactly zero
86  else me_rollNoiseBarrel_noise->Fill(std::min(noiseLevel, maxNoise));
87  }
88  }
89 
90  MEP me_noiseOccupancyEndcap_detId = dbe->get(subDir_+"/Occupancy/NoiseOccupancyEndcap_detId");
91  MEP me_rollAreaEndcap_detId = dbe->get(subDir_+"/Occupancy/RollAreaEndcap_detId");
92  if ( me_noiseOccupancyEndcap_detId and me_rollAreaEndcap_detId )
93  {
94  TH1* h_noiseOccupancyEndcap_detId = me_noiseOccupancyEndcap_detId->getTH1();
95  TH1* h_rollAreaEndcap_detId = me_rollAreaEndcap_detId->getTH1();
96 
97  for ( int bin = 1, nBin = h_noiseOccupancyEndcap_detId->GetNbinsX(); bin <= nBin; ++bin )
98  {
99  const double noiseCount = h_noiseOccupancyEndcap_detId->GetBinContent(bin);
100  const double area = h_rollAreaEndcap_detId->GetBinContent(bin);
101  const double noiseLevel = area > 0 ? noiseCount/area/nEvent : 0;
102  if ( noiseLevel == 0 ) me_rollNoiseEndcap_noise->Fill(-maxNoise/50); // Fill underflow bin if noise if exactly zero
103  else me_rollNoiseEndcap_noise->Fill(std::min(noiseLevel, maxNoise));
104  }
105  }
106 
107 }
int nEvent
Definition: myFastSimVal.cc:49
void Fill(long long x)
TH1 * getTH1(void) const
T min(T a, T b)
Definition: MathUtil.h:58

Member Data Documentation

std::string RPCRecHitValidClient::subDir_
private

Definition at line 30 of file RPCRecHitValidClient.h.