FWTrackingParticleProxyBuilderFullFramework.cc

Go to the documentation of this file.

#include "Fireworks/Core/interface/FWSimpleProxyBuilderTemplate.h"
#include "Fireworks/Core/interface/FWProxyBuilderBase.h"
#include "Fireworks/Core/interface/Context.h"
#include "Fireworks/Core/interface/FWEventItem.h"
#include "Fireworks/Core/interface/FWGeometry.h"
#include "SimDataFormats/TrackingAnalysis/interface/TrackingParticle.h"
#include "SimDataFormats/TrackingAnalysis/interface/TrackingParticleFwd.h"
#include "SimGeneral/TrackingAnalysis/interface/SimHitTPAssociationProducer.h"
#include "SimDataFormats/TrackingHit/interface/PSimHitContainer.h"
#include "DataFormats/TrackingRecHit/interface/TrackingRecHit.h"
 
#include "DataFormats/DetId/interface/DetId.h"
 
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Common/interface/EventBase.h"
#include "DataFormats/DetId/interface/DetId.h"
 
#include "Fireworks/Core/interface/FWProxyBuilderConfiguration.h"
#include "Fireworks/Core/interface/FWParameters.h"
 
#include "TEveTrack.h"
#include "TEveCompound.h"
#include "TEveManager.h"
#include "TEveBrowser.h"
#include "TEveTrackPropagator.h"
 
class FWTrackingParticleProxyBuilderFullFramework : public FWSimpleProxyBuilderTemplate<TrackingParticle> {
public:
  FWTrackingParticleProxyBuilderFullFramework(void) : m_assocList(nullptr) {}
  ~FWTrackingParticleProxyBuilderFullFramework(void) override {}
 
  // virtual void setItem(const FWEventItem* iItem) override;
 
  REGISTER_PROXYBUILDER_METHODS();
 
private:
  FWTrackingParticleProxyBuilderFullFramework(const FWTrackingParticleProxyBuilderFullFramework&) = delete;
  const FWTrackingParticleProxyBuilderFullFramework& operator=(const FWTrackingParticleProxyBuilderFullFramework&) =
      delete;
 
  using FWSimpleProxyBuilderTemplate<TrackingParticle>::build;
  void build(const FWEventItem* iItem, TEveElementList* product, const FWViewContext*) override;
 
  void build(const TrackingParticle& iData,
             unsigned int iIndex,
             TEveElement& oItemHolder,
             const FWViewContext*) override;
 
  edm::Handle<TrackingParticleCollection> tpch;
  const SimHitTPAssociationProducer::SimHitTPAssociationList* m_assocList;
};
 
//______________________________________________________________________________
 
/*
  void FWTrackingParticleProxyBuilderFullFramework::setItem(const FWEventItem* iItem)
  {
  printf("set item\n");
  FWProxyBuilderBase::setItem(iItem);
  }
*/
//______________________________________________________________________________
void FWTrackingParticleProxyBuilderFullFramework::build(const FWEventItem* iItem,
                                                        TEveElementList* product,
                                                        const FWViewContext*) {
  // setup event handles amd call function from parent class
 
  const edm::Event* event = (const edm::Event*)item()->getEvent();
  if (event) {
    // get collection handle
    edm::InputTag coltag(item()->moduleLabel(), item()->productInstanceLabel(), item()->processName());
    event->getByLabel(coltag, tpch);
 
    // AMT todo: check if there is any other way getting the list other than this
    //           ifnot, set proces name as a configurable parameter
    edm::Handle<SimHitTPAssociationProducer::SimHitTPAssociationList> simHitsTPAssoc;
    try {
      event->getByLabel("xxx", simHitsTPAssoc);
      m_assocList = &*simHitsTPAssoc;
    } catch (const std::exception& e) {
      std::cerr << " FWTrackingParticleProxyBuilderFullFramework::setItem() Can't get hits association list "
                << e.what() << std::endl;
    }
    /*
      // debug propagator
      gEve->GetBrowser()->MapWindow();
      gEve->AddToListTree(context().getTrackPropagator(), true);
      context().getTrackPropagator()->SetRnrReferences(true);
      */
  }
  FWSimpleProxyBuilder::build(iItem, product, nullptr);
}
//______________________________________________________________________________
void FWTrackingParticleProxyBuilderFullFramework::build(const TrackingParticle& iData,
                                                        unsigned int tpIdx,
                                                        TEveElement& comp,
                                                        const FWViewContext*) {
  TEveRecTrack t;
  t.fBeta = 1.0;
  t.fP = TEveVector(iData.px(), iData.py(), iData.pz());
  t.fV = TEveVector(iData.vx(), iData.vy(), iData.vz());
  t.fSign = iData.charge();
 
  TEveTrack* track = new TEveTrack(&t, context().getTrackPropagator());
  if (t.fSign == 0)
    track->SetLineStyle(7);
 
  track->MakeTrack();
  setupAddElement(track, &comp);
  // printf("add track %d \n", tpIdx);
 
  if (m_assocList) {
    TEvePointSet* pointSet = new TEvePointSet;
    setupAddElement(pointSet, &comp);
 
    const FWGeometry* geom = item()->getGeom();
    float local[3];
    float localDir[3];
    float global[3] = {0.0, 0.0, 0.0};
    float globalDir[3] = {0.0, 0.0, 0.0};
 
    TrackingParticleRef tpr(tpch, tpIdx);
    std::pair<TrackingParticleRef, TrackPSimHitRef> clusterTPpairWithDummyTP(tpr, TrackPSimHitRef());
    auto range = std::equal_range(m_assocList->begin(),
                                  m_assocList->end(),
                                  clusterTPpairWithDummyTP,
                                  SimHitTPAssociationProducer::simHitTPAssociationListGreater);
    // printf("TrackingParticle[%d] P(%.1f, %.1f, %.1f) matches %d hits\n", tpIdx,iData.px(), iData.py(), iData.pz() ,(int)(range.second-range.first ));
 
    std::vector<const PSimHit*> phits;
    for (auto ri = range.first; ri != range.second; ++ri)
      phits.push_back(ri->second.get());
 
    std::sort(phits.begin(), phits.end(), [](const PSimHit* a, const PSimHit* b) { return a->tof() < b->tof(); });
    for (auto phi = phits.begin(); phi != phits.end(); ++phi) {
      const PSimHit* phit = *phi;
 
      local[0] = phit->localPosition().x();
      local[1] = phit->localPosition().y();
      local[2] = phit->localPosition().z();
 
      localDir[0] = phit->momentumAtEntry().x();
      localDir[1] = phit->momentumAtEntry().y();
      localDir[2] = phit->momentumAtEntry().z();
 
      geom->localToGlobal(phit->detUnitId(), local, global);
      geom->localToGlobal(phit->detUnitId(), localDir, globalDir, false);
      pointSet->SetNextPoint(global[0], global[1], global[2]);
 
      //printf("localP = (%f, %f, %f) globalP = (%f, %f, %f), loss = %f, tof =%f\n", localDir[0], localDir[1], localDir[2],
      //       globalDir[0], globalDir[1], globalDir[2],
      //       phit->energyLoss(), phit->tof());
      track->AddPathMark(TEvePathMark(TEvePathMark::kReference,
                                      TEveVector(global[0], global[1], global[2]),
                                      TEveVector(globalDir[0], globalDir[1], globalDir[2])));
    }
  }
}
 
REGISTER_FWPROXYBUILDER(FWTrackingParticleProxyBuilderFullFramework,
                        TrackingParticle,
                        "TrackingParticleWithPSimHits",
                        FWViewType::kAll3DBits | FWViewType::kAllRPZBits);