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	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () 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_2018_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(), OrderedSet::t, HLT_2018_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 edm::ParameterSetDescription::setUnknown().

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().