CMS 3D CMS Logo

All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
PFRecoTauDiscriminationAgainstMuon.cc
Go to the documentation of this file.
1 /*
2  * class PFRecoTauDiscriminationAgainstMuon
3  * created : May 07 2008,
4  * revised : ,
5  * Authors : Sho Maruyama
6  */
7 
9 
12 
14 
15 
16 #include <string>
17 
18 using namespace reco;
19 
21  public:
23  discriminatorOption_ = iConfig.getParameter<std::string>("discriminatorOption");
24  a = iConfig.getParameter<double>("a");
25  b = iConfig.getParameter<double>("b");
26  c = iConfig.getParameter<double>("c");
27  }
28 
30 
31  double discriminate(const PFTauRef& pfTau);
32 
33  private:
34  std::string discriminatorOption_;
35  double a;
36  double b;
37  double c;
38 };
39 
41 {
42  bool decision = true;
43 
44  if((*thePFTauRef).hasMuonReference() ){
45 
46  MuonRef muonref = (*thePFTauRef).leadPFChargedHadrCand()->muonRef();
47  if (discriminatorOption_ == "noSegMatch") {
48  if ( muonref ->numberOfMatches() > 0 ) {
49  decision = false;
50  }
51  }
52  else if (discriminatorOption_ == "twoDCut") {
53  double seg = muon::segmentCompatibility(*muonref);
54  double calo= muonref->caloCompatibility();
55  double border = calo * a + seg * b +c;
56  if ( border > 0 ) {
57  decision = false;
58  }
59  }
60  else if (discriminatorOption_ == "merePresence") decision = false;
61  else if (discriminatorOption_ == "combined") { // testing purpose only
62  unsigned int muType = 0;
63  if(muonref->isGlobalMuon()) muType = 1;
64  else if(muonref->isCaloMuon()) muType = 2;
65  else if(muonref->isTrackerMuon()) muType = 3;
66  double muonEnergyFraction = (*thePFTauRef).pfTauTagInfoRef()->pfjetRef()->chargedMuEnergyFraction();
67  bool eta_veto = false;
68  bool phi_veto = false;
69  if(fabs(muonref->eta()) > 2.3 || (fabs(muonref->eta()) > 1.4 && fabs(muonref->eta()) < 1.6)) eta_veto = true;
70  if(muonref->phi() < 0.1 && muonref->phi() > -0.1) phi_veto = true;
71  if( muType != 1 || muonref ->numberOfMatches() > 0 || eta_veto || phi_veto || muonEnergyFraction > 0.9 ) decision = false; // as place holder
72  }
73  else if (discriminatorOption_ == "noAllArbitrated") { // One used in H->tautau 2010
75  decision = false;
76  }
77  else{
78  throw edm::Exception(edm::errors::UnimplementedFeature) << " Invalid Discriminator Option! Please check cfi file \n";
79  }
80  } // valid muon ref
81 
82  return (decision ? 1. : 0.);
83 }
84 
T getParameter(std::string const &) const
float segmentCompatibility(const reco::Muon &muon, reco::Muon::ArbitrationType arbitrationType=reco::Muon::SegmentAndTrackArbitration)
bool isGoodMuon(const reco::Muon &muon, SelectionType type, reco::Muon::ArbitrationType arbitrationType=reco::Muon::SegmentAndTrackArbitration)
main GoodMuon wrapper call
DEFINE_FWK_MODULE(CosmicTrackingParticleSelector)
double b
Definition: hdecay.h:120
PFRecoTauDiscriminationAgainstMuon(const edm::ParameterSet &iConfig)
double a
Definition: hdecay.h:121