#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_, HLT_25ns10e33_v2_cff::InputTag, lumiWeights_, and pileupSummaryInfoToken_.

                                                             :
   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_.

43 {return &event_;}

L1Analysis::L1AnalysisEvent::event_

L1Analysis::L1AnalysisEventDataFormat event_

Definition: L1AnalysisEvent.h:57

void L1Analysis::L1AnalysisEvent::Reset ( void )

inline

Definition at line 42 of file L1AnalysisEvent.h.

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

Referenced by L1EventTreeProducer::analyze().

42 {event_.Reset();}

L1Analysis::L1AnalysisEvent::event_

L1Analysis::L1AnalysisEventDataFormat event_

Definition: L1AnalysisEvent.h:57

L1Analysis::L1AnalysisEventDataFormat::Reset

void Reset()

Definition: L1AnalysisEventDataFormat.h:22

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

Definition at line 46 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().

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