#include <dEdxHitAnalyzer.h>

Inheritance diagram for dEdxHitAnalyzer:

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

	dEdxHitAnalyzer (const edm::ParameterSet &)

void	dqmBeginRun (const edm::Run &, const edm::EventSetup &) override

double	harmonic2 (const reco::DeDxHitInfo *dedxHits)

	~dEdxHitAnalyzer () 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

	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_

int	dEdxHarm2Bin

double	dEdxHarm2Max

double	dEdxHarm2Min

std::vector< std::string >	dEdxInputList_

std::vector< dEdxMEs >	dEdxMEsVector

int	dEdxNHitBin

double	dEdxNHitMax

double	dEdxNHitMin

int	dEdxPixelBin

double	dEdxPixelMax

double	dEdxPixelMin

int	dEdxStripBin

double	dEdxStripMax

double	dEdxStripMin

std::vector< edm::EDGetTokenT< reco::DeDxHitInfoAss > >	dEdxTokenList_

bool	doAllPlots_

bool	doDeDxPlots_

edm::ParameterSet	fullconf_

GenericTriggerEventFlag *	genTriggerEventFlag_

std::string	histname

std::string	MEFolderName

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 dEdxHitAnalyzer.h.

Constructor & Destructor Documentation

◆ dEdxHitAnalyzer()

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

explicit

Definition at line 22 of file dEdxHitAnalyzer.cc.

     : 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> >("deDxHitProducers");
   for (auto const& tag : dEdxInputList_) {
     dEdxTokenList_.push_back(consumes<reco::DeDxHitInfoAss>(edm::InputTag(tag)));
   }
  
   // parameters from the configuration
   MEFolderName = conf_.getParameter<std::string>("FolderName");
  
   dEdxNHitBin = conf_.getParameter<int>("dEdxNHitBin");
   dEdxNHitMin = conf_.getParameter<double>("dEdxNHitMin");
   dEdxNHitMax = conf_.getParameter<double>("dEdxNHitMax");
  
   dEdxStripBin = conf_.getParameter<int>("dEdxStripBin");
   dEdxStripMin = conf_.getParameter<double>("dEdxStripMin");
   dEdxStripMax = conf_.getParameter<double>("dEdxStripMax");
  
   dEdxPixelBin = conf_.getParameter<int>("dEdxPixelBin");
   dEdxPixelMin = conf_.getParameter<double>("dEdxPixelMin");
   dEdxPixelMax = conf_.getParameter<double>("dEdxPixelMax");
  
   dEdxHarm2Bin = conf_.getParameter<int>("dEdxHarm2Bin");
   dEdxHarm2Min = conf_.getParameter<double>("dEdxHarm2Min");
   dEdxHarm2Max = conf_.getParameter<double>("dEdxHarm2Max");
 }

References conf_, dEdxHarm2Bin, dEdxHarm2Max, dEdxHarm2Min, dEdxInputList_, dEdxNHitBin, dEdxNHitMax, dEdxNHitMin, dEdxPixelBin, dEdxPixelMax, dEdxPixelMin, dEdxStripBin, dEdxStripMax, dEdxStripMin, dEdxTokenList_, edm::ParameterSet::getParameter(), HLT_2018_cff::InputTag, MEFolderName, AlCaHLTBitMon_QueryRunRegistry::string, GlobalPosition_Frontier_DevDB_cff::tag, trackInputTag_, and trackToken_.

◆ ~dEdxHitAnalyzer()

dEdxHitAnalyzer::~dEdxHitAnalyzer ( )

override

Definition at line 57 of file dEdxHitAnalyzer.cc.

                                   {
   if (genTriggerEventFlag_)
     delete genTriggerEventFlag_;
 }

References genTriggerEventFlag_.

Member Function Documentation

◆ analyze()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 128 of file dEdxHitAnalyzer.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::DeDxHitInfoAss> dEdxObjectHandle;
       iEvent.getByToken(dEdxTokenList_[i], dEdxObjectHandle);
       if (!dEdxObjectHandle.isValid())
         continue;
  
       for (unsigned int t = 0; t < trackCollectionHandle->size(); t++) {
         reco::TrackRef track = reco::TrackRef(trackCollectionHandle, t);
  
         if (track->quality(reco::TrackBase::highPurity)) {
           const reco::DeDxHitInfo* dedxHits = nullptr;
           if (!track.isNull()) {
             reco::DeDxHitInfoRef dedxHitsRef = (*dEdxObjectHandle)[track];
             if (!dedxHitsRef.isNull())
               dedxHits = &(*dedxHitsRef);
           }
           if (!dedxHits)
             continue;
  
           for (unsigned int h = 0; h < dedxHits->size(); h++) {
             DetId detid(dedxHits->detId(h));
             if (detid.subdetId() >= 3)
               dEdxMEsVector[i].ME_StripHitDeDx->Fill(dedxHits->charge(h));
             if (detid.subdetId() < 3)
               dEdxMEsVector[i].ME_PixelHitDeDx->Fill(dedxHits->charge(h));
           }
           dEdxMEsVector[i].ME_NHitDeDx->Fill(dedxHits->size());
           dEdxMEsVector[i].ME_Harm2DeDx->Fill(harmonic2(dedxHits));
         }
       }
     }
   }
 }

References GenericTriggerEventFlag::accept(), reco::DeDxHitInfo::charge(), dEdxInputList_, dEdxMEsVector, dEdxTokenList_, reco::DeDxHitInfo::detId(), doAllPlots_, doDeDxPlots_, genTriggerEventFlag_, harmonic2(), reco::TrackBase::highPurity, mps_fire::i, iEvent, edm::Ref< C, T, F >::isNull(), edm::HandleBase::isValid(), GenericTriggerEventFlag::on(), reco::DeDxHitInfo::size(), OrderedSet::t, HLT_2018_cff::track, and trackToken_.

◆ bookHistograms()

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 70 of file dEdxHitAnalyzer.cc.

                                                                                                               {
   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 = "Strip_dEdxPerCluster_";
       dEdxMEsVector[i].ME_StripHitDeDx = ibooker.book1D(histname, histname, dEdxStripBin, dEdxStripMin, dEdxStripMax);
       dEdxMEsVector[i].ME_StripHitDeDx->setAxisTitle("dEdx of on-track strip cluster (ADC)");
       dEdxMEsVector[i].ME_StripHitDeDx->setAxisTitle("Number of Strip clusters", 2);
  
       histname = "Pixel_dEdxPerCluster_";
       dEdxMEsVector[i].ME_PixelHitDeDx = ibooker.book1D(histname, histname, dEdxPixelBin, dEdxPixelMin, dEdxPixelMax);
       dEdxMEsVector[i].ME_PixelHitDeDx->setAxisTitle("dEdx of on-track pixel cluster (ADC)");
       dEdxMEsVector[i].ME_PixelHitDeDx->setAxisTitle("Number of Pixel clusters", 2);
  
       histname = "NumberOfdEdxHitsPerTrack_";
       dEdxMEsVector[i].ME_NHitDeDx = ibooker.book1D(histname, histname, dEdxNHitBin, dEdxNHitMin, dEdxNHitMax);
       dEdxMEsVector[i].ME_NHitDeDx->setAxisTitle("Number of dEdxHits per Track");
       dEdxMEsVector[i].ME_NHitDeDx->setAxisTitle("Number of Tracks", 2);
  
       histname = "Harm2_dEdxPerTrack_";
       dEdxMEsVector[i].ME_Harm2DeDx = ibooker.book1D(histname, histname, dEdxHarm2Bin, dEdxHarm2Min, dEdxHarm2Max);
       dEdxMEsVector[i].ME_Harm2DeDx->setAxisTitle("Harmonic2 dEdx estimator for each Track");
       dEdxMEsVector[i].ME_Harm2DeDx->setAxisTitle("Number of Tracks", 2);
     }
   }
 }

References dqm::implementation::IBooker::book1D(), dEdxHarm2Bin, dEdxHarm2Max, dEdxHarm2Min, dEdxInputList_, dEdxMEsVector, dEdxNHitBin, dEdxNHitMax, dEdxNHitMin, dEdxPixelBin, dEdxPixelMax, dEdxPixelMin, dEdxStripBin, dEdxStripMax, dEdxStripMin, doAllPlots_, doDeDxPlots_, histname, mps_fire::i, MEFolderName, and dqm::implementation::NavigatorBase::setCurrentFolder().

◆ dqmBeginRun()

void dEdxHitAnalyzer::dqmBeginRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	iSetup
	)

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 64 of file dEdxHitAnalyzer.cc.

                                                                                  {
   // Initialize the GenericTriggerEventFlag
   if (genTriggerEventFlag_->on())
     genTriggerEventFlag_->initRun(iRun, iSetup);
 }

References genTriggerEventFlag_, GenericTriggerEventFlag::initRun(), and GenericTriggerEventFlag::on().

◆ fillDescriptions()

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

static

Definition at line 174 of file dEdxHitAnalyzer.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().

◆ harmonic2()

double dEdxHitAnalyzer::harmonic2 ( const reco::DeDxHitInfo * dedxHits )

Definition at line 104 of file dEdxHitAnalyzer.cc.

                                                                  {
   if (!dedxHits)
     return -1;
   std::vector<double> vect;
   for (unsigned int h = 0; h < dedxHits->size(); h++) {
     DetId detid(dedxHits->detId(h));
     double Norm = (detid.subdetId() < 3) ? 3.61e-06 : 3.61e-06 * 265;
     double ChargeOverPathlength = Norm * dedxHits->charge(h) / dedxHits->pathlength(h);
     vect.push_back(ChargeOverPathlength);  //save charge
   }
  
   int size = vect.size();
   if (size <= 0)
     return -1;
   double result = 0;
   double expo = -2;
   for (int i = 0; i < size; i++) {
     result += pow(vect[i], expo);
   }
   return pow(result / size, 1. / expo);
 }