#include <DTSegment2DSLPhiQuality.h>

Inheritance diagram for DTSegment2DSLPhiQuality:

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

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

void	endJob ()

void	endLuminosityBlock (edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &c)

virtual	~DTSegment2DSLPhiQuality ()
	Destructor. More...

Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

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

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

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

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
DQMStore *	dbe_

bool	debug

bool	doall

HEff2DHit *	h2DHitEff_SuperPhi

HRes2DHit *	h2DHitSuperPhi

bool	local

std::string	rootFileName

edm::InputTag	segment4DLabel

double	sigmaResAngle

double	sigmaResPos

edm::InputTag	simHitLabel

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

typedef WorkerT< EDAnalyzer >	WorkerType

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::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

CurrentProcessingContext const *	currentContext () const

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 reconstructed with both SL Phi and plot resolution comparing reconstructed and simulated quantities

Date:: 2010/09/13 09:49:18

Revision:: 1.7

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

Definition at line 32 of file DTSegment2DSLPhiQuality.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 42 of file DTSegment2DSLPhiQuality.cc.

References dbe_, debug, DTHitQualityUtils::debug, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), and cppFunctionSkipper::operator.

                                                                           {
    // 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");
   // the name of the 4D rec hit collection
   segment4DLabel = pset.getUntrackedParameter<InputTag>("segment4DLabel");
 
   //sigma resolution on position
   sigmaResPos = pset.getParameter<double>("sigmaResPos");
   //sigma resolution on angle
   sigmaResAngle = pset.getParameter<double>("sigmaResAngle");
   doall = pset.getUntrackedParameter<bool>("doall", false);
   local = pset.getUntrackedParameter<bool>("local", false);
 
   // Create the root file
   //theFile = new TFile(rootFileName.c_str(), "RECREATE");
   //theFile->cd();
   // ----------------------                 
   // 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();
   }
 
   // Book the histos
   h2DHitSuperPhi = new HRes2DHit ("SuperPhi",dbe_,doall,local);
   if(doall) h2DHitEff_SuperPhi = new HEff2DHit ("SuperPhi",dbe_);
 }

DTSegment2DSLPhiQuality::~DTSegment2DSLPhiQuality ( )

virtual

Destructor.

Definition at line 83 of file DTSegment2DSLPhiQuality.cc.

83 {

84 }

Member Function Documentation

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

virtual

Perform the real analysis.

Implements edm::EDAnalyzer.

Definition at line 106 of file DTSegment2DSLPhiQuality.cc.

References gather_cfg::cout, debug, PV3DBase< T, PVType, FrameType >::eta(), 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(), trackerHits::simHits, mathSSE::sqrt(), 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.getByLabel(simHitLabel, simHits); //FIXME: second string to be removed
 
   //Map simHits by chamber
   map<DTChamberId, PSimHitContainer > simHitsPerCh;
   for(PSimHitContainer::const_iterator simHit = simHits->begin();
       simHit != simHits->end(); simHit++){
     // Create the id of the chamber (the simHits in the DT known their wireId)
     DTChamberId chamberId = (((DTWireId(simHit->detUnitId())).layerId()).superlayerId()).chamberId();
     // Fill the map
     simHitsPerCh[chamberId].push_back(*simHit);
   }
 
   // Get the 4D rechits from the event
   Handle<DTRecSegment4DCollection> segment4Ds;
   event.getByLabel(segment4DLabel, segment4Ds);
 
   if(!segment4Ds.isValid()) {
     if(debug) cout << "[DTSegment2DSLPhiQuality]**Warning: no 4D Segments with label: " << segment4DLabel
    << " in this event, skipping!" << endl;
     return;
   }
     
   // Loop over all chambers containing a segment
   DTRecSegment4DCollection::id_iterator chamberId;
   for (chamberId = segment4Ds->id_begin();
        chamberId != segment4Ds->id_end();
        ++chamberId){
     
     //------------------------- simHits ---------------------------//
     //Get simHits of each chamber
     PSimHitContainer simHits =  simHitsPerCh[(*chamberId)];
        
     // 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 << "[DTSegment2DSLPhiQuality] Only " << nMuSimHit << " mu SimHit in this chamber, skipping!" << endl;
       continue; // If no or only one mu SimHit is found skip this chamber
     } 
     if(debug)
       cout << "=== Chamber " << (*chamberId) << " has " << nMuSimHit << " SimHits" << endl;
   
     //Find outer and inner mu SimHit to build a segment
     pair<const PSimHit*, const PSimHit*> inAndOutSimHit = DTHitQualityUtils::findMuSimSegment(muSimHitPerWire); 
 
     //Find direction and position of the sim Segment in Chamber RF
     pair<LocalVector, LocalPoint> dirAndPosSimSegm = DTHitQualityUtils::findMuSimSegmentDirAndPos(inAndOutSimHit,
                                                   (*chamberId),&(*dtGeom));
 
     LocalVector simSegmLocalDir = dirAndPosSimSegm.first;
     LocalPoint simSegmLocalPos = dirAndPosSimSegm.second;
     const DTChamber* chamber = dtGeom->chamber(*chamberId);
     GlobalPoint simSegmGlobalPos = chamber->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 lobal RF
     float etaSimSeg = simSegmGlobalPos.eta(); 
     float phiSimSeg = simSegmGlobalPos.phi();
 
     if(debug)
       cout<<"  Simulated segment:  local direction "<<simSegmLocalDir<<endl
       <<"                      local position  "<<simSegmLocalPos<<endl
       <<"                      angle           "<<angleSimSeg<<endl;
     
     //---------------------------- recHits --------------------------//
     // Get the range of rechit for the corresponding chamberId
     bool recHitFound = false;
     DTRecSegment4DCollection::range range = segment4Ds->get(*chamberId);
     int nsegm = distance(range.first, range.second);
     if(debug)
       cout << "   Chamber: " << *chamberId << " has " << nsegm
        << " 4D segments" << endl;
 
     if (nsegm!=0) {
       // Find the best RecHit: look for the 4D RecHit with the phi 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 chamberId
       for (DTRecSegment4DCollection::const_iterator segment4D = range.first;
            segment4D!=range.second;
            ++segment4D){
         // Check the dimension
         if((*segment4D).dimension() != 4) {
           if(debug) cout << "[DTSegment2DSLPhiQuality]***Error: This is not 4D segment!!!" << endl;
           continue;
         }
 
         //Get 2D superPhi segments from 4D segments
         const DTChamberRecSegment2D* phiSegment2D = (*segment4D).phiSegment();
         if((*phiSegment2D).dimension() != 2) {
           if(debug) cout << "[DTSegment2DQuality]***Error: This is not 2D segment!!!" << endl;
           abort();
         }
 
         // Segment Local Direction and position (in Chamber RF)
         LocalVector recSegDirection = (*phiSegment2D).localDirection();
 
         float recSegAlpha = DTHitQualityUtils::findSegmentAlphaAndBeta(recSegDirection).first;
         if(debug)
           cout << "  RecSegment direction: " << recSegDirection << endl
             << "             position : " <<  (*phiSegment2D).localPosition() << endl
             << "             alpha    : " << recSegAlpha << endl;
 
         if(fabs(recSegAlpha - angleSimSeg) < deltaAlpha) {
           deltaAlpha = fabs(recSegAlpha - angleSimSeg);
           bestRecHit = &(*phiSegment2D);
           bestRecHitFound = true;
         }
       }  // End of Loop over all 4D RecHits of this chambers
 
       if(bestRecHitFound) {
         // Best rechit direction and position in Chamber 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
         h2DHitSuperPhi->Fill(angleSimSeg,
                             angleBestRHit,
                             posSimSeg,
                             bestRecHitLocalPos.x(),
                             etaSimSeg,
                             phiSimSeg,
                             sqrt(bestRecHitLocalPosErr.xx()),
                             sqrt(bestRecHitLocalDirErr.xx())
                            );
       }
     } //end of if(nsegm!=0)
 
       // Fill Efficiency plot
     if(doall) {h2DHitEff_SuperPhi->Fill(etaSimSeg,
                             phiSimSeg,
                             posSimSeg,
                             angleSimSeg,
                     recHitFound);}
   } // End of loop over chambers
 }

void DTSegment2DSLPhiQuality::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 93 of file DTSegment2DSLPhiQuality.cc.

                                      {
   // Write the histos to file
   //theFile->cd();
   //h2DHitSuperPhi->Write();
 
   if(doall) h2DHitEff_SuperPhi->ComputeEfficiency();
   //h2DHitEff_SuperPhi->Write();
 
   //if ( rootFileName.size() != 0 && dbe_ ) dbe_->save(rootFileName); 
   //theFile->Close();
 } 

void DTSegment2DSLPhiQuality::endLuminosityBlock	(	edm::LuminosityBlock const &	lumiSeg,
		edm::EventSetup const &	c
	)

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 86 of file DTSegment2DSLPhiQuality.cc.

87 {

88 }

Member Data Documentation

DQMStore* DTSegment2DSLPhiQuality::dbe_

private

Definition at line 70 of file DTSegment2DSLPhiQuality.h.

bool DTSegment2DSLPhiQuality::debug

private

Definition at line 57 of file DTSegment2DSLPhiQuality.h.

bool DTSegment2DSLPhiQuality::doall

private

Definition at line 71 of file DTSegment2DSLPhiQuality.h.

HEff2DHit* DTSegment2DSLPhiQuality::h2DHitEff_SuperPhi

private

Definition at line 69 of file DTSegment2DSLPhiQuality.h.

HRes2DHit* DTSegment2DSLPhiQuality::h2DHitSuperPhi

private

Definition at line 68 of file DTSegment2DSLPhiQuality.h.

bool DTSegment2DSLPhiQuality::local

private

Definition at line 72 of file DTSegment2DSLPhiQuality.h.

std::string DTSegment2DSLPhiQuality::rootFileName

private

Definition at line 59 of file DTSegment2DSLPhiQuality.h.

edm::InputTag DTSegment2DSLPhiQuality::segment4DLabel

private

Definition at line 62 of file DTSegment2DSLPhiQuality.h.

double DTSegment2DSLPhiQuality::sigmaResAngle

private

Definition at line 66 of file DTSegment2DSLPhiQuality.h.

double DTSegment2DSLPhiQuality::sigmaResPos

private

Definition at line 64 of file DTSegment2DSLPhiQuality.h.

edm::InputTag DTSegment2DSLPhiQuality::simHitLabel

private

Definition at line 61 of file DTSegment2DSLPhiQuality.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation