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

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

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

	~DQMEDAnalyzer () override=default

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_

Detailed Description

Definition at line 26 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
   TP_minNLayersStub= conf_.getParameter< int >("TP_minNLayersStub"); //min number of layers with stubs in the tracking particle to consider matching                                                                  // consider matching
   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 63 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 71 of file OuterTrackerMonitorTrackingParticles.cc.

References patCaloMETCorrections_cff::A, TrackingParticle::charge(), funct::cos(), JetChargeProducer_cfi::exp, 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
	)

override

Definition at line 385 of file OuterTrackerMonitorTrackingParticles.cc.

References 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, MonitorElement::setAxisTitle(), 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