MuonToTrackingParticleAssociatorEDProducer.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    SimMuon/MCTruth
// Class:      MuonToTrackingParticleAssociatorEDProducer
// 
//
// Original Author:  Christopher Jones
//         Created:  Wed, 07 Jan 2015 21:30:14 GMT
//
//


// system include files
#include <memory>

// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/stream/EDProducer.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "SimMuon/MCTruth/interface/TrackerMuonHitExtractor.h"
#include "SimDataFormats/Associations/interface/MuonToTrackingParticleAssociator.h"
#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
#include "Geometry/Records/interface/TrackerTopologyRcd.h"
#include "MuonToTrackingParticleAssociatorByHitsImpl.h"

//
// class declaration
//
namespace {
  using TrackHitsCollection = MuonAssociatorByHitsHelper::TrackHitsCollection;

  class InputDumper {
  public:
    
    InputDumper(const edm::ParameterSet& conf, edm::ConsumesCollector && iC) :
      simtracksTag(conf.getParameter<edm::InputTag>("simtracksTag")),
      simtracksXFTag(conf.getParameter<edm::InputTag>("simtracksXFTag")),
      crossingframe(conf.getParameter<bool>("crossingframe")) 
    {
      if (crossingframe) {
        simtracksXFToken_ = iC.consumes<CrossingFrame<SimTrack> >(simtracksXFTag);
        simvertsXFToken_ = iC.consumes<CrossingFrame<SimVertex> >(simtracksXFTag);
      }
      else{
        simtracksToken_ = iC.consumes<edm::SimTrackContainer>(simtracksTag);
        simvertsToken_ = iC.consumes<edm::SimVertexContainer>(simtracksTag);
      }
    }

    void read(const edm::Event&);
    void dump(const TrackHitsCollection&, const TrackingParticleCollection&) const;

  private:
    edm::InputTag simtracksTag;
    edm::InputTag simtracksXFTag;
    edm::EDGetTokenT<CrossingFrame<SimTrack> > simtracksXFToken_;
    edm::EDGetTokenT<CrossingFrame<SimVertex> > simvertsXFToken_;
    edm::EDGetTokenT<edm::SimTrackContainer> simtracksToken_;
    edm::EDGetTokenT<edm::SimVertexContainer> simvertsToken_;

    edm::Handle<CrossingFrame<SimTrack>> simtracksXF_;
    edm::Handle<CrossingFrame<SimVertex>> simvertsXF_;
    edm::Handle<edm::SimTrackContainer> simtracks_;
    edm::Handle<edm::SimVertexContainer> simverts_;
    bool const crossingframe;
  };
  
  void
  InputDumper::read(const edm::Event& iEvent) {
    if(crossingframe) {
      iEvent.getByToken(simtracksXFToken_, simtracksXF_);
      iEvent.getByToken(simvertsXFToken_, simvertsXF_);
    }else {
      iEvent.getByToken(simtracksToken_, simtracks_);
      iEvent.getByToken(simvertsToken_, simverts_);
    }
  }

  void
  InputDumper::dump(const TrackHitsCollection& tC, const TrackingParticleCollection& tPC) const {
    using namespace std;
    // reco::Track collection
    edm::LogVerbatim("MuonToTrackingParticleAssociatorEDProducer")<<"\n"<<"reco::Track collection --- size = "<<tC.size();
 
    // TrackingParticle collection
    edm::LogVerbatim("MuonToTrackingParticleAssociatorEDProducer")<<"\n"<<"TrackingParticle collection --- size = "<<tPC.size();
    int j = 0;
    for(TrackingParticleCollection::const_iterator ITER=tPC.begin(); ITER!=tPC.end(); ITER++, j++) {
      edm::LogVerbatim("MuonToTrackingParticleAssociatorEDProducer")
        <<"TrackingParticle "<<j<<", q = "<<ITER->charge()<<", p = "<<ITER->p()
        <<", pT = "<<ITER->pt()<<", eta = "<<ITER->eta()<<", phi = "<<ITER->phi();
       
      edm::LogVerbatim("MuonToTrackingParticleAssociatorEDProducer")
        <<"\t pdg code = "<<ITER->pdgId()<<", made of "<<ITER->numberOfHits()<<" PSimHit"
        <<" (in "<<ITER->numberOfTrackerLayers()<<" layers)"
        <<" from "<<ITER->g4Tracks().size()<<" SimTrack:";
      for (TrackingParticle::g4t_iterator g4T=ITER->g4Track_begin(); g4T!=ITER->g4Track_end(); g4T++) {
        edm::LogVerbatim("MuonToTrackingParticleAssociatorEDProducer")
          <<"\t\t Id:"<<g4T->trackId()<<"/Evt:("<<g4T->eventId().event()<<","<<g4T->eventId().bunchCrossing()<<")";
      }    
    }
 
    if (crossingframe) {
      std:: unique_ptr<MixCollection<SimTrack> > SimTk( new MixCollection<SimTrack>(simtracksXF_.product()) );
      edm::LogVerbatim("MuonToTrackingParticleAssociatorEDProducer")
    <<"\n"<<"CrossingFrame<SimTrack> collection with InputTag = "<<simtracksXFTag<<" has size = "<<SimTk->size();
      int k = 0;
      for (MixCollection<SimTrack>::MixItr ITER=SimTk->begin(); ITER!=SimTk->end(); ITER++, k++) {
        edm::LogVerbatim("MuonToTrackingParticleAssociatorEDProducer")
          <<"SimTrack "<<k
          <<" - Id:"<<ITER->trackId()<<"/Evt:("<<ITER->eventId().event()<<","<<ITER->eventId().bunchCrossing()<<")"
          <<", q = "<<ITER->charge()<<", p = "<<ITER->momentum().P()
          <<", pT = "<<ITER->momentum().Pt()<<", eta = "<<ITER->momentum().Eta()<<", phi = "<<ITER->momentum().Phi()
      <<"\n\t pdgId = "<<ITER->type()<<", Vertex index = "<<ITER->vertIndex()
      <<", Gen Particle index = "<< (ITER->genpartIndex()>0 ? ITER->genpartIndex()-1 : ITER->genpartIndex()) <<endl;
      }

      std::unique_ptr<MixCollection<SimVertex> > SimVtx( new MixCollection<SimVertex>(simvertsXF_.product()) );
      edm::LogVerbatim("MuonToTrackingParticleAssociatorEDProducer")
    <<"\n"<<"CrossingFrame<SimVertex> collection with InputTag = "<<simtracksXFTag<<" has size = "<<SimVtx->size();
      int kv = 0;
      for (MixCollection<SimVertex>::MixItr VITER=SimVtx->begin(); VITER!=SimVtx->end(); VITER++, kv++){
        edm::LogVerbatim("MuonToTrackingParticleAssociatorEDProducer")
          <<"SimVertex "<<kv
      <<" - Id:"<<VITER->vertexId()<<", position = "<<VITER->position()<<", parent SimTrack Id = "<<VITER->parentIndex()
      <<", processType = "<<VITER->processType();
      }
    }
    else {
      const edm::SimTrackContainer simTC = *(simtracks_.product());
      edm::LogVerbatim("MuonToTrackingParticleAssociatorEDProducer")<<"\n"<<"SimTrack collection with InputTag = "<<simtracksTag
                                              <<" has size = "<<simTC.size()<<endl;
      int k = 0;
      for(edm::SimTrackContainer::const_iterator ITER=simTC.begin(); ITER!=simTC.end(); ITER++, k++){
        edm::LogVerbatim("MuonToTrackingParticleAssociatorEDProducer")
          <<"SimTrack "<<k
          <<" - Id:"<<ITER->trackId()<<"/Evt:("<<ITER->eventId().event()<<","<<ITER->eventId().bunchCrossing()<<")"
          <<", q = "<<ITER->charge()<<", p = "<<ITER->momentum().P()
          <<", pT = "<<ITER->momentum().Pt()<<", eta = "<<ITER->momentum().Eta()<<", phi = "<<ITER->momentum().Phi()
      <<"\n\t pdgId = "<<ITER->type()<<", Vertex index = "<<ITER->vertIndex()
      <<", Gen Particle index = "<< (ITER->genpartIndex()>0 ? ITER->genpartIndex()-1 : ITER->genpartIndex()) <<endl;
      }
      const edm::SimVertexContainer simVC = *(simverts_.product());
      edm::LogVerbatim("MuonToTrackingParticleAssociatorEDProducer")<<"\n"<<"SimVertex collection with InputTag = "<<"g4SimHits"
                                              <<" has size = "<<simVC.size()<<endl;
      int kv = 0;
      for (edm::SimVertexContainer::const_iterator VITER=simVC.begin(); VITER!=simVC.end(); VITER++, kv++){
        edm::LogVerbatim("MuonToTrackingParticleAssociatorEDProducer")
          <<"SimVertex "<<kv
      <<" - Id:"<<VITER->vertexId()<<", position = "<<VITER->position()<<", parent SimTrack Id = "<<VITER->parentIndex()
      <<", processType = "<<VITER->processType();
      }
    }

  }

}


class MuonToTrackingParticleAssociatorEDProducer : public edm::stream::EDProducer<> {
public:
  explicit MuonToTrackingParticleAssociatorEDProducer(const edm::ParameterSet&);
  ~MuonToTrackingParticleAssociatorEDProducer() override;
  
  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);
  
private:
  void produce(edm::Event&, const edm::EventSetup&) override;
  
  // ----------member data ---------------------------
  edm::ParameterSet const config_;
  MuonAssociatorByHitsHelper helper_;
  TrackerHitAssociator::Config trackerHitAssociatorConfig_;
  TrackerMuonHitExtractor hitExtractor_;

  std::unique_ptr<RPCHitAssociator> rpctruth_;
  std::unique_ptr<GEMHitAssociator> gemtruth_;
  std::unique_ptr<DTHitAssociator> dttruth_;
  std::unique_ptr<CSCHitAssociator> csctruth_;
  std::unique_ptr<TrackerHitAssociator> trackertruth_;
  std::unique_ptr<InputDumper> diagnostics_;
};

//
// constants, enums and typedefs
//


//
// static data member definitions
//

//
// constructors and destructor
//
MuonToTrackingParticleAssociatorEDProducer::MuonToTrackingParticleAssociatorEDProducer(const edm::ParameterSet& iConfig):
  config_(iConfig),
  helper_(iConfig),
  trackerHitAssociatorConfig_(iConfig,consumesCollector()),
  hitExtractor_(iConfig,consumesCollector())
{
   //register your products
   produces<reco::MuonToTrackingParticleAssociator>();

   //hack for consumes
   RPCHitAssociator rpctruth(iConfig,consumesCollector());
   GEMHitAssociator gemtruth(iConfig,consumesCollector());
   DTHitAssociator dttruth(iConfig,consumesCollector());
   CSCHitAssociator cscruth(iConfig,consumesCollector());

  if( iConfig.getUntrackedParameter<bool>("dumpInputCollections") ) {
    diagnostics_.reset( new InputDumper(iConfig, consumesCollector()) );
  }
}


MuonToTrackingParticleAssociatorEDProducer::~MuonToTrackingParticleAssociatorEDProducer()
{
 
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)

}


//
// member functions
//

// ------------ method called to produce the data  ------------
void
MuonToTrackingParticleAssociatorEDProducer::produce(edm::Event& iEvent, const edm::EventSetup& iSetup)
{
   using namespace edm;

   hitExtractor_.init(iEvent, iSetup);

   //Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHand;
   iSetup.get<TrackerTopologyRcd>().get(tTopoHand);
   const TrackerTopology *tTopo=tTopoHand.product();
   
   bool printRtS = true;
   
   //NOTE: This assumes that produce will not be called until the edm::Event used in the previous call
   // has been deleted. This is true for now. In the future, we may have to have the resources own
   // the memory.

   // Tracker hit association  
   trackertruth_.reset(new TrackerHitAssociator(iEvent, trackerHitAssociatorConfig_));
   // CSC hit association
   csctruth_.reset(new CSCHitAssociator(iEvent,iSetup,config_));
   // DT hit association
   printRtS = false;
   dttruth_.reset(new DTHitAssociator(iEvent,iSetup,config_,printRtS));
   // RPC hit association
   rpctruth_.reset( new RPCHitAssociator(iEvent,iSetup,config_) );
   // GEM hit association
   gemtruth_.reset( new GEMHitAssociator(iEvent,iSetup,config_) );
   
   MuonAssociatorByHitsHelper::Resources resources = {tTopo, trackertruth_.get(), csctruth_.get(), dttruth_.get(), rpctruth_.get(), gemtruth_.get()};

   if(diagnostics_) {
     diagnostics_->read(iEvent);
     resources.diagnostics_ = [this](const TrackHitsCollection& hC, const TrackingParticleCollection& pC) {
       diagnostics_->dump(hC,pC);
     };
   }

   std::unique_ptr<reco::MuonToTrackingParticleAssociatorBaseImpl> impl{ 
     new MuonToTrackingParticleAssociatorByHitsImpl(hitExtractor_,resources, &helper_) }; 
   std::unique_ptr<reco::MuonToTrackingParticleAssociator> toPut( new reco::MuonToTrackingParticleAssociator(std::move(impl)));
   iEvent.put(std::move(toPut));
}

 
// ------------ method fills 'descriptions' with the allowed parameters for the module  ------------
void
MuonToTrackingParticleAssociatorEDProducer::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
  //The following says we do not know what parameters are allowed so do no validation
  // Please change this to state exactly what you do use, even if it is no parameters
  edm::ParameterSetDescription desc;
  desc.setUnknown();
  descriptions.addDefault(desc);
}

//define this as a plug-in
DEFINE_FWK_MODULE(MuonToTrackingParticleAssociatorEDProducer);