43 , NumberOfTracks(
NULL)
44 , NumberOfMeanRecHitsPerTrack(
NULL)
45 , NumberOfMeanLayersPerTrack(
NULL)
47 , FractionOfGoodTracks(
NULL)
49 , NumberOfSeeds_lumiFlag(
NULL)
50 , NumberOfTrackCandidates(
NULL)
52 , NumberEventsOfVsLS(
NULL)
53 , NumberOfTracksVsLS(
NULL)
55 , GoodTracksFractionVsLS(
NULL)
59 , NumberOfGoodPVtxVsLS(
NULL)
60 , NumberEventsOfVsBX (
NULL)
61 , NumberOfTracksVsBX(
NULL)
62 , GoodTracksFractionVsBX(
NULL)
63 , NumberOfRecHitsPerTrackVsBX(
NULL)
64 , NumberOfGoodPVtxVsBX(
NULL)
65 , NumberOfTracksVsBXlumi(
NULL)
66 , NumberOfTracksVsGoodPVtx(
NULL)
67 , NumberOfTracksVsPUPVtx(
NULL)
68 , NumberEventsOfVsGoodPVtx(
NULL)
69 , GoodTracksFractionVsGoodPVtx(
NULL)
70 , NumberOfRecHitsPerTrackVsGoodPVtx(
NULL)
71 , NumberOfPVtxVsGoodPVtx(
NULL)
72 , NumberEventsOfVsLUMI(
NULL)
73 , NumberOfTracksVsLUMI(
NULL)
74 , GoodTracksFractionVsLUMI(
NULL)
75 , NumberOfRecHitsPerTrackVsLUMI(
NULL)
76 , NumberOfGoodPVtxVsLUMI(
NULL)
77 , NumberOfTracks_lumiFlag(
NULL)
80 , builderName ( conf_.getParameter<std::
string>(
"TTRHBuilder"))
81 , doTrackerSpecific_ ( conf_.getParameter<bool>(
"doTrackerSpecific") )
82 , doLumiAnalysis ( conf_.getParameter<bool>(
"doLumiAnalysis"))
83 , doProfilesVsLS_ ( conf_.getParameter<bool>(
"doProfilesVsLS"))
84 , doAllPlots ( conf_.getParameter<bool>(
"doAllPlots"))
85 , doGeneralPropertiesPlots_( conf_.getParameter<bool>(
"doGeneralPropertiesPlots"))
86 , doHitPropertiesPlots_ ( conf_.getParameter<bool>(
"doHitPropertiesPlots"))
87 , doPUmonitoring_ ( conf_.getParameter<bool>(
"doPUmonitoring") )
89 , numSelection_ (conf_.getParameter<std::
string>(
"numCut"))
90 , denSelection_ (conf_.getParameter<std::
string>(
"denCut"))
91 , pvNDOF_ ( conf_.getParameter<int> (
"pvNDOF") )
109 allTrackToken_ = consumes<reco::TrackCollection>(alltrackProducer);
112 seedToken_ = consumes<edm::View<TrajectorySeed> >(seedProducer);
120 if (alltrackProducer.
label()==trackProducer.label()) doFractionPlot_ =
false;
127 if ( doPlotsVsBXlumi_ )
135 std::vector<edm::InputTag> primaryVertexInputTags =
conf_.
getParameter<std::vector<edm::InputTag> >(
"primaryVertexInputTags");
136 std::vector<edm::InputTag> selPrimaryVertexInputTags =
conf_.
getParameter<std::vector<edm::InputTag> >(
"selPrimaryVertexInputTags");
137 std::vector<std::string> pvLabels =
conf_.
getParameter<std::vector<std::string> > (
"pvLabels");
139 if (primaryVertexInputTags.size()==pvLabels.size() and primaryVertexInputTags.size()==selPrimaryVertexInputTags.size()) {
140 for (
size_t i=0;
i<primaryVertexInputTags.size();
i++) {
182 edm::LogWarning(
"TrackingMonitor") <<
"Qualty Name is invalid, using no quality criterea by default";
217 StateName !=
"OuterSurface" &&
218 StateName !=
"InnerSurface" &&
219 StateName !=
"ImpactPoint" &&
220 StateName !=
"default" &&
224 edm::LogWarning(
"TrackingMonitor") <<
"State Name is invalid, using 'ImpactPoint' by default";
236 histname =
"NumberOfTracks_" + CategoryName;
242 histname =
"NumberOfMeanRecHitsPerTrack_" + CategoryName;
247 histname =
"NumberOfMeanLayersPerTrack_" + CategoryName;
253 histname =
"FractionOfGoodTracks_" + CategoryName;
264 histname =
"NumberOfTracks_lumiFlag_" + CategoryName;
279 histname =
"NumberOfTracksVsLS_"+ CategoryName;
285 histname =
"NumberOfRecHitsPerTrackVsLS_" + CategoryName;
291 histname =
"NumberEventsVsLS_" + CategoryName;
300 histname =
"NumberOfGoodPVtxVsLS_" + CategoryName;
307 histname =
"GoodTracksFractionVsLS_"+ CategoryName;
316 int BXBin = 3564;
double BXMin = 0.5;
double BXMax = 3564.5;
318 histname =
"NumberEventsVsBX_" + CategoryName;
324 histname =
"NumberOfTracksVsBX_"+ CategoryName;
330 histname =
"NumberOfRecHitsPerTrackVsBX_" + CategoryName;
336 histname =
"NumberOfGoodPVtxVsBX_" + CategoryName;
343 histname =
"GoodTracksFractionVsBX_"+ CategoryName;
369 histname =
"NumberOfTracksVsGoodPVtx";
375 histname =
"NumberOfTracksVsPUPVtx";
381 histname =
"NumberEventsVsGoodPVtx";
388 histname =
"GoodTracksFractionVsGoodPVtx";
395 histname =
"NumberOfRecHitsPerTrackVsGoodPVtx";
401 histname =
"NumberOfPVtxVsGoodPVtx";
421 histname =
"NumberOfTracksVsLUMI";
428 histname =
"GoodTracksFractionVsLUMI";
435 histname =
"NumberOfRecHitsPerTrackVsLUMI";
444 histname =
"NumberOfGoodPVtxVsLUMI";
456 int BXlumiBin = BXlumiParameters.
getParameter<
int>(
"BXlumiBin");
457 double BXlumiMin = BXlumiParameters.
getParameter<
double>(
"BXlumiMin");
458 double BXlumiMax = BXlumiParameters.
getParameter<
double>(
"BXlumiMax");
460 histname =
"NumberOfTracksVsBXlumi_"+ CategoryName;
488 histname =
"NumberOfSeeds_"+ seedProducer.
label() +
"_"+ CategoryName;
493 if ( doSeedLumiAnalysis_ ) {
495 histname =
"NumberOfSeeds_lumiFlag_"+ seedProducer.
label() +
"_"+ CategoryName;
508 std::vector<double> histoMin,histoMax;
509 std::vector<int> histoBin;
519 setMaxMinBin(histoMin,histoMax,histoBin,NClusStrMin,NClusStrMax,NClusStrBin,NClusPxMin,NClusPxMax,NClusPxBin);
524 TKNoSeedBin, TKNoSeedMin, TKNoSeedMax)));
538 histname =
"NumberOfTrackCandidates_"+ tcProducer.
label() +
"_"+ CategoryName;
553 if ( doSeedLumiAnalysis_ )
561 std::vector<double> histoMin,histoMax;
562 std::vector<int> histoBin;
577 NClusStrMin,NClusStrMax,NClusStrBin,
578 NClusPxMin, NClusPxMax, NClusPxBin);
587 histoBin[
i], histoMin[i], histoMax[i],
588 NTrk2DBin,NTrk2DMin,NTrk2DMax
592 title =
"# of Clusters in (Pixel+Strip) Detectors";
637 if ( lumiScalers.
isValid() && lumiScalers->size() ) {
638 LumiScalersCollection::const_iterator scalit = lumiScalers->begin();
639 lumi = scalit->instantLumi();
658 int numberOfTracks_den = 0;
661 if (allTrackHandle.
isValid()) {
662 for (reco::TrackCollection::const_iterator track = allTrackHandle->begin();
663 track!=allTrackHandle->end(); ++track) {
665 numberOfTracks_den++;
673 pv0 = &pvHandle->front();
677 || pv0->
ndof() <
pvNDOF_ || pv0->
z() > 24.) pv0 =
nullptr;
683 int numberOfTracks = trackHandle->size();
684 int numberOfTracks_num = 0;
685 int numberOfTracks_pv0 = 0;
689 int totalRecHits = 0, totalLayers = 0;
694 for (reco::TrackCollection::const_iterator track = trackCollection.begin();
695 track!=trackCollection.end(); ++track) {
700 numberOfTracks_num++;
710 totalRecHits += track->numberOfValidHits();
711 totalLayers += track->hitPattern().trackerLayersWithMeasurement();
719 if (numberOfTracks_den > 0) frac =
static_cast<double>(numberOfTracks_num)/static_cast<double>(numberOfTracks_den);
737 if( numberOfTracks > 0 ) {
738 double meanRecHits =
static_cast<double>(totalRecHits) / static_cast<double>(numberOfTracks);
739 double meanLayers =
static_cast<double>(totalLayers) / static_cast<double>(numberOfTracks);
783 for( TrackCandidateCollection::const_iterator cand = theTCCollection.begin(); cand != theTCCollection.end(); ++cand) {
788 edm::LogWarning(
"TrackingMonitor") <<
"No Track Candidates in the event. Not filling associated histograms";
812 std::vector<int> NClus;
832 for(
size_t i=0;
i < seedHandle->size(); ++
i) {
840 edm::LogWarning(
"TrackingMonitor") <<
"No Trajectory seeds in the event. Not filling associated histograms";
847 std::vector<int> NClus;
862 size_t totalNumGoodPV = 0;
865 for (reco::VertexCollection::const_iterator
pv = pvHandle->begin();
866 pv != pvHandle->end(); ++
pv) {
869 if (
pv->isFake() ||
pv->tracksSize()==0)
continue;
872 if (
pv->ndof() <
pvNDOF_ ||
pv->z() > 24.)
continue;
878 if (totalNumGoodPV>1)
NumberOfTracksVsPUPVtx->
Fill( totalNumGoodPV-1,
double(numberOfTracks-numberOfTracks_pv0)/
double(totalNumGoodPV-1) );
881 for (reco::TrackCollection::const_iterator track = trackCollection.begin();
882 track!=trackCollection.end(); ++track) {
915 void TrackingMonitor::setMaxMinBin(std::vector<double> &arrayMin, std::vector<double> &arrayMax, std::vector<int> &arrayBin,
double smin,
double smax,
int sbin,
double pmin,
double pmax,
int pbin)
927 arrayMin[
i]=smin; arrayMax[
i]=smax; arrayBin[
i]=sbin;}
937 int ncluster_strip=-1;
945 ncluster_pix = (*pixel_clusters).dataSize();
946 ncluster_strip = (*strip_clusters).dataSize();
956 arrayNclus[
i]=ncluster_strip ;}
T getParameter(std::string const &) const
edm::EDGetTokenT< edm::View< TrajectorySeed > > seedToken_
std::vector< MonitorElement * > NumberOfTrkVsClusters
MonitorElement * FractionOfGoodTracks
virtual void setMaxMinBin(std::vector< double > &, std::vector< double > &, std::vector< int > &, double, double, int, double, double, int)
MonitorElement * NumberOfTracksVsBX
TrackAnalyzer * theTrackAnalyzer
MonitorElement * NumberOfRecHitsPerTrackVsBX
MonitorElement * NumberOfRecHitsPerTrackVsLUMI
virtual void setNclus(const edm::Event &, std::vector< int > &)
edm::EDGetTokenT< reco::TrackCollection > allTrackToken_
MonitorElement * bookProfile(Args &&...args)
MonitorElement * NumberEventsOfVsGoodPVtx
GetLumi * theLumiDetails_
bool getByToken(EDGetToken token, Handle< PROD > &result) const
#define DEFINE_FWK_MODULE(type)
std::vector< TrackCandidate > TrackCandidateCollection
std::vector< Track > TrackCollection
collection of Tracks
int bunchCrossing() const
MonitorElement * GoodTracksFractionVsLUMI
virtual void analyze(const edm::Event &, const edm::EventSetup &)
Provides a code based selection for trigger and DCS information in order to have no failing filters i...
MonitorElement * NumberOfPVtxVsGoodPVtx
MonitorElement * NumberOfMeanRecHitsPerTrack
MonitorElement * NumberOfTracksVsGoodPVtx
LuminosityBlockNumber_t luminosityBlock() const
const Point & position() const
position
MonitorElement * GoodTracksFractionVsLS
void setLumi(const edm::Event &, const edm::EventSetup &iSetup)
double getValue(const edm::Event &)
virtual void analyze(const edm::Event &iEvent, const edm::EventSetup &iSetup, const TrajectorySeed &seed, const reco::BeamSpot &bs, const edm::ESHandle< MagneticField > &theMF, const edm::ESHandle< TransientTrackingRecHitBuilder > &theTTRHBuilder)
MonitorElement * NumberOfSeeds
MonitorElement * NumberOfGoodPVtxVsLS
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
MonitorElement * NumberOfGoodPVtxVsLUMI
ConsumesCollector consumesCollector()
Use a ConsumesCollector to gather consumes information from helper functions.
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
bool accept(const edm::Event &event, const edm::EventSetup &setup)
To be called from analyze/filter() methods.
void setBX(const edm::Event &)
virtual void initHisto(DQMStore::IBooker &ibooker)
MonitorElement * book1D(Args &&...args)
MonitorElement * NumberOfTracksVsLS
Abs< T >::type abs(const T &t)
double z() const
y coordinate
MonitorElement * NumberOfTracksVsPUPVtx
edm::EDGetTokenT< LumiScalersCollection > lumiscalersToken_
MonitorElement * NumberEventsOfVsBX
bool runTrackBuildingAnalyzerForSeed
std::vector< MonitorElement * > SeedsVsClusters
MonitorElement * NumberOfRecHitsPerTrackVsLS
MonitorElement * GoodTracksFractionVsGoodPVtx
virtual void endRun(const edm::Run &, const edm::EventSetup &)
MonitorElement * NumberOfRecHitsPerTrackVsGoodPVtx
std::vector< VertexMonitor * > theVertexMonitor
MonitorElement * NumberOfTracksVsBXlumi
edm::InputTag stripClusterInputTag_
bool doGeneralPropertiesPlots_
TrackBuildingAnalyzer * theTrackBuildingAnalyzer
void setCurrentFolder(const std::string &fullpath)
MonitorElement * book2D(Args &&...args)
virtual void analyze(const edm::Event &iEvent, const edm::EventSetup &iSetup, const reco::Track &track)
edm::EDGetTokenT< reco::TrackCollection > trackToken_
edm::EDGetTokenT< TrackCandidateCollection > trackCandidateToken_
MonitorElement * NumberOfMeanLayersPerTrack
StringCutObjectSelector< reco::Track, true > denSelection_
StringCutObjectSelector< reco::Track, true > numSelection_
GenericTriggerEventFlag * genTriggerEventFlag_
MonitorElement * GoodTracksFractionVsBX
edm::InputTag pixelClusterInputTag_
MonitorElement * NumberOfSeeds_lumiFlag
TrackingMonitor(const edm::ParameterSet &)
edm::ESHandle< TransientTrackingRecHitBuilder > theTTRHBuilder
virtual void beginJob(void)
MonitorElement * NumberOfTrackCandidates
virtual void initHisto(DQMStore::IBooker &ibooker, const edm::EventSetup &)
edm::EDGetTokenT< edmNew::DetSetVector< SiPixelCluster > > pixelClustersToken_
void Reset(std::vector< TH2F > &depth)
edm::EDGetTokenT< edmNew::DetSetVector< SiStripCluster > > stripClustersToken_
void initRun(const edm::Run &run, const edm::EventSetup &setup)
To be called from beginRun() methods.
MonitorElement * NumberOfGoodPVtxVsBX
void setLumiFlag(void)
this ME is meant to be stored for each luminosity section
MonitorElement * NumberEventsOfVsLUMI
void setNumberOfGoodVertices(const edm::Event &)
MonitorElement * NumberOfTracks_lumiFlag
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
void Reset(void)
reset ME (ie. contents, errors, etc)
edm::EDGetTokenT< reco::BeamSpot > bsSrcToken_
MonitorElement * NumberOfTracks
size_t tracksSize() const
number of tracks
edm::EDGetTokenT< reco::VertexCollection > pvSrcToken_
MonitorElement * NumberEventsOfVsLS
MonitorElement * NumberOfTracksVsLUMI
virtual void beginLuminosityBlock(const edm::LuminosityBlock &lumi, const edm::EventSetup &eSetup)
std::vector< std::string > ClusterLabels