#include <Phase2OTMonitorVectorHits.cc>

Inheritance diagram for Phase2OTMonitorVectorHits:

Classes
struct	VecHitME

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

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

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

	Phase2OTMonitorVectorHits (const edm::ParameterSet &)

	~Phase2OTMonitorVectorHits () override

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

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	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

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

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

Private Attributes
edm::ParameterSet	config_

const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	geomToken_

MonitorElement *	globalRZ_P_ = nullptr

MonitorElement *	globalRZ_S_ = nullptr

MonitorElement *	globalXY_P_ = nullptr

MonitorElement *	globalXY_S_ = nullptr

std::map< std::string, VecHitME >	layerMEs_

const MagneticField *	magField_ = nullptr

const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	magFieldToken_

MonitorElement *	numberVecHits_ = nullptr

const TrackerGeometry *	tkGeom_ = nullptr

const edm::EDGetTokenT< VectorHitCollection >	tokenVecHitsOT_

const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	topoToken_

const TrackerTopology *	tTopo_ = nullptr

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

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

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

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Description: DQM module to monitor Phase2 VectorHit

Definition at line 44 of file Phase2OTMonitorVectorHits.cc.

Constructor & Destructor Documentation

◆ Phase2OTMonitorVectorHits()

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

explicit

Definition at line 98 of file Phase2OTMonitorVectorHits.cc.

     : config_(iConfig),
       tokenVecHitsOT_(consumes<VectorHitCollection>(config_.getParameter<edm::InputTag>("vechitsSrc"))),
       geomToken_(esConsumes<TrackerGeometry, TrackerDigiGeometryRecord, edm::Transition::BeginRun>()),
       topoToken_(esConsumes<TrackerTopology, TrackerTopologyRcd, edm::Transition::BeginRun>()),
       magFieldToken_(esConsumes<MagneticField, IdealMagneticFieldRecord, edm::Transition::BeginRun>()) {
   edm::LogInfo("Phase2OTMonitorVectorHits") << ">>> Construct Phase2OTMonitorVectorHits ";
 }

◆ ~Phase2OTMonitorVectorHits()

Phase2OTMonitorVectorHits::~Phase2OTMonitorVectorHits ( )

override

Definition at line 110 of file Phase2OTMonitorVectorHits.cc.

                                                       {
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
   edm::LogInfo("Phase2OTMonitorVectorHits") << ">>> Destroy Phase2OTMonitorVectorHits ";
 }

Member Function Documentation

◆ analyze()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 126 of file Phase2OTMonitorVectorHits.cc.

References edmNew::DetSet< T >::begin(), PixelRecoUtilities::curvature(), edmNew::DetSet< T >::end(), PVValHelper::eta, PV3DBase< T, PVType, FrameType >::eta(), f, dqm::impl::MonitorElement::Fill(), TrackerGeometry::getDetectorType(), phase2tkutil::getOTHistoId(), globalRZ_P_, globalRZ_S_, globalXY_P_, globalXY_S_, iEvent, MagneticField::inTesla(), crabWrapper::key, layerMEs_, magField_, numberVecHits_, PV3DBase< T, PVType, FrameType >::perp(), TrackerGeometry::Ph2PSP, TrackerGeometry::Ph2SS, phi, PV3DBase< T, PVType, FrameType >::phi(), Validation_hcalonly_cfi::sign, AlCaHLTBitMon_QueryRunRegistry::string, tkGeom_, tokenVecHitsOT_, tTopo_, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                            {
   // Get the VecHits
   const auto& vechits = iEvent.getHandle(tokenVecHitsOT_);
   if (!vechits.isValid())
     return;
   std::map<std::string, unsigned int> nvechitLayerMapP;
   std::map<std::string, unsigned int> nvechitLayerMapS;
   unsigned long int nTotvechitsinevt = 0;
   // Loop over modules
   VectorHitCollection::const_iterator DSViter;
   for (DSViter = vechits->begin(); DSViter != vechits->end(); ++DSViter) {
     // Get the detector unit's id
     unsigned int rawid(DSViter->detId());
     DetId detId(rawid);
     TrackerGeometry::ModuleType mType = tkGeom_->getDetectorType(detId);
     if (mType != TrackerGeometry::ModuleType::Ph2PSP && mType != TrackerGeometry::ModuleType::Ph2SS)
       continue;
     std::string key = phase2tkutil::getOTHistoId(detId, tTopo_);
     nTotvechitsinevt += DSViter->size();
     if (mType == TrackerGeometry::ModuleType::Ph2PSP) {
       if (nvechitLayerMapP.find(key) == nvechitLayerMapP.end()) {
         nvechitLayerMapP.insert(std::make_pair(key, DSViter->size()));
       } else {
         nvechitLayerMapP[key] += DSViter->size();
       }
     } else if (mType == TrackerGeometry::ModuleType::Ph2SS) {
       if (nvechitLayerMapS.find(key) == nvechitLayerMapS.end()) {
         nvechitLayerMapS.insert(std::make_pair(key, DSViter->size()));
       } else {
         nvechitLayerMapS[key] += DSViter->size();
       }
     }
 
     edmNew::DetSet<VectorHit>::const_iterator vechitIt;
     for (vechitIt = DSViter->begin(); vechitIt != DSViter->end(); ++vechitIt) {
       const Global3DPoint globalPos = vechitIt->lowerGlobalPos();
       const LocalPoint lp = vechitIt->localPosition();
       const double gx = globalPos.x() * 10.;
       const double gy = globalPos.y() * 10.;
       const double gz = globalPos.z() * 10.;
       const double gr = globalPos.perp() * 10.;
       const Global3DVector globalVec = vechitIt->globalDirection();
       const float curvature = vechitIt->curvature();
       const float curverr = vechitIt->curvatureError();
       const float eta = globalVec.eta();
       const float phi = globalVec.phi();
       float QOverPT = vechitIt->transverseMomentum(magField_->inTesla(GlobalPoint(0., 0., 0.)).z());
       const int sign = QOverPT > 0. ? 1. : -1.;
       if (mType == TrackerGeometry::ModuleType::Ph2PSP) {
         globalXY_P_->Fill(gx, gy);
         globalRZ_P_->Fill(gz, gr);
         //layer wise histo
         layerMEs_[key].localPosXY_P->Fill(lp.x(), lp.y());
         layerMEs_[key].curvature_P->Fill(curvature);
         if (curvature != 0.f)
           layerMEs_[key].curvErr_P->Fill(curverr / curvature);
         else
           edm::LogError("Phase2OTMonitorVectorHits") << "VectorHit with curvature zero found";
         layerMEs_[key].phi_P->Fill(phi);
         layerMEs_[key].eta_P->Fill(eta);
         layerMEs_[key].pt_P->Fill(sign * QOverPT);
         layerMEs_[key].chi2_P->Fill(vechitIt->chi2());
         layerMEs_[key].curvatureVsEta_P->Fill(eta, curvature);
       } else if (mType == TrackerGeometry::ModuleType::Ph2SS) {
         globalXY_S_->Fill(gx, gy);
         globalRZ_S_->Fill(gz, gr);
         //layer wise histo
         layerMEs_[key].localPosXY_S->Fill(lp.x(), lp.y());
         layerMEs_[key].curvature_S->Fill(curvature);
         if (curvature != 0.f)
           layerMEs_[key].curvErr_S->Fill(curverr / curvature);
         else
           edm::LogError("Phase2OTMonitorVectorHits") << "VectorHit with curvature zero found";
         layerMEs_[key].phi_S->Fill(phi);
         layerMEs_[key].eta_S->Fill(eta);
         layerMEs_[key].pt_S->Fill(sign * QOverPT);
         layerMEs_[key].chi2_S->Fill(vechitIt->chi2());
         layerMEs_[key].curvatureVsEta_S->Fill(eta, curvature);
       }
     }
   }
   //fill nVecHits per event
   numberVecHits_->Fill(nTotvechitsinevt);
   //fill nVecHit counter per layer
   for (auto& lme : nvechitLayerMapP) {
     layerMEs_[lme.first].numberVecHits_P->Fill(lme.second);
   }
   for (auto& lme : nvechitLayerMapS) {
     layerMEs_[lme.first].numberVecHits_S->Fill(lme.second);
   }
 }

◆ bookHistograms()

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 220 of file Phase2OTMonitorVectorHits.cc.

References phase2tkutil::book1DFromPSet(), phase2tkutil::book2DFromPSet(), bookLayerHistos(), dqm::implementation::NavigatorBase::cd(), edm::ESWatcher< T >::check(), config_, TrackerGeometry::detUnits(), TrackerGeometry::getDetectorType(), edm::ParameterSet::getParameter(), globalRZ_P_, globalRZ_S_, globalXY_P_, globalXY_S_, numberVecHits_, TrackerGeometry::Ph2PSP, TrackerGeometry::Ph2SS, dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, and tkGeom_.

                                                                             {
   std::string top_folder = config_.getParameter<std::string>("TopFolderName");
   //std::stringstream folder_name;
 
   ibooker.cd();
   edm::LogInfo("Phase2OTMonitorVectorHits") << " Booking Histograms in : " << top_folder;
   ibooker.setCurrentFolder(top_folder);
 
   //Global histos for OT
   numberVecHits_ = phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("GlobalNVecHits"), ibooker);
 
   globalXY_P_ = phase2tkutil::book2DFromPSet(config_.getParameter<edm::ParameterSet>("GlobalPositionXY_P"), ibooker);
 
   globalRZ_P_ = phase2tkutil::book2DFromPSet(config_.getParameter<edm::ParameterSet>("GlobalPositionRZ_P"), ibooker);
 
   globalXY_S_ = phase2tkutil::book2DFromPSet(config_.getParameter<edm::ParameterSet>("GlobalPositionXY_S"), ibooker);
 
   globalRZ_S_ = phase2tkutil::book2DFromPSet(config_.getParameter<edm::ParameterSet>("GlobalPositionRZ_S"), ibooker);
 
   //Now book layer wise histos
   edm::ESWatcher<TrackerDigiGeometryRecord> theTkDigiGeomWatcher;
   if (theTkDigiGeomWatcher.check(iSetup)) {
     for (auto const& det_u : tkGeom_->detUnits()) {
       unsigned int detId_raw = det_u->geographicalId().rawId();
       //we only need the layerwise histos for the lower layer (?)
       if (tkGeom_->getDetectorType(detId_raw) != TrackerGeometry::ModuleType::Ph2PSP &&
           tkGeom_->getDetectorType(detId_raw) != TrackerGeometry::ModuleType::Ph2SS)
         continue;
       bookLayerHistos(ibooker, detId_raw, top_folder);
     }
   }
 }

◆ bookLayerHistos()

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

private

Definition at line 258 of file Phase2OTMonitorVectorHits.cc.

Referenced by bookHistograms().

                                                                                                                 {
   std::string key = phase2tkutil::getOTHistoId(det_id, tTopo_);
   if (layerMEs_.find(key) == layerMEs_.end()) {
     ibooker.cd();
     VecHitME local_histos;
     ibooker.setCurrentFolder(subdir + "/" + key);
     edm::LogInfo("Phase2OTMonitorVectorHits") << " Booking Histograms in : " << key;
     //either PSP or SS - ensured by the call from bookHisto
     if (tkGeom_->getDetectorType(det_id) == TrackerGeometry::ModuleType::Ph2PSP) {
       local_histos.numberVecHits_P =
           phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("NVecHitsLayer_P"), ibooker);
       local_histos.localPosXY_P =
           phase2tkutil::book2DFromPSet(config_.getParameter<edm::ParameterSet>("LocalPositionXY_P"), ibooker);
       local_histos.curvature_P =
           phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("Curvature"), ibooker);
       local_histos.curvErr_P =
           phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("CurvErr"), ibooker);
       local_histos.phi_P = phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("Phi"), ibooker);
       local_histos.eta_P = phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("Eta"), ibooker);
       local_histos.pt_P = phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("Pt"), ibooker);
       local_histos.chi2_P = phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("Chi2"), ibooker);
       local_histos.curvatureVsEta_P =
           phase2tkutil::bookProfile1DFromPSet(config_.getParameter<edm::ParameterSet>("CurvatureVsEta_P"), ibooker);
     } else {
       ibooker.setCurrentFolder(subdir + "/" + key);
       local_histos.numberVecHits_S =
           phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("NVecHitsLayer_P"), ibooker);
       local_histos.localPosXY_S =
           phase2tkutil::book2DFromPSet(config_.getParameter<edm::ParameterSet>("LocalPositionXY_S"), ibooker);
       local_histos.curvature_S =
           phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("Curvature"), ibooker);
       local_histos.curvErr_S =
           phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("CurvErr"), ibooker);
       local_histos.phi_S = phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("Phi"), ibooker);
       local_histos.eta_S = phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("Eta"), ibooker);
       local_histos.pt_S = phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("Pt"), ibooker);
       local_histos.chi2_S = phase2tkutil::book1DFromPSet(config_.getParameter<edm::ParameterSet>("Chi2"), ibooker);
       local_histos.curvatureVsEta_S =
           phase2tkutil::bookProfile1DFromPSet(config_.getParameter<edm::ParameterSet>("CurvatureVsEta_S"), ibooker);
     }
     layerMEs_.insert(std::make_pair(key, local_histos));
   }
 }

◆ dqmBeginRun()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 117 of file Phase2OTMonitorVectorHits.cc.

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

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

◆ fillDescriptions()

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

static

Definition at line 304 of file Phase2OTMonitorVectorHits.cc.

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

                                                                                            {
   // vechitMonitorOT
   edm::ParameterSetDescription desc;
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("name", "NumberOfVecHits");
     psd0.add<std::string>("title", ";Number of vechits per event;");
     psd0.add<double>("xmin", 0.0);
     psd0.add<bool>("switch", true);
     psd0.add<double>("xmax", 15000.0);
     psd0.add<int>("NxBins", 150);
     desc.add<edm::ParameterSetDescription>("GlobalNVecHits", psd0);
   }
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("name", "Global_VecHitPosition_XY_P");
     psd0.add<std::string>("title", "Global_VecHitPosition_XY_P;x [mm];y [mm];");
     psd0.add<int>("NxBins", 1250);
     psd0.add<double>("xmin", -1250.0);
     psd0.add<double>("xmax", 1250.0);
     psd0.add<int>("NyBins", 1250);
     psd0.add<double>("ymin", -1250.0);
     psd0.add<double>("ymax", 1250.0);
     psd0.add<bool>("switch", true);
     desc.add<edm::ParameterSetDescription>("GlobalPositionXY_P", psd0);
   }
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("name", "Global_VecHitPosition_XY_S");
     psd0.add<std::string>("title", "Global_VecHitPosition_XY_S;x [mm];y [mm];");
     psd0.add<int>("NxBins", 1250);
     psd0.add<double>("xmin", -1250.0);
     psd0.add<double>("xmax", 1250.0);
     psd0.add<int>("NyBins", 1250);
     psd0.add<double>("ymin", -1250.0);
     psd0.add<double>("ymax", 1250.0);
     psd0.add<bool>("switch", true);
     desc.add<edm::ParameterSetDescription>("GlobalPositionXY_S", psd0);
   }
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("name", "Global_VecHitPosition_RZ_P");
     psd0.add<std::string>("title", "Global_VecHitPosition_RZ_P;z [mm];r [mm]");
     psd0.add<int>("NxBins", 1500);
     psd0.add<double>("xmin", -3000.0);
     psd0.add<double>("xmax", 3000.0);
     psd0.add<int>("NyBins", 1250);
     psd0.add<double>("ymin", 0.0);
     psd0.add<double>("ymax", 1250.0);
     psd0.add<bool>("switch", true);
     desc.add<edm::ParameterSetDescription>("GlobalPositionRZ_P", psd0);
   }
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("name", "Global_VecHitPosition_RZ_S");
     psd0.add<std::string>("title", "Global_VecHitPosition_RZ_S;z [mm];r [mm]");
 
     psd0.add<int>("NxBins", 1500);
     psd0.add<double>("xmin", -3000.0);
     psd0.add<double>("xmax", 3000.0);
     psd0.add<int>("NyBins", 1250);
     psd0.add<double>("ymin", 0.0);
     psd0.add<double>("ymax", 1250.0);
     psd0.add<bool>("switch", true);
     desc.add<edm::ParameterSetDescription>("GlobalPositionRZ_S", psd0);
   }
   //Layer wise parameter
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("name", "NumberOfVecHitsLayerP");
     psd0.add<std::string>("title", ";Number of vector hits per event(macro pixel sensor);");
     psd0.add<double>("xmin", 0.0);
     psd0.add<double>("xmax", 5000.0);
     psd0.add<int>("NxBins", 100);
     psd0.add<bool>("switch", true);
     desc.add<edm::ParameterSetDescription>("NVecHitsLayer_P", psd0);
   }
 
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("name", "NumberOfVecHitsLayerS");
     psd0.add<std::string>("title", ";Number of vector hits per event(strip sensor);");
     psd0.add<double>("xmin", 0.0);
     psd0.add<double>("xmax", 5000.0);
     psd0.add<int>("NxBins", 100);
     psd0.add<bool>("switch", true);
     desc.add<edm::ParameterSetDescription>("NVecHitsLayer_S", psd0);
   }
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("name", "LocalPositionXY_P");
     psd0.add<std::string>("title", "LocalPositionXY_P;x ;y ;");
     psd0.add<int>("NxBins", 50);
     psd0.add<double>("xmin", -10.0);
     psd0.add<double>("xmax", 10.0);
     psd0.add<int>("NyBins", 50);
     psd0.add<double>("ymin", -10.0);
     psd0.add<double>("ymax", 10.0);
     psd0.add<bool>("switch", true);
     desc.add<edm::ParameterSetDescription>("LocalPositionXY_P", psd0);
   }
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("name", "LocalPositionXY_S");
     psd0.add<std::string>("title", "LocalPositionXY_S;x ;y ;");
     psd0.add<int>("NxBins", 50);
     psd0.add<double>("xmin", -10.0);
     psd0.add<double>("xmax", 10.0);
     psd0.add<int>("NyBins", 50);
     psd0.add<double>("ymin", -10.0);
     psd0.add<double>("ymax", 10.0);
     psd0.add<bool>("switch", true);
     desc.add<edm::ParameterSetDescription>("LocalPositionXY_S", psd0);
   }
 
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("name", "CurvatureOfVecHits");
     psd0.add<std::string>("title", ";VectorHit curvature;");
     psd0.add<double>("xmin", -0.05);
     psd0.add<bool>("switch", true);
     psd0.add<double>("xmax", 0.05);
     psd0.add<int>("NxBins", 200);
     desc.add<edm::ParameterSetDescription>("Curvature", psd0);
   }
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("name", "CurvatureErrOCurvature");
     psd0.add<std::string>("title", ";VectorHit #delta#rho/#rho;");
     psd0.add<double>("xmin", -0.05);
     psd0.add<bool>("switch", true);
     psd0.add<double>("xmax", 0.05);
     psd0.add<int>("NxBins", 500);
     desc.add<edm::ParameterSetDescription>("CurvErr", psd0);
   }
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("name", "PhiOfVecHits");
     psd0.add<std::string>("title", ";VectorHit #phi;");
     psd0.add<double>("xmin", -M_PI);
     psd0.add<bool>("switch", true);
     psd0.add<double>("xmax", M_PI);
     psd0.add<int>("NxBins", 30);
     desc.add<edm::ParameterSetDescription>("Phi", psd0);
   }
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("name", "EtaOfVecHits");
     psd0.add<std::string>("title", ";VectorHit #eta;");
     psd0.add<double>("xmin", -5.);
     psd0.add<bool>("switch", true);
     psd0.add<double>("xmax", 5.);
     psd0.add<int>("NxBins", 50);
     desc.add<edm::ParameterSetDescription>("Eta", psd0);
   }
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("name", "PtOfVecHits");
     psd0.add<std::string>("title", "VectorHit p_T;p_T ;");
     psd0.add<int>("NxBins", 100);
     psd0.add<double>("xmin", 0.);
     psd0.add<double>("xmax", 200.0);
     psd0.add<bool>("switch", true);
     desc.add<edm::ParameterSetDescription>("Pt", psd0);
   }
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("name", "Chi2OfVecHits");
     psd0.add<std::string>("title", "VectorHit chi squared; #chi^2;");
     psd0.add<int>("NxBins", 100);
     psd0.add<double>("xmin", 0.);
     psd0.add<double>("xmax", 0.000001);
     psd0.add<bool>("switch", true);
     desc.add<edm::ParameterSetDescription>("Chi2", psd0);
   }
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("name", "CurvatureVsEtaProf_P");
     psd0.add<std::string>("title", "Curvature vs #eta (macro-pixel);#eta ;curvature ;");
     psd0.add<int>("NxBins", 50);
     psd0.add<double>("xmin", -5.0);
     psd0.add<double>("xmax", 5.0);
     psd0.add<double>("ymin", -0.05);
     psd0.add<double>("ymax", 0.05);
     psd0.add<bool>("switch", true);
     desc.add<edm::ParameterSetDescription>("CurvatureVsEta_P", psd0);
   }
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("name", "CurvatureVsEtaProf_S");
     psd0.add<std::string>("title", "Curvature vs #eta (strip);#eta ;curvature ;");
     psd0.add<int>("NxBins", 25);
     psd0.add<double>("xmin", -5.0);
     psd0.add<double>("xmax", 5.0);
     psd0.add<double>("ymin", -0.05);
     psd0.add<double>("ymax", 0.05);
     psd0.add<bool>("switch", true);
     desc.add<edm::ParameterSetDescription>("CurvatureVsEta_S", psd0);
   }
   desc.add<std::string>("TopFolderName", "TrackerPhase2OTVectorHits/Accepted");
   desc.add<bool>("Verbosity", false);
   desc.add<edm::InputTag>("vechitsSrc", edm::InputTag("siPhase2VectorHits", "accepted"));
   descriptions.add("Phase2OTMonitorVectorHits", desc);
 }