---

VisDTRecSegment4D Class Reference

#include <VisReco/Analyzer/interface/VisDTRecSegment4D.h>

Inheritance diagram for VisDTRecSegment4D:

List of all members.

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
	VisDTRecSegment4D (const edm::ParameterSet &)
virtual	~VisDTRecSegment4D (void)
Private Attributes
edm::InputTag	inputTag_

---

Detailed Description

Definition at line 18 of file VisDTRecSegment4D.h.

---

Constructor & Destructor Documentation

VisDTRecSegment4D::VisDTRecSegment4D

(

const edm::ParameterSet &

iConfig

)

[explicit]

Definition at line 30 of file VisDTRecSegment4D.cc.

    : inputTag_ (iConfig.getParameter<edm::InputTag>("visDTRecSegment4DTag"))
{}

virtual VisDTRecSegment4D::~VisDTRecSegment4D

(

void

)

[inline, virtual]

Definition at line 26 of file VisDTRecSegment4D.h.

{}

---

Member Function Documentation

void VisDTRecSegment4D::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)			[virtual]

Implements edm::EDAnalyzer.

Definition at line 35 of file VisDTRecSegment4D.cc.

References IgCollection::addProperty(), collection, edmplugin::standard::config(), funct::cos(), IgCollection::create(), end, error, Exception, edm::EventSetup::get(), IgDataStorage::getCollection(), inputTag_, edm::InputTag::instance(), edm::Service< T >::isAvailable(), edm::Handle< T >::isValid(), edm::ESHandle< T >::isValid(), it, edm::InputTag::label(), PV3DBase< T, PVType, FrameType >::mag(), edm::InputTag::process(), PV3DBase< T, PVType, FrameType >::theta(), PV3DBase< T, PVType, FrameType >::x(), x, PV3DBase< T, PVType, FrameType >::y(), y, PV3DBase< T, PVType, FrameType >::z(), and z.

{
    edm::Service<IguanaService> config;
    if (! config.isAvailable ()) 
    {
        throw cms::Exception ("Configuration")
            << "VisDTRecSegment4D requires the IguanaService\n"
            "which is not present in the configuration file.\n"
            "You must add the service in the configuration file\n"
            "or remove the module that requires it";
    }
    
    edm::Handle<DTRecSegment4DCollection> collection;
    event.getByLabel (inputTag_, collection);

    edm::ESHandle<DTGeometry> geom;
    eventSetup.get<MuonGeometryRecord> ().get (geom);

    if (collection.isValid () && geom.isValid ())
    {       
        IgDataStorage *storage = config->storage ();
        IgCollection &segments = storage->getCollection ("DTRecSegment4D_V1");
        IgProperty DET_ID   = segments.addProperty ("detid", int (0)); 
        IgProperty POS_1    = segments.addProperty ("pos_1", IgV3d());
        IgProperty POS_2    = segments.addProperty ("pos_2", IgV3d());
        
        IgProperty SECTOR_ID = segments.addProperty("sectorId", int(0));
        IgProperty STATION_ID = segments.addProperty("stationId", int(0));
        IgProperty WHEEL_ID = segments.addProperty("wheelId", int(0));

        DTRecSegment4DCollection::const_iterator it = collection->begin ();
        DTRecSegment4DCollection::const_iterator end = collection->end ();
        for (; it != end; ++it) 
        {
            DTChamberId chId ((*it).geographicalId ().rawId ());
            const DTChamber *chamber = geom->chamber (chId);

            float halfHeight = chamber->surface ().bounds ().thickness () / 2.0;
            // float halfWidth = chamber->surface ().bounds ().width () / 2.0;
            LocalVector locDir = (*it).localDirection ();
            LocalPoint locPos = (*it).localPosition ();

            GlobalPoint posInner = chamber->toGlobal (locPos + locDir / locDir.mag () * halfHeight / cos (locDir.theta ()));              
            float x = posInner.x () / 100.0;  // cm -> m
            float y = posInner.y () / 100.0;  // cm -> m
            float z = posInner.z () / 100.0;  // cm -> m

            IgCollectionItem isegment = segments.create ();
            isegment[DET_ID] = static_cast<int> ((*it).geographicalId ().rawId ());
            isegment[POS_1]  = IgV3d(static_cast<double>(x), static_cast<double>(y), static_cast<double>(z));

            GlobalPoint posOuter = chamber->toGlobal (locPos + locDir / (-locDir.mag ()) * halfHeight / cos (locDir.theta ()));           
            x = posOuter.x () / 100.0;  // cm -> m
            y = posOuter.y () / 100.0;  // cm -> m
            z = posOuter.z () / 100.0;  // cm -> m
            isegment[POS_2]  = IgV3d(static_cast<double>(x), static_cast<double>(y), static_cast<double>(z));
            
            isegment[SECTOR_ID] = static_cast<int>(chId.sector());
            isegment[STATION_ID] = static_cast<int>(chId.station());
            isegment[WHEEL_ID] = static_cast<int>(chId.wheel());
        }
    }
    else 
    {
        // friendlyName:moduleLabel:instanceName:processName
        std::string error = "### Error: DTRecSegment4Ds "
                            + edm::TypeID (typeid (DTRecSegment4DCollection)).friendlyClassName () + ":" 
                            + inputTag_.label() + ":"
                            + inputTag_.instance() + ":" 
                            + inputTag_.process() + " are not found.";

        IgDataStorage *storage = config->storage ();
        IgCollection &collection = storage->getCollection ("Errors_V1");
        IgProperty ERROR_MSG = collection.addProperty ("Error", std::string ());
        IgCollectionItem item = collection.create ();
        item [ERROR_MSG] = error;
    }
}

---

Member Data Documentation

edm::InputTag VisDTRecSegment4D::inputTag_ [private]

Definition at line 31 of file VisDTRecSegment4D.h.

Referenced by analyze().

---

The documentation for this class was generated from the following files:

• VisReco/Analyzer/interface/VisDTRecSegment4D.h
• VisReco/Analyzer/src/VisDTRecSegment4D.cc

---