#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 &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 () noexcept(false) 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_

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

Constructor & Destructor Documentation

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

explicit

Definition at line 22 of file dEdxHitAnalyzer.cc.

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

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

dEdxHitAnalyzer::~dEdxHitAnalyzer ( )

override

Definition at line 58 of file dEdxHitAnalyzer.cc.

References genTriggerEventFlag_.

 { 
 
   if (genTriggerEventFlag_)      delete genTriggerEventFlag_;
 }

Member Function Documentation

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

override

Definition at line 134 of file dEdxHitAnalyzer.cc.

References GenericTriggerEventFlag::accept(), reco::DeDxHitInfo::charge(), dEdxInputList_, dEdxMEsVector, dEdxTokenList_, reco::DeDxHitInfo::detId(), doAllPlots_, doDeDxPlots_, HcalObjRepresent::Fill(), genTriggerEventFlag_, edm::Event::getByToken(), h, harmonic2(), reco::TrackBase::highPurity, mps_fire::i, edm::Ref< C, T, F >::isNull(), edm::HandleBase::isValid(), GenericTriggerEventFlag::on(), reco::DeDxHitInfo::size(), lumiQTWidget::t, HiIsolationCommonParameters_cff::track, 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::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));
             }
          }
       }
    }
 }

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 73 of file dEdxHitAnalyzer.cc.

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

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

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

override

Definition at line 66 of file dEdxHitAnalyzer.cc.

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

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

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

static

Definition at line 180 of file dEdxHitAnalyzer.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 dEdxHitAnalyzer::harmonic2 ( const reco::DeDxHitInfo * dedxHits )

Definition at line 111 of file dEdxHitAnalyzer.cc.

References reco::DeDxHitInfo::charge(), reco::DeDxHitInfo::detId(), MillePedeFileConverter_cfg::e, h, mps_fire::i, reco::DeDxHitInfo::pathlength(), funct::pow(), mps_fire::result, reco::DeDxHitInfo::size(), and findQualityFiles::size.

Referenced by analyze().

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