MCEfficiencyAnalyzer Class Reference

Inheritance diagram for MCEfficiencyAnalyzer:

Public Member Functions
virtual void	analyze (const edm::Event &event, const edm::EventSetup &setup) override
virtual void	endJob () override
	MCEfficiencyAnalyzer (const edm::ParameterSet &pset)
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Attributes
double	deltaRStacut_
double	etacut_
EDGetTokenT < GenParticleCollection >	genParticlesToken_
EDGetTokenT< CandMatchMap >	MuonsMapToken_
EDGetTokenT< CandidateCollection >	MuonsToken_
int	nMuMC
int	nMureco
int	nNotMuMatching
int	nSta
int	nTrk
double	ptcut_
EDGetTokenT< CandMatchMap >	StandAloneMapToken_
EDGetTokenT< CandidateCollection >	StandAloneToken_
EDGetTokenT< CandMatchMap >	TracksMapToken_
EDGetTokenT< CandidateCollection >	TracksToken_
vector< double >	v_
EDGetTokenT< CandidateCollection >	zMuMuToken_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
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 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 29 of file MC_Efficiency_Analyzer.cc.

Constructor & Destructor Documentation

MCEfficiencyAnalyzer::MCEfficiencyAnalyzer ( const edm::ParameterSet &  pset )

inline

Definition at line 32 of file MC_Efficiency_Analyzer.cc.

                                                    :
    zMuMuToken_( consumes<CandidateCollection>( pset.getParameter<InputTag>( "zMuMu" ) ) ),
    MuonsToken_( consumes<CandidateCollection>( pset.getParameter<InputTag>( "Muons" ) ) ),
    MuonsMapToken_( consumes<CandMatchMap>( pset.getParameter<InputTag>( "MuonsMap" ) ) ),
    TracksToken_( consumes<CandidateCollection>( pset.getParameter<InputTag>( "Tracks" ) ) ),
    TracksMapToken_( consumes<CandMatchMap>( pset.getParameter<InputTag>( "TracksMap" ) ) ),
    genParticlesToken_( consumes<GenParticleCollection>( pset.getParameter<InputTag>( "genParticles" ) ) ),
    StandAloneToken_( consumes<CandidateCollection>( pset.getParameter<InputTag>( "StandAlone" ) ) ),
    StandAloneMapToken_( consumes<CandMatchMap>( pset.getParameter<InputTag>( "StandAloneMap" ) ) ),
    etacut_( pset.getParameter<double>( "etacut" ) ),
    ptcut_( pset.getParameter<double>( "ptcut" ) ),
    deltaRStacut_( pset.getParameter<double>( "deltaRStacut" ) )

  {
    nMuMC =0; nMureco =0; nTrk=0; nSta=0; nNotMuMatching =0 ;
  }

Member Function Documentation

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

inlineoverridevirtual

Implements edm::EDAnalyzer.

Definition at line 49 of file MC_Efficiency_Analyzer.cc.

References funct::abs(), gather_cfg::cout, reco::Candidate::daughter(), deltaR(), alignCSCRings::e, reco::Candidate::eta(), genParticleCandidates2GenParticles_cfi::genParticles, i, gen::k, combine::key, ConfigFiles::l, n, reco::Candidate::numberOfDaughters(), reco::Candidate::pdgId(), reco::Candidate::pt(), reco::Candidate::status(), ntuplemaker::status, and ZMuMuAnalysisNtupler_cff::zMuMu.

  {

    Handle<CandidateCollection> zMuMu;
    event.getByToken(zMuMuToken_, zMuMu);


    Handle<CandidateCollection> Muons;
    event.getByToken(MuonsToken_, Muons);

    Handle<CandidateCollection> Tracks;
    event.getByToken(TracksToken_, Tracks);

    Handle<CandidateCollection> StandAlone;
    event.getByToken(StandAloneToken_, StandAlone);

    Handle<CandMatchMap> MuonsMap;
    event.getByToken(MuonsMapToken_, MuonsMap);

    Handle<CandMatchMap> TracksMap;
    event.getByToken(TracksMapToken_, TracksMap);

    Handle<CandMatchMap> StandAloneMap;
    event.getByToken(StandAloneMapToken_, StandAloneMap);

    Handle<GenParticleCollection> genParticles;
    event.getByToken(genParticlesToken_, genParticles);

    //Getting muons from Z MC
    for( unsigned int k = 0; k < genParticles->size(); k++ )
      {
    const Candidate & ZCand = (*genParticles)[ k ];
    int status = ZCand.status();

    if (ZCand.pdgId()==23&& status==3 )
      {
        // positive muons
        const Candidate  *  muCand1 =   ZCand.daughter(0);
        if (muCand1->status()==3)
          {
        for (unsigned int d =0 ; d< muCand1->numberOfDaughters(); d++)
          {
            const Candidate * muCandidate = muCand1->daughter(d);
            if (muCandidate->pdgId() == muCand1->pdgId() )
              {
            muCand1 = muCand1->daughter(d);
              }
          }
          }
        // negative muons
        const Candidate  *  muCand2 =   ZCand.daughter(1);
        if (muCand2->status()==3)
          {
        for (unsigned int e =0 ; e< muCand2->numberOfDaughters(); e++)
          {
            const Candidate * muCandidate = muCand2->daughter(e);
            if (muCandidate->pdgId() == muCand2->pdgId() )
              {
            muCand2 = muCand2->daughter(e);
              }
          }
          }

        double deltaR_Mu_Sta =0;
        int nMurecoTemp = nMureco;
        // getting mu matched
        CandMatchMap::const_iterator i;
        for(i = MuonsMap->begin(); i != MuonsMap->end(); i++ )
          {
        const Candidate/* & reco = * i -> key,*/ &  mc =  * i -> val;
        if ((muCand1 == &mc) && (mc.pt()>ptcut_) && (abs(mc.eta()<etacut_)))
          {
            nMuMC++;
            nMureco++;
            break;
          }
          }
        if (nMureco == nMurecoTemp ) // I.E. MU RECO NOT FOUND!!!!
          {
                int nTrkTemp = nTrk;
        int nStaTemp = nSta;

        // getting tracks matched and doing the same, CONTROLLING IF MU IS RECONSTRUCTED AS A TRACK AND NOT AS A MU
        CandMatchMap::const_iterator l;
        for(l = TracksMap->begin(); l != TracksMap->end(); l++ )
          {
            const Candidate /* & Trkreco = * l -> key, */ &  Trkmc =  * l -> val;
            if (( muCand1 == & Trkmc) && (Trkmc.pt()>ptcut_) && (abs(Trkmc.eta()<etacut_)))
              {
                nMuMC++;
            nTrk++;
            break;
              }
          }
        // the same for standalone
        CandMatchMap::const_iterator n;
        for(n = StandAloneMap->begin(); n != StandAloneMap->end(); n++ )
          {
            const Candidate & Stareco = * n -> key, &  Stamc =  * n -> val;
            if ((muCand1 == &Stamc ) && (Stamc.pt()>ptcut_) && (abs(Stamc.eta()<etacut_)))
              {
            nMuMC++;
            nSta++;
            deltaR_Mu_Sta = deltaR(Stareco, *muCand1);

            //  cout<<"Ho trovato un sta reco "<<endl;
            break;
              }
          }
        // controlling if sta and trk are reconstrucetd both, and if so the get the deltaR beetween muon MC and reco sta, to controll correlation to this happening
        if ((nSta == nStaTemp + 1) && (nTrk == nTrkTemp + 1 ) )
          {
            nNotMuMatching ++;
            if ((deltaR_Mu_Sta< deltaRStacut_))
              {
            v_.push_back(deltaR_Mu_Sta) ;
                   cout << "Not matching from trk and sta matched to MC mu, to reconstruct a recoMU" << endl;
              }
          }
          }
      }
      }
  }

virtual void MCEfficiencyAnalyzer::endJob ( void  )

inlineoverridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 174 of file MC_Efficiency_Analyzer.cc.

References gather_cfg::cout, i, and p2.

  {


    cout <<"--- nMuMC == "<<nMuMC<<endl;
    cout <<"--- nMureco == "<<nMureco<<endl;
    cout <<"--- nSta == "<<nSta<<endl;
    cout <<"--- nTrk == "<<nTrk<<endl;
    cout <<"--- nNotMuMatching from a trk and sta matched to a Mu MC == "<<nNotMuMatching<<endl;
    if (nMuMC!=0 )
      {
    cout<<" effMu == "<<(double) nMureco/nMuMC<<endl;
    cout<<" effTrk == "<< (double)(nTrk + nMureco) /nMuMC<<endl;
    cout<<" effSta == "<< (double)(nSta + nMureco) / nMuMC<<endl;
      }

vector< int >::const_iterator p2;
    for (unsigned int i =0 ; i < v_.size(); ++i )
      {
    cout<<" delta R Mu Sta == "<< v_[i]<<endl;
        }
  }

Member Data Documentation

double MCEfficiencyAnalyzer::deltaRStacut_

private

Definition at line 211 of file MC_Efficiency_Analyzer.cc.

double MCEfficiencyAnalyzer::etacut_

private

Definition at line 211 of file MC_Efficiency_Analyzer.cc.

EDGetTokenT<GenParticleCollection> MCEfficiencyAnalyzer::genParticlesToken_

private

Definition at line 208 of file MC_Efficiency_Analyzer.cc.

EDGetTokenT<CandMatchMap> MCEfficiencyAnalyzer::MuonsMapToken_

private

Definition at line 205 of file MC_Efficiency_Analyzer.cc.

EDGetTokenT<CandidateCollection> MCEfficiencyAnalyzer::MuonsToken_

private

Definition at line 204 of file MC_Efficiency_Analyzer.cc.

int MCEfficiencyAnalyzer::nMuMC

private

Definition at line 212 of file MC_Efficiency_Analyzer.cc.

int MCEfficiencyAnalyzer::nMureco

private

Definition at line 212 of file MC_Efficiency_Analyzer.cc.

int MCEfficiencyAnalyzer::nNotMuMatching

private

Definition at line 212 of file MC_Efficiency_Analyzer.cc.

int MCEfficiencyAnalyzer::nSta

private

Definition at line 212 of file MC_Efficiency_Analyzer.cc.

int MCEfficiencyAnalyzer::nTrk

private

Definition at line 212 of file MC_Efficiency_Analyzer.cc.

double MCEfficiencyAnalyzer::ptcut_

private

Definition at line 211 of file MC_Efficiency_Analyzer.cc.

EDGetTokenT<CandMatchMap> MCEfficiencyAnalyzer::StandAloneMapToken_

private

Definition at line 210 of file MC_Efficiency_Analyzer.cc.

EDGetTokenT<CandidateCollection> MCEfficiencyAnalyzer::StandAloneToken_

private

Definition at line 209 of file MC_Efficiency_Analyzer.cc.

EDGetTokenT<CandMatchMap> MCEfficiencyAnalyzer::TracksMapToken_

private

Definition at line 207 of file MC_Efficiency_Analyzer.cc.

EDGetTokenT<CandidateCollection> MCEfficiencyAnalyzer::TracksToken_

private

Definition at line 206 of file MC_Efficiency_Analyzer.cc.

vector<double> MCEfficiencyAnalyzer::v_

private

Definition at line 213 of file MC_Efficiency_Analyzer.cc.

EDGetTokenT<CandidateCollection> MCEfficiencyAnalyzer::zMuMuToken_

private

Definition at line 203 of file MC_Efficiency_Analyzer.cc.