zPdfUnc Class Reference

Inheritance diagram for zPdfUnc:

Public Member Functions
	zPdfUnc (const edm::ParameterSet &pset)
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 &event, const edm::EventSetup &setup) override
void	endJob () override

Private Attributes
double	accEtaMax_
double	accEtaMin_
double	accMassMax_
double	accMassMin_
double	accMassMinDenominator_
double	accPtMin_
double	angMax_
std::string	filename_
edm::EDGetTokenT< reco::GenParticleCollection >	genToken_
TH1F *	h_etaZMC_
TH1F *	h_mZMC_
TH1F *	h_nZ_
TH1F *	h_phiZMC_
TH1F *	h_ptZMC_
TH1F *	h_rapidityZMC_
TH1F *	h_thetaZMC_
TH1F *	h_weight_histo
TH1F *	hardeta
TH1F *	hardpt
bool	isMCatNLO_
double	massMax_
double	nAcc_
double	nAccReW_
unsigned int	nbinsAng_
unsigned int	nbinsMass_
unsigned int	nbinsPt_
double	nBothMuHasZHasGrandMa_
unsigned int	pdfmember_
edm::EDGetTokenT< std::vector< double > >	pdfweightsToken_
double	ptMax_
TH1F *	softeta
TH1F *	softpt

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
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 7 of file zPdfUnc.cc.

Constructor & Destructor Documentation

zPdfUnc::zPdfUnc

(

const edm::ParameterSet &

pset

)

Definition at line 49 of file zPdfUnc.cc.

References angMax_, gather_cfg::cout, h_etaZMC_, h_mZMC_, h_nZ_, h_phiZMC_, h_ptZMC_, h_rapidityZMC_, h_thetaZMC_, h_weight_histo, hardeta, hardpt, TFileDirectory::make(), massMax_, TFileService::mkdir(), nAcc_, nAccReW_, nbinsAng_, nbinsMass_, nbinsPt_, nBothMuHasZHasGrandMa_, ptMax_, softeta, and softpt.

    : genToken_(consumes<GenParticleCollection>(pset.getParameter<InputTag>("genParticles"))),
      pdfweightsToken_(consumes<std::vector<double> >(pset.getParameter<InputTag>("pdfweights"))),
      pdfmember_(pset.getUntrackedParameter<unsigned int>("pdfmember")),
      nbinsMass_(pset.getUntrackedParameter<unsigned int>("nbinsMass")),
      nbinsPt_(pset.getUntrackedParameter<unsigned int>("nbinsPt")),
      nbinsAng_(pset.getUntrackedParameter<unsigned int>("nbinsAng")),
      massMax_(pset.getUntrackedParameter<double>("massMax")),
      ptMax_(pset.getUntrackedParameter<double>("ptMax")),
      angMax_(pset.getUntrackedParameter<double>("angMax")),
      accPtMin_(pset.getUntrackedParameter<double>("accPtMin")),
      accMassMin_(pset.getUntrackedParameter<double>("accMassMin")),
      accMassMax_(pset.getUntrackedParameter<double>("accMassMax")),
      accEtaMin_(pset.getUntrackedParameter<double>("accEtaMin")),
      accEtaMax_(pset.getUntrackedParameter<double>("accEtaMax")),
      accMassMinDenominator_(pset.getUntrackedParameter<double>("accMassMinDenominator")),
      isMCatNLO_(pset.getUntrackedParameter<bool>("isMCatNLO")),
      filename_(pset.getUntrackedParameter<std::string>("outfilename")) {
  cout << ">>> Z Histogrammer constructor" << endl;
  Service<TFileService> fs;

  TFileDirectory ZMCHisto = fs->mkdir("ZMCHisto");
  h_nZ_ = ZMCHisto.make<TH1F>("ZNumber", "number of Z particles", 11, -0.5, 10.5);
  h_weight_histo = ZMCHisto.make<TH1F>("weight_histo", "weight_histo", 20, -10, 10);

  h_mZMC_ = ZMCHisto.make<TH1F>("ZMCMass", "Z MC mass (GeV/c^{2})", nbinsMass_, 0, massMax_);
  h_ptZMC_ = ZMCHisto.make<TH1F>("ZMCPt", "Z MC p_{t} (GeV/c)", nbinsPt_, 0, ptMax_);
  hardpt = ZMCHisto.make<TH1F>("hardpt", "hard muon p_{t} (GeV/c)", nbinsPt_, 0, ptMax_);
  softpt = ZMCHisto.make<TH1F>("softpt", "soft muon p_{t} (GeV/c)", nbinsPt_, 0, ptMax_);

  h_phiZMC_ = ZMCHisto.make<TH1F>("ZMCPhi", "Z MC #phi", nbinsAng_, -angMax_, angMax_);
  h_thetaZMC_ = ZMCHisto.make<TH1F>("ZMCTheta", "Z MC #theta", nbinsAng_, 0, angMax_);
  h_etaZMC_ = ZMCHisto.make<TH1F>("ZMCEta", "Z MC #eta", nbinsAng_, -angMax_, angMax_);
  h_rapidityZMC_ = ZMCHisto.make<TH1F>("ZMCRapidity", "Z MC y", nbinsAng_, -angMax_, angMax_);

  hardeta = ZMCHisto.make<TH1F>("hard muon eta", "hard muon #eta", nbinsAng_, -angMax_, angMax_);
  softeta = ZMCHisto.make<TH1F>("soft muon eta", "soft muon #eta", nbinsAng_, -angMax_, angMax_);
  nAcc_ = 0.;
  nAccReW_ = 0;
  nBothMuHasZHasGrandMa_ = 0;
}

Member Function Documentation

void zPdfUnc::analyze ( const edm::Event &  event, const edm::EventSetup &  setup  )

overrideprivate

Definition at line 91 of file zPdfUnc.cc.

References funct::abs(), accEtaMax_, accEtaMin_, accMassMax_, accMassMin_, accMassMinDenominator_, accPtMin_, gather_cfg::cout, HLT_2018_cff::eta1, HLT_2018_cff::eta2, relval_steps::gen(), genToken_, h_etaZMC_, h_mZMC_, h_phiZMC_, h_ptZMC_, h_rapidityZMC_, h_thetaZMC_, h_weight_histo, hardeta, hardpt, mps_fire::i, reco::CompositeRefCandidateT< D >::mother(), reco::Candidate::mother(), PDWG_BPHSkim_cff::muons, nAcc_, nAccReW_, nBothMuHasZHasGrandMa_, reco::CompositeRefCandidateT< D >::numberOfMothers(), reco::Candidate::numberOfMothers(), p4, pdfmember_, pdfweightsToken_, reco::Candidate::pdgId(), reco::LeafCandidate::pdgId(), HLT_2018_cff::pt1, HLT_2018_cff::pt2, softeta, softpt, mathSSE::sqrt(), reco::LeafCandidate::status(), mps_merge::weight, and HGCalRecHit_cfi::weights.

                                                                       {
  cout << ">>> Z Histogrammer analyze" << endl;

  Handle<GenParticleCollection> gen;
  Handle<vector<double> > pdfweights;

  event.getByToken(genToken_, gen);
  event.getByToken(pdfweightsToken_, pdfweights);

  // get weight and fill it to histogram
  std::vector<double> weights = (*pdfweights);

  unsigned int nmembers = weights.size();

  double weight = weights[pdfmember_];

  h_weight_histo->Fill(weight);

  cout << "found nmember: " << nmembers << endl;
  cout << "taken the member n " << pdfmember_ << " weight= " << weight << endl;

  std::vector<GenParticle> muons;

  for (unsigned int i = 0; i < gen->size(); ++i) {
    const GenParticle& muMC = (*gen)[i];
    // filling only muons coming form Z
    if (abs(muMC.pdgId()) == 13 && muMC.status() == 1 && muMC.numberOfMothers() > 0) {
      if (muMC.mother()->numberOfMothers() > 0) {
        cout << "I'm getting a muon \n"
             << "with "
             << "muMC.numberOfMothers()  " << muMC.numberOfMothers() << "\n the first mother has pdgId "
             << muMC.mother()->pdgId() << "with "
             << "muMC.mother()->numberOfMothers()  " << muMC.mother()->numberOfMothers()
             << "\n the first grandma has pdgId " << muMC.mother()->mother()->pdgId() << endl;
        if (muMC.mother()->mother()->pdgId() == 23)
          muons.push_back(muMC);
      }
    }
  }

  cout << "finally I selected " << muons.size() << " muons" << endl;

  // if there are at least two muons,
  // calculate invarant mass of first two and fill it into histogram

  double inv_mass = 0.0;
  double Zpt_ = 0.0;
  double Zeta_ = 0.0;
  double Ztheta_ = 0.0;
  double Zphi_ = 0.0;
  double Zrapidity_ = 0.0;

  if (muons.size() > 1) {
    if (muons[0].mother()->mother()->pdgId() == 23 && muons[1].mother()->mother()->pdgId() == 23)
      nBothMuHasZHasGrandMa_++;
    math::XYZTLorentzVector tot_momentum(muons[0].p4());
    math::XYZTLorentzVector mom2(muons[1].p4());
    tot_momentum += mom2;
    inv_mass = sqrt(tot_momentum.mass2());

    Zpt_ = tot_momentum.pt();
    Zeta_ = tot_momentum.eta();
    Ztheta_ = tot_momentum.theta();
    Zphi_ = tot_momentum.phi();
    Zrapidity_ = tot_momentum.Rapidity();

    // IMPORTANT: use the weight of the event ...

    double weight_sign = weight;
    //double weight_sign = 1.    ;
    // fill the denominator numbers only if mass>accMassMinDenominator and go on in that case
    if (inv_mass > accMassMinDenominator_) {
      h_mZMC_->Fill(inv_mass, weight_sign);
      h_ptZMC_->Fill(Zpt_, weight_sign);
      h_etaZMC_->Fill(Zeta_, weight_sign);
      h_thetaZMC_->Fill(Ztheta_, weight_sign);
      h_phiZMC_->Fill(Zphi_, weight_sign);
      h_rapidityZMC_->Fill(Zrapidity_, weight_sign);

      double pt1 = muons[0].pt();
      double pt2 = muons[1].pt();
      double eta1 = muons[0].eta();
      double eta2 = muons[1].eta();
      if (pt1 > pt2) {
        hardpt->Fill(pt1, weight_sign);
        softpt->Fill(pt2, weight_sign);
        hardeta->Fill(eta1, weight_sign);
        softeta->Fill(eta2, weight_sign);
      } else {
        hardpt->Fill(pt2, weight_sign);
        softpt->Fill(pt1, weight_sign);
        hardeta->Fill(eta2, weight_sign);
        softeta->Fill(eta1, weight_sign);
      }

      //evaluating the geometric acceptance
      if (pt1 >= accPtMin_ && pt2 >= accPtMin_ && fabs(eta1) >= accEtaMin_ && fabs(eta2) >= accEtaMin_ &&
          fabs(eta1) <= accEtaMax_ && fabs(eta2) <= accEtaMax_ && inv_mass >= accMassMin_ && inv_mass <= accMassMax_) {
        nAcc_++;
        nAccReW_ += weight;
      }
    }
  }
}

void zPdfUnc::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 196 of file zPdfUnc.cc.

References gather_cfg::cout, DEFINE_FWK_MODULE, runTheMatrix::err, filename_, h_mZMC_, nAcc_, nAccReW_, nbinsMass_, nBothMuHasZHasGrandMa_, and mathSSE::sqrt().

                     {
  cout << " number of events accepted :" << nAcc_ << endl;
  cout << " number of events accepted reweigthed :" << nAccReW_ << endl;
  double nev = h_mZMC_->GetEntries();
  double nev_weigthed = h_mZMC_->Integral();
  cout << " number of total events :" << nev << endl;
  cout << " number of total weighted events :" << nev_weigthed << endl;
  cout << " number of cases in which BothMuHasZHasGrandMa :" << nBothMuHasZHasGrandMa_ << endl;
  double eff = (double)nAcc_ / (double)h_mZMC_->GetEntries();
  double eff_rew = (double)nAccReW_ / (double)h_mZMC_->Integral(0, nbinsMass_ + 1);
  double err = sqrt(eff * (1. - eff) / (double)h_mZMC_->GetEntries());
  double err_rew = sqrt(eff_rew * (1. - eff_rew) / (double)h_mZMC_->Integral(0, nbinsMass_ + 1));
  cout << " geometric acceptance: " << eff << "+/-" << err << endl;
  cout << " geometric acceptance reweighted: " << eff_rew << "+/-" << err_rew << endl;

  ofstream myfile;
  myfile.open(filename_.c_str(), std::ios::app);
  myfile << eff << " " << eff_rew << " " << nev << " " << nev_weigthed << endl;
  myfile.close();
}

Member Data Documentation

double zPdfUnc::accEtaMax_

private

Definition at line 18 of file zPdfUnc.cc.

Referenced by analyze().

double zPdfUnc::accEtaMin_

private

Definition at line 18 of file zPdfUnc.cc.

Referenced by analyze().

double zPdfUnc::accMassMax_

private

Definition at line 18 of file zPdfUnc.cc.

Referenced by analyze().

double zPdfUnc::accMassMin_

private

Definition at line 18 of file zPdfUnc.cc.

Referenced by analyze().

double zPdfUnc::accMassMinDenominator_

private

Definition at line 18 of file zPdfUnc.cc.

Referenced by analyze().

double zPdfUnc::accPtMin_

private

Definition at line 18 of file zPdfUnc.cc.

Referenced by analyze().

double zPdfUnc::angMax_

private

Definition at line 17 of file zPdfUnc.cc.

Referenced by zPdfUnc().

std::string zPdfUnc::filename_

private

Definition at line 25 of file zPdfUnc.cc.

Referenced by endJob().

edm::EDGetTokenT<reco::GenParticleCollection> zPdfUnc::genToken_

private

Definition at line 14 of file zPdfUnc.cc.

Referenced by analyze().

TH1F * zPdfUnc::h_etaZMC_

private

Definition at line 20 of file zPdfUnc.cc.

Referenced by analyze(), and zPdfUnc().

TH1F* zPdfUnc::h_mZMC_

private

Definition at line 20 of file zPdfUnc.cc.

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

TH1F* zPdfUnc::h_nZ_

private

Definition at line 19 of file zPdfUnc.cc.

Referenced by zPdfUnc().

TH1F * zPdfUnc::h_phiZMC_

private

Definition at line 20 of file zPdfUnc.cc.

Referenced by analyze(), and zPdfUnc().

TH1F * zPdfUnc::h_ptZMC_

private

Definition at line 20 of file zPdfUnc.cc.

Referenced by analyze(), and zPdfUnc().

TH1F * zPdfUnc::h_rapidityZMC_

private

Definition at line 20 of file zPdfUnc.cc.

Referenced by analyze(), and zPdfUnc().

TH1F * zPdfUnc::h_thetaZMC_

private

Definition at line 20 of file zPdfUnc.cc.

Referenced by analyze(), and zPdfUnc().

TH1F* zPdfUnc::h_weight_histo

private

Definition at line 22 of file zPdfUnc.cc.

Referenced by analyze(), and zPdfUnc().

TH1F * zPdfUnc::hardeta

private

Definition at line 21 of file zPdfUnc.cc.

Referenced by analyze(), and zPdfUnc().

TH1F* zPdfUnc::hardpt

private

Definition at line 21 of file zPdfUnc.cc.

Referenced by analyze(), and zPdfUnc().

bool zPdfUnc::isMCatNLO_

private

Definition at line 23 of file zPdfUnc.cc.

double zPdfUnc::massMax_

private

Definition at line 17 of file zPdfUnc.cc.

Referenced by zPdfUnc().

double zPdfUnc::nAcc_

private

Definition at line 24 of file zPdfUnc.cc.

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

double zPdfUnc::nAccReW_

private

Definition at line 24 of file zPdfUnc.cc.

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

unsigned int zPdfUnc::nbinsAng_

private

Definition at line 16 of file zPdfUnc.cc.

Referenced by zPdfUnc().

unsigned int zPdfUnc::nbinsMass_

private

Definition at line 16 of file zPdfUnc.cc.

Referenced by endJob(), and zPdfUnc().

unsigned int zPdfUnc::nbinsPt_

private

Definition at line 16 of file zPdfUnc.cc.

Referenced by zPdfUnc().

double zPdfUnc::nBothMuHasZHasGrandMa_

private

Definition at line 24 of file zPdfUnc.cc.

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

unsigned int zPdfUnc::pdfmember_

private

Definition at line 16 of file zPdfUnc.cc.

Referenced by analyze().

edm::EDGetTokenT<std::vector<double> > zPdfUnc::pdfweightsToken_

private

Definition at line 15 of file zPdfUnc.cc.

Referenced by analyze().

double zPdfUnc::ptMax_

private

Definition at line 17 of file zPdfUnc.cc.

Referenced by zPdfUnc().

TH1F * zPdfUnc::softeta

private

Definition at line 21 of file zPdfUnc.cc.

Referenced by analyze(), and zPdfUnc().

TH1F * zPdfUnc::softpt

private

Definition at line 21 of file zPdfUnc.cc.

Referenced by analyze(), and zPdfUnc().