MCAcceptanceAnalyzer Class Reference

Inheritance diagram for MCAcceptanceAnalyzer:

Public Member Functions
	MCAcceptanceAnalyzer (const ParameterSet &cfg)
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 Event &, const EventSetup &) override
void	endJob () override

Private Attributes
long	nZToMuMu_
long	nZToMuMuMC_
long	nZToMuMuMCDen_
long	nZToMuMuMCMatched_
ZSelector	select_
ZSelector	select_OnlyMassCut_
long	selZToMuMu_
long	selZToMuMuMC_
long	selZToMuMuMCMatched_
InputTag	zToMuMu_
InputTag	zToMuMuMatched_
EDGetTokenT< GenParticleMatch >	zToMuMuMatchedToken_
InputTag	zToMuMuMC_
EDGetTokenT< CandidateView >	zToMuMuMCToken_
EDGetTokenT< CandidateView >	zToMuMuToken_

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 79 of file MCAcceptanceAnalyzer.cc.

Constructor & Destructor Documentation

MCAcceptanceAnalyzer::MCAcceptanceAnalyzer

(

const ParameterSet &

cfg

)

Definition at line 94 of file MCAcceptanceAnalyzer.cc.

    : zToMuMu_(cfg.getParameter<InputTag>("zToMuMu")),
      zToMuMuMC_(cfg.getParameter<InputTag>("zToMuMuMC")),
      zToMuMuMatched_(cfg.getParameter<InputTag>("zToMuMuMatched")),
      zToMuMuToken_(consumes<CandidateView>(zToMuMu_)),
      zToMuMuMCToken_(consumes<CandidateView>(zToMuMuMC_)),
      zToMuMuMatchedToken_(consumes<GenParticleMatch>(zToMuMuMatched_)),
      nZToMuMu_(0),
      selZToMuMu_(0),
      nZToMuMuMC_(0),
      selZToMuMuMC_(0),
      nZToMuMuMCMatched_(0),
      selZToMuMuMCMatched_(0),
      nZToMuMuMCDen_(0),
      select_(cfg.getParameter<double>("ptMin"),
              cfg.getParameter<double>("etaDau0Min"),
              cfg.getParameter<double>("etaDau0Max"),
              cfg.getParameter<double>("etaDau1Min"),
              cfg.getParameter<double>("etaDau1Max"),
              cfg.getParameter<double>("massMin"),
              cfg.getParameter<double>("massMax"),
              cfg.getParameter<double>("massMinZMC"),
              cfg.getParameter<double>("massMaxZMC")),
      select_OnlyMassCut_(-1,
                          -9999,
                          9999,
                          -9999,
                          9999,
                          0,
                          0,
                          cfg.getParameter<double>("massMinZMC"),
                          cfg.getParameter<double>("massMaxZMC")) {}

Member Function Documentation

void MCAcceptanceAnalyzer::analyze ( const Event &  evt, const EventSetup &    )

overrideprivate

Definition at line 127 of file MCAcceptanceAnalyzer.cc.

References gather_cfg::cout, HLTMuonOfflineAnalyzer_cfi::d0, reco::Candidate::daughter(), reco::Candidate::eta(), edm::Event::getByToken(), mps_fire::i, edm::Ref< C, T, F >::isNonnull(), visualization-live-secondInstance_cfg::m, EgHLTOffHistBins_cfi::mass, reco::Candidate::mass(), mcMuDaughter(), reco::Candidate::numberOfDaughters(), nZToMuMu_, nZToMuMuMC_, nZToMuMuMCDen_, nZToMuMuMCMatched_, reco::Candidate::p4(), reco::Candidate::pt(), edm::View< T >::refAt(), select_, select_OnlyMassCut_, selZToMuMu_, selZToMuMuMC_, selZToMuMuMCMatched_, edm::Association< C >::size(), edm::View< T >::size(), z, zToMuMuMatchedToken_, zToMuMuMCToken_, and zToMuMuToken_.

                                                                      {
  Handle<CandidateView> zToMuMu;
  evt.getByToken(zToMuMuToken_, zToMuMu);
  Handle<CandidateView> zToMuMuMC;
  evt.getByToken(zToMuMuMCToken_, zToMuMuMC);
  Handle<GenParticleMatch> zToMuMuMatched;
  evt.getByToken(zToMuMuMatchedToken_, zToMuMuMatched);
  //  long nZToMuMu = zToMuMu->size();
  long nZToMuMuMC = zToMuMuMC->size();
  long nZToMuMuMatched = zToMuMuMatched->size();

  // cout << ">>> " << zToMuMu_ << " has " << nZToMuMu << " entries" << endl;
  //cout << ">>> " << zToMuMuMC_ << " has " << nZToMuMuMC << " entries" << endl;
  //cout << ">>> " << zToMuMuMatched_ << " has " << nZToMuMuMatched << " entries" << endl;

  nZToMuMuMC_ += nZToMuMuMC;
  for (long i = 0; i < nZToMuMuMC; ++i) {
    const Candidate& z = (*zToMuMuMC)[i];
    if (select_(z) == 4)
      ++selZToMuMuMC_;
    if (select_OnlyMassCut_(z) == 4)
      ++nZToMuMuMCDen_;
  }

  for (long i = 0; i < nZToMuMuMatched; ++i) {
    const Candidate& z = (*zToMuMu)[i];
    CandidateBaseRef zRef = zToMuMu->refAt(i);
    GenParticleRef mcRef = (*zToMuMuMatched)[zRef];

    if (mcRef.isNonnull()) {  // z candidate matched to Z MC
      ++nZToMuMu_;
      ++nZToMuMuMCMatched_;

      int selectZ = select_(z);
      if (selectZ == 4)
        ++selZToMuMu_;

      int selectMC = select_(*mcRef);

      if (selectMC == 4)
        ++selZToMuMuMCMatched_;

      if (selectZ != selectMC) {
        cout << ">>> select reco: " << selectZ << ", select mc: " << selectMC << endl;
        if ((selectZ * selectMC) == 0)
          break;
        if (z.numberOfDaughters() > 1) {
          const Candidate *d0 = z.daughter(0), *d1 = z.daughter(1);
          if (mcRef->numberOfDaughters() > 1) {
            const Candidate *mcd0 = mcMuDaughter(mcRef->daughter(0)), *mcd1 = mcMuDaughter(mcRef->daughter(1));
            double m = z.mass(), mcm = (mcd0->p4() + mcd1->p4()).mass();
            cout << ">>> reco pt1, eta1: " << d0->pt() << ", " << d0->eta() << ", 2: " << d1->pt() << ", " << d1->eta()
                 << ", mass = " << m << endl;
            cout << ">>> mc   pt1, eta1: " << mcd0->pt() << ", " << mcd0->eta() << ", 2: " << mcd1->pt() << ", "
                 << mcd1->eta() << ", mass = " << mcm << endl;
          }
        }
      }
      // to avoid double counting
      if ((selectZ == 3) && (selectMC == 3))
        break;
    }
  }
}

void MCAcceptanceAnalyzer::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 192 of file MCAcceptanceAnalyzer.cc.

References gather_cfg::cout, DEFINE_FWK_MODULE, nZToMuMu_, nZToMuMuMC_, nZToMuMuMCDen_, nZToMuMuMCMatched_, selZToMuMu_, selZToMuMuMC_, selZToMuMuMCMatched_, mathSSE::sqrt(), zToMuMu_, and zToMuMuMC_.

                                  {
  double effZToMuMu = double(selZToMuMu_) / double(nZToMuMu_);
  double errZToMuMu = sqrt(effZToMuMu * (1. - effZToMuMu) / nZToMuMu_);
  double effZToMuMuMC = double(selZToMuMuMC_) / double(nZToMuMuMC_);
  double errZToMuMuMC = sqrt(effZToMuMuMC * (1. - effZToMuMuMC) / nZToMuMuMC_);
  double effZToMuMuMCDen = double(selZToMuMuMC_) / double(nZToMuMuMCDen_);
  double errZToMuMuMCDen = sqrt(effZToMuMuMCDen * (1. - effZToMuMuMCDen) / nZToMuMuMCDen_);
  double effZToMuMuMCMatched = double(selZToMuMuMCMatched_) / double(nZToMuMuMCMatched_);
  double errZToMuMuMCMatched = sqrt(effZToMuMuMCMatched * (1. - effZToMuMuMCMatched) / nZToMuMuMCMatched_);
  cout << ">>> " << zToMuMu_ << ": " << selZToMuMu_ << "/" << nZToMuMu_ << " = " << effZToMuMu << " +/- " << errZToMuMu
       << endl;
  cout << ">>> " << zToMuMuMC_ << " - matched: " << selZToMuMuMCMatched_ << "/" << nZToMuMuMCMatched_ << " = "
       << effZToMuMuMCMatched << " +/- " << errZToMuMuMCMatched << endl;
  cout << " if the two numbers above are the same we can neglete resolution effect and quote the acceptance as the "
          "number below.... "
       << endl;
  cout << "********* acceptance m>sampleMCMassCut (usually 20 or 40) ********  " << endl;
  cout << ">>> " << zToMuMuMC_ << ": " << selZToMuMuMC_ << "/" << nZToMuMuMC_ << " = " << effZToMuMuMC << " +/- "
       << errZToMuMuMC << endl;
  cout << "********* acceptance in the given mass range ********  " << endl;
  cout << ">>> " << zToMuMuMC_ << ": " << selZToMuMuMC_ << "/" << nZToMuMuMCDen_ << " = " << effZToMuMuMCDen << " +/- "
       << errZToMuMuMCDen << endl;
}

Member Data Documentation

long MCAcceptanceAnalyzer::nZToMuMu_

private

Definition at line 90 of file MCAcceptanceAnalyzer.cc.

Referenced by analyze(), and endJob().

long MCAcceptanceAnalyzer::nZToMuMuMC_

private

Definition at line 90 of file MCAcceptanceAnalyzer.cc.

Referenced by analyze(), and endJob().

long MCAcceptanceAnalyzer::nZToMuMuMCDen_

private

Definition at line 90 of file MCAcceptanceAnalyzer.cc.

Referenced by analyze(), and endJob().

long MCAcceptanceAnalyzer::nZToMuMuMCMatched_

private

Definition at line 90 of file MCAcceptanceAnalyzer.cc.

Referenced by analyze(), and endJob().

ZSelector MCAcceptanceAnalyzer::select_

private

Definition at line 91 of file MCAcceptanceAnalyzer.cc.

Referenced by analyze().

ZSelector MCAcceptanceAnalyzer::select_OnlyMassCut_

private

Definition at line 91 of file MCAcceptanceAnalyzer.cc.

Referenced by analyze().

long MCAcceptanceAnalyzer::selZToMuMu_

private

Definition at line 90 of file MCAcceptanceAnalyzer.cc.

Referenced by analyze(), and endJob().

long MCAcceptanceAnalyzer::selZToMuMuMC_

private

Definition at line 90 of file MCAcceptanceAnalyzer.cc.

Referenced by analyze(), and endJob().

long MCAcceptanceAnalyzer::selZToMuMuMCMatched_

private

Definition at line 90 of file MCAcceptanceAnalyzer.cc.

Referenced by analyze(), and endJob().

InputTag MCAcceptanceAnalyzer::zToMuMu_

private

Definition at line 86 of file MCAcceptanceAnalyzer.cc.

Referenced by endJob().

InputTag MCAcceptanceAnalyzer::zToMuMuMatched_

private

Definition at line 86 of file MCAcceptanceAnalyzer.cc.

EDGetTokenT<GenParticleMatch> MCAcceptanceAnalyzer::zToMuMuMatchedToken_

private

Definition at line 89 of file MCAcceptanceAnalyzer.cc.

Referenced by analyze().

InputTag MCAcceptanceAnalyzer::zToMuMuMC_

private

Definition at line 86 of file MCAcceptanceAnalyzer.cc.

Referenced by endJob().

EDGetTokenT<CandidateView> MCAcceptanceAnalyzer::zToMuMuMCToken_

private

Definition at line 88 of file MCAcceptanceAnalyzer.cc.

Referenced by analyze().

EDGetTokenT<CandidateView> MCAcceptanceAnalyzer::zToMuMuToken_

private

Definition at line 87 of file MCAcceptanceAnalyzer.cc.

Referenced by analyze().