dimuonHLTFilterCheck Class Reference

Inheritance diagram for dimuonHLTFilterCheck:

Public Member Functions
	dimuonHLTFilterCheck (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
vector< int >	nDimuonsByType (const Handle< CandidateView > d)
vector< int >	nMuonsByType (const Handle< CandidateView > d)

Private Attributes
EDGetTokenT< CandidateView >	anyDimuonToken
int	counterMatrix [5][5]
EDGetTokenT< CandidateView >	muonToken
EDGetTokenT< CandidateView >	tracksToken
EDGetTokenT< TriggerResults >	triggerResultsToken

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 42 of file dimuonHLTFilterCheck.cc.

Constructor & Destructor Documentation

dimuonHLTFilterCheck::dimuonHLTFilterCheck

(

const edm::ParameterSet &

pset

)

Definition at line 74 of file dimuonHLTFilterCheck.cc.

References counterMatrix, and mps_fire::i.

                                                                   :
  // trigger results
  triggerResultsToken(consumes<TriggerResults>(pset.getParameter<InputTag>("triggerResults"))),
  tracksToken(consumes<CandidateView>(pset.getParameter<InputTag>("tracks"))),
  muonToken(consumes<CandidateView>(pset.getParameter<InputTag>("muons"))),
  anyDimuonToken(consumes<CandidateView>(pset.getParameter<InputTag>("anyDimuon")))
{
  Service<TFileService> fs;


// general histograms

// general counters
  for (int i=0; i<5; i++) {
    for (int j=0;j<5;j++) {
      counterMatrix[i][j]=0;
    }
  }
}

Member Function Documentation

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

overrideprivate

Definition at line 94 of file dimuonHLTFilterCheck.cc.

References edm::HLTGlobalStatus::accept(), anyDimuonToken, counterMatrix, gather_cfg::cout, reco::Candidate::eta(), reco::Candidate::get(), mps_fire::i, reco::Candidate::isGlobalMuon(), extraflags_cff::muons, muons_cff::muonTable, muonToken, nDimuonsByType(), nMuonsByType(), edm::Handle< T >::product(), reco::Candidate::pt(), edm::View< T >::refAt(), edm::HLTGlobalStatus::size(), edm::View< T >::size(), edm::HLTGlobalStatus::state(), l1t::tracks, tracksToken, edm::TriggerNames::triggerName(), edm::TriggerNames::triggerNames(), triggerResults, triggerResultsToken, and edm::HLTGlobalStatus::wasrun().

                                                                              {

  Handle<TriggerResults> triggerResults;   // trigger response
  Handle<CandidateView> anyDimuon;         // any dimuon pair
  Handle<CandidateView> tracks;            // any track
  Handle<CandidateView> muons;             // any muon

  event.getByToken( triggerResultsToken, triggerResults );
  event.getByToken( anyDimuonToken, anyDimuon );
  event.getByToken( tracksToken, tracks );
  event.getByToken( muonToken, muons );

  const edm::TriggerNames & triggerNames = event.triggerNames(*triggerResults);

  // map of MU triggers of interest
  map<string,int> dimuonHLT_triggers;
  dimuonHLT_triggers.insert(make_pair("HLT1MuonPrescalePt3",0));
  dimuonHLT_triggers.insert(make_pair("HLT1MuonPrescalePt7x7",1));
  dimuonHLT_triggers.insert(make_pair("HLT1MuonIso",2));
  dimuonHLT_triggers.insert(make_pair("HLT1MuonNonIso15",3));
  dimuonHLT_triggers.insert(make_pair("HLT2MuonNonIso",4));
  // map of JET triggers of interest
  map<string,int> jetHLT_triggers;
  jetHLT_triggers.insert(make_pair("HLT1jetPE5",0));
  jetHLT_triggers.insert(make_pair("HLT1jetPE3",1));
  jetHLT_triggers.insert(make_pair("HLT1jetPE1",2));
  jetHLT_triggers.insert(make_pair("candHLT1jetPE7",3));

  bool trgMask[5];
  for (int i=0; i<5; i++) trgMask[i] = false;

  // table of possible dimuons
  string dimuonTableNames[10];
  dimuonTableNames[0] = "global-global";
  dimuonTableNames[1] = "global-trackerSta";
  dimuonTableNames[2] = "global-sta";
  dimuonTableNames[3] = "global-tracker";
  dimuonTableNames[4] = "trackerSta-trackerSta";
  dimuonTableNames[5] = "trackerSta-sta";
  dimuonTableNames[6] = "trackerSta-tracker";
  dimuonTableNames[7] = "sta-sta";
  dimuonTableNames[8] = "sta-tracker";
  dimuonTableNames[9] = "tracker-tracker";
  // table of possible muons
  string muonTableNames[10];
  muonTableNames[0] = "global";
  muonTableNames[1] = "trackerSta";
  muonTableNames[2] = "sta";
  muonTableNames[3] = "tracker";

  cout << "-------------------NEW event---------------------------" << endl;
  // check the dimuon reconstruction
  vector<int> dimuonTable;
  dimuonTable = nDimuonsByType(anyDimuon);
  // check the muon reconstruction
  vector<int> muonTable;
  muonTable = nMuonsByType(muons);

  if ( triggerResults.product()->wasrun() ){
      //      cout << "at least one path out of " << triggerResults.product()->size() << " ran? " << triggerResults.product()->wasrun() << endl;

    if ( triggerResults.product()->accept() )
      {
    //  cout << endl << "at least one path accepted? " << triggerResults.product()->accept() << endl;

    const unsigned int n_TriggerResults( triggerResults.product()->size() );
    for ( unsigned int itrig( 0 ); itrig < n_TriggerResults; ++itrig )
      {
        if ( triggerResults.product()->accept( itrig ) )
          {
        map<string,int>::iterator iterMuHLT = dimuonHLT_triggers.find(triggerNames.triggerName( itrig ));
        if (iterMuHLT != dimuonHLT_triggers.end()) {
          cout << "ecco la chiave Mu HLT " << (*iterMuHLT).second << endl;
          if (triggerResults.product()->state( itrig )==1) trgMask[(*iterMuHLT).second] = true;
          }   // end if key found
        map<string,int>::iterator iterjetHLT = jetHLT_triggers.find(triggerNames.triggerName( itrig ));
        if (iterjetHLT != jetHLT_triggers.end()) {
          cout << "ecco la chiave jet HLT " << (*iterjetHLT).second << endl;
          }   // end if key found

          }  // end if trigger accepted
      }   // end loop on triggerResults
      } // end if at least one triggerResult accepted
  }  // end if wasRun
  if ( muonTable[0]>1) {
    for(unsigned int i = 0; i < muons->size(); ++i) { //loop on candidates
      const Candidate & muCand = (*muons)[i]; //the candidate
      CandidateBaseRef muCandRef = muons->refAt(i);
      TrackRef muStaComponentRef = muCand.get<TrackRef,reco::StandAloneMuonTag>();  // standalone part of muon
      TrackRef muTrkComponentRef = muCand.get<TrackRef>();  // track part of muon
      if (muCandRef->isGlobalMuon()) {
    cout << "muCand : " << i << "  pt " << muCandRef->pt() << "  eta " << muCandRef->eta() << endl;
    cout << "muCandStaComponent : " << i << "  pt " << muStaComponentRef->pt() << "  eta " << muStaComponentRef->eta() << endl;
      }
    }
  }

  // fill counterMatrix
  for (int i=0; i<5; i++) {
    for (int j=0; j<5; j++) {
      if (trgMask[i] && trgMask[j]) counterMatrix[i][j]++;
    }
  }
}  // end analyze

void dimuonHLTFilterCheck::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 267 of file dimuonHLTFilterCheck.cc.

References counterMatrix, gather_cfg::cout, DEFINE_FWK_MODULE, and mps_fire::i.

                                  {

  cout << "------------------------------------  Counters  --------------------------------" << endl;
  for (int i=0; i<5; i++) {
    cout << "trg " << i << ": ";
    for (int j=0; j<5; j++) {
      cout << counterMatrix[i][j] << " ";
    }
    cout << endl;
  }
}

vector< int > dimuonHLTFilterCheck::nDimuonsByType ( const Handle< CandidateView >  d )

private

Definition at line 199 of file dimuonHLTFilterCheck.cc.

References reco::Candidate::daughter(), mps_fire::i, reco::Candidate::isGlobalMuon(), reco::Candidate::isStandAloneMuon(), reco::Candidate::isTrackerMuon(), edm::View< T >::refAt(), and edm::View< T >::size().

Referenced by analyze().

                                                                              {
  vector<int> n_;
  int nCat = 10;    // number of dimuon categories (0 = glb-glb, 1 = glb-trkSta, 2 = glb-sta, 3 = glb-trk, 4 = trkSta-trkSta, 5 = trkSta-sta, 6 = trkSta-trk, 7 = sta-sta, 8 = sta-trk, 9  trk-trk)
  // reset vector
  for (int i=0; i<nCat; i++) n_.push_back(0);
  for(unsigned int i = 0; i < d->size(); ++i) { //loop on candidates
    //      const Candidate & dCand = (*d)[i]; //the candidate
    bool mu0global = false;
    bool mu0trackerSta = false;
    bool mu0sta = false;
    bool mu0tracker = false;
    bool mu1global = false;
    bool mu1trackerSta = false;
    bool mu1sta = false;
    bool mu1tracker = false;
    CandidateBaseRef dCandRef = d->refAt(i);
    const Candidate * mu0 = dCandRef->daughter(0);
    const Candidate * mu1 = dCandRef->daughter(1);
    if (mu0->isGlobalMuon()) mu0global=true;
    if (! mu0->isGlobalMuon() && mu0->isStandAloneMuon() && mu0->isTrackerMuon()) mu0trackerSta=true;
    if (! mu0->isGlobalMuon() && mu0->isStandAloneMuon() && ! mu0->isTrackerMuon()) mu0sta=true;
    if (! mu0->isGlobalMuon() && ! mu0->isStandAloneMuon() && mu0->isTrackerMuon()) mu0tracker=true;
    if (mu1->isGlobalMuon()) mu1global=true;
    if (! mu1->isGlobalMuon() && mu1->isStandAloneMuon() && mu1->isTrackerMuon()) mu1trackerSta=true;
    if (! mu1->isGlobalMuon() && mu1->isStandAloneMuon() && ! mu1->isTrackerMuon()) mu1sta=true;
    if (! mu1->isGlobalMuon() && ! mu1->isStandAloneMuon() && mu1->isTrackerMuon()) mu1tracker=true;

    if (mu0global && mu1global) n_[0]++;
    if ( (mu0global && mu1trackerSta) || (mu1global && mu0trackerSta) ) n_[1]++;
    if ( (mu0global && mu1sta) || (mu1global && mu0sta) ) n_[2]++;
    if ( (mu0global && mu1tracker) || (mu1global && mu0tracker) ) n_[3]++;
    if (mu0trackerSta && mu1trackerSta) n_[4]++;
    if ( (mu0trackerSta && mu1sta) || (mu1trackerSta && mu0sta) ) n_[5]++;
    if ( (mu0trackerSta && mu1tracker) || (mu1trackerSta && mu0tracker) ) n_[6]++;
    if (mu0sta && mu1sta) n_[7]++;
    if ( (mu0sta && mu1tracker) || (mu1sta && mu0tracker) ) n_[8]++;
    if (mu0tracker && mu1tracker) n_[9]++;

  }
  return n_;
}

vector< int > dimuonHLTFilterCheck::nMuonsByType ( const Handle< CandidateView >  d )

private

Definition at line 241 of file dimuonHLTFilterCheck.cc.

References gather_cfg::cout, mps_fire::i, reco::Candidate::isGlobalMuon(), reco::Candidate::isStandAloneMuon(), reco::Candidate::isTrackerMuon(), edm::View< T >::refAt(), and edm::View< T >::size().

Referenced by analyze().

                                                                            {
  vector<int> n_;
  int nCat = 4;    // number of muon categories (0 = glb, 1 = trkSta, 2 = sta, 3 = trk)
  // reset vector
  for (int i=0; i<nCat; i++) n_.push_back(0);
  for(unsigned int i = 0; i < d->size(); ++i) { //loop on candidates
    //      const Candidate & dCand = (*d)[i]; //the candidate
    bool muglobal = false;
    bool mutrackerSta = false;
    bool musta = false;
    bool mutracker = false;
    CandidateBaseRef muCandRef = d->refAt(i);
    if (muCandRef->isGlobalMuon()) muglobal=true;
    if (! muCandRef->isGlobalMuon() && muCandRef->isStandAloneMuon() && muCandRef->isTrackerMuon()) mutrackerSta=true;
    if (! muCandRef->isGlobalMuon() && muCandRef->isStandAloneMuon() && ! muCandRef->isTrackerMuon()) musta=true;
    if (! muCandRef->isGlobalMuon() && ! muCandRef->isStandAloneMuon() && muCandRef->isTrackerMuon()) mutracker=true;
    cout << "muglobal " << muglobal << " mutrackserSta " << mutrackerSta << " must " << musta << " mutracker " << mutracker << endl;
    if (muglobal) n_[0]++;
    if (mutrackerSta) n_[1]++;
    if (musta) n_[2]++;
    if (mutracker) n_[3]++;

  }
  return n_;
}

Member Data Documentation

EDGetTokenT<CandidateView> dimuonHLTFilterCheck::anyDimuonToken

private

Definition at line 53 of file dimuonHLTFilterCheck.cc.

Referenced by analyze().

int dimuonHLTFilterCheck::counterMatrix[5][5]

private

Definition at line 59 of file dimuonHLTFilterCheck.cc.

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

EDGetTokenT<CandidateView> dimuonHLTFilterCheck::muonToken

private

Definition at line 52 of file dimuonHLTFilterCheck.cc.

Referenced by analyze().

EDGetTokenT<CandidateView> dimuonHLTFilterCheck::tracksToken

private

Definition at line 51 of file dimuonHLTFilterCheck.cc.

Referenced by analyze().

EDGetTokenT<TriggerResults> dimuonHLTFilterCheck::triggerResultsToken

private

Definition at line 50 of file dimuonHLTFilterCheck.cc.

Referenced by analyze().