CMS 3D CMS Logo
CSCValidation.h
Go to the documentation of this file.
1 #ifndef RecoLocalMuon_CSCValidation_H
2 #define RecoLocalMuon_CSCValidation_H
3 
18 #include "FWCore/Framework/interface/ConsumesCollector.h"
19 #include "FWCore/Framework/interface/MakerMacros.h"
20 #include "FWCore/Framework/interface/Frameworkfwd.h"
21 #include "FWCore/Framework/interface/one/EDAnalyzer.h"
22 #include "FWCore/Framework/interface/Event.h"
23 #include "DataFormats/Common/interface/Handle.h"
24 #include "FWCore/ParameterSet/interface/ParameterSet.h"
25 #include "FWCore/Framework/interface/EventSetup.h"
26 #include "FWCore/MessageLogger/interface/MessageLogger.h"
27 #include "FWCore/Utilities/interface/ESGetToken.h"
28 #include "DataFormats/CSCDigi/interface/CSCWireDigiCollection.h"
29 #include "DataFormats/CSCDigi/interface/CSCStripDigiCollection.h"
30 #include "DataFormats/CSCDigi/interface/CSCComparatorDigiCollection.h"
31 #include "DataFormats/CSCDigi/interface/CSCALCTDigiCollection.h"
32 #include "DataFormats/CSCDigi/interface/CSCCLCTDigiCollection.h"
33 #include "DataFormats/CSCDigi/interface/CSCCorrelatedLCTDigiCollection.h"
34 
35 #include "EventFilter/CSCRawToDigi/interface/CSCDCCEventData.h"
36 #include "EventFilter/CSCRawToDigi/interface/CSCDCCExaminer.h"
37 #include "EventFilter/CSCRawToDigi/interface/CSCEventData.h"
38 #include "EventFilter/CSCRawToDigi/interface/CSCCFEBData.h"
39 #include "EventFilter/CSCRawToDigi/interface/CSCALCTHeader.h"
40 #include "EventFilter/CSCRawToDigi/interface/CSCAnodeData.h"
41 #include "EventFilter/CSCRawToDigi/interface/CSCComparatorData.h"
42 #include "EventFilter/CSCRawToDigi/interface/CSCDDUEventData.h"
43 #include "EventFilter/CSCRawToDigi/interface/CSCTMBData.h"
44 #include "EventFilter/CSCRawToDigi/interface/CSCTMBHeader.h"
45 #include "EventFilter/CSCRawToDigi/interface/CSCRPCData.h"
46 #include "EventFilter/CSCRawToDigi/interface/CSCDCCExaminer.h"
47 #include "EventFilter/CSCRawToDigi/interface/CSCDCCEventData.h"
48 #include "EventFilter/CSCRawToDigi/interface/CSCCFEBData.h"
49 #include "EventFilter/CSCRawToDigi/interface/CSCCFEBTimeSlice.h"
50 #include "EventFilter/CSCRawToDigi/interface/CSCMonitorInterface.h"
51 #include "CondFormats/CSCObjects/interface/CSCCrateMap.h"
52 #include "CondFormats/DataRecord/interface/CSCCrateMapRcd.h"
53 #include "FWCore/ServiceRegistry/interface/Service.h"
54 //FEDRawData
55 #include "DataFormats/FEDRawData/interface/FEDRawData.h"
56 #include "DataFormats/FEDRawData/interface/FEDNumbering.h"
57 #include "DataFormats/FEDRawData/interface/FEDRawDataCollection.h"
58 
59 #include "DataFormats/CSCRecHit/interface/CSCSegmentCollection.h"
60 #include "DataFormats/MuonDetId/interface/RPCDetId.h"
61 
62 #include "DataFormats/L1GlobalMuonTrigger/interface/L1MuGMTReadoutRecord.h"
63 #include "DataFormats/L1GlobalMuonTrigger/interface/L1MuGMTReadoutCollection.h"
64 #include "DataFormats/Common/interface/TriggerResults.h"
65 #include "FWCore/Common/interface/TriggerNames.h"
66 
67 #include "Geometry/CSCGeometry/interface/CSCGeometry.h"
68 #include "Geometry/Records/interface/MuonGeometryRecord.h"
69 
70 #include "DataFormats/TrackReco/interface/Track.h"
71 #include "DataFormats/TrackReco/interface/TrackFwd.h"
72 #include "Geometry/CommonDetUnit/interface/GlobalTrackingGeometry.h"
73 #include "Geometry/Records/interface/GlobalTrackingGeometryRecord.h"
74 
75 #include "DataFormats/CSCRecHit/interface/CSCRecHit2DCollection.h"
76 #include "SimDataFormats/TrackingHit/interface/PSimHitContainer.h"
77 
78 #include "DataFormats/GeometryVector/interface/GlobalPoint.h"
79 #include "DataFormats/GeometryVector/interface/GlobalVector.h"
80 #include "DataFormats/GeometryVector/interface/LocalPoint.h"
81 #include "DataFormats/GeometryVector/interface/LocalVector.h"
82 #include "DataFormats/CLHEP/interface/AlgebraicObjects.h"
83 #include "DataFormats/MuonDetId/interface/CSCIndexer.h"
84 
85 #include "CondFormats/CSCObjects/interface/CSCDBGains.h"
86 #include "CondFormats/DataRecord/interface/CSCDBGainsRcd.h"
87 #include "CondFormats/CSCObjects/interface/CSCDBNoiseMatrix.h"
88 #include "CondFormats/DataRecord/interface/CSCDBNoiseMatrixRcd.h"
89 #include "CondFormats/CSCObjects/interface/CSCDBCrosstalk.h"
90 #include "CondFormats/DataRecord/interface/CSCDBCrosstalkRcd.h"
91 #include "CondFormats/CSCObjects/interface/CSCDBPedestals.h"
92 #include "CondFormats/DataRecord/interface/CSCDBPedestalsRcd.h"
93 
94 #include "CommonTools/Statistics/interface/ChiSquaredProbability.h"
95 
96 #include "RecoLocalMuon/CSCValidation/src/CSCValHists.h"
97 #include "TVector3.h"
98 #include "TH1F.h"
99 #include "TH2F.h"
100 #include "TFile.h"
101 #include "TString.h"
102 #include "TTree.h"
103 #include "TProfile2D.h"
104 
105 namespace {
106  class CSCValidation : public edm::one::EDAnalyzer<> {
107  public:
109  CSCValidation(const edm::ParameterSet &pset);
110 
112  ~CSCValidation() override;
113 
115  void analyze(edm::Event const &event, edm::EventSetup const &eventSetup) override;
116  void beginJob() override;
117  void endJob() override;
118 
119  // for noise module
120  struct ltrh {
121  bool operator()(const CSCRecHit2D &rh1, const CSCRecHit2D &rh2) const {
122  return ((rh1.localPosition()).x() - (rh2.localPosition()).x()) < 0;
123  }
124  };
125 
126  protected:
127  private:
128  // these are the "modules"
129  // if you would like to add code to CSCValidation, please do so by adding an
130  // extra module in the form of an additional private member function
131  void doOccupancies(edm::Handle<CSCStripDigiCollection> strips,
132  edm::Handle<CSCWireDigiCollection> wires,
133  edm::Handle<CSCRecHit2DCollection> recHits,
134  edm::Handle<CSCSegmentCollection> cscSegments);
135  void doStripDigis(edm::Handle<CSCStripDigiCollection> strips);
136  void doWireDigis(edm::Handle<CSCWireDigiCollection> wires);
137  void doRecHits(edm::Handle<CSCRecHit2DCollection> recHits, edm::ESHandle<CSCGeometry> cscGeom);
138  void doSimHits(edm::Handle<CSCRecHit2DCollection> recHits, edm::Handle<edm::PSimHitContainer> simHits);
139  void doPedestalNoise(edm::Handle<CSCStripDigiCollection> strips);
140  void doSegments(edm::Handle<CSCSegmentCollection> cscSegments, edm::ESHandle<CSCGeometry> cscGeom);
141  void doResolution(edm::Handle<CSCSegmentCollection> cscSegments, edm::ESHandle<CSCGeometry> cscGeom);
142  void doEfficiencies(edm::Handle<CSCWireDigiCollection> wires,
143  edm::Handle<CSCStripDigiCollection> strips,
144  edm::Handle<CSCRecHit2DCollection> recHits,
145  edm::Handle<CSCSegmentCollection> cscSegments,
146  edm::ESHandle<CSCGeometry> cscGeom);
147  void doGasGain(const CSCWireDigiCollection &, const CSCStripDigiCollection &, const CSCRecHit2DCollection &);
148  void doCalibrations(const edm::EventSetup &eventSetup);
149  void doAFEBTiming(const CSCWireDigiCollection &);
150  void doCompTiming(const CSCComparatorDigiCollection &);
151  void doADCTiming(const CSCRecHit2DCollection &);
152  void doNoiseHits(edm::Handle<CSCRecHit2DCollection> recHits,
153  edm::Handle<CSCSegmentCollection> cscSegments,
154  edm::ESHandle<CSCGeometry> cscGeom,
155  edm::Handle<CSCStripDigiCollection> strips);
156  bool doTrigger(edm::Handle<L1MuGMTReadoutCollection> pCollection);
157  void doStandalone(edm::Handle<reco::TrackCollection> saMuons);
158  void doTimeMonitoring(edm::Handle<CSCRecHit2DCollection> recHits,
159  edm::Handle<CSCSegmentCollection> cscSegments,
160  edm::Handle<CSCALCTDigiCollection> alcts,
161  edm::Handle<CSCCLCTDigiCollection> clcts,
162  edm::Handle<CSCCorrelatedLCTDigiCollection> correlatedlcts,
163  edm::Handle<L1MuGMTReadoutCollection> pCollection,
164  edm::ESHandle<CSCGeometry> cscGeom,
165  const edm::EventSetup &eventSetup,
166  const edm::Event &event);
167  bool doHLT(edm::Handle<edm::TriggerResults> hltResults);
168 
169  // some useful functions
170  bool filterEvents(edm::Handle<CSCRecHit2DCollection> recHits,
171  edm::Handle<CSCSegmentCollection> cscSegments,
172  edm::Handle<reco::TrackCollection> saMuons);
173  float fitX(const CLHEP::HepMatrix &sp, const CLHEP::HepMatrix &ep);
174  float getSignal(const CSCStripDigiCollection &stripdigis, CSCDetId idRH, int centerStrip);
175  float getthisSignal(const CSCStripDigiCollection &stripdigis, CSCDetId idRH, int centerStrip);
176  int getWidth(const CSCStripDigiCollection &stripdigis, CSCDetId idRH, int centerStrip);
177  void findNonAssociatedRecHits(edm::ESHandle<CSCGeometry> cscGeom, edm::Handle<CSCStripDigiCollection> strips);
178  int chamberSerial(CSCDetId id);
179  int ringSerial(CSCDetId id);
180 
181  // these functions handle Stoyan's efficiency code
182  void fillEfficiencyHistos(int bin, int flag);
183  void getEfficiency(float bin, float Norm, std::vector<float> &eff);
184  void histoEfficiency(TH1F *readHisto, TH1F *writeHisto);
185  double lineParametrization(double z1Position, double z2Position, double z1Direction) {
186  double parameterLine = (z2Position - z1Position) / z1Direction;
187  return parameterLine;
188  }
189  double extrapolate1D(double initPosition, double initDirection, double parameterOfTheLine) {
190  double extrapolatedPosition = initPosition + initDirection * parameterOfTheLine;
191  return extrapolatedPosition;
192  }
193  bool withinSensitiveRegion(LocalPoint localPos,
194  const std::array<const float, 4> &layerBounds,
195  int station,
196  int ring,
197  float shiftFromEdge,
198  float shiftFromDeadZone);
199 
200  // counters
201  int nEventsAnalyzed;
202  int rhTreeCount;
203  int segTreeCount;
204  bool firstEvent;
205  bool cleanEvent;
206 
207  //
208  //
209  // The root file for the histograms.
210  TFile *theFile;
211 
212  //
213  //
214  // input parameters for this module
215  bool makePlots;
216  bool makeComparisonPlots;
217  std::string refRootFile;
218  bool writeTreeToFile;
219  bool isSimulation;
220  std::string rootFileName;
221  bool detailedAnalysis;
222  bool useDigis;
223  bool makeHLTPlots;
224 
225  // filters
226  bool useQualityFilter;
227  bool useTriggerFilter;
228 
229  // quality filter parameters
230  double pMin;
231  double chisqMax;
232  int nCSCHitsMin, nCSCHitsMax;
233  double lengthMin, lengthMax;
234  double deltaPhiMax;
235  double polarMin, polarMax;
236 
237  edm::EDGetTokenT<FEDRawDataCollection> rd_token;
238  edm::EDGetTokenT<CSCWireDigiCollection> wd_token;
239  edm::EDGetTokenT<CSCStripDigiCollection> sd_token;
240  edm::EDGetTokenT<CSCComparatorDigiCollection> cd_token;
241  edm::EDGetTokenT<CSCALCTDigiCollection> al_token;
242  edm::EDGetTokenT<CSCCLCTDigiCollection> cl_token;
243  edm::EDGetTokenT<CSCCorrelatedLCTDigiCollection> co_token;
244  edm::EDGetTokenT<CSCRecHit2DCollection> rh_token;
245  edm::EDGetTokenT<CSCSegmentCollection> se_token;
246  edm::EDGetTokenT<L1MuGMTReadoutCollection> l1_token;
247  edm::EDGetTokenT<edm::TriggerResults> tr_token;
248  edm::EDGetTokenT<reco::TrackCollection> sa_token;
249  edm::EDGetTokenT<edm::PSimHitContainer> sh_token;
250  // geometry
251  edm::ESGetToken<CSCGeometry, MuonGeometryRecord> geomToken_;
252  edm::ESGetToken<CSCCrateMap, CSCCrateMapRcd> crateToken_;
253  // conditions data
254  edm::ESGetToken<CSCDBPedestals, CSCDBPedestalsRcd> pedestalsToken_;
255  edm::ESGetToken<CSCDBGains, CSCDBGainsRcd> gainsToken_;
256  edm::ESGetToken<CSCDBNoiseMatrix, CSCDBNoiseMatrixRcd> noiseToken_;
257  edm::ESGetToken<CSCDBCrosstalk, CSCDBCrosstalkRcd> crosstalkToken_;
258 
259  // module on/off switches
260  bool makeOccupancyPlots;
261  bool makeTriggerPlots;
262  bool makeStripPlots;
263  bool makeWirePlots;
264  bool makeRecHitPlots;
265  bool makeSimHitPlots;
266  bool makeSegmentPlots;
267  bool makeResolutionPlots;
268  bool makePedNoisePlots;
269  bool makeEfficiencyPlots;
270  bool makeGasGainPlots;
271  bool makeAFEBTimingPlots;
272  bool makeCompTimingPlots;
273  bool makeADCTimingPlots;
274  bool makeRHNoisePlots;
275  bool makeCalibPlots;
276  bool makeStandalonePlots;
277  bool makeTimeMonitorPlots;
278 
279  // The histo managing object
280  CSCValHists *histos;
281 
282  // tmp histos for Efficiency
283  TH1F *hSSTE;
284  TH1F *hRHSTE;
285  TH1F *hSEff;
286  TH1F *hRHEff;
287  TH2F *hSSTE2;
288  TH2F *hRHSTE2;
289  TH2F *hStripSTE2;
290  TH2F *hWireSTE2;
291  TH2F *hSEff2;
292  TH2F *hRHEff2;
293  TH2F *hStripEff2;
294  TH2F *hWireEff2;
295  TH2F *hEffDenominator;
296  TH2F *hSensitiveAreaEvt;
297 
298  TH1F *hSSTETight;
299  TH1F *hRHSTETight;
300  TH2F *hSSTE2Tight;
301  TH2F *hRHSTE2Tight;
302  TH2F *hStripSTE2Tight;
303  TH2F *hWireSTE2Tight;
304  TH2F *hEffDenominatorTight;
305 
306  // occupancy
307  TH2I *hOWires;
308  TH2I *hOStrips;
309  TH2I *hORecHits;
310  TH2I *hOSegments;
311 
313  std::vector<int> nmbhvsegm;
314  std::map<int, std::vector<int> > m_wire_hvsegm;
315  std::map<int, int> m_single_wire_layer;
316 
317  //maps to store the DetId and associated RecHits
318  std::multimap<CSCDetId, CSCRecHit2D> AllRechits;
319  std::multimap<CSCDetId, CSCRecHit2D> SegRechits;
320  std::multimap<CSCDetId, CSCRecHit2D> NonAssociatedRechits;
321  std::map<CSCRecHit2D, float, ltrh> distRHmap;
322 
323  int typeIndex(CSCDetId id) {
324  // linearlized index bases on endcap, station, and ring
325  int index = 0;
326  if (id.station() == 1) {
327  index = id.ring() + 1;
328  if (id.ring() == 4)
329  index = 1;
330  } else
331  index = id.station() * 2 + id.ring();
332  if (id.endcap() == 1)
333  index = index + 10;
334  if (id.endcap() == 2)
335  index = 11 - index;
336  return index;
337  }
338  };
339 } // namespace
340 #endif
SiStripBadComponentsDQMServiceTemplate_cfg.ep
ep
Definition: SiStripBadComponentsDQMServiceTemplate_cfg.py:86
FastTrackerRecHitMaskProducer_cfi.recHits
recHits
Definition: FastTrackerRecHitMaskProducer_cfi.py:8
muonDTDigis_cfi.pset
pset
Definition: muonDTDigis_cfi.py:27
dtChamberEfficiency_cfi.cscSegments
cscSegments
Definition: dtChamberEfficiency_cfi.py:15
dtNoiseAnalysis_cfi.detailedAnalysis
detailedAnalysis
Definition: dtNoiseAnalysis_cfi.py:8
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
edm::one::EDAnalyzerBase::analyze
virtual void analyze(Event const &, EventSetup const &)=0
amptDefault_cfi.firstEvent
firstEvent
Definition: amptDefault_cfi.py:8
ALCARECOTkAlMinBias_cff.pMin
pMin
GeV.
Definition: ALCARECOTkAlMinBias_cff.py:35
CSCRecHit2D::localPosition
LocalPoint localPosition() const override
Definition: CSCRecHit2D.h:56
portablecollection::typeIndex
constexpr std::size_t typeIndex
Definition: PortableCollectionCommon.h:99
ecalTB2006H4_GenSimDigiReco_cfg.doSimHits
doSimHits
Definition: ecalTB2006H4_GenSimDigiReco_cfg.py:252
combine.histos
histos
Definition: combine.py:4
JetHT_cfg.hltResults
hltResults
Definition: JetHT_cfg.py:259
DigiDM_cff.strips
strips
#turn off noise in all subdetectors simHcalUnsuppressedDigis.doNoise = False mix.digitizers.hcal.doNoise = False simEcalUnsuppressedDigis.doNoise = False mix.digitizers.ecal.doNoise = False simEcalUnsuppressedDigis.doESNoise = False simSiPixelDigis.AddNoise = False mix.digitizers.pixel.AddNoise = False simSiStripDigis.Noise = False mix.digitizers.strip.AddNoise = False
Definition: DigiDM_cff.py:32
DigiDM_cff.wires
wires
Definition: DigiDM_cff.py:33
options_cfi.eventSetup
eventSetup
Definition: options_cfi.py:11
makePlots
void makePlots(Int_t nFiles, TString *files, TString *names, TString misalignment, Double_t *values, Double_t *phases, TString directory, Bool_t matrix[xsize][ysize])
FastTrackerRecHitCombiner_cfi.simHits
simHits
Definition: FastTrackerRecHitCombiner_cfi.py:5
event
Definition: event.py:1
CSCSkim_cfi.rootFileName
rootFileName
Definition: CSCSkim_cfi.py:9
Generated for CMSSW Reference Manual by doxygen 1.8.14