#include <DTLocalTriggerEfficiencyTest.h>

Inheritance diagram for DTLocalTriggerEfficiencyTest:

Public Member Functions
	DTLocalTriggerEfficiencyTest (const edm::ParameterSet &ps)
	Constructor. More...

	~DTLocalTriggerEfficiencyTest () override
	Destructor. More...

Public Member Functions inherited from DTLocalTriggerBaseTest
	DTLocalTriggerBaseTest ()
	Constructor. More...

	~DTLocalTriggerBaseTest () override
	Destructor. More...

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

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

void	beginRun (edm::Run const &, edm::EventSetup const &) override

	DQMEDHarvester ()

void	endJob () final

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

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

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

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

	~DQMEDHarvester () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns, edm::one::SharedResources >
	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 () 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)

Protected Member Functions
void	beginRun (const edm::Run &r, const edm::EventSetup &c) override
	BeginRun. More...

void	bookChambHistos (DQMStore::IBooker &, DTChamberId chambId, std::string htype)
	Book the new MEs (for each chamber) More...

void	Bookings (DQMStore::IBooker &, DQMStore::IGetter &)

void	makeEfficiencyME (TH1D numerator, TH1D denominator, MonitorElement *result)
	Compute efficiency plots. More...

void	makeEfficiencyME2D (TH2F numerator, TH2F denominator, MonitorElement *result)
	Compute 2D efficiency plots. More...

void	runClientDiagnostic (DQMStore::IBooker &, DQMStore::IGetter &) override
	DQM Client Diagnostic. More...

Protected Member Functions inherited from DTLocalTriggerBaseTest
void	beginRun (edm::Run const &run, edm::EventSetup const &context) override
	BeginRun. More...

void	bookCmsHistos (DQMStore::IBooker &, std::string hTag, std::string folder="", bool isGlb=false)
	Book the new MEs (CMS summary) More...

void	bookSectorHistos (DQMStore::IBooker &, int wheel, int sector, std::string hTag, std::string folder="")
	Book the new MEs (for each sector) More...

void	bookWheelHistos (DQMStore::IBooker &, int wheel, std::string hTag, std::string folder="")
	Book the new MEs (for each wheel) More...

std::string	category ()
	Get message logger name. More...

void	dqmEndJob (DQMStore::IBooker &, DQMStore::IGetter &) override

void	dqmEndLuminosityBlock (DQMStore::IBooker &, DQMStore::IGetter &, edm::LuminosityBlock const &, edm::EventSetup const &) override
	Perform client diagnostic in online. More...

void	endRun (edm::Run const &run, edm::EventSetup const &context) override
	Perform client diagnostic in offline. More...

std::string	fullName (std::string htype)
	Create fullname from histo partial name. More...

template<class T >
T *	getHisto (MonitorElement *me)
	Convert ME to Histogram fo type T. More...

std::string	getMEName (std::string histoTag, std::string subfolder, const DTChamberId &chambid)
	Get the ME name (by chamber) More...

std::string	getMEName (std::string histoTag, std::string subfolder, int wh)
	Get the ME name (by wheel) More...

std::pair< float, float >	phiRange (const DTChamberId &id)
	Calculate phi range for histograms. More...

void	setConfig (const edm::ParameterSet &ps, std::string name)
	Set configuration variables. More...

std::string &	topFolder (bool isTM)
	Get top folder name. More...

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
const int	wheelArrayShift = 3

Protected Attributes inherited from DTLocalTriggerBaseTest
std::string	baseFolderDDU

std::string	baseFolderTM

std::map< std::string, MonitorElement * >	cmsME

std::string	hwSource

std::vector< std::string >	hwSources

edm::ESHandle< DTGeometry >	muonGeom

int	nevents

unsigned int	nLumiSegs

edm::ParameterSet	parameters

int	prescaleFactor

int	run

bool	runOnline

std::map< int, std::map< std::string, MonitorElement * > >	secME

std::string	sourceFolder

std::string	testName

std::string	trigSource

std::vector< std::string >	trigSources

std::map< int, std::map< std::string, MonitorElement * > >	whME

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

edm::EDPutTokenT< DQMToken >	runToken_

Private Attributes
bool	bookingdone

std::map< uint32_t, std::map< std::string, MonitorElement * > >	chambME

DTTrigGeomUtils *	trigGeomUtils

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)

Detailed Description

DQM Test Client

Author
C. Battilana S. Marcellini - INFN Bologna

threadsafe version (//-) oct/nov 2014 - WATWanAbdullah ncpp-um-my

Definition at line 22 of file DTLocalTriggerEfficiencyTest.h.

Constructor & Destructor Documentation

DTLocalTriggerEfficiencyTest::DTLocalTriggerEfficiencyTest ( const edm::ParameterSet & ps )

Constructor.

Definition at line 38 of file DTLocalTriggerEfficiencyTest.cc.

References DTLocalTriggerBaseTest::baseFolderDDU, DTLocalTriggerBaseTest::baseFolderTM, bookingdone, and DTLocalTriggerBaseTest::setConfig().

                                                                                     : trigGeomUtils(nullptr) {
 
   setConfig(ps,"DTLocalTriggerEfficiency");
   baseFolderTM = "DT/03-LocalTrigger-TM/";
   baseFolderDDU = "DT/04-LocalTrigger-DDU/";
 
   bookingdone = false;
 
 }

DTLocalTriggerEfficiencyTest::~DTLocalTriggerEfficiencyTest ( )

override

Destructor.

Definition at line 49 of file DTLocalTriggerEfficiencyTest.cc.

References trigGeomUtils.

                                                            {
 
   if (trigGeomUtils) {
     delete trigGeomUtils;
   }
   
 }

Member Function Documentation

void DTLocalTriggerEfficiencyTest::beginRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)

overrideprotected

BeginRun.

Definition at line 58 of file DTLocalTriggerEfficiencyTest.cc.

References DTLocalTriggerBaseTest::beginRun(), DTLocalTriggerBaseTest::muonGeom, and trigGeomUtils.

                                                                                   {
 
   DTLocalTriggerBaseTest::beginRun(r,c);
   trigGeomUtils = new DTTrigGeomUtils(muonGeom);
 }

void DTLocalTriggerEfficiencyTest::bookChambHistos	(	DQMStore::IBooker &	,
		DTChamberId	chambId,
		std::string	htype
	)

protected

Book the new MEs (for each chamber)

Definition at line 312 of file DTLocalTriggerEfficiencyTest.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DTLocalTriggerBaseTest::category(), chambME, DTLocalTriggerBaseTest::fullName(), DTLocalTriggerBaseTest::hwSource, LogTrace, SiStripPI::max, min(), pileupCalc::nbins, DTTrigGeomUtils::phiRange(), DetId::rawId(), DTChamberId::sector(), DQMStore::IBooker::setCurrentFolder(), DTChamberId::station(), relativeConstraints::station, DTLocalTriggerBaseTest::testName, DTTrigGeomUtils::thetaRange(), DTLocalTriggerBaseTest::topFolder(), trigGeomUtils, DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

Referenced by Bookings(), and runClientDiagnostic().

                                                                                                               {
   
   stringstream wheel; wheel << chambId.wheel();
   stringstream station; station << chambId.station();   
   stringstream sector; sector << chambId.sector();
 
   string fullType  = fullName(htype);
   bool isTM = hwSource=="TM" ;
   string HistoName = fullType + "_W" + wheel.str() + "_Sec" + sector.str() + "_St" + station.str();
 
   ibooker.setCurrentFolder(topFolder(isTM) + "Wheel" + wheel.str() +
             "/Sector" + sector.str() +
             "/Station" + station.str() + "/Segment");
 
   LogTrace(category()) << "[" << testName << "Test]: booking " + topFolder(isTM) + "Wheel" << wheel.str() 
                <<"/Sector" << sector.str() << "/Station" << station.str() << "/Segment/" << HistoName;
 
   
   uint32_t indexChId = chambId.rawId();
   if (htype.find("TrigEffAnglePhi") == 0){
     chambME[indexChId][fullType] = ibooker.book1D(HistoName.c_str(),"Trigger efficiency vs angle of incidence (Phi)",16,-40.,40.);
   }
   else if (htype.find("TrigEffAngleHHHLPhi") == 0){
     chambME[indexChId][fullType] = ibooker.book1D(HistoName.c_str(),"Trigger efficiency (HH/HL) vs angle of incidence (Phi)",16,-40.,40.);
   }
   else if (htype.find("TrigEffAngleTheta") == 0){
     chambME[indexChId][fullType] = ibooker.book1D(HistoName.c_str(),"Trigger efficiency vs angle of incidence (Theta)",16,-40.,40.);
   }
   else if (htype.find("TrigEffAngleHTheta") == 0){
     chambME[indexChId][fullType] = ibooker.book1D(HistoName.c_str(),"Trigger efficiency (H) vs angle of incidence (Theta)",16,-40.,40.);
   }
   else if (htype.find("TrigEffPosPhi") == 0 ){
     float min,max;
     int nbins;
     trigGeomUtils->phiRange(chambId,min,max,nbins);
     chambME[indexChId][fullType] = ibooker.book1D(HistoName.c_str(),"Trigger efficiency vs position (Phi)",nbins,min,max);
   }
   else if (htype.find("TrigEffPosvsAnglePhi") == 0 ){
     float min,max;
     int nbins;
     trigGeomUtils->phiRange(chambId,min,max,nbins);
     chambME[indexChId][fullType] = ibooker.book2D(HistoName.c_str(),"Trigger efficiency position vs angle (Phi)",16,-40.,40.,nbins,min,max);
   }
   else if (htype.find("TrigEffPosvsAngleHHHLPhi") == 0 ){
     float min,max;
     int nbins;
     trigGeomUtils->phiRange(chambId,min,max,nbins);
     chambME[indexChId][fullType] = ibooker.book2D(HistoName.c_str(),"Trigger efficiency (HH/HL) pos vs angle (Phi)",16,-40.,40.,nbins,min,max);
   }
   else if (htype.find("TrigEffPosHHHLPhi") == 0 ){
     float min,max;
     int nbins;
     trigGeomUtils->phiRange(chambId,min,max,nbins);
     chambME[indexChId][fullType] = ibooker.book1D(HistoName.c_str(),"Trigger efficiency (HH/HL) vs position (Phi)",nbins,min,max);
   }
   else if (htype.find("TrigEffPosTheta") == 0){
     float min,max;
     int nbins;
     trigGeomUtils->thetaRange(chambId,min,max,nbins);
     chambME[indexChId][fullType] = ibooker.book1D(HistoName.c_str(),"Trigger efficiency vs position (Theta)",nbins,min,max);
   }
   else if (htype.find("TrigEffPosHTheta") == 0){
     float min,max;
     int nbins;
     trigGeomUtils->thetaRange(chambId,min,max,nbins);
     chambME[indexChId][fullType] = ibooker.book1D(HistoName.c_str(),"Trigger efficiency (H) vs position (Theta)",nbins,min,max);
   }
   else if (htype.find("TrigEffPosvsAngleTheta") == 0 ){
     float min,max;
     int nbins;
     trigGeomUtils->thetaRange(chambId,min,max,nbins);
     chambME[indexChId][fullType] = ibooker.book2D(HistoName.c_str(),"Trigger efficiency pos vs angle (Theta)",16,-40.,40.,nbins,min,max);
   }
   else if (htype.find("TrigEffPosvsAngleHTheta") == 0 ){
     float min,max;
     int nbins;
     trigGeomUtils->thetaRange(chambId,min,max,nbins);
     chambME[indexChId][fullType] = ibooker.book2D(HistoName.c_str(),"Trigger efficiency (H) pos vs angle (Theta)",16,-40.,40.,nbins,min,max);
   }
 
 }

void DTLocalTriggerEfficiencyTest::Bookings	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter
	)

protected

Definition at line 66 of file DTLocalTriggerEfficiencyTest.cc.

References bookChambHistos(), bookingdone, DTLocalTriggerBaseTest::bookSectorHistos(), DTLocalTriggerBaseTest::bookWheelHistos(), DTLocalTriggerBaseTest::hwSource, DTLocalTriggerBaseTest::hwSources, trackingPlots::stat, DTLocalTriggerBaseTest::trigSource, and DTLocalTriggerBaseTest::trigSources.

Referenced by runClientDiagnostic().

                                                                                                 {
 
   vector<string>::const_iterator iTr   = trigSources.begin();
   vector<string>::const_iterator trEnd = trigSources.end();
   vector<string>::const_iterator iHw   = hwSources.begin();
   vector<string>::const_iterator hwEnd = hwSources.end();
 
 
   //Booking
   if(parameters.getUntrackedParameter<bool>("staticBooking", true)){
     for (; iTr != trEnd; ++iTr){
       trigSource = (*iTr);
       for (; iHw != hwEnd; ++iHw){
     hwSource = (*iHw);
     // Loop over the TriggerUnits
     for (int wh=-2; wh<=2; ++wh){
       for (int sect=1; sect<=12; ++sect){
         for (int stat=1; stat<=4; ++stat){
           DTChamberId chId(wh,stat,sect);
           bookChambHistos(ibooker,chId,"TrigEffPosvsAnglePhi");
           bookChambHistos(ibooker,chId,"TrigEffPosvsAngleHHHLPhi");
           bookChambHistos(ibooker,chId,"TrigEffPosPhi");
           bookChambHistos(ibooker,chId,"TrigEffPosHHHLPhi");
           bookChambHistos(ibooker,chId,"TrigEffAnglePhi");
           bookChambHistos(ibooker,chId,"TrigEffAngleHHHLPhi");
           if (stat<=3) {
         bookChambHistos(ibooker,chId,"TrigEffPosvsAngleTheta");
         bookChambHistos(ibooker,chId,"TrigEffPosvsAngleHTheta");
         bookChambHistos(ibooker,chId,"TrigEffPosTheta");
         bookChambHistos(ibooker,chId,"TrigEffPosHTheta");
         bookChambHistos(ibooker,chId,"TrigEffAngleTheta");
         bookChambHistos(ibooker,chId,"TrigEffAngleHTheta");
           }
         }
 
         bookSectorHistos(ibooker,wh,sect,"TrigEffPhi");  
         bookSectorHistos(ibooker,wh,sect,"TrigEffTheta");  
       }
  
       bookWheelHistos(ibooker,wh,"TrigEffPhi");  
       bookWheelHistos(ibooker,wh,"TrigEffHHHLPhi");  
       bookWheelHistos(ibooker,wh,"TrigEffTheta");  
       bookWheelHistos(ibooker,wh,"TrigEffHTheta");  
     }
       }
     }
   }
 
   bookingdone = true; 
 }

void DTLocalTriggerEfficiencyTest::makeEfficiencyME	(	TH1D *	numerator,
		TH1D *	denominator,
		MonitorElement *	result
	)

protected

Compute efficiency plots.

Definition at line 262 of file DTLocalTriggerEfficiencyTest.cc.

References stringResolutionProvider_cfi::bin, PostProcessor_cff::efficiency, relativeConstraints::error, MonitorElement::getTH1F(), pileupCalc::nbins, and mathSSE::sqrt().

Referenced by runClientDiagnostic().

                                                                                                              {
   
   TH1F* efficiency = result->getTH1F();
   efficiency->Divide(numerator,denominator,1,1,"");
   
   int nbins = efficiency->GetNbinsX();
   for (int bin=1; bin<=nbins; ++bin){
     float error = 0;
     float bineff = efficiency->GetBinContent(bin);
 
     if (denominator->GetBinContent(bin)){
       error = sqrt(bineff*(1-bineff)/denominator->GetBinContent(bin));
     }
     else {
       error = 1;
       efficiency->SetBinContent(bin,1.);
     }
  
     efficiency->SetBinError(bin,error);
   }
 
 }

void DTLocalTriggerEfficiencyTest::makeEfficiencyME2D	(	TH2F *	numerator,
		TH2F *	denominator,
		MonitorElement *	result
	)

protected

Compute 2D efficiency plots.

Definition at line 286 of file DTLocalTriggerEfficiencyTest.cc.

References PostProcessor_cff::efficiency, relativeConstraints::error, MonitorElement::getTH2F(), and mathSSE::sqrt().

Referenced by runClientDiagnostic().

                                                                                                                {
   
   TH2F* efficiency = result->getTH2F();
   efficiency->Divide(numerator,denominator,1,1,"");
   
   int nbinsx = efficiency->GetNbinsX();
   int nbinsy = efficiency->GetNbinsY();
   for (int binx=1; binx<=nbinsx; ++binx){
     for (int biny=1; biny<=nbinsy; ++biny){
       float error = 0;
       float bineff = efficiency->GetBinContent(binx,biny);
 
       if (denominator->GetBinContent(binx,biny)){
     error = sqrt(bineff*(1-bineff)/denominator->GetBinContent(binx,biny));
       }
       else {
     error = 1;
     efficiency->SetBinContent(binx,biny,0.);
       }
  
       efficiency->SetBinError(binx,biny,error);
     }
   }
 
 }    

void DTLocalTriggerEfficiencyTest::runClientDiagnostic	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter
	)

overrideprotectedvirtual

DQM Client Diagnostic.

Implements DTLocalTriggerBaseTest.

Definition at line 117 of file DTLocalTriggerEfficiencyTest.cc.

References bookChambHistos(), bookingdone, Bookings(), DTLocalTriggerBaseTest::bookSectorHistos(), DTLocalTriggerBaseTest::bookWheelHistos(), chambME, end, spr::find(), objects.autophobj::float, DTLocalTriggerBaseTest::fullName(), DQMStore::IGetter::get(), DTLocalTriggerBaseTest::getMEName(), DTLocalTriggerBaseTest::hwSource, DTLocalTriggerBaseTest::hwSources, makeEfficiencyME(), makeEfficiencyME2D(), DetId::rawId(), DTLocalTriggerBaseTest::secME, MonitorElement::setBinContent(), mathSSE::sqrt(), trackingPlots::stat, reco::tau::disc::TrackAngle(), DTLocalTriggerBaseTest::trigSource, DTLocalTriggerBaseTest::trigSources, wheelArrayShift, and DTLocalTriggerBaseTest::whME.

                                                                                              {
  if (!bookingdone) Bookings(ibooker,igetter);
 
   // Loop over Trig & Hw sources
   for (vector<string>::const_iterator iTr = trigSources.begin(); iTr != trigSources.end(); ++iTr){
     trigSource = (*iTr);
     for (vector<string>::const_iterator iHw = hwSources.begin(); iHw != hwSources.end(); ++iHw){
       hwSource = (*iHw);
       // Loop over the TriggerUnits
       for (int stat=1; stat<=4; ++stat){
     for (int wh=-2; wh<=2; ++wh){
       for (int sect=1; sect<=12; ++sect){
         DTChamberId chId(wh,stat,sect);
         int sector_id = (wh+wheelArrayShift)+(sect-1)*5;
         uint32_t indexCh = chId.rawId();
 
 
         TH2F * TrackPosvsAngle            = getHisto<TH2F>(igetter.get(getMEName("TrackPosvsAngle","Segment", chId)));
         TH2F * TrackPosvsAngleandTrig     = getHisto<TH2F>(igetter.get(getMEName("TrackPosvsAngleandTrig","Segment", chId)));
         TH2F * TrackPosvsAngleandTrigHHHL = getHisto<TH2F>(igetter.get(getMEName("TrackPosvsAngleandTrigHHHL","Segment", chId)));
         
         if (TrackPosvsAngle && TrackPosvsAngleandTrig && TrackPosvsAngleandTrigHHHL && TrackPosvsAngle->GetEntries()>1) {
           
           if( chambME[indexCh].find(fullName("TrigEffAnglePhi")) == chambME[indexCh].end()){
 
         bookChambHistos(ibooker,chId,"TrigEffPosvsAnglePhi");
         bookChambHistos(ibooker,chId,"TrigEffPosvsAngleHHHLPhi");
         bookChambHistos(ibooker,chId,"TrigEffPosPhi");
         bookChambHistos(ibooker,chId,"TrigEffPosHHHLPhi");
         bookChambHistos(ibooker,chId,"TrigEffAnglePhi");
         bookChambHistos(ibooker,chId,"TrigEffAngleHHHLPhi");
           }
           if( secME[sector_id].find(fullName("TrigEffPhi")) == secME[sector_id].end() ){
 
         bookSectorHistos(ibooker,wh,sect,"TrigEffPhi");  
           }
           if( whME[wh].find(fullName("TrigEffPhi")) == whME[wh].end() ){
 
         bookWheelHistos(ibooker,wh,"TrigEffPhi");  
         bookWheelHistos(ibooker,wh,"TrigEffHHHLPhi");  
           }
 
           std::map<std::string,MonitorElement*> *innerME = &(secME[sector_id]);
           TH1D* TrackPos               = TrackPosvsAngle->ProjectionY();
           TH1D* TrackAngle             = TrackPosvsAngle->ProjectionX();
           TH1D* TrackPosandTrig        = TrackPosvsAngleandTrig->ProjectionY();
           TH1D* TrackAngleandTrig      = TrackPosvsAngleandTrig->ProjectionX();
           TH1D* TrackPosandTrigHHHL    = TrackPosvsAngleandTrigHHHL->ProjectionY();
           TH1D* TrackAngleandTrigHHHL  = TrackPosvsAngleandTrigHHHL->ProjectionX();
           float binEff     = float(TrackPosandTrig->GetEntries())/TrackPos->GetEntries();
           float binEffHHHL = float(TrackPosandTrigHHHL->GetEntries())/TrackPos->GetEntries();
           float binErr     = sqrt(binEff*(1-binEff)/TrackPos->GetEntries());
           float binErrHHHL = sqrt(binEffHHHL*(1-binEffHHHL)/TrackPos->GetEntries());
       
           MonitorElement* globalEff = innerME->find(fullName("TrigEffPhi"))->second;
           globalEff->setBinContent(stat,binEff);
           globalEff->setBinError(stat,binErr);
 
           innerME = &(whME[wh]);
           globalEff = innerME->find(fullName("TrigEffPhi"))->second;
           globalEff->setBinContent(sect,stat,binEff);
           globalEff->setBinError(sect,stat,binErr);
           globalEff = innerME->find(fullName("TrigEffHHHLPhi"))->second;
           globalEff->setBinContent(sect,stat,binEffHHHL);
           globalEff->setBinError(sect,stat,binErrHHHL);
       
       
           innerME = &(chambME[indexCh]);
           makeEfficiencyME(TrackPosandTrig,TrackPos,innerME->find(fullName("TrigEffPosPhi"))->second);
           makeEfficiencyME(TrackPosandTrigHHHL,TrackPos,innerME->find(fullName("TrigEffPosHHHLPhi"))->second);
           makeEfficiencyME(TrackAngleandTrig,TrackAngle,innerME->find(fullName("TrigEffAnglePhi"))->second);
           makeEfficiencyME(TrackAngleandTrigHHHL,TrackAngle,innerME->find(fullName("TrigEffAngleHHHLPhi"))->second);
           makeEfficiencyME2D(TrackPosvsAngleandTrig,TrackPosvsAngle,innerME->find(fullName("TrigEffPosvsAnglePhi"))->second);
           makeEfficiencyME2D(TrackPosvsAngleandTrigHHHL,TrackPosvsAngle,innerME->find(fullName("TrigEffPosvsAngleHHHLPhi"))->second);
          
         }
     
         // Perform Efficiency analysis (Theta+Segments)  CB FIXME -> no TM theta qual info
         TH2F * TrackThetaPosvsAngle            = getHisto<TH2F>(igetter.get(getMEName("TrackThetaPosvsAngle","Segment", chId)));
         TH2F * TrackThetaPosvsAngleandTrig     = getHisto<TH2F>(igetter.get(getMEName("TrackThetaPosvsAngleandTrig","Segment", chId)));
         TH2F * TrackThetaPosvsAngleandTrigH    = getHisto<TH2F>(igetter.get(getMEName("TrackThetaPosvsAngleandTrigH","Segment", chId)));
         
         if (TrackThetaPosvsAngle && TrackThetaPosvsAngleandTrig && TrackThetaPosvsAngleandTrigH && TrackThetaPosvsAngle->GetEntries()>1) {
           
           if( chambME[indexCh].find(fullName("TrigEffAngleTheta")) == chambME[indexCh].end()){
 
         bookChambHistos(ibooker,chId,"TrigEffPosvsAngleTheta");
         bookChambHistos(ibooker,chId,"TrigEffPosvsAngleHTheta");
         bookChambHistos(ibooker,chId,"TrigEffPosTheta");
         bookChambHistos(ibooker,chId,"TrigEffPosHTheta");
         bookChambHistos(ibooker,chId,"TrigEffAngleTheta");
         bookChambHistos(ibooker,chId,"TrigEffAngleHTheta");
           }
           if( secME[sector_id].find(fullName("TrigEffTheta")) == secME[sector_id].end() ){ 
         bookSectorHistos(ibooker,wh,sect,"TrigEffTheta");  
           }
           if( whME[wh].find(fullName("TrigEffTheta")) == whME[wh].end() ){
         bookWheelHistos(ibooker,wh,"TrigEffTheta");  
         bookWheelHistos(ibooker,wh,"TrigEffHTheta");  
           }
 
 
           std::map<std::string,MonitorElement*> *innerME = &(secME[sector_id]);
           TH1D* TrackThetaPos               = TrackThetaPosvsAngle->ProjectionY();
           TH1D* TrackThetaAngle             = TrackThetaPosvsAngle->ProjectionX();
           TH1D* TrackThetaPosandTrig        = TrackThetaPosvsAngleandTrig->ProjectionY();
           TH1D* TrackThetaAngleandTrig      = TrackThetaPosvsAngleandTrig->ProjectionX();
           TH1D* TrackThetaPosandTrigH       = TrackThetaPosvsAngleandTrigH->ProjectionY();
           TH1D* TrackThetaAngleandTrigH     = TrackThetaPosvsAngleandTrigH->ProjectionX();
           float binEff  = float(TrackThetaPosandTrig->GetEntries())/TrackThetaPos->GetEntries();
           float binErr  = sqrt(binEff*(1-binEff)/TrackThetaPos->GetEntries());
           float binEffH = float(TrackThetaPosandTrigH->GetEntries())/TrackThetaPos->GetEntries();
           float binErrH = sqrt(binEffH*(1-binEffH)/TrackThetaPos->GetEntries());
       
           MonitorElement* globalEff = innerME->find(fullName("TrigEffTheta"))->second;
           globalEff->setBinContent(stat,binEff);
           globalEff->setBinError(stat,binErr);
 
           innerME = &(whME[wh]);
           globalEff = innerME->find(fullName("TrigEffTheta"))->second;
           globalEff->setBinContent(sect,stat,binEff);
           globalEff->setBinError(sect,stat,binErr);
           globalEff = innerME->find(fullName("TrigEffHTheta"))->second;
           globalEff->setBinContent(sect,stat,binEffH);
           globalEff->setBinError(sect,stat,binErrH);
       
           innerME = &(chambME[indexCh]);
           makeEfficiencyME(TrackThetaPosandTrig,TrackThetaPos,innerME->find(fullName("TrigEffPosTheta"))->second);
           makeEfficiencyME(TrackThetaPosandTrigH,TrackThetaPos,innerME->find(fullName("TrigEffPosHTheta"))->second);
           makeEfficiencyME(TrackThetaAngleandTrig,TrackThetaAngle,innerME->find(fullName("TrigEffAngleTheta"))->second);
           makeEfficiencyME(TrackThetaAngleandTrigH,TrackThetaAngle,innerME->find(fullName("TrigEffAngleHTheta"))->second);
           makeEfficiencyME2D(TrackThetaPosvsAngleandTrig,TrackThetaPosvsAngle,innerME->find(fullName("TrigEffPosvsAngleTheta"))->second);
           makeEfficiencyME2D(TrackThetaPosvsAngleandTrigH,TrackThetaPosvsAngle,innerME->find(fullName("TrigEffPosvsAngleHTheta"))->second);      
         }
 
       }
     }
       }
     }
   } 
 
 }