Phase2ITValidateRecHitBase Class Reference

#include <Phase2ITValidateRecHitBase.h>

Inheritance diagram for Phase2ITValidateRecHitBase:

Classes
struct	RecHitME

Public Member Functions
void	bookHistograms (DQMStore::IBooker &ibooker, edm::Run const &iRun, edm::EventSetup const &iSetup) override
void	dqmBeginRun (const edm::Run &, const edm::EventSetup &) override
	Phase2ITValidateRecHitBase (const edm::ParameterSet &)
	~Phase2ITValidateRecHitBase () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
virtual void	analyze (edm::Event const &, edm::EventSetup const &)
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
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Static Public Member Functions
static void	fillPSetDescription (edm::ParameterSetDescription &desc)
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 &)

Protected Member Functions
void	bookLayerHistos (DQMStore::IBooker &ibooker, unsigned int det_id, std::string &subdir)
void	fillRechitHistos (const PSimHit *simhitClosest, const SiPixelRecHit *rechit, const std::map< unsigned int, SimTrack > &selectedSimTrackMap, std::map< std::string, unsigned int > &nrechitLayerMap_primary)
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Protected Attributes
edm::ParameterSet	config_
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	geomToken_
std::map< std::string, RecHitME >	layerMEs_
const TrackerGeometry *	tkGeom_ = nullptr
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	topoToken_
const TrackerTopology *	tTopo_ = nullptr
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

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 >
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Description: Base Class for Phase2 Validation

Definition at line 41 of file Phase2ITValidateRecHitBase.h.

Constructor & Destructor Documentation

Phase2ITValidateRecHitBase()

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

explicit

Definition at line 18 of file Phase2ITValidateRecHitBase.cc.

    : config_(iConfig),
      geomToken_(esConsumes<TrackerGeometry, TrackerDigiGeometryRecord, edm::Transition::BeginRun>()),
      topoToken_(esConsumes<TrackerTopology, TrackerTopologyRcd, edm::Transition::BeginRun>()) {
  edm::LogInfo("Phase2ITValidateRecHitBase") << ">>> Construct Phase2ITValidateRecHitBase ";
}

~Phase2ITValidateRecHitBase()

Phase2ITValidateRecHitBase::~Phase2ITValidateRecHitBase ( )

overridedefault

Member Function Documentation

bookHistograms()

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 32 of file Phase2ITValidateRecHitBase.cc.

References bookLayerHistos(), edm::ESWatcher< T >::check(), config_, TrackerGeometry::detUnits(), edm::ParameterSet::getParameter(), GeomDetEnumerators::P2PXB, GeomDetEnumerators::P2PXEC, AlCaHLTBitMon_QueryRunRegistry::string, and tkGeom_.

                                                                             {
  std::string top_folder = config_.getParameter<std::string>("TopFolderName");
  edm::LogInfo("Phase2ITValidateRecHitBase") << " Booking Histograms in : " << top_folder;
  edm::ESWatcher<TrackerDigiGeometryRecord> theTkDigiGeomWatcher;
  if (theTkDigiGeomWatcher.check(iSetup)) {
    for (auto const& det_u : tkGeom_->detUnits()) {
      //Always check TrackerNumberingBuilder before changing this part
      if (!(det_u->subDetector() == GeomDetEnumerators::SubDetector::P2PXB ||
            det_u->subDetector() == GeomDetEnumerators::SubDetector::P2PXEC))
        continue;
      unsigned int detId_raw = det_u->geographicalId().rawId();
      bookLayerHistos(ibooker, detId_raw, top_folder);
    }
  }
}

bookLayerHistos()

void Phase2ITValidateRecHitBase::bookLayerHistos ( DQMStore::IBooker &  ibooker, unsigned int  det_id, std::string &  subdir  )

protected

Definition at line 51 of file Phase2ITValidateRecHitBase.cc.

References phase2tkutil::book1DFromPSet(), phase2tkutil::book2DFromPSet(), phase2tkutil::bookProfile1DFromPSet(), dqm::implementation::NavigatorBase::cd(), config_, Phase2ITValidateRecHitBase::RecHitME::deltaX, Phase2ITValidateRecHitBase::RecHitME::deltaX_clsizex, Phase2ITValidateRecHitBase::RecHitME::deltaX_clsizey, Phase2ITValidateRecHitBase::RecHitME::deltaX_eta, Phase2ITValidateRecHitBase::RecHitME::deltaX_phi, Phase2ITValidateRecHitBase::RecHitME::deltaX_primary, Phase2ITValidateRecHitBase::RecHitME::deltaY, Phase2ITValidateRecHitBase::RecHitME::deltaY_clsizex, Phase2ITValidateRecHitBase::RecHitME::deltaY_clsizey, Phase2ITValidateRecHitBase::RecHitME::deltaY_eta, Phase2ITValidateRecHitBase::RecHitME::deltaY_phi, Phase2ITValidateRecHitBase::RecHitME::deltaY_primary, Phase2ITValidateRecHitBase::RecHitME::deltaYvsdeltaX, phase2tkutil::getITHistoId(), edm::ParameterSet::getParameter(), crabWrapper::key, layerMEs_, Phase2ITValidateRecHitBase::RecHitME::numberRecHitsprimary, Phase2ITValidateRecHitBase::RecHitME::pullX, Phase2ITValidateRecHitBase::RecHitME::pullX_eta, Phase2ITValidateRecHitBase::RecHitME::pullX_primary, Phase2ITValidateRecHitBase::RecHitME::pullY, Phase2ITValidateRecHitBase::RecHitME::pullY_eta, Phase2ITValidateRecHitBase::RecHitME::pullY_primary, dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, mkLumiAveragedPlots::subdir, and tTopo_.

Referenced by bookHistograms().

                                                                                                                 {
  ibooker.cd();
  std::string key = phase2tkutil::getITHistoId(det_id, tTopo_);
  if (key.empty())
    return;
  if (layerMEs_.find(key) == layerMEs_.end()) {
    ibooker.cd();
    RecHitME local_histos;
    ibooker.setCurrentFolder(subdir + "/" + key);
    edm::LogInfo("Phase2ITValidateRecHit") << " Booking Histograms in : " << (subdir + "/" + key);

    local_histos.deltaX = phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("DeltaX"), ibooker);

    local_histos.deltaY = phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("DeltaY"), ibooker);

    local_histos.pullX = phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("PullX"), ibooker);

    local_histos.pullY = phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("PullY"), ibooker);

    local_histos.deltaX_eta =
        phase2tkutil::book2DFromPSet(config_.getParameter<edm::ParameterSet>("DeltaX_eta"), ibooker);

    local_histos.deltaX_phi =
        phase2tkutil::book2DFromPSet(config_.getParameter<edm::ParameterSet>("DeltaX_phi"), ibooker);

    local_histos.deltaY_eta =
        phase2tkutil::book2DFromPSet(config_.getParameter<edm::ParameterSet>("DeltaY_eta"), ibooker);

    local_histos.deltaY_phi =
        phase2tkutil::book2DFromPSet(config_.getParameter<edm::ParameterSet>("DeltaY_phi"), ibooker);

    local_histos.deltaX_clsizex =
        phase2tkutil::bookProfile1DFromPSet(config_.getParameter<edm::ParameterSet>("DeltaX_clsizex"), ibooker);

    local_histos.deltaX_clsizey =
        phase2tkutil::bookProfile1DFromPSet(config_.getParameter<edm::ParameterSet>("DeltaX_clsizey"), ibooker);

    local_histos.deltaY_clsizex =
        phase2tkutil::bookProfile1DFromPSet(config_.getParameter<edm::ParameterSet>("DeltaY_clsizex"), ibooker);

    local_histos.deltaY_clsizey =
        phase2tkutil::bookProfile1DFromPSet(config_.getParameter<edm::ParameterSet>("DeltaY_clsizey"), ibooker);

    local_histos.deltaYvsdeltaX =
        phase2tkutil::book2DFromPSet(config_.getParameter<edm::ParameterSet>("DeltaY_vs_DeltaX"), ibooker);

    local_histos.pullX_eta =
        phase2tkutil::bookProfile1DFromPSet(config_.getParameter<edm::ParameterSet>("PullX_eta"), ibooker);

    local_histos.pullY_eta =
        phase2tkutil::bookProfile1DFromPSet(config_.getParameter<edm::ParameterSet>("PullY_eta"), ibooker);
    ibooker.setCurrentFolder(subdir + "/" + key + "/PrimarySimHits");
    //all histos for Primary particles
    local_histos.numberRecHitsprimary =
        phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("nRecHits_primary"), ibooker);

    local_histos.deltaX_primary =
        phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("DeltaX_primary"), ibooker);

    local_histos.deltaY_primary =
        phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("DeltaY_primary"), ibooker);

    local_histos.pullX_primary =
        phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("PullX_primary"), ibooker);

    local_histos.pullY_primary =
        phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("PullY_primary"), ibooker);

    layerMEs_.emplace(key, local_histos);
  }
}

dqmBeginRun()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 25 of file Phase2ITValidateRecHitBase.cc.

References geomToken_, edm::EventSetup::getData(), tkGeom_, topoToken_, and tTopo_.

                                                                                            {
  tkGeom_ = &iSetup.getData(geomToken_);
  tTopo_ = &iSetup.getData(topoToken_);
}

fillPSetDescription()

void Phase2ITValidateRecHitBase::fillPSetDescription ( edm::ParameterSetDescription &  desc )

static

Definition at line 179 of file Phase2ITValidateRecHitBase.cc.

References edm::ParameterSetDescription::add(), submitPVResolutionJobs::desc, M_PI, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by Phase2ITValidateRecHit::fillDescriptions(), and Phase2ITValidateTrackingRecHit::fillDescriptions().

                                                                                     {
  edm::ParameterSetDescription psd0;
  psd0.add<std::string>("name", "Delta_X");
  psd0.add<std::string>("title", "Delta_X;RecHit resolution X coordinate [#mum]");
  psd0.add<bool>("switch", true);
  psd0.add<double>("xmin", -100.0);
  psd0.add<double>("xmax", 100.0);
  psd0.add<int>("NxBins", 100);
  desc.add<edm::ParameterSetDescription>("DeltaX", psd0);

  edm::ParameterSetDescription psd1;
  psd1.add<std::string>("name", "Delta_Y");
  psd1.add<std::string>("title", "Delta_Y;RecHit resolution Y coordinate [#mum];");
  psd1.add<bool>("switch", true);
  psd1.add<double>("xmin", -100.0);
  psd1.add<double>("xmax", 100.0);
  psd1.add<int>("NxBins", 100);
  desc.add<edm::ParameterSetDescription>("DeltaY", psd1);

  edm::ParameterSetDescription psd2;
  psd2.add<std::string>("name", "Pull_X");
  psd2.add<std::string>("title", "Pull_X;pull x;");
  psd2.add<double>("xmin", -4.0);
  psd2.add<bool>("switch", true);
  psd2.add<double>("xmax", 4.0);
  psd2.add<int>("NxBins", 100);
  desc.add<edm::ParameterSetDescription>("PullX", psd2);

  edm::ParameterSetDescription psd3;
  psd3.add<std::string>("name", "Pull_Y");
  psd3.add<std::string>("title", "Pull_Y;pull y;");
  psd3.add<double>("xmin", -4.0);
  psd3.add<bool>("switch", true);
  psd3.add<double>("xmax", 4.0);
  psd3.add<int>("NxBins", 100);
  desc.add<edm::ParameterSetDescription>("PullY", psd3);

  edm::ParameterSetDescription psd4;
  psd4.add<std::string>("name", "Delta_X_vs_Eta");
  psd4.add<std::string>("title", "Delta_X_vs_Eta;|#eta|;#Delta x [#mum]");
  psd4.add<int>("NyBins", 100);
  psd4.add<double>("ymin", -100.0);
  psd4.add<double>("ymax", 100.0);
  psd4.add<int>("NxBins", 41);
  psd4.add<bool>("switch", true);
  psd4.add<double>("xmax", 4.1);
  psd4.add<double>("xmin", 0.);
  desc.add<edm::ParameterSetDescription>("DeltaX_eta", psd4);

  edm::ParameterSetDescription psd4_y;
  psd4_y.add<std::string>("name", "Delta_X_vs_Phi");
  ;
  psd4_y.add<std::string>("title", "Delta_X_vs_Phi;#phi;#Delta x [#mum]");
  psd4_y.add<int>("NyBins", 100);
  psd4_y.add<double>("ymin", -100.0);
  psd4_y.add<double>("ymax", 100.0);
  psd4_y.add<int>("NxBins", 36);
  psd4_y.add<bool>("switch", true);
  psd4_y.add<double>("xmax", M_PI);
  psd4_y.add<double>("xmin", -M_PI);
  desc.add<edm::ParameterSetDescription>("DeltaX_phi", psd4_y);

  edm::ParameterSetDescription psd5;
  psd5.add<std::string>("name", "Delta_Y_vs_Eta");
  psd5.add<std::string>("title", "Delta_Y_vs_Eta;|#eta|;#Delta y [#mum]");
  psd5.add<int>("NyBins", 100);
  psd5.add<double>("ymin", -100.0);
  psd5.add<double>("ymax", 100.0);
  psd5.add<int>("NxBins", 41);
  psd5.add<bool>("switch", true);
  psd5.add<double>("xmax", 4.1);
  psd5.add<double>("xmin", 0.);
  desc.add<edm::ParameterSetDescription>("DeltaY_eta", psd5);

  edm::ParameterSetDescription psd5_y;
  psd5_y.add<std::string>("name", "Delta_Y_vs_Phi");
  psd5_y.add<std::string>("title", "Delta_Y_vs_Phi;#phi;#Delta y [#mum]");
  psd5_y.add<int>("NyBins", 100);
  psd5_y.add<double>("ymin", -100.0);
  psd5_y.add<double>("ymax", 100.0);
  psd5_y.add<int>("NxBins", 36);
  psd5_y.add<bool>("switch", true);
  psd5_y.add<double>("xmax", M_PI);
  psd5_y.add<double>("xmin", -M_PI);
  desc.add<edm::ParameterSetDescription>("DeltaY_phi", psd5_y);

  edm::ParameterSetDescription psd6;
  psd6.add<std::string>("name", "Delta_X_vs_ClusterSizeX");
  psd6.add<std::string>("title", ";Cluster size X;#Delta x [#mum]");
  psd6.add<double>("ymin", -100.0);
  psd6.add<double>("ymax", 100.0);
  psd6.add<int>("NxBins", 21);
  psd6.add<bool>("switch", true);
  psd6.add<double>("xmax", 20.5);
  psd6.add<double>("xmin", -0.5);
  desc.add<edm::ParameterSetDescription>("DeltaX_clsizex", psd6);

  edm::ParameterSetDescription psd7;
  psd7.add<std::string>("name", "Delta_X_vs_ClusterSizeY");
  psd7.add<std::string>("title", ";Cluster size Y;#Delta x [#mum]");
  psd7.add<double>("ymin", -100.0);
  psd7.add<double>("ymax", 100.0);
  psd7.add<int>("NxBins", 21);
  psd7.add<bool>("switch", true);
  psd7.add<double>("xmax", 20.5);
  psd7.add<double>("xmin", -0.5);
  desc.add<edm::ParameterSetDescription>("DeltaX_clsizey", psd7);

  edm::ParameterSetDescription psd8;
  psd8.add<std::string>("name", "Delta_Y_vs_ClusterSizeX");
  psd8.add<std::string>("title", ";Cluster size X;#Delta y [#mum]");
  psd8.add<double>("ymin", -100.0);
  psd8.add<double>("ymax", 100.0);
  psd8.add<int>("NxBins", 21);
  psd8.add<bool>("switch", true);
  psd8.add<double>("xmax", 20.5);
  psd8.add<double>("xmin", -0.5);
  desc.add<edm::ParameterSetDescription>("DeltaY_clsizex", psd8);

  edm::ParameterSetDescription psd9;
  psd9.add<std::string>("name", "Delta_Y_vs_ClusterSizeY");
  psd9.add<std::string>("title", ";Cluster size Y;#Delta y [#mum]");
  psd9.add<double>("ymin", -100.0);
  psd9.add<double>("ymax", 100.0);
  psd9.add<int>("NxBins", 21);
  psd9.add<bool>("switch", true);
  psd9.add<double>("xmax", 20.5);
  psd9.add<double>("xmin", -0.5);
  desc.add<edm::ParameterSetDescription>("DeltaY_clsizey", psd9);

  edm::ParameterSetDescription psd10;
  psd10.add<std::string>("name", "Delta_Y_vs_DeltaX");
  psd10.add<std::string>("title", ";#Delta x[#mum];#Delta y[#mum]");
  psd10.add<bool>("switch", true);
  psd10.add<double>("ymin", -100.0);
  psd10.add<double>("ymax", 100.0);
  psd10.add<int>("NyBins", 100);
  psd10.add<double>("xmax", 100.);
  psd10.add<double>("xmin", -100.);
  psd10.add<int>("NxBins", 100);
  desc.add<edm::ParameterSetDescription>("DeltaY_vs_DeltaX", psd10);

  edm::ParameterSetDescription psd11;
  psd11.add<std::string>("name", "Pull_X_vs_Eta");
  psd11.add<std::string>("title", "Pull_X_vs_Eta;#eta;pull x");
  psd11.add<double>("ymax", 4.0);
  psd11.add<int>("NxBins", 82);
  psd11.add<bool>("switch", true);
  psd11.add<double>("xmax", 4.1);
  psd11.add<double>("xmin", -4.1);
  psd11.add<double>("ymin", -4.0);
  desc.add<edm::ParameterSetDescription>("PullX_eta", psd11);

  edm::ParameterSetDescription psd12;
  psd12.add<std::string>("name", "Pull_Y_vs_Eta");
  psd12.add<std::string>("title", "Pull_Y_vs_Eta;#eta;pull y");
  psd12.add<double>("ymax", 4.0);
  psd12.add<int>("NxBins", 82);
  psd12.add<bool>("switch", true);
  psd12.add<double>("xmax", 4.1);
  psd12.add<double>("xmin", -4.1);
  psd12.add<double>("ymin", -4.0);
  desc.add<edm::ParameterSetDescription>("PullY_eta", psd12);

  //simhits primary

  edm::ParameterSetDescription psd13;
  psd13.add<std::string>("name", "Number_RecHits_matched_PrimarySimTrack");
  psd13.add<std::string>("title", "Number of RecHits matched to primary SimTrack;;");
  psd13.add<double>("xmin", 0.0);
  psd13.add<bool>("switch", true);
  psd13.add<double>("xmax", 0.0);
  psd13.add<int>("NxBins", 100);
  desc.add<edm::ParameterSetDescription>("nRecHits_primary", psd13);

  edm::ParameterSetDescription psd14;
  psd14.add<std::string>("name", "Delta_X_SimHitPrimary");
  psd14.add<std::string>("title", "Delta_X_SimHitPrimary;#delta x [#mum];");
  psd14.add<double>("xmin", -100.0);
  psd14.add<bool>("switch", true);
  psd14.add<double>("xmax", 100.0);
  psd14.add<int>("NxBins", 100);
  desc.add<edm::ParameterSetDescription>("DeltaX_primary", psd14);

  edm::ParameterSetDescription psd15;
  psd15.add<std::string>("name", "Delta_Y_SimHitPrimary");
  psd15.add<std::string>("title", "Delta_Y_SimHitPrimary;#Delta y [#mum];");
  psd15.add<double>("xmin", -100.0);
  psd15.add<bool>("switch", true);
  psd15.add<double>("xmax", 100.0);
  psd15.add<int>("NxBins", 100);
  desc.add<edm::ParameterSetDescription>("DeltaY_primary", psd15);

  edm::ParameterSetDescription psd16;
  psd16.add<std::string>("name", "Pull_X_SimHitPrimary");
  psd16.add<std::string>("title", "Pull_X_SimHitPrimary;pull x;");
  psd16.add<double>("ymax", 4.0);
  psd16.add<int>("NxBins", 82);
  psd16.add<bool>("switch", true);
  psd16.add<double>("xmax", 4.1);
  psd16.add<double>("xmin", -4.1);
  psd16.add<double>("ymin", -4.0);
  desc.add<edm::ParameterSetDescription>("PullX_primary", psd16);

  edm::ParameterSetDescription psd17;
  psd17.add<std::string>("name", "Pull_Y_SimHitPrimary");
  psd17.add<std::string>("title", "Pull_Y_SimHitPrimary;pull y;");
  psd17.add<double>("ymax", 4.0);
  psd17.add<int>("NxBins", 82);
  psd17.add<bool>("switch", true);
  psd17.add<double>("xmax", 4.1);
  psd17.add<double>("xmin", -4.1);
  psd17.add<double>("ymin", -4.0);
  desc.add<edm::ParameterSetDescription>("PullY_primary", psd17);
}

fillRechitHistos()

void Phase2ITValidateRecHitBase::fillRechitHistos ( const PSimHit *  simhitClosest, const SiPixelRecHit *  rechit, const std::map< unsigned int, SimTrack > &  selectedSimTrackMap, std::map< std::string, unsigned int > &  nrechitLayerMap_primary  )

protected

Definition at line 123 of file Phase2ITValidateRecHitBase.cc.

References funct::abs(), SiPixelRecHit::cluster(), phase2tkutil::cmtomicron, PVValHelper::dx, PVValHelper::dy, PVValHelper::eta, TrackingRecHit::geographicalId(), phase2tkutil::getITHistoId(), TrackerGeometry::idToDetUnit(), phase2tkutil::isPrimary(), crabWrapper::key, layerMEs_, PSimHit::localPosition(), BaseTrackerRecHit::localPosition(), BaseTrackerRecHit::localPositionError(), phi, mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, tkGeom_, PSimHit::trackId(), tTopo_, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by Phase2ITValidateRecHit::fillITHistos(), and Phase2ITValidateTrackingRecHit::fillITHistos().

                                                                                                            {
  auto id = rechit->geographicalId();
  std::string key = phase2tkutil::getITHistoId(id.rawId(), tTopo_);
  const GeomDetUnit* geomDetunit(tkGeom_->idToDetUnit(id));
  if (!geomDetunit)
    return;

  LocalPoint lp = rechit->localPosition();
  auto simTrackIt(selectedSimTrackMap.find(simhitClosest->trackId()));
  bool isPrimary = false;
  //check if simhit is primary
  if (simTrackIt != selectedSimTrackMap.end())
    isPrimary = phase2tkutil::isPrimary(simTrackIt->second, simhitClosest);

  Local3DPoint simlp(simhitClosest->localPosition());
  const LocalError& lperr = rechit->localPositionError();
  double dx = phase2tkutil::cmtomicron * (lp.x() - simlp.x());
  double dy = phase2tkutil::cmtomicron * (lp.y() - simlp.y());
  double pullx = 999.;
  double pully = 999.;
  if (lperr.xx())
    pullx = (lp.x() - simlp.x()) / std::sqrt(lperr.xx());
  if (lperr.yy())
    pully = (lp.y() - simlp.y()) / std::sqrt(lperr.yy());
  float eta = geomDetunit->surface().toGlobal(lp).eta();
  float phi = geomDetunit->surface().toGlobal(lp).phi();
  layerMEs_[key].deltaX->Fill(dx);
  layerMEs_[key].deltaY->Fill(dy);
  layerMEs_[key].pullX->Fill(pullx);
  layerMEs_[key].pullY->Fill(pully);

  layerMEs_[key].deltaX_eta->Fill(std::abs(eta), dx);
  layerMEs_[key].deltaY_eta->Fill(std::abs(eta), dy);
  layerMEs_[key].deltaX_phi->Fill(phi, dx);
  layerMEs_[key].deltaY_phi->Fill(phi, dy);

  layerMEs_[key].deltaX_clsizex->Fill(rechit->cluster()->sizeX(), dx);
  layerMEs_[key].deltaX_clsizey->Fill(rechit->cluster()->sizeY(), dx);
  layerMEs_[key].deltaY_clsizex->Fill(rechit->cluster()->sizeX(), dy);
  layerMEs_[key].deltaY_clsizey->Fill(rechit->cluster()->sizeY(), dy);
  layerMEs_[key].deltaYvsdeltaX->Fill(dx, dy);
  layerMEs_[key].pullX_eta->Fill(eta, pullx);
  layerMEs_[key].pullY_eta->Fill(eta, pully);
  if (isPrimary) {
    layerMEs_[key].deltaX_primary->Fill(dx);
    layerMEs_[key].deltaY_primary->Fill(dy);
    layerMEs_[key].pullX_primary->Fill(pullx);
    layerMEs_[key].pullY_primary->Fill(pully);
  } else {
    nrechitLayerMap_primary[key]--;
  }
}

Member Data Documentation

config_

edm::ParameterSet Phase2ITValidateRecHitBase::config_

protected

Definition at line 56 of file Phase2ITValidateRecHitBase.h.

Referenced by bookHistograms(), and bookLayerHistos().

geomToken_

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> Phase2ITValidateRecHitBase::geomToken_

protected

Definition at line 57 of file Phase2ITValidateRecHitBase.h.

Referenced by dqmBeginRun().

layerMEs_

std::map<std::string, RecHitME> Phase2ITValidateRecHitBase::layerMEs_

protected

Definition at line 85 of file Phase2ITValidateRecHitBase.h.

Referenced by bookLayerHistos(), Phase2ITValidateRecHit::fillITHistos(), Phase2ITValidateTrackingRecHit::fillITHistos(), and fillRechitHistos().

tkGeom_

const TrackerGeometry* Phase2ITValidateRecHitBase::tkGeom_ = nullptr

protected

Definition at line 59 of file Phase2ITValidateRecHitBase.h.

Referenced by bookHistograms(), dqmBeginRun(), Phase2ITValidateTrackingRecHit::fillITHistos(), and fillRechitHistos().

topoToken_

const edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> Phase2ITValidateRecHitBase::topoToken_

protected

Definition at line 58 of file Phase2ITValidateRecHitBase.h.

Referenced by dqmBeginRun().

tTopo_

const TrackerTopology* Phase2ITValidateRecHitBase::tTopo_ = nullptr

protected

Definition at line 60 of file Phase2ITValidateRecHitBase.h.

Referenced by bookLayerHistos(), dqmBeginRun(), Phase2ITValidateRecHit::fillITHistos(), Phase2ITValidateTrackingRecHit::fillITHistos(), and fillRechitHistos().