L1Analysis::L1AnalysisEvent Class Reference

#include <L1AnalysisEvent.h>

Public Member Functions
L1AnalysisEventDataFormat *	getData ()
	L1AnalysisEvent (std::string puMCFile, std::string puMCHist, std::string puDataFile, std::string puDataHist, bool useAvgVtx, double maxWeight, edm::ConsumesCollector &&)
void	Reset ()
void	Set (const edm::Event &e, const edm::EDGetTokenT< edm::TriggerResults > &hlt_)
	~L1AnalysisEvent ()

Private Attributes
bool	doPUWeights_
L1Analysis::L1AnalysisEventDataFormat	event_
bool	fillHLT_
edm::LumiReWeighting	lumiWeights_
double	maxAllowedWeight_
edm::EDGetTokenT< std::vector< PileupSummaryInfo > >	pileupSummaryInfoToken_
bool	useAvgVtx_

Detailed Description

Definition at line 27 of file L1AnalysisEvent.h.

Constructor & Destructor Documentation

L1AnalysisEvent()

L1Analysis::L1AnalysisEvent::L1AnalysisEvent	(	std::string	puMCFile,
		std::string	puMCHist,
		std::string	puDataFile,
		std::string	puDataHist,
		bool	useAvgVtx,
		double	maxWeight,
		edm::ConsumesCollector &&	iConsumes
	)

Definition at line 10 of file L1AnalysisEvent.cc.

References visDQMUpload::buf, doPUWeights_, HLT_2022v12_cff::InputTag, lumiWeights_, pileupSummaryInfoToken_, l1EventTree_cfi::puDataFile, l1EventTree_cfi::puDataHist, l1EventTree_cfi::puMCFile, l1EventTree_cfi::puMCHist, and edm_modernize_messagelogger::stat.

    : fillHLT_(true), doPUWeights_(false), useAvgVtx_(useAvgVtx), maxAllowedWeight_(maxWeight), lumiWeights_() {
  pileupSummaryInfoToken_ = iConsumes.consumes<std::vector<PileupSummaryInfo>>(edm::InputTag("addPileupInfo"));
  // check PU files exists, and reweight if they do
  struct stat buf;
  if ((stat(puMCFile.c_str(), &buf) != -1) && (stat(puDataFile.c_str(), &buf) != -1)) {
    lumiWeights_ = edm::LumiReWeighting(puMCFile, puDataFile, puMCHist, puDataHist);
    doPUWeights_ = true;
  } else {
    edm::LogWarning("L1Prompt") << "No PU reweighting inputs - not going to calculate weights" << std::endl;
  }
}

~L1AnalysisEvent()

L1Analysis::L1AnalysisEvent::~L1AnalysisEvent

(

)

Definition at line 29 of file L1AnalysisEvent.cc.

{}

Member Function Documentation

getData()

L1AnalysisEventDataFormat* L1Analysis::L1AnalysisEvent::getData ( void  )

inline

Definition at line 41 of file L1AnalysisEvent.h.

References event_.

{ return &event_; }

Reset()

void L1Analysis::L1AnalysisEvent::Reset ( )

inline

Definition at line 40 of file L1AnalysisEvent.h.

References event_, and L1Analysis::L1AnalysisEventDataFormat::Reset().

{ event_.Reset(); }

Set()

void L1Analysis::L1AnalysisEvent::Set	(	const edm::Event &	e,
		const edm::EDGetTokenT< edm::TriggerResults > &	hlt_
	)

Definition at line 31 of file L1AnalysisEvent.cc.

References simKBmtfDigis_cfi::bx, MillePedeFileConverter_cfg::e, HLTBitAnalyser_cfi::hltresults, edm::EDGetTokenT< T >::isUninitialized(), edm::HandleBase::isValid(), edm::TriggerNames::triggerName(), cscTnPEfficiencyTask_cfi::trigName, and mps_merge::weight.

                                                                                                    {
  event_.run = e.id().run();
  event_.event = e.id().event();
  event_.time = e.time().value();
  event_.bx = e.bunchCrossing();  //overwritten by EVM info until fixed by fw
  event_.lumi = e.luminosityBlock();
  event_.orbit = e.orbitNumber();  //overwritten by EVM info until fixed by fw

  if (!hlt_.isUninitialized()) {
    edm::Handle<edm::TriggerResults> hltresults;
    e.getByToken(hlt_, hltresults);
    const edm::TriggerNames& TrigNames_ = e.triggerNames(*hltresults);
    const int ntrigs = hltresults->size();

    for (int itr = 0; itr < ntrigs; itr++) {
      TString trigName = TrigNames_.triggerName(itr);
      if (!hltresults->accept(itr))
        continue;
      event_.hlt.push_back(trigName);
    }
  }

  // do PU re-weighting for MC only
  double weight = 1.;

  if (doPUWeights_ && (!e.eventAuxiliary().isRealData())) {
    edm::Handle<std::vector<PileupSummaryInfo>> puInfo;
    e.getByLabel(edm::InputTag("addPileupInfo"), puInfo);

    if (puInfo.isValid()) {
      std::vector<PileupSummaryInfo>::const_iterator pvi;

      float npv = -1;
      for (pvi = puInfo->begin(); pvi != puInfo->end(); ++pvi) {
        int bx = pvi->getBunchCrossing();

        if (bx == 0) {
          npv = useAvgVtx_ ? pvi->getTrueNumInteractions() : pvi->getPU_NumInteractions();
          continue;
        }
      }

      weight = lumiWeights_.weight(npv);
      if (maxAllowedWeight_ > 0. && weight > maxAllowedWeight_)
        weight = maxAllowedWeight_;
    }
  }

  if (!e.eventAuxiliary().isRealData()) {
    edm::Handle<std::vector<PileupSummaryInfo>> puInfo;
    e.getByToken(pileupSummaryInfoToken_, puInfo);
    if (puInfo.isValid()) {
      for (std::vector<PileupSummaryInfo>::const_iterator pvi = puInfo->begin(); pvi != puInfo->end(); pvi++) {
        int bx = pvi->getBunchCrossing();
        if (bx == 0) {
          event_.nPV = pvi->getPU_NumInteractions();
          event_.nPV_True = pvi->getTrueNumInteractions();
        }
      }
    }
  }

  event_.puWeight = weight;
}

Member Data Documentation

doPUWeights_

bool L1Analysis::L1AnalysisEvent::doPUWeights_

private

Definition at line 47 of file L1AnalysisEvent.h.

Referenced by L1AnalysisEvent().

event_

L1Analysis::L1AnalysisEventDataFormat L1Analysis::L1AnalysisEvent::event_

private

Definition at line 55 of file L1AnalysisEvent.h.

Referenced by getData(), and Reset().

fillHLT_

bool L1Analysis::L1AnalysisEvent::fillHLT_

private

Definition at line 46 of file L1AnalysisEvent.h.

lumiWeights_

edm::LumiReWeighting L1Analysis::L1AnalysisEvent::lumiWeights_

private

Definition at line 52 of file L1AnalysisEvent.h.

Referenced by L1AnalysisEvent().

maxAllowedWeight_

double L1Analysis::L1AnalysisEvent::maxAllowedWeight_

private

Definition at line 50 of file L1AnalysisEvent.h.

pileupSummaryInfoToken_

edm::EDGetTokenT<std::vector<PileupSummaryInfo> > L1Analysis::L1AnalysisEvent::pileupSummaryInfoToken_

private

Definition at line 54 of file L1AnalysisEvent.h.

Referenced by L1AnalysisEvent().

useAvgVtx_

bool L1Analysis::L1AnalysisEvent::useAvgVtx_

private

Definition at line 49 of file L1AnalysisEvent.h.