#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

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_FULL_cff::InputTag, lumiWeights_, pileupSummaryInfoToken_, 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;
   }
 }

L1Analysis::L1AnalysisEvent::~L1AnalysisEvent ( )

Definition at line 29 of file L1AnalysisEvent.cc.

29 {}

Member Function Documentation

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

inline

Definition at line 41 of file L1AnalysisEvent.h.

References event_.

41 { return &event_; }

L1Analysis::L1AnalysisEvent::event_

L1Analysis::L1AnalysisEventDataFormat event_

Definition: L1AnalysisEvent.h:55

void L1Analysis::L1AnalysisEvent::Reset ( )

inline

Definition at line 40 of file L1AnalysisEvent.h.

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

Referenced by L1EventTreeProducer::analyze().

40 { event_.Reset(); }

L1Analysis::L1AnalysisEvent::event_

L1Analysis::L1AnalysisEventDataFormat event_

Definition: L1AnalysisEvent.h:55

L1Analysis::L1AnalysisEventDataFormat::Reset

void Reset()

Definition: L1AnalysisEventDataFormat.h:20

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

Definition at line 31 of file L1AnalysisEvent.cc.

References edm::EventBase::bunchCrossing(), makePileupJSON::bx, edm::EventID::event(), event_(), edm::Event::eventAuxiliary(), edm::Event::getByLabel(), edm::Event::getByToken(), edm::EventBase::id(), edm::EventAuxiliary::isRealData(), edm::EDGetTokenT< T >::isUninitialized(), edm::HandleBase::isValid(), edm::EventBase::luminosityBlock(), edm::EventBase::orbitNumber(), edm::EventID::run(), edm::EventBase::time(), edm::TriggerNames::triggerName(), edm::Event::triggerNames(), edm::Timestamp::value(), and histoStyle::weight.

Referenced by L1EventTreeProducer::analyze().

                                                                                                     {
   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;
 }