CMS 3D CMS Logo
List of all members | Public Member Functions | Private Attributes
WMuNuSelector Class Reference
Inheritance diagram for WMuNuSelector:
edm::EDFilter edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

void beginJob () override
 
void endJob () override
 
bool filter (edm::Event &, const edm::EventSetup &) override
 
void init_histograms ()
 
 WMuNuSelector (const edm::ParameterSet &)
 
- Public Member Functions inherited from edm::EDFilter
 EDFilter ()
 
SerialTaskQueue * globalLuminosityBlocksQueue ()
 
SerialTaskQueue * globalRunsQueue ()
 
ModuleDescription const & moduleDescription () const
 
 ~EDFilter () override
 
- Public Member Functions inherited from edm::ProducerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
std::vector< edm::ProductResolverIndex > const & indiciesForPutProducts (BranchType iBranchType) const
 
 ProducerBase ()
 
std::vector< edm::ProductResolverIndex > const & putTokenIndexToProductResolverIndex () const
 
std::vector< bool > const & recordProvenanceList () const
 
void registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
 
std::function< void(BranchDescription const &)> registrationCallback () const
 used by the fwk to register list of products More...
 
void resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
 
TypeLabelList const & typeLabelList () const
 used by the fwk to register the list of products of this module More...
 
 ~ProducerBase () noexcept(false) 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 &&)=default
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
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 & operator= (EDConsumerBase &&)=default
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
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 Attributes

double acopCut_
 
edm::EDGetTokenT< reco::BeamSpot > beamSpotToken_
 
double dxyCut_
 
double eJetMin_
 
double etaCut_
 
std::map< std::string, TH1D * > h1_
 
std::map< std::string, TH2D * > h2_
 
bool isAlsoTrackerMuon_
 
bool isCombinedIso_
 
double isoCut03_
 
bool isRelativeIso_
 
edm::EDGetTokenT< edm::View< reco::Jet > > jetToken_
 
double metMax_
 
double metMin_
 
double mtMax_
 
double mtMin_
 
edm::EDGetTokenT< edm::View< reco::Muon > > muonToken_
 
const std::string muonTrig_
 
double nacop
 
double nall
 
double ncharge
 
double nid
 
double niso
 
int nJetMax_
 
double nkin
 
double nmass
 
double normalizedChi2Cut_
 
double npresel
 
double nsel
 
double ntrig
 
bool plotHistograms_
 
double ptCut_
 
double ptThrForZ1_
 
double ptThrForZ2_
 
int selectByCharge_
 
int trackerHitsCut_
 
edm::EDGetTokenT< edm::TriggerResults > trigToken_
 
edm::EDGetTokenT< reco::WMuNuCandidateCollection > WMuNuCollectionToken_
 

Additional Inherited Members

- Public Types inherited from edm::EDFilter
typedef EDFilter ModuleType
 
- Public Types inherited from edm::ProducerBase
using ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex > >
 
typedef ProductRegistryHelper::TypeLabelList TypeLabelList
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::EDFilter
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::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType > produces ()
 declare what type of product will make and with which optional label More...
 
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType > produces ()
 
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType > produces ()
 
BranchAliasSetter produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
 
template<BranchType B>
BranchAliasSetter produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
 
template<Transition B>
BranchAliasSetter produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
 
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType > produces (std::string instanceName)
 
template<class ProductType >
BranchAliasSetterT< ProductType > produces (std::string instanceName)
 
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType > produces (std::string instanceName)
 
ProducesCollector producesCollector ()
 
- Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (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)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 

Detailed Description

Definition at line 35 of file WMuNuSelector.cc.

Constructor & Destructor Documentation

◆ WMuNuSelector()

WMuNuSelector::WMuNuSelector ( const edm::ParameterSet &  cfg)

Definition at line 102 of file WMuNuSelector.cc.

103  : // Fast selection (no histograms)
104  plotHistograms_(cfg.getUntrackedParameter<bool>("plotHistograms", true)),
105 
106  // Input collections
107  trigToken_(consumes<TriggerResults>(
108  cfg.getUntrackedParameter<edm::InputTag>("TrigTag", edm::InputTag("TriggerResults::HLT")))),
109  muonToken_(consumes<View<Muon> >(cfg.getUntrackedParameter<edm::InputTag>("MuonTag", edm::InputTag("muons")))),
110  jetToken_(
111  consumes<View<Jet> >(cfg.getUntrackedParameter<edm::InputTag>("JetTag", edm::InputTag("sisCone5CaloJets")))),
112  beamSpotToken_(consumes<reco::BeamSpot>(edm::InputTag("offlineBeamSpot"))),
113  WMuNuCollectionToken_(consumes<reco::WMuNuCandidateCollection>(
114  cfg.getUntrackedParameter<edm::InputTag>("WMuNuCollectionTag", edm::InputTag("WMuNus")))),
115 
116  // Preselection cuts
117  muonTrig_(cfg.getUntrackedParameter<std::string>("MuonTrig", "HLT_Mu9")),
118  ptThrForZ1_(cfg.getUntrackedParameter<double>("PtThrForZ1", 20.)),
119  ptThrForZ2_(cfg.getUntrackedParameter<double>("PtThrForZ2", 10.)),
120  eJetMin_(cfg.getUntrackedParameter<double>("EJetMin", 999999.)),
121  nJetMax_(cfg.getUntrackedParameter<int>("NJetMax", 999999)),
122 
123  // Main cuts
124  ptCut_(cfg.getUntrackedParameter<double>("PtCut", 25.)),
125  etaCut_(cfg.getUntrackedParameter<double>("EtaCut", 2.1)),
126  isRelativeIso_(cfg.getUntrackedParameter<bool>("IsRelativeIso", true)),
127  isCombinedIso_(cfg.getUntrackedParameter<bool>("IsCombinedIso", false)),
128  isoCut03_(cfg.getUntrackedParameter<double>("IsoCut03", 0.1)),
129  mtMin_(cfg.getUntrackedParameter<double>("MtMin", 50.)),
130  mtMax_(cfg.getUntrackedParameter<double>("MtMax", 200.)),
131  metMin_(cfg.getUntrackedParameter<double>("MetMin", -999999.)),
132  metMax_(cfg.getUntrackedParameter<double>("MetMax", 999999.)),
133  acopCut_(cfg.getUntrackedParameter<double>("AcopCut", 2.)),
134 
135  // Muon quality cuts
136  dxyCut_(cfg.getUntrackedParameter<double>("DxyCut", 0.2)),
137  normalizedChi2Cut_(cfg.getUntrackedParameter<double>("NormalizedChi2Cut", 10.)),
138  trackerHitsCut_(cfg.getUntrackedParameter<int>("TrackerHitsCut", 11)),
139  isAlsoTrackerMuon_(cfg.getUntrackedParameter<bool>("IsAlsoTrackerMuon", true)),
140 
141  // W+/W- Selection
142  selectByCharge_(cfg.getUntrackedParameter<int>("SelectByCharge", 0)) {}

Member Function Documentation

◆ beginJob()

void WMuNuSelector::beginJob ( void  )
overridevirtual

Reimplemented from edm::EDFilter.

Definition at line 144 of file WMuNuSelector.cc.

144  {
145  nall = 0;
146  ntrig = 0;
147  npresel = 0;
148  ncharge = 0;
149  nkin = 0;
150  nid = 0;
151  nacop = 0;
152  niso = 0;
153  nsel = 0;
154 
155  if (plotHistograms_) {
156  edm::Service<TFileService> fs;
157  h1_["hNWCand_sel"] =
158  fs->make<TH1D>("NWCand_sel", "Nb. of WCandidates passing pre-selection (ordered by pt)", 10, 0., 10.);
159  h1_["hPtMu_sel"] = fs->make<TH1D>("ptMu_sel", "Pt mu", 100, 0., 100.);
160  h1_["hEtaMu_sel"] = fs->make<TH1D>("etaMu_sel", "Eta mu", 50, -2.5, 2.5);
161  h1_["hd0_sel"] = fs->make<TH1D>("d0_sel", "Impact parameter", 1000, -1., 1.);
162  h1_["hNHits_sel"] = fs->make<TH1D>("NumberOfValidHits_sel", "Number of Hits in Silicon", 100, 0., 100.);
163  h1_["hNormChi2_sel"] = fs->make<TH1D>("NormChi2_sel", "Chi2/ndof of global track", 1000, 0., 50.);
164  h1_["hTracker_sel"] = fs->make<TH1D>("isTrackerMuon_sel", "is Tracker Muon?", 2, 0., 2.);
165  h1_["hMET_sel"] = fs->make<TH1D>("MET_sel", "Missing Transverse Energy (GeV)", 300, 0., 300.);
166  h1_["hTMass_sel"] = fs->make<TH1D>("TMass_sel", "Rec. Transverse Mass (GeV)", 300, 0., 300.);
167  h1_["hAcop_sel"] = fs->make<TH1D>("Acop_sel", "Mu-MET acoplanarity", 50, 0., M_PI);
168  h1_["hPtSum_sel"] = fs->make<TH1D>("ptSum_sel", "Track Isolation, Sum pT (GeV)", 200, 0., 100.);
169  h1_["hPtSumN_sel"] = fs->make<TH1D>("ptSumN_sel", "Track Isolation, Sum pT/pT", 1000, 0., 10.);
170  h1_["hCal_sel"] = fs->make<TH1D>("Cal_sel", "Calorimetric isolation, HCAL+ECAL (GeV)", 200, 0., 100.);
171  h1_["hIsoTotN_sel"] = fs->make<TH1D>("isoTotN_sel", "(Sum pT + Cal)/pT", 1000, 0., 10.);
172  h1_["hIsoTot_sel"] = fs->make<TH1D>("isoTot_sel", "(Sum pT + Cal)", 200, 0., 100.);
173  h2_["hTMass_PtSum_inclusive"] = fs->make<TH2D>(
174  "TMass_PtSum_inclusive", "Rec. Transverse Mass (GeV) vs Sum pT (GeV)", 200, 0., 100., 300, 0., 300.);
175  h2_["hTMass_PtSumNorm_inclusive"] = fs->make<TH2D>(
176  "TMass_PtSumNorm_inclusive", "Rec. Transverse Mass (GeV) vs Sum Pt / Pt", 1000, 0, 10, 300, 0, 300);
177  h2_["hTMass_TotIsoNorm_inclusive"] = fs->make<TH2D>(
178  "TMass_TotIsoNorm_inclusive", "Rec. Transverse Mass (GeV) vs (Sum Pt + Cal)/Pt", 10000, 0, 10, 200, 0, 200);
179  h2_["hMET_PtSum_inclusive"] = fs->make<TH2D>(
180  "MET_PtSum_inclusive", "Missing Transverse Energy (GeV) vs Sum Pt (GeV)", 200, 0., 100., 300, 0., 300.);
181  h2_["hMET_PtSumNorm_inclusive"] = fs->make<TH2D>(
182  "MET_PtSumNorm_inclusive", "Missing Transverse Energy (GeV) vs Sum Pt/Pt", 1000, 0, 10, 300, 0, 300);
183  h2_["hMET_TotIsoNorm_inclusive"] = fs->make<TH2D>(
184  "MET_TotIsoNorm_inclusive", "Missing Transverse Energy (GeV) vs (SumPt + Cal)/Pt", 10000, 0, 10, 200, 0, 200);
185  }
186 }

References h1_, h2_, M_PI, TFileService::make(), nacop, nall, ncharge, nid, niso, nkin, npresel, nsel, ntrig, and plotHistograms_.

◆ endJob()

void WMuNuSelector::endJob ( void  )
overridevirtual

Reimplemented from edm::EDFilter.

Definition at line 188 of file WMuNuSelector.cc.

188  {
189  double all = nall;
190  double epresel = npresel / all;
191  double etrig = ntrig / all;
192  double ekin = nkin / all;
193  double eid = nid / all;
194  double eacop = nacop / all;
195  double eiso = niso / all;
196  double esel = nsel / all;
197 
198  LogVerbatim("") << "\n>>>>>> W SELECTION SUMMARY BEGIN >>>>>>>>>>>>>>>";
199  LogVerbatim("") << "Total number of events analyzed: " << nall << " [events]";
200  LogVerbatim("") << "Total number of events triggered: " << ntrig << " [events]";
201  LogVerbatim("") << "Total number of events pre-selected: " << npresel << " [events]";
202  LogVerbatim("") << "Total number of events after kinematic cuts: " << nkin << " [events]";
203  LogVerbatim("") << "Total number of events after Muon ID cuts: " << nid << " [events]";
204  LogVerbatim("") << "Total number of events after Acop cut: " << nacop << " [events]";
205  LogVerbatim("") << "Total number of events after iso cut: " << niso << " [events]";
206  LogVerbatim("") << "Total number of events selected: " << nsel << " [events]";
207  LogVerbatim("") << "Efficiencies:";
208  LogVerbatim("") << "Trigger Efficiency: "
209  << "(" << setprecision(4) << etrig * 100. << " +/- " << setprecision(2)
210  << sqrt(etrig * (1 - etrig) / all) * 100. << ")%";
211  LogVerbatim("") << "Pre-Selection Efficiency: "
212  << "(" << setprecision(4) << epresel * 100. << " +/- " << setprecision(2)
213  << sqrt(epresel * (1 - epresel) / all) * 100. << ")%";
214  LogVerbatim("") << "Pt, Eta Selection Efficiency: "
215  << "(" << setprecision(4) << ekin * 100. << " +/- " << setprecision(2)
216  << sqrt(ekin * (1 - ekin) / all) * 100. << ")%";
217  LogVerbatim("") << "MuonID Efficiency: "
218  << "(" << setprecision(4) << eid * 100. << " +/- " << setprecision(2)
219  << sqrt(eid * (1 - eid) / all) * 100. << ")%";
220  LogVerbatim("") << "Acop Efficiency: "
221  << "(" << setprecision(4) << eacop * 100. << " +/- " << setprecision(2)
222  << sqrt(eacop * (1 - eacop) / all) * 100. << ")%";
223  LogVerbatim("") << "Iso Efficiency: "
224  << "(" << setprecision(4) << eiso * 100. << " +/- " << setprecision(2)
225  << sqrt(eiso * (1 - eiso) / all) * 100. << ")%";
226  LogVerbatim("") << "Selection Efficiency: "
227  << "(" << setprecision(4) << esel * 100. << " +/- " << setprecision(2)
228  << sqrt(esel * (1 - esel) / nall) * 100. << ")%";
229 
230  if (fabs(selectByCharge_) == 1) {
231  esel = nsel / ncharge;
232  LogVerbatim("") << "\n>>>>>> W+(-) SELECTION >>>>>>>>>>>>>>>";
233  LogVerbatim("") << "Total number of W+(-) events pre-selected: " << ncharge << " [events]";
234  LogVerbatim("") << "Total number of events selected: " << nsel << " [events]";
235  LogVerbatim("") << "Selection Efficiency only W+(-): "
236  << "(" << setprecision(4) << esel * 100. << " +/- " << setprecision(2)
237  << sqrt(esel * (1 - esel) / ncharge) * 100. << ")%";
238  }
239  LogVerbatim("") << ">>>>>> W SELECTION SUMMARY END >>>>>>>>>>>>>>>\n";
240 }

References python.cmstools::all(), runTauDisplay::eid, nacop, nall, ncharge, nid, niso, nkin, npresel, nsel, ntrig, selectByCharge_, and mathSSE::sqrt().

◆ filter()

bool WMuNuSelector::filter ( edm::Event &  ev,
const edm::EventSetup &   
)
overridevirtual

Implements edm::EDFilter.

Definition at line 242 of file WMuNuSelector.cc.

242  {
243  nall++;
244 
245  // Repeat Pre-Selection Cuts just in case...
246  // Muon collection
247  Handle<View<Muon> > muonCollection;
248  if (!ev.getByToken(muonToken_, muonCollection)) {
249  LogError("") << ">>> Muon collection does not exist !!!";
250  return false;
251  }
252  unsigned int muonCollectionSize = muonCollection->size();
253 
254  // Trigger
255  Handle<TriggerResults> triggerResults;
256  if (!ev.getByToken(trigToken_, triggerResults)) {
257  LogError("") << ">>> TRIGGER collection does not exist !!!";
258  return false;
259  }
260  const edm::TriggerNames& triggerNames = ev.triggerNames(*triggerResults);
261  bool trigger_fired = false;
262  int itrig1 = triggerNames.triggerIndex(muonTrig_);
263  if (triggerResults->accept(itrig1))
264  trigger_fired = true;
265  LogTrace("") << ">>> Trigger bit: " << trigger_fired << " (" << muonTrig_ << ")";
266 
267  // Loop to reject/control Z->mumu is done separately
268  unsigned int nmuonsForZ1 = 0;
269  unsigned int nmuonsForZ2 = 0;
270  for (unsigned int i = 0; i < muonCollectionSize; i++) {
271  const Muon& mu = muonCollection->at(i);
272  if (!mu.isGlobalMuon())
273  continue;
274  double pt = mu.pt();
275  if (pt > ptThrForZ1_)
276  nmuonsForZ1++;
277  if (pt > ptThrForZ2_)
278  nmuonsForZ2++;
279  }
280  LogTrace("") << "> Z rejection: muons above " << ptThrForZ1_ << " [GeV]: " << nmuonsForZ1;
281  LogTrace("") << "> Z rejection: muons above " << ptThrForZ2_ << " [GeV]: " << nmuonsForZ2;
282 
283  // Jet collection
284  Handle<View<Jet> > jetCollection;
285  if (!ev.getByToken(jetToken_, jetCollection)) {
286  LogError("") << ">>> JET collection does not exist !!!";
287  return false;
288  }
289  unsigned int jetCollectionSize = jetCollection->size();
290  int njets = 0;
291  for (unsigned int i = 0; i < jetCollectionSize; i++) {
292  const Jet& jet = jetCollection->at(i);
293  if (jet.et() > eJetMin_)
294  njets++;
295  }
296  LogTrace("") << ">>> Total number of jets: " << jetCollectionSize;
297  LogTrace("") << ">>> Number of jets above " << eJetMin_ << " [GeV]: " << njets;
298 
299  // Beam spot
300  Handle<reco::BeamSpot> beamSpotHandle;
301  if (!ev.getByToken(beamSpotToken_, beamSpotHandle)) {
302  LogTrace("") << ">>> No beam spot found !!!";
303  return false;
304  }
305 
306  // Get WMuNu candidates from file:
307 
308  Handle<reco::WMuNuCandidateCollection> WMuNuCollection;
309  if (!ev.getByToken(WMuNuCollectionToken_, WMuNuCollection)) {
310  LogTrace("") << ">>> WMuNu not found !!!";
311  return false;
312  }
313 
314  if (WMuNuCollection->empty()) {
315  LogTrace("") << "No WMuNu Candidates in the Event!";
316  return false;
317  }
318  if (WMuNuCollection->size() > 1) {
319  LogTrace("") << "This event contains more than one W Candidate";
320  }
321 
322  // W->mu nu selection criteria
323 
324  LogTrace("") << "> WMuNu Candidate with: ";
325  const WMuNuCandidate& WMuNu = WMuNuCollection->at(0);
326  // WMuNuCandidates are ordered by Pt!
327  // The Inclusive Selection WMuNu Candidate is the first one
328 
329  const reco::Muon& mu = WMuNu.getMuon();
330  const reco::MET& met = WMuNu.getNeutrino();
331  if (plotHistograms_) {
332  h1_["hNWCand_sel"]->Fill(WMuNuCollection->size());
333  }
334 
335  // Preselection cuts:
336 
337  if (!trigger_fired) {
338  LogTrace("") << "Event did not fire the Trigger";
339  return false;
340  }
341  ntrig++;
342 
343  if (nmuonsForZ1 >= 1 && nmuonsForZ2 >= 2) {
344  LogTrace("") << "Z Candidate!!";
345  return false;
346  }
347  if (njets > nJetMax_) {
348  LogTrace("") << "NJets > threshold";
349  return false;
350  }
351 
352  npresel++;
353 
354  // Select Ws by charge:
355 
356  if (selectByCharge_ * WMuNu.charge() == -1) {
357  ncharge++;
358  }
359 
360  // W->mu nu selection criteria
361 
362  if (!mu.isGlobalMuon())
363  return false;
364 
365  reco::TrackRef gm = mu.globalTrack();
366  //reco::TrackRef tk = mu.innerTrack();
367 
368  // Pt,eta cuts
369  double pt = mu.pt();
370  double eta = mu.eta();
371  LogTrace("") << "\t... Muon pt, eta: " << pt << " [GeV], " << eta;
372  if (plotHistograms_) {
373  h1_["hPtMu_sel"]->Fill(pt);
374  }
375  if (pt < ptCut_)
376  return false;
377  if (plotHistograms_) {
378  h1_["hEtaMu_sel"]->Fill(eta);
379  }
380  if (fabs(eta) > etaCut_)
381  return false;
382 
383  nkin++;
384 
385  // d0, chi2, nhits quality cuts
386  double dxy = gm->dxy(beamSpotHandle->position());
387  double normalizedChi2 = gm->normalizedChi2();
388  LogTrace("") << "Im here" << endl;
389  double trackerHits = gm->hitPattern().numberOfValidTrackerHits();
390  LogTrace("") << "\t... Muon dxy, normalizedChi2, trackerHits, isTrackerMuon?: " << dxy << " [cm], " << normalizedChi2
391  << ", " << trackerHits << ", " << mu.isTrackerMuon();
392 
393  if (plotHistograms_) {
394  h1_["hd0_sel"]->Fill(dxy);
395  }
396  if (!muon::isGoodMuon(mu, muon::GlobalMuonPromptTight))
397  return false;
398  if (plotHistograms_) {
399  h1_["hNormChi2_sel"]->Fill(normalizedChi2);
400  }
401  if (normalizedChi2 > normalizedChi2Cut_)
402  return false;
403  if (plotHistograms_) {
404  h1_["hNHits_sel"]->Fill(trackerHits);
405  }
406  if (trackerHits < trackerHitsCut_)
407  return false;
408  if (plotHistograms_) {
409  h1_["hTracker_sel"]->Fill(mu.isTrackerMuon());
410  }
411  if (!mu.isTrackerMuon())
412  return false;
413 
414  nid++;
415 
416  // Acoplanarity cuts
417  double acop = WMuNu.acop();
418  LogTrace("") << "\t... acoplanarity: " << acop;
419 
420  // Isolation cuts
421  double SumPt = mu.isolationR03().sumPt;
422  double isovar = SumPt;
423  double Cal = mu.isolationR03().emEt + mu.isolationR03().hadEt;
424  if (isCombinedIso_)
425  isovar += Cal;
426  if (plotHistograms_) {
427  h1_["hPtSum_sel"]->Fill(SumPt);
428  h1_["hPtSumN_sel"]->Fill(SumPt / pt);
429  h1_["hCal_sel"]->Fill(Cal);
430  h1_["hIsoTot_sel"]->Fill((SumPt + Cal));
431  h1_["hIsoTotN_sel"]->Fill((SumPt + Cal) / pt);
432  }
433 
434  if (isRelativeIso_)
435  isovar /= pt;
436  bool iso = (isovar <= isoCut03_);
437  LogTrace("") << "\t... isolation value" << isovar << ", isolated? " << iso;
438 
439  double met_et = met.pt();
440  LogTrace("") << "\t... Met pt: " << WMuNu.getNeutrino().pt() << "[GeV]";
441 
442  double massT = WMuNu.massT();
443  double w_et = WMuNu.eT();
444 
445  LogTrace("") << "\t... W mass, W_et, W_px, W_py: " << massT << ", " << w_et << ", " << WMuNu.px() << ", "
446  << WMuNu.py() << " [GeV]";
447 
448  // Plot 2D Histograms before final cuts
449  if (plotHistograms_ && acop < acopCut_) {
450  h2_["hTMass_PtSum_inclusive"]->Fill(SumPt, massT);
451  h2_["hTMass_PtSumNorm_inclusive"]->Fill(SumPt / pt, massT);
452  h2_["hTMass_TotIsoNorm_inclusive"]->Fill((SumPt + Cal) / pt, massT);
453  h2_["hMET_PtSum_inclusive"]->Fill(SumPt, met_et);
454  h2_["hMET_PtSumNorm_inclusive"]->Fill(SumPt / pt, met_et);
455  h2_["hMET_TotIsoNorm_inclusive"]->Fill((SumPt + Cal) / pt, met_et);
456  }
457 
458  if (!iso)
459  return false;
460 
461  niso++;
462 
463  if (plotHistograms_) {
464  h1_["hAcop_sel"]->Fill(acop);
465  }
466  if (acop >= acopCut_)
467  return false;
468 
469  nacop++;
470 
471  if (plotHistograms_) {
472  h1_["hMET_sel"]->Fill(met_et);
473  h1_["hTMass_sel"]->Fill(massT);
474  }
475 
476  if (massT <= mtMin_ || massT >= mtMax_)
477  return false;
478  if (met_et <= metMin_ || met_et >= metMax_)
479  return false;
480 
481  LogTrace("") << ">>>> Event ACCEPTED";
482  nsel++;
483 
484  // (To be continued ;-) )
485 
486  return true;
487 }

References reco::WMuNuCandidate::acop(), acopCut_, beamSpotToken_, reco::LeafCandidate::charge(), PVValHelper::dxy, eJetMin_, reco::WMuNuCandidate::eT(), PVValHelper::eta, etaCut_, ev, reco::WMuNuCandidate::getMuon(), reco::WMuNuCandidate::getNeutrino(), muon::GlobalMuonPromptTight, h1_, h2_, mps_fire::i, isCombinedIso_, muon::isGoodMuon(), isoCut03_, isRelativeIso_, metsig::jet, jetfilter_cfi::jetCollection, jetToken_, LogTrace, reco::WMuNuCandidate::massT(), BTaggingMonitor_cfi::met, metMax_, mtMax_, amptDefaultParameters_cff::mu, pdwgLeptonRecoSkim_cfi::muonCollection, muonToken_, muonTrig_, nacop, nall, ncharge, nid, niso, nJetMax_, BTaggingMonitoring_cff::njets, nkin, TrackCollections2monitor_cff::normalizedChi2, normalizedChi2Cut_, npresel, nsel, ntrig, plotHistograms_, reco::BeamSpot::position(), DiDispStaMuonMonitor_cfi::pt, reco::LeafCandidate::pt(), ptCut_, ptThrForZ1_, ptThrForZ2_, reco::LeafCandidate::px(), reco::LeafCandidate::py(), selectByCharge_, trackerHitsCut_, L1TEGammaOffline_cfi::triggerNames, triggerResults, trigToken_, and WMuNuCollectionToken_.

◆ init_histograms()

void WMuNuSelector::init_histograms ( )

Member Data Documentation

◆ acopCut_

double WMuNuSelector::acopCut_
private

Definition at line 64 of file WMuNuSelector.cc.

Referenced by filter().

◆ beamSpotToken_

edm::EDGetTokenT<reco::BeamSpot> WMuNuSelector::beamSpotToken_
private

Definition at line 48 of file WMuNuSelector.cc.

Referenced by filter().

◆ dxyCut_

double WMuNuSelector::dxyCut_
private

Definition at line 66 of file WMuNuSelector.cc.

◆ eJetMin_

double WMuNuSelector::eJetMin_
private

Definition at line 53 of file WMuNuSelector.cc.

Referenced by filter().

◆ etaCut_

double WMuNuSelector::etaCut_
private

Definition at line 56 of file WMuNuSelector.cc.

Referenced by filter().

◆ h1_

std::map<std::string, TH1D*> WMuNuSelector::h1_
private

Definition at line 79 of file WMuNuSelector.cc.

Referenced by beginJob(), and filter().

◆ h2_

std::map<std::string, TH2D*> WMuNuSelector::h2_
private

Definition at line 80 of file WMuNuSelector.cc.

Referenced by beginJob(), and filter().

◆ isAlsoTrackerMuon_

bool WMuNuSelector::isAlsoTrackerMuon_
private

Definition at line 69 of file WMuNuSelector.cc.

◆ isCombinedIso_

bool WMuNuSelector::isCombinedIso_
private

Definition at line 58 of file WMuNuSelector.cc.

Referenced by filter().

◆ isoCut03_

double WMuNuSelector::isoCut03_
private

Definition at line 59 of file WMuNuSelector.cc.

Referenced by filter().

◆ isRelativeIso_

bool WMuNuSelector::isRelativeIso_
private

Definition at line 57 of file WMuNuSelector.cc.

Referenced by filter().

◆ jetToken_

edm::EDGetTokenT<edm::View<reco::Jet> > WMuNuSelector::jetToken_
private

Definition at line 47 of file WMuNuSelector.cc.

Referenced by filter().

◆ metMax_

double WMuNuSelector::metMax_
private

Definition at line 63 of file WMuNuSelector.cc.

Referenced by filter().

◆ metMin_

double WMuNuSelector::metMin_
private

Definition at line 62 of file WMuNuSelector.cc.

◆ mtMax_

double WMuNuSelector::mtMax_
private

Definition at line 61 of file WMuNuSelector.cc.

Referenced by filter().

◆ mtMin_

double WMuNuSelector::mtMin_
private

Definition at line 60 of file WMuNuSelector.cc.

◆ muonToken_

edm::EDGetTokenT<edm::View<reco::Muon> > WMuNuSelector::muonToken_
private

Definition at line 46 of file WMuNuSelector.cc.

Referenced by filter().

◆ muonTrig_

const std::string WMuNuSelector::muonTrig_
private

Definition at line 50 of file WMuNuSelector.cc.

Referenced by filter().

◆ nacop

double WMuNuSelector::nacop
private

Definition at line 77 of file WMuNuSelector.cc.

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

◆ nall

double WMuNuSelector::nall
private

Definition at line 73 of file WMuNuSelector.cc.

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

◆ ncharge

double WMuNuSelector::ncharge
private

Definition at line 76 of file WMuNuSelector.cc.

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

◆ nid

double WMuNuSelector::nid
private

Definition at line 77 of file WMuNuSelector.cc.

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

◆ niso

double WMuNuSelector::niso
private

Definition at line 77 of file WMuNuSelector.cc.

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

◆ nJetMax_

int WMuNuSelector::nJetMax_
private

Definition at line 54 of file WMuNuSelector.cc.

Referenced by filter().

◆ nkin

double WMuNuSelector::nkin
private

Definition at line 77 of file WMuNuSelector.cc.

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

◆ nmass

double WMuNuSelector::nmass
private

Definition at line 77 of file WMuNuSelector.cc.

◆ normalizedChi2Cut_

double WMuNuSelector::normalizedChi2Cut_
private

Definition at line 67 of file WMuNuSelector.cc.

Referenced by filter().

◆ npresel

double WMuNuSelector::npresel
private

Definition at line 74 of file WMuNuSelector.cc.

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

◆ nsel

double WMuNuSelector::nsel
private

Definition at line 75 of file WMuNuSelector.cc.

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

◆ ntrig

double WMuNuSelector::ntrig
private

Definition at line 74 of file WMuNuSelector.cc.

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

◆ plotHistograms_

bool WMuNuSelector::plotHistograms_
private

Definition at line 44 of file WMuNuSelector.cc.

Referenced by beginJob(), and filter().

◆ ptCut_

double WMuNuSelector::ptCut_
private

Definition at line 55 of file WMuNuSelector.cc.

Referenced by badGlobalMuonTaggerFWLite.BadGlobalMuonTagger::badMuons(), and filter().

◆ ptThrForZ1_

double WMuNuSelector::ptThrForZ1_
private

Definition at line 51 of file WMuNuSelector.cc.

Referenced by filter().

◆ ptThrForZ2_

double WMuNuSelector::ptThrForZ2_
private

Definition at line 52 of file WMuNuSelector.cc.

Referenced by filter().

◆ selectByCharge_

int WMuNuSelector::selectByCharge_
private

Definition at line 71 of file WMuNuSelector.cc.

Referenced by endJob(), and filter().

◆ trackerHitsCut_

int WMuNuSelector::trackerHitsCut_
private

Definition at line 68 of file WMuNuSelector.cc.

Referenced by filter().

◆ trigToken_

edm::EDGetTokenT<edm::TriggerResults> WMuNuSelector::trigToken_
private

Definition at line 45 of file WMuNuSelector.cc.

Referenced by filter().

◆ WMuNuCollectionToken_

edm::EDGetTokenT<reco::WMuNuCandidateCollection> WMuNuSelector::WMuNuCollectionToken_
private

Definition at line 49 of file WMuNuSelector.cc.

Referenced by filter().

WMuNuSelector::dxyCut_
double dxyCut_
Definition: WMuNuSelector.cc:66
WMuNuSelector::metMin_
double metMin_
Definition: WMuNuSelector.cc:62
mps_fire.i
i
Definition: mps_fire.py:355
WMuNuSelector::nall
double nall
Definition: WMuNuSelector.cc:73
WMuNuSelector::trackerHitsCut_
int trackerHitsCut_
Definition: WMuNuSelector.cc:68
DiDispStaMuonMonitor_cfi.pt
pt
Definition: DiDispStaMuonMonitor_cfi.py:39
amptDefaultParameters_cff.mu
mu
Definition: amptDefaultParameters_cff.py:16
trackerHits
Definition: trackerHits.py:1
reco::WMuNuCandidate::acop
double acop() const
Definition: WMuNuCandidate.cc:42
muon::GlobalMuonPromptTight
Definition: MuonSelectors.h:26
triggerResults
static const std::string triggerResults
Definition: EdmProvDump.cc:45
WMuNuSelector::nkin
double nkin
Definition: WMuNuSelector.cc:77
reco::LeafCandidate::pt
double pt() const final
transverse momentum
Definition: LeafCandidate.h:146
WMuNuSelector::selectByCharge_
int selectByCharge_
Definition: WMuNuSelector.cc:71
WMuNuSelector::isAlsoTrackerMuon_
bool isAlsoTrackerMuon_
Definition: WMuNuSelector.cc:69
edm::Handle
Definition: AssociativeIterator.h:50
python.cmstools.all
def all(container)
workaround iterator generators for ROOT classes
Definition: cmstools.py:26
Muon
Definition: Muon.py:1
reco::Muon
Definition: Muon.h:27
edm::Ref< TrackCollection >
BTaggingMonitor_cfi.met
met
Definition: BTaggingMonitor_cfi.py:84
WMuNuSelector::ptThrForZ2_
double ptThrForZ2_
Definition: WMuNuSelector.cc:52
reco::MET
Definition: MET.h:41
WMuNuSelector::muonTrig_
const std::string muonTrig_
Definition: WMuNuSelector.cc:50
muon::isGoodMuon
bool isGoodMuon(const reco::Muon &muon, SelectionType type, reco::Muon::ArbitrationType arbitrationType=reco::Muon::SegmentAndTrackArbitration)
main GoodMuon wrapper call
Definition: MuonSelectors.cc:664
WMuNuSelector::WMuNuCollectionToken_
edm::EDGetTokenT< reco::WMuNuCandidateCollection > WMuNuCollectionToken_
Definition: WMuNuSelector.cc:49
reco::WMuNuCandidate::massT
double massT() const
Definition: WMuNuCandidate.cc:31
Jet
Definition: Jet.py:1
reco::LeafCandidate::py
double py() const final
y coordinate of momentum vector
Definition: LeafCandidate.h:142
L1TEGammaOffline_cfi.triggerNames
triggerNames
Definition: L1TEGammaOffline_cfi.py:40
PVValHelper::eta
Definition: PVValidationHelpers.h:69
reco::WMuNuCandidate::getMuon
const reco::Muon & getMuon() const
Definition: WMuNuCandidate.h:38
WMuNuSelector::h2_
std::map< std::string, TH2D * > h2_
Definition: WMuNuSelector.cc:80
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
WMuNuSelector::ncharge
double ncharge
Definition: WMuNuSelector.cc:76
WMuNuSelector::metMax_
double metMax_
Definition: WMuNuSelector.cc:63
reco::BeamSpot::position
const Point & position() const
position
Definition: BeamSpot.h:59
jetfilter_cfi.jetCollection
jetCollection
Definition: jetfilter_cfi.py:4
WMuNuSelector::isCombinedIso_
bool isCombinedIso_
Definition: WMuNuSelector.cc:58
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
reco::WMuNuCandidate::getNeutrino
const reco::MET & getNeutrino() const
Definition: WMuNuCandidate.h:39
edm::View
Definition: CaloClusterFwd.h:14
edm::LogError
Definition: MessageLogger.h:183
WMuNuSelector::mtMax_
double mtMax_
Definition: WMuNuSelector.cc:61
WMuNuSelector::isRelativeIso_
bool isRelativeIso_
Definition: WMuNuSelector.cc:57
WMuNuSelector::muonToken_
edm::EDGetTokenT< edm::View< reco::Muon > > muonToken_
Definition: WMuNuSelector.cc:46
edm::Service< TFileService >
runTauDisplay.eid
eid
Definition: runTauDisplay.py:298
M_PI
#define M_PI
Definition: BXVectorInputProducer.cc:50
edm::LogVerbatim
Definition: MessageLogger.h:297
WMuNuSelector::ptCut_
double ptCut_
Definition: WMuNuSelector.cc:55
reco::LeafCandidate::charge
int charge() const final
electric charge
Definition: LeafCandidate.h:106
WMuNuSelector::nsel
double nsel
Definition: WMuNuSelector.cc:75
WMuNuSelector::nacop
double nacop
Definition: WMuNuSelector.cc:77
TrackCollections2monitor_cff.normalizedChi2
normalizedChi2
Definition: TrackCollections2monitor_cff.py:247
pdwgLeptonRecoSkim_cfi.muonCollection
muonCollection
Definition: pdwgLeptonRecoSkim_cfi.py:7
WMuNuSelector::mtMin_
double mtMin_
Definition: WMuNuSelector.cc:60
WMuNuSelector::nid
double nid
Definition: WMuNuSelector.cc:77
looper.cfg
cfg
Definition: looper.py:297
reco::WMuNuCandidate
Definition: WMuNuCandidate.h:17
PVValHelper::dxy
Definition: PVValidationHelpers.h:47
BTaggingMonitoring_cff.njets
njets
Definition: BTaggingMonitoring_cff.py:10
metsig::jet
Definition: SignAlgoResolutions.h:47
reco::WMuNuCandidate::eT
double eT() const
Definition: WMuNuCandidate.cc:25
ev
bool ev
Definition: Hydjet2Hadronizer.cc:95
edm::TriggerNames
Definition: TriggerNames.h:55
WMuNuSelector::h1_
std::map< std::string, TH1D * > h1_
Definition: WMuNuSelector.cc:79
WMuNuSelector::jetToken_
edm::EDGetTokenT< edm::View< reco::Jet > > jetToken_
Definition: WMuNuSelector.cc:47
WMuNuSelector::plotHistograms_
bool plotHistograms_
Definition: WMuNuSelector.cc:44
WMuNuSelector::eJetMin_
double eJetMin_
Definition: WMuNuSelector.cc:53
WMuNuSelector::trigToken_
edm::EDGetTokenT< edm::TriggerResults > trigToken_
Definition: WMuNuSelector.cc:45
WMuNuSelector::nJetMax_
int nJetMax_
Definition: WMuNuSelector.cc:54
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:671
WMuNuSelector::beamSpotToken_
edm::EDGetTokenT< reco::BeamSpot > beamSpotToken_
Definition: WMuNuSelector.cc:48
edm::EDConsumerBase::consumes
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
Definition: EDConsumerBase.h:126
WMuNuSelector::etaCut_
double etaCut_
Definition: WMuNuSelector.cc:56
WMuNuSelector::npresel
double npresel
Definition: WMuNuSelector.cc:74
WMuNuSelector::acopCut_
double acopCut_
Definition: WMuNuSelector.cc:64
reco::LeafCandidate::px
double px() const final
x coordinate of momentum vector
Definition: LeafCandidate.h:140
WMuNuSelector::ptThrForZ1_
double ptThrForZ1_
Definition: WMuNuSelector.cc:51
edm::InputTag
Definition: InputTag.h:15
WMuNuSelector::normalizedChi2Cut_
double normalizedChi2Cut_
Definition: WMuNuSelector.cc:67
TFileService::make
T * make(const Args &... args) const
make new ROOT object
Definition: TFileService.h:64
WMuNuSelector::isoCut03_
double isoCut03_
Definition: WMuNuSelector.cc:59
WMuNuSelector::niso
double niso
Definition: WMuNuSelector.cc:77
WMuNuSelector::ntrig
double ntrig
Definition: WMuNuSelector.cc:74
Generated for CMSSW Reference Manual by doxygen 1.8.16