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)
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_
bool	useAvgVtx_

Detailed Description

Definition at line 25 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
	)

Definition at line 10 of file L1AnalysisEvent.cc.

References doPUWeights_, and lumiWeights_.

                                           :
  fillHLT_(true),
  doPUWeights_(false),
  useAvgVtx_(useAvgVtx),
  maxAllowedWeight_(maxWeight),
  lumiWeights_()
{
 
  // 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 39 of file L1AnalysisEvent.cc.

{

}

Member Function Documentation

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

inline

Definition at line 39 of file L1AnalysisEvent.h.

References event_.

{return &event_;}

void L1Analysis::L1AnalysisEvent::Reset ( void  )

inline

Definition at line 38 of file L1AnalysisEvent.h.

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

Referenced by L1EventTreeProducer::analyze(), and L1NtupleProducer::analyzeEvent().

{event_.Reset();}

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

Definition at line 44 of file L1AnalysisEvent.cc.

References edm::EventBase::bunchCrossing(), 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(), and L1NtupleProducer::analyzeEvent().

{

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

  }

  event_.puWeight = weight;




}

Member Data Documentation

bool L1Analysis::L1AnalysisEvent::doPUWeights_

private

Definition at line 45 of file L1AnalysisEvent.h.

Referenced by L1AnalysisEvent().

L1Analysis::L1AnalysisEventDataFormat L1Analysis::L1AnalysisEvent::event_

private

Definition at line 52 of file L1AnalysisEvent.h.

Referenced by getData(), and Reset().

bool L1Analysis::L1AnalysisEvent::fillHLT_

private

Definition at line 44 of file L1AnalysisEvent.h.

edm::LumiReWeighting L1Analysis::L1AnalysisEvent::lumiWeights_

private

Definition at line 50 of file L1AnalysisEvent.h.

Referenced by L1AnalysisEvent().

double L1Analysis::L1AnalysisEvent::maxAllowedWeight_

private

Definition at line 48 of file L1AnalysisEvent.h.

bool L1Analysis::L1AnalysisEvent::useAvgVtx_

private

Definition at line 47 of file L1AnalysisEvent.h.