dEdxAnalyzer Class Reference

#include <dEdxAnalyzer.h>

Inheritance diagram for dEdxAnalyzer:

Classes
struct	dEdxMEs

Public Member Functions
void	analyze (const edm::Event &iEvent, const edm::EventSetup &iSetup) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	dEdxAnalyzer (const edm::ParameterSet &)
double	mass (double P, double I)
	~dEdxAnalyzer () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)
static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Private Attributes
std::vector< std::string >	AlgoNames
edm::ParameterSet	conf_
double	dEdxC
std::vector< std::string >	dEdxInputList_
double	dEdxK
std::vector< dEdxMEs >	dEdxMEsVector
std::vector< edm::EDGetTokenT< reco::DeDxDataValueMap > >	dEdxTokenList_
bool	doAllPlots_
bool	doDeDxPlots_
DQMStore *	dqmStore_
edm::ParameterSet	fullconf_
GenericTriggerEventFlag *	genTriggerEventFlag_
double	HighPtThreshold
double	HIPdEdxMin
std::string	histname
double	TrackHitMin
edm::InputTag	trackInputTag_
std::string	TrackName
edm::EDGetTokenT< reco::TrackCollection >	trackToken_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Monitoring source for general quantities related to track dEdx.

Definition at line 34 of file dEdxAnalyzer.h.

Constructor & Destructor Documentation

dEdxAnalyzer()

dEdxAnalyzer::dEdxAnalyzer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 22 of file dEdxAnalyzer.cc.

    : dqmStore_(edm::Service<DQMStore>().operator->()),
      fullconf_(iConfig),
      conf_(fullconf_.getParameter<edm::ParameterSet>("dEdxParameters")),
      doAllPlots_(conf_.getParameter<bool>("doAllPlots")),
      doDeDxPlots_(conf_.getParameter<bool>("doDeDxPlots")),
      genTriggerEventFlag_(new GenericTriggerEventFlag(
          conf_.getParameter<edm::ParameterSet>("genericTriggerEventPSet"), consumesCollector(), *this)) {
  trackInputTag_ = edm::InputTag(conf_.getParameter<std::string>("TracksForDeDx"));
  trackToken_ = consumes<reco::TrackCollection>(trackInputTag_);
 
  dEdxInputList_ = conf_.getParameter<std::vector<std::string> >("deDxProducers");
  for (auto const& tag : dEdxInputList_) {
    dEdxTokenList_.push_back(consumes<reco::DeDxDataValueMap>(edm::InputTag(tag)));
  }
}

References conf_, dEdxInputList_, dEdxTokenList_, edm::ParameterSet::getParameter(), HLT_FULL_cff::InputTag, AlCaHLTBitMon_QueryRunRegistry::string, GlobalPosition_Frontier_DevDB_cff::tag, trackInputTag_, and trackToken_.

~dEdxAnalyzer()

dEdxAnalyzer::~dEdxAnalyzer ( )

override

Definition at line 39 of file dEdxAnalyzer.cc.

                            {
  if (genTriggerEventFlag_)
    delete genTriggerEventFlag_;
}

References genTriggerEventFlag_.

Member Function Documentation

analyze()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 146 of file dEdxAnalyzer.cc.

                                                                              {
  // Filter out events if Trigger Filtering is requested
  if (genTriggerEventFlag_->on() && !genTriggerEventFlag_->accept(iEvent, iSetup))
    return;
 
  if (doDeDxPlots_ || doAllPlots_) {
    edm::Handle<reco::TrackCollection> trackCollectionHandle;
    iEvent.getByToken(trackToken_, trackCollectionHandle);
    if (!trackCollectionHandle.isValid())
      return;
 
    for (unsigned int i = 0; i < dEdxInputList_.size(); i++) {
      edm::Handle<reco::DeDxDataValueMap> dEdxObjectHandle;
      iEvent.getByToken(dEdxTokenList_[i], dEdxObjectHandle);
      if (!dEdxObjectHandle.isValid())
        continue;
      const edm::ValueMap<reco::DeDxData> dEdxColl = *dEdxObjectHandle.product();
 
      for (unsigned int t = 0; t < trackCollectionHandle->size(); t++) {
        reco::TrackRef track = reco::TrackRef(trackCollectionHandle, t);
 
        if (track->quality(reco::TrackBase::highPurity)) {
          //MIPs
          if (track->pt() >= 5.0 && track->numberOfValidHits() > TrackHitMin) {
            dEdxMEsVector[i].ME_MipDeDx->Fill(dEdxColl[track].dEdx());
            dEdxMEsVector[i].ME_MipDeDxNHits->Fill(dEdxColl[track].numberOfMeasurements());
            if (dEdxColl[track].numberOfMeasurements() != 0)
              dEdxMEsVector[i].ME_MipDeDxNSatHits->Fill((1.0 * dEdxColl[track].numberOfSaturatedMeasurements()) /
                                                        dEdxColl[track].numberOfMeasurements());
            dEdxMEsVector[i].ME_MipDeDxMass->Fill(mass(track->p(), dEdxColl[track].dEdx()));
 
            if (track->pt() >= HighPtThreshold) {
              dEdxMEsVector[i].ME_MipHighPtDeDx->Fill(dEdxColl[track].dEdx());
              dEdxMEsVector[i].ME_MipHighPtDeDxNHits->Fill(dEdxColl[track].numberOfMeasurements());
            }
 
            //HighlyIonizing particles
          } else if (track->pt() < 2 && dEdxColl[track].dEdx() > HIPdEdxMin) {
            dEdxMEsVector[i].ME_HipDeDxMass->Fill(mass(track->p(), dEdxColl[track].dEdx()));
          }
        }
      }
    }
  }
}

References GenericTriggerEventFlag::accept(), plot_hgcal_utils::dEdx, dEdxInputList_, dEdxMEsVector, dEdxTokenList_, doAllPlots_, doDeDxPlots_, genTriggerEventFlag_, HighPtThreshold, reco::TrackBase::highPurity, HIPdEdxMin, mps_fire::i, iEvent, edm::HandleBase::isValid(), mass(), GenericTriggerEventFlag::on(), edm::Handle< T >::product(), submitPVValidationJobs::t, HLT_FULL_cff::track, TrackHitMin, and trackToken_.

bookHistograms()

void dEdxAnalyzer::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  iRun, edm::EventSetup const &  iSetup  )

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 55 of file dEdxAnalyzer.cc.

                                                                                                           {
  // Initialize the GenericTriggerEventFlag
  if (genTriggerEventFlag_->on())
    genTriggerEventFlag_->initRun(iRun, iSetup);
 
  // parameters from the configuration
  std::string MEFolderName = conf_.getParameter<std::string>("FolderName");
 
  // get binning from the configuration
  TrackHitMin = conf_.getParameter<double>("TrackHitMin");
  HIPdEdxMin = conf_.getParameter<double>("HIPdEdxMin");
  HighPtThreshold = conf_.getParameter<double>("HighPtThreshold");
 
  dEdxK = conf_.getParameter<double>("dEdxK");
  dEdxC = conf_.getParameter<double>("dEdxC");
 
  int dEdxNHitBin = conf_.getParameter<int>("dEdxNHitBin");
  double dEdxNHitMin = conf_.getParameter<double>("dEdxNHitMin");
  double dEdxNHitMax = conf_.getParameter<double>("dEdxNHitMax");
 
  int dEdxBin = conf_.getParameter<int>("dEdxBin");
  double dEdxMin = conf_.getParameter<double>("dEdxMin");
  double dEdxMax = conf_.getParameter<double>("dEdxMax");
 
  int dEdxHIPmassBin = conf_.getParameter<int>("dEdxHIPmassBin");
  double dEdxHIPmassMin = conf_.getParameter<double>("dEdxHIPmassMin");
  double dEdxHIPmassMax = conf_.getParameter<double>("dEdxHIPmassMax");
 
  int dEdxMIPmassBin = conf_.getParameter<int>("dEdxMIPmassBin");
  double dEdxMIPmassMin = conf_.getParameter<double>("dEdxMIPmassMin");
  double dEdxMIPmassMax = conf_.getParameter<double>("dEdxMIPmassMax");
 
  ibooker.setCurrentFolder(MEFolderName);
 
  // book the Hit Property histograms
  // ---------------------------------------------------------------------------------//
 
  if (doDeDxPlots_ || doAllPlots_) {
    for (unsigned int i = 0; i < dEdxInputList_.size(); i++) {
      ibooker.setCurrentFolder(MEFolderName + "/" + dEdxInputList_[i]);
      dEdxMEsVector.push_back(dEdxMEs());
 
      histname = "MIP_dEdxPerTrack_";
      dEdxMEsVector[i].ME_MipDeDx = ibooker.book1D(histname, histname, dEdxBin, dEdxMin, dEdxMax);
      dEdxMEsVector[i].ME_MipDeDx->setAxisTitle("dEdx of each MIP Track (MeV/cm)");
      dEdxMEsVector[i].ME_MipDeDx->setAxisTitle("Number of Tracks", 2);
 
      histname = "MIP_NumberOfdEdxHitsPerTrack_";
      dEdxMEsVector[i].ME_MipDeDxNHits = ibooker.book1D(histname, histname, dEdxNHitBin, dEdxNHitMin, dEdxNHitMax);
      dEdxMEsVector[i].ME_MipDeDxNHits->setAxisTitle("Number of dEdxHits of each MIP Track");
      dEdxMEsVector[i].ME_MipDeDxNHits->setAxisTitle("Number of Tracks", 2);
 
      histname = "MIP_FractionOfSaturateddEdxHitsPerTrack_";
      dEdxMEsVector[i].ME_MipDeDxNSatHits = ibooker.book1D(histname, histname, 2 * dEdxNHitBin, 0, 1);
      dEdxMEsVector[i].ME_MipDeDxNSatHits->setAxisTitle("Fraction of Saturated dEdxHits of each MIP Track");
      dEdxMEsVector[i].ME_MipDeDxNSatHits->setAxisTitle("Number of Tracks", 2);
 
      histname = "MIP_MassPerTrack_";
      dEdxMEsVector[i].ME_MipDeDxMass =
          ibooker.book1D(histname, histname, dEdxMIPmassBin, dEdxMIPmassMin, dEdxMIPmassMax);
      dEdxMEsVector[i].ME_MipDeDxMass->setAxisTitle("dEdx Mass of each MIP Track (GeV/c^{2})");
      dEdxMEsVector[i].ME_MipDeDxMass->setAxisTitle("Number of Tracks", 2);
 
      histname = "HIP_MassPerTrack_";
      dEdxMEsVector[i].ME_HipDeDxMass =
          ibooker.book1D(histname, histname, dEdxHIPmassBin, dEdxHIPmassMin, dEdxHIPmassMax);
      dEdxMEsVector[i].ME_HipDeDxMass->setAxisTitle("dEdx Mass of each HIP Track (GeV/c^{2})");
      dEdxMEsVector[i].ME_HipDeDxMass->setAxisTitle("Number of Tracks", 2);
 
      histname = "MIPOfHighPt_dEdxPerTrack_";
      dEdxMEsVector[i].ME_MipHighPtDeDx = ibooker.book1D(histname, histname, dEdxBin, dEdxMin, dEdxMax);
      dEdxMEsVector[i].ME_MipHighPtDeDx->setAxisTitle("dEdx of each MIP (of High pT) Track (MeV/cm)");
      dEdxMEsVector[i].ME_MipHighPtDeDx->setAxisTitle("Number of Tracks", 2);
 
      histname = "MIPOfHighPt_NumberOfdEdxHitsPerTrack_";
      dEdxMEsVector[i].ME_MipHighPtDeDxNHits =
          ibooker.book1D(histname, histname, dEdxNHitBin, dEdxNHitMin, dEdxNHitMax);
      dEdxMEsVector[i].ME_MipHighPtDeDxNHits->setAxisTitle("Number of dEdxHits of each MIP (of High pT) Track");
      dEdxMEsVector[i].ME_MipHighPtDeDxNHits->setAxisTitle("Number of Tracks", 2);
    }
  }
}

References dqm::implementation::IBooker::book1D(), conf_, dEdxAnalyzer_cfi::dEdxBin, dEdxC, dEdxAnalyzer_cfi::dEdxHIPmassBin, dEdxAnalyzer_cfi::dEdxHIPmassMax, dEdxAnalyzer_cfi::dEdxHIPmassMin, dEdxInputList_, dEdxK, dEdxAnalyzer_cfi::dEdxMax, dEdxMEsVector, dEdxAnalyzer_cfi::dEdxMin, dEdxAnalyzer_cfi::dEdxMIPmassBin, dEdxAnalyzer_cfi::dEdxMIPmassMax, dEdxAnalyzer_cfi::dEdxMIPmassMin, dEdxAnalyzer_cfi::dEdxNHitBin, dEdxAnalyzer_cfi::dEdxNHitMax, dEdxAnalyzer_cfi::dEdxNHitMin, doAllPlots_, doDeDxPlots_, genTriggerEventFlag_, edm::ParameterSet::getParameter(), HighPtThreshold, HIPdEdxMin, histname, mps_fire::i, GenericTriggerEventFlag::initRun(), GenericTriggerEventFlag::on(), dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, and TrackHitMin.

fillDescriptions()

void dEdxAnalyzer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 193 of file dEdxAnalyzer.cc.

                                                                              {
  //The following says we do not know what parameters are allowed so do no validation
  // Please change this to state exactly what you do use, even if it is no parameters
  edm::ParameterSetDescription desc;
  desc.setUnknown();
  descriptions.addDefault(desc);
}

References edm::ConfigurationDescriptions::addDefault(), and submitPVResolutionJobs::desc.

mass()

double dEdxAnalyzer::mass	(	double	P,
		double	I
	)

Definition at line 138 of file dEdxAnalyzer.cc.

                                            {
  if (I - dEdxC < 0)
    return -1;
  return sqrt((I - dEdxC) / dEdxK) * P;
}

References dEdxC, dEdxK, Exhume::I, and mathSSE::sqrt().

Referenced by Particle.Particle::__str__(), DiObject.DiMuon::__str__(), and analyze().

Member Data Documentation

AlgoNames

std::vector<std::string> dEdxAnalyzer::AlgoNames

private

Definition at line 84 of file dEdxAnalyzer.h.

conf_

edm::ParameterSet dEdxAnalyzer::conf_

private

Definition at line 50 of file dEdxAnalyzer.h.

Referenced by bookHistograms(), and dEdxAnalyzer().

dEdxC

double dEdxAnalyzer::dEdxC

private

Definition at line 75 of file dEdxAnalyzer.h.

Referenced by bookHistograms(), and mass().

dEdxInputList_

std::vector<std::string> dEdxAnalyzer::dEdxInputList_

private

Definition at line 80 of file dEdxAnalyzer.h.

Referenced by analyze(), bookHistograms(), and dEdxAnalyzer().

dEdxK

double dEdxAnalyzer::dEdxK

private

Definition at line 75 of file dEdxAnalyzer.h.

Referenced by bookHistograms(), and mass().

dEdxMEsVector

std::vector<dEdxMEs> dEdxAnalyzer::dEdxMEsVector

private

Definition at line 85 of file dEdxAnalyzer.h.

Referenced by analyze(), and bookHistograms().

dEdxTokenList_

std::vector<edm::EDGetTokenT<reco::DeDxDataValueMap> > dEdxAnalyzer::dEdxTokenList_

private

Definition at line 81 of file dEdxAnalyzer.h.

Referenced by analyze(), and dEdxAnalyzer().

doAllPlots_

bool dEdxAnalyzer::doAllPlots_

private

Definition at line 52 of file dEdxAnalyzer.h.

Referenced by analyze(), and bookHistograms().

doDeDxPlots_

bool dEdxAnalyzer::doDeDxPlots_

private

Definition at line 53 of file dEdxAnalyzer.h.

Referenced by analyze(), and bookHistograms().

dqmStore_

DQMStore* dEdxAnalyzer::dqmStore_

private

Definition at line 48 of file dEdxAnalyzer.h.

fullconf_

edm::ParameterSet dEdxAnalyzer::fullconf_

private

Definition at line 49 of file dEdxAnalyzer.h.

genTriggerEventFlag_

GenericTriggerEventFlag* dEdxAnalyzer::genTriggerEventFlag_

private

Definition at line 88 of file dEdxAnalyzer.h.

Referenced by analyze(), bookHistograms(), and ~dEdxAnalyzer().

HighPtThreshold

double dEdxAnalyzer::HighPtThreshold

private

Definition at line 74 of file dEdxAnalyzer.h.

Referenced by analyze(), and bookHistograms().

HIPdEdxMin

double dEdxAnalyzer::HIPdEdxMin

private

Definition at line 74 of file dEdxAnalyzer.h.

Referenced by analyze(), and bookHistograms().

histname

std::string dEdxAnalyzer::histname

private

Definition at line 86 of file dEdxAnalyzer.h.

Referenced by bookHistograms().

TrackHitMin

double dEdxAnalyzer::TrackHitMin

private

Definition at line 74 of file dEdxAnalyzer.h.

Referenced by analyze(), and bookHistograms().

trackInputTag_

edm::InputTag dEdxAnalyzer::trackInputTag_

private

Definition at line 77 of file dEdxAnalyzer.h.

Referenced by dEdxAnalyzer().

TrackName

std::string dEdxAnalyzer::TrackName

private

Definition at line 83 of file dEdxAnalyzer.h.

trackToken_

edm::EDGetTokenT<reco::TrackCollection> dEdxAnalyzer::trackToken_

private

Definition at line 78 of file dEdxAnalyzer.h.

Referenced by analyze(), and dEdxAnalyzer().