Inheritance diagram for ResolutionCreator:

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

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::SharedResources >
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const noexcept

bool	wantsStreamRuns () const noexcept

	~EDAnalyzerBase () 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

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

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices 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::EDGetTokenT< std::vector< pat::Electron > >	electronsToken_

std::vector< double >	etabinVals_

int	etanrbins

TF1 *	fResEtaBin [10][20]

TF1 *	fResPtEtaBin [10][20][20]

edm::EDGetTokenT< TtGenEvent >	genEvtToken_

TH1F *	hResEtaBin [10][20]

TH1F *	hResPtEtaBin [10][20][20]

edm::EDGetTokenT< std::vector< pat::Jet > >	jetsToken_

std::string	labelName_

int	matchingAlgo_

double	maxDist_

edm::EDGetTokenT< std::vector< pat::MET > >	metsToken_

double	minDR_

edm::EDGetTokenT< std::vector< pat::Muon > >	muonsToken_

int	nrFilled

std::string	objectType_

std::vector< double >	pTbinVals_

int	ptnrbins

edm::EDGetTokenT< std::vector< pat::Tau > >	tausToken_

TTree *	tResVar

bool	useDeltaR_

bool	useMaxDist_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

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

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 38 of file ResolutionCreator.cc.

Constructor & Destructor Documentation

◆ ResolutionCreator()

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

explicit

Definition at line 75 of file ResolutionCreator.cc.

References electronsToken_, etabinVals_, genEvtToken_, edm::ParameterSet::getParameter(), ProducerED_cfi::InputTag, jetsToken_, labelName_, metsToken_, minDR_, muonsToken_, nrFilled, objectType_, pTbinVals_, AlCaHLTBitMon_QueryRunRegistry::string, and tausToken_.

                                                                    {
   usesResource("TFileService");
 
   // input parameters
   genEvtToken_ = consumes<TtGenEvent>(edm::InputTag("genEvt"));
   objectType_ = iConfig.getParameter<std::string>("object");
   labelName_ = iConfig.getParameter<std::string>("label");
   if (objectType_ == "electron")
     electronsToken_ = consumes<std::vector<pat::Electron> >(edm::InputTag(labelName_));
   else if (objectType_ == "muon")
     muonsToken_ = consumes<std::vector<pat::Muon> >(edm::InputTag(labelName_));
   else if (objectType_ == "lJets" || objectType_ == "bJets")
     jetsToken_ = consumes<std::vector<pat::Jet> >(edm::InputTag(labelName_));
   else if (objectType_ == "met")
     metsToken_ = consumes<std::vector<pat::MET> >(edm::InputTag(labelName_));
   else if (objectType_ == "tau")
     tausToken_ = consumes<std::vector<pat::Tau> >(edm::InputTag(labelName_));
   if (objectType_ != "met") {
     etabinVals_ = iConfig.getParameter<std::vector<double> >("etabinValues");
   }
   pTbinVals_ = iConfig.getParameter<std::vector<double> >("pTbinValues");
   minDR_ = iConfig.getParameter<double>("minMatchingDR");
 
   nrFilled = 0;
 }

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 106 of file ResolutionCreator.cc.

References funct::abs(), b, bgen, electronAnalyzer_cfi::DeltaR, MillePedeFileConverter_cfg::e, pwdgSkimBPark_cfi::electrons, electronsToken_, reco::LeafCandidate::energy(), etabinVals_, etanrbins, TtGenEvtProducer_cfi::genEvt, genEvtToken_, hResPtEtaBin, iEvent, dqmiolumiharvest::j, PDWG_EXODelayedJetMET_cff::jets, jetsToken_, visualization-live-secondInstance_cfg::m, singleTopDQM_cfi::mets, metsToken_, minDR_, reco::CompositeRefCandidateT< D >::mother(), reco::Candidate::mother(), DiMuonV_cfg::muons, muonsToken_, nrFilled, reco::CompositeRefCandidateT< D >::numberOfMothers(), reco::Candidate::numberOfMothers(), objectType_, AlCaHLTBitMon_ParallelJobs::p, reco::LeafCandidate::p4(), reco::Candidate::pdgId(), reco::LeafCandidate::pdgId(), hitfit::phidiff(), DiDispStaMuonMonitor_cfi::pt, pTbinVals_, multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, reco::LeafCandidate::status(), metsig::tau, Tau3MuMonitor_cff::taus, and tausToken_.

                                                                                    {
   // Get the gen and cal object fourvector
   std::vector<reco::Particle *> p4gen, p4rec;
 
   edm::Handle<TtGenEvent> genEvt;
   iEvent.getByToken(genEvtToken_, genEvt);
 
   if (genEvt->particles().size() < 10)
     return;
 
   if (objectType_ == "electron") {
     edm::Handle<std::vector<pat::Electron> >
         electrons;  //to calculate the ResolutionCreator for the electrons, i used the TopElectron instead of the AOD information
     iEvent.getByToken(electronsToken_, electrons);
     for (size_t e = 0; e < electrons->size(); e++) {
       for (size_t p = 0; p < genEvt->particles().size(); p++) {
         if ((std::abs(genEvt->particles()[p].pdgId()) == 11) &&
             (ROOT::Math::VectorUtil::DeltaR(genEvt->particles()[p].p4(), (*electrons)[e].p4()) < minDR_)) {
           //p4gen.push_back(new reco::Particle(genEvt->particles()[p]));
           //p4rec.push_back(new reco::Particle((pat::Electron)((*electrons)[e])));
         }
       }
     }
   } else if (objectType_ == "muon") {
     edm::Handle<std::vector<pat::Muon> > muons;
     iEvent.getByToken(muonsToken_, muons);
     for (size_t m = 0; m < muons->size(); m++) {
       for (size_t p = 0; p < genEvt->particles().size(); p++) {
         if ((std::abs(genEvt->particles()[p].pdgId()) == 13) &&
             (ROOT::Math::VectorUtil::DeltaR(genEvt->particles()[p].p4(), (*muons)[m].p4()) < minDR_)) {
           //p4gen.push_back(new reco::Particle(genEvt->particles()[p]));
           //p4rec.push_back(new reco::Particle((pat::Muon)((*muons)[m])));
         }
       }
     }
   } else if (objectType_ == "lJets") {
     edm::Handle<std::vector<pat::Jet> > jets;
     iEvent.getByToken(jetsToken_, jets);
     if (jets->size() >= 4) {
       for (unsigned int j = 0; j < 4; j++) {
         for (size_t p = 0; p < genEvt->particles().size(); p++) {
           if ((std::abs(genEvt->particles()[p].pdgId()) < 5) &&
               (ROOT::Math::VectorUtil::DeltaR(genEvt->particles()[p].p4(), (*jets)[j].p4()) < minDR_)) {
             //p4gen.push_back(new reco::Particle(genEvt->particles()[p]));
             //p4rec.push_back(new reco::Particle((pat::Jet)(*jets)[j]));
           }
         }
       }
     }
   } else if (objectType_ == "bJets") {
     edm::Handle<std::vector<pat::Jet> > jets;
     iEvent.getByToken(jetsToken_, jets);
     if (jets->size() >= 4) {
       for (unsigned int j = 0; j < 4; j++) {
         for (size_t p = 0; p < genEvt->particles().size(); p++) {
           if ((std::abs(genEvt->particles()[p].pdgId()) == 5) &&
               (ROOT::Math::VectorUtil::DeltaR(genEvt->particles()[p].p4(), (*jets)[j].p4()) < minDR_)) {
             //p4gen.push_back(new reco::Particle(genEvt->particles()[p]));
             //p4rec.push_back(new reco::Particle((pat::Jet)(*jets)[j]));
           }
         }
       }
     }
   } else if (objectType_ == "met") {
     edm::Handle<std::vector<pat::MET> > mets;
     iEvent.getByToken(metsToken_, mets);
     if (!mets->empty()) {
       if (genEvt->isSemiLeptonic() && genEvt->singleNeutrino() != nullptr &&
           ROOT::Math::VectorUtil::DeltaR(genEvt->singleNeutrino()->p4(), (*mets)[0].p4()) < minDR_) {
         //p4gen.push_back(new reco::Particle(0,genEvt->singleNeutrino()->p4(),math::XYZPoint()));
         //p4rec.push_back(new reco::Particle((pat::MET)((*mets)[0])));
       }
     }
   } else if (objectType_ == "tau") {
     edm::Handle<std::vector<pat::Tau> > taus;
     iEvent.getByToken(tausToken_, taus);
     for (std::vector<pat::Tau>::const_iterator tau = taus->begin(); tau != taus->end(); ++tau) {
       // find the tau (if any) that matches a MC tau from W
       reco::GenParticle genLepton = *(tau->genLepton());
       if (std::abs(genLepton.pdgId()) == 15 && genLepton.status() == 2 && genLepton.numberOfMothers() > 0 &&
           std::abs(genLepton.mother(0)->pdgId()) == 15 && genLepton.mother(0)->numberOfMothers() > 0 &&
           std::abs(genLepton.mother(0)->mother(0)->pdgId()) == 24 &&
           ROOT::Math::VectorUtil::DeltaR(genLepton.p4(), tau->p4()) < minDR_) {
       }
       //p4gen.push_back(new reco::Particle(genLepton));
       //p4rec.push_back(new reco::Particle(*tau));
     }
   }
   // Fill the object's value
   for (unsigned m = 0; m < p4gen.size(); m++) {
     double Egen = p4gen[m]->energy();
     double Thetagen = p4gen[m]->theta();
     double Phigen = p4gen[m]->phi();
     double Etgen = p4gen[m]->et();
     double Etagen = p4gen[m]->eta();
     double Ecal = p4rec[m]->energy();
     double Thetacal = p4rec[m]->theta();
     double Phical = p4rec[m]->phi();
     double Etcal = p4rec[m]->et();
     double Etacal = p4rec[m]->eta();
     double phidiff = Phical - Phigen;
     if (phidiff > 3.14159)
       phidiff = 2. * 3.14159 - phidiff;
     if (phidiff < -3.14159)
       phidiff = -phidiff - 2. * 3.14159;
 
     // find eta and et bin
     int etabin = 0;
     if (etanrbins > 1) {
       for (unsigned int b = 0; b < etabinVals_.size() - 1; b++) {
         if (fabs(Etacal) > etabinVals_[b])
           etabin = b;
       }
     }
 
     int ptbin = 0;
     for (unsigned int b = 0; b < pTbinVals_.size() - 1; b++) {
       if (p4rec[m]->pt() > pTbinVals_[b])
         ptbin = b;
     }
 
     // calculate the resolution on "a", "b", "c" & "d" according to the definition (CMS-NOTE-2006-023):
     // p = a*|p_meas|*u_1 + b*u_2 + c*u_3
     // E(fit) = E_meas * d
     //
     // with u_1 = p/|p_meas|
     //      u_3 = (u_z x u_1)/|u_z x u_1|
     //      u_2 = (u_1 x u_3)/|u_1 x u_3|
     //
     // The initial parameters values are chosen like (a, b, c, d) = (1., 0., 0., 1.)
 
     // 1/ calculate the unitary vectors of the basis u_1, u_2, u_3
     ROOT::Math::SVector<double, 3> pcalvec(p4rec[m]->px(), p4rec[m]->py(), p4rec[m]->pz());
     ROOT::Math::SVector<double, 3> pgenvec(p4gen[m]->px(), p4gen[m]->py(), p4gen[m]->pz());
 
     ROOT::Math::SVector<double, 3> u_z(0, 0, 1);
     ROOT::Math::SVector<double, 3> u_1 = ROOT::Math::Unit(pcalvec);
     ROOT::Math::SVector<double, 3> u_3 = ROOT::Math::Cross(u_z, u_1) / ROOT::Math::Mag(ROOT::Math::Cross(u_z, u_1));
     ROOT::Math::SVector<double, 3> u_2 = ROOT::Math::Cross(u_1, u_3) / ROOT::Math::Mag(ROOT::Math::Cross(u_1, u_3));
     double acal = 1.;
     double bcal = 0.;
     double ccal = 0.;
     double dcal = 1.;
     double agen = ROOT::Math::Dot(pgenvec, u_1) / ROOT::Math::Mag(pcalvec);
     double bgen = ROOT::Math::Dot(pgenvec, u_2);
     double cgen = ROOT::Math::Dot(pgenvec, u_3);
     double dgen = Egen / Ecal;
 
     //fill histograms
     ++nrFilled;
     hResPtEtaBin[0][etabin][ptbin]->Fill(acal - agen);
     hResPtEtaBin[1][etabin][ptbin]->Fill(bcal - bgen);
     hResPtEtaBin[2][etabin][ptbin]->Fill(ccal - cgen);
     hResPtEtaBin[3][etabin][ptbin]->Fill(dcal - dgen);
     hResPtEtaBin[4][etabin][ptbin]->Fill(Thetacal - Thetagen);
     hResPtEtaBin[5][etabin][ptbin]->Fill(phidiff);
     hResPtEtaBin[6][etabin][ptbin]->Fill(Etcal - Etgen);
     hResPtEtaBin[7][etabin][ptbin]->Fill(Etacal - Etagen);
 
     delete p4gen[m];
     delete p4rec[m];
   }
 }

◆ beginJob()

void ResolutionCreator::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 271 of file ResolutionCreator.cc.

References b, edm::errors::Configuration, defaults_cfi::etabins, etabinVals_, etanrbins, Exception, fResEtaBin, fResPtEtaBin, make_classfiles::fs, hResEtaBin, hResPtEtaBin, objectType_, LepHTMonitor_cff::ptbins, pTbinVals_, ptnrbins, and tResVar.

                                  {
   edm::Service<TFileService> fs;
   if (!fs)
     throw edm::Exception(edm::errors::Configuration, "TFileService missing from configuration!");
 
   // input constants
   TString resObsName[8] = {"_ares", "_bres", "_cres", "_dres", "_thres", "_phres", "_etres", "_etares"};
   int resObsNrBins = 120;
   if ((objectType_ == "muon") || (objectType_ == "electron"))
     resObsNrBins = 80;
   std::vector<double> resObsMin, resObsMax;
   if (objectType_ == "electron") {
     resObsMin.push_back(-0.15);
     resObsMin.push_back(-0.2);
     resObsMin.push_back(-0.1);
     resObsMin.push_back(-0.15);
     resObsMin.push_back(-0.0012);
     resObsMin.push_back(-0.009);
     resObsMin.push_back(-16);
     resObsMin.push_back(-0.0012);
     resObsMax.push_back(0.15);
     resObsMax.push_back(0.2);
     resObsMax.push_back(0.1);
     resObsMax.push_back(0.15);
     resObsMax.push_back(0.0012);
     resObsMax.push_back(0.009);
     resObsMax.push_back(16);
     resObsMax.push_back(0.0012);
   } else if (objectType_ == "muon") {
     resObsMin.push_back(-0.15);
     resObsMin.push_back(-0.1);
     resObsMin.push_back(-0.05);
     resObsMin.push_back(-0.15);
     resObsMin.push_back(-0.004);
     resObsMin.push_back(-0.003);
     resObsMin.push_back(-8);
     resObsMin.push_back(-0.004);
     resObsMax.push_back(0.15);
     resObsMax.push_back(0.1);
     resObsMax.push_back(0.05);
     resObsMax.push_back(0.15);
     resObsMax.push_back(0.004);
     resObsMax.push_back(0.003);
     resObsMax.push_back(8);
     resObsMax.push_back(0.004);
   } else if (objectType_ == "tau") {
     resObsMin.push_back(-1.);
     resObsMin.push_back(-10.);
     resObsMin.push_back(-10);
     resObsMin.push_back(-1.);
     resObsMin.push_back(-0.1);
     resObsMin.push_back(-0.1);
     resObsMin.push_back(-80);
     resObsMin.push_back(-0.1);
     resObsMax.push_back(1.);
     resObsMax.push_back(10.);
     resObsMax.push_back(10);
     resObsMax.push_back(1.);
     resObsMax.push_back(0.1);
     resObsMax.push_back(0.1);
     resObsMax.push_back(50);
     resObsMax.push_back(0.1);
   } else if (objectType_ == "lJets" || objectType_ == "bJets") {
     resObsMin.push_back(-1.);
     resObsMin.push_back(-10.);
     resObsMin.push_back(-10.);
     resObsMin.push_back(-1.);
     resObsMin.push_back(-0.4);
     resObsMin.push_back(-0.6);
     resObsMin.push_back(-80);
     resObsMin.push_back(-0.6);
     resObsMax.push_back(1.);
     resObsMax.push_back(10.);
     resObsMax.push_back(10.);
     resObsMax.push_back(1.);
     resObsMax.push_back(0.4);
     resObsMax.push_back(0.6);
     resObsMax.push_back(80);
     resObsMax.push_back(0.6);
   } else {
     resObsMin.push_back(-2.);
     resObsMin.push_back(-150.);
     resObsMin.push_back(-150.);
     resObsMin.push_back(-2.);
     resObsMin.push_back(-6);
     resObsMin.push_back(-6);
     resObsMin.push_back(-180);
     resObsMin.push_back(-6);
     resObsMax.push_back(3.);
     resObsMax.push_back(150.);
     resObsMax.push_back(150.);
     resObsMax.push_back(3.);
     resObsMax.push_back(6);
     resObsMax.push_back(6);
     resObsMax.push_back(180);
     resObsMax.push_back(6);
   }
 
   const char *resObsVsPtFit[8] = {"[0]+[1]*exp(-[2]*x)",
                                   "[0]+[1]*exp(-[2]*x)",
                                   "[0]+[1]*exp(-[2]*x)",
                                   "[0]+[1]*exp(-[2]*x)",
                                   "[0]+[1]*exp(-[2]*x)",
                                   "[0]+[1]*exp(-[2]*x)",
                                   "pol1",
                                   "[0]+[1]*exp(-[2]*x)"};
 
   ptnrbins = pTbinVals_.size() - 1;
   double *ptbins = new double[pTbinVals_.size()];
   for (unsigned int b = 0; b < pTbinVals_.size(); b++)
     ptbins[b] = pTbinVals_[b];
   double *etabins;
   if (objectType_ != "met") {
     etanrbins = etabinVals_.size() - 1;
     etabins = new double[etabinVals_.size()];
     for (unsigned int b = 0; b < etabinVals_.size(); b++)
       etabins[b] = etabinVals_[b];
   } else {
     etanrbins = 1;
     etabins = new double[2];
     etabins[0] = 0;
     etabins[1] = 5.;
   }
 
   //define the histograms booked
   for (Int_t ro = 0; ro < 8; ro++) {
     for (Int_t etab = 0; etab < etanrbins; etab++) {
       for (Int_t ptb = 0; ptb < ptnrbins; ptb++) {
         TString obsName = objectType_;
         obsName += resObsName[ro];
         obsName += "_etabin";
         obsName += etab;
         obsName += "_ptbin";
         obsName += ptb;
         hResPtEtaBin[ro][etab][ptb] = fs->make<TH1F>(obsName, obsName, resObsNrBins, resObsMin[ro], resObsMax[ro]);
         fResPtEtaBin[ro][etab][ptb] = fs->make<TF1>("F_" + obsName, "gaus");
       }
       TString obsName2 = objectType_;
       obsName2 += resObsName[ro];
       obsName2 += "_etabin";
       obsName2 += etab;
       hResEtaBin[ro][etab] = fs->make<TH1F>(obsName2, obsName2, ptnrbins, ptbins);
       fResEtaBin[ro][etab] =
           fs->make<TF1>("F_" + obsName2, resObsVsPtFit[ro], pTbinVals_[0], pTbinVals_[pTbinVals_.size() - 1]);
     }
   }
   tResVar = fs->make<TTree>("tResVar", "Resolution tree");
 
   delete[] etabins;
   delete[] ptbins;
 }

◆ endJob()

void ResolutionCreator::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 424 of file ResolutionCreator.cc.

References relativeConstraints::error, PVValHelper::eta, etabinVals_, etanrbins, fResEtaBin, fResPtEtaBin, hResEtaBin, hResPtEtaBin, createfilelist::int, nrFilled, objectType_, DiDispStaMuonMonitor_cfi::pt, pTbinVals_, ptnrbins, isotrackApplyRegressor::range, tResVar, and relativeConstraints::value.

                                {
   TString resObsName2[8] = {"a", "b", "c", "d", "theta", "phi", "et", "eta"};
   Int_t ro = 0;
   Double_t pt = 0.;
   Double_t eta = 0.;
   Double_t value, error;
 
   tResVar->Branch("Pt", &pt, "Pt/D");
   tResVar->Branch("Eta", &eta, "Eta/D");
   tResVar->Branch("ro", &ro, "ro/I");
   tResVar->Branch("value", &value, "value/D");
   tResVar->Branch("error", &error, "error/D");
 
   for (ro = 0; ro < 8; ro++) {
     for (int etab = 0; etab < etanrbins; etab++) {
       //CD set eta at the center of the bin
       eta = etanrbins > 1 ? (etabinVals_[etab] + etabinVals_[etab + 1]) / 2. : 2.5;
       for (int ptb = 0; ptb < ptnrbins; ptb++) {
         //CD set pt at the center of the bin
         pt = (pTbinVals_[ptb] + pTbinVals_[ptb + 1]) / 2.;
         double maxcontent = 0.;
         int maxbin = 0;
         for (int nb = 1; nb < hResPtEtaBin[ro][etab][ptb]->GetNbinsX(); nb++) {
           if (hResPtEtaBin[ro][etab][ptb]->GetBinContent(nb) > maxcontent) {
             maxcontent = hResPtEtaBin[ro][etab][ptb]->GetBinContent(nb);
             maxbin = nb;
           }
         }
         int range = (int)(hResPtEtaBin[ro][etab][ptb]->GetNbinsX() / 6);  //in order that ~1/3 of X-axis range is fitted
         fResPtEtaBin[ro][etab][ptb]->SetRange(hResPtEtaBin[ro][etab][ptb]->GetBinCenter(maxbin - range),
                                               hResPtEtaBin[ro][etab][ptb]->GetBinCenter(maxbin + range));
         fResPtEtaBin[ro][etab][ptb]->SetParameters(hResPtEtaBin[ro][etab][ptb]->GetMaximum(),
                                                    hResPtEtaBin[ro][etab][ptb]->GetMean(),
                                                    hResPtEtaBin[ro][etab][ptb]->GetRMS());
         hResPtEtaBin[ro][etab][ptb]->Fit(fResPtEtaBin[ro][etab][ptb]->GetName(), "RQ SERIAL");
         hResEtaBin[ro][etab]->SetBinContent(ptb + 1, fResPtEtaBin[ro][etab][ptb]->GetParameter(2));
         hResEtaBin[ro][etab]->SetBinError(ptb + 1, fResPtEtaBin[ro][etab][ptb]->GetParError(2));
         //CD: Fill the tree
         value = fResPtEtaBin[ro][etab][ptb]->GetParameter(2);  //parameter value
         error = fResPtEtaBin[ro][etab][ptb]->GetParError(2);   //parameter error
         tResVar->Fill();
       }
       //CD: add a fake entry in pt=0 for the NN training
       // for that, use a linear extrapolation.
       pt = 0.;
       value = ((pTbinVals_[0] + pTbinVals_[1]) / 2.) *
                   (fResPtEtaBin[ro][etab][0]->GetParameter(2) - fResPtEtaBin[ro][etab][1]->GetParameter(2)) /
                   ((pTbinVals_[2] - pTbinVals_[0]) / 2.) +
               fResPtEtaBin[ro][etab][0]->GetParameter(2);
       error = fResPtEtaBin[ro][etab][0]->GetParError(2) + fResPtEtaBin[ro][etab][1]->GetParError(2);
       tResVar->Fill();
       // standard fit
       hResEtaBin[ro][etab]->Fit(fResEtaBin[ro][etab]->GetName(), "RQ SERIAL");
     }
   }
   if (objectType_ == "lJets" && nrFilled == 0)
     edm::LogProblem("SummaryError") << "No plots filled for light jets \n";
   if (objectType_ == "bJets" && nrFilled == 0)
     edm::LogProblem("SummaryError") << "No plots filled for bjets \n";
   if (objectType_ == "muon" && nrFilled == 0)
     edm::LogProblem("SummaryError") << "No plots filled for muons \n";
   if (objectType_ == "electron" && nrFilled == 0)
     edm::LogProblem("SummaryError") << "No plots filled for electrons \n";
   if (objectType_ == "tau" && nrFilled == 0)
     edm::LogProblem("SummaryError") << "No plots filled for taus \n";
   if (objectType_ == "met" && nrFilled == 0)
     edm::LogProblem("SummaryError") << "No plots filled for met \n";
 
   edm::LogVerbatim("MainResults") << " \n\n";
   edm::LogVerbatim("MainResults") << " ----------------------------------------------";
   edm::LogVerbatim("MainResults") << " ----------------------------------------------";
   edm::LogVerbatim("MainResults") << " Resolutions on " << objectType_ << " with nrfilled: " << nrFilled;
   edm::LogVerbatim("MainResults") << " ----------------------------------------------";
   edm::LogVerbatim("MainResults") << " ----------------------------------------------";
   if (nrFilled != 0 && objectType_ != "met") {
     for (ro = 0; ro < 8; ro++) {
       edm::LogVerbatim("MainResults") << "-------------------- ";
       edm::LogVerbatim("MainResults") << "\n Resolutions on " << resObsName2[ro] << "\n";
       edm::LogVerbatim("MainResults") << "-------------------- ";
       for (int etab = 0; etab < etanrbins; etab++) {
         if (nrFilled != 0 && ro != 6) {
           if (etab == 0) {
             edm::LogVerbatim("MainResults")
                 << "if(fabs(eta)<" << etabinVals_[etab + 1] << ") res = " << fResEtaBin[ro][etab]->GetParameter(0)
                 << "+" << fResEtaBin[ro][etab]->GetParameter(1) << "*exp(-(" << fResEtaBin[ro][etab]->GetParameter(2)
                 << "*pt));";
           } else {
             edm::LogVerbatim("MainResults")
                 << "else if(fabs(eta)<" << etabinVals_[etab + 1] << ") res = " << fResEtaBin[ro][etab]->GetParameter(0)
                 << "+" << fResEtaBin[ro][etab]->GetParameter(1) << "*exp(-(" << fResEtaBin[ro][etab]->GetParameter(2)
                 << "*pt));";
           }
         } else if (nrFilled != 0 && ro == 6) {
           if (etab == 0) {
             edm::LogVerbatim("MainResults")
                 << "if(fabs(eta)<" << etabinVals_[etab + 1] << ") res = " << fResEtaBin[ro][etab]->GetParameter(0)
                 << "+" << fResEtaBin[ro][etab]->GetParameter(1) << "*pt;";
           } else {
             edm::LogVerbatim("MainResults")
                 << "else if(fabs(eta)<" << etabinVals_[etab + 1] << ") res = " << fResEtaBin[ro][etab]->GetParameter(0)
                 << "+" << fResEtaBin[ro][etab]->GetParameter(1) << "*pt;";
           }
         }
       }
     }
   } else if (nrFilled != 0 && objectType_ == "met") {
     for (ro = 0; ro < 8; ro++) {
       edm::LogVerbatim("MainResults") << "-------------------- ";
       edm::LogVerbatim("MainResults") << "\n Resolutions on " << resObsName2[ro] << "\n";
       edm::LogVerbatim("MainResults") << "-------------------- ";
       if (nrFilled != 0 && ro != 6) {
         edm::LogVerbatim("MainResults") << "res = " << fResEtaBin[ro][0]->GetParameter(0) << "+"
                                         << fResEtaBin[ro][0]->GetParameter(1) << "*exp(-("
                                         << fResEtaBin[ro][0]->GetParameter(2) << "*pt));";
       } else if (nrFilled != 0 && ro == 6) {
         edm::LogVerbatim("MainResults") << "res = " << fResEtaBin[ro][0]->GetParameter(0) << "+"
                                         << fResEtaBin[ro][0]->GetParameter(1) << "*pt;";
       }
     }
   }
 }

Member Data Documentation

◆ electronsToken_

edm::EDGetTokenT<std::vector<pat::Electron> > ResolutionCreator::electronsToken_

private

Definition at line 50 of file ResolutionCreator.cc.

Referenced by analyze(), and ResolutionCreator().

◆ etabinVals_

std::vector<double> ResolutionCreator::etabinVals_

private

Definition at line 55 of file ResolutionCreator.cc.

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

◆ etanrbins

int ResolutionCreator::etanrbins

private

Definition at line 61 of file ResolutionCreator.cc.

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

◆ fResEtaBin

TF1* ResolutionCreator::fResEtaBin[10][20]

private

Definition at line 66 of file ResolutionCreator.cc.

Referenced by beginJob(), and endJob().

◆ fResPtEtaBin

TF1* ResolutionCreator::fResPtEtaBin[10][20][20]

private

Definition at line 65 of file ResolutionCreator.cc.

Referenced by beginJob(), and endJob().

◆ genEvtToken_

edm::EDGetTokenT<TtGenEvent> ResolutionCreator::genEvtToken_

private

Definition at line 48 of file ResolutionCreator.cc.

Referenced by analyze(), and ResolutionCreator().

◆ hResEtaBin

TH1F* ResolutionCreator::hResEtaBin[10][20]

private

Definition at line 68 of file ResolutionCreator.cc.

Referenced by beginJob(), and endJob().

◆ hResPtEtaBin

TH1F* ResolutionCreator::hResPtEtaBin[10][20][20]

private

Definition at line 67 of file ResolutionCreator.cc.

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

◆ jetsToken_

edm::EDGetTokenT<std::vector<pat::Jet> > ResolutionCreator::jetsToken_

private

Definition at line 52 of file ResolutionCreator.cc.

Referenced by analyze(), and ResolutionCreator().

◆ labelName_

std::string ResolutionCreator::labelName_

private

Definition at line 49 of file ResolutionCreator.cc.

Referenced by ResolutionCreator().

◆ matchingAlgo_

int ResolutionCreator::matchingAlgo_

private

Definition at line 57 of file ResolutionCreator.cc.

◆ maxDist_

double ResolutionCreator::maxDist_

private

Definition at line 60 of file ResolutionCreator.cc.

◆ metsToken_

edm::EDGetTokenT<std::vector<pat::MET> > ResolutionCreator::metsToken_

private

Definition at line 53 of file ResolutionCreator.cc.

Referenced by analyze(), and ResolutionCreator().

◆ minDR_

double ResolutionCreator::minDR_

private

Definition at line 56 of file ResolutionCreator.cc.

Referenced by analyze(), and ResolutionCreator().

◆ muonsToken_

edm::EDGetTokenT<std::vector<pat::Muon> > ResolutionCreator::muonsToken_

private

Definition at line 51 of file ResolutionCreator.cc.

Referenced by analyze(), and ResolutionCreator().

◆ nrFilled

int ResolutionCreator::nrFilled

private

Definition at line 62 of file ResolutionCreator.cc.

Referenced by analyze(), endJob(), and ResolutionCreator().

◆ objectType_

std::string ResolutionCreator::objectType_

private

Definition at line 49 of file ResolutionCreator.cc.

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

◆ pTbinVals_

std::vector<double> ResolutionCreator::pTbinVals_

private

Definition at line 55 of file ResolutionCreator.cc.

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

◆ ptnrbins

int ResolutionCreator::ptnrbins

private

Definition at line 61 of file ResolutionCreator.cc.

Referenced by beginJob(), and endJob().

◆ tausToken_

edm::EDGetTokenT<std::vector<pat::Tau> > ResolutionCreator::tausToken_

private

Definition at line 54 of file ResolutionCreator.cc.

Referenced by analyze(), and ResolutionCreator().

◆ tResVar

TTree* ResolutionCreator::tResVar

private

Definition at line 69 of file ResolutionCreator.cc.

Referenced by beginJob(), and endJob().

◆ useDeltaR_

bool ResolutionCreator::useDeltaR_

private

Definition at line 59 of file ResolutionCreator.cc.

◆ useMaxDist_

bool ResolutionCreator::useMaxDist_

private

Definition at line 58 of file ResolutionCreator.cc.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ ResolutionCreator()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ endJob()

Member Data Documentation

◆ electronsToken_

◆ etabinVals_

◆ etanrbins

◆ fResEtaBin

◆ fResPtEtaBin

◆ genEvtToken_

◆ hResEtaBin

◆ hResPtEtaBin

◆ jetsToken_

◆ labelName_

◆ matchingAlgo_

◆ maxDist_

◆ metsToken_

◆ minDR_

◆ muonsToken_

◆ nrFilled

◆ objectType_

◆ pTbinVals_

◆ ptnrbins

◆ tausToken_

◆ tResVar

◆ useDeltaR_

◆ useMaxDist_