CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
List of all members | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes | Static Private Attributes
SplitVertexResolution Class Reference

#include <Alignment/OfflineValidation/plugins/SplitVertexResolution.cc>

Inheritance diagram for SplitVertexResolution:
edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::SharedResources > edm::one::EDAnalyzerBase edm::EDConsumerBase

Public Member Functions

 SplitVertexResolution (const edm::ParameterSet &)
 
 ~SplitVertexResolution () override
 
- Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::SharedResources >
 EDAnalyzer ()=default
 
 EDAnalyzer (const EDAnalyzer &)=delete
 
SerialTaskQueue * globalLuminosityBlocksQueue () final
 
SerialTaskQueue * globalRunsQueue () final
 
const EDAnalyzer & operator= (const EDAnalyzer &)=delete
 
bool wantsGlobalLuminosityBlocks () const final
 
bool wantsGlobalRuns () const final
 
bool wantsInputProcessBlocks () const final
 
bool wantsProcessBlocks () const final
 
- Public Member Functions inherited from edm::one::EDAnalyzerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzerBase ()
 
ModuleDescription const & moduleDescription () const
 
bool wantsStreamLuminosityBlocks () const
 
bool wantsStreamRuns () const
 
 ~EDAnalyzerBase () override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ESProxyIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
std::vector< ESProxyIndex > const & esGetTokenIndicesVector (edm::Transition iTrans) const
 
std::vector< ESRecordIndex >
const & 		esGetTokenRecordIndicesVector (edm::Transition iTrans) const
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector
< ProductResolverIndexAndSkipBit >
const & 		itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const * > *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBase & operator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 
static bool mysorter (reco::Track i, reco::Track j)
 
- Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 

Private Member Functions

void analyze (const edm::Event &, const edm::EventSetup &) override
 
virtual void beginEvent () final
 
void beginJob () override
 
void beginRun (edm::Run const &iEvent, edm::EventSetup const &) override
 
std::vector< TH1F * > bookResidualsHistogram (TFileDirectory dir, unsigned int theNOfBins, TString resType, TString varType)
 
template<std::size_t SIZE>
bool checkBinOrdering (std::array< float, SIZE > &bins)
 
void endJob () override
 
void endRun (edm::Run const &, edm::EventSetup const &) override
 
void fillTrendPlotByIndex (TH1F *trendPlot, std::vector< TH1F * > &h, PVValHelper::estimator fitPar_)
 
statmode::fitParams fitResiduals (TH1 *hist, bool singleTime=false)
 
statmode::fitParams fitResiduals_v0 (TH1 *hist)
 
std::pair< long long, long long > getRunTime (const edm::EventSetup &iSetup) const
 

Private Attributes

const int compressionSettings_
 
bool debug_
 
std::mt19937 engine_
 
pvEvent event_
 
TH1D * evtsByTrigger_
 
TH1F * h_avgSumPt
 
TH1F * h_diffX
 
TH1F * h_diffY
 
TH1F * h_diffZ
 
TH1F * h_errX
 
TH1F * h_errY
 
TH1F * h_errZ
 
TH1F * h_lumiFromConfig
 
TH1F * h_minWTrks1
 
TH1F * h_minWTrks2
 
TH1I * h_nFinalVertices
 
TH1I * h_nNonFakeVertices
 
TH1I * h_nOfflineVertices
 
TH1F * h_ntrks
 
TH1I * h_nVertices
 
TH1F * h_OrigVertexErrX
 
TH1F * h_OrigVertexErrY
 
TH1F * h_OrigVertexErrZ
 
TH1F * h_profileBinnings
 
TH1F * h_pullX
 
std::vector< TH1F * > h_pullX_Ntracks_
 
std::vector< TH1F * > h_pullX_Nvtx_
 
std::vector< TH1F * > h_pullX_sumPt_
 
TH1F * h_pullY
 
std::vector< TH1F * > h_pullY_Ntracks_
 
std::vector< TH1F * > h_pullY_Nvtx_
 
std::vector< TH1F * > h_pullY_sumPt_
 
TH1F * h_pullZ
 
std::vector< TH1F * > h_pullZ_Ntracks_
 
std::vector< TH1F * > h_pullZ_Nvtx_
 
std::vector< TH1F * > h_pullZ_sumPt_
 
TH1F * h_PVCL_subVtx1
 
TH1F * h_PVCL_subVtx2
 
std::vector< TH1F * > h_resolX_Ntracks_
 
std::vector< TH1F * > h_resolX_Nvtx_
 
std::vector< TH1F * > h_resolX_sumPt_
 
std::vector< TH1F * > h_resolY_Ntracks_
 
std::vector< TH1F * > h_resolY_Nvtx_
 
std::vector< TH1F * > h_resolY_sumPt_
 
std::vector< TH1F * > h_resolZ_Ntracks_
 
std::vector< TH1F * > h_resolZ_Nvtx_
 
std::vector< TH1F * > h_resolZ_sumPt_
 
TH1I * h_runEndTimes
 
TH1I * h_runFromConfig
 
TH1F * h_runNumber
 
TH1I * h_runStartTimes
 
TH1F * h_sumPt
 
TH1F * h_sumPt1
 
TH1F * h_sumPt2
 
TH1F * h_wTrks1
 
TH1F * h_wTrks2
 
int ievt
 
double intLumi_
 
int itrks
 
double minVtxNdf_
 
double minVtxWgt_
 
std::array< float, nTrackBins_+1 > myNTrack_bins_
 
std::array< float, nVtxBins_+1 > myNVtx_bins_
 
std::array< float, nPtBins_+1 > mypT_bins_
 
const double nVisTrackBins_
 
const double nVisVtxBins_
 
edm::Service< TFileService > outfile_
 
TH1F * p_pullX_vsNtracks
 
TH1F * p_pullX_vsNvtx
 
TH1F * p_pullX_vsSumPt
 
TH1F * p_pullY_vsNtracks
 
TH1F * p_pullY_vsNvtx
 
TH1F * p_pullY_vsSumPt
 
TH1F * p_pullZ_vsNtracks
 
TH1F * p_pullZ_vsNvtx
 
TH1F * p_pullZ_vsSumPt
 
TH1F * p_resolX_vsNtracks
 
TH1F * p_resolX_vsNvtx
 
TH1F * p_resolX_vsSumPt
 
TH1F * p_resolY_vsNtracks
 
TH1F * p_resolY_vsNvtx
 
TH1F * p_resolY_vsSumPt
 
TH1F * p_resolZ_vsNtracks
 
TH1F * p_resolZ_vsNvtx
 
TH1F * p_resolZ_vsSumPt
 
edm::InputTag pvsTag_
 
edm::EDGetTokenT
< reco::VertexCollection > 		pvsToken_
 
bool runControl_
 
std::vector< unsigned int > runControlNumbers_
 
edm::ESGetToken< RunInfo,
RunInfoRcd > 		runInfoToken_
 
std::map< unsigned int,
std::pair< long long, long
long > > 		runNumbersTimesLog_
 
bool storeNtuple_
 
const float sumpTEndScale_
 
const float sumpTStartScale_
 
TH1D * tksByTrigger_
 
edm::InputTag tracksTag_
 
edm::EDGetTokenT
< reco::TrackCollection > 		tracksToken_
 
edm::ESGetToken
< TransientTrackBuilder,
TransientTrackRecord > 		transientTrackBuilderToken_
 
TTree * tree_
 
std::map< std::string,
std::pair< int, int > > 		triggerMap_
 
edm::InputTag triggerResultsTag_ = edm::InputTag("TriggerResults", "", "HLT")
 
edm::EDGetTokenT
< edm::TriggerResults > 		triggerResultsToken_
 

Static Private Attributes

static constexpr double cmToUm = 10000.
 
static const int nPtBins_ = 30
 
static const int nTrackBins_ = 120
 
static const int nVtxBins_ = 60
 

Additional Inherited Members

- Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B > consumes (edm::InputTag tag) noexcept
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<Transition Tr = Transition::Event>
constexpr auto esConsumes ()
 
template<Transition Tr = Transition::Event>
auto esConsumes (ESInputTag tag)
 
template<Transition Tr = Transition::Event>
ESGetTokenGeneric esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
 Used with EventSetupRecord::doGet. More...
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
void resetItemsToGetFrom (BranchType iType)
 

Detailed Description

Definition at line 76 of file SplitVertexResolution.cc.

Constructor & Destructor Documentation

SplitVertexResolution::SplitVertexResolution ( const edm::ParameterSet &  iConfig)
explicit

Definition at line 270 of file SplitVertexResolution.cc.

References filterCSVwithJSON::copy, debug_, PVValHelper::generateBins(), edm::ParameterSet::getUntrackedParameter(), TFileService::kSharedResource, myNTrack_bins_, myNVtx_bins_, mypT_bins_, nPtBins_, nTrackBins_, nVtxBins_, runControlNumbers_, AlCaHLTBitMon_QueryRunRegistry::string, sumpTEndScale_, and sumpTStartScale_.

271  : compressionSettings_(iConfig.getUntrackedParameter<int>("compressionSettings", -1)),
272  storeNtuple_(iConfig.getParameter<bool>("storeNtuple")),
273  intLumi_(iConfig.getUntrackedParameter<double>("intLumi", 0.)),
274  debug_(iConfig.getUntrackedParameter<bool>("Debug", false)),
275  pvsTag_(iConfig.getParameter<edm::InputTag>("vtxCollection")),
276  pvsToken_(consumes<reco::VertexCollection>(pvsTag_)),
277  tracksTag_(iConfig.getParameter<edm::InputTag>("trackCollection")),
278  tracksToken_(consumes<reco::TrackCollection>(tracksTag_)),
279  triggerResultsToken_(consumes<edm::TriggerResults>(triggerResultsTag_)),
280  transientTrackBuilderToken_(
281  esConsumes<TransientTrackBuilder, TransientTrackRecord>(edm::ESInputTag("", "TransientTrackBuilder"))),
282  runInfoToken_(esConsumes<RunInfo, RunInfoRcd, edm::Transition::BeginRun>()),
283  minVtxNdf_(iConfig.getUntrackedParameter<double>("minVertexNdf")),
284  minVtxWgt_(iConfig.getUntrackedParameter<double>("minVertexMeanWeight")),
285  runControl_(iConfig.getUntrackedParameter<bool>("runControl", false)),
286  // binning
287  sumpTStartScale_(iConfig.getUntrackedParameter<double>("sumpTStartScale", 1.)),
288  sumpTEndScale_(iConfig.getUntrackedParameter<double>("sumpTEndScale", 1e3)),
289  nVisTrackBins_(iConfig.getUntrackedParameter<double>("nTrackBins", 120.)),
290  nVisVtxBins_(iConfig.getUntrackedParameter<double>("nVtxBins", 60.)) {
291  usesResource(TFileService::kSharedResource);
292 
293  std::vector<unsigned int> defaultRuns;
294  defaultRuns.push_back(0);
295  runControlNumbers_ = iConfig.getUntrackedParameter<std::vector<unsigned int>>("runControlNumber", defaultRuns);
296 
297  mypT_bins_ = PVValHelper::makeLogBins<float, nPtBins_>(sumpTStartScale_, sumpTEndScale_);
298 
299  // IMPORTANT
300  // first argument is start, second argument is the range so it's [-0.5;nTracksBins_-0.5]
301  std::vector<float> vect = PVValHelper::generateBins(nTrackBins_ + 1, -0.5, nTrackBins_);
302  std::copy(vect.begin(), vect.begin() + nTrackBins_ + 1, myNTrack_bins_.begin());
303 
304  vect.clear();
305 
306  // IMPORTANT
307  // first argument is start, second argument is the range so it's [-0.5;nVtxBins_-0.5]
308  vect = PVValHelper::generateBins(nVtxBins_ + 1, -0.5, nVtxBins_);
309  std::copy(vect.begin(), vect.begin() + nVtxBins_ + 1, myNVtx_bins_.begin());
310 
311  if (debug_) {
312  std::string toOutput = "";
313  for (const auto& pTbin : mypT_bins_ | boost::adaptors::indexed(1)) {
314  if (pTbin.index() != 1) {
315  toOutput += " ";
316  }
317  toOutput += fmt::sprintf("%.2f", pTbin.value());
318  if (pTbin.index() != nPtBins_ + 1) {
319  toOutput += ",";
320  }
321  }
322  edm::LogVerbatim("SplitVertexResolution") << "sum(pT) bins = [" << toOutput << "] \n";
323 
324  toOutput.clear();
325  for (const auto& tkbin : myNTrack_bins_ | boost::adaptors::indexed(1)) {
326  if (tkbin.index() != 1) {
327  toOutput += " ";
328  }
329  toOutput += fmt::sprintf("%.1f", tkbin.value());
330  if (tkbin.index() != nTrackBins_ + 1) {
331  toOutput += ",";
332  }
333  }
334  edm::LogVerbatim("SplitVertexResolution") << "n. track bins = [" << toOutput << "] \n";
335 
336  toOutput.clear();
337  for (const auto& vtxbin : myNVtx_bins_ | boost::adaptors::indexed(1)) {
338  if (vtxbin.index() != 1) {
339  toOutput += " ";
340  }
341  toOutput += fmt::sprintf("%.1f", vtxbin.value());
342  if (vtxbin.index() != nVtxBins_ + 1) {
343  toOutput += ",";
344  }
345  }
346  edm::LogVerbatim("SplitVertexResolution") << "n. vertices bins = [" << toOutput << "] \n";
347  }
348 }
TFileService::kSharedResource
static const std::string kSharedResource
Definition: TFileService.h:76
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
PVValHelper::generateBins
std::vector< float > generateBins(int n, float start, float range)
Definition: PVValidationHelpers.cc:161
SplitVertexResolution::pvsToken_
edm::EDGetTokenT< reco::VertexCollection > pvsToken_
Definition: SplitVertexResolution.cc:120
SplitVertexResolution::runInfoToken_
edm::ESGetToken< RunInfo, RunInfoRcd > runInfoToken_
Definition: SplitVertexResolution.cc:130
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
SplitVertexResolution::mypT_bins_
std::array< float, nPtBins_+1 > mypT_bins_
Definition: SplitVertexResolution.cc:257
SplitVertexResolution::myNVtx_bins_
std::array< float, nVtxBins_+1 > myNVtx_bins_
Definition: SplitVertexResolution.cc:265
SplitVertexResolution::triggerResultsToken_
edm::EDGetTokenT< edm::TriggerResults > triggerResultsToken_
Definition: SplitVertexResolution.cc:126
SplitVertexResolution::myNTrack_bins_
std::array< float, nTrackBins_+1 > myNTrack_bins_
Definition: SplitVertexResolution.cc:261
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
SplitVertexResolution::transientTrackBuilderToken_
edm::ESGetToken< TransientTrackBuilder, TransientTrackRecord > transientTrackBuilderToken_
Definition: SplitVertexResolution.cc:129
SplitVertexResolution::tracksToken_
edm::EDGetTokenT< reco::TrackCollection > tracksToken_
Definition: SplitVertexResolution.cc:123
SplitVertexResolution::runControlNumbers_
std::vector< unsigned int > runControlNumbers_
Definition: SplitVertexResolution.cc:138
SplitVertexResolution::~SplitVertexResolution ( )
overridedefault

Member Function Documentation

void SplitVertexResolution::analyze ( const edm::Event &  iEvent,
const edm::EventSetup &  iSetup 
)
overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 353 of file SplitVertexResolution.cc.

References accept(), muonTagProbeFilters_cff::allTracks, beginEvent(), pvCand::CL_subVtx1, pvCand::CL_subVtx2, cmToUm, counter, debug_, TransientVertex::degreesOfFreedom(), engine_, edm::EventID::event(), event_, edm::Event::eventAuxiliary(), pvEvent::eventNumber, PVValHelper::fillByIndex(), edm::Event::getByToken(), edm::EventSetup::getData(), h_avgSumPt, h_diffX, h_diffY, h_diffZ, h_errX, h_errY, h_errZ, h_minWTrks1, h_minWTrks2, h_nFinalVertices, h_nNonFakeVertices, h_nOfflineVertices, h_ntrks, h_nVertices, h_OrigVertexErrX, h_OrigVertexErrY, h_OrigVertexErrZ, h_pullX, h_pullX_Ntracks_, h_pullX_Nvtx_, h_pullX_sumPt_, h_pullY, h_pullY_Ntracks_, h_pullY_Nvtx_, h_pullY_sumPt_, h_pullZ, h_pullZ_Ntracks_, h_pullZ_Nvtx_, h_pullZ_sumPt_, h_PVCL_subVtx1, h_PVCL_subVtx2, h_resolX_Ntracks_, h_resolX_Nvtx_, h_resolX_sumPt_, h_resolY_Ntracks_, h_resolY_Nvtx_, h_resolY_sumPt_, h_resolZ_Ntracks_, h_resolZ_Nvtx_, h_resolZ_sumPt_, h_runNumber, h_sumPt, h_sumPt1, h_sumPt2, h_wTrks1, h_wTrks2, edm::EventBase::id(), ievt, pvCand::ipos, reco::Vertex::isFake(), TransientVertex::isValid(), itrks, edm::EventID::luminosityBlock(), pvEvent::luminosityBlockNumber, minVtxNdf_, minVtxWgt_, pvCand::minW_subVtx1, pvCand::minW_subVtx2, myNTrack_bins_, myNVtx_bins_, mypT_bins_, mysorter(), pvCand::n_subVtx1, pvCand::n_subVtx2, reco::Vertex::ndof(), nPtBins_, nTrackBins_, reco::Vertex::nTracks(), pvCand::nTrks, pvEvent::nVtx, nVtxBins_, TransientVertex::originalTracks(), funct::pow(), pvEvent::pvs, pvsToken_, edm::EventID::run(), edm::EventAuxiliary::run(), runControl_, runControlNumbers_, pvEvent::runNumber, edm::TriggerNames::size(), mathSSE::sqrt(), storeNtuple_, AlCaHLTBitMon_QueryRunRegistry::string, pvCand::sumPt_subVtx1, pvCand::sumPt_subVtx2, TransientVertex::totalChiSquared(), tracks, reco::Vertex::tracks_begin(), reco::Vertex::tracks_end(), reco::Vertex::tracksSize(), tracksToken_, TransientVertex::trackWeight(), transientTrackBuilderToken_, tree_, triggerMap_, edm::Event::triggerNames(), triggerResultsToken_, parallelization::uint(), AdaptiveVertexFitter::vertex(), beam_dqm_sourceclient-live_cfg::vertices, reco::Vertex::x(), pvCand::x_origVtx, pvCand::x_subVtx1, pvCand::x_subVtx2, pvCand::xErr_origVtx, pvCand::xErr_subVtx1, pvCand::xErr_subVtx2, reco::Vertex::xError(), reco::Vertex::y(), pvCand::y_origVtx, pvCand::y_subVtx1, pvCand::y_subVtx2, pvCand::yErr_origVtx, pvCand::yErr_subVtx1, pvCand::yErr_subVtx2, reco::Vertex::yError(), reco::Vertex::z(), pvCand::z_origVtx, pvCand::z_subVtx1, pvCand::z_subVtx2, pvCand::zErr_origVtx, pvCand::zErr_subVtx1, pvCand::zErr_subVtx2, and reco::Vertex::zError().

353  {
354  using namespace edm;
355 
356  // deterministic seed from the event number
357  // should not bias the result as the event number is already
358  // assigned randomly-enough
359  engine_.seed(iEvent.id().event() + (iEvent.id().luminosityBlock() << 10) + (iEvent.id().run() << 20));
360 
361  // first check if the event passes the run control
362  bool passesRunControl = false;
363 
364  if (runControl_) {
365  for (const auto& runControlNumber : runControlNumbers_) {
366  if (iEvent.eventAuxiliary().run() == runControlNumber) {
367  if (debug_) {
368  edm::LogInfo("SplitVertexResolution")
369  << " run number: " << iEvent.eventAuxiliary().run() << " keeping run:" << runControlNumber;
370  }
371  passesRunControl = true;
372  break;
373  }
374  }
375  if (!passesRunControl)
376  return;
377  }
378 
379  // Fill general info
380  h_runNumber->Fill(iEvent.id().run());
381 
382  ievt++;
383  edm::Handle<edm::TriggerResults> hltresults;
384  iEvent.getByToken(triggerResultsToken_, hltresults);
385 
386  const edm::TriggerNames& triggerNames_ = iEvent.triggerNames(*hltresults);
387  int ntrigs = hltresults->size();
388  //const std::vector<std::string>& triggernames = triggerNames_.triggerNames();
389 
390  beginEvent();
391 
392  // Fill general info
393  event_.runNumber = iEvent.id().run();
394  event_.luminosityBlockNumber = iEvent.id().luminosityBlock();
395  event_.eventNumber = iEvent.id().event();
396 
397  TransientTrackBuilder const& theB = iSetup.getData(transientTrackBuilderToken_);
398 
399  edm::Handle<reco::VertexCollection> vertices;
400  iEvent.getByToken(pvsToken_, vertices);
401  const reco::VertexCollection pvtx = *(vertices.product());
402 
403  event_.nVtx = pvtx.size();
404  int nOfflineVtx = pvtx.size();
405  h_nOfflineVertices->Fill(nOfflineVtx);
406 
407  edm::Handle<reco::TrackCollection> tracks;
408  iEvent.getByToken(tracksToken_, tracks);
409  itrks += tracks.product()->size();
410 
411  for (int itrig = 0; itrig != ntrigs; ++itrig) {
412  const std::string& trigName = triggerNames_.triggerName(itrig);
413  bool accept = hltresults->accept(itrig);
414  if (accept == 1) {
415  triggerMap_[trigName].first += 1;
416  triggerMap_[trigName].second += tracks.product()->size();
417  // triggerInfo.push_back(pair <string, int> (trigName, accept));
418  }
419  }
420 
421  int counter = 0;
422  int noFakecounter = 0;
423  int goodcounter = 0;
424 
425  for (auto pvIt = pvtx.cbegin(); pvIt != pvtx.cend(); ++pvIt) {
426  reco::Vertex iPV = *pvIt;
427  counter++;
428  if (iPV.isFake())
429  continue;
430  noFakecounter++;
431 
432  // vertex selection as in bs code
433  if (iPV.ndof() < minVtxNdf_ || (iPV.ndof() + 3.) / iPV.tracksSize() < 2 * minVtxWgt_)
434  continue;
435 
436  goodcounter++;
437  reco::TrackCollection allTracks;
438  reco::TrackCollection groupOne, groupTwo;
439  for (auto trki = iPV.tracks_begin(); trki != iPV.tracks_end(); ++trki) {
440  if (trki->isNonnull()) {
441  reco::TrackRef trk_now(tracks, (*trki).key());
442  allTracks.push_back(*trk_now);
443  }
444  }
445 
446  if (goodcounter > 1)
447  continue;
448 
449  // order with decreasing pt
450  std::sort(allTracks.begin(), allTracks.end(), mysorter);
451 
452  int ntrks = allTracks.size();
453  h_ntrks->Fill(ntrks);
454 
455  // discard lowest pt track
456  uint even_ntrks;
457  ntrks % 2 == 0 ? even_ntrks = ntrks : even_ntrks = ntrks - 1;
458 
459  // split into two sets equally populated
460  for (uint tracksIt = 0; tracksIt < even_ntrks; tracksIt = tracksIt + 2) {
461  reco::Track firstTrk = allTracks.at(tracksIt);
462  reco::Track secondTrk = allTracks.at(tracksIt + 1);
463  auto dis = std::uniform_int_distribution<>(0, 1); // [0, 1]
464 
465  if (dis(engine_) > 0.5) {
466  groupOne.push_back(firstTrk);
467  groupTwo.push_back(secondTrk);
468  } else {
469  groupOne.push_back(secondTrk);
470  groupTwo.push_back(firstTrk);
471  }
472  }
473 
474  if (!(groupOne.size() >= 2 && groupTwo.size() >= 2))
475  continue;
476 
477  h_OrigVertexErrX->Fill(iPV.xError() * cmToUm);
478  h_OrigVertexErrY->Fill(iPV.yError() * cmToUm);
479  h_OrigVertexErrZ->Fill(iPV.zError() * cmToUm);
480 
481  float sumPt = 0, sumPt1 = 0, sumPt2 = 0, avgSumPt = 0;
482 
483  // refit the two sets of tracks
484  std::vector<reco::TransientTrack> groupOne_ttks;
485  groupOne_ttks.clear();
486  for (auto itrk = groupOne.cbegin(); itrk != groupOne.cend(); itrk++) {
487  reco::TransientTrack tmpTransientTrack = theB.build(*itrk);
488  groupOne_ttks.push_back(tmpTransientTrack);
489  sumPt1 += itrk->pt();
490  sumPt += itrk->pt();
491  }
492 
493  AdaptiveVertexFitter pvFitter;
494  TransientVertex pvOne = pvFitter.vertex(groupOne_ttks);
495  if (!pvOne.isValid())
496  continue;
497 
498  reco::Vertex onePV = pvOne;
499 
500  std::vector<reco::TransientTrack> groupTwo_ttks;
501  groupTwo_ttks.clear();
502  for (auto itrk = groupTwo.cbegin(); itrk != groupTwo.cend(); itrk++) {
503  reco::TransientTrack tmpTransientTrack = theB.build(*itrk);
504  groupTwo_ttks.push_back(tmpTransientTrack);
505  sumPt2 += itrk->pt();
506  sumPt += itrk->pt();
507  }
508 
509  // average sumPt
510  avgSumPt = (sumPt1 + sumPt2) / 2.;
511  h_avgSumPt->Fill(avgSumPt);
512 
513  TransientVertex pvTwo = pvFitter.vertex(groupTwo_ttks);
514  if (!pvTwo.isValid())
515  continue;
516 
517  reco::Vertex twoPV = pvTwo;
518 
519  float theminW1 = 1.;
520  float theminW2 = 1.;
521  for (auto otrk = pvOne.originalTracks().cbegin(); otrk != pvOne.originalTracks().cend(); ++otrk) {
522  h_wTrks1->Fill(pvOne.trackWeight(*otrk));
523  if (pvOne.trackWeight(*otrk) < theminW1) {
524  theminW1 = pvOne.trackWeight(*otrk);
525  }
526  }
527  for (auto otrk = pvTwo.originalTracks().cbegin(); otrk != pvTwo.originalTracks().end(); ++otrk) {
528  h_wTrks2->Fill(pvTwo.trackWeight(*otrk));
529  if (pvTwo.trackWeight(*otrk) < theminW2) {
530  theminW2 = pvTwo.trackWeight(*otrk);
531  }
532  }
533 
534  h_sumPt->Fill(sumPt);
535 
536  int half_trks = twoPV.nTracks();
537 
538  const double invSqrt2 = 1. / std::sqrt(2.);
539 
540  double deltaX = (twoPV.x() - onePV.x());
541  double deltaY = (twoPV.y() - onePV.y());
542  double deltaZ = (twoPV.z() - onePV.z());
543 
544  double resX = deltaX * invSqrt2;
545  double resY = deltaY * invSqrt2;
546  double resZ = deltaZ * invSqrt2;
547 
548  h_diffX->Fill(resX * cmToUm);
549  h_diffY->Fill(resY * cmToUm);
550  h_diffZ->Fill(resZ * cmToUm);
551 
552  double errX = sqrt(pow(twoPV.xError(), 2) + pow(onePV.xError(), 2));
553  double errY = sqrt(pow(twoPV.yError(), 2) + pow(onePV.yError(), 2));
554  double errZ = sqrt(pow(twoPV.zError(), 2) + pow(onePV.zError(), 2));
555 
556  h_errX->Fill(errX * cmToUm);
557  h_errY->Fill(errY * cmToUm);
558  h_errZ->Fill(errZ * cmToUm);
559 
560  h_pullX->Fill(deltaX / errX);
561  h_pullY->Fill(deltaY / errY);
562  h_pullZ->Fill(deltaZ / errZ);
563 
564  // filling the pT-binned distributions
565 
566  for (int ipTBin = 0; ipTBin < nPtBins_; ipTBin++) {
567  float pTF = mypT_bins_[ipTBin];
568  float pTL = mypT_bins_[ipTBin + 1];
569 
570  if (avgSumPt >= pTF && avgSumPt < pTL) {
571  PVValHelper::fillByIndex(h_resolX_sumPt_, ipTBin, resX * cmToUm, "1");
572  PVValHelper::fillByIndex(h_resolY_sumPt_, ipTBin, resY * cmToUm, "2");
573  PVValHelper::fillByIndex(h_resolZ_sumPt_, ipTBin, resZ * cmToUm, "3");
574 
575  PVValHelper::fillByIndex(h_pullX_sumPt_, ipTBin, deltaX / errX, "4");
576  PVValHelper::fillByIndex(h_pullY_sumPt_, ipTBin, deltaY / errY, "5");
577  PVValHelper::fillByIndex(h_pullZ_sumPt_, ipTBin, deltaZ / errZ, "6");
578  }
579  }
580 
581  // filling the track multeplicity binned distributions
582 
583  for (int inTrackBin = 0; inTrackBin < nTrackBins_; inTrackBin++) {
584  float nTrackF = myNTrack_bins_[inTrackBin];
585  float nTrackL = myNTrack_bins_[inTrackBin + 1];
586  if (ntrks >= nTrackF && ntrks < nTrackL) {
587  //if (ntrks == inTrackBin) {
588  PVValHelper::fillByIndex(h_resolX_Ntracks_, inTrackBin, resX * cmToUm, "7");
589  PVValHelper::fillByIndex(h_resolY_Ntracks_, inTrackBin, resY * cmToUm, "8");
590  PVValHelper::fillByIndex(h_resolZ_Ntracks_, inTrackBin, resZ * cmToUm, "9");
591 
592  PVValHelper::fillByIndex(h_pullX_Ntracks_, inTrackBin, deltaX / errX, "10");
593  PVValHelper::fillByIndex(h_pullY_Ntracks_, inTrackBin, deltaY / errY, "11");
594  PVValHelper::fillByIndex(h_pullZ_Ntracks_, inTrackBin, deltaZ / errZ, "12");
595  }
596  }
597 
598  // filling the vertex multeplicity binned distributions
599 
600  for (int inVtxBin = 0; inVtxBin < nVtxBins_; inVtxBin++) {
601  float nVtxF = myNVtx_bins_[inVtxBin];
602  float nVtxL = myNVtx_bins_[inVtxBin + 1];
603  if (nOfflineVtx >= nVtxF && nOfflineVtx < nVtxL) {
604  //if (nOfflineVtx == inVtxBin) {
605  PVValHelper::fillByIndex(h_resolX_Nvtx_, inVtxBin, deltaX * cmToUm, "7");
606  PVValHelper::fillByIndex(h_resolY_Nvtx_, inVtxBin, deltaY * cmToUm, "8");
607  PVValHelper::fillByIndex(h_resolZ_Nvtx_, inVtxBin, deltaZ * cmToUm, "9");
608 
609  PVValHelper::fillByIndex(h_pullX_Nvtx_, inVtxBin, deltaX / errX, "10");
610  PVValHelper::fillByIndex(h_pullY_Nvtx_, inVtxBin, deltaY / errY, "11");
611  PVValHelper::fillByIndex(h_pullZ_Nvtx_, inVtxBin, deltaZ / errZ, "12");
612  }
613  }
614 
615  h_sumPt1->Fill(sumPt1);
616  h_sumPt2->Fill(sumPt2);
617 
618  h_minWTrks1->Fill(theminW1);
619  h_minWTrks2->Fill(theminW2);
620 
621  h_PVCL_subVtx1->Fill(TMath::Prob(pvOne.totalChiSquared(), (int)(pvOne.degreesOfFreedom())));
622  h_PVCL_subVtx2->Fill(TMath::Prob(pvTwo.totalChiSquared(), (int)(pvTwo.degreesOfFreedom())));
623 
624  // fill ntuples
625  pvCand thePV;
626  thePV.ipos = counter;
627  thePV.nTrks = ntrks;
628 
629  thePV.x_origVtx = iPV.x();
630  thePV.y_origVtx = iPV.y();
631  thePV.z_origVtx = iPV.z();
632 
633  thePV.xErr_origVtx = iPV.xError();
634  thePV.yErr_origVtx = iPV.yError();
635  thePV.zErr_origVtx = iPV.zError();
636 
637  thePV.n_subVtx1 = half_trks;
638  thePV.x_subVtx1 = onePV.x();
639  thePV.y_subVtx1 = onePV.y();
640  thePV.z_subVtx1 = onePV.z();
641 
642  thePV.xErr_subVtx1 = onePV.xError();
643  thePV.yErr_subVtx1 = onePV.yError();
644  thePV.zErr_subVtx1 = onePV.zError();
645  thePV.sumPt_subVtx1 = sumPt1;
646 
647  thePV.n_subVtx2 = half_trks;
648  thePV.x_subVtx2 = twoPV.x();
649  thePV.y_subVtx2 = twoPV.y();
650  thePV.z_subVtx2 = twoPV.z();
651 
652  thePV.xErr_subVtx2 = twoPV.xError();
653  thePV.yErr_subVtx2 = twoPV.yError();
654  thePV.zErr_subVtx2 = twoPV.zError();
655  thePV.sumPt_subVtx2 = sumPt2;
656 
657  thePV.CL_subVtx1 = TMath::Prob(pvOne.totalChiSquared(), (int)(pvOne.degreesOfFreedom()));
658  thePV.CL_subVtx2 = TMath::Prob(pvTwo.totalChiSquared(), (int)(pvTwo.degreesOfFreedom()));
659 
660  thePV.minW_subVtx1 = theminW1;
661  thePV.minW_subVtx2 = theminW2;
662 
663  event_.pvs.push_back(thePV);
664 
665  } // loop on the vertices
666 
667  // fill the histogram of vertices per event
668  h_nVertices->Fill(counter);
669  h_nNonFakeVertices->Fill(noFakecounter);
670  h_nFinalVertices->Fill(goodcounter);
671 
672  if (storeNtuple_) {
673  tree_->Fill();
674  }
675 }
edm::TriggerNames::size
std::size_t size() const
Definition: TriggerNames.cc:59
edm::EventID::run
RunNumber_t run() const
Definition: EventID.h:38
edm::EventID::event
EventNumber_t event() const
Definition: EventID.h:40
pvEvent::eventNumber
int eventNumber
Definition: pvTree.h:58
pvCand::z_origVtx
float z_origVtx
Definition: pvTree.h:16
PVValHelper::fillByIndex
void fillByIndex(std::vector< TH1F * > &h, unsigned int index, double x, std::string tag="")
Definition: PVValidationHelpers.cc:42
SplitVertexResolution::h_resolZ_Nvtx_
std::vector< TH1F * > h_resolZ_Nvtx_
Definition: SplitVertexResolution.cc:207
SplitVertexResolution::h_pullY_sumPt_
std::vector< TH1F * > h_pullY_sumPt_
Definition: SplitVertexResolution.cc:223
pvCand::sumPt_subVtx2
float sumPt_subVtx2
Definition: pvTree.h:40
pvEvent::nVtx
int nVtx
Definition: pvTree.h:60
SplitVertexResolution::h_pullZ_sumPt_
std::vector< TH1F * > h_pullZ_sumPt_
Definition: SplitVertexResolution.cc:224
pvCand::y_subVtx1
float y_subVtx1
Definition: pvTree.h:24
SplitVertexResolution::h_resolX_sumPt_
std::vector< TH1F * > h_resolX_sumPt_
Definition: SplitVertexResolution.cc:197
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:539
reco::Vertex::zError
double zError() const
error on z
Definition: Vertex.h:141
pvCand::zErr_subVtx1
float zErr_subVtx1
Definition: pvTree.h:29
reco::Vertex::y
double y() const
y coordinate
Definition: Vertex.h:131
SplitVertexResolution::pvsToken_
edm::EDGetTokenT< reco::VertexCollection > pvsToken_
Definition: SplitVertexResolution.cc:120
TransientVertex::totalChiSquared
float totalChiSquared() const
Definition: TransientVertex.h:188
pvCand::y_origVtx
float y_origVtx
Definition: pvTree.h:15
reco::TrackCollection
std::vector< Track > TrackCollection
collection of Tracks
Definition: TrackFwd.h:14
edm::EventAuxiliary::run
RunNumber_t run() const
Definition: EventAuxiliary.h:73
reco::VertexCollection
std::vector< Vertex > VertexCollection
collection of Vertex objects
Definition: VertexFwd.h:9
tracks
auto const & tracks
cannot be loose
Definition: CAHitNtupletGeneratorKernelsImpl.h:148
pvEvent::runNumber
int runNumber
Definition: pvTree.h:56
beam_dqm_sourceclient-live_cfg.vertices
tuple vertices
Definition: beam_dqm_sourceclient-live_cfg.py:130
edm::Event::eventAuxiliary
EventAuxiliary const & eventAuxiliary() const override
Definition: Event.h:95
pvCand::minW_subVtx1
float minW_subVtx1
Definition: pvTree.h:45
pvCand::x_origVtx
float x_origVtx
Definition: pvTree.h:14
AdaptiveVertexFitter::vertex
CachingVertex< 5 > vertex(const std::vector< reco::TransientTrack > &) const override
Definition: AdaptiveVertexFitter.cc:158
accept
bool accept(const edm::Event &event, const edm::TriggerResults &triggerTable, const std::string &triggerPath)
Definition: TopDQMHelpers.h:31
edm::EventID::luminosityBlock
LuminosityBlockNumber_t luminosityBlock() const
Definition: EventID.h:39
SplitVertexResolution::h_pullY_Ntracks_
std::vector< TH1F * > h_pullY_Ntracks_
Definition: SplitVertexResolution.cc:227
pvCand::CL_subVtx2
float CL_subVtx2
Definition: pvTree.h:43
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
pvCand
Definition: pvTree.h:9
pvCand::y_subVtx2
float y_subVtx2
Definition: pvTree.h:34
edm::Event::triggerNames
edm::TriggerNames const & triggerNames(edm::TriggerResults const &triggerResults) const override
Definition: Event.cc:275
edm::EventSetup::getData
bool getData(T &iHolder) const
Definition: EventSetup.h:122
pvCand::x_subVtx1
float x_subVtx1
Definition: pvTree.h:23
pvCand::minW_subVtx2
float minW_subVtx2
Definition: pvTree.h:46
SplitVertexResolution::mypT_bins_
std::array< float, nPtBins_+1 > mypT_bins_
Definition: SplitVertexResolution.cc:257
pvCand::n_subVtx2
int n_subVtx2
Definition: pvTree.h:32
SplitVertexResolution::cmToUm
static constexpr double cmToUm
Definition: SplitVertexResolution.cc:140
TransientVertex::originalTracks
std::vector< reco::TransientTrack > const & originalTracks() const
Definition: TransientVertex.h:200
TransientVertex::degreesOfFreedom
float degreesOfFreedom() const
Definition: TransientVertex.h:190
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
pvEvent::pvs
std::vector< pvCand > pvs
Definition: pvTree.h:62
pvCand::zErr_origVtx
float zErr_origVtx
Definition: pvTree.h:20
reco::Vertex::tracks_end
trackRef_iterator tracks_end() const
last iterator over tracks
Definition: Vertex.h:110
reco::Vertex::z
double z() const
z coordinate
Definition: Vertex.h:133
pvCand::sumPt_subVtx1
float sumPt_subVtx1
Definition: pvTree.h:30
pvCand::x_subVtx2
float x_subVtx2
Definition: pvTree.h:33
pvCand::n_subVtx1
int n_subVtx1
Definition: pvTree.h:22
reco::Vertex::tracks_begin
trackRef_iterator tracks_begin() const
first iterator over tracks
Definition: Vertex.h:108
SplitVertexResolution::myNVtx_bins_
std::array< float, nVtxBins_+1 > myNVtx_bins_
Definition: SplitVertexResolution.cc:265
SplitVertexResolution::h_resolX_Nvtx_
std::vector< TH1F * > h_resolX_Nvtx_
Definition: SplitVertexResolution.cc:205
SplitVertexResolution::triggerResultsToken_
edm::EDGetTokenT< edm::TriggerResults > triggerResultsToken_
Definition: SplitVertexResolution.cc:126
pvCand::xErr_subVtx1
float xErr_subVtx1
Definition: pvTree.h:27
TransientVertex::trackWeight
float trackWeight(const reco::TransientTrack &track) const
Definition: TransientVertex.cc:241
reco::Vertex::tracksSize
size_t tracksSize() const
number of tracks
Definition: Vertex.h:112
SplitVertexResolution::triggerMap_
std::map< std::string, std::pair< int, int > > triggerMap_
Definition: SplitVertexResolution.cc:267
reco::Vertex::ndof
double ndof() const
Definition: Vertex.h:123
pvCand::nTrks
int nTrks
Definition: pvTree.h:11
pvCand::ipos
int ipos
Definition: pvTree.h:12
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
reco::Vertex::x
double x() const
x coordinate
Definition: Vertex.h:129
SplitVertexResolution::mysorter
static bool mysorter(reco::Track i, reco::Track j)
Definition: SplitVertexResolution.cc:82
SplitVertexResolution::myNTrack_bins_
std::array< float, nTrackBins_+1 > myNTrack_bins_
Definition: SplitVertexResolution.cc:261
SplitVertexResolution::h_pullY_Nvtx_
std::vector< TH1F * > h_pullY_Nvtx_
Definition: SplitVertexResolution.cc:231
reco::Vertex::xError
double xError() const
error on x
Definition: Vertex.h:137
SplitVertexResolution::h_resolZ_Ntracks_
std::vector< TH1F * > h_resolZ_Ntracks_
Definition: SplitVertexResolution.cc:203
reco::Vertex::isFake
bool isFake() const
Definition: Vertex.h:76
pvCand::z_subVtx1
float z_subVtx1
Definition: pvTree.h:25
SplitVertexResolution::h_pullX_Nvtx_
std::vector< TH1F * > h_pullX_Nvtx_
Definition: SplitVertexResolution.cc:230
SplitVertexResolution::beginEvent
virtual void beginEvent() final
Definition: SplitVertexResolution.cc:677
pvCand::CL_subVtx1
float CL_subVtx1
Definition: pvTree.h:42
reco::Vertex::nTracks
unsigned int nTracks(float minWeight=0.5) const
Returns the number of tracks in the vertex with weight above minWeight.
Definition: Vertex.cc:131
pvEvent::luminosityBlockNumber
int luminosityBlockNumber
Definition: pvTree.h:57
pvCand::yErr_subVtx2
float yErr_subVtx2
Definition: pvTree.h:38
pvCand::xErr_subVtx2
float xErr_subVtx2
Definition: pvTree.h:37
SplitVertexResolution::h_pullZ_Nvtx_
std::vector< TH1F * > h_pullZ_Nvtx_
Definition: SplitVertexResolution.cc:232
muonTagProbeFilters_cff.allTracks
tuple allTracks
Definition: muonTagProbeFilters_cff.py:22
SplitVertexResolution::h_resolY_Nvtx_
std::vector< TH1F * > h_resolY_Nvtx_
Definition: SplitVertexResolution.cc:206
pvCand::z_subVtx2
float z_subVtx2
Definition: pvTree.h:35
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:59
pvCand::yErr_origVtx
float yErr_origVtx
Definition: pvTree.h:19
counter
static std::atomic< unsigned int > counter
Definition: SharedResourceNames.cc:18
SplitVertexResolution::transientTrackBuilderToken_
edm::ESGetToken< TransientTrackBuilder, TransientTrackRecord > transientTrackBuilderToken_
Definition: SplitVertexResolution.cc:129
SplitVertexResolution::h_pullX_sumPt_
std::vector< TH1F * > h_pullX_sumPt_
Definition: SplitVertexResolution.cc:222
pvCand::zErr_subVtx2
float zErr_subVtx2
Definition: pvTree.h:39
pvCand::xErr_origVtx
float xErr_origVtx
Definition: pvTree.h:18
pvCand::yErr_subVtx1
float yErr_subVtx1
Definition: pvTree.h:28
SplitVertexResolution::h_pullX_Ntracks_
std::vector< TH1F * > h_pullX_Ntracks_
Definition: SplitVertexResolution.cc:226
SplitVertexResolution::tracksToken_
edm::EDGetTokenT< reco::TrackCollection > tracksToken_
Definition: SplitVertexResolution.cc:123
SplitVertexResolution::h_pullZ_Ntracks_
std::vector< TH1F * > h_pullZ_Ntracks_
Definition: SplitVertexResolution.cc:228
SplitVertexResolution::h_resolY_sumPt_
std::vector< TH1F * > h_resolY_sumPt_
Definition: SplitVertexResolution.cc:198
TransientVertex::isValid
bool isValid() const
Definition: TransientVertex.h:195
SplitVertexResolution::runControlNumbers_
std::vector< unsigned int > runControlNumbers_
Definition: SplitVertexResolution.cc:138
SplitVertexResolution::h_resolX_Ntracks_
std::vector< TH1F * > h_resolX_Ntracks_
Definition: SplitVertexResolution.cc:201
SplitVertexResolution::h_resolY_Ntracks_
std::vector< TH1F * > h_resolY_Ntracks_
Definition: SplitVertexResolution.cc:202
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:29
reco::Vertex::yError
double yError() const
error on y
Definition: Vertex.h:139
SplitVertexResolution::h_resolZ_sumPt_
std::vector< TH1F * > h_resolZ_sumPt_
Definition: SplitVertexResolution.cc:199
void SplitVertexResolution::beginEvent ( )
finalprivatevirtual

Definition at line 677 of file SplitVertexResolution.cc.

References event_, pvEvent::nVtx, and pvEvent::pvs.

Referenced by analyze().

677  {
678  event_.pvs.clear();
679  event_.nVtx = -1;
680 }
pvEvent::nVtx
int nVtx
Definition: pvTree.h:60
pvEvent::pvs
std::vector< pvCand > pvs
Definition: pvTree.h:62
void SplitVertexResolution::beginJob ( )
overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 700 of file SplitVertexResolution.cc.

References bookResidualsHistogram(), checkBinOrdering(), compressionSettings_, event_, TFileService::file(), h_avgSumPt, h_diffX, h_diffY, h_diffZ, h_errX, h_errY, h_errZ, h_lumiFromConfig, h_minWTrks1, h_minWTrks2, h_nFinalVertices, h_nNonFakeVertices, h_nOfflineVertices, h_ntrks, h_nVertices, h_OrigVertexErrX, h_OrigVertexErrY, h_OrigVertexErrZ, h_profileBinnings, h_pullX, h_pullX_Ntracks_, h_pullX_Nvtx_, h_pullX_sumPt_, h_pullY, h_pullY_Ntracks_, h_pullY_Nvtx_, h_pullY_sumPt_, h_pullZ, h_pullZ_Ntracks_, h_pullZ_Nvtx_, h_pullZ_sumPt_, h_PVCL_subVtx1, h_PVCL_subVtx2, h_resolX_Ntracks_, h_resolX_Nvtx_, h_resolX_sumPt_, h_resolY_Ntracks_, h_resolY_Nvtx_, h_resolY_sumPt_, h_resolZ_Ntracks_, h_resolZ_Nvtx_, h_resolZ_sumPt_, h_runFromConfig, h_runNumber, h_sumPt, h_sumPt1, h_sumPt2, h_wTrks1, h_wTrks2, ievt, intLumi_, itrks, TFileDirectory::make(), TFileService::make(), TFileService::mkdir(), myNTrack_bins_, myNVtx_bins_, mypT_bins_, nPtBins_, nTrackBins_, nVisTrackBins_, nVisVtxBins_, nVtxBins_, outfile_, p_pullX_vsNtracks, p_pullX_vsNvtx, p_pullX_vsSumPt, p_pullY_vsNtracks, p_pullY_vsNvtx, p_pullY_vsSumPt, p_pullZ_vsNtracks, p_pullZ_vsNvtx, p_pullZ_vsSumPt, p_resolX_vsNtracks, p_resolX_vsNvtx, p_resolX_vsSumPt, p_resolY_vsNtracks, p_resolY_vsNvtx, p_resolY_vsSumPt, p_resolZ_vsNtracks, p_resolZ_vsNvtx, p_resolZ_vsSumPt, submitPVValidationJobs::run, runControlNumbers_, sumpTEndScale_, sumpTStartScale_, and tree_.

700  {
701  ievt = 0;
702  itrks = 0;
703 
704  if (compressionSettings_ > 0) {
705  outfile_->file().SetCompressionSettings(compressionSettings_);
706  }
707 
708  // binning histogram
709  TFileDirectory BinningFeatures = outfile_->mkdir("BinningFeatures");
710  // sumpT
711  h_profileBinnings = BinningFeatures.make<TH1F>("h_profileBinnings", "profile binnings", 4, -0.5, 3.5);
712  h_profileBinnings->GetXaxis()->SetBinLabel(1, "#sum p_{T} min");
713  h_profileBinnings->GetXaxis()->SetBinLabel(2, "#sum p_{T} max");
714  // n. tracks
715  h_profileBinnings->GetXaxis()->SetBinLabel(3, "n. track bins");
716  // n. vertices
717  h_profileBinnings->GetXaxis()->SetBinLabel(4, "n. vertices bins");
718 
719  // fill the binning book-keeping
720  h_profileBinnings->SetBinContent(1, sumpTStartScale_);
721  h_profileBinnings->SetBinContent(2, sumpTEndScale_);
722  h_profileBinnings->SetBinContent(3, nVisTrackBins_);
723  h_profileBinnings->SetBinContent(4, nVisVtxBins_);
724 
725  // luminosity histo
726  TFileDirectory EventFeatures = outfile_->mkdir("EventFeatures");
727  h_lumiFromConfig =
728  EventFeatures.make<TH1F>("h_lumiFromConfig", "luminosity from config;;luminosity of present run", 1, -0.5, 0.5);
729  h_lumiFromConfig->SetBinContent(1, intLumi_);
730 
731  h_runFromConfig = EventFeatures.make<TH1I>("h_runFromConfig",
732  "run number from config;;run number (from configuration)",
733  runControlNumbers_.size(),
734  0.,
735  runControlNumbers_.size());
736 
737  for (const auto& run : runControlNumbers_ | boost::adaptors::indexed(1)) {
738  h_runFromConfig->SetBinContent(run.index(), run.value());
739  }
740 
741  // resolutions
742 
743  if (!checkBinOrdering(mypT_bins_)) {
744  edm::LogError("SplitVertexResolution") << " Warning - the vector of pT bins is not ordered " << std::endl;
745  }
746 
747  if (!checkBinOrdering(myNTrack_bins_)) {
748  edm::LogError("SplitVertexResolution") << " Warning -the vector of n. tracks bins is not ordered " << std::endl;
749  }
750 
751  if (!checkBinOrdering(myNVtx_bins_)) {
752  edm::LogError("SplitVertexResolution") << " Warning -the vector of n. vertices bins is not ordered " << std::endl;
753  }
754 
755  TFileDirectory xResolSumPt = outfile_->mkdir("xResolSumPt");
756  h_resolX_sumPt_ = bookResidualsHistogram(xResolSumPt, nPtBins_, "resolX", "sumPt");
757 
758  TFileDirectory yResolSumPt = outfile_->mkdir("yResolSumPt");
759  h_resolY_sumPt_ = bookResidualsHistogram(yResolSumPt, nPtBins_, "resolY", "sumPt");
760 
761  TFileDirectory zResolSumPt = outfile_->mkdir("zResolSumPt");
762  h_resolZ_sumPt_ = bookResidualsHistogram(zResolSumPt, nPtBins_, "resolZ", "sumPt");
763 
764  TFileDirectory xResolNtracks_ = outfile_->mkdir("xResolNtracks");
765  h_resolX_Ntracks_ = bookResidualsHistogram(xResolNtracks_, nTrackBins_, "resolX", "Ntracks");
766 
767  TFileDirectory yResolNtracks_ = outfile_->mkdir("yResolNtracks");
768  h_resolY_Ntracks_ = bookResidualsHistogram(yResolNtracks_, nTrackBins_, "resolY", "Ntracks");
769 
770  TFileDirectory zResolNtracks_ = outfile_->mkdir("zResolNtracks");
771  h_resolZ_Ntracks_ = bookResidualsHistogram(zResolNtracks_, nTrackBins_, "resolZ", "Ntracks");
772 
773  TFileDirectory xResolNvtx_ = outfile_->mkdir("xResolNvtx");
774  h_resolX_Nvtx_ = bookResidualsHistogram(xResolNvtx_, nVtxBins_, "resolX", "Nvtx");
775 
776  TFileDirectory yResolNvtx_ = outfile_->mkdir("yResolNvtx");
777  h_resolY_Nvtx_ = bookResidualsHistogram(yResolNvtx_, nVtxBins_, "resolY", "Nvtx");
778 
779  TFileDirectory zResolNvtx_ = outfile_->mkdir("zResolNvtx");
780  h_resolZ_Nvtx_ = bookResidualsHistogram(zResolNvtx_, nVtxBins_, "resolZ", "Nvtx");
781 
782  // pulls
783 
784  TFileDirectory xPullSumPt = outfile_->mkdir("xPullSumPt");
785  h_pullX_sumPt_ = bookResidualsHistogram(xPullSumPt, nPtBins_, "pullX", "sumPt");
786 
787  TFileDirectory yPullSumPt = outfile_->mkdir("yPullSumPt");
788  h_pullY_sumPt_ = bookResidualsHistogram(yPullSumPt, nPtBins_, "pullY", "sumPt");
789 
790  TFileDirectory zPullSumPt = outfile_->mkdir("zPullSumPt");
791  h_pullZ_sumPt_ = bookResidualsHistogram(zPullSumPt, nPtBins_, "pullZ", "sumPt");
792 
793  TFileDirectory xPullNtracks_ = outfile_->mkdir("xPullNtracks");
794  h_pullX_Ntracks_ = bookResidualsHistogram(xPullNtracks_, nTrackBins_, "pullX", "Ntracks");
795 
796  TFileDirectory yPullNtracks_ = outfile_->mkdir("yPullNtracks");
797  h_pullY_Ntracks_ = bookResidualsHistogram(yPullNtracks_, nTrackBins_, "pullY", "Ntracks");
798 
799  TFileDirectory zPullNtracks_ = outfile_->mkdir("zPullNtracks");
800  h_pullZ_Ntracks_ = bookResidualsHistogram(zPullNtracks_, nTrackBins_, "pullZ", "Ntracks");
801 
802  TFileDirectory xPullNvtx_ = outfile_->mkdir("xPullNvtx");
803  h_pullX_Nvtx_ = bookResidualsHistogram(xPullNvtx_, nVtxBins_, "pullX", "Nvtx");
804 
805  TFileDirectory yPullNvtx_ = outfile_->mkdir("yPullNvtx");
806  h_pullY_Nvtx_ = bookResidualsHistogram(yPullNvtx_, nVtxBins_, "pullY", "Nvtx");
807 
808  TFileDirectory zPullNvtx_ = outfile_->mkdir("zPullNvtx");
809  h_pullZ_Nvtx_ = bookResidualsHistogram(zPullNvtx_, nVtxBins_, "pullZ", "Nvtx");
810 
811  // control plots
812  h_runNumber = outfile_->make<TH1F>("h_runNumber", "run number;run number;n_{events}", 100000, 250000., 350000.);
813 
814  h_nOfflineVertices = outfile_->make<TH1I>("h_nOfflineVertices", "n. of vertices;n. vertices; events", 100, 0, 100);
815  h_nVertices = outfile_->make<TH1I>("h_nVertices", "n. of vertices;n. vertices; events", 100, 0, 100);
816  h_nNonFakeVertices =
817  outfile_->make<TH1I>("h_nRealVertices", "n. of non-fake vertices;n. vertices; events", 100, 0, 100);
818  h_nFinalVertices =
819  outfile_->make<TH1I>("h_nSelectedVertices", "n. of selected vertices vertices;n. vertices; events", 100, 0, 100);
820 
821  h_diffX = outfile_->make<TH1F>(
822  "h_diffX", "x-coordinate vertex resolution;vertex resolution (x) [#mum];vertices", 100, -300, 300.);
823  h_diffY = outfile_->make<TH1F>(
824  "h_diffY", "y-coordinate vertex resolution;vertex resolution (y) [#mum];vertices", 100, -300, 300.);
825  h_diffZ = outfile_->make<TH1F>(
826  "h_diffZ", "z-coordinate vertex resolution;vertex resolution (z) [#mum];vertices", 100, -500, 500.);
827 
828  h_OrigVertexErrX = outfile_->make<TH1F>(
829  "h_OrigVertexErrX", "x-coordinate vertex error;vertex error (x) [#mum];vertices", 300, 0., 300.);
830  h_OrigVertexErrY = outfile_->make<TH1F>(
831  "h_OrigVertexErrY", "y-coordinate vertex error;vertex error (y) [#mum];vertices", 300, 0., 300.);
832  h_OrigVertexErrZ = outfile_->make<TH1F>(
833  "h_OrigVertexErrZ", "z-coordinate vertex error;vertex error (z) [#mum];vertices", 500, 0., 500.);
834 
835  h_errX = outfile_->make<TH1F>(
836  "h_errX", "x-coordinate vertex resolution error;vertex resoltuion error (x) [#mum];vertices", 300, 0., 300.);
837  h_errY = outfile_->make<TH1F>(
838  "h_errY", "y-coordinate vertex resolution error;vertex resolution error (y) [#mum];vertices", 300, 0., 300.);
839  h_errZ = outfile_->make<TH1F>(
840  "h_errZ", "z-coordinate vertex resolution error;vertex resolution error (z) [#mum];vertices", 500, 0., 500.);
841 
842  h_pullX = outfile_->make<TH1F>("h_pullX", "x-coordinate vertex pull;vertex pull (x);vertices", 500, -10, 10.);
843  h_pullY = outfile_->make<TH1F>("h_pullY", "y-coordinate vertex pull;vertex pull (y);vertices", 500, -10, 10.);
844  h_pullZ = outfile_->make<TH1F>("h_pullZ", "z-coordinate vertex pull;vertex pull (z);vertices", 500, -10, 10.);
845 
846  h_ntrks = outfile_->make<TH1F>("h_ntrks",
847  "number of tracks in vertex;vertex multeplicity;vertices",
848  myNTrack_bins_.size() - 1,
849  myNTrack_bins_.data());
850 
851  h_sumPt = outfile_->make<TH1F>(
852  "h_sumPt", "#Sigma p_{T};#sum p_{T} [GeV];vertices", mypT_bins_.size() - 1, mypT_bins_.data());
853 
854  h_avgSumPt = outfile_->make<TH1F>(
855  "h_avgSumPt", "#LT #Sigma p_{T} #GT;#LT #sum p_{T} #GT [GeV];vertices", mypT_bins_.size() - 1, mypT_bins_.data());
856 
857  h_sumPt1 = outfile_->make<TH1F>("h_sumPt1",
858  "#Sigma p_{T} sub-vertex 1;#sum p_{T} sub-vertex 1 [GeV];subvertices",
859  mypT_bins_.size() - 1,
860  mypT_bins_.data());
861  h_sumPt2 = outfile_->make<TH1F>("h_sumPt2",
862  "#Sigma p_{T} sub-vertex 2;#sum p_{T} sub-vertex 2 [GeV];subvertices",
863  mypT_bins_.size() - 1,
864  mypT_bins_.data());
865 
866  h_wTrks1 = outfile_->make<TH1F>("h_wTrks1", "weight of track for sub-vertex 1;track weight;subvertices", 500, 0., 1.);
867  h_wTrks2 = outfile_->make<TH1F>("h_wTrks2", "weithg of track for sub-vertex 2;track weight;subvertices", 500, 0., 1.);
868 
869  h_minWTrks1 = outfile_->make<TH1F>(
870  "h_minWTrks1", "minimum weight of track for sub-vertex 1;minimum track weight;subvertices", 500, 0., 1.);
871  h_minWTrks2 = outfile_->make<TH1F>(
872  "h_minWTrks2", "minimum weithg of track for sub-vertex 2;minimum track weight;subvertices", 500, 0., 1.);
873 
874  h_PVCL_subVtx1 =
875  outfile_->make<TH1F>("h_PVCL_subVtx1",
876  "#chi^{2} probability for sub-vertex 1;Prob(#chi^{2},ndof) for sub-vertex 1;subvertices",
877  100,
878  0.,
879  1);
880  h_PVCL_subVtx2 =
881  outfile_->make<TH1F>("h_PVCL_subVtx2",
882  "#chi^{2} probability for sub-vertex 2;Prob(#chi^{2},ndof) for sub-vertex 2;subvertices",
883  100,
884  0.,
885  1);
886 
887  // resolutions
888 
889  p_resolX_vsSumPt = outfile_->make<TH1F>("p_resolX_vsSumPt",
890  "x-resolution vs #Sigma p_{T};#sum p_{T} [GeV]; x vertex resolution [#mum]",
891  mypT_bins_.size() - 1,
892  mypT_bins_.data());
893  p_resolY_vsSumPt = outfile_->make<TH1F>("p_resolY_vsSumPt",
894  "y-resolution vs #Sigma p_{T};#sum p_{T} [GeV]; y vertex resolution [#mum]",
895  mypT_bins_.size() - 1,
896  mypT_bins_.data());
897  p_resolZ_vsSumPt = outfile_->make<TH1F>("p_resolZ_vsSumPt",
898  "z-resolution vs #Sigma p_{T};#sum p_{T} [GeV]; z vertex resolution [#mum]",
899  mypT_bins_.size() - 1,
900  mypT_bins_.data());
901 
902  p_resolX_vsNtracks = outfile_->make<TH1F>("p_resolX_vsNtracks",
903  "x-resolution vs n_{tracks};n_{tracks}; x vertex resolution [#mum]",
904  myNTrack_bins_.size() - 1,
905  myNTrack_bins_.data());
906  p_resolY_vsNtracks = outfile_->make<TH1F>("p_resolY_vsNtracks",
907  "y-resolution vs n_{tracks};n_{tracks}; y vertex resolution [#mum]",
908  myNTrack_bins_.size() - 1,
909  myNTrack_bins_.data());
910  p_resolZ_vsNtracks = outfile_->make<TH1F>("p_resolZ_vsNtracks",
911  "z-resolution vs n_{tracks};n_{tracks}; z vertex resolution [#mum]",
912  myNTrack_bins_.size() - 1,
913  myNTrack_bins_.data());
914 
915  p_resolX_vsNvtx = outfile_->make<TH1F>("p_resolX_vsNvtx",
916  "x-resolution vs n_{vertices};n_{vertices}; x vertex resolution [#mum]",
917  myNVtx_bins_.size() - 1,
918  myNVtx_bins_.data());
919  p_resolY_vsNvtx = outfile_->make<TH1F>("p_resolY_vsNvtx",
920  "y-resolution vs n_{vertices};n_{vertices}; y vertex resolution [#mum]",
921  myNVtx_bins_.size() - 1,
922  myNVtx_bins_.data());
923  p_resolZ_vsNvtx = outfile_->make<TH1F>("p_resolZ_vsNvtx",
924  "z-resolution vs n_{vertices};n_{vertices}; z vertex resolution [#mum]",
925  myNVtx_bins_.size() - 1,
926  myNVtx_bins_.data());
927 
928  // pulls
929 
930  p_pullX_vsSumPt = outfile_->make<TH1F>("p_pullX_vsSumPt",
931  "x-pull vs #Sigma p_{T};#sum p_{T} [GeV]; x vertex pull",
932  mypT_bins_.size() - 1,
933  mypT_bins_.data());
934  p_pullY_vsSumPt = outfile_->make<TH1F>("p_pullY_vsSumPt",
935  "y-pull vs #Sigma p_{T};#sum p_{T} [GeV]; y vertex pull",
936  mypT_bins_.size() - 1,
937  mypT_bins_.data());
938  p_pullZ_vsSumPt = outfile_->make<TH1F>("p_pullZ_vsSumPt",
939  "z-pull vs #Sigma p_{T};#sum p_{T} [GeV]; z vertex pull",
940  mypT_bins_.size() - 1,
941  mypT_bins_.data());
942 
943  p_pullX_vsNtracks = outfile_->make<TH1F>("p_pullX_vsNtracks",
944  "x-pull vs n_{tracks};n_{tracks}; x vertex pull",
945  myNTrack_bins_.size() - 1,
946  myNTrack_bins_.data());
947  p_pullY_vsNtracks = outfile_->make<TH1F>("p_pullY_vsNtracks",
948  "y-pull vs n_{tracks};n_{tracks}; y vertex pull",
949  myNTrack_bins_.size() - 1,
950  myNTrack_bins_.data());
951  p_pullZ_vsNtracks = outfile_->make<TH1F>("p_pullZ_vsNtracks",
952  "z-pull vs n_{tracks};n_{tracks}; z vertex pull",
953  myNTrack_bins_.size() - 1,
954  myNTrack_bins_.data());
955 
956  p_pullX_vsNvtx = outfile_->make<TH1F>("p_pullX_vsNvtx",
957  "x-pull vs n_{vertices};n_{vertices}; x vertex pull",
958  myNVtx_bins_.size() - 1,
959  myNVtx_bins_.data());
960  p_pullY_vsNvtx = outfile_->make<TH1F>("p_pullY_vsNvtx",
961  "y-pull vs n_{vertices};n_{vertices}; y vertex pull",
962  myNVtx_bins_.size() - 1,
963  myNVtx_bins_.data());
964  p_pullZ_vsNvtx = outfile_->make<TH1F>("p_pullZ_vsNvtx",
965  "z-pull vs n_{vertices};n_{vertices}; z vertex pull",
966  myNVtx_bins_.size() - 1,
967  myNVtx_bins_.data());
968 
969  tree_ = outfile_->make<TTree>("pvTree", "pvTree");
970  tree_->Branch("event", &event_, 64000, 2);
971 }
SplitVertexResolution::h_resolZ_Nvtx_
std::vector< TH1F * > h_resolZ_Nvtx_
Definition: SplitVertexResolution.cc:207
SplitVertexResolution::h_pullY_sumPt_
std::vector< TH1F * > h_pullY_sumPt_
Definition: SplitVertexResolution.cc:223
SplitVertexResolution::h_pullZ_sumPt_
std::vector< TH1F * > h_pullZ_sumPt_
Definition: SplitVertexResolution.cc:224
SplitVertexResolution::h_resolX_sumPt_
std::vector< TH1F * > h_resolX_sumPt_
Definition: SplitVertexResolution.cc:197
TFileService::make
T * make(const Args &...args) const
make new ROOT object
Definition: TFileService.h:64
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:123
SplitVertexResolution::h_pullY_Ntracks_
std::vector< TH1F * > h_pullY_Ntracks_
Definition: SplitVertexResolution.cc:227
SplitVertexResolution::mypT_bins_
std::array< float, nPtBins_+1 > mypT_bins_
Definition: SplitVertexResolution.cc:257
TFileDirectory::make
T * make(const Args &...args) const
make new ROOT object
Definition: TFileDirectory.h:53
SplitVertexResolution::myNVtx_bins_
std::array< float, nVtxBins_+1 > myNVtx_bins_
Definition: SplitVertexResolution.cc:265
SplitVertexResolution::h_resolX_Nvtx_
std::vector< TH1F * > h_resolX_Nvtx_
Definition: SplitVertexResolution.cc:205
SplitVertexResolution::outfile_
edm::Service< TFileService > outfile_
Definition: SplitVertexResolution.cc:142
TFileService::file
TFile & file() const
return opened TFile
Definition: TFileService.h:37
SplitVertexResolution::myNTrack_bins_
std::array< float, nTrackBins_+1 > myNTrack_bins_
Definition: SplitVertexResolution.cc:261
SplitVertexResolution::h_pullY_Nvtx_
std::vector< TH1F * > h_pullY_Nvtx_
Definition: SplitVertexResolution.cc:231
SplitVertexResolution::h_resolZ_Ntracks_
std::vector< TH1F * > h_resolZ_Ntracks_
Definition: SplitVertexResolution.cc:203
TFileService::mkdir
TFileDirectory mkdir(const std::string &dir, const std::string &descr="")
create a new subdirectory
Definition: TFileService.h:69
SplitVertexResolution::h_pullX_Nvtx_
std::vector< TH1F * > h_pullX_Nvtx_
Definition: SplitVertexResolution.cc:230
SplitVertexResolution::h_pullZ_Nvtx_
std::vector< TH1F * > h_pullZ_Nvtx_
Definition: SplitVertexResolution.cc:232
SplitVertexResolution::bookResidualsHistogram
std::vector< TH1F * > bookResidualsHistogram(TFileDirectory dir, unsigned int theNOfBins, TString resType, TString varType)
Definition: SplitVertexResolution.cc:976
SplitVertexResolution::h_resolY_Nvtx_
std::vector< TH1F * > h_resolY_Nvtx_
Definition: SplitVertexResolution.cc:206
SplitVertexResolution::h_pullX_sumPt_
std::vector< TH1F * > h_pullX_sumPt_
Definition: SplitVertexResolution.cc:222
SplitVertexResolution::h_pullX_Ntracks_
std::vector< TH1F * > h_pullX_Ntracks_
Definition: SplitVertexResolution.cc:226
SplitVertexResolution::checkBinOrdering
bool checkBinOrdering(std::array< float, SIZE > &bins)
Definition: SplitVertexResolution.cc:1260
SplitVertexResolution::h_pullZ_Ntracks_
std::vector< TH1F * > h_pullZ_Ntracks_
Definition: SplitVertexResolution.cc:228
SplitVertexResolution::h_resolY_sumPt_
std::vector< TH1F * > h_resolY_sumPt_
Definition: SplitVertexResolution.cc:198
SplitVertexResolution::runControlNumbers_
std::vector< unsigned int > runControlNumbers_
Definition: SplitVertexResolution.cc:138
SplitVertexResolution::h_resolX_Ntracks_
std::vector< TH1F * > h_resolX_Ntracks_
Definition: SplitVertexResolution.cc:201
SplitVertexResolution::h_resolY_Ntracks_
std::vector< TH1F * > h_resolY_Ntracks_
Definition: SplitVertexResolution.cc:202
SplitVertexResolution::h_resolZ_sumPt_
std::vector< TH1F * > h_resolZ_sumPt_
Definition: SplitVertexResolution.cc:199
void SplitVertexResolution::beginRun ( edm::Run const &  iEvent,
edm::EventSetup const &  iSetup 
)
overrideprivate

Definition at line 682 of file SplitVertexResolution.cc.

References debug_, getRunTime(), edm::RunBase::run(), and runNumbersTimesLog_.

682  {
683  unsigned int RunNumber_ = run.run();
684 
685  if (!runNumbersTimesLog_.count(RunNumber_)) {
686  auto times = getRunTime(iSetup);
687 
688  if (debug_) {
689  const time_t start_time = times.first / 1.0e+6;
690  edm::LogInfo("SplitVertexResolution")
691  << RunNumber_ << " has start time: " << times.first << " - " << times.second << std::endl;
692  edm::LogInfo("SplitVertexResolution")
693  << "human readable time: " << std::asctime(std::gmtime(&start_time)) << std::endl;
694  }
695  runNumbersTimesLog_[RunNumber_] = times;
696  }
697 }
SplitVertexResolution::getRunTime
std::pair< long long, long long > getRunTime(const edm::EventSetup &iSetup) const
Definition: SplitVertexResolution.cc:1114
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
SplitVertexResolution::runNumbersTimesLog_
std::map< unsigned int, std::pair< long long, long long > > runNumbersTimesLog_
Definition: SplitVertexResolution.cc:147
std::vector< TH1F * > SplitVertexResolution::bookResidualsHistogram ( TFileDirectory  dir,
unsigned int  theNOfBins,
TString  resType,
TString  varType 
)
private

Definition at line 976 of file SplitVertexResolution.cc.

References down, h, mps_fire::i, TFileDirectory::make(), and up.

Referenced by beginJob().

979  {
980  TH1F::SetDefaultSumw2(kTRUE);
981 
982  double up = 500.;
983  double down = -500.;
984 
985  if (resType.Contains("pull")) {
986  up *= 0.01;
987  down *= 0.01;
988  }
989 
990  std::vector<TH1F*> h;
991  h.reserve(theNOfBins);
992 
993  const char* auxResType = (resType.ReplaceAll("_", "")).Data();
994 
995  for (unsigned int i = 0; i < theNOfBins; i++) {
996  TH1F* htemp = dir.make<TH1F>(Form("histo_%s_%s_plot%i", resType.Data(), varType.Data(), i),
997  Form("%s vs %s - bin %i;%s;vertices", auxResType, varType.Data(), i, auxResType),
998  250,
999  down,
1000  up);
1001  h.push_back(htemp);
1002  }
1003 
1004  return h;
1005 }
up
Definition: BitonicSort.h:7
mps_fire.i
i
Definition: mps_fire.py:428
dumpRecoGeometry_cfg.varType
tuple varType
Definition: dumpRecoGeometry_cfg.py:8
TFileDirectory::make
T * make(const Args &...args) const
make new ROOT object
Definition: TFileDirectory.h:53
h
The Signals That Services Can Subscribe To This is based on ActivityRegistry h
Helper function to determine trigger accepts.
Definition: Activities.doc:4
template<std::size_t SIZE>
bool SplitVertexResolution::checkBinOrdering ( std::array< float, SIZE > &  bins)
private

Definition at line 1260 of file SplitVertexResolution.cc.

References newFWLiteAna::bin, and mps_fire::i.

Referenced by beginJob().

1262 {
1263  int i = 1;
1264 
1265  if (std::is_sorted(bins.begin(), bins.end())) {
1266  return true;
1267  } else {
1268  for (const auto& bin : bins) {
1269  edm::LogInfo("SplitVertexResolution") << "bin: " << i << " : " << bin << std::endl;
1270  i++;
1271  }
1272  edm::LogInfo("SplitVertexResolution") << "--------------------------------" << std::endl;
1273  return false;
1274  }
1275 }
mps_fire.i
i
Definition: mps_fire.py:428
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
void SplitVertexResolution::endJob ( void  )
overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 1008 of file SplitVertexResolution.cc.

References submitPVResolutionJobs::count, evtsByTrigger_, fillTrendPlotByIndex(), first, h_pullX_Ntracks_, h_pullX_Nvtx_, h_pullX_sumPt_, h_pullY_Ntracks_, h_pullY_Nvtx_, h_pullY_sumPt_, h_pullZ_Ntracks_, h_pullZ_Nvtx_, h_pullZ_sumPt_, h_resolX_Ntracks_, h_resolX_Nvtx_, h_resolX_sumPt_, h_resolY_Ntracks_, h_resolY_Nvtx_, h_resolY_sumPt_, h_resolZ_Ntracks_, h_resolZ_Nvtx_, h_resolZ_sumPt_, h_runEndTimes, h_runStartTimes, mps_fire::i, ievt, itrks, TFileDirectory::make(), TFileService::make(), TFileService::mkdir(), outfile_, p_pullX_vsNtracks, p_pullX_vsNvtx, p_pullX_vsSumPt, p_pullY_vsNtracks, p_pullY_vsNvtx, p_pullY_vsSumPt, p_pullZ_vsNtracks, p_pullZ_vsNvtx, p_pullZ_vsSumPt, p_resolX_vsNtracks, p_resolX_vsNvtx, p_resolX_vsSumPt, p_resolY_vsNtracks, p_resolY_vsNvtx, p_resolY_vsSumPt, p_resolZ_vsNtracks, p_resolZ_vsNvtx, p_resolZ_vsSumPt, submitPVValidationJobs::run, runNumbersTimesLog_, edm::second(), AlCaHLTBitMon_QueryRunRegistry::string, tksByTrigger_, cond::impl::to_string(), triggerMap_, and PVValHelper::WIDTH.

Referenced by o2olib.O2ORunMgr::executeJob().

1008  {
1009  edm::LogVerbatim("SplitVertexResolution") << "*******************************" << std::endl;
1010  edm::LogVerbatim("SplitVertexResolution") << "Events run in total: " << ievt << std::endl;
1011  edm::LogVerbatim("SplitVertexResolution") << "n. tracks: " << itrks << std::endl;
1012  edm::LogVerbatim("SplitVertexResolution") << "*******************************" << std::endl;
1013 
1014  int nFiringTriggers = triggerMap_.size();
1015  edm::LogVerbatim("SplitVertexResolution") << "firing triggers: " << nFiringTriggers << std::endl;
1016  edm::LogVerbatim("SplitVertexResolution") << "*******************************" << std::endl;
1017 
1018  tksByTrigger_ = outfile_->make<TH1D>(
1019  "tksByTrigger", "tracks by HLT path;;% of # traks", nFiringTriggers, -0.5, nFiringTriggers - 0.5);
1020  evtsByTrigger_ = outfile_->make<TH1D>(
1021  "evtsByTrigger", "events by HLT path;;% of # events", nFiringTriggers, -0.5, nFiringTriggers - 0.5);
1022 
1023  int i = 0;
1024  for (std::map<std::string, std::pair<int, int>>::iterator it = triggerMap_.begin(); it != triggerMap_.end(); ++it) {
1025  i++;
1026 
1027  double trkpercent = ((it->second).second) * 100. / double(itrks);
1028  double evtpercent = ((it->second).first) * 100. / double(ievt);
1029 
1030  edm::LogVerbatim("SplitVertexResolution")
1031  << "HLT path: " << std::setw(60) << std::left << it->first << " | events firing: " << std::right << std::setw(8)
1032  << (it->second).first << " (" << std::setw(8) << std::fixed << std::setprecision(4) << evtpercent << "%)"
1033  << " | tracks collected: " << std::setw(8) << (it->second).second << " (" << std::setw(8) << std::fixed
1034  << std::setprecision(4) << trkpercent << "%)";
1035 
1036  tksByTrigger_->SetBinContent(i, trkpercent);
1037  tksByTrigger_->GetXaxis()->SetBinLabel(i, (it->first).c_str());
1038 
1039  evtsByTrigger_->SetBinContent(i, evtpercent);
1040  evtsByTrigger_->GetXaxis()->SetBinLabel(i, (it->first).c_str());
1041  }
1042 
1043  TFileDirectory RunFeatures = outfile_->mkdir("RunFeatures");
1044  h_runStartTimes = RunFeatures.make<TH1I>(
1045  "runStartTimes", "run start times", runNumbersTimesLog_.size(), 0, runNumbersTimesLog_.size());
1046  h_runEndTimes =
1047  RunFeatures.make<TH1I>("runEndTimes", "run end times", runNumbersTimesLog_.size(), 0, runNumbersTimesLog_.size());
1048 
1049  unsigned int count = 1;
1050  for (const auto& run : runNumbersTimesLog_) {
1051  // strip down the microseconds
1052  h_runStartTimes->SetBinContent(count, run.second.first / 10e6);
1053  h_runStartTimes->GetXaxis()->SetBinLabel(count, (std::to_string(run.first)).c_str());
1054 
1055  h_runEndTimes->SetBinContent(count, run.second.second / 10e6);
1056  h_runEndTimes->GetXaxis()->SetBinLabel(count, (std::to_string(run.first)).c_str());
1057 
1058  count++;
1059  }
1060 
1061  // resolutions
1062 
1063  fillTrendPlotByIndex(p_resolX_vsSumPt, h_resolX_sumPt_, PVValHelper::WIDTH);
1064  fillTrendPlotByIndex(p_resolY_vsSumPt, h_resolY_sumPt_, PVValHelper::WIDTH);
1065  fillTrendPlotByIndex(p_resolZ_vsSumPt, h_resolZ_sumPt_, PVValHelper::WIDTH);
1066 
1067  fillTrendPlotByIndex(p_resolX_vsNtracks, h_resolX_Ntracks_, PVValHelper::WIDTH);
1068  fillTrendPlotByIndex(p_resolY_vsNtracks, h_resolY_Ntracks_, PVValHelper::WIDTH);
1069  fillTrendPlotByIndex(p_resolZ_vsNtracks, h_resolZ_Ntracks_, PVValHelper::WIDTH);
1070 
1071  fillTrendPlotByIndex(p_resolX_vsNvtx, h_resolX_Nvtx_, PVValHelper::WIDTH);
1072  fillTrendPlotByIndex(p_resolY_vsNvtx, h_resolY_Nvtx_, PVValHelper::WIDTH);
1073  fillTrendPlotByIndex(p_resolZ_vsNvtx, h_resolZ_Nvtx_, PVValHelper::WIDTH);
1074 
1075  // pulls
1076 
1077  fillTrendPlotByIndex(p_pullX_vsSumPt, h_pullX_sumPt_, PVValHelper::WIDTH);
1078  fillTrendPlotByIndex(p_pullY_vsSumPt, h_pullY_sumPt_, PVValHelper::WIDTH);
1079  fillTrendPlotByIndex(p_pullZ_vsSumPt, h_pullZ_sumPt_, PVValHelper::WIDTH);
1080 
1081  fillTrendPlotByIndex(p_pullX_vsNtracks, h_pullX_Ntracks_, PVValHelper::WIDTH);
1082  fillTrendPlotByIndex(p_pullY_vsNtracks, h_pullY_Ntracks_, PVValHelper::WIDTH);
1083  fillTrendPlotByIndex(p_pullZ_vsNtracks, h_pullZ_Ntracks_, PVValHelper::WIDTH);
1084 
1085  fillTrendPlotByIndex(p_pullX_vsNvtx, h_pullX_Nvtx_, PVValHelper::WIDTH);
1086  fillTrendPlotByIndex(p_pullY_vsNvtx, h_pullY_Nvtx_, PVValHelper::WIDTH);
1087  fillTrendPlotByIndex(p_pullZ_vsNvtx, h_pullZ_Nvtx_, PVValHelper::WIDTH);
1088 }
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
mps_fire.i
i
Definition: mps_fire.py:428
SplitVertexResolution::h_resolZ_Nvtx_
std::vector< TH1F * > h_resolZ_Nvtx_
Definition: SplitVertexResolution.cc:207
SplitVertexResolution::h_pullY_sumPt_
std::vector< TH1F * > h_pullY_sumPt_
Definition: SplitVertexResolution.cc:223
SplitVertexResolution::h_pullZ_sumPt_
std::vector< TH1F * > h_pullZ_sumPt_
Definition: SplitVertexResolution.cc:224
SplitVertexResolution::h_resolX_sumPt_
std::vector< TH1F * > h_resolX_sumPt_
Definition: SplitVertexResolution.cc:197
cond::impl::to_string
std::string to_string(const V &value)
Definition: OMSAccess.h:71
TFileService::make
T * make(const Args &...args) const
make new ROOT object
Definition: TFileService.h:64
SplitVertexResolution::fillTrendPlotByIndex
void fillTrendPlotByIndex(TH1F *trendPlot, std::vector< TH1F * > &h, PVValHelper::estimator fitPar_)
Definition: SplitVertexResolution.cc:1126
SplitVertexResolution::h_pullY_Ntracks_
std::vector< TH1F * > h_pullY_Ntracks_
Definition: SplitVertexResolution.cc:227
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
edm::second
U second(std::pair< T, U > const &p)
Definition: ParameterSet.cc:222
TFileDirectory::make
T * make(const Args &...args) const
make new ROOT object
Definition: TFileDirectory.h:53
SplitVertexResolution::h_resolX_Nvtx_
std::vector< TH1F * > h_resolX_Nvtx_
Definition: SplitVertexResolution.cc:205
SplitVertexResolution::outfile_
edm::Service< TFileService > outfile_
Definition: SplitVertexResolution.cc:142
SplitVertexResolution::triggerMap_
std::map< std::string, std::pair< int, int > > triggerMap_
Definition: SplitVertexResolution.cc:267
SplitVertexResolution::runNumbersTimesLog_
std::map< unsigned int, std::pair< long long, long long > > runNumbersTimesLog_
Definition: SplitVertexResolution.cc:147
SplitVertexResolution::h_pullY_Nvtx_
std::vector< TH1F * > h_pullY_Nvtx_
Definition: SplitVertexResolution.cc:231
SplitVertexResolution::h_resolZ_Ntracks_
std::vector< TH1F * > h_resolZ_Ntracks_
Definition: SplitVertexResolution.cc:203
TFileService::mkdir
TFileDirectory mkdir(const std::string &dir, const std::string &descr="")
create a new subdirectory
Definition: TFileService.h:69
SplitVertexResolution::h_pullX_Nvtx_
std::vector< TH1F * > h_pullX_Nvtx_
Definition: SplitVertexResolution.cc:230
SplitVertexResolution::h_pullZ_Nvtx_
std::vector< TH1F * > h_pullZ_Nvtx_
Definition: SplitVertexResolution.cc:232
SplitVertexResolution::h_resolY_Nvtx_
std::vector< TH1F * > h_resolY_Nvtx_
Definition: SplitVertexResolution.cc:206
SplitVertexResolution::h_pullX_sumPt_
std::vector< TH1F * > h_pullX_sumPt_
Definition: SplitVertexResolution.cc:222
SplitVertexResolution::h_pullX_Ntracks_
std::vector< TH1F * > h_pullX_Ntracks_
Definition: SplitVertexResolution.cc:226
SplitVertexResolution::h_pullZ_Ntracks_
std::vector< TH1F * > h_pullZ_Ntracks_
Definition: SplitVertexResolution.cc:228
SplitVertexResolution::h_resolY_sumPt_
std::vector< TH1F * > h_resolY_sumPt_
Definition: SplitVertexResolution.cc:198
SplitVertexResolution::h_resolX_Ntracks_
std::vector< TH1F * > h_resolX_Ntracks_
Definition: SplitVertexResolution.cc:201
SplitVertexResolution::h_resolY_Ntracks_
std::vector< TH1F * > h_resolY_Ntracks_
Definition: SplitVertexResolution.cc:202
SplitVertexResolution::h_resolZ_sumPt_
std::vector< TH1F * > h_resolZ_sumPt_
Definition: SplitVertexResolution.cc:199
void SplitVertexResolution::endRun ( edm::Run const &  ,
edm::EventSetup const &   
)
inlineoverrideprivate

Definition at line 90 of file SplitVertexResolution.cc.

90 {};
void SplitVertexResolution::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions)
static

Definition at line 1091 of file SplitVertexResolution.cc.

References edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addUntracked(), edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, HLT_FULL_cff::InputTag, and edm::ParameterSetDescription::setComment().

1091  {
1092  edm::ParameterSetDescription desc;
1093  desc.setComment(
1094  "Validates alignment payloads by evaluating the resulting Primary Vertex Resolution via vertex splitting method");
1095 
1096  desc.addUntracked<int>("compressionSettings", -1);
1097  desc.add<bool>("storeNtuple", false);
1098  desc.addUntracked<double>("intLumi", 0.);
1099  desc.addUntracked<bool>("Debug", false);
1100  desc.add<edm::InputTag>("vtxCollection", edm::InputTag("offlinePrimaryVertices"));
1101  desc.add<edm::InputTag>("trackCollection", edm::InputTag("generalTracks"));
1102  desc.addUntracked<double>("minVertexNdf", 10);
1103  desc.addUntracked<double>("minVertexMeanWeight", 0.5);
1104  desc.addUntracked<bool>("runControl", false);
1105  desc.addUntracked<std::vector<unsigned int>>("runControlNumber", {});
1106  desc.addUntracked<double>("sumpTStartScale", 1.);
1107  desc.addUntracked<double>("sumpTEndScale", 1e3);
1108  desc.addUntracked<double>("nTrackBins", 120.);
1109  desc.addUntracked<double>("nVtxBins", 60.);
1110  descriptions.addWithDefaultLabel(desc);
1111 }
edm::ConfigurationDescriptions::addWithDefaultLabel
void addWithDefaultLabel(ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:87
edm::ParameterSetDescription::addUntracked
ParameterDescriptionBase * addUntracked(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:100
edm::ParameterSetDescription::setComment
void setComment(std::string const &value)
Definition: ParameterSetDescription.cc:33
edm::ParameterSetDescription::add
ParameterDescriptionBase * add(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:95
HLT_FULL_cff.InputTag
tuple InputTag
Definition: HLT_FULL_cff.py:71427
void SplitVertexResolution::fillTrendPlotByIndex ( TH1F *  trendPlot,
std::vector< TH1F * > &  h,
PVValHelper::estimator  fitPar_ 
)
private

Definition at line 1126 of file SplitVertexResolution.cc.

References newFWLiteAna::bin, HLT_FULL_cff::distance, Measurement1D::error(), fitResiduals(), PVValHelper::getMAD(), PVValHelper::getMedian(), PVValHelper::MAD, PVValHelper::MEAN, PVValHelper::MEDIAN, Measurement1D::value(), and PVValHelper::WIDTH.

Referenced by endJob().

1128 {
1129  for (auto iterator = h.begin(); iterator != h.end(); iterator++) {
1130  unsigned int bin = std::distance(h.begin(), iterator) + 1;
1131  statmode::fitParams myFit = fitResiduals((*iterator));
1132 
1133  switch (fitPar_) {
1134  case PVValHelper::MEAN: {
1135  float mean_ = myFit.first.value();
1136  float meanErr_ = myFit.first.error();
1137  trendPlot->SetBinContent(bin, mean_);
1138  trendPlot->SetBinError(bin, meanErr_);
1139  break;
1140  }
1141  case PVValHelper::WIDTH: {
1142  float width_ = myFit.second.value();
1143  float widthErr_ = myFit.second.error();
1144  trendPlot->SetBinContent(bin, width_);
1145  trendPlot->SetBinError(bin, widthErr_);
1146  break;
1147  }
1148  case PVValHelper::MEDIAN: {
1149  float median_ = PVValHelper::getMedian((*iterator)).value();
1150  float medianErr_ = PVValHelper::getMedian((*iterator)).error();
1151  trendPlot->SetBinContent(bin, median_);
1152  trendPlot->SetBinError(bin, medianErr_);
1153  break;
1154  }
1155  case PVValHelper::MAD: {
1156  float mad_ = PVValHelper::getMAD((*iterator)).value();
1157  float madErr_ = PVValHelper::getMAD((*iterator)).error();
1158  trendPlot->SetBinContent(bin, mad_);
1159  trendPlot->SetBinError(bin, madErr_);
1160  break;
1161  }
1162  default:
1163  edm::LogWarning("SplitVertexResolution")
1164  << "fillTrendPlotByIndex() " << fitPar_ << " unknown estimator!" << std::endl;
1165  break;
1166  }
1167  }
1168 }
PVValHelper::getMedian
Measurement1D getMedian(TH1F *histo)
Definition: PVValidationHelpers.cc:180
HLT_FULL_cff.distance
tuple distance
Definition: HLT_FULL_cff.py:6880
Measurement1D::error
double error() const
Definition: Measurement1D.h:27
PVValHelper::getMAD
Measurement1D getMAD(TH1F *histo)
Definition: PVValidationHelpers.cc:196
statmode::fitParams
std::pair< Measurement1D, Measurement1D > fitParams
Definition: SplitVertexResolution.cc:69
SplitVertexResolution::fitResiduals
statmode::fitParams fitResiduals(TH1 *hist, bool singleTime=false)
Definition: SplitVertexResolution.cc:1171
Measurement1D::value
double value() const
Definition: Measurement1D.h:25
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:122
h
The Signals That Services Can Subscribe To This is based on ActivityRegistry h
Helper function to determine trigger accepts.
Definition: Activities.doc:4
statmode::fitParams SplitVertexResolution::fitResiduals ( TH1 *  hist,
bool  singleTime = false 
)
private

Definition at line 1171 of file SplitVertexResolution.cc.

References cms::cuda::func, edm::isNotFinite(), LogDebug, SiStripPI::max, SiStripPI::mean, SiStripPI::min, and mps_fire::result.

Referenced by fillTrendPlotByIndex().

1173 {
1174  if (hist->GetEntries() < 10) {
1175  LogDebug("SplitVertexResolution") << "hist name: " << hist->GetName() << " has less than 10 entries" << std::endl;
1176  return std::make_pair(Measurement1D(0., 0.), Measurement1D(0., 0.));
1177  }
1178 
1179  float maxHist = hist->GetXaxis()->GetXmax();
1180  float minHist = hist->GetXaxis()->GetXmin();
1181  float mean = hist->GetMean();
1182  float sigma = hist->GetRMS();
1183 
1184  if (edm::isNotFinite(mean) || edm::isNotFinite(sigma)) {
1185  mean = 0;
1186  //sigma= - hist->GetXaxis()->GetBinLowEdge(1) + hist->GetXaxis()->GetBinLowEdge(hist->GetNbinsX()+1);
1187  sigma = -minHist + maxHist;
1188  edm::LogWarning("SplitVertexResolution")
1189  << "FitPVResiduals::fitResiduals(): histogram" << hist->GetName() << " mean or sigma are NaN!!" << std::endl;
1190  }
1191 
1192  TF1 func("tmp", "gaus", mean - 2. * sigma, mean + 2. * sigma);
1193  if (0 == hist->Fit(&func, "QNR")) { // N: do not blow up file by storing fit!
1194  mean = func.GetParameter(1);
1195  sigma = func.GetParameter(2);
1196 
1197  if (!singleTime) {
1198  // second fit: three sigma of first fit around mean of first fit
1199  func.SetRange(std::max(mean - 3 * sigma, minHist), std::min(mean + 3 * sigma, maxHist));
1200  // I: integral gives more correct results if binning is too wide
1201  // L: Likelihood can treat empty bins correctly (if hist not weighted...)
1202  if (0 == hist->Fit(&func, "Q0LR")) {
1203  if (hist->GetFunction(func.GetName())) { // Take care that it is later on drawn:
1204  hist->GetFunction(func.GetName())->ResetBit(TF1::kNotDraw);
1205  }
1206  }
1207  }
1208  }
1209 
1210  float res_mean = func.GetParameter(1);
1211  float res_width = func.GetParameter(2);
1212 
1213  float res_mean_err = func.GetParError(1);
1214  float res_width_err = func.GetParError(2);
1215 
1216  Measurement1D resultM(res_mean, res_mean_err);
1217  Measurement1D resultW(res_width, res_width_err);
1218 
1219  statmode::fitParams result = std::make_pair(resultM, resultW);
1220  return result;
1221 }
edm::isNotFinite
constexpr bool isNotFinite(T x)
Definition: isFinite.h:9
cms::cuda::func
uint32_t T const *__restrict__ uint32_t const *__restrict__ int32_t int Histo::index_type cudaStream_t Func __host__ __device__ V int Func func
Definition: HistoContainer.h:73
statmode::fitParams
std::pair< Measurement1D, Measurement1D > fitParams
Definition: SplitVertexResolution.cc:69
mps_fire.result
tuple result
Definition: mps_fire.py:311
gpuVertexFinder::hist
__shared__ Hist hist
Definition: gpuClusterTracksDBSCAN.h:48
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:122
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:233
statmode::fitParams SplitVertexResolution::fitResiduals_v0 ( TH1 *  hist)
private

Definition at line 1224 of file SplitVertexResolution.cc.

References cms::cuda::func, SiStripPI::mean, and mps_fire::result.

1226 {
1227  float mean = hist->GetMean();
1228  float sigma = hist->GetRMS();
1229 
1230  TF1 func("tmp", "gaus", mean - 1.5 * sigma, mean + 1.5 * sigma);
1231  if (0 == hist->Fit(&func, "QNR")) { // N: do not blow up file by storing fit!
1232  mean = func.GetParameter(1);
1233  sigma = func.GetParameter(2);
1234  // second fit: three sigma of first fit around mean of first fit
1235  func.SetRange(mean - 2 * sigma, mean + 2 * sigma);
1236  // I: integral gives more correct results if binning is too wide
1237  // L: Likelihood can treat empty bins correctly (if hist not weighted...)
1238  if (0 == hist->Fit(&func, "Q0LR")) {
1239  if (hist->GetFunction(func.GetName())) { // Take care that it is later on drawn:
1240  hist->GetFunction(func.GetName())->ResetBit(TF1::kNotDraw);
1241  }
1242  }
1243  }
1244 
1245  float res_mean = func.GetParameter(1);
1246  float res_width = func.GetParameter(2);
1247 
1248  float res_mean_err = func.GetParError(1);
1249  float res_width_err = func.GetParError(2);
1250 
1251  Measurement1D resultM(res_mean, res_mean_err);
1252  Measurement1D resultW(res_width, res_width_err);
1253 
1254  statmode::fitParams result = std::make_pair(resultM, resultW);
1255  return result;
1256 }
cms::cuda::func
uint32_t T const *__restrict__ uint32_t const *__restrict__ int32_t int Histo::index_type cudaStream_t Func __host__ __device__ V int Func func
Definition: HistoContainer.h:73
statmode::fitParams
std::pair< Measurement1D, Measurement1D > fitParams
Definition: SplitVertexResolution.cc:69
mps_fire.result
tuple result
Definition: mps_fire.py:311
gpuVertexFinder::hist
__shared__ Hist hist
Definition: gpuClusterTracksDBSCAN.h:48
std::pair< long long, long long > SplitVertexResolution::getRunTime ( const edm::EventSetup &  iSetup) const
private

Definition at line 1114 of file SplitVertexResolution.cc.

References debug_, submitPVValidationJobs::runInfo, and runInfoToken_.

Referenced by beginRun().

1116 {
1117  const auto& runInfo = iSetup.getData(runInfoToken_);
1118  if (debug_) {
1119  edm::LogInfo("SplitVertexResolution")
1120  << "start time: " << runInfo.m_start_time_str << " - stop time: " << runInfo.m_stop_time_str << std::endl;
1121  }
1122  return std::make_pair(runInfo.m_start_time_ll, runInfo.m_stop_time_ll);
1123 }
SplitVertexResolution::runInfoToken_
edm::ESGetToken< RunInfo, RunInfoRcd > runInfoToken_
Definition: SplitVertexResolution.cc:130
edm::EventSetup::getData
bool getData(T &iHolder) const
Definition: EventSetup.h:122
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
static bool SplitVertexResolution::mysorter ( reco::Track  i,
reco::Track  j 
)
inlinestatic

Definition at line 82 of file SplitVertexResolution.cc.

References reco::TrackBase::pt().

Referenced by analyze().

82 { return (i.pt() > j.pt()); }
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:637

Member Data Documentation

constexpr double SplitVertexResolution::cmToUm = 10000.
staticprivate

Definition at line 140 of file SplitVertexResolution.cc.

Referenced by analyze().

const int SplitVertexResolution::compressionSettings_
private

Definition at line 110 of file SplitVertexResolution.cc.

Referenced by beginJob().

bool SplitVertexResolution::debug_
private

Definition at line 117 of file SplitVertexResolution.cc.

Referenced by analyze(), beginRun(), getRunTime(), and SplitVertexResolution().

std::mt19937 SplitVertexResolution::engine_
private

Definition at line 248 of file SplitVertexResolution.cc.

Referenced by analyze().

pvEvent SplitVertexResolution::event_
private

Definition at line 250 of file SplitVertexResolution.cc.

Referenced by analyze(), beginEvent(), and beginJob().

TH1D* SplitVertexResolution::evtsByTrigger_
private

Definition at line 193 of file SplitVertexResolution.cc.

Referenced by endJob().

TH1F* SplitVertexResolution::h_avgSumPt
private

Definition at line 169 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

TH1F* SplitVertexResolution::h_diffX
private

Definition at line 151 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

TH1F* SplitVertexResolution::h_diffY
private

Definition at line 152 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

TH1F* SplitVertexResolution::h_diffZ
private

Definition at line 153 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

TH1F* SplitVertexResolution::h_errX
private

Definition at line 159 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

TH1F* SplitVertexResolution::h_errY
private

Definition at line 160 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

TH1F* SplitVertexResolution::h_errZ
private

Definition at line 161 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

TH1F* SplitVertexResolution::h_lumiFromConfig
private

Definition at line 144 of file SplitVertexResolution.cc.

Referenced by beginJob().

TH1F* SplitVertexResolution::h_minWTrks1
private

Definition at line 177 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

TH1F* SplitVertexResolution::h_minWTrks2
private

Definition at line 178 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

TH1I* SplitVertexResolution::h_nFinalVertices
private

Definition at line 188 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

TH1I* SplitVertexResolution::h_nNonFakeVertices
private

Definition at line 187 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

TH1I* SplitVertexResolution::h_nOfflineVertices
private

Definition at line 185 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

TH1F* SplitVertexResolution::h_ntrks
private

Definition at line 167 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

TH1I* SplitVertexResolution::h_nVertices
private

Definition at line 186 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

TH1F* SplitVertexResolution::h_OrigVertexErrX
private

Definition at line 155 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

TH1F* SplitVertexResolution::h_OrigVertexErrY
private

Definition at line 156 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

TH1F* SplitVertexResolution::h_OrigVertexErrZ
private

Definition at line 157 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

TH1F* SplitVertexResolution::h_profileBinnings
private

Definition at line 246 of file SplitVertexResolution.cc.

Referenced by beginJob().

TH1F* SplitVertexResolution::h_pullX
private

Definition at line 163 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

std::vector<TH1F*> SplitVertexResolution::h_pullX_Ntracks_
private

Definition at line 226 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and endJob().

std::vector<TH1F*> SplitVertexResolution::h_pullX_Nvtx_
private

Definition at line 230 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and endJob().

std::vector<TH1F*> SplitVertexResolution::h_pullX_sumPt_
private

Definition at line 222 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and endJob().

TH1F* SplitVertexResolution::h_pullY
private

Definition at line 164 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

std::vector<TH1F*> SplitVertexResolution::h_pullY_Ntracks_
private

Definition at line 227 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and endJob().

std::vector<TH1F*> SplitVertexResolution::h_pullY_Nvtx_
private

Definition at line 231 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and endJob().

std::vector<TH1F*> SplitVertexResolution::h_pullY_sumPt_
private

Definition at line 223 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and endJob().

TH1F* SplitVertexResolution::h_pullZ
private

Definition at line 165 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

std::vector<TH1F*> SplitVertexResolution::h_pullZ_Ntracks_
private

Definition at line 228 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and endJob().

std::vector<TH1F*> SplitVertexResolution::h_pullZ_Nvtx_
private

Definition at line 232 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and endJob().

std::vector<TH1F*> SplitVertexResolution::h_pullZ_sumPt_
private

Definition at line 224 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and endJob().

TH1F* SplitVertexResolution::h_PVCL_subVtx1
private

Definition at line 180 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

TH1F* SplitVertexResolution::h_PVCL_subVtx2
private

Definition at line 181 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

std::vector<TH1F*> SplitVertexResolution::h_resolX_Ntracks_
private

Definition at line 201 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and endJob().

std::vector<TH1F*> SplitVertexResolution::h_resolX_Nvtx_
private

Definition at line 205 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and endJob().

std::vector<TH1F*> SplitVertexResolution::h_resolX_sumPt_
private

Definition at line 197 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and endJob().

std::vector<TH1F*> SplitVertexResolution::h_resolY_Ntracks_
private

Definition at line 202 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and endJob().

std::vector<TH1F*> SplitVertexResolution::h_resolY_Nvtx_
private

Definition at line 206 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and endJob().

std::vector<TH1F*> SplitVertexResolution::h_resolY_sumPt_
private

Definition at line 198 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and endJob().

std::vector<TH1F*> SplitVertexResolution::h_resolZ_Ntracks_
private

Definition at line 203 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and endJob().

std::vector<TH1F*> SplitVertexResolution::h_resolZ_Nvtx_
private

Definition at line 207 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and endJob().

std::vector<TH1F*> SplitVertexResolution::h_resolZ_sumPt_
private

Definition at line 199 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and endJob().

TH1I* SplitVertexResolution::h_runEndTimes
private

Definition at line 149 of file SplitVertexResolution.cc.

Referenced by endJob().

TH1I* SplitVertexResolution::h_runFromConfig
private

Definition at line 145 of file SplitVertexResolution.cc.

Referenced by beginJob().

TH1F* SplitVertexResolution::h_runNumber
private

Definition at line 183 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

TH1I* SplitVertexResolution::h_runStartTimes
private

Definition at line 148 of file SplitVertexResolution.cc.

Referenced by endJob().

TH1F* SplitVertexResolution::h_sumPt
private

Definition at line 168 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

TH1F* SplitVertexResolution::h_sumPt1
private

Definition at line 171 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

TH1F* SplitVertexResolution::h_sumPt2
private

Definition at line 172 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

TH1F* SplitVertexResolution::h_wTrks1
private

Definition at line 174 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

TH1F* SplitVertexResolution::h_wTrks2
private

Definition at line 175 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

int SplitVertexResolution::ievt
private

Definition at line 106 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and endJob().

double SplitVertexResolution::intLumi_
private

Definition at line 116 of file SplitVertexResolution.cc.

Referenced by beginJob().

int SplitVertexResolution::itrks
private

Definition at line 107 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and endJob().

double SplitVertexResolution::minVtxNdf_
private

Definition at line 132 of file SplitVertexResolution.cc.

Referenced by analyze().

double SplitVertexResolution::minVtxWgt_
private

Definition at line 133 of file SplitVertexResolution.cc.

Referenced by analyze().

std::array<float, nTrackBins_ + 1> SplitVertexResolution::myNTrack_bins_
private

Definition at line 261 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and SplitVertexResolution().

std::array<float, nVtxBins_ + 1> SplitVertexResolution::myNVtx_bins_
private

Definition at line 265 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and SplitVertexResolution().

std::array<float, nPtBins_ + 1> SplitVertexResolution::mypT_bins_
private

Definition at line 257 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and SplitVertexResolution().

const int SplitVertexResolution::nPtBins_ = 30
staticprivate

Definition at line 256 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and SplitVertexResolution().

const int SplitVertexResolution::nTrackBins_ = 120
staticprivate

Definition at line 259 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and SplitVertexResolution().

const double SplitVertexResolution::nVisTrackBins_
private

Definition at line 260 of file SplitVertexResolution.cc.

Referenced by beginJob().

const double SplitVertexResolution::nVisVtxBins_
private

Definition at line 264 of file SplitVertexResolution.cc.

Referenced by beginJob().

const int SplitVertexResolution::nVtxBins_ = 60
staticprivate

Definition at line 263 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and SplitVertexResolution().

edm::Service<TFileService> SplitVertexResolution::outfile_
private

Definition at line 142 of file SplitVertexResolution.cc.

Referenced by beginJob(), and endJob().

TH1F* SplitVertexResolution::p_pullX_vsNtracks
private

Definition at line 238 of file SplitVertexResolution.cc.

Referenced by beginJob(), and endJob().

TH1F* SplitVertexResolution::p_pullX_vsNvtx
private

Definition at line 242 of file SplitVertexResolution.cc.

Referenced by beginJob(), and endJob().

TH1F* SplitVertexResolution::p_pullX_vsSumPt
private

Definition at line 234 of file SplitVertexResolution.cc.

Referenced by beginJob(), and endJob().

TH1F* SplitVertexResolution::p_pullY_vsNtracks
private

Definition at line 239 of file SplitVertexResolution.cc.

Referenced by beginJob(), and endJob().

TH1F* SplitVertexResolution::p_pullY_vsNvtx
private

Definition at line 243 of file SplitVertexResolution.cc.

Referenced by beginJob(), and endJob().

TH1F* SplitVertexResolution::p_pullY_vsSumPt
private

Definition at line 235 of file SplitVertexResolution.cc.

Referenced by beginJob(), and endJob().

TH1F* SplitVertexResolution::p_pullZ_vsNtracks
private

Definition at line 240 of file SplitVertexResolution.cc.

Referenced by beginJob(), and endJob().

TH1F* SplitVertexResolution::p_pullZ_vsNvtx
private

Definition at line 244 of file SplitVertexResolution.cc.

Referenced by beginJob(), and endJob().

TH1F* SplitVertexResolution::p_pullZ_vsSumPt
private

Definition at line 236 of file SplitVertexResolution.cc.

Referenced by beginJob(), and endJob().

TH1F* SplitVertexResolution::p_resolX_vsNtracks
private

Definition at line 213 of file SplitVertexResolution.cc.

Referenced by beginJob(), and endJob().

TH1F* SplitVertexResolution::p_resolX_vsNvtx
private

Definition at line 217 of file SplitVertexResolution.cc.

Referenced by beginJob(), and endJob().

TH1F* SplitVertexResolution::p_resolX_vsSumPt
private

Definition at line 209 of file SplitVertexResolution.cc.

Referenced by beginJob(), and endJob().

TH1F* SplitVertexResolution::p_resolY_vsNtracks
private

Definition at line 214 of file SplitVertexResolution.cc.

Referenced by beginJob(), and endJob().

TH1F* SplitVertexResolution::p_resolY_vsNvtx
private

Definition at line 218 of file SplitVertexResolution.cc.

Referenced by beginJob(), and endJob().

TH1F* SplitVertexResolution::p_resolY_vsSumPt
private

Definition at line 210 of file SplitVertexResolution.cc.

Referenced by beginJob(), and endJob().

TH1F* SplitVertexResolution::p_resolZ_vsNtracks
private

Definition at line 215 of file SplitVertexResolution.cc.

Referenced by beginJob(), and endJob().

TH1F* SplitVertexResolution::p_resolZ_vsNvtx
private

Definition at line 219 of file SplitVertexResolution.cc.

Referenced by beginJob(), and endJob().

TH1F* SplitVertexResolution::p_resolZ_vsSumPt
private

Definition at line 211 of file SplitVertexResolution.cc.

Referenced by beginJob(), and endJob().

edm::InputTag SplitVertexResolution::pvsTag_
private

Definition at line 119 of file SplitVertexResolution.cc.

edm::EDGetTokenT<reco::VertexCollection> SplitVertexResolution::pvsToken_
private

Definition at line 120 of file SplitVertexResolution.cc.

Referenced by analyze().

bool SplitVertexResolution::runControl_
private

Definition at line 137 of file SplitVertexResolution.cc.

Referenced by analyze().

std::vector<unsigned int> SplitVertexResolution::runControlNumbers_
private

Definition at line 138 of file SplitVertexResolution.cc.

Referenced by analyze(), beginJob(), and SplitVertexResolution().

edm::ESGetToken<RunInfo, RunInfoRcd> SplitVertexResolution::runInfoToken_
private

Definition at line 130 of file SplitVertexResolution.cc.

Referenced by getRunTime().

std::map<unsigned int, std::pair<long long, long long> > SplitVertexResolution::runNumbersTimesLog_
private

Definition at line 147 of file SplitVertexResolution.cc.

Referenced by beginRun(), and endJob().

bool SplitVertexResolution::storeNtuple_
private

Definition at line 113 of file SplitVertexResolution.cc.

Referenced by analyze().

const float SplitVertexResolution::sumpTEndScale_
private

Definition at line 255 of file SplitVertexResolution.cc.

Referenced by beginJob(), and SplitVertexResolution().

const float SplitVertexResolution::sumpTStartScale_
private

Definition at line 254 of file SplitVertexResolution.cc.

Referenced by beginJob(), and SplitVertexResolution().

TH1D* SplitVertexResolution::tksByTrigger_
private

Definition at line 192 of file SplitVertexResolution.cc.

Referenced by endJob().

edm::InputTag SplitVertexResolution::tracksTag_
private

Definition at line 122 of file SplitVertexResolution.cc.

edm::EDGetTokenT<reco::TrackCollection> SplitVertexResolution::tracksToken_
private

Definition at line 123 of file SplitVertexResolution.cc.

Referenced by analyze().

edm::ESGetToken<TransientTrackBuilder, TransientTrackRecord> SplitVertexResolution::transientTrackBuilderToken_
private

Definition at line 129 of file SplitVertexResolution.cc.

Referenced by analyze().

TTree* SplitVertexResolution::tree_
private

Definition at line 251 of file SplitVertexResolution.cc.

Referenced by analyze(), and beginJob().

std::map<std::string, std::pair<int, int> > SplitVertexResolution::triggerMap_
private

Definition at line 267 of file SplitVertexResolution.cc.

Referenced by analyze(), and endJob().

edm::InputTag SplitVertexResolution::triggerResultsTag_ = edm::InputTag("TriggerResults", "", "HLT")
private

Definition at line 125 of file SplitVertexResolution.cc.

edm::EDGetTokenT<edm::TriggerResults> SplitVertexResolution::triggerResultsToken_
private

Definition at line 126 of file SplitVertexResolution.cc.

Referenced by analyze().

Generated for CMSSW Reference Manual by doxygen 1.8.5