#include <LumiReWeighting.h>

Public Member Functions
	LumiReWeighting ()

	LumiReWeighting (const std::vector< float > &MC_distr, const std::vector< float > &Lumi_distr, const edm::InputTag &PileupSumInfoInputTag=edm::InputTag("addPileupInfo"))

	LumiReWeighting (std::string generatedFile, std::string dataFile, std::string GenHistName="pileup", std::string DataHistName="pileup", const edm::InputTag &PileupSumInfoInputTag=edm::InputTag("addPileupInfo"))

double	weight (const edm::EventBase &e)

double	weight (float npv)

double	weight (int npv)

double	weightOOT (const edm::EventBase &e)

Protected Attributes
std::shared_ptr< TH1 >	Data_distr_

std::shared_ptr< TFile >	dataFile_

std::string	dataFileName_

std::string	DataHistName_

bool	FirstWarning_

std::shared_ptr< TFile >	generatedFile_

std::string	generatedFileName_

std::string	GenHistName_

std::shared_ptr< TH1 >	MC_distr_

int	OldLumiSection_

edm::InputTag	pileupSumInfoTag_

std::shared_ptr< TH1 >	weights_

Detailed Description

Definition at line 105 of file LumiReWeighting.h.

Constructor & Destructor Documentation

◆ LumiReWeighting() [1/3]

LumiReWeighting::LumiReWeighting	(	std::string	generatedFile,
		std::string	dataFile,
		std::string	GenHistName = `"pileup"`,
		std::string	DataHistName = `"pileup"`,
		const edm::InputTag &	PileupSumInfoInputTag = `edm::InputTag("addPileupInfo")`
	)

Definition at line 39 of file LumiReWeighting.cc.

     : generatedFileName_(generatedFile),
       dataFileName_(dataFile),
       GenHistName_(GenHistName),
       DataHistName_(DataHistName),
       pileupSumInfoTag_(PileupSumInfoInputTag) {
   generatedFile_ = std::make_shared<TFile>(generatedFileName_.c_str());  //MC distribution
   dataFile_ = std::make_shared<TFile>(dataFileName_.c_str());            //Data distribution
  
   Data_distr_ = std::shared_ptr<TH1>((static_cast<TH1*>(dataFile_->Get(DataHistName_.c_str())->Clone())));
   MC_distr_ = std::shared_ptr<TH1>((static_cast<TH1*>(generatedFile_->Get(GenHistName_.c_str())->Clone())));
  
   // MC * data/MC = data, so the weights are data/MC:
  
   // normalize both histograms first
  
   Data_distr_->Scale(1.0 / Data_distr_->Integral());
   MC_distr_->Scale(1.0 / MC_distr_->Integral());
  
   weights_ = std::shared_ptr<TH1>(static_cast<TH1*>(Data_distr_->Clone()));
  
   weights_->SetName("lumiWeights");
  
   TH1* den = dynamic_cast<TH1*>(MC_distr_->Clone());
  
   //den->Scale(1.0/ den->Integral());
  
   weights_->Divide(den);  // so now the average weight should be 1.0
  
   std::cout << " Lumi/Pileup Reweighting: Computed Weights per In-Time Nint " << std::endl;
  
   int NBins = weights_->GetNbinsX();
  
   for (int ibin = 1; ibin < NBins + 1; ++ibin) {
     std::cout << "   " << ibin - 1 << " " << weights_->GetBinContent(ibin) << std::endl;
   }
  
   FirstWarning_ = true;
   OldLumiSection_ = -1;
 }

References gather_cfg::cout, Data_distr_, dataFile_, dataFileName_, DataHistName_, FirstWarning_, generatedFile_, generatedFileName_, GenHistName_, MC_distr_, SiStripSourceConfigP5_cff::NBins, OldLumiSection_, and weights_.

◆ LumiReWeighting() [2/3]

LumiReWeighting::LumiReWeighting	(	const std::vector< float > &	MC_distr,
		const std::vector< float > &	Lumi_distr,
		const edm::InputTag &	PileupSumInfoInputTag = `edm::InputTag("addPileupInfo")`
	)

Definition at line 84 of file LumiReWeighting.cc.

     : pileupSumInfoTag_(PileupSumInfoInputTag) {
   // no histograms for input: use vectors
  
   // now, make histograms out of them:
  
   // first, check they are the same size...
  
   if (MC_distr.size() != Lumi_distr.size()) {
     std::cerr << "ERROR: LumiReWeighting: input vectors have different sizes. Quitting... \n";
     return;
   }
  
   Int_t NBins = MC_distr.size();
  
   MC_distr_ = std::shared_ptr<TH1>(new TH1F("MC_distr", "MC dist", NBins, -0.5, float(NBins) - 0.5));
   Data_distr_ = std::shared_ptr<TH1>(new TH1F("Data_distr", "Data dist", NBins, -0.5, float(NBins) - 0.5));
  
   weights_ = std::shared_ptr<TH1>(new TH1F("luminumer", "luminumer", NBins, -0.5, float(NBins) - 0.5));
   TH1* den = new TH1F("lumidenom", "lumidenom", NBins, -0.5, float(NBins) - 0.5);
  
   for (int ibin = 1; ibin < NBins + 1; ++ibin) {
     weights_->SetBinContent(ibin, Lumi_distr[ibin - 1]);
     Data_distr_->SetBinContent(ibin, Lumi_distr[ibin - 1]);
     den->SetBinContent(ibin, MC_distr[ibin - 1]);
     MC_distr_->SetBinContent(ibin, MC_distr[ibin - 1]);
   }
  
   // check integrals, make sure things are normalized
  
   float deltaH = weights_->Integral();
   if (fabs(1.0 - deltaH) > 0.02) {  //*OOPS*...
     weights_->Scale(1.0 / weights_->Integral());
     Data_distr_->Scale(1.0 / Data_distr_->Integral());
   }
   float deltaMC = den->Integral();
   if (fabs(1.0 - deltaMC) > 0.02) {
     den->Scale(1.0 / den->Integral());
     MC_distr_->Scale(1.0 / MC_distr_->Integral());
   }
  
   weights_->Divide(den);  // so now the average weight should be 1.0
  
   std::cout << " Lumi/Pileup Reweighting: Computed Weights per In-Time Nint " << std::endl;
  
   for (int ibin = 1; ibin < NBins + 1; ++ibin) {
     std::cout << "   " << ibin - 1 << " " << weights_->GetBinContent(ibin) << std::endl;
   }
  
   FirstWarning_ = true;
   OldLumiSection_ = -1;
 }

References beam_dqm_sourceclient-live_cfg::cerr, gather_cfg::cout, Data_distr_, FirstWarning_, MC_distr_, SiStripSourceConfigP5_cff::NBins, OldLumiSection_, and weights_.

◆ LumiReWeighting() [3/3]

edm::LumiReWeighting::LumiReWeighting ( )

inline

Definition at line 117 of file LumiReWeighting.h.

117 {};

Member Function Documentation

◆ weight() [1/3]

double LumiReWeighting::weight ( const edm::EventBase & e )

Definition at line 152 of file LumiReWeighting.cc.

                                                     {
   if (FirstWarning_) {
     e.processHistory();
     //    SetFirstFalse();
     FirstWarning_ = false;
   }
   Handle<std::vector<PileupSummaryInfo> > PupInfo;
   e.getByLabel(pileupSumInfoTag_, PupInfo);
  
   std::vector<PileupSummaryInfo>::const_iterator PVI;
  
   int npv = -1;
   for (PVI = PupInfo->begin(); PVI != PupInfo->end(); ++PVI) {
     int BX = PVI->getBunchCrossing();
  
     if (BX == 0) {
       npv = PVI->getPU_NumInteractions();
       continue;
     }
   }
  
   if (npv < 0)
     std::cerr << " no in-time beam crossing found\n! ";
  
   return weight(npv);
 }

References L1TStage2uGTEmulatorClient_cff::BX, beam_dqm_sourceclient-live_cfg::cerr, MillePedeFileConverter_cfg::e, FirstWarning_, pileupSumInfoTag_, and weight().

◆ weight() [2/3]

double LumiReWeighting::weight ( float npv )

Definition at line 144 of file LumiReWeighting.cc.

                                         {
   int bin = weights_->GetXaxis()->FindBin(npv);
   return weights_->GetBinContent(bin);
 }

References newFWLiteAna::bin, and weights_.

◆ weight() [3/3]

double LumiReWeighting::weight ( int npv )

Definition at line 139 of file LumiReWeighting.cc.

                                       {
   int bin = weights_->GetXaxis()->FindBin(npv);
   return weights_->GetBinContent(bin);
 }

References newFWLiteAna::bin, and weights_.

Referenced by weight().

◆ weightOOT()

double LumiReWeighting::weightOOT ( const edm::EventBase & e )

Definition at line 183 of file LumiReWeighting.cc.

                                                        {
   //int Run = e.run();
   int LumiSection = e.luminosityBlock();
   // do some caching here, attempt to catch file boundaries
  
   if (LumiSection != OldLumiSection_) {
     e.processHistory();  // keep the function call
     OldLumiSection_ = LumiSection;
   }
   // find the pileup summary information
  
   Handle<std::vector<PileupSummaryInfo> > PupInfo;
   e.getByLabel(pileupSumInfoTag_, PupInfo);
  
   std::vector<PileupSummaryInfo>::const_iterator PVI;
  
   int npv = -1;
   int npv50ns = -1;
  
   for (PVI = PupInfo->begin(); PVI != PupInfo->end(); ++PVI) {
     int BX = PVI->getBunchCrossing();
  
     if (BX == 0) {
       npv = PVI->getPU_NumInteractions();
     }
  
     if (BX == 1) {
       npv50ns = PVI->getPU_NumInteractions();
     }
   }
  
   // Note: for the "uncorrelated" out-of-time pileup, reweighting is only done on the 50ns
   // "late" bunch (BX=+1), since that is basically the only one that matters in terms of
   // energy deposition.
  
   if (npv < 0) {
     std::cerr << " no in-time beam crossing found\n! ";
     std::cerr << " Returning event weight=0\n! ";
     return 0.;
   }
   if (npv50ns < 0) {
     std::cerr << " no out-of-time beam crossing found\n! ";
     std::cerr << " Returning event weight=0\n! ";
     return 0.;
   }
  
   int bin = weights_->GetXaxis()->FindBin(npv);
  
   double inTimeWeight = weights_->GetBinContent(bin);
  
   double TotalWeight = 1.0;
  
   TotalWeight = inTimeWeight;
  
   return TotalWeight;
 }

References newFWLiteAna::bin, L1TStage2uGTEmulatorClient_cff::BX, beam_dqm_sourceclient-live_cfg::cerr, MillePedeFileConverter_cfg::e, OldLumiSection_, pileupSumInfoTag_, and weights_.

Member Data Documentation

◆ Data_distr_

std::shared_ptr<TH1> edm::LumiReWeighting::Data_distr_

protected

Definition at line 140 of file LumiReWeighting.h.

Referenced by LumiReWeighting().

◆ dataFile_

std::shared_ptr<TFile> edm::LumiReWeighting::dataFile_

protected

Definition at line 134 of file LumiReWeighting.h.

Referenced by LumiReWeighting().

◆ dataFileName_

std::string edm::LumiReWeighting::dataFileName_

protected

Definition at line 129 of file LumiReWeighting.h.

Referenced by LumiReWeighting().

◆ DataHistName_

std::string edm::LumiReWeighting::DataHistName_

protected

Definition at line 131 of file LumiReWeighting.h.

Referenced by LumiReWeighting().

◆ FirstWarning_

bool edm::LumiReWeighting::FirstWarning_

protected

Definition at line 143 of file LumiReWeighting.h.

Referenced by LumiReWeighting(), and weight().

◆ generatedFile_

std::shared_ptr<TFile> edm::LumiReWeighting::generatedFile_

protected

Definition at line 133 of file LumiReWeighting.h.

Referenced by LumiReWeighting().

◆ generatedFileName_

std::string edm::LumiReWeighting::generatedFileName_

protected

Definition at line 128 of file LumiReWeighting.h.

Referenced by LumiReWeighting().

◆ GenHistName_

std::string edm::LumiReWeighting::GenHistName_

protected

Definition at line 130 of file LumiReWeighting.h.

Referenced by LumiReWeighting().

◆ MC_distr_

std::shared_ptr<TH1> edm::LumiReWeighting::MC_distr_

protected

Definition at line 139 of file LumiReWeighting.h.

Referenced by LumiReWeighting().

◆ OldLumiSection_

int edm::LumiReWeighting::OldLumiSection_

protected

Definition at line 142 of file LumiReWeighting.h.

Referenced by LumiReWeighting(), and weightOOT().

◆ pileupSumInfoTag_

edm::InputTag edm::LumiReWeighting::pileupSumInfoTag_

protected

Definition at line 132 of file LumiReWeighting.h.

Referenced by weight(), and weightOOT().

◆ weights_

std::shared_ptr<TH1> edm::LumiReWeighting::weights_

protected

Definition at line 135 of file LumiReWeighting.h.

Referenced by LumiReWeighting(), weight(), and weightOOT().

Public Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ LumiReWeighting() [1/3]

◆ LumiReWeighting() [2/3]

◆ LumiReWeighting() [3/3]

Member Function Documentation

◆ weight() [1/3]

◆ weight() [2/3]

◆ weight() [3/3]

◆ weightOOT()

Member Data Documentation

◆ Data_distr_

◆ dataFile_

◆ dataFileName_

◆ DataHistName_

◆ FirstWarning_

◆ generatedFile_

◆ generatedFileName_

◆ GenHistName_

◆ MC_distr_

◆ OldLumiSection_

◆ pileupSumInfoTag_

◆ weights_