|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
|
#ifndef DPGAnalysisSkims_CSCSkim_H
#define DPGAnalysisSkims_CSCSkim_H
/** \class CSCSkim
*
* This simple program selects minimal CSC events for output.
*
* Michael Schmitt, Northwestern University, July 2008
*/
#include <memory>
#include <iostream>
#include <vector>
#include <map>
#include <string>
#include <iomanip>
#include <fstream>
#include "FWCore/Framework/interface/ConsumesCollector.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/one/EDFilter.h"
#include "FWCore/Framework/interface/Event.h"
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "DataFormats/MuonDetId/interface/CSCDetId.h"
#include "DataFormats/CSCDigi/interface/CSCWireDigi.h"
#include "DataFormats/CSCDigi/interface/CSCWireDigiCollection.h"
#include "DataFormats/CSCDigi/interface/CSCStripDigi.h"
#include "DataFormats/CSCDigi/interface/CSCStripDigiCollection.h"
#include <DataFormats/CSCRecHit/interface/CSCRecHit2D.h>
#include "DataFormats/CSCRecHit/interface/CSCRecHit2DCollection.h"
#include <DataFormats/CSCRecHit/interface/CSCSegmentCollection.h>
#include "DataFormats/DTRecHit/interface/DTRecSegment4D.h"
#include "DataFormats/DTRecHit/interface/DTRecSegment4DCollection.h"
#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackReco/interface/TrackFwd.h"
#include "Geometry/Records/interface/GlobalTrackingGeometryRecord.h"
#include "Geometry/CommonDetUnit/interface/GlobalTrackingGeometry.h"
#include "Geometry/DTGeometry/interface/DTGeometry.h"
#include "Geometry/CSCGeometry/interface/CSCGeometry.h"
#include "Geometry/CSCGeometry/interface/CSCChamber.h"
#include "Geometry/CSCGeometry/interface/CSCLayer.h"
#include "Geometry/CSCGeometry/interface/CSCLayerGeometry.h"
#include "Geometry/Records/interface/MuonGeometryRecord.h"
#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "DataFormats/GeometryVector/interface/GlobalVector.h"
#include "DataFormats/GeometryVector/interface/LocalPoint.h"
#include "DataFormats/GeometryVector/interface/LocalVector.h"
#include "DataFormats/CLHEP/interface/AlgebraicObjects.h"
// #include "SimDataFormats/TrackingHit/interface/PSimHitContainer.h"
#include "TVector3.h"
#include "TH1F.h"
#include "TH2F.h"
#include "TFile.h"
#include "TString.h"
#include "TTree.h"
class CSCSkim : public edm::one::EDFilter<> {
public:
// Constructor
explicit CSCSkim(const edm::ParameterSet &pset);
// Destructor
~CSCSkim() override;
// Analysis routines
void beginJob() override;
bool filter(edm::Event &event, const edm::EventSetup &eventSetup) override;
void endJob() override;
protected:
private:
// main skimming routine
bool doCSCSkimming(edm::Handle<CSCRecHit2DCollection> cscRecHits, edm::Handle<CSCSegmentCollection> cscSegments);
// extra skimming routine for alignment studies
bool doOverlapSkimming(edm::Handle<CSCSegmentCollection> cscSegments);
// skimming routine for messy events
bool doMessyEventSkimming(edm::Handle<CSCRecHit2DCollection> cscRecHits,
edm::Handle<CSCSegmentCollection> cscSegments);
// select events with DIGIs in a certain chamber
bool doCertainChamberSelection(edm::Handle<CSCWireDigiCollection> wires, edm::Handle<CSCStripDigiCollection> strips);
// select events which might probe the DT-CSC overlap region
bool doDTOverlap(edm::Handle<CSCSegmentCollection> cscSegments);
// select muons which go through inner part of stations 1,2,3,4
bool doHaloLike(edm::Handle<CSCSegmentCollection> cscSegments);
// select events with long stand-alone (cosmic) muons
bool doLongSATrack(edm::Handle<reco::TrackCollection> saTracks);
// select events suitable for magnetic field studies - they have a track in the tracker
bool doBFieldStudySelection(edm::Handle<reco::TrackCollection> saTracks,
edm::Handle<reco::TrackCollection> Tracks,
edm::Handle<reco::MuonCollection> gMuons);
// some useful functions
int chamberSerial(int kE, int kS, int kR, int kCh);
// counters
int nEventsAnalyzed;
int nEventsSelected;
int nEventsChambersBothSides;
int nEventsOverlappingChambers;
int nEventsMessy;
int nEventsCertainChamber;
int nEventsDTOverlap;
int nEventsHaloLike;
int nEventsLongSATrack;
int nEventsForBFieldStudies;
// run and event number
int iRun;
int iEvent;
// The root file for the histograms.
TFile *theHistogramFile;
// file names
std::string outputFileName;
std::string histogramFileName;
// es token names
const edm::ESGetToken<CSCGeometry, MuonGeometryRecord> m_CSCGeomToken;
// token names
edm::EDGetTokenT<CSCWireDigiCollection> wds_token;
edm::EDGetTokenT<CSCStripDigiCollection> sds_token;
edm::EDGetTokenT<CSCWireDigiCollection> wdr_token;
edm::EDGetTokenT<CSCStripDigiCollection> sdr_token;
edm::EDGetTokenT<CSCRecHit2DCollection> rh_token;
edm::EDGetTokenT<CSCSegmentCollection> seg_token;
edm::EDGetTokenT<reco::TrackCollection> sam_token;
edm::EDGetTokenT<reco::TrackCollection> trk_token;
edm::EDGetTokenT<reco::MuonCollection> glm_token;
// parameters for the selection
bool isSimulation;
int typeOfSkim;
int nLayersWithHitsMinimum;
int minimumHitChambers;
int minimumSegments;
bool demandChambersBothSides;
bool makeHistograms;
bool makeHistogramsForMessyEvents;
int whichEndcap;
int whichStation;
int whichRing;
int whichChamber;
// parameters for B-field study skim
float pMin;
float zLengthMin;
int nCSCHitsMin;
float zInnerMax;
int nTrHitsMin;
float zLengthTrMin;
float rExtMax;
float redChiSqMax;
int nValidHitsMin;
// histograms for skimming module
TH1F *hxnRecHits;
TH1F *hxnSegments;
TH1F *hxnHitChambers;
TH1F *hxnRecHitsSel;
TH1F *mevnRecHits0;
TH1F *mevnChambers0;
TH1F *mevnSegments0;
TH1F *mevnRecHits1;
TH1F *mevnChambers1;
TH1F *mevnSegments1;
TH1F *xxP, *xxnValidHits, *xxnCSCHits, *xxredChiSq, *xxnTrackerHits;
};
#endif
|