#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	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

	dEdxAnalyzer (const edm::ParameterSet &)

void	endJob () override

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

double	mass (double P, double I)

	~dEdxAnalyzer () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

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

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) override

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

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

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::vector< ModuleDescription const * > &modules, 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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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 edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

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)

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

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Monitoring source for general quantities related to track dEdx.

Definition at line 35 of file dEdxAnalyzer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 22 of file dEdxAnalyzer.cc.

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

   : 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) ) );
   }
 }

dEdxAnalyzer::~dEdxAnalyzer ( )

override

Definition at line 40 of file dEdxAnalyzer.cc.

References genTriggerEventFlag_.

 { 
 
   if (genTriggerEventFlag_)      delete genTriggerEventFlag_;
 }

Member Function Documentation

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

override

Definition at line 164 of file dEdxAnalyzer.cc.

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

 {
 
   // 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()));
               }
             }
          }
       }
    }
 }

void dEdxAnalyzer::beginLuminosityBlock	(	edm::LuminosityBlock const &	,
		edm::EventSetup const &
	)

override

Definition at line 215 of file dEdxAnalyzer.cc.

216 {

217 }

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 70 of file dEdxAnalyzer.cc.

References DQMStore::IBooker::book1D(), conf_, dEdxC, dEdxInputList_, dEdxK, dEdxMEsVector, doAllPlots_, doDeDxPlots_, genTriggerEventFlag_, edm::ParameterSet::getParameter(), HighPtThreshold, HIPdEdxMin, histname, mps_fire::i, GenericTriggerEventFlag::initRun(), GenericTriggerEventFlag::on(), DQMStore::IBooker::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, and TrackHitMin.

 {
 
   // 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);
        }
     }
 
 }

void dEdxAnalyzer::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDProducerBase.

Definition at line 48 of file dEdxAnalyzer.cc.

References conf_, dqmStore_, edm::ParameterSet::getParameter(), postProcessorL1Gen_cff::outputFileName, DQMStore::save(), DQMStore::showDirStructure(), and AlCaHLTBitMon_QueryRunRegistry::string.

 {
     bool outputMEsInRootFile   = conf_.getParameter<bool>("OutputMEsInRootFile");
     std::string outputFileName = conf_.getParameter<std::string>("OutputFileName");
     if(outputMEsInRootFile)
     {
       dqmStore_->showDirStructure();
       dqmStore_->save(outputFileName);
     }
 }

void dEdxAnalyzer::endLuminosityBlock	(	edm::LuminosityBlock const &	,
		edm::EventSetup const &
	)

override

Definition at line 221 of file dEdxAnalyzer.cc.

222 {

223 }

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

static

Definition at line 227 of file dEdxAnalyzer.cc.

References edm::ConfigurationDescriptions::addDefault(), DEFINE_FWK_MODULE, and edm::ParameterSetDescription::setUnknown().

                                                                          {
   //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);
 }

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

Definition at line 157 of file dEdxAnalyzer.cc.

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

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

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

Member Data Documentation

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

private

Definition at line 92 of file dEdxAnalyzer.h.

edm::ParameterSet dEdxAnalyzer::conf_

private

Definition at line 56 of file dEdxAnalyzer.h.

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

double dEdxAnalyzer::dEdxC

private

Definition at line 83 of file dEdxAnalyzer.h.

Referenced by bookHistograms(), and mass().

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

private

Definition at line 88 of file dEdxAnalyzer.h.

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

double dEdxAnalyzer::dEdxK

private

Definition at line 83 of file dEdxAnalyzer.h.

Referenced by bookHistograms(), and mass().

std::vector< dEdxMEs > dEdxAnalyzer::dEdxMEsVector

private

Definition at line 93 of file dEdxAnalyzer.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 89 of file dEdxAnalyzer.h.

Referenced by analyze(), and dEdxAnalyzer().

bool dEdxAnalyzer::doAllPlots_

private

Definition at line 58 of file dEdxAnalyzer.h.

Referenced by analyze(), and bookHistograms().

bool dEdxAnalyzer::doDeDxPlots_

private

Definition at line 59 of file dEdxAnalyzer.h.

Referenced by analyze(), and bookHistograms().

DQMStore* dEdxAnalyzer::dqmStore_

private

Definition at line 54 of file dEdxAnalyzer.h.

Referenced by endJob().

edm::ParameterSet dEdxAnalyzer::fullconf_

private

Definition at line 55 of file dEdxAnalyzer.h.

GenericTriggerEventFlag* dEdxAnalyzer::genTriggerEventFlag_

private

Definition at line 96 of file dEdxAnalyzer.h.

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

double dEdxAnalyzer::HighPtThreshold

private

Definition at line 82 of file dEdxAnalyzer.h.

Referenced by analyze(), and bookHistograms().

double dEdxAnalyzer::HIPdEdxMin

private

Definition at line 82 of file dEdxAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::string dEdxAnalyzer::histname

private

Definition at line 94 of file dEdxAnalyzer.h.

Referenced by bookHistograms().

double dEdxAnalyzer::TrackHitMin

private

Definition at line 82 of file dEdxAnalyzer.h.

Referenced by analyze(), and bookHistograms().

edm::InputTag dEdxAnalyzer::trackInputTag_

private

Definition at line 85 of file dEdxAnalyzer.h.

Referenced by dEdxAnalyzer().

std::string dEdxAnalyzer::TrackName

private

Definition at line 91 of file dEdxAnalyzer.h.

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

private

Definition at line 86 of file dEdxAnalyzer.h.

Referenced by analyze(), and dEdxAnalyzer().

Classes

Public Member Functions

Static Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation