#include <OuterTrackerMonitorTTTrack.h>

Inheritance diagram for OuterTrackerMonitorTTTrack:

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

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

	OuterTrackerMonitorTTTrack (const edm::ParameterSet &)

	~OuterTrackerMonitorTTTrack () override

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

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override

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

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

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

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

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

void	endRun (edm::Run const &run, edm::EventSetup const &setup) override

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

	~DQMEDAnalyzer () override=default

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

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

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

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Public Attributes
MonitorElement *	Track_Eta_NStubs = 0

MonitorElement *	Track_HQ_Chi2 = 0

MonitorElement *	Track_HQ_Chi2_Probability = 0

MonitorElement *	Track_HQ_Chi2Red = 0

MonitorElement *	Track_HQ_Chi2Red_Eta = 0

MonitorElement *	Track_HQ_Chi2Red_NStubs = 0

MonitorElement *	Track_HQ_D0 = 0

MonitorElement *	Track_HQ_Eta = 0

MonitorElement *	Track_HQ_Eta_BarrelStubs = 0

MonitorElement *	Track_HQ_Eta_ECStubs = 0

MonitorElement *	Track_HQ_N = 0
	High-quality TTTracks (NStubs >=5, chi2/dof<10) More...

MonitorElement *	Track_HQ_Phi = 0

MonitorElement *	Track_HQ_Pt = 0

MonitorElement *	Track_HQ_VtxZ = 0

MonitorElement *	Track_LQ_Chi2 = 0

MonitorElement *	Track_LQ_Chi2_Probability = 0

MonitorElement *	Track_LQ_Chi2Red = 0

MonitorElement *	Track_LQ_Chi2Red_Eta = 0

MonitorElement *	Track_LQ_Chi2Red_NStubs = 0

MonitorElement *	Track_LQ_D0 = 0

MonitorElement *	Track_LQ_Eta = 0

MonitorElement *	Track_LQ_Eta_BarrelStubs = 0

MonitorElement *	Track_LQ_Eta_ECStubs = 0

MonitorElement *	Track_LQ_N = 0
	Low-quality TTTracks (All tracks) More...

MonitorElement *	Track_LQ_Phi = 0

MonitorElement *	Track_LQ_Pt = 0

MonitorElement *	Track_LQ_VtxZ = 0

MonitorElement *	Track_NStubs = 0

Private Attributes
edm::ParameterSet	conf_

double	HQChi2dof_

unsigned int	HQNStubs_

std::string	topFolderName_

edm::EDGetTokenT< std::vector< TTTrack< Ref_Phase2TrackerDigi_ > > >	ttTrackToken_

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

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

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 22 of file OuterTrackerMonitorTTTrack.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 39 of file OuterTrackerMonitorTTTrack.cc.

References conf_, edm::ParameterSet::getParameter(), HQChi2dof_, HQNStubs_, AlCaHLTBitMon_QueryRunRegistry::string, topFolderName_, and ttTrackToken_.

 :conf_(iConfig)
 {
   topFolderName_ = conf_.getParameter<std::string>("TopFolderName");
   ttTrackToken_ = consumes< std::vector < TTTrack< Ref_Phase2TrackerDigi_ > > > (conf_.getParameter<edm::InputTag>("TTTracksTag"));
   HQNStubs_        = conf_.getParameter<int>("HQNStubs");
   HQChi2dof_       = conf_.getParameter<double>("HQChi2dof");
 }

OuterTrackerMonitorTTTrack::~OuterTrackerMonitorTTTrack ( )

override

Definition at line 48 of file OuterTrackerMonitorTTTrack.cc.

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

Member Function Documentation

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

override

Track Trigger Tracks

Loop over TTTracks

Make the pointer

Definition at line 57 of file OuterTrackerMonitorTTTrack.cc.

References ChiSquaredProbability(), MonitorElement::Fill(), edm::Event::getByToken(), HQChi2dof_, HQNStubs_, edm::HandleBase::isValid(), mathSSE::sqrt(), StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, Track_Eta_NStubs, Track_HQ_Chi2, Track_HQ_Chi2_Probability, Track_HQ_Chi2Red, Track_HQ_Chi2Red_Eta, Track_HQ_Chi2Red_NStubs, Track_HQ_D0, Track_HQ_Eta, Track_HQ_Eta_BarrelStubs, Track_HQ_Eta_ECStubs, Track_HQ_N, Track_HQ_Phi, Track_HQ_Pt, Track_HQ_VtxZ, Track_LQ_Chi2, Track_LQ_Chi2_Probability, Track_LQ_Chi2Red, Track_LQ_Chi2Red_Eta, Track_LQ_Chi2Red_NStubs, Track_LQ_D0, Track_LQ_Eta, Track_LQ_Eta_BarrelStubs, Track_LQ_Eta_ECStubs, Track_LQ_N, Track_LQ_Phi, Track_LQ_Pt, Track_LQ_VtxZ, Track_NStubs, reco::btau::trackChi2, align::Tracker, and ttTrackToken_.

 {
   // L1 Primaries
   edm::Handle< std::vector< TTTrack< Ref_Phase2TrackerDigi_ > > > TTTrackHandle;
   iEvent.getByToken(ttTrackToken_, TTTrackHandle);
 
   unsigned int numHQTracks = 0;
   unsigned int numLQTracks = 0;
 
   // Adding protection
   if ( !TTTrackHandle.isValid() ) return;
 
   unsigned int tkCnt = 0;
   for(auto iterTTTrack: *TTTrackHandle){
     edm::Ptr< TTTrack< Ref_Phase2TrackerDigi_ > > tempTrackPtr(TTTrackHandle, tkCnt++); 
 
     unsigned int nStubs = tempTrackPtr->getStubRefs().size();
     int nBarrelStubs = 0;
     int nECStubs = 0;
 
 
     double track_d0 = sqrt(tempTrackPtr->getPOCA().x()*tempTrackPtr->getPOCA().x() + tempTrackPtr->getPOCA().y()*tempTrackPtr->getPOCA().y());
     double trackChi2 = tempTrackPtr->getChi2();
     double trackChi2R = tempTrackPtr->getChi2Red();
 
     Track_NStubs->Fill(nStubs);
     Track_Eta_NStubs->Fill(tempTrackPtr->getMomentum().eta(), nStubs);
 
     std::vector< edm::Ref< edmNew::DetSetVector< TTStub< Ref_Phase2TrackerDigi_ > >, TTStub< Ref_Phase2TrackerDigi_ > > >  theStubs = iterTTTrack.getStubRefs() ;
     for (auto istub: theStubs) {
       bool inTIB = false;
       bool inTOB = false;
       bool inTEC = false;
       bool inTID = false;
 
       // Verify if stubs are in EC or barrel
       DetId detId(istub->getDetId());
       if (detId.det() == DetId::Detector::Tracker) {
 
         if      (detId.subdetId() == StripSubdetector::TIB) inTIB=true;
         else if (detId.subdetId() == StripSubdetector::TOB) inTOB=true;
         else if (detId.subdetId() == StripSubdetector::TEC) inTEC=true;
         else if (detId.subdetId() == StripSubdetector::TID) inTID=true;
       }
       if (inTIB) nBarrelStubs++;
       else if (inTOB) nBarrelStubs++;
       else if (inTEC) nECStubs++;
       else if (inTID) nECStubs++;
     } //end loop over stubs
 
     //HQ tracks: >=5 stubs and chi2/dof < 1
     if (nStubs >= HQNStubs_ && trackChi2R <= HQChi2dof_) {
       numHQTracks++;
 
       Track_HQ_Pt->Fill(tempTrackPtr->getMomentum().perp());
       Track_HQ_Eta->Fill(tempTrackPtr->getMomentum().eta());
       Track_HQ_Phi->Fill(tempTrackPtr->getMomentum().phi());
       Track_HQ_VtxZ->Fill(tempTrackPtr->getPOCA().z());
       Track_HQ_Chi2->Fill(trackChi2);
       Track_HQ_Chi2Red->Fill(trackChi2R);
       Track_HQ_D0->Fill(track_d0);
       Track_HQ_Chi2Red_NStubs->Fill(nStubs, trackChi2R);
       Track_HQ_Chi2Red_Eta->Fill(tempTrackPtr->getMomentum().eta(), trackChi2R);
       Track_HQ_Eta_BarrelStubs->Fill(tempTrackPtr->getMomentum().eta(), nBarrelStubs);
       Track_HQ_Eta_ECStubs->Fill(tempTrackPtr->getMomentum().eta(), nECStubs);
       Track_HQ_Chi2_Probability->Fill(ChiSquaredProbability(trackChi2, nStubs));
     }
 
     //LQ: now defined as all tracks (including HQ tracks)
     numLQTracks++;
 
     Track_LQ_Pt->Fill(tempTrackPtr->getMomentum().perp());
     Track_LQ_Eta->Fill(tempTrackPtr->getMomentum().eta());
     Track_LQ_Phi->Fill(tempTrackPtr->getMomentum().phi());
     Track_LQ_VtxZ->Fill(tempTrackPtr->getPOCA().z());
     Track_LQ_Chi2->Fill(trackChi2);
     Track_LQ_Chi2Red->Fill(trackChi2R);
     Track_LQ_D0->Fill(track_d0);
     Track_LQ_Chi2Red_NStubs->Fill(nStubs, trackChi2R);
     Track_LQ_Chi2Red_Eta->Fill(tempTrackPtr->getMomentum().eta(), trackChi2R);
     Track_LQ_Eta_BarrelStubs->Fill(tempTrackPtr->getMomentum().eta(), nBarrelStubs);
     Track_LQ_Eta_ECStubs->Fill(tempTrackPtr->getMomentum().eta(),  nECStubs);
     Track_LQ_Chi2_Probability->Fill(ChiSquaredProbability(trackChi2, nStubs));
 
   } // End of loop over TTTracks
 
   Track_HQ_N->Fill(numHQTracks);
   Track_LQ_N->Fill(numLQTracks);
 } // end of method

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

overridevirtual

Low-quality tracks (All tracks, including HQ tracks)

High-quality tracks (Now defined as >=5 stubs and chi2/dof < 10)

Implements DQMEDAnalyzer.

Definition at line 151 of file OuterTrackerMonitorTTTrack.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), conf_, DEFINE_FWK_MODULE, edm::ParameterSet::getParameter(), MonitorElement::setAxisTitle(), DQMStore::IBooker::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, topFolderName_, Track_Eta_NStubs, Track_HQ_Chi2, Track_HQ_Chi2_Probability, Track_HQ_Chi2Red, Track_HQ_Chi2Red_Eta, Track_HQ_Chi2Red_NStubs, Track_HQ_D0, Track_HQ_Eta, Track_HQ_Eta_BarrelStubs, Track_HQ_Eta_ECStubs, Track_HQ_N, Track_HQ_Phi, Track_HQ_Pt, Track_HQ_VtxZ, Track_LQ_Chi2, Track_LQ_Chi2_Probability, Track_LQ_Chi2Red, Track_LQ_Chi2Red_Eta, Track_LQ_Chi2Red_NStubs, Track_LQ_D0, Track_LQ_Eta, Track_LQ_Eta_BarrelStubs, Track_LQ_Eta_ECStubs, Track_LQ_N, Track_LQ_Phi, Track_LQ_Pt, Track_LQ_VtxZ, and Track_NStubs.

                                                                                                                       {
   std::string HistoName;
   iBooker.setCurrentFolder(topFolderName_+"/Tracks/");
 
   //Number of stubs
   edm::ParameterSet psTrack_NStubs =  conf_.getParameter<edm::ParameterSet>("TH1_NStubs");
   HistoName = "Track_NStubs";
   Track_NStubs = iBooker.book1D(HistoName, HistoName,
       psTrack_NStubs.getParameter<int32_t>("Nbinsx"),
       psTrack_NStubs.getParameter<double>("xmin"),
       psTrack_NStubs.getParameter<double>("xmax"));
   Track_NStubs->setAxisTitle("# L1 Stubs per L1 Track", 1);
   Track_NStubs->setAxisTitle("# L1 Tracks", 2);
 
   edm::ParameterSet psTrack_Eta_NStubs =  conf_.getParameter<edm::ParameterSet>("TH2_Track_Eta_NStubs");
   HistoName = "Track_Eta_NStubs";
   Track_Eta_NStubs = iBooker.book2D(HistoName, HistoName,
       psTrack_Eta_NStubs.getParameter<int32_t>("Nbinsx"),
       psTrack_Eta_NStubs.getParameter<double>("xmin"),
       psTrack_Eta_NStubs.getParameter<double>("xmax"),
       psTrack_Eta_NStubs.getParameter<int32_t>("Nbinsy"),
       psTrack_Eta_NStubs.getParameter<double>("ymin"),
       psTrack_Eta_NStubs.getParameter<double>("ymax"));
   Track_Eta_NStubs->setAxisTitle("#eta", 1);
   Track_Eta_NStubs->setAxisTitle("# L1 Stubs", 2);
 
   iBooker.setCurrentFolder(topFolderName_+"/Tracks/LQ");
   // Nb of L1Tracks
   HistoName = "Track_LQ_N";
   edm::ParameterSet psTrack_N =  conf_.getParameter<edm::ParameterSet>("TH1_NTracks");
   Track_LQ_N = iBooker.book1D(HistoName, HistoName,
       psTrack_N.getParameter<int32_t>("Nbinsx"),
       psTrack_N.getParameter<double>("xmin"),
       psTrack_N.getParameter<double>("xmax"));
   Track_LQ_N->setAxisTitle("# L1 Tracks", 1);
   Track_LQ_N->setAxisTitle("# Events", 2);
 
   //Pt of the tracks
   edm::ParameterSet psTrack_Pt =  conf_.getParameter<edm::ParameterSet>("TH1_Track_Pt");
   HistoName = "Track_LQ_Pt";
   Track_LQ_Pt = iBooker.book1D(HistoName, HistoName,
       psTrack_Pt.getParameter<int32_t>("Nbinsx"),
       psTrack_Pt.getParameter<double>("xmin"),
       psTrack_Pt.getParameter<double>("xmax"));
   Track_LQ_Pt->setAxisTitle("p_{T} [GeV]", 1);
   Track_LQ_Pt->setAxisTitle("# L1 Tracks", 2);
 
   //Phi
   edm::ParameterSet psTrack_Phi =  conf_.getParameter<edm::ParameterSet>("TH1_Track_Phi");
   HistoName = "Track_LQ_Phi";
   Track_LQ_Phi = iBooker.book1D(HistoName, HistoName,
       psTrack_Phi.getParameter<int32_t>("Nbinsx"),
       psTrack_Phi.getParameter<double>("xmin"),
       psTrack_Phi.getParameter<double>("xmax"));
   Track_LQ_Phi->setAxisTitle("#phi", 1);
   Track_LQ_Phi->setAxisTitle("# L1 Tracks", 2);
 
   //D0
   edm::ParameterSet psTrack_D0 =  conf_.getParameter<edm::ParameterSet>("TH1_Track_D0");
   HistoName = "Track_LQ_D0";
   Track_LQ_D0 = iBooker.book1D(HistoName, HistoName,
       psTrack_D0.getParameter<int32_t>("Nbinsx"),
       psTrack_D0.getParameter<double>("xmin"),
       psTrack_D0.getParameter<double>("xmax"));
   Track_LQ_D0->setAxisTitle("Track D0", 1);
   Track_LQ_D0->setAxisTitle("# L1 Tracks", 2);
 
   //Eta
   edm::ParameterSet psTrack_Eta =  conf_.getParameter<edm::ParameterSet>("TH1_Track_Eta");
   HistoName = "Track_LQ_Eta";
   Track_LQ_Eta = iBooker.book1D(HistoName, HistoName,
       psTrack_Eta.getParameter<int32_t>("Nbinsx"),
       psTrack_Eta.getParameter<double>("xmin"),
       psTrack_Eta.getParameter<double>("xmax"));
   Track_LQ_Eta->setAxisTitle("#eta", 1);
   Track_LQ_Eta->setAxisTitle("# L1 Tracks", 2);
 
   //VtxZ
   edm::ParameterSet psTrack_VtxZ =  conf_.getParameter<edm::ParameterSet>("TH1_Track_VtxZ");
   HistoName = "Track_LQ_VtxZ";
   Track_LQ_VtxZ = iBooker.book1D(HistoName, HistoName,
       psTrack_VtxZ.getParameter<int32_t>("Nbinsx"),
       psTrack_VtxZ.getParameter<double>("xmin"),
       psTrack_VtxZ.getParameter<double>("xmax"));
   Track_LQ_VtxZ->setAxisTitle("L1 Track vertex position z [cm]", 1);
   Track_LQ_VtxZ->setAxisTitle("# L1 Tracks", 2);
 
   //chi2
   edm::ParameterSet psTrack_Chi2 =  conf_.getParameter<edm::ParameterSet>("TH1_Track_Chi2");
   HistoName = "Track_LQ_Chi2";
   Track_LQ_Chi2 = iBooker.book1D(HistoName, HistoName,
       psTrack_Chi2.getParameter<int32_t>("Nbinsx"),
       psTrack_Chi2.getParameter<double>("xmin"),
       psTrack_Chi2.getParameter<double>("xmax"));
   Track_LQ_Chi2->setAxisTitle("L1 Track #chi^{2}", 1);
   Track_LQ_Chi2->setAxisTitle("# L1 Tracks", 2);
 
   //chi2Red
   edm::ParameterSet psTrack_Chi2Red =  conf_.getParameter<edm::ParameterSet>("TH1_Track_Chi2R");
   HistoName = "Track_LQ_Chi2Red";
   Track_LQ_Chi2Red = iBooker.book1D(HistoName, HistoName,
       psTrack_Chi2Red.getParameter<int32_t>("Nbinsx"),
       psTrack_Chi2Red.getParameter<double>("xmin"),
       psTrack_Chi2Red.getParameter<double>("xmax"));
   Track_LQ_Chi2Red->setAxisTitle("L1 Track #chi^{2}/ndf", 1);
   Track_LQ_Chi2Red->setAxisTitle("# L1 Tracks", 2);
 
   //Chi2 prob
   edm::ParameterSet psTrack_Chi2_Probability =  conf_.getParameter<edm::ParameterSet>("TH1_Track_Chi2_Probability");
   HistoName = "Track_LQ_Chi2_Probability";
   Track_LQ_Chi2_Probability = iBooker.book1D(HistoName, HistoName,
       psTrack_Chi2_Probability.getParameter<int32_t>("Nbinsx"),
       psTrack_Chi2_Probability.getParameter<double>("xmin"),
       psTrack_Chi2_Probability.getParameter<double>("xmax"));
   Track_LQ_Chi2_Probability->setAxisTitle("#chi^{2} probability", 1);
   Track_LQ_Chi2_Probability->setAxisTitle("# L1 Tracks", 2);
 
   //Reduced chi2 vs #stubs
   edm::ParameterSet psTrack_Chi2Red_NStubs =  conf_.getParameter<edm::ParameterSet>("TH2_Track_Chi2R_NStubs");
   HistoName = "Track_LQ_Chi2Red_NStubs";
   Track_LQ_Chi2Red_NStubs = iBooker.book2D(HistoName, HistoName,
       psTrack_Chi2Red_NStubs.getParameter<int32_t>("Nbinsx"),
       psTrack_Chi2Red_NStubs.getParameter<double>("xmin"),
       psTrack_Chi2Red_NStubs.getParameter<double>("xmax"),
       psTrack_Chi2Red_NStubs.getParameter<int32_t>("Nbinsy"),
       psTrack_Chi2Red_NStubs.getParameter<double>("ymin"),
       psTrack_Chi2Red_NStubs.getParameter<double>("ymax"));
   Track_LQ_Chi2Red_NStubs->setAxisTitle("# L1 Stubs", 1);
   Track_LQ_Chi2Red_NStubs->setAxisTitle("L1 Track #chi^{2}/ndf", 2);
 
   //chi2/dof vs eta
   edm::ParameterSet psTrack_Chi2R_Eta =  conf_.getParameter<edm::ParameterSet>("TH2_Track_Chi2R_Eta");
   HistoName = "Track_LQ_Chi2Red_Eta";
   Track_LQ_Chi2Red_Eta = iBooker.book2D(HistoName, HistoName,
       psTrack_Chi2R_Eta.getParameter<int32_t>("Nbinsx"),
       psTrack_Chi2R_Eta.getParameter<double>("xmin"),
       psTrack_Chi2R_Eta.getParameter<double>("xmax"),
       psTrack_Chi2R_Eta.getParameter<int32_t>("Nbinsy"),
       psTrack_Chi2R_Eta.getParameter<double>("ymin"),
       psTrack_Chi2R_Eta.getParameter<double>("ymax"));
   Track_LQ_Chi2Red_Eta->setAxisTitle("#eta", 1);
   Track_LQ_Chi2Red_Eta->setAxisTitle("L1 Track #chi^{2}/ndf", 2);
 
   //Eta vs #stubs in barrel
   HistoName = "Track_LQ_Eta_BarrelStubs";
   Track_LQ_Eta_BarrelStubs = iBooker.book2D(HistoName, HistoName,
       psTrack_Eta_NStubs.getParameter<int32_t>("Nbinsx"),
       psTrack_Eta_NStubs.getParameter<double>("xmin"),
       psTrack_Eta_NStubs.getParameter<double>("xmax"),
       psTrack_Eta_NStubs.getParameter<int32_t>("Nbinsy"),
       psTrack_Eta_NStubs.getParameter<double>("ymin"),
       psTrack_Eta_NStubs.getParameter<double>("ymax"));
   Track_LQ_Eta_BarrelStubs->setAxisTitle("#eta", 1);
   Track_LQ_Eta_BarrelStubs->setAxisTitle("# L1 Barrel Stubs", 2);
 
   //Eta vs #stubs in EC
   HistoName = "Track_LQ_Eta_ECStubs";
   Track_LQ_Eta_ECStubs = iBooker.book2D(HistoName, HistoName,
       psTrack_Eta_NStubs.getParameter<int32_t>("Nbinsx"),
       psTrack_Eta_NStubs.getParameter<double>("xmin"),
       psTrack_Eta_NStubs.getParameter<double>("xmax"),
       psTrack_Eta_NStubs.getParameter<int32_t>("Nbinsy"),
       psTrack_Eta_NStubs.getParameter<double>("ymin"),
       psTrack_Eta_NStubs.getParameter<double>("ymax"));
   Track_LQ_Eta_ECStubs->setAxisTitle("#eta", 1);
   Track_LQ_Eta_ECStubs->setAxisTitle("# L1 EC Stubs", 2);
 
   iBooker.setCurrentFolder(topFolderName_+"/Tracks/HQ");
   // Nb of L1Tracks
   HistoName = "Track_HQ_N";
   Track_HQ_N = iBooker.book1D(HistoName, HistoName,
       psTrack_N.getParameter<int32_t>("Nbinsx"),
       psTrack_N.getParameter<double>("xmin"),
       psTrack_N.getParameter<double>("xmax"));
   Track_HQ_N->setAxisTitle("# L1 Tracks", 1);
   Track_HQ_N->setAxisTitle("# Events", 2);
 
   //Pt of the tracks
   HistoName = "Track_HQ_Pt";
   Track_HQ_Pt = iBooker.book1D(HistoName, HistoName,
       psTrack_Pt.getParameter<int32_t>("Nbinsx"),
       psTrack_Pt.getParameter<double>("xmin"),
       psTrack_Pt.getParameter<double>("xmax"));
   Track_HQ_Pt->setAxisTitle("p_{T} [GeV]", 1);
   Track_HQ_Pt->setAxisTitle("# L1 Tracks", 2);
 
   //Phi
   HistoName = "Track_HQ_Phi";
   Track_HQ_Phi = iBooker.book1D(HistoName, HistoName,
       psTrack_Phi.getParameter<int32_t>("Nbinsx"),
       psTrack_Phi.getParameter<double>("xmin"),
       psTrack_Phi.getParameter<double>("xmax"));
   Track_HQ_Phi->setAxisTitle("#phi", 1);
   Track_HQ_Phi->setAxisTitle("# L1 Tracks", 2);
 
   //D0
   HistoName = "Track_HQ_D0";
   Track_HQ_D0 = iBooker.book1D(HistoName, HistoName,
       psTrack_D0.getParameter<int32_t>("Nbinsx"),
       psTrack_D0.getParameter<double>("xmin"),
       psTrack_D0.getParameter<double>("xmax"));
   Track_HQ_D0->setAxisTitle("Track D0", 1);
   Track_HQ_D0->setAxisTitle("# L1 Tracks", 2);
 
   //Eta
   HistoName = "Track_HQ_Eta";
   Track_HQ_Eta = iBooker.book1D(HistoName, HistoName,
       psTrack_Eta.getParameter<int32_t>("Nbinsx"),
       psTrack_Eta.getParameter<double>("xmin"),
       psTrack_Eta.getParameter<double>("xmax"));
   Track_HQ_Eta->setAxisTitle("#eta", 1);
   Track_HQ_Eta->setAxisTitle("# L1 Tracks", 2);
 
   //VtxZ
   HistoName = "Track_HQ_VtxZ";
   Track_HQ_VtxZ = iBooker.book1D(HistoName, HistoName,
       psTrack_VtxZ.getParameter<int32_t>("Nbinsx"),
       psTrack_VtxZ.getParameter<double>("xmin"),
       psTrack_VtxZ.getParameter<double>("xmax"));
   Track_HQ_VtxZ->setAxisTitle("L1 Track vertex position z [cm]", 1);
   Track_HQ_VtxZ->setAxisTitle("# L1 Tracks", 2);
 
   //chi2
   HistoName = "Track_HQ_Chi2";
   Track_HQ_Chi2 = iBooker.book1D(HistoName, HistoName,
       psTrack_Chi2.getParameter<int32_t>("Nbinsx"),
       psTrack_Chi2.getParameter<double>("xmin"),
       psTrack_Chi2.getParameter<double>("xmax"));
   Track_HQ_Chi2->setAxisTitle("L1 Track #chi^{2}", 1);
   Track_HQ_Chi2->setAxisTitle("# L1 Tracks", 2);
 
   //chi2Red
   HistoName = "Track_HQ_Chi2Red";
   Track_HQ_Chi2Red = iBooker.book1D(HistoName, HistoName,
       psTrack_Chi2Red.getParameter<int32_t>("Nbinsx"),
       psTrack_Chi2Red.getParameter<double>("xmin"),
       psTrack_Chi2Red.getParameter<double>("xmax"));
   Track_HQ_Chi2Red->setAxisTitle("L1 Track #chi^{2}/ndf", 1);
   Track_HQ_Chi2Red->setAxisTitle("# L1 Tracks", 2);
 
   //Chi2 prob
   HistoName = "Track_HQ_Chi2_Probability";
   Track_HQ_Chi2_Probability = iBooker.book1D(HistoName, HistoName,
       psTrack_Chi2_Probability.getParameter<int32_t>("Nbinsx"),
       psTrack_Chi2_Probability.getParameter<double>("xmin"),
       psTrack_Chi2_Probability.getParameter<double>("xmax"));
   Track_HQ_Chi2_Probability->setAxisTitle("#chi^{2} probability", 1);
   Track_HQ_Chi2_Probability->setAxisTitle("# L1 Tracks", 2);
 
   //Reduced chi2 vs #stubs
   HistoName = "Track_HQ_Chi2Red_NStubs";
   Track_HQ_Chi2Red_NStubs = iBooker.book2D(HistoName, HistoName,
       psTrack_Chi2Red_NStubs.getParameter<int32_t>("Nbinsx"),
       psTrack_Chi2Red_NStubs.getParameter<double>("xmin"),
       psTrack_Chi2Red_NStubs.getParameter<double>("xmax"),
       psTrack_Chi2Red_NStubs.getParameter<int32_t>("Nbinsy"),
       psTrack_Chi2Red_NStubs.getParameter<double>("ymin"),
       psTrack_Chi2Red_NStubs.getParameter<double>("ymax"));
   Track_HQ_Chi2Red_NStubs->setAxisTitle("# L1 Stubs", 1);
   Track_HQ_Chi2Red_NStubs->setAxisTitle("L1 Track #chi^{2}/ndf", 2);
 
   //chi2/dof vs eta
   HistoName = "Track_HQ_Chi2Red_Eta";
   Track_HQ_Chi2Red_Eta = iBooker.book2D(HistoName, HistoName,
       psTrack_Chi2R_Eta.getParameter<int32_t>("Nbinsx"),
       psTrack_Chi2R_Eta.getParameter<double>("xmin"),
       psTrack_Chi2R_Eta.getParameter<double>("xmax"),
       psTrack_Chi2R_Eta.getParameter<int32_t>("Nbinsy"),
       psTrack_Chi2R_Eta.getParameter<double>("ymin"),
       psTrack_Chi2R_Eta.getParameter<double>("ymax"));
   Track_HQ_Chi2Red_Eta->setAxisTitle("#eta", 1);
   Track_HQ_Chi2Red_Eta->setAxisTitle("L1 Track #chi^{2}/ndf", 2);
 
   //eta vs #stubs in barrel
   HistoName = "Track_HQ_Eta_BarrelStubs";
   Track_HQ_Eta_BarrelStubs = iBooker.book2D(HistoName, HistoName,
       psTrack_Eta_NStubs.getParameter<int32_t>("Nbinsx"),
       psTrack_Eta_NStubs.getParameter<double>("xmin"),
       psTrack_Eta_NStubs.getParameter<double>("xmax"),
       psTrack_Eta_NStubs.getParameter<int32_t>("Nbinsy"),
       psTrack_Eta_NStubs.getParameter<double>("ymin"),
       psTrack_Eta_NStubs.getParameter<double>("ymax"));
   Track_HQ_Eta_BarrelStubs->setAxisTitle("#eta", 1);
   Track_HQ_Eta_BarrelStubs->setAxisTitle("# L1 Barrel Stubs", 2);
 
   //eta vs #stubs in EC
   HistoName = "Track_HQ_Eta_ECStubs";
   Track_HQ_Eta_ECStubs = iBooker.book2D(HistoName, HistoName,
       psTrack_Eta_NStubs.getParameter<int32_t>("Nbinsx"),
       psTrack_Eta_NStubs.getParameter<double>("xmin"),
       psTrack_Eta_NStubs.getParameter<double>("xmax"),
       psTrack_Eta_NStubs.getParameter<int32_t>("Nbinsy"),
       psTrack_Eta_NStubs.getParameter<double>("ymin"),
       psTrack_Eta_NStubs.getParameter<double>("ymax"));
   Track_HQ_Eta_ECStubs->setAxisTitle("#eta", 1);
   Track_HQ_Eta_ECStubs->setAxisTitle("# L1 EC Stubs", 2);
 
 }//end of method