Public Member Functions
double	getP (const reco::Candidate::LorentzVector &genMuonPlusP4_beforeRad, const reco::Candidate::LorentzVector &genMuonPlusP4_afterRad, const reco::Candidate::LorentzVector &genMuonMinusP4_beforeRad, const reco::Candidate::LorentzVector &genMuonMinusP4_afterRad)

bool	isWithinBounds (double muPlusEn, double muMinusEn)

	lutEntryType (TFile &inputFile, const std::string &directory, double minMuPlusEn, double maxMuPlusEn, double minMuMinusEn, double maxMuMinusEn)

	~lutEntryType ()

Public Attributes
TH2 *	lut_

double	maxMuMinusEn_

double	maxMuPlusEn_

double	minMuMinusEn_

double	minMuPlusEn_

Detailed Description

Definition at line 47 of file MuonRadiationCorrWeightProducer.h.

Constructor & Destructor Documentation

MuonRadiationCorrWeightProducer::lutEntryType::lutEntryType	(	TFile &	inputFile,
		const std::string &	directory,
		double	minMuPlusEn,
		double	maxMuPlusEn,
		double	minMuMinusEn,
		double	maxMuMinusEn
	)

inline

Definition at line 49 of file MuonRadiationCorrWeightProducer.h.

References python.multivaluedict::append(), createBeamHaloJobs::directory, Exception, lumiPlot::lut, lut_, maxMuMinusEn_, maxMuPlusEn_, minMuMinusEn_, minMuPlusEn_, and AlCaHLTBitMon_QueryRunRegistry::string.

       : minMuPlusEn_(minMuPlusEn),
     maxMuPlusEn_(maxMuPlusEn),
     minMuMinusEn_(minMuMinusEn),
     maxMuMinusEn_(maxMuMinusEn),
     lut_(0)
     {
       TString lutName = directory;
       if ( lutName.Length() > 0 && !lutName.EndsWith("/") ) lutName.Append("/");
       lutName.Append(directory); // CV: repeat directory, as it is done by GenMuonRadCorrAnalyzer
       if ( lutName.Length() > 0 && !lutName.EndsWith("/") ) lutName.Append("/");
       lutName.Append("genMuonRadCorr");
       if      ( minMuPlusEn_  < 0. && maxMuPlusEn_  < 0. ) lutName.Append("");
       else if (                       maxMuPlusEn_  < 0. ) lutName.Append(Form("_muPlusEnGt%1.0f", minMuPlusEn_));
       else if ( minMuPlusEn_  < 0.                       ) lutName.Append(Form("_muPlusEnLt%1.0f", maxMuPlusEn_));
       else                                                 lutName.Append(Form("_muPlusEn%1.0fto%1.0f", minMuPlusEn_, maxMuPlusEn_));
       if      ( minMuMinusEn_ < 0. && maxMuMinusEn_ < 0. ) lutName.Append("");
       else if (                       maxMuMinusEn_ < 0. ) lutName.Append(Form("_muMinusEnGt%1.0f", minMuMinusEn_));
       else if ( minMuMinusEn_ < 0.                       ) lutName.Append(Form("_muMinusEnLt%1.0f", maxMuMinusEn_));
       else                                                 lutName.Append(Form("_muMinusEn%1.0fto%1.0f", minMuMinusEn_, maxMuMinusEn_));
       TH2* lut = dynamic_cast<TH2*>(inputFile.Get(lutName.Data()));
       if ( !lut ) 
     throw cms::Exception("MuonRadiationCorrWeightProducer") 
       << " Failed to load LUT = " << lutName << " from file = " << inputFile.GetName() << " !!\n";
       lut_ = (TH2*)lut->Clone(std::string(lut->GetName()).append("_cloned").data());
       if ( !lut_->GetSumw2N() ) lut_->Sumw2();
       if ( lut_->Integral() != 0. ) lut_->Scale(1./lut_->Integral());
     }

MuonRadiationCorrWeightProducer::lutEntryType::~lutEntryType ( )

inline

Definition at line 77 of file MuonRadiationCorrWeightProducer.h.

References lut_.

     {
       delete lut_;
     }

Member Function Documentation

double MuonRadiationCorrWeightProducer::lutEntryType::getP	(	const reco::Candidate::LorentzVector &	genMuonPlusP4_beforeRad,
		const reco::Candidate::LorentzVector &	genMuonPlusP4_afterRad,
		const reco::Candidate::LorentzVector &	genMuonMinusP4_beforeRad,
		const reco::Candidate::LorentzVector &	genMuonMinusP4_afterRad
	)

inline

Definition at line 90 of file MuonRadiationCorrWeightProducer.h.

References lut_, x, and y.

     {
       TAxis* xAxis = lut_->GetXaxis();
       double x = (genMuonPlusP4_beforeRad.E() - genMuonPlusP4_afterRad.E())/genMuonPlusP4_beforeRad.E();    
       int binX = xAxis->FindBin(x);
       if ( binX < 1 ) binX = 1;
       int numBinsX = xAxis->GetNbins();
       if ( binX > numBinsX ) binX = numBinsX;
       
       TAxis* yAxis = lut_->GetYaxis();
       double y = (genMuonMinusP4_beforeRad.E() - genMuonMinusP4_afterRad.E())/genMuonMinusP4_beforeRad.E();
       int binY = yAxis->FindBin(y);
       if ( binY < 1 ) binY = 1;
       int numBinsY = yAxis->GetNbins();
       if ( binY > numBinsY ) binY = numBinsY;
       
       return lut_->GetBinContent(binX, binY);
     }

bool MuonRadiationCorrWeightProducer::lutEntryType::isWithinBounds	(	double	muPlusEn,
		double	muMinusEn
	)

inline

Definition at line 82 of file MuonRadiationCorrWeightProducer.h.

References maxMuMinusEn_, maxMuPlusEn_, minMuMinusEn_, and minMuPlusEn_.

     {
       if ( (minMuPlusEn_  <= 0. || muPlusEn  > minMuPlusEn_  ) &&
        (maxMuPlusEn_  <= 0. || muPlusEn  < maxMuPlusEn_  ) &&
        (minMuMinusEn_ <= 0. || muMinusEn > minMuMinusEn_ ) &&
        (maxMuMinusEn_ <= 0. || muMinusEn < maxMuMinusEn_ ) ) return true;
       else return false;
     }