#include <DTSegment2DQuality.h>

Inheritance diagram for DTSegment2DQuality:

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &eventSetup)
	Perform the real analysis. More...

virtual void	beginRun (const edm::Run &iRun, const edm::EventSetup &setup)

	DTSegment2DQuality (const edm::ParameterSet &pset)
	Constructor. More...

void	endJob ()

virtual	~DTSegment2DQuality ()
	Destructor. More...

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

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

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

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
DQMStore *	dbe_

bool	debug

HEff2DHit *	h2DHitEff_RPhi

HEff2DHit *	h2DHitEff_RZ

HEff2DHit *	h2DHitEff_RZ_W0

HEff2DHit *	h2DHitEff_RZ_W1

HEff2DHit *	h2DHitEff_RZ_W2

HRes2DHit *	h2DHitRPhi

HRes2DHit *	h2DHitRZ

HRes2DHit *	h2DHitRZ_W0

HRes2DHit *	h2DHitRZ_W1

HRes2DHit *	h2DHitRZ_W2

std::string	rootFileName

edm::InputTag	segment2DLabel

edm::EDGetTokenT< DTRecSegment2DCollection >	segment2DToken_

double	sigmaResAngle

double	sigmaResPos

edm::InputTag	simHitLabel

edm::EDGetTokenT< edm::PSimHitContainer >	simHitToken_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

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)

Detailed Description

Basic analyzer class which accesses 2D DTSegments and plot resolution comparing reconstructed and simulated quantities

Author: S. Bolognesi and G. Cerminara - INFN Torino

Definition at line 32 of file DTSegment2DQuality.h.

Constructor & Destructor Documentation

DTSegment2DQuality::DTSegment2DQuality ( const edm::ParameterSet & pset )

Constructor.

Definition at line 37 of file DTSegment2DQuality.cc.

References gather_cfg::cout, debug, DTHitQualityUtils::debug, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), DTRecHitQuality_cfi::segment2DLabel, DTRecHitQuality_cfi::sigmaResAngle, DTRecHitQuality_cfi::sigmaResPos, and AddPileupSummary_cfi::simHitLabel.

                                                                 {
    // Get the debug parameter for verbose output
   debug = pset.getUntrackedParameter<bool>("debug");
   DTHitQualityUtils::debug = debug;
   // the name of the simhit collection
   simHitLabel = pset.getUntrackedParameter<InputTag>("simHitLabel");
   simHitToken_ = consumes<PSimHitContainer>(pset.getUntrackedParameter<InputTag>("simHitLabel"));
   // the name of the 2D rec hit collection
   segment2DLabel = pset.getUntrackedParameter<InputTag>("segment2DLabel");
   segment2DToken_ = consumes<DTRecSegment2DCollection>(pset.getUntrackedParameter<InputTag>("segment2DLabel"));
 
   //sigma resolution on position
   sigmaResPos = pset.getParameter<double>("sigmaResPos");
   //sigma resolution on angle
   sigmaResAngle = pset.getParameter<double>("sigmaResAngle");
 
   if(debug)
     cout << "[DTSegment2DQuality] Constructor called " << endl;
 }

DTSegment2DQuality::~DTSegment2DQuality ( )

virtual

Destructor.

Definition at line 90 of file DTSegment2DQuality.cc.

                                        {
 
 }

Member Function Documentation

void DTSegment2DQuality::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)

Perform the real analysis.

Definition at line 121 of file DTSegment2DQuality.cc.

References funct::abs(), gather_cfg::cout, debug, SoftLeptonByDistance_cfi::distance, PV3DBase< T, PVType, FrameType >::eta(), HRes2DHit::Fill(), HEff2DHit::Fill(), DTHitQualityUtils::findMuSimSegment(), DTHitQualityUtils::findMuSimSegmentDirAndPos(), DTHitQualityUtils::findSegmentAlphaAndBeta(), edm::EventSetup::get(), edm::HandleBase::isValid(), DTRecSegment2D::localDirection(), DTRecSegment2D::localDirectionError(), DTRecSegment2D::localPosition(), DTRecSegment2D::localPositionError(), DTHitQualityUtils::mapMuSimHitsPerWire(), DTHitQualityUtils::mapSimHitsPerWire(), PV3DBase< T, PVType, FrameType >::phi(), DTRecHitQuality_cfi::segment2DLabel, DTRecHitQuality_cfi::sigmaResAngle, DTRecHitQuality_cfi::sigmaResPos, rpcPointValidation_cfi::simHit, trackerHits::simHits, mathSSE::sqrt(), DTGeometry::superLayer(), GeomDet::toGlobal(), PV3DBase< T, PVType, FrameType >::x(), and LocalError::xx().

                                                                                  {
   //theFile->cd();
 
   // Get the DT Geometry
   ESHandle<DTGeometry> dtGeom;
   eventSetup.get<MuonGeometryRecord>().get(dtGeom);
 
   // Get the SimHit collection from the event
   edm::Handle<PSimHitContainer> simHits;
   event.getByToken(simHitToken_, simHits); //FIXME: second string to be removed
 
   //Map simHits by sl
   map<DTSuperLayerId, PSimHitContainer > simHitsPerSl;
   for(PSimHitContainer::const_iterator simHit = simHits->begin();
       simHit != simHits->end(); simHit++){
     // Create the id of the sl (the simHits in the DT known their wireId)
     DTSuperLayerId slId = ((DTWireId(simHit->detUnitId())).layerId()).superlayerId();
     // Fill the map
     simHitsPerSl[slId].push_back(*simHit);
   }
 
   // Get the 2D rechits from the event
   Handle<DTRecSegment2DCollection> segment2Ds;
   event.getByToken(segment2DToken_, segment2Ds);
 
   if(!segment2Ds.isValid()) {
     if(debug) cout << "[DTSegment2DQuality]**Warning: no 2DSegments with label: " << segment2DLabel
            << " in this event, skipping!" << endl;
     return;
   }    
 
   // Loop over all superlayers containing a segment
   DTRecSegment2DCollection::id_iterator slId;
   for (slId = segment2Ds->id_begin();
        slId != segment2Ds->id_end();
        ++slId){
 
       //------------------------- simHits ---------------------------//
     //Get simHits of each superlayer
     PSimHitContainer simHits =  simHitsPerSl[(*slId)];
        
     // Map simhits per wire
     map<DTWireId, PSimHitContainer > simHitsPerWire = DTHitQualityUtils::mapSimHitsPerWire(simHits);
     map<DTWireId, const PSimHit*> muSimHitPerWire = DTHitQualityUtils::mapMuSimHitsPerWire(simHitsPerWire);
     int nMuSimHit = muSimHitPerWire.size();
     if(nMuSimHit == 0 || nMuSimHit == 1) {
       if(debug && nMuSimHit == 1)
         cout << "[DTSegment2DQuality] Only " << nMuSimHit << " mu SimHit in this SL, skipping!" << endl;
       continue; // If no or only one mu SimHit is found skip this SL
     } 
     if(debug)
       cout << "=== SL " << (*slId) << " has " << nMuSimHit << " SimHits" << endl;
   
     //Find outer and inner mu SimHit to build a segment
     pair<const PSimHit*, const PSimHit*> inAndOutSimHit = DTHitQualityUtils::findMuSimSegment(muSimHitPerWire); 
     //Check that outermost and innermost SimHit are not the same
     if(inAndOutSimHit.first ==inAndOutSimHit.second ) {
       cout << "[DTHitQualityUtils]***Warning: outermost and innermost SimHit are the same!" << endl;
       continue;
     }
 
     //Find direction and position of the sim Segment in SL RF
     pair<LocalVector, LocalPoint> dirAndPosSimSegm = DTHitQualityUtils::findMuSimSegmentDirAndPos(inAndOutSimHit,
                                                   (*slId),&(*dtGeom));
 
     LocalVector simSegmLocalDir = dirAndPosSimSegm.first;
     LocalPoint simSegmLocalPos = dirAndPosSimSegm.second;
     if(debug)
       cout<<"  Simulated segment:  local direction "<<simSegmLocalDir<<endl
         <<"                      local position  "<<simSegmLocalPos<<endl;
     const DTSuperLayer* superLayer = dtGeom->superLayer(*slId);
     GlobalPoint simSegmGlobalPos = superLayer->toGlobal(simSegmLocalPos);
 
     //Atan(x/z) angle and x position in SL RF
     float angleSimSeg = DTHitQualityUtils::findSegmentAlphaAndBeta(simSegmLocalDir).first;
     float posSimSeg = simSegmLocalPos.x(); 
     //Position (in eta,phi coordinates) in the global RF
     float etaSimSeg = simSegmGlobalPos.eta(); 
     float phiSimSeg = simSegmGlobalPos.phi();
 
 
     //---------------------------- recHits --------------------------//
     // Get the range of rechit for the corresponding slId
     bool recHitFound = false;
     DTRecSegment2DCollection::range range = segment2Ds->get(*slId);
     int nsegm = distance(range.first, range.second);
     if(debug)
       cout << "   Sl: " << *slId << " has " << nsegm
        << " 2D segments" << endl;
 
     if (nsegm!=0) {
       // Find the best RecHit: look for the 2D RecHit with the angle closest
       // to that of segment made of SimHits. 
       // RecHits must have delta alpha and delta position within 5 sigma of
       // the residual distribution (we are looking for residuals of segments
       // usefull to the track fit) for efficency purpose
       const DTRecSegment2D* bestRecHit = 0;
       bool bestRecHitFound = false;
       double deltaAlpha = 99999;
   
       // Loop over the recHits of this slId
       for (DTRecSegment2DCollection::const_iterator segment2D = range.first;
        segment2D!=range.second;
        ++segment2D){
     // Check the dimension
     if((*segment2D).dimension() != 2) {
       if(debug) cout << "[DTSegment2DQuality]***Error: This is not 2D segment!!!" << endl;
       abort();
     }
     // Segment Local Direction and position (in SL RF)
     LocalVector recSegDirection = (*segment2D).localDirection();
     LocalPoint recSegPosition = (*segment2D).localPosition();
       
     float recSegAlpha = DTHitQualityUtils::findSegmentAlphaAndBeta(recSegDirection).first;
     if(debug)
       cout << "  RecSegment direction: " << recSegDirection << endl
          << "             position : " << recSegPosition << endl
          << "             alpha    : " << recSegAlpha << endl;
       
     if(fabs(recSegAlpha - angleSimSeg) < deltaAlpha) {
       deltaAlpha = fabs(recSegAlpha - angleSimSeg);
       bestRecHit = &(*segment2D);
       bestRecHitFound = true;
     }
       }  // End of Loop over all 2D RecHits
 
       if(bestRecHitFound) {
     // Best rechit direction and position in SL RF
     LocalPoint bestRecHitLocalPos = bestRecHit->localPosition();
     LocalVector bestRecHitLocalDir = bestRecHit->localDirection();
 
     LocalError bestRecHitLocalPosErr = bestRecHit->localPositionError();
     LocalError bestRecHitLocalDirErr = bestRecHit->localDirectionError();
 
     float angleBestRHit = DTHitQualityUtils::findSegmentAlphaAndBeta(bestRecHitLocalDir).first;
 
     if(fabs(angleBestRHit - angleSimSeg) < 5*sigmaResAngle &&
        fabs(bestRecHitLocalPos.x() - posSimSeg) < 5*sigmaResPos) {
       recHitFound = true;
     }
 
     // Fill Residual histos
     HRes2DHit *hRes = 0;
     if((*slId).superlayer() == 1 || (*slId).superlayer() == 3) { //RPhi SL
       hRes = h2DHitRPhi;
     } else if((*slId).superlayer() == 2) { //RZ SL
       h2DHitRZ->Fill(angleSimSeg,
                     angleBestRHit,
                     posSimSeg,
                     bestRecHitLocalPos.x(),
                     etaSimSeg,
                     phiSimSeg,
                     sqrt(bestRecHitLocalPosErr.xx()),
                     sqrt(bestRecHitLocalDirErr.xx())
                     );
       if(abs((*slId).wheel()) == 0)
         hRes = h2DHitRZ_W0;
       else if(abs((*slId).wheel()) == 1)
         hRes = h2DHitRZ_W1;
       else if(abs((*slId).wheel()) == 2)
         hRes = h2DHitRZ_W2;
     }
     hRes->Fill(angleSimSeg,
                angleBestRHit,
                posSimSeg,
                bestRecHitLocalPos.x(),
                etaSimSeg,
                phiSimSeg,
                sqrt(bestRecHitLocalPosErr.xx()),
                sqrt(bestRecHitLocalDirErr.xx())
                );
       }
     } //end of if(nsegm!=0)
 
       // Fill Efficiency plot
     HEff2DHit *hEff = 0;
     if((*slId).superlayer() == 1 || (*slId).superlayer() == 3) { //RPhi SL
       hEff = h2DHitEff_RPhi;
     } else if((*slId).superlayer() == 2) { //RZ SL
       h2DHitEff_RZ->Fill(etaSimSeg, phiSimSeg, posSimSeg, angleSimSeg, recHitFound);
       if(abs((*slId).wheel()) == 0)
     hEff = h2DHitEff_RZ_W0;
       else if(abs((*slId).wheel()) == 1)
     hEff = h2DHitEff_RZ_W1;
       else if(abs((*slId).wheel()) == 2)
     hEff = h2DHitEff_RZ_W2;
     }
     hEff->Fill(etaSimSeg, phiSimSeg, posSimSeg, angleSimSeg, recHitFound);
   } // End of loop over superlayers
 }

void DTSegment2DQuality::beginRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	setup
	)

virtual

Definition at line 57 of file DTSegment2DQuality.cc.

References gather_cfg::cout, dbe_, debug, and Utilities::operator.

                                                                                 {
 
   // ----------------------                 
   // get hold of back-end interface 
   dbe_ = 0;
   dbe_ = Service<DQMStore>().operator->();
   if ( dbe_ ) {
     if (debug) {
       dbe_->setVerbose(1);
     } else {
       dbe_->setVerbose(0);
     }
   }
   if ( dbe_ ) {
     if ( debug ) dbe_->showDirStructure();
   }
 
   h2DHitRPhi = new HRes2DHit ("RPhi",dbe_,true,true);
   h2DHitRZ= new HRes2DHit ("RZ",dbe_,true,true);
   h2DHitRZ_W0= new HRes2DHit ("RZ_W0",dbe_,true,true);
   h2DHitRZ_W1= new HRes2DHit ("RZ_W1",dbe_,true,true);
   h2DHitRZ_W2= new HRes2DHit ("RZ_W2",dbe_,true,true);
 
   h2DHitEff_RPhi= new HEff2DHit ("RPhi",dbe_);
   h2DHitEff_RZ= new HEff2DHit ("RZ",dbe_);
   h2DHitEff_RZ_W0= new HEff2DHit ("RZ_W0",dbe_);
   h2DHitEff_RZ_W1= new HEff2DHit ("RZ_W1",dbe_);
   h2DHitEff_RZ_W2= new HEff2DHit ("RZ_W2",dbe_);
   if(debug)
     cout << "[DTSegment2DQuality] hitsos created " << endl;
 }

void DTSegment2DQuality::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 94 of file DTSegment2DQuality.cc.

                                 {
   // Write the histos to file
   //theFile->cd();
 
   //h2DHitRPhi->Write();
   //h2DHitRZ->Write();
   //h2DHitRZ_W0->Write();
   //h2DHitRZ_W1->Write();
   //h2DHitRZ_W2->Write();
 
   h2DHitEff_RPhi->ComputeEfficiency();
   h2DHitEff_RZ->ComputeEfficiency();
   h2DHitEff_RZ_W0->ComputeEfficiency();
   h2DHitEff_RZ_W1->ComputeEfficiency();
   h2DHitEff_RZ_W2->ComputeEfficiency();
 
   //h2DHitEff_RPhi->Write();
   //h2DHitEff_RZ->Write();
   //h2DHitEff_RZ_W0->Write();
   //h2DHitEff_RZ_W1->Write();
   //h2DHitEff_RZ_W2->Write();
   //if ( rootFileName.size() != 0 && dbe_ ) dbe_->save(rootFileName); 
 
   //theFile->Close();
 } 

Member Data Documentation

DQMStore* DTSegment2DQuality::dbe_

private

Definition at line 82 of file DTSegment2DQuality.h.

bool DTSegment2DQuality::debug

private

Definition at line 57 of file DTSegment2DQuality.h.

HEff2DHit* DTSegment2DQuality::h2DHitEff_RPhi

private

Definition at line 77 of file DTSegment2DQuality.h.

HEff2DHit* DTSegment2DQuality::h2DHitEff_RZ

private

Definition at line 78 of file DTSegment2DQuality.h.

HEff2DHit* DTSegment2DQuality::h2DHitEff_RZ_W0

private

Definition at line 79 of file DTSegment2DQuality.h.

HEff2DHit* DTSegment2DQuality::h2DHitEff_RZ_W1

private

Definition at line 80 of file DTSegment2DQuality.h.

HEff2DHit* DTSegment2DQuality::h2DHitEff_RZ_W2

private

Definition at line 81 of file DTSegment2DQuality.h.

HRes2DHit* DTSegment2DQuality::h2DHitRPhi

private

Definition at line 71 of file DTSegment2DQuality.h.

HRes2DHit* DTSegment2DQuality::h2DHitRZ

private

Definition at line 72 of file DTSegment2DQuality.h.

HRes2DHit* DTSegment2DQuality::h2DHitRZ_W0

private

Definition at line 73 of file DTSegment2DQuality.h.

HRes2DHit* DTSegment2DQuality::h2DHitRZ_W1

private

Definition at line 74 of file DTSegment2DQuality.h.

HRes2DHit* DTSegment2DQuality::h2DHitRZ_W2

private

Definition at line 75 of file DTSegment2DQuality.h.

std::string DTSegment2DQuality::rootFileName

private

Definition at line 59 of file DTSegment2DQuality.h.

edm::InputTag DTSegment2DQuality::segment2DLabel

private

Definition at line 62 of file DTSegment2DQuality.h.

edm::EDGetTokenT<DTRecSegment2DCollection> DTSegment2DQuality::segment2DToken_

private

Definition at line 64 of file DTSegment2DQuality.h.

double DTSegment2DQuality::sigmaResAngle

private

Definition at line 69 of file DTSegment2DQuality.h.

double DTSegment2DQuality::sigmaResPos

private

Definition at line 67 of file DTSegment2DQuality.h.

edm::InputTag DTSegment2DQuality::simHitLabel

private

Definition at line 61 of file DTSegment2DQuality.h.

edm::EDGetTokenT<edm::PSimHitContainer> DTSegment2DQuality::simHitToken_

private

Definition at line 63 of file DTSegment2DQuality.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation