42 theSeedsCollectionLabel_ = consumes<TrajectorySeedCollection>(
parameters.getParameter<
InputTag>(
"SeedCollection"));
44 seedHitBin =
parameters.getParameter<
int>(
"RecHitBin");
45 seedHitMin =
parameters.getParameter<
double>(
"RecHitMin");
46 seedHitMax =
parameters.getParameter<
double>(
"RecHitMax");
47 PhiBin =
parameters.getParameter<
int>(
"PhiBin");
50 EtaBin =
parameters.getParameter<
int>(
"EtaBin");
53 ThetaBin =
parameters.getParameter<
int>(
"ThetaBin");
54 ThetaMin =
parameters.getParameter<
double>(
"ThetaMin");
55 ThetaMax =
parameters.getParameter<
double>(
"ThetaMax");
56 seedPtBin =
parameters.getParameter<
int>(
"seedPtBin");
57 seedPtMin =
parameters.getParameter<
double>(
"seedPtMin");
58 seedPtMax =
parameters.getParameter<
double>(
"seedPtMax");
59 seedPxyzBin =
parameters.getParameter<
int>(
"seedPxyzBin");
60 seedPxyzMin =
parameters.getParameter<
double>(
"seedPxyzMin");
61 seedPxyzMax =
parameters.getParameter<
double>(
"seedPxyzMax");
62 pErrBin =
parameters.getParameter<
int>(
"pErrBin");
63 pErrMin =
parameters.getParameter<
double>(
"pErrMin");
64 pErrMax =
parameters.getParameter<
double>(
"pErrMax");
65 pxyzErrBin =
parameters.getParameter<
int>(
"pxyzErrBin");
66 pxyzErrMin =
parameters.getParameter<
double>(
"pxyzErrMin");
67 pxyzErrMax =
parameters.getParameter<
double>(
"pxyzErrMax");
68 phiErrBin =
parameters.getParameter<
int>(
"phiErrBin");
69 phiErrMin =
parameters.getParameter<
double>(
"phiErrMin");
70 phiErrMax =
parameters.getParameter<
double>(
"phiErrMax");
71 etaErrBin =
parameters.getParameter<
int>(
"etaErrBin");
72 etaErrMin =
parameters.getParameter<
double>(
"etaErrMin");
73 etaErrMax =
parameters.getParameter<
double>(
"etaErrMax");
83 string histname =
"NumberOfRecHitsPerSeed_";
84 NumberOfRecHitsPerSeed = ibooker.
book1D(histname,
"Number of seed recHits", seedHitBin, seedHitMin, seedHitMax);
86 histname =
"seedPhi_";
90 histname =
"seedEta_";
93 histname =
"seedTheta_";
94 seedTheta = ibooker.
book1D(histname,
"Seed #theta", ThetaBin, ThetaMin, ThetaMax);
98 seedPt = ibooker.
book1D(histname,
"Seed p_{t}", seedPtBin, seedPtMin, seedPtMax);
101 histname =
"seedPx_";
102 seedPx = ibooker.
book1D(histname,
"Seed p_{x}", seedPxyzBin, seedPxyzMin, seedPxyzMax);
104 histname =
"seedPy_";
105 seedPy = ibooker.
book1D(histname,
"Seed p_{y}", seedPxyzBin, seedPxyzMin, seedPxyzMax);
107 histname =
"seedPz_";
108 seedPz = ibooker.
book1D(histname,
"Seed p_{z}", seedPxyzBin, seedPxyzMin, seedPxyzMax);
111 histname =
"seedPtErrOverPt_";
112 seedPtErr = ibooker.
book1D(histname,
"Seed p_{t}Err/p_{t}", pErrBin, pErrMin, pErrMax);
113 histname =
"seedPtErrOverPtVsPhi_";
114 seedPtErrVsPhi = ibooker.
book2D(histname,
"Seed p_{t}Err/p_{t} vs #phi", PhiBin,
PhiMin,
PhiMax, pErrBin, pErrMin, pErrMax);
116 histname =
"seedPtErrOverPtVsEta_";
117 seedPtErrVsEta = ibooker.
book2D(histname,
"Seed p_{t}Err/p_{t} vs #eta", EtaBin,
EtaMin,
EtaMax, pErrBin, pErrMin, pErrMax);
118 histname =
"seedPtErrOverPtVsPt_";
119 seedPtErrVsPt = ibooker.
book2D(histname,
"Seed p_{t}Err/p_{t} vs p_{t}", seedPtBin/5, seedPtMin, seedPtMax, pErrBin, pErrMin, pErrMax);
121 histname =
"seedPErrOverP_";
122 seedPErr = ibooker.
book1D(histname,
"Seed pErr/p", pErrBin, pErrMin, pErrMax);
123 histname =
"seedPErrOverPVsPhi_";
124 seedPErrVsPhi = ibooker.
book2D(histname,
"Seed pErr/p vs #phi", PhiBin,
PhiMin,
PhiMax, pErrBin, pErrMin, pErrMax);
126 histname =
"seedPErrOverPVsEta_";
127 seedPErrVsEta = ibooker.
book2D(histname,
"Seed pErr/p vs #eta", EtaBin,
EtaMin,
EtaMax, pErrBin, pErrMin, pErrMax);
128 histname =
"seedPErrOverPVsPt_";
129 seedPErrVsPt = ibooker.
book2D(histname,
"Seed pErr/p vs p_{t}", seedPtBin/5, seedPtMin, seedPtMax, pErrBin, pErrMin, pErrMax);
132 histname =
"seedPxErrOverPx_";
133 seedPxErr = ibooker.
book1D(histname,
"Seed p_{x}Err/p_{x}", pxyzErrBin, pxyzErrMin, pxyzErrMax);
134 histname =
"seedPyErrOverPy_";
135 seedPyErr = ibooker.
book1D(histname,
"Seed p_{y}Err/p_{y}", pxyzErrBin, pxyzErrMin, pxyzErrMax);
136 histname =
"seedPzErrOverPz_";
137 seedPzErr = ibooker.
book1D(histname,
"Seed p_{z}Err/p_{z}", pxyzErrBin, pxyzErrMin, pxyzErrMax);
140 histname =
"seedPhiErr_";
141 seedPhiErr = ibooker.
book1D(histname,
"Seed #phi error", phiErrBin, phiErrMin, phiErrMax);
143 histname =
"seedEtaErr_";
144 seedEtaErr = ibooker.
book1D(histname,
"Seed #eta error", etaErrBin, etaErrMin, etaErrMax);
149 theService->update(iSetup);
153 iEvent.
getByToken(theSeedsCollectionLabel_, seeds);
158 for(TrajectorySeedCollection::const_iterator
seed = seeds->begin();
seed != seeds->end(); ++
seed){
166 const GeomDet* gdet = theService->trackingGeometry()->idToDet( seedDetId );
175 NumberOfRecHitsPerSeed->Fill(
seed->nHits());
bool getByToken(EDGetToken token, Handle< PROD > &result) const
const std::string metname
const CurvilinearTrajectoryError & curvilinearError() const
Sin< T >::type sin(const T &t)
ROOT::Math::SMatrix< double, 6, 6, ROOT::Math::MatRepSym< double, 6 > > AlgebraicSymMatrix66
Geom::Phi< T > phi() const
const CartesianTrajectoryError cartesianError() const
const double EtaMax[kNumberCalorimeter]
Geom::Theta< T > theta() const
MonitorElement * book1D(Args &&...args)
Abs< T >::type abs(const T &t)
unsigned int detId() const
const double EtaMin[kNumberCalorimeter]
const AlgebraicSymMatrix66 & matrix() const
void setCurrentFolder(const std::string &fullpath)
virtual ~MuonSeedsAnalyzer()
Destructor.
MonitorElement * book2D(Args &&...args)
void analyze(const edm::Event &, const edm::EventSetup &) override
const AlgebraicSymMatrix55 & matrix() const
GlobalVector globalMomentum() const
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
Power< A, B >::type pow(const A &a, const B &b)
MuonSeedsAnalyzer(const edm::ParameterSet &)
Constructor.
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override