CMS 3D CMS Logo

List of all members | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes
FilterOutLowPt Class Reference
Inheritance diagram for FilterOutLowPt:
edm::stream::EDFilter<> edm::stream::EDFilterBase edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

 FilterOutLowPt (const edm::ParameterSet &)
 
 ~FilterOutLowPt () override
 
- Public Member Functions inherited from edm::stream::EDFilter<>
 EDFilter ()=default
 
bool hasAbilityToProduceInLumis () const final
 
bool hasAbilityToProduceInRuns () const final
 
- Public Member Functions inherited from edm::stream::EDFilterBase
 EDFilterBase ()
 
ModuleDescription const & moduleDescription () const
 
 ~EDFilterBase () 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
 
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)
 
 ~ProducerBase () noexcept(false) override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfoconsumesInfo () 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
 
EDConsumerBaseoperator= (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)
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &)
 
- Static Public Member Functions inherited from edm::stream::EDFilterBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 

Private Member Functions

virtual void beginJob ()
 
virtual void endJob ()
 
bool filter (edm::Event &, const edm::EventSetup &) override
 

Private Attributes

reco::TrackBase::TrackQuality _trackQuality
 
const bool applyfilter
 
const bool debugOn
 
std::map< unsigned int, std::pair< int, int > > eventsInRun_
 
const unsigned int numtrack
 
double passes
 
const double ptmin
 
const bool runControl_
 
const std::vector< unsigned int > runControlNumbers_
 
edm::EDGetTokenT< reco::TrackCollectiontheTrackCollectionToken
 
const double thresh
 
double trials
 

Additional Inherited Members

- Public Types inherited from edm::stream::EDFilter<>
typedef CacheContexts< T... > CacheTypes
 
typedef CacheTypes::GlobalCache GlobalCache
 
typedef AbilityChecker< T... > HasAbility
 
typedef CacheTypes::LuminosityBlockCache LuminosityBlockCache
 
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCacheLuminosityBlockContext
 
typedef CacheTypes::LuminosityBlockSummaryCache LuminosityBlockSummaryCache
 
typedef CacheTypes::RunCache RunCache
 
typedef RunContextT< RunCache, GlobalCacheRunContext
 
typedef CacheTypes::RunSummaryCache RunSummaryCache
 
- Public Types inherited from edm::stream::EDFilterBase
typedef EDFilterAdaptorBase 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
 
- 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 30 of file FilterOutLowPt.cc.

Constructor & Destructor Documentation

FilterOutLowPt::FilterOutLowPt ( const edm::ParameterSet iConfig)
explicit

Definition at line 61 of file FilterOutLowPt.cc.

61  :
62  applyfilter(iConfig.getUntrackedParameter<bool>("applyfilter")),
63  debugOn(iConfig.getUntrackedParameter<bool>("debugOn")),
64  runControl_(iConfig.getUntrackedParameter<bool>("runControl")),
65  ptmin(iConfig.getUntrackedParameter<double>("ptmin")),
66  thresh(iConfig.getUntrackedParameter<int>("thresh")),
67  numtrack(iConfig.getUntrackedParameter<unsigned int>("numtrack")),
68  runControlNumbers_(iConfig.getUntrackedParameter<std::vector<unsigned int> >("runControlNumber")),
69  theTrackCollectionToken(consumes<reco::TrackCollection>(iConfig.getUntrackedParameter<edm::InputTag>("src")))
70 {
71 }
edm::EDGetTokenT< reco::TrackCollection > theTrackCollectionToken
T getUntrackedParameter(std::string const &, T const &) const
const double thresh
const unsigned int numtrack
const double ptmin
const bool debugOn
const bool applyfilter
const bool runControl_
const std::vector< unsigned int > runControlNumbers_
FilterOutLowPt::~FilterOutLowPt ( )
override

Definition at line 73 of file FilterOutLowPt.cc.

74 {
75 }

Member Function Documentation

void FilterOutLowPt::beginJob ( )
privatevirtual

Definition at line 77 of file FilterOutLowPt.cc.

References passes, and trials.

77  {
78  trials=0;
79  passes=0;
80 }
void FilterOutLowPt::endJob ( )
privatevirtual

Definition at line 160 of file FilterOutLowPt.cc.

References eventsInRun_, plotBeamSpotDB::first, passes, mathSSE::sqrt(), and trials.

Referenced by o2olib.O2ORunMgr::executeJob().

160  {
161 
162  double eff = passes/trials;
163  double eff_err = std::sqrt((eff*(1-eff))/trials);
164 
165  edm::LogVerbatim("FilterOutLowPt")
166  <<"###################################### \n"
167  <<"# FilterOutLowPt::endJob() report \n"
168  <<"# Number of analyzed events: "<<trials<<" \n"
169  <<"# Number of accpeted events: "<<passes<<" \n"
170  <<"# Efficiency: "<< eff*100 << " +/- " << eff_err*100 << " %\n"
171  <<"######################################";
172 
173  edm::LogVerbatim("FilterOutLowPt")<<"###################################### \n"
174  <<"# Filter Summary events accepted by run";
175  for (auto & it : eventsInRun_)
176  edm::LogVerbatim("FilterOutLowPt")<<"# run:" << it.first << " => events tested: " << (it.second).first << " | events passed: " << (it.second).second;
177  edm::LogVerbatim("FilterOutLowPt")<<"###################################### \n";
178 }
std::map< unsigned int, std::pair< int, int > > eventsInRun_
T sqrt(T t)
Definition: SSEVec.h:18
void FilterOutLowPt::fillDescriptions ( edm::ConfigurationDescriptions descriptions)
static

Definition at line 182 of file FilterOutLowPt.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::addUntracked(), DEFINE_FWK_MODULE, and edm::ParameterSetDescription::setComment().

182  {
183 
184  std::vector<unsigned int> defaultRuns;
185  defaultRuns.push_back(0);
186 
188  desc.setComment("Filters out soft events without at least `numtrack` tracks with pT> `ptmin`");
189  desc.addUntracked<bool>("applyfilter",true);
190  desc.addUntracked<bool>("debugOn",false);
191  desc.addUntracked<bool>("runControl",false);
192  desc.addUntracked<unsigned int>("numtrack",0);
193  desc.addUntracked<double>("ptmin",3.);
194  desc.addUntracked<int>("thresh",1);
195  desc.addUntracked<edm::InputTag>("src",edm::InputTag("generalTracks"));
196  desc.addUntracked<std::vector<unsigned int> >("runControlNumber",defaultRuns);
197  descriptions.add("filterOutLowPt", desc);
198 }
ParameterDescriptionBase * addUntracked(U const &iLabel, T const &value)
void setComment(std::string const &value)
void add(std::string const &label, ParameterSetDescription const &psetDescription)
bool FilterOutLowPt::filter ( edm::Event iEvent,
const edm::EventSetup iSetup 
)
overrideprivate

Definition at line 82 of file FilterOutLowPt.cc.

References _trackQuality, cms::dd::accepted(), applyfilter, edm::EventBase::bunchCrossing(), debugOn, edm::EventID::event(), edm::Event::eventAuxiliary(), eventsInRun_, dedxEstimators_cff::fraction, edm::Event::getByToken(), edm::EventBase::id(), edm::EventBase::luminosityBlock(), numtrack, passes, edm::Handle< T >::product(), ptmin, reco::TrackBase::qualityByName(), edm::EventID::run(), edm::EventAuxiliary::run(), runControl_, runControlNumbers_, theTrackCollectionToken, thresh, and trials.

Referenced by Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::filter(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::find(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::setDataAccessor(), and Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::switchCenterView().

83 {
84 
85  bool passesRunControl = false;
86 
87  if(runControl_){
88  if (debugOn){
89  for(unsigned int runControlNumber : runControlNumbers_){
90  edm::LogInfo("FilterOutLowPt")<<"run number:" <<iEvent.id().run()<<" keeping runs:"<<runControlNumber<<std::endl;
91  }
92  }
93 
94  for(unsigned int runControlNumber : runControlNumbers_){
95  if(iEvent.eventAuxiliary().run() == runControlNumber){
96  if (debugOn){
97  edm::LogInfo("FilterOutLowPt")<<"run number"<< runControlNumber << " match!"<<std::endl;
98  }
99  passesRunControl = true;
100  break;
101  }
102  }
103  if (!passesRunControl) return false;
104  }
105 
106  trials++;
107 
108  bool accepted = false;
109  float fraction = 0;
110  // get GeneralTracks collection
111 
113  iEvent.getByToken(theTrackCollectionToken,tkRef);
114  const reco::TrackCollection* tkColl = tkRef.product();
115 
116  int numhighpurity=0;
118 
119  if(tkColl->size()>numtrack){
120  reco::TrackCollection::const_iterator itk = tkColl->begin();
121  reco::TrackCollection::const_iterator itk_e = tkColl->end();
122  for(;itk!=itk_e;++itk){
123  if( itk->quality(_trackQuality) &&
124  (itk->pt() >= ptmin)
125  ) numhighpurity++;
126  }
127  fraction = numhighpurity; //(float)tkColl->size();
128  if(fraction>=thresh) accepted=true;
129  }
130 
131  if (debugOn) {
132  int ievt = iEvent.id().event();
133  int irun = iEvent.id().run();
134  int ils = iEvent.luminosityBlock();
135  int bx = iEvent.bunchCrossing();
136 
137  edm::LogInfo("FilterOutLowPt")<<" Run " << irun << " Event " << ievt << " Lumi Block " << ils << " Bunch Crossing " << bx << " Fraction " << fraction << " NTracks " << tkColl->size() << " Accepted " << accepted << std::endl;
138 
139  }
140 
141  // count the trials and passes
142  unsigned int iRun = iEvent.id().run();
143  if (eventsInRun_.count(iRun)>0){
144  eventsInRun_[iRun].first+=1;
145  if(accepted) eventsInRun_[iRun].second+=1;
146  } else {
147  std::pair<int,int> mypass = std::make_pair(1,0);
148  if(accepted) mypass.second = 1;
149  eventsInRun_[iRun]= mypass;
150  }
151 
152  if (applyfilter){
153  if(accepted) passes++;
154  return accepted;
155  } else
156  return true;
157 
158 }
RunNumber_t run() const
Definition: EventID.h:39
edm::EDGetTokenT< reco::TrackCollection > theTrackCollectionToken
EventNumber_t event() const
Definition: EventID.h:41
EventAuxiliary const & eventAuxiliary() const override
Definition: Event.h:92
const double thresh
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:517
const unsigned int numtrack
std::vector< Track > TrackCollection
collection of Tracks
Definition: TrackFwd.h:15
int bunchCrossing() const
Definition: EventBase.h:64
RunNumber_t run() const
edm::LuminosityBlockNumber_t luminosityBlock() const
Definition: EventBase.h:61
std::map< unsigned int, std::pair< int, int > > eventsInRun_
const double ptmin
reco::TrackBase::TrackQuality _trackQuality
const bool debugOn
static TrackQuality qualityByName(const std::string &name)
Definition: TrackBase.cc:134
T const * product() const
Definition: Handle.h:74
bool accepted(std::vector< std::string_view > const &, std::string_view)
const bool applyfilter
edm::EventID id() const
Definition: EventBase.h:59
const bool runControl_
const std::vector< unsigned int > runControlNumbers_

Member Data Documentation

reco::TrackBase::TrackQuality FilterOutLowPt::_trackQuality
private

Definition at line 56 of file FilterOutLowPt.cc.

Referenced by filter().

const bool FilterOutLowPt::applyfilter
private

Definition at line 43 of file FilterOutLowPt.cc.

Referenced by filter().

const bool FilterOutLowPt::debugOn
private

Definition at line 44 of file FilterOutLowPt.cc.

Referenced by filter().

std::map<unsigned int,std::pair<int,int> > FilterOutLowPt::eventsInRun_
private

Definition at line 54 of file FilterOutLowPt.cc.

Referenced by endJob(), and filter().

const unsigned int FilterOutLowPt::numtrack
private

Definition at line 48 of file FilterOutLowPt.cc.

Referenced by filter().

double FilterOutLowPt::passes
private

Definition at line 51 of file FilterOutLowPt.cc.

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

const double FilterOutLowPt::ptmin
private

Definition at line 46 of file FilterOutLowPt.cc.

Referenced by filter().

const bool FilterOutLowPt::runControl_
private

Definition at line 45 of file FilterOutLowPt.cc.

Referenced by filter().

const std::vector<unsigned int> FilterOutLowPt::runControlNumbers_
private

Definition at line 53 of file FilterOutLowPt.cc.

Referenced by filter().

edm::EDGetTokenT<reco::TrackCollection> FilterOutLowPt::theTrackCollectionToken
private

Definition at line 57 of file FilterOutLowPt.cc.

Referenced by filter().

const double FilterOutLowPt::thresh
private

Definition at line 47 of file FilterOutLowPt.cc.

Referenced by filter().

double FilterOutLowPt::trials
private

Definition at line 50 of file FilterOutLowPt.cc.

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