#include <DQM/TrackerMonitorTrack/src/TkAlCaRecoMonitor.cc>

Inheritance diagram for TkAlCaRecoMonitor:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

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

	TkAlCaRecoMonitor (const edm::ParameterSet &)

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

Private Member Functions
void	fillHitmaps (const reco::Track &track, const TrackerGeometry &geometry)

void	fillRawIdMap (const TrackerGeometry &geometry)

Private Attributes
MonitorElement *	AlCaRecoTrackEfficiency_

std::map< int, int >	binByRawId_

edm::ParameterSet	conf_

double	daughterMass_

bool	fillInvariantMass_

bool	fillRawIdMap_

MonitorElement *	Hits_perDetId_

MonitorElement *	Hits_XvsY_

MonitorElement *	Hits_ZvsR_

MonitorElement *	invariantMass_

edm::EDGetTokenT< reco::CaloJetCollection >	jetCollection_

MonitorElement *	jetPt_

double	maxJetPt_

MonitorElement *	minJetDeltaR_

MonitorElement *	minTrackDeltaR_

edm::EDGetTokenT< reco::TrackCollection >	referenceTrackProducer_

bool	runsOnReco_

MonitorElement *	sumCharge_

MonitorElement *	TrackCurvature_

edm::EDGetTokenT< reco::TrackCollection >	trackProducer_

MonitorElement *	TrackPtNegative_

MonitorElement *	TrackPtPositive_

MonitorElement *	TrackQuality_

bool	useSignedR_

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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 special quantities related to Tracker Alignment AlCaReco Production.

Definition at line 32 of file TkAlCaRecoMonitor.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 27 of file TkAlCaRecoMonitor.cc.

References conf_, edm::ParameterSet::getParameter(), jetCollection_, referenceTrackProducer_, and trackProducer_.

                                                                    {
   conf_ = iConfig;
   trackProducer_ = consumes<reco::TrackCollection>(conf_.getParameter<edm::InputTag>("TrackProducer"));
   referenceTrackProducer_ = consumes<reco::TrackCollection>(conf_.getParameter<edm::InputTag>("ReferenceTrackProducer"));
   jetCollection_ = mayConsume<reco::CaloJetCollection>(conf_.getParameter<edm::InputTag>("CaloJetCollection"));
 }

TkAlCaRecoMonitor::~TkAlCaRecoMonitor ( )

override

Definition at line 34 of file TkAlCaRecoMonitor.cc.

34 { }

Member Function Documentation

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

override

Definition at line 185 of file TkAlCaRecoMonitor.cc.

References AlCaRecoTrackEfficiency_, TtFullHadDaughter::B, binByRawId_, daughterMass_, deltaR(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, MillePedeFileConverter_cfg::e, MonitorElement::Fill(), fillHitmaps(), fillInvariantMass_, fillRawIdMap(), fillRawIdMap_, geometry, edm::EventSetup::get(), edm::Event::getByToken(), mps_fire::i, MagneticField::inTesla(), invariantMass_, edm::ESHandleBase::isValid(), edm::HandleBase::isValid(), jetCollection_, jetPt_, fwrapper::jets, seedCreatorFromRegionConsecutiveHitsEDProducer_cff::magneticField, maxJetPt_, ALCARECOTkAlMuonIsolated_cff::minJetDeltaR, minJetDeltaR_, minTrackDeltaR_, reco::TrackBase::qualitySize, referenceTrackProducer_, runsOnReco_, mathSSE::sqrt(), sumCharge_, HiIsolationCommonParameters_cff::track, findElectronsInSiStrips_cfi::trackCollection, TrackCurvature_, trackProducer_, TrackPtNegative_, TrackPtPositive_, TrackQuality_, and PV3DBase< T, PVType, FrameType >::z().

                                                                                    {
 
   edm::Handle<reco::TrackCollection> trackCollection;
   iEvent.getByToken(trackProducer_, trackCollection);
   if (!trackCollection.isValid()){
     edm::LogError("Alignment")<<"invalid trackcollection encountered!";
     return;
   }
 
   edm::Handle<reco::TrackCollection> referenceTrackCollection;
   iEvent.getByToken(referenceTrackProducer_, referenceTrackCollection);
   if (!trackCollection.isValid()){
     edm::LogError("Alignment")<<"invalid reference track-collection encountered!";
     return;
   }
 
   edm::ESHandle<TrackerGeometry> geometry;
   iSetup.get<TrackerDigiGeometryRecord>().get(geometry);
   if(! geometry.isValid()){
     edm::LogError("Alignment")<<"invalid geometry found in event setup!";
   }
 
   edm::ESHandle<MagneticField> magneticField;
   iSetup.get<IdealMagneticFieldRecord>().get(magneticField);
   if (!magneticField.isValid()){
     edm::LogError("Alignment")<<"invalid magnetic field configuration encountered!";
     return;
   }
 
   edm::Handle<reco::CaloJetCollection> jets;
   if(runsOnReco_){
     iEvent.getByToken(jetCollection_, jets);
     if(! jets.isValid()){
       edm::LogError("Alignment")<<"no jet collection found in event!";
     }
   }
   // fill only once - not yet in beginJob since no access to geometry
   if (fillRawIdMap_ && binByRawId_.empty()) this->fillRawIdMap(*geometry);
 
   AlCaRecoTrackEfficiency_->Fill( static_cast<double>((*trackCollection).size()) / (*referenceTrackCollection).size() );
   
   double sumOfCharges = 0;
   for( reco::TrackCollection::const_iterator track = (*trackCollection).begin(); track < (*trackCollection).end(); ++track ){
     double dR = 0;  
     if(runsOnReco_){
       double minJetDeltaR = 10; // some number > 2pi
       for(reco::CaloJetCollection::const_iterator itJet = jets->begin(); itJet != jets->end() ; ++itJet) {
     jetPt_->Fill( (*itJet).pt() );
     dR = deltaR((*track),(*itJet));
     if((*itJet).pt() > maxJetPt_ && dR < minJetDeltaR)
       minJetDeltaR = dR; 
     
     //edm::LogInfo("Alignment") <<">  isolated: "<< isolated << " jetPt "<< (*itJet).pt() <<" deltaR: "<< deltaR(*(*it),(*itJet)) ;
       }
       minJetDeltaR_->Fill( minJetDeltaR );
     }
     
     double minTrackDeltaR = 10; // some number > 2pi
     for( reco::TrackCollection::const_iterator track2 = (*trackCollection).begin(); track2 < (*trackCollection).end(); ++track2 ){
       dR = deltaR((*track),(*track2));
       if(dR < minTrackDeltaR && dR > 1e-6)
     minTrackDeltaR = dR;
     }
 
     for ( int i = 0; i<reco::TrackBase::qualitySize; ++i){
       if( (*track).quality( reco::TrackBase::TrackQuality(i) ) ){
     TrackQuality_->Fill( i );
       }
     }
     
     GlobalPoint gPoint((*track).vx(), (*track).vy(), (*track).vz());
     double B = magneticField->inTesla(gPoint).z();
     double curv = -(*track).charge()*0.002998*B/(*track).pt();
     TrackCurvature_->Fill( curv );
 
     if( (*track).charge() > 0 )
       TrackPtPositive_->Fill( (*track).pt() );
     if( (*track).charge() < 0 )
       TrackPtNegative_->Fill( (*track).pt() );
 
     minTrackDeltaR_->Fill( minTrackDeltaR );
     fillHitmaps( *track, *geometry );
     sumOfCharges += (*track).charge();
   }
 
   sumCharge_->Fill( sumOfCharges );
 
   if(fillInvariantMass_){
     if((*trackCollection).size() == 2){
       TLorentzVector track0((*trackCollection).at(0).px(),(*trackCollection).at(0).py(),(*trackCollection).at(0).pz(),
                 sqrt(((*trackCollection).at(0).p()*(*trackCollection).at(0).p())+daughterMass_*daughterMass_));
       TLorentzVector track1((*trackCollection).at(1).px(),(*trackCollection).at(1).py(),(*trackCollection).at(1).pz(),
                 sqrt(((*trackCollection).at(1).p()*(*trackCollection).at(1).p())+daughterMass_*daughterMass_));
       TLorentzVector mother = track0+track1;
       
       invariantMass_->Fill( mother.M() );
     }else{
       edm::LogInfo("Alignment")<<"wrong number of tracks trackcollection encountered: "<<(*trackCollection).size();
     }
   }
 }

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 36 of file TkAlCaRecoMonitor.cc.

References AlCaRecoTrackEfficiency_, DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), conf_, daughterMass_, fillInvariantMass_, fillRawIdMap_, edm::ParameterSet::getParameter(), MonitorElement::getTH1(), Hits_perDetId_, Hits_XvsY_, Hits_ZvsR_, mps_fire::i, invariantMass_, jetPt_, CaloRecoTauProducer_cfi::JetPtMin, edm::EDConsumerBase::labelsForToken(), maxJetPt_, minJetDeltaR_, minTrackDeltaR_, edm::ProductLabels::module, reco::TrackBase::qualityName(), reco::TrackBase::qualitySize, conversionPostprocessing_cfi::rBin, referenceTrackProducer_, conversionPostprocessing_cfi::rMax, conversionPostprocessing_cfi::rMin, runsOnReco_, MonitorElement::setAxisTitle(), DQMStore::IBooker::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, sumCharge_, TrackCurvature_, trackProducer_, TrackPtNegative_, TrackPtPositive_, TrackQuality_, useSignedR_, conversionPostprocessing_cfi::zBin, and conversionPostprocessing_cfi::zMax.

                                                                                                       {
 
   std::string histname;  //for naming the histograms according to algorithm used
 
   std::string AlgoName     = conf_.getParameter<std::string>("AlgoName");
   std::string MEFolderName = conf_.getParameter<std::string>("FolderName"); 
 
   daughterMass_ = conf_.getParameter<double>("daughterMass");
 
   maxJetPt_ = conf_.getParameter<double>("maxJetPt");
 
   iBooker.setCurrentFolder(MEFolderName+"/TkAlignmentSpecific");
   fillInvariantMass_ = conf_.getParameter<bool>("fillInvariantMass");
   runsOnReco_ = conf_.getParameter<bool>("runsOnReco");
   useSignedR_ = conf_.getParameter<bool>("useSignedR");
   fillRawIdMap_ = conf_.getParameter<bool>("fillRawIdMap");
 
   //    
   unsigned int MassBin = conf_.getParameter<unsigned int>("MassBin");
   double MassMin = conf_.getParameter<double>("MassMin");
   double MassMax = conf_.getParameter<double>("MassMax");
 
   if(fillInvariantMass_){
     histname = "InvariantMass_";
     invariantMass_ = iBooker.book1D(histname+AlgoName, histname+AlgoName, MassBin, MassMin, MassMax);
     invariantMass_->setAxisTitle("invariant Mass / GeV");
   }else{
     invariantMass_ = nullptr;
   }
 
   unsigned int TrackPtPositiveBin = conf_.getParameter<unsigned int>("TrackPtBin");
   double TrackPtPositiveMin = conf_.getParameter<double>("TrackPtMin");
   double TrackPtPositiveMax = conf_.getParameter<double>("TrackPtMax");
 
   histname = "TrackPtPositive_";
   TrackPtPositive_ = iBooker.book1D(histname+AlgoName, histname+AlgoName, TrackPtPositiveBin, TrackPtPositiveMin, TrackPtPositiveMax);
   TrackPtPositive_->setAxisTitle("p_{T} of tracks charge > 0");
 
   unsigned int TrackPtNegativeBin = conf_.getParameter<unsigned int>("TrackPtBin");
   double TrackPtNegativeMin = conf_.getParameter<double>("TrackPtMin");
   double TrackPtNegativeMax = conf_.getParameter<double>("TrackPtMax");
 
   histname = "TrackPtNegative_";
   TrackPtNegative_ = iBooker.book1D(histname+AlgoName, histname+AlgoName, TrackPtNegativeBin, TrackPtNegativeMin, TrackPtNegativeMax);
   TrackPtNegative_->setAxisTitle("p_{T} of tracks charge < 0");
 
   histname = "TrackQuality_";
   TrackQuality_ = iBooker.book1D(histname+AlgoName, histname+AlgoName, 
                  reco::TrackBase::qualitySize, -0.5, reco::TrackBase::qualitySize-0.5);
   TrackQuality_->setAxisTitle("quality");
   for ( int i = 0; i<reco::TrackBase::qualitySize; ++i){
     TrackQuality_->getTH1()->GetXaxis()->SetBinLabel(i+1,
                 reco::TrackBase::qualityName( reco::TrackBase::TrackQuality(i) ).c_str());
   } 
 
   unsigned int SumChargeBin = conf_.getParameter<unsigned int>("SumChargeBin");
   double SumChargeMin = conf_.getParameter<double>("SumChargeMin");
   double SumChargeMax = conf_.getParameter<double>("SumChargeMax");
 
   histname = "SumCharge_";
   sumCharge_ = iBooker.book1D(histname+AlgoName, histname+AlgoName, SumChargeBin, SumChargeMin, SumChargeMax);
   sumCharge_->setAxisTitle("#SigmaCharge");
 
   unsigned int TrackCurvatureBin = conf_.getParameter<unsigned int>("TrackCurvatureBin");
   double TrackCurvatureMin = conf_.getParameter<double>("TrackCurvatureMin");
   double TrackCurvatureMax = conf_.getParameter<double>("TrackCurvatureMax");
   
   histname = "TrackCurvature_";
   TrackCurvature_ = iBooker.book1D(histname+AlgoName, histname+AlgoName, TrackCurvatureBin, TrackCurvatureMin, TrackCurvatureMax);
   TrackCurvature_->setAxisTitle("#kappa track");
 
   if( runsOnReco_ ){
 
     unsigned int    JetPtBin = conf_.getParameter<unsigned int>("JetPtBin");
     double JetPtMin = conf_.getParameter<double>("JetPtMin");
     double JetPtMax = conf_.getParameter<double>("JetPtMax");
 
     histname = "JetPt_";
     jetPt_ = iBooker.book1D(histname+AlgoName, histname+AlgoName, JetPtBin, JetPtMin, JetPtMax);
     jetPt_->setAxisTitle("jet p_{T} / GeV");
 
     unsigned int    MinJetDeltaRBin = conf_.getParameter<unsigned int>("MinJetDeltaRBin");
     double MinJetDeltaRMin = conf_.getParameter<double>("MinJetDeltaRMin");
     double MinJetDeltaRMax = conf_.getParameter<double>("MinJetDeltaRMax");
 
     histname = "MinJetDeltaR_";
     minJetDeltaR_ = iBooker.book1D(histname+AlgoName, histname+AlgoName, MinJetDeltaRBin, MinJetDeltaRMin, MinJetDeltaRMax);
     minJetDeltaR_->setAxisTitle("minimal Jet #DeltaR / rad");
   }else{
     jetPt_ = nullptr;
     minJetDeltaR_ = nullptr;
   }
 
   unsigned int    MinTrackDeltaRBin = conf_.getParameter<unsigned int>("MinTrackDeltaRBin");
   double MinTrackDeltaRMin = conf_.getParameter<double>("MinTrackDeltaRMin");
   double MinTrackDeltaRMax = conf_.getParameter<double>("MinTrackDeltaRMax");
 
   histname = "MinTrackDeltaR_";
   minTrackDeltaR_ = iBooker.book1D(histname+AlgoName, histname+AlgoName, MinTrackDeltaRBin, MinTrackDeltaRMin, MinTrackDeltaRMax);
   minTrackDeltaR_->setAxisTitle("minimal Track #DeltaR / rad");
 
   unsigned int TrackEfficiencyBin = conf_.getParameter<unsigned int>("TrackEfficiencyBin");
   double TrackEfficiencyMin = conf_.getParameter<double>("TrackEfficiencyMin");
   double TrackEfficiencyMax = conf_.getParameter<double>("TrackEfficiencyMax");
 
   histname = "AlCaRecoTrackEfficiency_";
   AlCaRecoTrackEfficiency_ = iBooker.book1D(histname+AlgoName, histname+AlgoName, TrackEfficiencyBin, TrackEfficiencyMin, TrackEfficiencyMax);
   Labels l_tp, l_rtp;
   labelsForToken(referenceTrackProducer_, l_rtp);
   labelsForToken(trackProducer_, l_tp);
   AlCaRecoTrackEfficiency_->setAxisTitle("n("+std::string(l_tp.module)+") / n("+std::string(l_rtp.module)+")");
 
   int zBin =  conf_.getParameter<unsigned int>("HitMapsZBin"); //300
   double zMax = conf_.getParameter<double>("HitMapZMax"); //300.0; //cm
   
   int rBin = conf_.getParameter<unsigned int>("HitMapsRBin");//120;
   double rMax = conf_.getParameter<double>("HitMapRMax"); //120.0; //cm
 
   histname = "Hits_ZvsR_";
   double rMin = 0.0;
   if( useSignedR_ )
     rMin = -rMax;
 
   Hits_ZvsR_ = iBooker.book2D(histname+AlgoName, histname+AlgoName, zBin, -zMax, zMax, rBin, rMin, rMax);
 
   histname = "Hits_XvsY_";
   Hits_XvsY_ = iBooker.book2D(histname+AlgoName, histname+AlgoName, rBin, -rMax, rMax, rBin, -rMax, rMax);
 
   if( fillRawIdMap_){
     histname = "Hits_perDetId_";
     
     //leads to differences in axsis between samples??
     //int nModules = binByRawId_.size();
     //Hits_perDetId_ = iBooker.book1D(histname+AlgoName, histname+AlgoName, nModules, static_cast<double>(nModules) -0.5, static_cast<double>(nModules) -0.5);
     Hits_perDetId_ = iBooker.book1D(histname+AlgoName, histname+AlgoName, 16601,-0.5, 16600.5 );
     Hits_perDetId_->setAxisTitle("rawId Bins");
 
     //  std::stringstream binLabel;
     //  for( std::map<int,int>::iterator it = binByRawId_.begin(); it != binByRawId_.end(); ++it ){
     //    binLabel.str() = "";
     //    binLabel << (*it).first;
     //    Hits_perDetId_->getTH1()->GetXaxis()->SetBinLabel( (*it).second +1, binLabel.str().c_str());
     //  }
   }
 }

void TkAlCaRecoMonitor::fillHitmaps	(	const reco::Track &	track,
		const TrackerGeometry &	geometry
	)

private

Definition at line 287 of file TkAlCaRecoMonitor.cc.

References binByRawId_, MonitorElement::Fill(), fillRawIdMap_, TrackingRecHit::geographicalId(), Hits_perDetId_, Hits_XvsY_, Hits_ZvsR_, TrackerGeometry::idToDet(), alignCSCRings::r, reco::Track::recHitsBegin(), reco::Track::recHitsEnd(), mathSSE::sqrt(), and useSignedR_.

Referenced by analyze().

 {
   for (trackingRecHit_iterator iHit = track.recHitsBegin(); iHit != track.recHitsEnd(); ++iHit) {    
     if( (*iHit)->isValid() ){
       const TrackingRecHit *hit = (*iHit);
       const DetId geoId(hit->geographicalId());
       const GeomDet* gd = geometry.idToDet(geoId);
       // since 2_1_X local hit positions are transient. taking center of the hit module for now.
       // The alternative would be the coarse estimation or a refit.
       //const GlobalPoint globP( gd->toGlobal( hit->localPosition() ) );
       const GlobalPoint globP( gd->toGlobal( Local3DPoint(0.,0.,0.) ) );
       double r = sqrt( globP.x()*globP.x() + globP.y()*globP.y() );
       if( useSignedR_ )
     r*= globP.y() / fabs( globP.y() );
 
       Hits_ZvsR_->Fill( globP.z(), r );
       Hits_XvsY_->Fill( globP.x(), globP.y() );
 
       if( fillRawIdMap_)
      Hits_perDetId_->Fill( binByRawId_[ geoId.rawId() ]);  
     }
   }
 }

void TkAlCaRecoMonitor::fillRawIdMap ( const TrackerGeometry & geometry )

private

Definition at line 311 of file TkAlCaRecoMonitor.cc.

References binByRawId_, DEFINE_FWK_MODULE, TrackerGeometry::detUnitIds(), and mps_fire::i.

Referenced by analyze().

 {
   std::vector<int> sortedRawIds;
   for (std::vector<DetId>::const_iterator iDetId = geometry.detUnitIds().begin();
        iDetId != geometry.detUnitIds().end(); ++iDetId) {
     sortedRawIds.push_back((*iDetId).rawId());
   }
   std::sort(sortedRawIds.begin(), sortedRawIds.end());
   
   int i = 0;
   for (std::vector<int>::iterator iRawId = sortedRawIds.begin();
        iRawId != sortedRawIds.end(); ++iRawId){
     binByRawId_[*iRawId] = i;
     ++i;
   }
 }