#include <ME0SegmentsValidation.h>

Inheritance diagram for ME0SegmentsValidation:

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

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

	ME0SegmentsValidation (const edm::ParameterSet &)

	~ME0SegmentsValidation ()

Public Member Functions inherited from ME0BaseValidation
MonitorElement *	BookHistXY (DQMStore::IBooker &, const char name, const char label, unsigned int region_num, unsigned int layer_num=99)

MonitorElement *	BookHistZR (DQMStore::IBooker &, const char name, const char label, unsigned int region_num, unsigned int layer_num=99)

	ME0BaseValidation (const edm::ParameterSet &ps)

virtual	~ME0BaseValidation ()

Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final

virtual void	beginStream (edm::StreamID id) final

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

	DQMEDAnalyzer (void)

virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

uint32_t	streamId () const

Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

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

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

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
edm::EDGetToken	InputTagToken_Segments

MonitorElement *	me0_rh_xy_Muon [2][6]

MonitorElement *	me0_segment_chi2

MonitorElement *	me0_segment_EtaRH

MonitorElement *	me0_segment_numRH

MonitorElement *	me0_segment_PhiRH

MonitorElement *	me0_segment_time

MonitorElement *	me0_segment_timeErr

MonitorElement *	me0_specRH_DeltaX [2][6]

MonitorElement *	me0_specRH_DeltaY [2][6]

MonitorElement *	me0_specRH_PullX [2][6]

MonitorElement *	me0_specRH_PullY [2][6]

MonitorElement *	me0_specRH_xy [2][6]

MonitorElement *	me0_specRH_zr [2]

Int_t	npart

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T...>	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T...>	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT < LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr < dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)

static std::shared_ptr < dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)

static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const , dqmDetails::NoCache )

static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const , dqmDetails::NoCache )

Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
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 ME0BaseValidation
edm::EDGetToken	InputTagToken_

std::vector< std::string >	layerLabel

int	nBinXY_

std::vector< double >	nBinZR_

std::vector< double >	RangeZR_

std::vector< std::string >	regionLabel

Detailed Description

Definition at line 12 of file ME0SegmentsValidation.h.

Constructor & Destructor Documentation

ME0SegmentsValidation::ME0SegmentsValidation ( const edm::ParameterSet & cfg )

explicit

Definition at line 5 of file ME0SegmentsValidation.cc.

References edm::ParameterSet::getParameter(), and InputTagToken_Segments.

                                                                       :  ME0BaseValidation(cfg)
 {
   InputTagToken_Segments = consumes<ME0SegmentCollection>(cfg.getParameter<edm::InputTag>("segmentInputLabel"));
 }

ME0SegmentsValidation::~ME0SegmentsValidation ( )

Definition at line 62 of file ME0SegmentsValidation.cc.

62 {

63 }

Member Function Documentation

void ME0SegmentsValidation::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	iSetup
	)

overridevirtual

Implements ME0BaseValidation.

Definition at line 66 of file ME0SegmentsValidation.cc.

References ME0Geometry::chamber(), beam_dqm_sourceclient-live_cfg::chi2, ME0Geometry::etaPartition(), MonitorElement::Fill(), edm::EventSetup::get(), edm::Event::getByToken(), InputTagToken_Segments, edm::HandleBase::isValid(), me0_segment_chi2, me0_segment_numRH, me0_segment_time, me0_segment_timeErr, me0_specRH_DeltaX, me0_specRH_DeltaY, me0_specRH_PullX, me0_specRH_PullY, me0_specRH_xy, me0_specRH_zr, ndof, HLT_25ns10e33_v2_cff::region, reco::return(), mathSSE::sqrt(), GeomDet::toGlobal(), x, PV3DBase< T, PVType, FrameType >::x(), and y.

 {
  edm::ESHandle<ME0Geometry> hGeom;
  iSetup.get<MuonGeometryRecord>().get(hGeom);
  const ME0Geometry* ME0Geometry_ =( &*hGeom);
 
   edm::Handle<ME0SegmentCollection> ME0Segments;
   e.getByToken(InputTagToken_Segments, ME0Segments);
 
   if (!ME0Segments.isValid() ) {
     edm::LogError("ME0SegmentsValidation") << "Cannot get ME0RecHits/ME0Segments by Token InputTagToken";
     return ;
   }
 
 
 
  for (auto me0s = ME0Segments->begin(); me0s != ME0Segments->end(); me0s++) {
 
    // The ME0 Ensamble DetId refers to layer = 1
    ME0DetId id = me0s->me0DetId();
    //std::cout <<" Original ME0DetID "<<id<<std::endl;
    auto chamber = ME0Geometry_->chamber(id);
    //std::cout <<"Global Segment Position "<< chamber->toGlobal(me0s->localPosition())<<std::endl;
    auto segLP = me0s->localPosition();
    auto segLD = me0s->localDirection();
    //std::cout <<" Global Direction theta = "<<segLD.theta()<<" phi="<<segLD.phi()<<std::endl;
    auto me0rhs = me0s->specificRecHits();
    //std::cout <<"ME0 Ensamble Det Id "<<id<<" Number of RecHits "<<me0rhs.size()<<std::endl;
 
 //   Int_t detId = id;
 //   Float_t localX = segLP.x();
 //   Float_t localY = segLP.y();
 //   Float_t localZ = segLP.z();
 //   Float_t dirTheta = segLD.theta();
 //   Float_t dirPhi = segLD.phi();
    Short_t numberRH = me0rhs.size();
    Float_t chi2 = (Float_t) me0s->chi2();
    Float_t ndof = me0s->degreesOfFreedom();
    Double_t time = me0s->time();   
    Double_t timeErr = me0s->timeErr();   
 
    Float_t reducedChi2 = chi2/ndof;
 
    me0_segment_chi2->Fill(reducedChi2);
    me0_segment_numRH->Fill(numberRH);
    me0_segment_time->Fill(time);
    me0_segment_timeErr->Fill(timeErr);
 
    for (auto rh = me0rhs.begin(); rh!= me0rhs.end(); rh++)
    {
 
      auto me0id = rh->me0Id();
      auto rhr = ME0Geometry_->etaPartition(me0id);
      auto rhLP = rh->localPosition();
      auto erhLEP = rh->localPositionError();
      auto rhGP = rhr->toGlobal(rhLP);
      auto rhLPSegm = chamber->toLocal(rhGP);
      float xe = segLP.x()+segLD.x()*rhLPSegm.z()/segLD.z();
      float ye = segLP.y()+segLD.y()*rhLPSegm.z()/segLD.z();
      float ze = rhLPSegm.z();
      LocalPoint extrPoint(xe,ye,ze); // in segment rest frame
      auto extSegm = rhr->toLocal(chamber->toGlobal(extrPoint)); // in layer restframe
 
 //     Int_t detId = me0id;
 
      Short_t region = me0id.region();
 //     Short_t station = 0;
 //     Short_t ring = 0;
      Short_t layer = me0id.layer();
 //     Short_t chamber = me0id.chamber();
 //     Short_t chamber = me0id.chamber();
 
      Float_t x = rhLP.x();
      Float_t xErr = erhLEP.xx();
      Float_t y = rhLP.y();
      Float_t yErr = erhLEP.yy();
 
      Float_t globalR = rhGP.perp();
      Float_t globalX = rhGP.x();
      Float_t globalY = rhGP.y();
      Float_t globalZ = rhGP.z();
 //     Float_t globalEta = rhGP.eta();
 //     Float_t globalPhi = rhGP.phi();
 
      Float_t xExt = extSegm.x();
      Float_t yExt = extSegm.y();
 
      Float_t pull_x = (x - xExt) / sqrt(xErr);
      Float_t pull_y = (y - yExt) / sqrt(yErr);
 
      int region_num=0 ;
      if ( region ==-1 ) region_num = 0 ;
      else if ( region==1) region_num = 1;
      int layer_num = layer-1;
 
      me0_specRH_xy[region_num][layer_num]->Fill(globalX,globalY);
      me0_specRH_zr[region_num]->Fill(globalZ,globalR);
 
      me0_specRH_DeltaX[region_num][layer_num]->Fill(x - xExt);
      me0_specRH_DeltaY[region_num][layer_num]->Fill(y - yExt);
      me0_specRH_PullX[region_num][layer_num]->Fill(pull_x);
      me0_specRH_PullY[region_num][layer_num]->Fill(pull_y);
 
 
   }
 }
 
 }

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 10 of file ME0SegmentsValidation.cc.

References DQMStore::IBooker::book1D(), ME0BaseValidation::BookHistXY(), ME0BaseValidation::BookHistZR(), ME0BaseValidation::layerLabel, LogDebug, me0_segment_chi2, me0_segment_numRH, me0_segment_time, me0_segment_timeErr, me0_specRH_DeltaX, me0_specRH_DeltaY, me0_specRH_PullX, me0_specRH_PullY, me0_specRH_xy, me0_specRH_zr, ME0BaseValidation::regionLabel, DQMStore::IBooker::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                                        {
    //edm::ESHandle<ME0Geometry> hGeom;
    //iSetup.get<MuonGeometryRecord>().get(hGeom);
    //const ME0Geometry* ME0Geometry_ =( &*hGeom);
 
   LogDebug("MuonME0SegmentsValidation")<<"Info : Loading Geometry information\n";
   ibooker.setCurrentFolder("MuonME0RecHitsV/ME0SegmentsTask");
 
   unsigned int nregion  = 2;
 
   edm::LogInfo("MuonME0SegmentsValidation")<<"+++ Info : # of region : "<<nregion<<std::endl;
 
   LogDebug("MuonME0SegmentsValidation")<<"+++ Info : finish to get geometry information from ES.\n";
 
   me0_segment_chi2    = ibooker.book1D("me0_seg_ReducedChi2","#chi^{2}/ndof; #chi^{2}/ndof; # Segments",100,0,5);
   me0_segment_numRH   = ibooker.book1D("me0_seg_NumberRH","Number of fitted RecHits; # RecHits; entries",11,-0.5,10.5);
   //me0_segment_EtaRH   = ibooker.book1D("me0_specRH_globalEta","Fitted RecHits Eta Distribution; #eta; entries",200,-4.0,4.0);
   //me0_segment_PhiRH   = ibooker.book1D("me0_specRH_globalPhi","Fitted RecHits Phi Distribution; #eta; entries",18,-3.14,3.14);
   me0_segment_time    = ibooker.book1D("me0_seg_time","Segment Timing; ns; entries",40,15,25);
   me0_segment_timeErr = ibooker.book1D("me0_seg_timErr","Segment Timing Error; ns; entries",50,0,0.5);
 
   for( unsigned int region_num = 0 ; region_num < nregion ; region_num++ ) {
       me0_specRH_zr[region_num] = BookHistZR(ibooker,"me0_specRH_tot","Segment RecHits",region_num);
       for( unsigned int layer_num = 0 ; layer_num < 6 ; layer_num++) {
           //me0_strip_dg_zr[region_num][layer_num] = BookHistZR(ibooker,"me0_strip_dg","SimHit",region_num,layer_num);
           me0_specRH_xy[region_num][layer_num] = BookHistXY(ibooker,"me0_specRH","Segment RecHits",region_num,layer_num);
           //me0_rh_xy_Muon[region_num][layer_num] = BookHistXY(ibooker,"me0_rh","RecHit Muon",region_num,layer_num);
 
           std::string histo_name_DeltaX = std::string("me0_specRH_DeltaX_r")+regionLabel[region_num]+"_l"+layerLabel[layer_num];
           std::string histo_name_DeltaY = std::string("me0_specRH_DeltaY_r")+regionLabel[region_num]+"_l"+layerLabel[layer_num];
           std::string histo_label_DeltaX = "Segment RecHits Delta X : region"+regionLabel[region_num]+
                 " layer "+layerLabel[layer_num]+" "+" ; x_{SimHit} - x_{Segment RecHits} ; entries";
           std::string histo_label_DeltaY = "Segment RecHits Delta Y : region"+regionLabel[region_num]+
                 " layer "+layerLabel[layer_num]+" "+" ; y_{SimHit} - y_{Segment RecHit} ; entries";
 
           me0_specRH_DeltaX[region_num][layer_num] = ibooker.book1D(histo_name_DeltaX.c_str(), histo_label_DeltaX.c_str(),100,-10,10);
           me0_specRH_DeltaY[region_num][layer_num] = ibooker.book1D(histo_name_DeltaY.c_str(), histo_label_DeltaY.c_str(),100,-10,10);
 
           std::string histo_name_PullX = std::string("me0_specRH_PullX_r")+regionLabel[region_num]+"_l"+layerLabel[layer_num];
           std::string histo_name_PullY = std::string("me0_specRH_PullY_r")+regionLabel[region_num]+"_l"+layerLabel[layer_num];
           std::string histo_label_PullX = "Segment RecHits Pull X : region"+regionLabel[region_num]+
                 " layer "+layerLabel[layer_num]+" "+" ; #frac{x_{SimHit} - x_{Segment RecHit}}{#sigma_{x,RecHit}} ; entries";
           std::string histo_label_PullY = "Segment RecHits Pull Y : region"+regionLabel[region_num]+
                 " layer "+layerLabel[layer_num]+" "+" ; #frac{y_{SimHit} - y_{Segment RecHit}}{#sigma_{y,RecHit}} ; entries";
 
           me0_specRH_PullX[region_num][layer_num] = ibooker.book1D(histo_name_PullX.c_str(), histo_label_DeltaX.c_str(),100,-10,10);
           me0_specRH_PullY[region_num][layer_num] = ibooker.book1D(histo_name_PullY.c_str(), histo_label_DeltaY.c_str(),100,-10,10);
       }
   }
 }