Inheritance diagram for JetCorrectorDemo:

Public Types
typedef reco::Particle::LorentzVector	LorentzVector

Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Public Member Functions
	JetCorrectorDemo (const edm::ParameterSet &)

	~JetCorrectorDemo () override

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	beginJob () override

void	endJob () override

Private Attributes
edm::Service< TFileService >	fs

bool	mDebug

double	mEtaMax

double	mEtaMin

TH2F *	mJECvsEta

TH2F *	mJECvsPt

edm::EDGetTokenT< reco::JetCorrector >	mJetCorrector

int	mNGraphPoints

int	mNHistoPoints

std::string	mPayloadName

double	mPtMax

double	mPtMin

TRandom *	mRandom

TGraph *	mUncCorPt [100]

std::string	mUncertaintyFile

std::string	mUncertaintyTag

TGraph *	mUncEta [100]

bool	mUseCondDB

std::vector< double >	mVEta

TGraphErrors *	mVGraphCorPt [100]

TGraphErrors *	mVGraphEta [100]

TGraphErrors *	mVGraphPt [100]

std::vector< double >	mVPt

double	veta [100][1000]

double	vex_eta [100][1000]

double	vex_pt [100][1000]

double	vjec_eta [100][1000]

double	vjec_pt [100][1000]

double	vjecUnc_eta [100][1000]

double	vjecUnc_pt [100][1000]

double	vpt [100][1000]

double	vptcor [100][1000]

double	vUnc_eta [100][1000]

double	vUnc_pt [100][1000]

Additional Inherited Members
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 28 of file JetCorrectorDemo.cc.

Member Typedef Documentation

typedef reco::Particle::LorentzVector JetCorrectorDemo::LorentzVector

Definition at line 32 of file JetCorrectorDemo.cc.

Constructor & Destructor Documentation

JetCorrectorDemo::JetCorrectorDemo ( const edm::ParameterSet & iConfig )

explicit

Definition at line 58 of file JetCorrectorDemo.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), mDebug, mEtaMax, mEtaMin, mJetCorrector, mNGraphPoints, mNHistoPoints, mPayloadName, mPtMax, mPtMin, mUncertaintyFile, mUncertaintyTag, mUseCondDB, mVEta, mVPt, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                  {
   mJetCorrector = consumes<reco::JetCorrector>(iConfig.getParameter<edm::InputTag>("JetCorrector"));
   mPayloadName = iConfig.getParameter<std::string>("PayloadName");
   mUncertaintyTag = iConfig.getParameter<std::string>("UncertaintyTag");
   mUncertaintyFile = iConfig.getParameter<std::string>("UncertaintyFile");
   mNHistoPoints = iConfig.getParameter<int>("NHistoPoints");
   mNGraphPoints = iConfig.getParameter<int>("NGraphPoints");
   mEtaMin = iConfig.getParameter<double>("EtaMin");
   mEtaMax = iConfig.getParameter<double>("EtaMax");
   mPtMin = iConfig.getParameter<double>("PtMin");
   mPtMax = iConfig.getParameter<double>("PtMax");
   mVEta = iConfig.getParameter<std::vector<double> >("VEta");
   mVPt = iConfig.getParameter<std::vector<double> >("VPt");
   mDebug = iConfig.getUntrackedParameter<bool>("Debug", false);
   mUseCondDB = iConfig.getUntrackedParameter<bool>("UseCondDB", false);
 }

JetCorrectorDemo::~JetCorrectorDemo ( )

override

Definition at line 75 of file JetCorrectorDemo.cc.

75 {}

Member Function Documentation

void JetCorrectorDemo::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 77 of file JetCorrectorDemo.cc.

References reco::JetCorrector::correction(), mitigatedMETSequence_cff::corrector, gather_cfg::cout, HLT_2018_cff::dEta, MillePedeFileConverter_cfg::e, PVValHelper::eta, edm::FileInPath::fullPath(), edm::EventSetup::get(), edm::Event::getByToken(), JetCorrectionUncertainty::getUncertainty(), mps_fire::i, LEDCalibrationChannels::ieta, jetMETDQMOfflineSource_cff::jec, mDebug, mEtaMax, mEtaMin, mJECvsEta, mJECvsPt, mJetCorrector, mNGraphPoints, mNHistoPoints, mPayloadName, mPtMax, mPtMin, mRandom, mUncertaintyFile, mUncertaintyTag, mUseCondDB, mVEta, mVPt, funct::pow(), jets_cff::rawPt, JetCorrectionUncertainty::setJetEta(), JetCorrectionUncertainty::setJetPt(), createJobs::tmp, veta, vex_eta, vex_pt, vjec_eta, vjec_pt, vjecUnc_eta, vjecUnc_pt, vpt, vptcor, vUnc_eta, and vUnc_pt.

                                                                                   {
   edm::Handle<reco::JetCorrector> corrector;
   iEvent.getByToken(mJetCorrector, corrector);
   JetCorrectionUncertainty *jecUnc(nullptr);
   if (!mUncertaintyTag.empty()) {
     if (mUseCondDB) {
       edm::ESHandle<JetCorrectorParametersCollection> JetCorParColl;
       iSetup.get<JetCorrectionsRecord>().get(mPayloadName, JetCorParColl);
       JetCorrectorParameters const &JetCorPar = (*JetCorParColl)[mUncertaintyTag];
       jecUnc = new JetCorrectionUncertainty(JetCorPar);
       std::cout << "Configured Uncertainty from CondDB" << std::endl;
     } else {
       edm::FileInPath fip("CondFormats/JetMETObjects/data/" + mUncertaintyFile + ".txt");
       jecUnc = new JetCorrectionUncertainty(fip.fullPath());
     }
   }
   double jec, rawPt, corPt, eta, unc;
   TLorentzVector P4;
   double dEta = (mEtaMax - mEtaMin) / mNGraphPoints;
   for (int i = 0; i < mNHistoPoints; i++) {
     rawPt = mRandom->Uniform(mPtMin, mPtMax);
     eta = mRandom->Uniform(mEtaMin, mEtaMax);
     P4.SetPtEtaPhiE(rawPt, eta, 0, 0);
     LorentzVector rawP4(P4.Px(), P4.Py(), P4.Pz(), P4.E());
     jec = corrector->correction(rawP4);
     mJECvsEta->Fill(eta, jec);
     mJECvsPt->Fill(rawPt, jec);
   }
   //--------- Pt Graphs ------------------
   for (unsigned ieta = 0; ieta < mVEta.size(); ieta++) {
     double rPt = pow((3500. / TMath::CosH(mVEta[ieta])) / mPtMin, 1. / mNGraphPoints);
     for (int i = 0; i < mNGraphPoints; i++) {
       rawPt = mPtMin * pow(rPt, i);
       eta = mVEta[ieta];
       vpt[ieta][i] = rawPt;
       P4.SetPtEtaPhiE(rawPt, eta, 0, 0);
       LorentzVector rawP4(P4.Px(), P4.Py(), P4.Pz(), P4.E());
       jec = corrector->correction(rawP4);  // the jec is a function of the raw pt
       vjec_eta[ieta][i] = jec;
       vptcor[ieta][i] = rawPt * jec;
       unc = 0.0;
       if (!mUncertaintyTag.empty()) {
         jecUnc->setJetEta(eta);
         jecUnc->setJetPt(rawPt * jec);  // the uncertainty is a function of the corrected pt
         unc = jecUnc->getUncertainty(true);
       }
       vjecUnc_eta[ieta][i] = unc * jec;
       vUnc_eta[ieta][i] = unc;
       vex_eta[ieta][i] = 0.0;
       if (mDebug)
         std::cout << rawPt << " " << eta << " " << jec << " " << rawPt * jec << " " << unc << std::endl;
     }
   }
   //--------- Eta Graphs -------------------
   for (unsigned ipt = 0; ipt < mVPt.size(); ipt++) {
     for (int i = 0; i < mNGraphPoints; i++) {
       eta = mEtaMin + i * dEta;
       corPt = mVPt[ipt];
       veta[ipt][i] = eta;
       //---------- find the raw pt -----------
       double e = 1.0;
       int nLoop(0);
       rawPt = corPt;
       while (e > 0.0001 && nLoop < 10) {
         P4.SetPtEtaPhiE(rawPt, eta, 0, 0);
         LorentzVector rawP4(P4.Px(), P4.Py(), P4.Pz(), P4.E());
         jec = corrector->correction(rawP4);
         double tmp = rawPt * jec;
         e = fabs(tmp - corPt) / corPt;
         if (jec > 0)
           rawPt = corPt / jec;
         nLoop++;
       }
       //--------- calculate the jec for the rawPt --------
       P4.SetPtEtaPhiE(rawPt, eta, 0, 0);
       LorentzVector rawP4(P4.Px(), P4.Py(), P4.Pz(), P4.E());
       jec = corrector->correction(rawP4);  // the jec is a function of the raw pt
       vjec_pt[ipt][i] = jec;
       unc = 0.0;
       if (!mUncertaintyTag.empty()) {
         jecUnc->setJetEta(eta);
         jecUnc->setJetPt(corPt);  // the uncertainty is a function of the corrected pt
         unc = jecUnc->getUncertainty(true);
       }
       vjecUnc_pt[ipt][i] = unc * jec;
       vUnc_pt[ipt][i] = unc;
       vex_pt[ipt][i] = 0.0;
       if (mDebug)
         std::cout << rawPt << " " << eta << " " << jec << " " << rawPt * jec << " " << unc << std::endl;
     }
   }
 }

void JetCorrectorDemo::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 170 of file JetCorrectorDemo.cc.

References Exception, fs, TFileService::make(), mEtaMax, mEtaMin, mJECvsEta, mJECvsPt, mNGraphPoints, mPtMax, mPtMin, mRandom, mVEta, and mVPt.

                                 {
   if (mNGraphPoints > 1000)
     throw cms::Exception("JetCorrectorDemo", "Too many graph points !!! Maximum is 1000 !!!");
   if (mVEta.size() > 100)
     throw cms::Exception("JetCorrectorDemo", "Too many eta values !!! Maximum is 100 !!!");
   if (mVPt.size() > 100)
     throw cms::Exception("JetCorrectorDemo", "Too many pt values !!! Maximum is 100 !!!");
   mJECvsEta = fs->make<TH2F>("JECvsEta", "JECvsEta", 200, mEtaMin, mEtaMax, 100, 0, 5);
   mJECvsPt = fs->make<TH2F>("JECvsPt", "JECvsPt", 200, mPtMin, mPtMax, 100, 0, 5);
   mRandom = new TRandom();
   mRandom->SetSeed(0);
 }

void JetCorrectorDemo::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 183 of file JetCorrectorDemo.cc.

References DEFINE_FWK_MODULE, fs, LEDCalibrationChannels::ieta, TFileService::make(), mNGraphPoints, mUncCorPt, mUncEta, mVEta, mVGraphCorPt, mVGraphEta, mVGraphPt, mVPt, Skims_PA_cff::name, veta, vex_eta, vex_pt, vjec_eta, vjec_pt, vjecUnc_eta, vjecUnc_pt, vpt, vptcor, vUnc_eta, and vUnc_pt.

                               {
   char name[1000];
   for (unsigned ipt = 0; ipt < mVPt.size(); ipt++) {
     mVGraphEta[ipt] = fs->make<TGraphErrors>(mNGraphPoints, veta[ipt], vjec_pt[ipt], vex_pt[ipt], vjecUnc_pt[ipt]);
     sprintf(name, "JEC_vs_Eta_CorPt%1.1f", mVPt[ipt]);
     mVGraphEta[ipt]->SetName(name);
     mUncEta[ipt] = fs->make<TGraph>(mNGraphPoints, veta[ipt], vUnc_pt[ipt]);
     sprintf(name, "UNC_vs_Eta_CorPt%1.1f", mVPt[ipt]);
     mUncEta[ipt]->SetName(name);
   }
   for (unsigned ieta = 0; ieta < mVEta.size(); ieta++) {
     mVGraphPt[ieta] =
         fs->make<TGraphErrors>(mNGraphPoints, vpt[ieta], vjec_eta[ieta], vex_eta[ieta], vjecUnc_eta[ieta]);
     sprintf(name, "JEC_vs_RawPt_eta%1.1f", mVEta[ieta]);
     mVGraphPt[ieta]->SetName(name);
     mVGraphCorPt[ieta] =
         fs->make<TGraphErrors>(mNGraphPoints, vptcor[ieta], vjec_eta[ieta], vex_eta[ieta], vjecUnc_eta[ieta]);
     sprintf(name, "JEC_vs_CorPt_eta%1.1f", mVEta[ieta]);
     mVGraphCorPt[ieta]->SetName(name);
     mUncCorPt[ieta] = fs->make<TGraph>(mNGraphPoints, vptcor[ieta], vUnc_eta[ieta]);
     sprintf(name, "UNC_vs_CorPt_eta%1.1f", mVEta[ieta]);
     mUncCorPt[ieta]->SetName(name);
   }
 }

Member Data Documentation

edm::Service<TFileService> JetCorrectorDemo::fs

private

Definition at line 48 of file JetCorrectorDemo.cc.

Referenced by beginJob(), and endJob().

bool JetCorrectorDemo::mDebug

private

Definition at line 41 of file JetCorrectorDemo.cc.

Referenced by analyze(), and JetCorrectorDemo().

double JetCorrectorDemo::mEtaMax

private

Definition at line 43 of file JetCorrectorDemo.cc.

Referenced by analyze(), beginJob(), and JetCorrectorDemo().

double JetCorrectorDemo::mEtaMin

private

Definition at line 43 of file JetCorrectorDemo.cc.

Referenced by analyze(), beginJob(), and JetCorrectorDemo().

TH2F* JetCorrectorDemo::mJECvsEta

private

Definition at line 49 of file JetCorrectorDemo.cc.

Referenced by analyze(), and beginJob().

TH2F * JetCorrectorDemo::mJECvsPt

private

Definition at line 49 of file JetCorrectorDemo.cc.

Referenced by analyze(), and beginJob().

edm::EDGetTokenT<reco::JetCorrector> JetCorrectorDemo::mJetCorrector

private

Definition at line 39 of file JetCorrectorDemo.cc.

Referenced by analyze(), and JetCorrectorDemo().

int JetCorrectorDemo::mNGraphPoints

private

Definition at line 42 of file JetCorrectorDemo.cc.

Referenced by analyze(), beginJob(), endJob(), and JetCorrectorDemo().

int JetCorrectorDemo::mNHistoPoints

private

Definition at line 42 of file JetCorrectorDemo.cc.

Referenced by analyze(), and JetCorrectorDemo().

std::string JetCorrectorDemo::mPayloadName

private

Definition at line 40 of file JetCorrectorDemo.cc.

Referenced by analyze(), and JetCorrectorDemo().

double JetCorrectorDemo::mPtMax

private

Definition at line 43 of file JetCorrectorDemo.cc.

Referenced by analyze(), beginJob(), and JetCorrectorDemo().

double JetCorrectorDemo::mPtMin

private

Definition at line 43 of file JetCorrectorDemo.cc.

Referenced by analyze(), beginJob(), and JetCorrectorDemo().

TRandom* JetCorrectorDemo::mRandom

private

Definition at line 52 of file JetCorrectorDemo.cc.

Referenced by analyze(), and beginJob().

TGraph * JetCorrectorDemo::mUncCorPt[100]

private

Definition at line 51 of file JetCorrectorDemo.cc.

Referenced by endJob().

std::string JetCorrectorDemo::mUncertaintyFile

private

Definition at line 40 of file JetCorrectorDemo.cc.

Referenced by analyze(), and JetCorrectorDemo().

std::string JetCorrectorDemo::mUncertaintyTag

private

Definition at line 40 of file JetCorrectorDemo.cc.

Referenced by analyze(), and JetCorrectorDemo().

TGraph* JetCorrectorDemo::mUncEta[100]

private

Definition at line 51 of file JetCorrectorDemo.cc.

Referenced by endJob().

bool JetCorrectorDemo::mUseCondDB

private

Definition at line 41 of file JetCorrectorDemo.cc.

Referenced by analyze(), and JetCorrectorDemo().

std::vector<double> JetCorrectorDemo::mVEta

private

Definition at line 44 of file JetCorrectorDemo.cc.

Referenced by analyze(), beginJob(), endJob(), and JetCorrectorDemo().

TGraphErrors * JetCorrectorDemo::mVGraphCorPt[100]

private

Definition at line 50 of file JetCorrectorDemo.cc.

Referenced by endJob().

TGraphErrors* JetCorrectorDemo::mVGraphEta[100]

private

Definition at line 50 of file JetCorrectorDemo.cc.

Referenced by endJob().

TGraphErrors * JetCorrectorDemo::mVGraphPt[100]

private

Definition at line 50 of file JetCorrectorDemo.cc.

Referenced by endJob().

std::vector<double> JetCorrectorDemo::mVPt

private

Definition at line 44 of file JetCorrectorDemo.cc.

Referenced by analyze(), beginJob(), endJob(), and JetCorrectorDemo().

double JetCorrectorDemo::veta[100][1000]

private

Definition at line 45 of file JetCorrectorDemo.cc.

Referenced by analyze(), and endJob().

double JetCorrectorDemo::vex_eta[100][1000]

private

Definition at line 46 of file JetCorrectorDemo.cc.

Referenced by analyze(), and endJob().

double JetCorrectorDemo::vex_pt[100][1000]

private

Definition at line 46 of file JetCorrectorDemo.cc.

Referenced by analyze(), and endJob().

double JetCorrectorDemo::vjec_eta[100][1000]

private

Definition at line 45 of file JetCorrectorDemo.cc.

Referenced by analyze(), and endJob().

double JetCorrectorDemo::vjec_pt[100][1000]

private

Definition at line 45 of file JetCorrectorDemo.cc.

Referenced by analyze(), and endJob().

double JetCorrectorDemo::vjecUnc_eta[100][1000]

private

Definition at line 46 of file JetCorrectorDemo.cc.

Referenced by analyze(), and endJob().

double JetCorrectorDemo::vjecUnc_pt[100][1000]

private

Definition at line 46 of file JetCorrectorDemo.cc.

Referenced by analyze(), and endJob().

double JetCorrectorDemo::vpt[100][1000]

private

Definition at line 45 of file JetCorrectorDemo.cc.

Referenced by analyze(), and endJob().

double JetCorrectorDemo::vptcor[100][1000]

private

Definition at line 45 of file JetCorrectorDemo.cc.

Referenced by analyze(), and endJob().

double JetCorrectorDemo::vUnc_eta[100][1000]

private

Definition at line 46 of file JetCorrectorDemo.cc.

Referenced by analyze(), and endJob().

double JetCorrectorDemo::vUnc_pt[100][1000]

private

Definition at line 46 of file JetCorrectorDemo.cc.

Referenced by analyze(), and endJob().

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation