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 28 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 doPUWeights_, lumiWeights_, pileupSummaryInfoToken_, and trackingPlots::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 41 of file L1AnalysisEvent.cc.

{

}

Member Function Documentation

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

inline

Definition at line 43 of file L1AnalysisEvent.h.

References event_.

{return &event_;}

void L1Analysis::L1AnalysisEvent::Reset ( )

inline

Definition at line 42 of file L1AnalysisEvent.h.

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

Referenced by L1EventTreeProducer::analyze().

{event_.Reset();}

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

Definition at line 46 of file L1AnalysisEvent.cc.

References edm::HLTGlobalStatus::accept(), edm::EventBase::bunchCrossing(), L1Analysis::L1AnalysisEventDataFormat::bx, doPUWeights_, L1Analysis::L1AnalysisEventDataFormat::event, edm::EventID::event(), event_, edm::Event::eventAuxiliary(), edm::Event::getByLabel(), edm::Event::getByToken(), L1Analysis::L1AnalysisEventDataFormat::hlt, edm::EventBase::id(), edm::EventAuxiliary::isRealData(), edm::EDGetTokenT< T >::isUninitialized(), edm::HandleBase::isValid(), L1Analysis::L1AnalysisEventDataFormat::lumi, edm::EventBase::luminosityBlock(), lumiWeights_, maxAllowedWeight_, L1Analysis::L1AnalysisEventDataFormat::nPV, L1Analysis::L1AnalysisEventDataFormat::nPV_True, L1Analysis::L1AnalysisEventDataFormat::orbit, edm::EventBase::orbitNumber(), pileupSummaryInfoToken_, L1Analysis::L1AnalysisEventDataFormat::puWeight, L1Analysis::L1AnalysisEventDataFormat::run, edm::EventID::run(), edm::HLTGlobalStatus::size(), L1Analysis::L1AnalysisEventDataFormat::time, edm::EventBase::time(), edm::TriggerNames::triggerName(), edm::Event::triggerNames(), useAvgVtx_, edm::Timestamp::value(), mps_merge::weight, and edm::LumiReWeighting::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;

}

Member Data Documentation

bool L1Analysis::L1AnalysisEvent::doPUWeights_

private

Definition at line 49 of file L1AnalysisEvent.h.

Referenced by L1AnalysisEvent(), and Set().

L1Analysis::L1AnalysisEventDataFormat L1Analysis::L1AnalysisEvent::event_

private

Definition at line 57 of file L1AnalysisEvent.h.

Referenced by getData(), Reset(), and Set().

bool L1Analysis::L1AnalysisEvent::fillHLT_

private

Definition at line 48 of file L1AnalysisEvent.h.

edm::LumiReWeighting L1Analysis::L1AnalysisEvent::lumiWeights_

private

Definition at line 54 of file L1AnalysisEvent.h.

Referenced by L1AnalysisEvent(), and Set().

double L1Analysis::L1AnalysisEvent::maxAllowedWeight_

private

Definition at line 52 of file L1AnalysisEvent.h.

Referenced by Set().

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

private

Definition at line 56 of file L1AnalysisEvent.h.

Referenced by L1AnalysisEvent(), and Set().

bool L1Analysis::L1AnalysisEvent::useAvgVtx_

private

Definition at line 51 of file L1AnalysisEvent.h.

Referenced by Set().