CaloCleaner.h

Go to the documentation of this file.

#ifndef TauAnalysis_MCEmbeddingTools_CaloCleaner_H
#define TauAnalysis_MCEmbeddingTools_CaloCleaner_H

#include "FWCore/Framework/interface/stream/EDProducer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/PatCandidates/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonEnergy.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "TrackingTools/TrackAssociator/interface/TrackAssociatorParameters.h"
#include "TrackingTools/TrackAssociator/interface/TrackDetectorAssociator.h"
#include "TrackingTools/Records/interface/TrackingComponentsRecord.h"

#include "DataFormats/Common/interface/SortedCollection.h"
#include "DataFormats/EcalRecHit/interface/EcalRecHit.h"

#include <string>
#include <iostream>
#include <map>

template <typename T>
class CaloCleaner : public edm::stream::EDProducer<> {
public:
  explicit CaloCleaner(const edm::ParameterSet&);
  ~CaloCleaner() override;

private:
  void produce(edm::Event&, const edm::EventSetup&) override;

  typedef edm::SortedCollection<T> RecHitCollection;

  const edm::EDGetTokenT<edm::View<pat::Muon> > mu_input_;

  std::map<std::string, edm::EDGetTokenT<RecHitCollection> > inputs_;

  TrackDetectorAssociator trackAssociator_;
  TrackAssociatorParameters parameters_;

  bool is_preshower_;
  void fill_correction_map(TrackDetMatchInfo*, std::map<uint32_t, float>*);
};

template <typename T>
CaloCleaner<T>::CaloCleaner(const edm::ParameterSet& iConfig)
    : mu_input_(consumes<edm::View<pat::Muon> >(iConfig.getParameter<edm::InputTag>("MuonCollection"))) {
  std::vector<edm::InputTag> inCollections = iConfig.getParameter<std::vector<edm::InputTag> >("oldCollection");
  for (auto inCollection : inCollections) {
    inputs_[inCollection.instance()] = consumes<RecHitCollection>(inCollection);
    produces<RecHitCollection>(inCollection.instance());
  }

  is_preshower_ = iConfig.getUntrackedParameter<bool>("is_preshower", false);
  edm::ParameterSet parameters = iConfig.getParameter<edm::ParameterSet>("TrackAssociatorParameters");
  edm::ConsumesCollector iC = consumesCollector();
  parameters_.loadParameters(parameters, iC);
  //trackAssociator_.useDefaultPropagator();
}

template <typename T>
CaloCleaner<T>::~CaloCleaner() {
  // nothing to be done yet...
}

template <typename T>
void CaloCleaner<T>::produce(edm::Event& iEvent, const edm::EventSetup& iSetup) {
  edm::ESHandle<Propagator> propagator;
  iSetup.get<TrackingComponentsRecord>().get("SteppingHelixPropagatorAny", propagator);
  trackAssociator_.setPropagator(propagator.product());

  edm::Handle<edm::View<pat::Muon> > muonHandle;
  iEvent.getByToken(mu_input_, muonHandle);
  edm::View<pat::Muon> muons = *muonHandle;

  std::map<uint32_t, float> correction_map;

  // Fill the correction map
  for (edm::View<pat::Muon>::const_iterator iMuon = muons.begin(); iMuon != muons.end(); ++iMuon) {
    // get the basic informaiton like fill reco mouon does
    //     RecoMuon/MuonIdentification/plugins/MuonIdProducer.cc
    const reco::Track* track = nullptr;
    if (iMuon->track().isNonnull())
      track = iMuon->track().get();
    else if (iMuon->standAloneMuon().isNonnull())
      track = iMuon->standAloneMuon().get();
    else
      throw cms::Exception("FatalError")
          << "Failed to fill muon id information for a muon with undefined references to tracks";
    TrackDetMatchInfo info =
        trackAssociator_.associate(iEvent, iSetup, *track, parameters_, TrackDetectorAssociator::Any);
    fill_correction_map(&info, &correction_map);
  }

  // Copy the old collection and correct if necessary
  for (auto input_ : inputs_) {
    std::unique_ptr<RecHitCollection> recHitCollection_output(new RecHitCollection());
    edm::Handle<RecHitCollection> recHitCollection;
    // iEvent.getByToken(input_.second[0], recHitCollection);
    iEvent.getByToken(input_.second, recHitCollection);
    for (typename RecHitCollection::const_iterator recHit = recHitCollection->begin();
         recHit != recHitCollection->end();
         ++recHit) {
      if (correction_map[recHit->detid().rawId()] > 0) {
        float new_energy = recHit->energy() - correction_map[recHit->detid().rawId()];
        if (new_energy <= 0)
          continue;  // Do not save empty Hits
        T newRecHit(*recHit);
        newRecHit.setEnergy(new_energy);
        recHitCollection_output->push_back(newRecHit);
      } else {
        recHitCollection_output->push_back(*recHit);
      }
      /* For the inveted collection   
     if (correction_map[recHit->detid().rawId()] > 0){
    float new_energy =   correction_map[recHit->detid().rawId()];
    if (new_energy < 0) new_energy =0;
    T newRecHit(*recHit);
    newRecHit.setEnergy(new_energy); 
    recHitCollection_output->push_back(newRecHit);
      }*/
    }
    // Save the new collection
    recHitCollection_output->sort();
    iEvent.put(std::move(recHitCollection_output), input_.first);
  }
}
#endif