#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 one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

Static Public Member Functions
static void	fillDescriptions (edm::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_

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(), protons_cff::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
	)

override

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