25 , doTrackerSpecific_ ( conf_.getParameter<bool>(
"doTrackerSpecific") )
26 , doAllPlots_ ( conf_.getParameter<bool>(
"doAllPlots") )
27 , doBSPlots_ ( conf_.getParameter<bool>(
"doBeamSpotPlots") )
28 , doPVPlots_ ( conf_.getParameter<bool>(
"doPrimaryVertexPlots") )
29 , doDCAPlots_ ( conf_.getParameter<bool>(
"doDCAPlots") )
30 , doGeneralPropertiesPlots_ ( conf_.getParameter<bool>(
"doGeneralPropertiesPlots") )
31 , doMeasurementStatePlots_ ( conf_.getParameter<bool>(
"doMeasurementStatePlots") )
32 , doHitPropertiesPlots_ ( conf_.getParameter<bool>(
"doHitPropertiesPlots") )
33 , doRecHitVsPhiVsEtaPerTrack_ ( conf_.getParameter<bool>(
"doRecHitVsPhiVsEtaPerTrack") )
34 , doLayersVsPhiVsEtaPerTrack_ ( conf_.getParameter<bool>(
"doLayersVsPhiVsEtaPerTrack") )
35 , doRecHitsPerTrackProfile_ ( conf_.getParameter<bool>(
"doRecHitsPerTrackProfile") )
36 , doThetaPlots_ ( conf_.getParameter<bool>(
"doThetaPlots") )
37 , doTrackPxPyPlots_ ( conf_.getParameter<bool>(
"doTrackPxPyPlots") )
38 , doDCAwrtPVPlots_ ( conf_.getParameter<bool>(
"doDCAwrtPVPlots") )
39 , doDCAwrt000Plots_ ( conf_.getParameter<bool>(
"doDCAwrt000Plots") )
40 , doLumiAnalysis_ ( conf_.getParameter<bool>(
"doLumiAnalysis") )
41 , doTestPlots_ ( conf_.getParameter<bool>(
"doTestPlots") )
42 , doHIPlots_ ( conf_.getParameter<bool>(
"doHIPlots") )
43 , doSIPPlots_ ( conf_.getParameter<bool>(
"doSIPPlots") )
44 , doEffFromHitPatternVsPU_ ( conf_.getParameter<bool>(
"doEffFromHitPatternVsPU") )
45 , doEffFromHitPatternVsBX_ ( conf_.getParameter<bool>(
"doEffFromHitPatternVsBX") )
46 , pvNDOF_ ( conf_.getParameter<int> (
"pvNDOF") )
47 , qualityString_ ( conf_.getParameter<std::
string>(
"qualityString"))
163 if (StateName ==
"All") {
168 StateName !=
"OuterSurface" &&
169 StateName !=
"InnerSurface" &&
170 StateName !=
"ImpactPoint" &&
171 StateName !=
"default"
190 int NBINS[] = { 50, 3564 };
191 float MAX[] = { 50.5, 3564.5 };
196 for (
int i=0;
i<monQuantity::END;
i++) {
209 const char * dets[] = {
"None",
"PXB",
"PXF",
"TIB",
"TID",
"TOB",
"TEC"};
215 const char * hit_category[] = {
"valid",
"missing",
"inactive",
"bad",
"total"};
219 for (
unsigned int det = 1; det <
sizeof(dets)/
sizeof(
char*); ++det ) {
220 for (
unsigned int sub_det = 1;
221 sub_det <= trackerGeometry->numberOfLayers(det); ++sub_det) {
222 for (
unsigned int cat = 0;
223 cat <
sizeof(hit_category)/
sizeof(
char *); ++
cat) {
224 memset(title, 0,
sizeof(title));
225 snprintf(title,
sizeof(title),
"Hits%s_%s_%s_Subdet%d", suffix.c_str(), hit_category[
cat], dets[det], sub_det);
229 Key(det, sub_det, mon),
230 ibooker.
book1D(title, title, nbins, 0.5, max)));
234 Key(det, sub_det, mon),
235 ibooker.
book1D(title, title, nbins, 0.5, max)));
239 Key(det, sub_det, mon),
240 ibooker.
book1D(title, title, nbins, 0.5, max)));
244 Key(det, sub_det, mon),
245 ibooker.
book1D(title, title, nbins, 0.5, max)));
249 Key(det, sub_det, mon),
250 ibooker.
book1D(title, title, nbins, 0.5, max)));
253 LogDebug(
"TrackAnalyzer") <<
"Invalid hit category used " <<
cat <<
" ignored\n";
315 histname =
"NumberOfRecHitsPerTrack_";
320 histname =
"NumberOfValidRecHitsPerTrack_";
325 histname =
"NumberOfLostRecHitsPerTrack_";
330 histname =
"NumberOfMissingInnerRecHitsPerTrack_";
335 histname =
"NumberOfMissingOuterRecHitsPerTrack_";
340 histname =
"ValidFractionPerTrack_";
349 histname =
"NumberOfValidRecHitVsPhiVsEtaPerTrack_";
351 EtaBin, EtaMin, EtaMax, PhiBin, PhiMin, PhiMax, 0, 40.,
"");
355 histname =
"NumberOfLostRecHitVsPhiVsEtaPerTrack_";
357 EtaBin, EtaMin, EtaMax, PhiBin, PhiMin, PhiMax, 0, 5.,
"");
362 histname =
"NumberMIRecHitVsPhiVsEtaPerTrack_";
364 EtaBin, EtaMin, EtaMax, PhiBin, PhiMin, PhiMax, 0, 15.,
"");
368 histname =
"NumberMORecHitVsPhiVsEtaPerTrack_";
370 EtaBin, EtaMin, EtaMax, PhiBin, PhiMin, PhiMax, 0, 15.,
"");
374 histname =
"ValidFractionVsPhiVsEtaPerTrack_";
376 EtaBin, EtaMin, EtaMax, PhiBin, PhiMin, PhiMax, 0, 2.,
"");
383 std::string layerTypeName[4] = {
"",
"Off",
"3D",
"Missing"};
384 for (
int i=0;
i<4; ++
i) {
385 histname =
"NumberOf"+ layerTypeName[
i] +
"LayersPerTrack_";
391 for (
int i=0;
i<4; ++
i) {
392 histname =
"NumberOf"+ layerTypeName[
i] +
"LayersVsPhiVsEtaPerTrack_";
394 EtaBin, EtaMin, EtaMax, PhiBin, PhiMin, PhiMax, 0, 40.,
"");
433 Chi2->setAxisTitle(
"Track #chi^{2}" ,1);
434 Chi2->setAxisTitle(
"Number of Tracks",2);
497 histname =
"xPointOfClosestApproach_";
502 histname =
"yPointOfClosestApproach_";
507 histname =
"zPointOfClosestApproach_";
512 histname =
"xPointOfClosestApproachToPV_";
517 histname =
"yPointOfClosestApproachToPV_";
522 histname =
"zPointOfClosestApproachToPV_";
546 std::vector<std::string> StopReasonName = {
"UNINITIALIZED",
"MAX_HITS",
"MAX_LOST_HITS",
"MAX_CONSECUTIVE_LOST_HITS",
"LOST_HIT_FRACTION",
"MIN_PT",
"CHARGE_SIGNIFICANCE",
"LOOPER",
"MAX_CCC_LOST_HITS",
"NO_SEGMENTS_FOR_VALID_LAYERS",
"NOT_STOPPED" };
563 for (
size_t ibin=0; ibin<StopReasonName.size(); ibin++) {
599 histname =
"NumberOfRecHitsPerTrack_lumiFlag_";
653 histname =
"DistanceOfClosestApproachToBS_";
658 histname =
"DistanceOfClosestApproachToBSVsPhi_";
664 histname =
"xPointOfClosestApproachVsZ0wrt000_";
669 histname =
"yPointOfClosestApproachVsZ0wrt000_";
674 histname =
"xPointOfClosestApproachVsZ0wrtBS_";
679 histname =
"yPointOfClosestApproachVsZ0wrtBS_";
684 histname =
"zPointOfClosestApproachVsPhi_";
714 histname =
"DistanceOfClosestApproachToPV_";
719 histname =
"DistanceOfClosestApproachToPVVsPhi_";
725 histname =
"xPointOfClosestApproachVsZ0wrtPV_";
730 histname =
"yPointOfClosestApproachVsZ0wrtPV_";
748 histname =
"TESTDistanceOfClosestApproachToBS_";
753 histname =
"TESTDistanceOfClosestApproachToBSVsPhi_";
787 histname =
"DistanceOfClosestApproachVsTheta_";
793 histname =
"DistanceOfClosestApproachVsEta_";
799 histname =
"DistanceOfClosestApproach_";
804 histname =
"DistanceOfClosestApproachVsPhi_";
814 const double sipBins = 200;
815 const double sipMin = -20;
816 const double sipMax = 20;
862 if (recoPrimaryVerticesHandle.
isValid())
863 if (recoPrimaryVerticesHandle->size() > 0)
864 for (
auto v : *recoPrimaryVerticesHandle)
889 double chi2prob = TMath::Prob(track.
chi2(),(int)track.
ndof());
956 double stop = track.
stopReason() > max ? double(max-1) :
static_cast<double>(track.
stopReason());
994 if (recoPrimaryVerticesHandle.
isValid() && recoPrimaryVerticesHandle->size() > 0) {
1004 double longDCAsig = 0, transDCAsig = 0;
1007 if(zerr2 > 0) longDCAsig = track.
dz(pv.
position())/zerr2;
1008 if(xyerr2 > 0) transDCAsig = track.
dxy(pv.
position())/xyerr2;
1045 if(ip3d.first)
sip3dToPV->
Fill(ip3d.second.value() / ip3d.second.error());
1046 if(ip2d.first)
sip2dToPV->
Fill(ip2d.second.value() / ip2d.second.error());
1071 if (StateName ==
"All") {
1076 StateName !=
"OuterSurface" &&
1077 StateName !=
"InnerSurface" &&
1078 StateName !=
"ImpactPoint" &&
1079 StateName !=
"default"
1097 for (
int i=0;
i<monQuantity::END;
i++) {
1101 if (track.
pt() > 1.0 && track.
dxy() < 0.1 and monitoring > 0) {
1110 bool valid =
hp.validHitFilter(
pattern);
1112 bool inactive =
hp.inactiveHitFilter(
pattern);
1114 hit_type = valid ? 0 :
1116 ( inactive ? 2 : 3));
1118 LogDebug(
"TrackAnalyzer") <<
"Invalid combination of detector and subdetector: ("
1121 <<
"): ignoring it.\n";
1140 LogDebug(
"TrackAnalyzer") <<
"Invalid hit category used " << hit_type <<
" ignored\n";
1246 std::string histTag = (sname ==
"default") ? CategoryName : sname +
"_" + CategoryName;
1254 histname =
"Chi2oNDFVsTheta_" + histTag;
1259 histname =
"Chi2oNDFVsPhi_" + histTag;
1264 histname =
"Chi2oNDFVsEta_" + histTag;
1269 histname =
"Chi2ProbVsPhi_" + histTag;
1274 histname =
"Chi2ProbVsEta_" + histTag;
1320 histname =
"TrackEtaPhi_" + histTag;
1325 histname =
"TrackEtaPhiInner_" + histTag;
1330 histname =
"TrackEtaPhiOuter_" + histTag;
1338 histname =
"TrackTheta_" + histTag;
1348 histname =
"TrackPErrOverP_" + histTag;
1353 histname =
"TrackPtErrOverPt_" + histTag;
1358 histname =
"TrackPtErrOverPtVsEta_" + histTag;
1364 histname =
"TrackPxErrOverPx_" + histTag;
1369 histname =
"TrackPyErrOverPy_" + histTag;
1374 histname =
"TrackPzErrOverPz_" + histTag;
1379 histname =
"TrackPhiErr_" + histTag;
1384 histname =
"TrackEtaErr_" + histTag;
1391 histname =
"NumberOfRecHitsPerTrackVsPhi_" + histTag;
1397 histname =
"NumberOfRecHitsPerTrackVsTheta_" + histTag;
1402 histname =
"NumberOfRecHitsPerTrackVsEta_" + histTag;
1407 histname =
"NumberOfValidRecHitsPerTrackVsPhi_" + histTag;
1413 histname =
"NumberOfValidRecHitsPerTrackVsEta_" + histTag;
1419 histname =
"NumberOfLayersPerTrackVsPhi_" + histTag;
1425 histname =
"NumberOfLayersPerTrackVsTheta_" + histTag;
1430 histname =
"NumberOfLayersPerTrackVsEta_" + histTag;
1436 histname =
"Chi2oNDFVsTheta_" + histTag;
1442 histname =
"Chi2oNDFVsPhi_" + histTag;
1447 histname =
"Chi2oNDFVsEta_" + histTag;
1452 histname =
"Chi2ProbVsPhi_" + histTag;
1457 histname =
"Chi2ProbVsEta_" + histTag;
1475 double pxerror, pyerror, pzerror, pterror, perror, phierror,
etaerror;
1483 if (sname ==
"default") {
1491 theta = track.
theta();
1538 std::map<std::string, TkParameterMEs>::iterator iPos =
TkParameterMEMap.find(sname);
1596 double chi2prob = TMath::Prob(track.
chi2(),(int)track.
ndof());
1643 for (
auto det : subdetectors ) {
1660 histname =
"NumberOfRecHitsPerTrack_" + det +
"_" + CategoryName;
1665 histname =
"NumberOfRecHitsPerTrackVsPhi_" + det +
"_" + CategoryName;
1670 histname =
"NumberOfRecHitsPerTrackVsEta_" + det +
"_" + CategoryName;
1675 histname =
"NumberOfLayersPerTrack_" + det +
"_" + CategoryName;
1680 histname =
"NumberOfLayersPerTrackVsPhi_" + det +
"_" + CategoryName;
1685 histname =
"NumberOfLayersPerTrackVsEta_" + det +
"_" + CategoryName;
1702 double phi = track.
phi();
1703 double eta = track.
eta();
1708 int nValidLayers = 0;
1709 int nValidRecHits = 0;
1710 int substr = it->second.detectorId;
1741 it->second.NumberOfRecHitsPerTrack ->
Fill(nValidRecHits);
1742 it->second.NumberOfRecHitsPerTrackVsPhi ->
Fill(phi, nValidRecHits);
1743 it->second.NumberOfRecHitsPerTrackVsEta ->
Fill(eta, nValidRecHits);
1745 it->second.NumberOfLayersPerTrack ->
Fill(nValidLayers);
1746 it->second.NumberOfLayersPerTrackVsPhi ->
Fill(phi, nValidLayers);
1747 it->second.NumberOfLayersPerTrackVsEta ->
Fill(eta, nValidLayers);
MonitorElement * NumberOfRecHitsPerTrackVsPhi
MonitorElement * stoppingSourceVSphi
MonitorElement * NumberOfValidRecHitsPerTrackVsPhi
MonitorElement * Chi2ProbVsPhi
double p() const
momentum vector magnitude
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
MonitorElement * NumberOfLayersPerTrackVsPhi
int stripTOBLayersWithMeasurement() const
T getParameter(std::string const &) const
MonitorElement * TrackEtaPhi
MonitorElement * DistanceOfClosestApproachToPVVsPhi
const Point & referencePoint() const
Reference point on the track.
MonitorElement * NumberOfRecHitVsPhiVsEtaPerTrack
MonitorElement * ValidFractionVsPhiVsEtaPerTrack
MonitorElement * TrackPErr
double d0Error() const
error on d0
MonitorElement * dNhitdPt_HighPurity
void fillHistosForTrackerSpecific(const reco::Track &track)
MonitorElement * NumberOfRecHitsPerTrackVsEta
std::unordered_map< Key, MonitorElement *, KeyHasher > hits_valid_
void bookHistosForLScertification(DQMStore::IBooker &ibooker)
MonitorElement * stoppingSourceVSeta
TrackCharge charge() const
double validFraction() const
fraction of valid hits on the track
double d0() const
dxy parameter in perigee convention (d0 = -dxy)
std::unordered_map< Key, MonitorElement *, KeyHasher > hits_bad_
MonitorElement * NumberOfValidRecHitsPerTrack
int stripTIBLayersWithMeasurement() const
double normalizedChi2() const
chi-squared divided by n.d.o.f. (or chi-squared * 1e6 if n.d.o.f. is zero)
MonitorElement * bookProfile(Args &&...args)
bool getByToken(EDGetToken token, Handle< PROD > &result) const
double zError() const
error on z
MonitorElement * Chi2oNDF_lumiFlag
double theta() const
polar angle
double dxyError() const
error on dxy
MonitorElement * xPointOfClosestApproachToPV
const CurvilinearTrajectoryError & curvilinearError() const
MonitorElement * TrackEta
MonitorElement * TrackPtErr
MonitorElement * sip3dToPV
Sin< T >::type sin(const T &t)
MonitorElement * NumberOfLostRecHitsPerTrack
MonitorElement * TrackPyErr
uint8_t stopReason() const
ROOT::Math::SMatrix< double, 6, 6, ROOT::Math::MatRepSym< double, 6 > > AlgebraicSymMatrix66
Geom::Phi< T > phi() const
const CartesianTrajectoryError cartesianError() const
Geom::Theta< T > theta() const
bool doGeneralPropertiesPlots_
double etaError() const
error on eta
MonitorElement * Chi2oNDFVsPhi
int bunchCrossing() const
MonitorElement * Chi2ProbVsPhi
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)
MonitorElement * DistanceOfClosestApproachVsPhi
double phi() const
azimuthal angle of momentum vector
unsigned short numberOfLostHits() const
number of cases where track crossed a layer without getting a hit.
std::vector< Vertex > VertexCollection
collection of Vertex objects
bool doEffFromHitPatternVsBX_
MonitorElement * sipDzToBS
std::map< std::string, TkParameterMEs > TkParameterMEMap
MonitorElement * NumberOfLayersPerTrack[4]
MonitorElement * NhitVsPhi_HighPurity
MonitorElement * NumberOfRecHitsPerTrackVsTheta
MonitorElement * algorithm
double px() const
x coordinate of momentum vector
int pixelLayersWithMeasurement() const
bool doHitPropertiesPlots_
MonitorElement * DistanceOfClosestApproach
const math::XYZPoint & outerPosition() const
position of the outermost hit
bool doRecHitVsPhiVsEtaPerTrack_
MonitorElement * NumberOfLayersPerTrackVsTheta
MonitorElement * NumberOfRecHitsPerTrack_lumiFlag
std::map< std::string, TkRecHitsPerSubDetMEs > TkRecHitsPerSubDetMEMap
int trackerLayersWithMeasurement() const
const Point & position() const
position
MonitorElement * NumberOfValidRecHitVsPhiVsEtaPerTrack
int pixelEndcapLayersWithMeasurement() const
int numberOfValidStripTOBHits() const
bool doMeasurementStatePlots_
TrajectoryStateOnSurface innermostMeasurementState() const
MonitorElement * NumberOfRecHitsPerTrack
MonitorElement * Ptdist_HighPurity
int numberOfLostTrackerHits(HitCategory category) const
Geom::Theta< T > theta() const
const math::XYZPoint & innerPosition() const
position of the innermost hit
MonitorElement * NumberOfLayersVsPhiVsEtaPerTrack[4]
TrackAlgorithm algo() const
void disableSoftReset(MonitorElement *me)
MonitorElement * TESTDistanceOfClosestApproachToBSVsPhi
MonitorElement * NumberOfValidRecHitsPerTrackVsEta
std::unordered_map< Key, MonitorElement *, KeyHasher > hits_total_
MonitorElement * Chi2oNDFVsEta
MonitorElement * TransDCASig
MonitorElement * NumberOfMIRecHitVsPhiVsEtaPerTrack
MonitorElement * NumberOfRecHitsPerTrackVsPhi
MonitorElement * xPointOfClosestApproachVsZ0wrt000
double eta() const
pseudorapidity of momentum vector
int numberOfValidStripLayersWithMonoAndStereo(uint16_t stripdet, uint16_t layer) const
std::string qualityString_
int trackerLayersTotallyOffOrBad() const
int numberOfValidPixelBarrelHits() const
bool doEffFromHitPatternVsPU_
MonitorElement * NumberOfRecHitsPerTrackVsTheta
MonitorElement * NumberOfLostRecHitVsPhiVsEtaPerTrack
MonitorElement * xPointOfClosestApproach
MonitorElement * Chi2oNDFVsTheta
double chi2() const
chi-squared of the fit
MonitorElement * NumberOfMIRecHitsPerTrack
MonitorElement * bookProfile2D(Args &&...args)
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
double ndof() const
number of degrees of freedom of the fit
int stripTIDLayersWithMeasurement() const
MonitorElement * NumberOfValidRecHitsPerTrackVsEta
MonitorElement * NumberOfRecHitsPerTrackVsPhi
double pt() const
track transverse momentum
MonitorElement * dNdPt_HighPurity
void setBX(const edm::Event &)
double ptError() const
error on Pt (set to 1000 TeV if charge==0 for safety)
MonitorElement * xPointOfClosestApproachVsZ0wrtBS
double phiError() const
error on phi
MonitorElement * NumberOfValidRecHitsPerTrackVsPhi
MonitorElement * dNdPhi_HighPurity
MonitorElement * book1D(Args &&...args)
int numberOfValidStripTIDHits() const
MonitorElement * yPointOfClosestApproachToPV
const double EtaMin[kNumberCalorimeter]
MonitorElement * NumberOfRecHitsPerTrackVsEta
unsigned short numberOfValidHits() const
number of valid hits found
int numberOfValidStripTECHits() const
MonitorElement * sip2dToPV
void softReset(MonitorElement *me)
MonitorElement * Chi2oNDFVsTheta
TrajectoryStateOnSurface outermostMeasurementState() const
std::unordered_map< Key, MonitorElement *, KeyHasher > hits_inactive_
MonitorElement * sipDzToPV
MonitorElement * sipDxyToBS
MonitorElement * TrackEtaErr
void bookHistosForHitProperties(DQMStore::IBooker &ibooker)
MonitorElement * TrackPxErr
MonitorElement * DistanceOfClosestApproachVsEta
MonitorElement * zPointOfClosestApproachVsPhi
const AlgebraicSymMatrix66 & matrix() const
double pz() const
z coordinate of momentum vector
MonitorElement * TrackPzErr
double dz() const
dz parameter (= dsz/cos(lambda)). This is the track z0 w.r.t (0,0,0) only if the refPoint is close to...
MonitorElement * NumberOfLayersPerTrackVsEta
double dzError() const
error on dz
MonitorElement * dNdEta_HighPurity
MonitorElement * NumberOfRecHitsPerTrack
MonitorElement * TrackPhiErr
double vz() const
z coordinate of the reference point on track
void bookHistosForTrackerSpecific(DQMStore::IBooker &ibooker)
MonitorElement * stoppingSource
void fillHistosForState(const edm::EventSetup &iSetup, const reco::Track &track, std::string sname)
TrackAlgorithm originalAlgo() const
TrajectoryStateOnSurface TSOS
static TrackQuality qualityByName(const std::string &name)
double xError() const
error on x
MonitorElement * NumberOfLayersPerTrackVsEta
void setCurrentFolder(const std::string &fullpath)
edm::EDGetTokenT< reco::BeamSpot > beamSpotToken_
MonitorElement * sipDxyToPV
MonitorElement * book2D(Args &&...args)
virtual void analyze(const edm::Event &iEvent, const edm::EventSetup &iSetup, const reco::Track &track)
const HitPattern & hitPattern() const
Access the hit pattern, indicating in which Tracker layers the track has hits.
MonitorElement * xPointOfClosestApproachVsZ0wrtPV
MonitorElement * yPointOfClosestApproachVsZ0wrtBS
void fillHistosForEfficiencyFromHitPatter(const reco::Track &track, const std::string suffix, const unsigned int monitoring)
std::unordered_map< Key, MonitorElement *, KeyHasher > hits_missing_
int pixelBarrelLayersWithMeasurement() const
MonitorElement * DistanceOfClosestApproachVsTheta
void bookHistosForBeamSpot(DQMStore::IBooker &ibooker)
MonitorElement * Chi2ProbVsEta
int numberOfValidStripTIBHits() const
MonitorElement * TrackPhi
static const std::string algoNames[]
bool quality(const TrackQuality) const
Track quality.
int numberOfValidPixelEndcapHits() const
MonitorElement * zPointOfClosestApproachToPV
int trackerLayersWithoutMeasurement(HitCategory category) const
MonitorElement * DistanceOfClosestApproachToPV
double vy() const
y coordinate of the reference point on track
const AlgebraicSymMatrix55 & matrix() const
GlobalVector globalMomentum() const
MonitorElement * NumberOfLayersPerTrack
MonitorElement * DistanceOfClosestApproachToBSVsPhi
bool doLayersVsPhiVsEtaPerTrack_
MonitorElement * ValidFractionPerTrack
MonitorElement * NhitVsEta_HighPurity
void doSoftReset(DQMStore *dqmStore_)
int getNbinsX(void) const
get # of bins in X-axis
int stripTECLayersWithMeasurement() const
MonitorElement * yPointOfClosestApproachVsZ0wrt000
unsigned int good_vertices_
MonitorElement * TrackPtErrVsEta
MonitorElement * TESTDistanceOfClosestApproachToBS
virtual void initHisto(DQMStore::IBooker &ibooker, const edm::EventSetup &)
MonitorElement * NumberOfRecHitsPerTrackVsEta
int charge() const
track electric charge
const Point & position() const
position
TrajectoryStateOnSurface impactPointState() const
void bookHistosForState(std::string sname, DQMStore::IBooker &ibooker)
MonitorElement * yPointOfClosestApproach
void Reset(std::vector< TH2F > &depth)
MonitorElement * Chi2Prob
MonitorElement * TrackEtaPhiOuter
MonitorElement * yPointOfClosestApproachVsZ0wrtPV
void setNumberOfGoodVertices(const edm::Event &)
void bookHistosForEfficiencyFromHitPatter(DQMStore::IBooker &ibooker, const edm::EventSetup &iSetup, const std::string suffix)
double dxy() const
dxy parameter. (This is the transverse impact parameter w.r.t. to (0,0,0) ONLY if refPoint is close t...
MonitorElement * TrackTheta
MonitorElement * zPointOfClosestApproach
TrackAnalyzer(const edm::ParameterSet &)
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
MonitorElement * Chi2ProbVsEta
MonitorElement * TrackEtaPhiInner
MonitorElement * DistanceOfClosestApproachToBS
std::string monName[monQuantity::END]
Power< A, B >::type pow(const A &a, const B &b)
MonitorElement * Chi2oNDF
double yError() const
error on y
MonitorElement * Chi2oNDFVsEta
double py() const
y coordinate of momentum vector
double vx() const
x coordinate of the reference point on track
void undoSoftReset(DQMStore *dqmStore_)
int numberOfHits(HitCategory category) const
MonitorElement * NumberOfMORecHitVsPhiVsEtaPerTrack
MonitorElement * NumberOfMORecHitsPerTrack
MonitorElement * NumberOfLayersPerTrackVsPhi
MonitorElement * LongDCASig
edm::EDGetTokenT< reco::VertexCollection > pvToken_
MonitorElement * Chi2oNDFVsPhi