#include <OuterTrackerMonitorTrackingParticles.h>

Inheritance diagram for OuterTrackerMonitorTrackingParticles:

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

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

	OuterTrackerMonitorTrackingParticles (const edm::ParameterSet &)

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

virtual void	dqmBeginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

virtual void	dqmEndLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

virtual void	dqmEndRun (edm::Run const &, edm::EventSetup const &)

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &, edm::EventSetup const &) final

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, edm::Accumulator >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () 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

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

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)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Public Attributes
MonitorElement *	match_tp_d0 = 0

MonitorElement *	match_tp_eta = 0

MonitorElement *	match_tp_pt = 0

MonitorElement *	match_tp_pt_zoom = 0

MonitorElement *	match_tp_VtxR = 0

MonitorElement *	match_tp_VtxZ = 0

MonitorElement *	res_eta = 0

MonitorElement *	res_pt = 0

MonitorElement *	res_ptRel = 0

MonitorElement *	resd0_eta0p7to1 = 0

MonitorElement *	resd0_eta0to0p7 = 0

MonitorElement *	resd0_eta1p2to1p6 = 0

MonitorElement *	resd0_eta1p6to2 = 0

MonitorElement *	resd0_eta1to1p2 = 0

MonitorElement *	resd0_eta2to2p4 = 0

MonitorElement *	reseta_eta0p7to1 = 0

MonitorElement *	reseta_eta0to0p7 = 0

MonitorElement *	reseta_eta1p2to1p6 = 0

MonitorElement *	reseta_eta1p6to2 = 0

MonitorElement *	reseta_eta1to1p2 = 0

MonitorElement *	reseta_eta2to2p4 = 0

MonitorElement *	resphi_eta0p7to1 = 0

MonitorElement *	resphi_eta0to0p7 = 0

MonitorElement *	resphi_eta1p2to1p6 = 0

MonitorElement *	resphi_eta1p6to2 = 0

MonitorElement *	resphi_eta1to1p2 = 0

MonitorElement *	resphi_eta2to2p4 = 0

MonitorElement *	respt_eta0p7to1_pt2to3 = 0

MonitorElement *	respt_eta0p7to1_pt3to8 = 0

MonitorElement *	respt_eta0p7to1_pt8toInf = 0

MonitorElement *	respt_eta0to0p7_pt2to3 = 0

MonitorElement *	respt_eta0to0p7_pt3to8 = 0

MonitorElement *	respt_eta0to0p7_pt8toInf = 0

MonitorElement *	respt_eta1p2to1p6_pt2to3 = 0

MonitorElement *	respt_eta1p2to1p6_pt3to8 = 0

MonitorElement *	respt_eta1p2to1p6_pt8toInf = 0

MonitorElement *	respt_eta1p6to2_pt2to3 = 0

MonitorElement *	respt_eta1p6to2_pt3to8 = 0

MonitorElement *	respt_eta1p6to2_pt8toInf = 0

MonitorElement *	respt_eta1to1p2_pt2to3 = 0

MonitorElement *	respt_eta1to1p2_pt3to8 = 0

MonitorElement *	respt_eta1to1p2_pt8toInf = 0

MonitorElement *	respt_eta2to2p4_pt2to3 = 0

MonitorElement *	respt_eta2to2p4_pt3to8 = 0

MonitorElement *	respt_eta2to2p4_pt8toInf = 0

MonitorElement *	resVtxZ_eta0p7to1 = 0

MonitorElement *	resVtxZ_eta0to0p7 = 0

MonitorElement *	resVtxZ_eta1p2to1p6 = 0

MonitorElement *	resVtxZ_eta1p6to2 = 0

MonitorElement *	resVtxZ_eta1to1p2 = 0

MonitorElement *	resVtxZ_eta2to2p4 = 0

MonitorElement *	tp_d0 = 0

MonitorElement *	tp_eta = 0

MonitorElement *	tp_pt = 0

MonitorElement *	tp_pt_zoom = 0

MonitorElement *	tp_VtxR = 0

MonitorElement *	tp_VtxZ = 0

MonitorElement *	Track_MatchedChi2 = 0

MonitorElement *	Track_MatchedChi2Red = 0

MonitorElement *	trackParts_Eta = 0

MonitorElement *	trackParts_Phi = 0

MonitorElement *	trackParts_Pt = 0

Private Attributes
edm::ParameterSet	conf_

double	L1Tk_maxChi2

double	L1Tk_maxChi2dof

int	L1Tk_minNStub

int	L1Tk_nPar

std::string	topFolderName_

double	TP_maxEta

double	TP_maxPt

double	TP_maxVtxZ

int	TP_minNLayersStub

int	TP_minNStub

double	TP_minPt

int	TP_select_eventid

edm::EDGetTokenT< std::vector< TrackingParticle > >	trackingParticleToken_

edm::EDGetTokenT< TTClusterAssociationMap< Ref_Phase2TrackerDigi_ > >	ttClusterMCTruthToken_

edm::EDGetTokenT< TTStubAssociationMap< Ref_Phase2TrackerDigi_ > >	ttStubMCTruthToken_

edm::EDGetTokenT< TTTrackAssociationMap< Ref_Phase2TrackerDigi_ > >	ttTrackMCTruthToken_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

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::ProducerBase
ProducesCollector	producesCollector ()

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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

Definition at line 24 of file OuterTrackerMonitorTrackingParticles.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 39 of file OuterTrackerMonitorTrackingParticles.cc.

References conf_, edm::ParameterSet::getParameter(), L1Tk_maxChi2, L1Tk_maxChi2dof, L1Tk_minNStub, L1Tk_nPar, AlCaHLTBitMon_QueryRunRegistry::string, topFolderName_, TP_maxEta, TP_maxPt, TP_maxVtxZ, TP_minNLayersStub, TP_minNStub, TP_minPt, TP_select_eventid, trackingParticleToken_, ttClusterMCTruthToken_, ttStubMCTruthToken_, and ttTrackMCTruthToken_.

     : conf_(iConfig) {
   topFolderName_ = conf_.getParameter<std::string>("TopFolderName");
   trackingParticleToken_ =
       consumes<std::vector<TrackingParticle>>(conf_.getParameter<edm::InputTag>("trackingParticleToken"));
   ttStubMCTruthToken_ =
       consumes<TTStubAssociationMap<Ref_Phase2TrackerDigi_>>(conf_.getParameter<edm::InputTag>("MCTruthStubInputTag"));
   ttClusterMCTruthToken_ = consumes<TTClusterAssociationMap<Ref_Phase2TrackerDigi_>>(
       conf_.getParameter<edm::InputTag>("MCTruthClusterInputTag"));
   ttTrackMCTruthToken_ = consumes<TTTrackAssociationMap<Ref_Phase2TrackerDigi_>>(
       conf_.getParameter<edm::InputTag>("MCTruthTrackInputTag"));
   L1Tk_nPar = conf_.getParameter<int>("L1Tk_nPar");                 // 4 or 5(d0) track parameters
   L1Tk_minNStub = conf_.getParameter<int>("L1Tk_minNStub");         // min number of stubs in the track
   L1Tk_maxChi2 = conf_.getParameter<double>("L1Tk_maxChi2");        // maximum chi2 of the track
   L1Tk_maxChi2dof = conf_.getParameter<double>("L1Tk_maxChi2dof");  // maximum chi2/dof of the track
   TP_minNStub = conf_.getParameter<int>("TP_minNStub");             // min number of stubs in the tracking particle to
   //min number of layers with stubs in the tracking particle to consider matching
   TP_minNLayersStub = conf_.getParameter<int>("TP_minNLayersStub");
   TP_minPt = conf_.getParameter<double>("TP_minPt");                 // min pT to consider matching
   TP_maxPt = conf_.getParameter<double>("TP_maxPt");                 // max pT to consider matching
   TP_maxEta = conf_.getParameter<double>("TP_maxEta");               // max eta to consider matching
   TP_maxVtxZ = conf_.getParameter<double>("TP_maxVtxZ");             // max vertZ (or z0) to consider matching
   TP_select_eventid = conf_.getParameter<int>("TP_select_eventid");  // PI or not
 }

OuterTrackerMonitorTrackingParticles::~OuterTrackerMonitorTrackingParticles ( )

override

Definition at line 64 of file OuterTrackerMonitorTrackingParticles.cc.

                                                                             {
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
 }

Member Function Documentation

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

override

Definition at line 72 of file OuterTrackerMonitorTrackingParticles.cc.

References MaterialEffects_cfi::A, TrackingParticle::charge(), funct::cos(), JetChargeProducer_cfi::exp, dqm::impl::MonitorElement::Fill(), edm::EventSetup::get(), edm::Event::getByToken(), createfilelist::int, isum, edm::HandleBase::isValid(), L1Tk_maxChi2, L1Tk_maxChi2dof, L1Tk_minNStub, L1Tk_nPar, TrackerTopology::layer(), match_tp_d0, match_tp_eta, match_tp_pt, match_tp_pt_zoom, match_tp_VtxR, match_tp_VtxZ, TrackingParticle::p4(), TrackingParticle::pdgId(), pi, edm::ESHandle< T >::product(), res_eta, res_pt, res_ptRel, resd0_eta0p7to1, resd0_eta0to0p7, resd0_eta1p2to1p6, resd0_eta1p6to2, resd0_eta1to1p2, resd0_eta2to2p4, reseta_eta0p7to1, reseta_eta0to0p7, reseta_eta1p2to1p6, reseta_eta1p6to2, reseta_eta1to1p2, reseta_eta2to2p4, resphi_eta0p7to1, resphi_eta0to0p7, resphi_eta1p2to1p6, resphi_eta1p6to2, resphi_eta1to1p2, resphi_eta2to2p4, respt_eta0p7to1_pt2to3, respt_eta0p7to1_pt3to8, respt_eta0p7to1_pt8toInf, respt_eta0to0p7_pt2to3, respt_eta0to0p7_pt3to8, respt_eta0to0p7_pt8toInf, respt_eta1p2to1p6_pt2to3, respt_eta1p2to1p6_pt3to8, respt_eta1p2to1p6_pt8toInf, respt_eta1p6to2_pt2to3, respt_eta1p6to2_pt3to8, respt_eta1p6to2_pt8toInf, respt_eta1to1p2_pt2to3, respt_eta1to1p2_pt3to8, respt_eta1to1p2_pt8toInf, respt_eta2to2p4_pt2to3, respt_eta2to2p4_pt3to8, respt_eta2to2p4_pt8toInf, resVtxZ_eta0p7to1, resVtxZ_eta0to0p7, resVtxZ_eta1p2to1p6, resVtxZ_eta1p6to2, resVtxZ_eta1to1p2, resVtxZ_eta2to2p4, funct::sin(), mathSSE::sqrt(), funct::tan(), StripSubdetector::TID, StripSubdetector::TOB, tp_d0, tp_eta, TP_maxEta, TP_maxPt, TP_maxVtxZ, TP_minNLayersStub, TP_minNStub, TP_minPt, tp_pt, tp_pt_zoom, TP_select_eventid, tp_VtxR, tp_VtxZ, Track_MatchedChi2, Track_MatchedChi2Red, trackingParticleToken_, trackParts_Eta, trackParts_Phi, trackParts_Pt, ttClusterMCTruthToken_, ttStubMCTruthToken_, and ttTrackMCTruthToken_.

                                                                                                       {
   // Tracking Particles
   edm::Handle<std::vector<TrackingParticle>> trackingParticleHandle;
   iEvent.getByToken(trackingParticleToken_, trackingParticleHandle);
 
   // Truth Association Maps
   edm::Handle<TTTrackAssociationMap<Ref_Phase2TrackerDigi_>> MCTruthTTTrackHandle;
   iEvent.getByToken(ttTrackMCTruthToken_, MCTruthTTTrackHandle);
   edm::Handle<TTClusterAssociationMap<Ref_Phase2TrackerDigi_>> MCTruthTTClusterHandle;
   iEvent.getByToken(ttClusterMCTruthToken_, MCTruthTTClusterHandle);
   edm::Handle<TTStubAssociationMap<Ref_Phase2TrackerDigi_>> MCTruthTTStubHandle;
   iEvent.getByToken(ttStubMCTruthToken_, MCTruthTTStubHandle);
 
   // Geometries
   edm::ESHandle<TrackerTopology> tTopoHandle;
   iSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
   const TrackerTopology *const tTopo = tTopoHandle.product();
 
   // Loop over tracking particles
   int this_tp = 0;
   for (auto iterTP : *trackingParticleHandle) {
     edm::Ptr<TrackingParticle> tp_ptr(trackingParticleHandle, this_tp);
     this_tp++;
 
     int tmp_eventid = iterTP.eventId().event();
     float tmp_tp_pt = iterTP.pt();
     float tmp_tp_phi = iterTP.phi();
     float tmp_tp_eta = iterTP.eta();
 
     //Calculate nLayers variable
     std::vector<edm::Ref<edmNew::DetSetVector<TTStub<Ref_Phase2TrackerDigi_>>, TTStub<Ref_Phase2TrackerDigi_>>>
         theStubRefs = MCTruthTTStubHandle->findTTStubRefs(tp_ptr);
 
     int hasStubInLayer[11] = {0};
     for (unsigned int is = 0; is < theStubRefs.size(); is++) {
       DetId detid(theStubRefs.at(is)->getDetId());
       int layer = -1;
       if (detid.subdetId() == StripSubdetector::TOB)
         layer = static_cast<int>(tTopo->layer(detid)) - 1;  //fill in array as entries 0-5
       else if (detid.subdetId() == StripSubdetector::TID)
         layer = static_cast<int>(tTopo->layer(detid)) + 5;  //fill in array as entries 6-10
 
       //treat genuine stubs separately (==2 is genuine, ==1 is not)
       if (MCTruthTTStubHandle->findTrackingParticlePtr(theStubRefs.at(is)).isNull() && hasStubInLayer[layer] < 2)
         hasStubInLayer[layer] = 1;
       else
         hasStubInLayer[layer] = 2;
     }
 
     int nStubLayerTP = 0;
     int nStubLayerTP_g = 0;
     for (int isum = 0; isum < 11; isum++) {
       if (hasStubInLayer[isum] >= 1)
         nStubLayerTP += 1;
       else if (hasStubInLayer[isum] == 2)
         nStubLayerTP_g += 1;
     }
 
     if (std::fabs(tmp_tp_eta) > TP_maxEta)
       continue;
     // Fill the 1D distribution plots for tracking particles, to monitor change in stub definition
     if (tmp_tp_pt > TP_minPt && nStubLayerTP >= TP_minNLayersStub) {
       trackParts_Pt->Fill(tmp_tp_pt);
       trackParts_Eta->Fill(tmp_tp_eta);
       trackParts_Phi->Fill(tmp_tp_phi);
     }
 
     if (TP_select_eventid == 0 && tmp_eventid != 0)
       continue;  //only care about tracking particles from the primary interaction for efficiency/resolution
     int nStubTP = -1;
     if (MCTruthTTStubHandle.isValid()) {
       std::vector<edm::Ref<edmNew::DetSetVector<TTStub<Ref_Phase2TrackerDigi_>>, TTStub<Ref_Phase2TrackerDigi_>>>
           theStubRefs = MCTruthTTStubHandle->findTTStubRefs(tp_ptr);
       nStubTP = (int)theStubRefs.size();
     }
     if (MCTruthTTClusterHandle.isValid() && MCTruthTTClusterHandle->findTTClusterRefs(tp_ptr).empty())
       continue;
 
     float tmp_tp_vz = iterTP.vz();
     float tmp_tp_vx = iterTP.vx();
     float tmp_tp_vy = iterTP.vy();
     float tmp_tp_charge = tp_ptr->charge();
     int tmp_tp_pdgid = iterTP.pdgId();
 
     // ----------------------------------------------------------------------------------------------
     // calculate d0 and VtxZ propagated back to the IP, pass if greater than max
     // VtxZ
     float tmp_tp_t = tan(2.0 * atan(1.0) - 2.0 * atan(exp(-tmp_tp_eta)));
     float delx = -tmp_tp_vx;
     float dely = -tmp_tp_vy;
     float K = 0.01 * 0.5696 / tmp_tp_pt * tmp_tp_charge;  // curvature correction
     float A = 1. / (2. * K);
     float tmp_tp_x0p = delx - A * sin(tmp_tp_phi);
     float tmp_tp_y0p = dely + A * cos(tmp_tp_phi);
     float tmp_tp_rp = sqrt(tmp_tp_x0p * tmp_tp_x0p + tmp_tp_y0p * tmp_tp_y0p);
     static double pi = 4.0 * atan(1.0);
     float delphi = tmp_tp_phi - atan2(-K * tmp_tp_x0p, K * tmp_tp_y0p);
     if (delphi < -pi)
       delphi += 2.0 * pi;
     if (delphi > pi)
       delphi -= 2.0 * pi;
 
     float tmp_tp_VtxZ = tmp_tp_vz + tmp_tp_t * delphi / (2.0 * K);
     float tmp_tp_VtxR = sqrt(tmp_tp_vx * tmp_tp_vx + tmp_tp_vy * tmp_tp_vy);
     float tmp_tp_d0 = tmp_tp_charge * tmp_tp_rp - (1. / (2. * K));
 
     // simpler formula for d0, in cases where the charge is zero:
     // https://github.com/cms-sw/cmssw/blob/master/DataFormats/TrackReco/interface/TrackBase.h
     float other_d0 = -tmp_tp_vx * sin(tmp_tp_phi) + tmp_tp_vy * cos(tmp_tp_phi);
     tmp_tp_d0 = tmp_tp_d0 * (-1);  // fix d0 sign
     if (K == 0) {
       tmp_tp_d0 = other_d0;
       tmp_tp_VtxZ = tmp_tp_vz;
     }
     if (std::fabs(tmp_tp_VtxZ) > TP_maxVtxZ)
       continue;
 
     // To make efficiency plots where the denominator has NO stub cuts
     if (tmp_tp_VtxR < 1.0) {
       tp_pt->Fill(tmp_tp_pt);  //pT effic, no cut on pT, but VtxR cut
       if (tmp_tp_pt <= 10)
         tp_pt_zoom->Fill(tmp_tp_pt);  //pT effic, no cut on pT, but VtxR cut
     }
     if (tmp_tp_pt < TP_minPt)
       continue;
     tp_VtxR->Fill(tmp_tp_VtxR);  // VtxR efficiency has no cut on VtxR
     if (tmp_tp_VtxR > 1.0)
       continue;
     tp_eta->Fill(tmp_tp_eta);
     tp_d0->Fill(tmp_tp_d0);
     tp_VtxZ->Fill(tmp_tp_VtxZ);
 
     if (nStubTP < TP_minNStub || nStubLayerTP < TP_minNLayersStub)
       continue;  //nStub cut not included in denominator of efficiency plots
 
     // ----------------------------------------------------------------------------------------------
     // look for L1 tracks matched to the tracking particle
     int tp_nMatch = 0;
     int i_track = -1;
     float i_chi2dof = 99999;
     if (MCTruthTTTrackHandle.isValid()) {
       std::vector<edm::Ptr<TTTrack<Ref_Phase2TrackerDigi_>>> matchedTracks =
           MCTruthTTTrackHandle->findTTTrackPtrs(tp_ptr);
 
       // ----------------------------------------------------------------------------------------------
       // loop over matched L1 tracks
       // here, "match" means tracks that can be associated to a TrackingParticle
       // with at least one hit of at least one of its clusters
       // https://twiki.cern.ch/twiki/bin/viewauth/CMS/SLHCTrackerTriggerSWTools#MC_truth_for_TTTrack
       int trkCounter = 0;
       for (auto thisTrack : matchedTracks) {
         if (!MCTruthTTTrackHandle->isGenuine(thisTrack))
           continue;
         // ----------------------------------------------------------------------------------------------
         // further require L1 track to be (loosely) genuine, that there is only
         // one TP matched to the track
         // + have >= L1Tk_minNStub stubs for it to be a valid match
         int tmp_trk_nstub = thisTrack->getStubRefs().size();
         if (tmp_trk_nstub < L1Tk_minNStub)
           continue;
         float dmatch_pt = 999;
         float dmatch_eta = 999;
         float dmatch_phi = 999;
         int match_id = 999;
 
         edm::Ptr<TrackingParticle> my_tp = MCTruthTTTrackHandle->findTrackingParticlePtr(thisTrack);
         dmatch_pt = std::fabs(my_tp->p4().pt() - tmp_tp_pt);
         dmatch_eta = std::fabs(my_tp->p4().eta() - tmp_tp_eta);
         dmatch_phi = std::fabs(my_tp->p4().phi() - tmp_tp_phi);
         match_id = my_tp->pdgId();
         float tmp_trk_chi2dof = (thisTrack->getChi2(L1Tk_nPar)) / (2 * tmp_trk_nstub - L1Tk_nPar);
 
         // ensure that track is uniquely matched to the TP we are looking at!
         if (dmatch_pt < 0.1 && dmatch_eta < 0.1 && dmatch_phi < 0.1 && tmp_tp_pdgid == match_id) {
           tp_nMatch++;
           if (i_track < 0 || tmp_trk_chi2dof < i_chi2dof) {
             i_track = trkCounter;
             i_chi2dof = tmp_trk_chi2dof;
           }
         }
         trkCounter++;
       }  // end loop over matched L1 tracks
 
       if (tp_nMatch < 1)
         continue;
       // Get information on the matched tracks
       float tmp_matchtrk_pt = -999;
       float tmp_matchtrk_eta = -999;
       float tmp_matchtrk_phi = -999;
       float tmp_matchtrk_VtxZ = -999;
       float tmp_matchtrk_chi2 = -999;
       float tmp_matchtrk_chi2dof = -999;
       int tmp_matchTrk_nStub = -999;
       float tmp_matchtrk_d0 = -999;
 
       tmp_matchtrk_pt = matchedTracks[i_track]->getMomentum(L1Tk_nPar).perp();
       tmp_matchtrk_eta = matchedTracks[i_track]->getMomentum(L1Tk_nPar).eta();
       tmp_matchtrk_phi = matchedTracks[i_track]->getMomentum(L1Tk_nPar).phi();
       tmp_matchtrk_VtxZ = matchedTracks[i_track]->getPOCA(L1Tk_nPar).z();
       tmp_matchtrk_chi2 = matchedTracks[i_track]->getChi2(L1Tk_nPar);
       tmp_matchtrk_chi2dof = matchedTracks[i_track]->getChi2Red(L1Tk_nPar);
       tmp_matchTrk_nStub = (int)matchedTracks[i_track]->getStubRefs().size();
 
       Track_MatchedChi2->Fill(tmp_matchtrk_chi2);
       Track_MatchedChi2Red->Fill(tmp_matchtrk_chi2dof);
 
       //for d0
       float tmp_matchtrk_x0 = matchedTracks[i_track]->getPOCA(L1Tk_nPar).x();
       float tmp_matchtrk_y0 = matchedTracks[i_track]->getPOCA(L1Tk_nPar).y();
       tmp_matchtrk_d0 = -tmp_matchtrk_x0 * sin(tmp_matchtrk_phi) + tmp_matchtrk_y0 * cos(tmp_matchtrk_phi);
 
       //Add cuts for the matched tracks, numerator
       if (tmp_tp_pt > TP_maxPt)
         continue;
       if (tmp_matchTrk_nStub < L1Tk_minNStub || tmp_matchtrk_chi2 > L1Tk_maxChi2 ||
           tmp_matchtrk_chi2dof > L1Tk_maxChi2dof)
         continue;
 
       // fill matched track histograms (if passes all criteria)
       match_tp_pt->Fill(tmp_tp_pt);
       if (tmp_tp_pt > 0 && tmp_tp_pt <= 10)
         match_tp_pt_zoom->Fill(tmp_tp_pt);
       match_tp_eta->Fill(tmp_tp_eta);
       match_tp_d0->Fill(tmp_tp_d0);
       match_tp_VtxR->Fill(tmp_tp_VtxR);
       match_tp_VtxZ->Fill(tmp_tp_VtxZ);
 
       // Eta and pT histograms for resolution
       float pt_diff = tmp_matchtrk_pt - tmp_tp_pt;
       float pt_res = pt_diff / tmp_tp_pt;
       float eta_res = tmp_matchtrk_eta - tmp_tp_eta;
       float phi_res = tmp_matchtrk_phi - tmp_tp_phi;
       float VtxZ_res = tmp_matchtrk_VtxZ - tmp_tp_VtxZ;
       float d0_res = tmp_matchtrk_d0 - tmp_tp_d0;
 
       // fill total resolution histograms
       res_pt->Fill(pt_diff);
       res_ptRel->Fill(pt_res);
       res_eta->Fill(eta_res);
 
       // Fill resolution plots for different abs(eta) bins:
       // (0, 0.7), (0.7, 1.0), (1.0, 1.2), (1.2, 1.6), (1.6, 2.0), (2.0, 2.4)
       if (std::fabs(tmp_tp_eta) >= 0 && std::fabs(tmp_tp_eta) < 0.7) {
         reseta_eta0to0p7->Fill(eta_res);
         resphi_eta0to0p7->Fill(phi_res);
         resVtxZ_eta0to0p7->Fill(VtxZ_res);
         resd0_eta0to0p7->Fill(d0_res);
         if (tmp_tp_pt >= 2 && tmp_tp_pt < 3)
           respt_eta0to0p7_pt2to3->Fill(pt_res);
         else if (tmp_tp_pt >= 3 && tmp_tp_pt < 8)
           respt_eta0to0p7_pt3to8->Fill(pt_res);
         else if (tmp_tp_pt >= 8)
           respt_eta0to0p7_pt8toInf->Fill(pt_res);
       } else if (std::fabs(tmp_tp_eta) >= 0.7 && std::fabs(tmp_tp_eta) < 1.0) {
         reseta_eta0p7to1->Fill(eta_res);
         resphi_eta0p7to1->Fill(phi_res);
         resVtxZ_eta0p7to1->Fill(VtxZ_res);
         resd0_eta0p7to1->Fill(d0_res);
         if (tmp_tp_pt >= 2 && tmp_tp_pt < 3)
           respt_eta0p7to1_pt2to3->Fill(pt_res);
         else if (tmp_tp_pt >= 3 && tmp_tp_pt < 8)
           respt_eta0p7to1_pt3to8->Fill(pt_res);
         else if (tmp_tp_pt >= 8)
           respt_eta0p7to1_pt8toInf->Fill(pt_res);
       } else if (std::fabs(tmp_tp_eta) >= 1.0 && std::fabs(tmp_tp_eta) < 1.2) {
         reseta_eta1to1p2->Fill(eta_res);
         resphi_eta1to1p2->Fill(phi_res);
         resVtxZ_eta1to1p2->Fill(VtxZ_res);
         resd0_eta1to1p2->Fill(d0_res);
         if (tmp_tp_pt >= 2 && tmp_tp_pt < 3)
           respt_eta1to1p2_pt2to3->Fill(pt_res);
         else if (tmp_tp_pt >= 3 && tmp_tp_pt < 8)
           respt_eta1to1p2_pt3to8->Fill(pt_res);
         else if (tmp_tp_pt >= 8)
           respt_eta1to1p2_pt8toInf->Fill(pt_res);
       } else if (std::fabs(tmp_tp_eta) >= 1.2 && std::fabs(tmp_tp_eta) < 1.6) {
         reseta_eta1p2to1p6->Fill(eta_res);
         resphi_eta1p2to1p6->Fill(phi_res);
         resVtxZ_eta1p2to1p6->Fill(VtxZ_res);
         resd0_eta1p2to1p6->Fill(d0_res);
         if (tmp_tp_pt >= 2 && tmp_tp_pt < 3)
           respt_eta1p2to1p6_pt2to3->Fill(pt_res);
         else if (tmp_tp_pt >= 3 && tmp_tp_pt < 8)
           respt_eta1p2to1p6_pt3to8->Fill(pt_res);
         else if (tmp_tp_pt >= 8)
           respt_eta1p2to1p6_pt8toInf->Fill(pt_res);
       } else if (std::fabs(tmp_tp_eta) >= 1.6 && std::fabs(tmp_tp_eta) < 2.0) {
         reseta_eta1p6to2->Fill(eta_res);
         resphi_eta1p6to2->Fill(phi_res);
         resVtxZ_eta1p6to2->Fill(VtxZ_res);
         resd0_eta1p6to2->Fill(d0_res);
         if (tmp_tp_pt >= 2 && tmp_tp_pt < 3)
           respt_eta1p6to2_pt2to3->Fill(pt_res);
         else if (tmp_tp_pt >= 3 && tmp_tp_pt < 8)
           respt_eta1p6to2_pt3to8->Fill(pt_res);
         else if (tmp_tp_pt >= 8)
           respt_eta1p6to2_pt8toInf->Fill(pt_res);
       } else if (std::fabs(tmp_tp_eta) >= 2.0 && std::fabs(tmp_tp_eta) <= 2.4) {
         reseta_eta2to2p4->Fill(eta_res);
         resphi_eta2to2p4->Fill(phi_res);
         resVtxZ_eta2to2p4->Fill(VtxZ_res);
         resd0_eta2to2p4->Fill(d0_res);
         if (tmp_tp_pt >= 2 && tmp_tp_pt < 3)
           respt_eta2to2p4_pt2to3->Fill(pt_res);
         else if (tmp_tp_pt >= 3 && tmp_tp_pt < 8)
           respt_eta2to2p4_pt3to8->Fill(pt_res);
         else if (tmp_tp_pt >= 8)
           respt_eta2to2p4_pt8toInf->Fill(pt_res);
       }
     }  //if MC TTTrack handle is valid
   }    //end loop over tracking particles
 }  // end of method

void OuterTrackerMonitorTrackingParticles::bookHistograms	(	DQMStore::IBooker &	iBooker,
		edm::Run const &	run,
		edm::EventSetup const &	es
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 387 of file OuterTrackerMonitorTrackingParticles.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book1D(), conf_, DEFINE_FWK_MODULE, edm::ParameterSet::getParameter(), match_tp_d0, match_tp_eta, match_tp_pt, match_tp_pt_zoom, match_tp_VtxR, match_tp_VtxZ, res_eta, res_pt, res_ptRel, resd0_eta0p7to1, resd0_eta0to0p7, resd0_eta1p2to1p6, resd0_eta1p6to2, resd0_eta1to1p2, resd0_eta2to2p4, reseta_eta0p7to1, reseta_eta0to0p7, reseta_eta1p2to1p6, reseta_eta1p6to2, reseta_eta1to1p2, reseta_eta2to2p4, resphi_eta0p7to1, resphi_eta0to0p7, resphi_eta1p2to1p6, resphi_eta1p6to2, resphi_eta1to1p2, resphi_eta2to2p4, respt_eta0p7to1_pt2to3, respt_eta0p7to1_pt3to8, respt_eta0p7to1_pt8toInf, respt_eta0to0p7_pt2to3, respt_eta0to0p7_pt3to8, respt_eta0to0p7_pt8toInf, respt_eta1p2to1p6_pt2to3, respt_eta1p2to1p6_pt3to8, respt_eta1p2to1p6_pt8toInf, respt_eta1p6to2_pt2to3, respt_eta1p6to2_pt3to8, respt_eta1p6to2_pt8toInf, respt_eta1to1p2_pt2to3, respt_eta1to1p2_pt3to8, respt_eta1to1p2_pt8toInf, respt_eta2to2p4_pt2to3, respt_eta2to2p4_pt3to8, respt_eta2to2p4_pt8toInf, resVtxZ_eta0p7to1, resVtxZ_eta0to0p7, resVtxZ_eta1p2to1p6, resVtxZ_eta1p6to2, resVtxZ_eta1to1p2, resVtxZ_eta2to2p4, dqm::impl::MonitorElement::setAxisTitle(), dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, topFolderName_, tp_d0, tp_eta, tp_pt, tp_pt_zoom, tp_VtxR, tp_VtxZ, Track_MatchedChi2, Track_MatchedChi2Red, trackParts_Eta, trackParts_Phi, and trackParts_Pt.

                                                                                    {
   // Histogram setup and definitions
   std::string HistoName;
   iBooker.setCurrentFolder(topFolderName_ + "/trackParticles");
 
   // 1D: pT
   edm::ParameterSet psTrackParts_Pt = conf_.getParameter<edm::ParameterSet>("TH1TrackParts_Pt");
   HistoName = "trackParts_Pt";
   trackParts_Pt = iBooker.book1D(HistoName,
                                  HistoName,
                                  psTrackParts_Pt.getParameter<int32_t>("Nbinsx"),
                                  psTrackParts_Pt.getParameter<double>("xmin"),
                                  psTrackParts_Pt.getParameter<double>("xmax"));
   trackParts_Pt->setAxisTitle("p_{T} [GeV]", 1);
   trackParts_Pt->setAxisTitle("# tracking particles", 2);
 
   // 1D: eta
   edm::ParameterSet psTrackParts_Eta = conf_.getParameter<edm::ParameterSet>("TH1TrackParts_Eta");
   HistoName = "trackParts_Eta";
   trackParts_Eta = iBooker.book1D(HistoName,
                                   HistoName,
                                   psTrackParts_Eta.getParameter<int32_t>("Nbinsx"),
                                   psTrackParts_Eta.getParameter<double>("xmin"),
                                   psTrackParts_Eta.getParameter<double>("xmax"));
   trackParts_Eta->setAxisTitle("#eta", 1);
   trackParts_Eta->setAxisTitle("# tracking particles", 2);
 
   // 1D: phi
   edm::ParameterSet psTrackParts_Phi = conf_.getParameter<edm::ParameterSet>("TH1TrackParts_Phi");
   HistoName = "trackParts_Phi";
   trackParts_Phi = iBooker.book1D(HistoName,
                                   HistoName,
                                   psTrackParts_Phi.getParameter<int32_t>("Nbinsx"),
                                   psTrackParts_Phi.getParameter<double>("xmin"),
                                   psTrackParts_Phi.getParameter<double>("xmax"));
   trackParts_Phi->setAxisTitle("#phi", 1);
   trackParts_Phi->setAxisTitle("# tracking particles", 2);
 
   // For tracks correctly matched to truth-level track
   iBooker.setCurrentFolder(topFolderName_ + "/Tracks");
   // chi2
   edm::ParameterSet psTrack_Chi2 = conf_.getParameter<edm::ParameterSet>("TH1_Track_Chi2");
   HistoName = "Track_MatchedChi2";
   Track_MatchedChi2 = iBooker.book1D(HistoName,
                                      HistoName,
                                      psTrack_Chi2.getParameter<int32_t>("Nbinsx"),
                                      psTrack_Chi2.getParameter<double>("xmin"),
                                      psTrack_Chi2.getParameter<double>("xmax"));
   Track_MatchedChi2->setAxisTitle("L1 Track #chi^{2}", 1);
   Track_MatchedChi2->setAxisTitle("# Correctly Matched L1 Tracks", 2);
 
   // chi2Red
   edm::ParameterSet psTrack_Chi2Red = conf_.getParameter<edm::ParameterSet>("TH1_Track_Chi2R");
   HistoName = "Track_MatchedChi2Red";
   Track_MatchedChi2Red = iBooker.book1D(HistoName,
                                         HistoName,
                                         psTrack_Chi2Red.getParameter<int32_t>("Nbinsx"),
                                         psTrack_Chi2Red.getParameter<double>("xmin"),
                                         psTrack_Chi2Red.getParameter<double>("xmax"));
   Track_MatchedChi2Red->setAxisTitle("L1 Track #chi^{2}/ndf", 1);
   Track_MatchedChi2Red->setAxisTitle("# Correctly Matched L1 Tracks", 2);
 
   // 1D plots for efficiency
   iBooker.setCurrentFolder(topFolderName_ + "/Tracks/Efficiency");
   // pT
   edm::ParameterSet psEffic_pt = conf_.getParameter<edm::ParameterSet>("TH1Effic_pt");
   HistoName = "tp_pt";
   tp_pt = iBooker.book1D(HistoName,
                          HistoName,
                          psEffic_pt.getParameter<int32_t>("Nbinsx"),
                          psEffic_pt.getParameter<double>("xmin"),
                          psEffic_pt.getParameter<double>("xmax"));
   tp_pt->setAxisTitle("p_{T} [GeV]", 1);
   tp_pt->setAxisTitle("# tracking particles", 2);
 
   // Matched TP's pT
   HistoName = "match_tp_pt";
   match_tp_pt = iBooker.book1D(HistoName,
                                HistoName,
                                psEffic_pt.getParameter<int32_t>("Nbinsx"),
                                psEffic_pt.getParameter<double>("xmin"),
                                psEffic_pt.getParameter<double>("xmax"));
   match_tp_pt->setAxisTitle("p_{T} [GeV]", 1);
   match_tp_pt->setAxisTitle("# matched tracking particles", 2);
 
   // pT zoom (0-10 GeV)
   edm::ParameterSet psEffic_pt_zoom = conf_.getParameter<edm::ParameterSet>("TH1Effic_pt_zoom");
   HistoName = "tp_pt_zoom";
   tp_pt_zoom = iBooker.book1D(HistoName,
                               HistoName,
                               psEffic_pt_zoom.getParameter<int32_t>("Nbinsx"),
                               psEffic_pt_zoom.getParameter<double>("xmin"),
                               psEffic_pt_zoom.getParameter<double>("xmax"));
   tp_pt_zoom->setAxisTitle("p_{T} [GeV]", 1);
   tp_pt_zoom->setAxisTitle("# tracking particles", 2);
 
   // Matched pT zoom (0-10 GeV)
   HistoName = "match_tp_pt_zoom";
   match_tp_pt_zoom = iBooker.book1D(HistoName,
                                     HistoName,
                                     psEffic_pt_zoom.getParameter<int32_t>("Nbinsx"),
                                     psEffic_pt_zoom.getParameter<double>("xmin"),
                                     psEffic_pt_zoom.getParameter<double>("xmax"));
   match_tp_pt_zoom->setAxisTitle("p_{T} [GeV]", 1);
   match_tp_pt_zoom->setAxisTitle("# matched tracking particles", 2);
 
   // eta
   edm::ParameterSet psEffic_eta = conf_.getParameter<edm::ParameterSet>("TH1Effic_eta");
   HistoName = "tp_eta";
   tp_eta = iBooker.book1D(HistoName,
                           HistoName,
                           psEffic_eta.getParameter<int32_t>("Nbinsx"),
                           psEffic_eta.getParameter<double>("xmin"),
                           psEffic_eta.getParameter<double>("xmax"));
   tp_eta->setAxisTitle("#eta", 1);
   tp_eta->setAxisTitle("# tracking particles", 2);
 
   // Matched eta
   HistoName = "match_tp_eta";
   match_tp_eta = iBooker.book1D(HistoName,
                                 HistoName,
                                 psEffic_eta.getParameter<int32_t>("Nbinsx"),
                                 psEffic_eta.getParameter<double>("xmin"),
                                 psEffic_eta.getParameter<double>("xmax"));
   match_tp_eta->setAxisTitle("#eta", 1);
   match_tp_eta->setAxisTitle("# matched tracking particles", 2);
 
   // d0
   edm::ParameterSet psEffic_d0 = conf_.getParameter<edm::ParameterSet>("TH1Effic_d0");
   HistoName = "tp_d0";
   tp_d0 = iBooker.book1D(HistoName,
                          HistoName,
                          psEffic_d0.getParameter<int32_t>("Nbinsx"),
                          psEffic_d0.getParameter<double>("xmin"),
                          psEffic_d0.getParameter<double>("xmax"));
   tp_d0->setAxisTitle("d_{0} [cm]", 1);
   tp_d0->setAxisTitle("# tracking particles", 2);
 
   // Matched d0
   HistoName = "match_tp_d0";
   match_tp_d0 = iBooker.book1D(HistoName,
                                HistoName,
                                psEffic_d0.getParameter<int32_t>("Nbinsx"),
                                psEffic_d0.getParameter<double>("xmin"),
                                psEffic_d0.getParameter<double>("xmax"));
   match_tp_d0->setAxisTitle("d_{0} [cm]", 1);
   match_tp_d0->setAxisTitle("# matched tracking particles", 2);
 
   // VtxR (also known as vxy)
   edm::ParameterSet psEffic_VtxR = conf_.getParameter<edm::ParameterSet>("TH1Effic_VtxR");
   HistoName = "tp_VtxR";
   tp_VtxR = iBooker.book1D(HistoName,
                            HistoName,
                            psEffic_VtxR.getParameter<int32_t>("Nbinsx"),
                            psEffic_VtxR.getParameter<double>("xmin"),
                            psEffic_VtxR.getParameter<double>("xmax"));
   tp_VtxR->setAxisTitle("d_{xy} [cm]", 1);
   tp_VtxR->setAxisTitle("# tracking particles", 2);
 
   // Matched VtxR
   HistoName = "match_tp_VtxR";
   match_tp_VtxR = iBooker.book1D(HistoName,
                                  HistoName,
                                  psEffic_VtxR.getParameter<int32_t>("Nbinsx"),
                                  psEffic_VtxR.getParameter<double>("xmin"),
                                  psEffic_VtxR.getParameter<double>("xmax"));
   match_tp_VtxR->setAxisTitle("d_{xy} [cm]", 1);
   match_tp_VtxR->setAxisTitle("# matched tracking particles", 2);
 
   // VtxZ
   edm::ParameterSet psEffic_VtxZ = conf_.getParameter<edm::ParameterSet>("TH1Effic_VtxZ");
   HistoName = "tp_VtxZ";
   tp_VtxZ = iBooker.book1D(HistoName,
                            HistoName,
                            psEffic_VtxZ.getParameter<int32_t>("Nbinsx"),
                            psEffic_VtxZ.getParameter<double>("xmin"),
                            psEffic_VtxZ.getParameter<double>("xmax"));
   tp_VtxZ->setAxisTitle("z_{0} [cm]", 1);
   tp_VtxZ->setAxisTitle("# tracking particles", 2);
 
   // Matched d0
   HistoName = "match_tp_VtxZ";
   match_tp_VtxZ = iBooker.book1D(HistoName,
                                  HistoName,
                                  psEffic_VtxZ.getParameter<int32_t>("Nbinsx"),
                                  psEffic_VtxZ.getParameter<double>("xmin"),
                                  psEffic_VtxZ.getParameter<double>("xmax"));
   match_tp_VtxZ->setAxisTitle("z_{0} [cm]", 1);
   match_tp_VtxZ->setAxisTitle("# matched tracking particles", 2);
 
   // 1D plots for resolution
   iBooker.setCurrentFolder(topFolderName_ + "/Tracks/Resolution");
   // full pT
   edm::ParameterSet psRes_pt = conf_.getParameter<edm::ParameterSet>("TH1Res_pt");
   HistoName = "res_pt";
   res_pt = iBooker.book1D(HistoName,
                           HistoName,
                           psRes_pt.getParameter<int32_t>("Nbinsx"),
                           psRes_pt.getParameter<double>("xmin"),
                           psRes_pt.getParameter<double>("xmax"));
   res_pt->setAxisTitle("p_{T} [GeV]", 1);
   res_pt->setAxisTitle("# tracking particles", 2);
 
   // Full eta
   edm::ParameterSet psRes_eta = conf_.getParameter<edm::ParameterSet>("TH1Res_eta");
   HistoName = "res_eta";
   res_eta = iBooker.book1D(HistoName,
                            HistoName,
                            psRes_eta.getParameter<int32_t>("Nbinsx"),
                            psRes_eta.getParameter<double>("xmin"),
                            psRes_eta.getParameter<double>("xmax"));
   res_eta->setAxisTitle("#eta", 1);
   res_eta->setAxisTitle("# tracking particles", 2);
 
   // Relative pT
   edm::ParameterSet psRes_ptRel = conf_.getParameter<edm::ParameterSet>("TH1Res_ptRel");
   HistoName = "res_ptRel";
   res_ptRel = iBooker.book1D(HistoName,
                              HistoName,
                              psRes_ptRel.getParameter<int32_t>("Nbinsx"),
                              psRes_ptRel.getParameter<double>("xmin"),
                              psRes_ptRel.getParameter<double>("xmax"));
   res_ptRel->setAxisTitle("Relative p_{T} [GeV]", 1);
   res_ptRel->setAxisTitle("# tracking particles", 2);
 
   // Eta parts (for resolution)
   // Eta 1 (0 to 0.7)
   HistoName = "reseta_eta0to0p7";
   reseta_eta0to0p7 = iBooker.book1D(HistoName,
                                     HistoName,
                                     psRes_eta.getParameter<int32_t>("Nbinsx"),
                                     psRes_eta.getParameter<double>("xmin"),
                                     psRes_eta.getParameter<double>("xmax"));
   reseta_eta0to0p7->setAxisTitle("#eta_{trk} - #eta_{tp}", 1);
   reseta_eta0to0p7->setAxisTitle("# tracking particles", 2);
 
   // Eta 2 (0.7 to 1.0)
   HistoName = "reseta_eta0p7to1";
   reseta_eta0p7to1 = iBooker.book1D(HistoName,
                                     HistoName,
                                     psRes_eta.getParameter<int32_t>("Nbinsx"),
                                     psRes_eta.getParameter<double>("xmin"),
                                     psRes_eta.getParameter<double>("xmax"));
   reseta_eta0p7to1->setAxisTitle("#eta_{trk} - #eta_{tp}", 1);
   reseta_eta0p7to1->setAxisTitle("# tracking particles", 2);
 
   // Eta 3 (1.0 to 1.2)
   HistoName = "reseta_eta1to1p2";
   reseta_eta1to1p2 = iBooker.book1D(HistoName,
                                     HistoName,
                                     psRes_eta.getParameter<int32_t>("Nbinsx"),
                                     psRes_eta.getParameter<double>("xmin"),
                                     psRes_eta.getParameter<double>("xmax"));
   reseta_eta1to1p2->setAxisTitle("#eta_{trk} - #eta_{tp}", 1);
   reseta_eta1to1p2->setAxisTitle("# tracking particles", 2);
 
   // Eta 4 (1.2 to 1.6)
   HistoName = "reseta_eta1p2to1p6";
   reseta_eta1p2to1p6 = iBooker.book1D(HistoName,
                                       HistoName,
                                       psRes_eta.getParameter<int32_t>("Nbinsx"),
                                       psRes_eta.getParameter<double>("xmin"),
                                       psRes_eta.getParameter<double>("xmax"));
   reseta_eta1p2to1p6->setAxisTitle("#eta_{trk} - #eta_{tp}", 1);
   reseta_eta1p2to1p6->setAxisTitle("# tracking particles", 2);
 
   // Eta 5 (1.6 to 2.0)
   HistoName = "reseta_eta1p6to2";
   reseta_eta1p6to2 = iBooker.book1D(HistoName,
                                     HistoName,
                                     psRes_eta.getParameter<int32_t>("Nbinsx"),
                                     psRes_eta.getParameter<double>("xmin"),
                                     psRes_eta.getParameter<double>("xmax"));
   reseta_eta1p6to2->setAxisTitle("#eta_{trk} - #eta_{tp}", 1);
   reseta_eta1p6to2->setAxisTitle("# tracking particles", 2);
 
   // Eta 6 (2.0 to 2.4)
   HistoName = "reseta_eta2to2p4";
   reseta_eta2to2p4 = iBooker.book1D(HistoName,
                                     HistoName,
                                     psRes_eta.getParameter<int32_t>("Nbinsx"),
                                     psRes_eta.getParameter<double>("xmin"),
                                     psRes_eta.getParameter<double>("xmax"));
   reseta_eta2to2p4->setAxisTitle("#eta_{trk} - #eta_{tp}", 1);
   reseta_eta2to2p4->setAxisTitle("# tracking particles", 2);
 
   // pT parts for resolution (pT res vs eta)
   // pT a (2 to 3 GeV)
   // Eta 1 (0 to 0.7)
   HistoName = "respt_eta0to0p7_pt2to3";
   respt_eta0to0p7_pt2to3 = iBooker.book1D(HistoName,
                                           HistoName,
                                           psRes_pt.getParameter<int32_t>("Nbinsx"),
                                           psRes_pt.getParameter<double>("xmin"),
                                           psRes_pt.getParameter<double>("xmax"));
   respt_eta0to0p7_pt2to3->setAxisTitle("(p_{T}(trk) - p_{T}(tp))/p_{T}(tp)", 1);
   respt_eta0to0p7_pt2to3->setAxisTitle("# tracking particles", 2);
 
   // Eta 2 (0.7 to 1.0)
   HistoName = "respt_eta0p7to1_pt2to3";
   respt_eta0p7to1_pt2to3 = iBooker.book1D(HistoName,
                                           HistoName,
                                           psRes_pt.getParameter<int32_t>("Nbinsx"),
                                           psRes_pt.getParameter<double>("xmin"),
                                           psRes_pt.getParameter<double>("xmax"));
   respt_eta0p7to1_pt2to3->setAxisTitle("(p_{T}(trk) - p_{T}(tp))/p_{T}(tp)", 1);
   respt_eta0p7to1_pt2to3->setAxisTitle("# tracking particles", 2);
 
   // Eta 3 (1.0 to 1.2)
   HistoName = "respt_eta1to1p2_pt2to3";
   respt_eta1to1p2_pt2to3 = iBooker.book1D(HistoName,
                                           HistoName,
                                           psRes_pt.getParameter<int32_t>("Nbinsx"),
                                           psRes_pt.getParameter<double>("xmin"),
                                           psRes_pt.getParameter<double>("xmax"));
   respt_eta1to1p2_pt2to3->setAxisTitle("(p_{T}(trk) - p_{T}(tp))/p_{T}(tp)", 1);
   respt_eta1to1p2_pt2to3->setAxisTitle("# tracking particles", 2);
 
   // Eta 4 (1.2 to 1.6)
   HistoName = "respt_eta1p2to1p6_pt2to3";
   respt_eta1p2to1p6_pt2to3 = iBooker.book1D(HistoName,
                                             HistoName,
                                             psRes_pt.getParameter<int32_t>("Nbinsx"),
                                             psRes_pt.getParameter<double>("xmin"),
                                             psRes_pt.getParameter<double>("xmax"));
   respt_eta1p2to1p6_pt2to3->setAxisTitle("(p_{T}(trk) - p_{T}(tp))/p_{T}(tp)", 1);
   respt_eta1p2to1p6_pt2to3->setAxisTitle("# tracking particles", 2);
 
   // Eta 5 (1.6 to 2.0)
   HistoName = "respt_eta1p6to2_pt2to3";
   respt_eta1p6to2_pt2to3 = iBooker.book1D(HistoName,
                                           HistoName,
                                           psRes_pt.getParameter<int32_t>("Nbinsx"),
                                           psRes_pt.getParameter<double>("xmin"),
                                           psRes_pt.getParameter<double>("xmax"));
   respt_eta1p6to2_pt2to3->setAxisTitle("(p_{T}(trk) - p_{T}(tp))/p_{T}(tp)", 1);
   respt_eta1p6to2_pt2to3->setAxisTitle("# tracking particles", 2);
 
   // Eta 6 (2.0 to 2.4)
   HistoName = "respt_eta2to2p4_pt2to3";
   respt_eta2to2p4_pt2to3 = iBooker.book1D(HistoName,
                                           HistoName,
                                           psRes_pt.getParameter<int32_t>("Nbinsx"),
                                           psRes_pt.getParameter<double>("xmin"),
                                           psRes_pt.getParameter<double>("xmax"));
   respt_eta2to2p4_pt2to3->setAxisTitle("(p_{T}(trk) - p_{T}(tp))/p_{T}(tp)", 1);
   respt_eta2to2p4_pt2to3->setAxisTitle("# tracking particles", 2);
 
   // pT b (3 to 8 GeV)
   // Eta 1 (0 to 0.7)
   HistoName = "respt_eta0to0p7_pt3to8";
   respt_eta0to0p7_pt3to8 = iBooker.book1D(HistoName,
                                           HistoName,
                                           psRes_pt.getParameter<int32_t>("Nbinsx"),
                                           psRes_pt.getParameter<double>("xmin"),
                                           psRes_pt.getParameter<double>("xmax"));
   respt_eta0to0p7_pt3to8->setAxisTitle("(p_{T}(trk) - p_{T}(tp))/p_{T}(tp)", 1);
   respt_eta0to0p7_pt3to8->setAxisTitle("# tracking particles", 2);
 
   // Eta 2 (0.7 to 1.0)
   HistoName = "respt_eta0p7to1_pt3to8";
   respt_eta0p7to1_pt3to8 = iBooker.book1D(HistoName,
                                           HistoName,
                                           psRes_pt.getParameter<int32_t>("Nbinsx"),
                                           psRes_pt.getParameter<double>("xmin"),
                                           psRes_pt.getParameter<double>("xmax"));
   respt_eta0p7to1_pt3to8->setAxisTitle("(p_{T}(trk) - p_{T}(tp))/p_{T}(tp)", 1);
   respt_eta0p7to1_pt3to8->setAxisTitle("# tracking particles", 2);
 
   // Eta 3 (1.0 to 1.2)
   HistoName = "respt_eta1to1p2_pt3to8";
   respt_eta1to1p2_pt3to8 = iBooker.book1D(HistoName,
                                           HistoName,
                                           psRes_pt.getParameter<int32_t>("Nbinsx"),
                                           psRes_pt.getParameter<double>("xmin"),
                                           psRes_pt.getParameter<double>("xmax"));
   respt_eta1to1p2_pt3to8->setAxisTitle("(p_{T}(trk) - p_{T}(tp))/p_{T}(tp)", 1);
   respt_eta1to1p2_pt3to8->setAxisTitle("# tracking particles", 2);
 
   // Eta 4 (1.2 to 1.6)
   HistoName = "respt_eta1p2to1p6_pt3to8";
   respt_eta1p2to1p6_pt3to8 = iBooker.book1D(HistoName,
                                             HistoName,
                                             psRes_pt.getParameter<int32_t>("Nbinsx"),
                                             psRes_pt.getParameter<double>("xmin"),
                                             psRes_pt.getParameter<double>("xmax"));
   respt_eta1p2to1p6_pt3to8->setAxisTitle("(p_{T}(trk) - p_{T}(tp))/p_{T}(tp)", 1);
   respt_eta1p2to1p6_pt3to8->setAxisTitle("# tracking particles", 2);
 
   // Eta 5 (1.6 to 2.0)
   HistoName = "respt_eta1p6to2_pt3to8";
   respt_eta1p6to2_pt3to8 = iBooker.book1D(HistoName,
                                           HistoName,
                                           psRes_pt.getParameter<int32_t>("Nbinsx"),
                                           psRes_pt.getParameter<double>("xmin"),
                                           psRes_pt.getParameter<double>("xmax"));
   respt_eta1p6to2_pt3to8->setAxisTitle("(p_{T}(trk) - p_{T}(tp))/p_{T}(tp)", 1);
   respt_eta1p6to2_pt3to8->setAxisTitle("# tracking particles", 2);
 
   // Eta 6 (2.0 to 2.4)
   HistoName = "respt_eta2to2p4_pt3to8";
   respt_eta2to2p4_pt3to8 = iBooker.book1D(HistoName,
                                           HistoName,
                                           psRes_pt.getParameter<int32_t>("Nbinsx"),
                                           psRes_pt.getParameter<double>("xmin"),
                                           psRes_pt.getParameter<double>("xmax"));
   respt_eta2to2p4_pt3to8->setAxisTitle("(p_{T}(trk) - p_{T}(tp))/p_{T}(tp)", 1);
   respt_eta2to2p4_pt3to8->setAxisTitle("# tracking particles", 2);
 
   // pT c (>8 GeV)
   // Eta 1 (0 to 0.7)
   HistoName = "respt_eta0to0p7_pt8toInf";
   respt_eta0to0p7_pt8toInf = iBooker.book1D(HistoName,
                                             HistoName,
                                             psRes_pt.getParameter<int32_t>("Nbinsx"),
                                             psRes_pt.getParameter<double>("xmin"),
                                             psRes_pt.getParameter<double>("xmax"));
   respt_eta0to0p7_pt8toInf->setAxisTitle("(p_{T}(trk) - p_{T}(tp))/p_{T}(tp)", 1);
   respt_eta0to0p7_pt8toInf->setAxisTitle("# tracking particles", 2);
 
   // Eta 2 (0.7 to 1.0)
   HistoName = "respt_eta0p7to1_pt8toInf";
   respt_eta0p7to1_pt8toInf = iBooker.book1D(HistoName,
                                             HistoName,
                                             psRes_pt.getParameter<int32_t>("Nbinsx"),
                                             psRes_pt.getParameter<double>("xmin"),
                                             psRes_pt.getParameter<double>("xmax"));
   respt_eta0p7to1_pt8toInf->setAxisTitle("(p_{T}(trk) - p_{T}(tp))/p_{T}(tp)", 1);
   respt_eta0p7to1_pt8toInf->setAxisTitle("# tracking particles", 2);
 
   // Eta 3 (1.0 to 1.2)
   HistoName = "respt_eta1to1p2_pt8toInf";
   respt_eta1to1p2_pt8toInf = iBooker.book1D(HistoName,
                                             HistoName,
                                             psRes_pt.getParameter<int32_t>("Nbinsx"),
                                             psRes_pt.getParameter<double>("xmin"),
                                             psRes_pt.getParameter<double>("xmax"));
   respt_eta1to1p2_pt8toInf->setAxisTitle("(p_{T}(trk) - p_{T}(tp))/p_{T}(tp)", 1);
   respt_eta1to1p2_pt8toInf->setAxisTitle("# tracking particles", 2);
 
   // Eta 4 (1.2 to 1.6)
   HistoName = "respt_eta1p2to1p6_pt8toInf";
   respt_eta1p2to1p6_pt8toInf = iBooker.book1D(HistoName,
                                               HistoName,
                                               psRes_pt.getParameter<int32_t>("Nbinsx"),
                                               psRes_pt.getParameter<double>("xmin"),
                                               psRes_pt.getParameter<double>("xmax"));
   respt_eta1p2to1p6_pt8toInf->setAxisTitle("(p_{T}(trk) - p_{T}(tp))/p_{T}(tp)", 1);
   respt_eta1p2to1p6_pt8toInf->setAxisTitle("# tracking particles", 2);
 
   // Eta 5 (1.6 to 2.0)
   HistoName = "respt_eta1p6to2_pt8toInf";
   respt_eta1p6to2_pt8toInf = iBooker.book1D(HistoName,
                                             HistoName,
                                             psRes_pt.getParameter<int32_t>("Nbinsx"),
                                             psRes_pt.getParameter<double>("xmin"),
                                             psRes_pt.getParameter<double>("xmax"));
   respt_eta1p6to2_pt8toInf->setAxisTitle("(p_{T}(trk) - p_{T}(tp))/p_{T}(tp)", 1);
   respt_eta1p6to2_pt8toInf->setAxisTitle("# tracking particles", 2);
 
   // Eta 6 (2.0 to 2.4)
   HistoName = "respt_eta2to2p4_pt8toInf";
   respt_eta2to2p4_pt8toInf = iBooker.book1D(HistoName,
                                             HistoName,
                                             psRes_pt.getParameter<int32_t>("Nbinsx"),
                                             psRes_pt.getParameter<double>("xmin"),
                                             psRes_pt.getParameter<double>("xmax"));
   respt_eta2to2p4_pt8toInf->setAxisTitle("(p_{T}(trk) - p_{T}(tp))/p_{T}(tp)", 1);
   respt_eta2to2p4_pt8toInf->setAxisTitle("# tracking particles", 2);
 
   // Phi parts (for resolution)
   // Eta 1 (0 to 0.7)
   edm::ParameterSet psRes_phi = conf_.getParameter<edm::ParameterSet>("TH1Res_phi");
   HistoName = "resphi_eta0to0p7";
   resphi_eta0to0p7 = iBooker.book1D(HistoName,
                                     HistoName,
                                     psRes_phi.getParameter<int32_t>("Nbinsx"),
                                     psRes_phi.getParameter<double>("xmin"),
                                     psRes_phi.getParameter<double>("xmax"));
   resphi_eta0to0p7->setAxisTitle("#phi_{trk} - #phi_{tp}", 1);
   resphi_eta0to0p7->setAxisTitle("# tracking particles", 2);
 
   // Eta 2 (0.7 to 1.0)
   HistoName = "resphi_eta0p7to1";
   resphi_eta0p7to1 = iBooker.book1D(HistoName,
                                     HistoName,
                                     psRes_phi.getParameter<int32_t>("Nbinsx"),
                                     psRes_phi.getParameter<double>("xmin"),
                                     psRes_phi.getParameter<double>("xmax"));
   resphi_eta0p7to1->setAxisTitle("#phi_{trk} - #phi_{tp}", 1);
   resphi_eta0p7to1->setAxisTitle("# tracking particles", 2);
 
   // Eta 3 (1.0 to 1.2)
   HistoName = "resphi_eta1to1p2";
   resphi_eta1to1p2 = iBooker.book1D(HistoName,
                                     HistoName,
                                     psRes_phi.getParameter<int32_t>("Nbinsx"),
                                     psRes_phi.getParameter<double>("xmin"),
                                     psRes_phi.getParameter<double>("xmax"));
   resphi_eta1to1p2->setAxisTitle("#phi_{trk} - #phi_{tp}", 1);
   resphi_eta1to1p2->setAxisTitle("# tracking particles", 2);
 
   // Eta 4 (1.2 to 1.6)
   HistoName = "resphi_eta1p2to1p6";
   resphi_eta1p2to1p6 = iBooker.book1D(HistoName,
                                       HistoName,
                                       psRes_phi.getParameter<int32_t>("Nbinsx"),
                                       psRes_phi.getParameter<double>("xmin"),
                                       psRes_phi.getParameter<double>("xmax"));
   resphi_eta1p2to1p6->setAxisTitle("#phi_{trk} - #phi_{tp}", 1);
   resphi_eta1p2to1p6->setAxisTitle("# tracking particles", 2);
 
   // Eta 5 (1.6 to 2.0)
   HistoName = "resphi_eta1p6to2";
   resphi_eta1p6to2 = iBooker.book1D(HistoName,
                                     HistoName,
                                     psRes_phi.getParameter<int32_t>("Nbinsx"),
                                     psRes_phi.getParameter<double>("xmin"),
                                     psRes_phi.getParameter<double>("xmax"));
   resphi_eta1p6to2->setAxisTitle("#phi_{trk} - #phi_{tp}", 1);
   resphi_eta1p6to2->setAxisTitle("# tracking particles", 2);
 
   // Eta 6 (2.0 to 2.4)
   HistoName = "resphi_eta2to2p4";
   resphi_eta2to2p4 = iBooker.book1D(HistoName,
                                     HistoName,
                                     psRes_phi.getParameter<int32_t>("Nbinsx"),
                                     psRes_phi.getParameter<double>("xmin"),
                                     psRes_phi.getParameter<double>("xmax"));
   resphi_eta2to2p4->setAxisTitle("#phi_{trk} - #phi_{tp}", 1);
   resphi_eta2to2p4->setAxisTitle("# tracking particles", 2);
 
   // VtxZ parts (for resolution)
   // Eta 1 (0 to 0.7)
   edm::ParameterSet psRes_VtxZ = conf_.getParameter<edm::ParameterSet>("TH1Res_VtxZ");
   HistoName = "resVtxZ_eta0to0p7";
   resVtxZ_eta0to0p7 = iBooker.book1D(HistoName,
                                      HistoName,
                                      psRes_VtxZ.getParameter<int32_t>("Nbinsx"),
                                      psRes_VtxZ.getParameter<double>("xmin"),
                                      psRes_VtxZ.getParameter<double>("xmax"));
   resVtxZ_eta0to0p7->setAxisTitle("VtxZ_{trk} - VtxZ_{tp} [cm]", 1);
   resVtxZ_eta0to0p7->setAxisTitle("# tracking particles", 2);
 
   // Eta 2 (0.7 to 1.0)
   HistoName = "resVtxZ_eta0p7to1";
   resVtxZ_eta0p7to1 = iBooker.book1D(HistoName,
                                      HistoName,
                                      psRes_VtxZ.getParameter<int32_t>("Nbinsx"),
                                      psRes_VtxZ.getParameter<double>("xmin"),
                                      psRes_VtxZ.getParameter<double>("xmax"));
   resVtxZ_eta0p7to1->setAxisTitle("VtxZ_{trk} - VtxZ_{tp} [cm]", 1);
   resVtxZ_eta0p7to1->setAxisTitle("# tracking particles", 2);
 
   // Eta 3 (1.0 to 1.2)
   HistoName = "resVtxZ_eta1to1p2";
   resVtxZ_eta1to1p2 = iBooker.book1D(HistoName,
                                      HistoName,
                                      psRes_VtxZ.getParameter<int32_t>("Nbinsx"),
                                      psRes_VtxZ.getParameter<double>("xmin"),
                                      psRes_VtxZ.getParameter<double>("xmax"));
   resVtxZ_eta1to1p2->setAxisTitle("VtxZ_{trk} - VtxZ_{tp} [cm]", 1);
   resVtxZ_eta1to1p2->setAxisTitle("# tracking particles", 2);
 
   // Eta 4 (1.2 to 1.6)
   HistoName = "resVtxZ_eta1p2to1p6";
   resVtxZ_eta1p2to1p6 = iBooker.book1D(HistoName,
                                        HistoName,
                                        psRes_VtxZ.getParameter<int32_t>("Nbinsx"),
                                        psRes_VtxZ.getParameter<double>("xmin"),
                                        psRes_VtxZ.getParameter<double>("xmax"));
   resVtxZ_eta1p2to1p6->setAxisTitle("VtxZ_{trk} - VtxZ_{tp} [cm]", 1);
   resVtxZ_eta1p2to1p6->setAxisTitle("# tracking particles", 2);
 
   // Eta 5 (1.6 to 2.0)
   HistoName = "resVtxZ_eta1p6to2";
   resVtxZ_eta1p6to2 = iBooker.book1D(HistoName,
                                      HistoName,
                                      psRes_VtxZ.getParameter<int32_t>("Nbinsx"),
                                      psRes_VtxZ.getParameter<double>("xmin"),
                                      psRes_VtxZ.getParameter<double>("xmax"));
   resVtxZ_eta1p6to2->setAxisTitle("VtxZ_{trk} - VtxZ_{tp} [cm]", 1);
   resVtxZ_eta1p6to2->setAxisTitle("# tracking particles", 2);
 
   // Eta 6 (2.0 to 2.4)
   HistoName = "resVtxZ_eta2to2p4";
   resVtxZ_eta2to2p4 = iBooker.book1D(HistoName,
                                      HistoName,
                                      psRes_VtxZ.getParameter<int32_t>("Nbinsx"),
                                      psRes_VtxZ.getParameter<double>("xmin"),
                                      psRes_VtxZ.getParameter<double>("xmax"));
   resVtxZ_eta2to2p4->setAxisTitle("VtxZ_{trk} - VtxZ_{tp} [cm]", 1);
   resVtxZ_eta2to2p4->setAxisTitle("# tracking particles", 2);
 
   // d0 parts (for resolution)
   // Eta 1 (0 to 0.7)
   edm::ParameterSet psRes_d0 = conf_.getParameter<edm::ParameterSet>("TH1Res_d0");
   HistoName = "resd0_eta0to0p7";
   resd0_eta0to0p7 = iBooker.book1D(HistoName,
                                    HistoName,
                                    psRes_d0.getParameter<int32_t>("Nbinsx"),
                                    psRes_d0.getParameter<double>("xmin"),
                                    psRes_d0.getParameter<double>("xmax"));
   resd0_eta0to0p7->setAxisTitle("d0_{trk} - d0_{tp} [cm]", 1);
   resd0_eta0to0p7->setAxisTitle("# tracking particles", 2);
 
   // Eta 2 (0.7 to 1.0)
   HistoName = "resd0_eta0p7to1";
   resd0_eta0p7to1 = iBooker.book1D(HistoName,
                                    HistoName,
                                    psRes_d0.getParameter<int32_t>("Nbinsx"),
                                    psRes_d0.getParameter<double>("xmin"),
                                    psRes_d0.getParameter<double>("xmax"));
   resd0_eta0p7to1->setAxisTitle("d0_{trk} - d0_{tp} [cm]", 1);
   resd0_eta0p7to1->setAxisTitle("# tracking particles", 2);
 
   // Eta 3 (1.0 to 1.2)
   HistoName = "resd0_eta1to1p2";
   resd0_eta1to1p2 = iBooker.book1D(HistoName,
                                    HistoName,
                                    psRes_d0.getParameter<int32_t>("Nbinsx"),
                                    psRes_d0.getParameter<double>("xmin"),
                                    psRes_d0.getParameter<double>("xmax"));
   resd0_eta1to1p2->setAxisTitle("d0_{trk} - d0_{tp} [cm]", 1);
   resd0_eta1to1p2->setAxisTitle("# tracking particles", 2);
 
   // Eta 4 (1.2 to 1.6)
   HistoName = "resd0_eta1p2to1p6";
   resd0_eta1p2to1p6 = iBooker.book1D(HistoName,
                                      HistoName,
                                      psRes_d0.getParameter<int32_t>("Nbinsx"),
                                      psRes_d0.getParameter<double>("xmin"),
                                      psRes_d0.getParameter<double>("xmax"));
   resd0_eta1p2to1p6->setAxisTitle("d0_{trk} - d0_{tp} [cm]", 1);
   resd0_eta1p2to1p6->setAxisTitle("# tracking particles", 2);
 
   // Eta 5 (1.6 to 2.0)
   HistoName = "resd0_eta1p6to2";
   resd0_eta1p6to2 = iBooker.book1D(HistoName,
                                    HistoName,
                                    psRes_d0.getParameter<int32_t>("Nbinsx"),
                                    psRes_d0.getParameter<double>("xmin"),
                                    psRes_d0.getParameter<double>("xmax"));
   resd0_eta1p6to2->setAxisTitle("d0_{trk} - d0_{tp} [cm]", 1);
   resd0_eta1p6to2->setAxisTitle("# tracking particles", 2);
 
   // Eta 6 (2.0 to 2.4)
   HistoName = "resd0_eta2to2p4";
   resd0_eta2to2p4 = iBooker.book1D(HistoName,
                                    HistoName,
                                    psRes_d0.getParameter<int32_t>("Nbinsx"),
                                    psRes_d0.getParameter<double>("xmin"),
                                    psRes_d0.getParameter<double>("xmax"));
   resd0_eta2to2p4->setAxisTitle("d0_{trk} - d0_{tp} [cm]", 1);
   resd0_eta2to2p4->setAxisTitle("# tracking particles", 2);
 
 }  // end of method