PatternOptimizerBase.h

Go to the documentation of this file.

/*
 * PatternOptimizerBase.h
 *
 *  Created on: Oct 17, 2018
 *      Author: kbunkow
 */

#ifndef L1T_OmtfP1_PATTERNOPTIMIZERBASE_H_
#define L1T_OmtfP1_PATTERNOPTIMIZERBASE_H_

#include "L1Trigger/L1TMuonOverlapPhase1/interface/Tools/EmulationObserverBase.h"
#include "L1Trigger/L1TMuonOverlapPhase1/interface/Omtf/GoldenPatternWithStat.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "SimDataFormats/Track/interface/SimTrackContainer.h"

#include "TH1I.h"
#include "TH2I.h"
#include "TFile.h"

#include <functional>

class PatternOptimizerBase : public EmulationObserverBase {
public:
  static double vxMuRate(double pt_GeV);
  static double vxIntegMuRate(double pt_GeV, double dpt, double etaFrom, double etaTo);

  PatternOptimizerBase(const edm::ParameterSet& edmCfg,
                       const OMTFConfiguration* omtfConfig,
                       GoldenPatternVec<GoldenPatternWithStat>& gps);

  ~PatternOptimizerBase() override;

  void observeEventEnd(const edm::Event& iEvent,
                       std::unique_ptr<l1t::RegionalMuonCandBxCollection>& finalCandidates) override;

  void endJob() override;

protected:
  void savePatternsInRoot(std::string rootFileName);

  void printPatterns();

  virtual double getEventRateWeight(double pt) { return 1; }

  virtual void saveHists(TFile& outfile) {}

  GoldenPatternVec<GoldenPatternWithStat>& goldenPatterns;

  TH1I* simMuPt = nullptr;
  TH1I* simMuFoundByOmtfPt = nullptr;
  TH1I* simMuEta = nullptr;

  TH1I* candEta = nullptr;

  TH1F* simMuPtSpectrum = nullptr;
  TH2I* simMuPtVsDispl = nullptr;
  TH2I* simMuPtVsRho = nullptr;

  bool writeLayerStat = false;
};

#endif /* L1T_OmtfP1_PATTERNOPTIMIZERBASE_H_ */