---

VisTrack Class Reference

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

Inheritance diagram for VisTrack:

List of all members.

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

---

Detailed Description

Definition at line 18 of file VisTrack.h.

---

Constructor & Destructor Documentation

VisTrack::VisTrack

(

const edm::ParameterSet &

iConfig

)

[explicit]

Definition at line 42 of file VisTrack.cc.

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

virtual VisTrack::~VisTrack

(

)

[inline, virtual]

Definition at line 26 of file VisTrack.h.

{}

---

Member Function Documentation

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

Implements edm::EDAnalyzer.

Definition at line 47 of file VisTrack.cc.

References IgCollection::addProperty(), IgAssociationSet::associate(), collection, edmplugin::standard::config(), IgCollection::create(), error, ETA, Exception, edm::EventSetup::get(), IgDataStorage::getAssociationSet(), IgDataStorage::getCollection(), inputTag_, edm::InputTag::instance(), edm::Service< T >::isAvailable(), edm::Handle< T >::isValid(), edm::ESHandle< T >::isValid(), it, edm::InputTag::label(), VisLocalPosition::localPosition(), P, PHI, edm::InputTag::process(), edm::ESHandle< T >::product(), PT, track, and tracks.

{
    edm::Service<IguanaService> config;
    if (! config.isAvailable ()) 
    {
        throw cms::Exception ("Configuration")
            << "VisTrack 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<reco::TrackCollection> collection;
    event.getByLabel (inputTag_, collection);

    edm::ESHandle<GlobalTrackingGeometry> geometry;
    eventSetup.get<GlobalTrackingGeometryRecord> ().get (geometry);

    edm::ESHandle<MagneticField> field;    
    eventSetup.get<IdealMagneticFieldRecord> ().get (field);

    if (collection.isValid () && geometry.isValid () && field.isValid ())
    {       
        IgDataStorage *storage = config->storage ();
        IgCollection &tracks = storage->getCollection ("Tracks_V1");

        IgProperty VTX = tracks.addProperty ("pos", IgV3d());
        IgProperty P   = tracks.addProperty ("dir", IgV3d());
        IgProperty PT  = tracks.addProperty ("pt", static_cast<double>(0.0)); 
        IgProperty PHI = tracks.addProperty ("phi", static_cast<double>(0.0));
        IgProperty ETA = tracks.addProperty ("eta", static_cast<double>(0.0));
        IgProperty CHARGE = tracks.addProperty ("charge", static_cast<double>(0.0));

        // Create a collection for track extras.
        IgCollection &extras = storage->getCollection ("Extras_V1");
        IgProperty IPOS = extras.addProperty ("pos_1", IgV3d());
        IgProperty IP   = extras.addProperty ("dir_1", IgV3d());
        IgProperty OPOS = extras.addProperty ("pos_2", IgV3d());
        IgProperty OP   = extras.addProperty ("dir_2", IgV3d());
 
        IgAssociationSet &trackExtras = storage->getAssociationSet ("TrackExtras_V1");

        // Create a collection for tracking rec hits.
        IgCollection &hits = storage->getCollection ("Hits_V1");
        IgProperty HIT_POS = hits.addProperty ("pos", IgV3d());

        IgAssociationSet &trackHits = storage->getAssociationSet ("TrackHits_V1");
 
        for (reco::TrackCollection::const_iterator track = collection->begin (), trackEnd = collection->end ();
             track != trackEnd; ++track)
        {
            IgCollectionItem item = tracks.create ();
            item[VTX] = IgV3d(static_cast<double>((*track).referencePoint().x()/100.),
                              static_cast<double>((*track).referencePoint().y()/100.),
                              static_cast<double>((*track).referencePoint().z()/100.));
            item[P]   = IgV3d(static_cast<double>((*track).px()),
                              static_cast<double>((*track).py()),
                              static_cast<double>((*track).pz()));
            item[PT]  = static_cast<double>((*track).pt ());
            item[PHI] = static_cast<double>((*track).phi ());
            item[ETA] = static_cast<double>((*track).eta ());
            item[CHARGE] = static_cast<double>((*track).charge ());
            
            if ((*track).innerOk () && (*track).outerOk ())
            {
                IgCollectionItem eitem = extras.create ();
                eitem[IPOS]  = IgV3d(static_cast<double>((*track).innerPosition().x()/100.0),
                                     static_cast<double>((*track).innerPosition().y()/100.0),
                                     static_cast<double>((*track).innerPosition().z()/100.0));
                eitem[IP]    = IgV3d(static_cast<double>((*track).innerMomentum().x()),
                                     static_cast<double>((*track).innerMomentum().y()),
                                     static_cast<double>((*track).innerMomentum().z()));
                eitem[OPOS]  = IgV3d(static_cast<double>((*track).outerPosition().x()/100.0),
                                     static_cast<double>((*track).outerPosition().y()/100.0),
                                     static_cast<double>((*track).outerPosition().z()/100.0));
                eitem[OP]     = IgV3d(static_cast<double>((*track).outerMomentum().x()),
                                      static_cast<double>((*track).outerMomentum ().y ()),
                                      static_cast<double>((*track).outerMomentum ().z ()));

                trackExtras.associate (item, eitem);
            }
            
            if ((*track).extra ().isNonnull ())
            {
                for (trackingRecHit_iterator it = track->recHitsBegin ();  it != track->recHitsEnd (); ++it) 
                {
                    if ((*it)->isValid () && !(*it)->geographicalId ().null ())
                    {
                        IgCollectionItem hit = hits.create ();
                        LocalPoint point = VisLocalPosition::localPosition(&(**it), geometry.product());
                        hit[HIT_POS] = IgV3d(static_cast<double>(geometry->idToDet((*it)->geographicalId())->surface().toGlobal(point).x()/100.0),
                                             static_cast<double>(geometry->idToDet((*it)->geographicalId())->surface().toGlobal(point).y()/100.0),
                                             static_cast<double>(geometry->idToDet((*it)->geographicalId())->surface().toGlobal(point).z()/100.0));

                        trackHits.associate (item, hit);
                    }
                }
            }
        }
    }
    else 
    {
        // friendlyName:moduleLabel:instanceName:processName
        std::string error = "### Error: Tracks "
                            + edm::TypeID (typeid (reco::TrackCollection)).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 VisTrack::inputTag_ [private]

Definition at line 31 of file VisTrack.h.

Referenced by analyze().

---

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

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

---