CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
Functions
egHLT::trigTools Namespace Reference

Functions

template<class T , class U >
void fillHLTposition (T &obj, U &hltData, const std::vector< std::string > &filters, const trigger::TriggerEvent *trigEvt, const std::string &hltTag)
 
void filterInactiveTightLooseTriggers (std::vector< std::string > &namesToFilter, const std::vector< std::string > &activeFilters)
 
void filterInactiveTriggers (std::vector< std::string > &namesToFilter, std::vector< std::string > &activeFilters)
 
void getActiveFilters (const HLTConfigProvider &hltConfig, std::vector< std::string > &activeFilters, std::vector< std::string > &activeEleFilters, std::vector< std::string > &activeEle2LegFilters, std::vector< std::string > &activePhoFilters, std::vector< std::string > &activePho2LegFilters)
 
float getEtThresFromName (const std::string &trigName)
 
TrigCodes::TrigBitSet getFiltersPassed (const std::vector< std::pair< std::string, int > > &filters, const trigger::TriggerEvent *trigEvt, const std::string &hltTag)
 
std::string getL1SeedFilterOfPath (const HLTConfigProvider &hltConfig, const std::string &path)
 
int getMinNrObjsRequiredByFilter (const std::string &filterName)
 
float getSecondEtThresFromName (const std::string &trigName)
 
template<class T >
void setFiltersObjPasses (std::vector< T > &objs, const std::vector< std::string > &filters, const std::vector< std::pair< std::string, std::string > > &l1PreAndSeedFilters, const TrigCodes::TrigBitSet &evtTrigBits, const trigger::TriggerEvent *trigEvt, const std::string &hltTag)
 
void translateFiltersToPathNames (const HLTConfigProvider &hltConfig, const std::vector< std::string > &filters, std::vector< std::string > &paths)
 

Function Documentation

template<class T , class U >
void egHLT::trigTools::fillHLTposition ( T &  obj,
U &  hltData,
const std::vector< std::string > &  filters,
const trigger::TriggerEvent *  trigEvt,
const std::string &  hltTag 
)

Definition at line 102 of file EgHLTTrigTools.h.

References reco::deltaR(), alcaDQMUpload::encode(), trigger::TriggerEvent::filterIndex(), trigger::TriggerEvent::filterKeys(), and trigger::TriggerEvent::getObjects().

Referenced by egHLT::OffHelper::fillHLTData(), and egHLT::OffHelper::fillHLTDataPho().

107 {
108  std::vector<TrigCodes::TrigBitSet> partTrigBits(1);
109  const double maxDeltaR=0.1;
110  for(size_t filterNrInVec=0;filterNrInVec<filters.size();filterNrInVec++){
111  size_t filterNrInEvt = trigEvt->filterIndex(edm::InputTag(filters[filterNrInVec],"",hltTag).encode());
112  //const TrigCodes::TrigBitSet filterCode = TrigCodes::getCode(filters[filterNrInVec].c_str());
113  if(filterNrInEvt<trigEvt->sizeFilters()){ //filter found in event, something passes it
114  const trigger::Keys& trigKeys = trigEvt->filterKeys(filterNrInEvt); //trigger::Keys is actually a vector<uint16_t> holding the position of trigger objects in the trigger collection passing the filter
115  const trigger::TriggerObjectCollection & trigObjColl(trigEvt->getObjects());
116  for(trigger::Keys::const_iterator keyIt=trigKeys.begin();keyIt!=trigKeys.end();++keyIt){
117  float trigObjEta = trigObjColl[*keyIt].eta();
118  float trigObjPhi = trigObjColl[*keyIt].phi();
119  float trigObjE = trigObjColl[*keyIt].energy();
120  if (reco::deltaR(particle.superCluster()->eta(),particle.superCluster()->phi(),trigObjEta,trigObjPhi) < maxDeltaR){
121  hltData.HLTeta=trigObjEta;
122  hltData.HLTphi=trigObjPhi;
123  hltData.HLTenergy=trigObjE;
124  }//end dR<maxDeltaR trig obj match test
125  }//end loop over all objects passing filter`
126  }//end check if filter is present
127  }//end check if path has fired in the event
128 }//end loop over all filters
trigger::TriggerEvent::filterKeys
const Keys & filterKeys(trigger::size_type index) const
Definition: TriggerEvent.h:111
trigger::TriggerEvent::filterIndex
trigger::size_type filterIndex(const edm::InputTag &filterTag) const
find index of filter in data-member vector from filter tag
Definition: TriggerEvent.h:123
filters
std::vector< TPRegexp > filters
Definition: eve_filter.cc:25
trigger::TriggerEvent::getObjects
const TriggerObjectCollection & getObjects() const
Definition: TriggerEvent.h:98
reco::deltaR
auto deltaR(const T1 &t1, const T2 &t2) -> decltype(t1.eta())
Definition: deltaR.h:30
trigger::TriggerObjectCollection
std::vector< TriggerObject > TriggerObjectCollection
collection of trigger physics objects (e.g., all isolated muons)
Definition: TriggerObject.h:81
patRefSel_triggerSelection_cff.hltTag
tuple hltTag
Definition: patRefSel_triggerSelection_cff.py:4
trigger::Keys
std::vector< size_type > Keys
Definition: TriggerTypeDefs.h:20
void egHLT::trigTools::filterInactiveTightLooseTriggers ( std::vector< std::string > &  namesToFilter,
const std::vector< std::string > &  activeFilters 
)

Definition at line 202 of file EgHLTTrigTools.cc.

References cscdqm::h::names, split, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by EgHLTOfflineClient::beginRun(), and EgHLTOfflineSource::filterTriggers().

203 {
204  //tempory vector to store the filtered results
205  std::vector<std::string> filteredNames;
206 
207  for(size_t inputFilterNr=0;inputFilterNr<namesToFilter.size();inputFilterNr++){
208  std::vector<std::string> names;
209  boost::split(names,namesToFilter[inputFilterNr],boost::is_any_of(std::string(":")));
210  if(names.size()!=2) continue; //format incorrect, reject it
211  if(std::binary_search(activeFilters.begin(),activeFilters.end(),names[0]) &&
212  std::binary_search(activeFilters.begin(),activeFilters.end(),names[1])){ //both filters are valid
213  filteredNames.push_back(namesToFilter[inputFilterNr]);
214  }
215  }
216 
217  namesToFilter.swap(filteredNames);
218 }
cscdqm::h::names
static const HistoName names[]
Definition: CSCDQM_HistoDef.h:828
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
split
double split
Definition: MVATrainer.cc:139
void egHLT::trigTools::filterInactiveTriggers ( std::vector< std::string > &  namesToFilter,
std::vector< std::string > &  activeFilters 
)

Definition at line 185 of file EgHLTTrigTools.cc.

Referenced by EgHLTOfflineSummaryClient::beginRun(), EgHLTOfflineClient::beginRun(), and EgHLTOfflineSource::filterTriggers().

186 {
187  //tempory vector to store the filtered results
188  std::vector<std::string> filteredNames;
189  /*
190  for(size_t inputFilterNr=0;inputFilterNr<namesToFilter.size();inputFilterNr++){
191  if(std::binary_search(activeFilters.begin(),activeFilters.end(),namesToFilter[inputFilterNr])){
192  filteredNames.push_back(namesToFilter[inputFilterNr]);
193  }//std::cout<<filteredNames[inputFilterNr]<<std::endl;
194  }
195  */
196  //namesToFilter.swap(activeFilters);
197  filteredNames = activeFilters;
198  namesToFilter.swap(filteredNames);
199 }
void egHLT::trigTools::getActiveFilters ( const HLTConfigProvider &  hltConfig,
std::vector< std::string > &  activeFilters,
std::vector< std::string > &  activeEleFilters,
std::vector< std::string > &  activeEle2LegFilters,
std::vector< std::string > &  activePhoFilters,
std::vector< std::string > &  activePho2LegFilters 
)

Definition at line 101 of file EgHLTTrigTools.cc.

References alcazmumu_cfi::filter, filters, spr::find(), getMinNrObjsRequiredByFilter(), EgammaValidation_cff::pathName, HLTConfigProvider::saveTagsModules(), HLTConfigProvider::size(), python.multivaluedict::sort(), AlCaHLTBitMon_QueryRunRegistry::string, and HLTConfigProvider::triggerName().

Referenced by EgHLTOfflineSummaryClient::beginRun(), EgHLTOfflineClient::beginRun(), and EgHLTOfflineSource::filterTriggers().

102 {
103 
104  activeFilters.clear();
105  activeEleFilters.clear();
106  activeEle2LegFilters.clear();
107  activePhoFilters.clear();
108  activePho2LegFilters.clear();
109 
110  for(size_t pathNr=0;pathNr<hltConfig.size();pathNr++){
111  const std::string& pathName = hltConfig.triggerName(pathNr);
112  if(pathName.find("HLT_")==0){ //hlt path as they all start with HLT_XXXX
113  if((pathName.find("Photon")==4 || pathName.find("Ele")==4 || pathName.find("EG")==4 || pathName.find("Activity")==4 || pathName.find("Physics")==4)// e/g paths, pho or ele always come first
114  && (pathName.find("Jet")==pathName.npos && pathName.find("Muon")==pathName.npos
115  && pathName.find("Tau")==pathName.npos && pathName.find("HT")==pathName.npos
116  && pathName.find("MR")==pathName.npos && pathName.find("LEITI")==pathName.npos
117  && pathName.find("Jpsi")==pathName.npos && pathName.find("Ups")==pathName.npos ) ){//veto x-triggers
118  //std::string lastFilter;
119  const std::vector<std::string>& filters = hltConfig.saveTagsModules(pathNr);
120 
121  //std::cout<<"Number of prescale sets: "<<hltConfig.prescaleSize()<<std::endl;
122  //std::cout<<std::endl<<"Path Name: "<<pathName<<" Prescale: "<<hltConfig.prescaleValue(1,pathName)<<std::endl;
123 
124  if(!filters.empty()){//std::cout<<"Path Name: "<<pathName<<std::endl;
125  //if(filters.back()=="hltBoolEnd" && filters.size()>=2){
126  for(size_t filter=0;filter<filters.size();filter++){
127  if(filters[filter].find("Filter")!=filters[filter].npos){//keep only modules that contain the word "Filter"
128  //std::cout<<" Module Name: "<<filters[filter]<<" filter#: "<<int(filter)<<"/"<<filters.size()<<" ncandcut: "<<trigTools::getMinNrObjsRequiredByFilter(filters[filter])<<std::endl;
129  if(//keep only the last filter and the last one with ncandcut==1 (for di-object triggers)
130  (filter<filters.size()-1 && trigTools::getMinNrObjsRequiredByFilter(filters[filter])==1 && trigTools::getMinNrObjsRequiredByFilter(filters[filter+1])==2)
131  || (filter<filters.size()-1 && trigTools::getMinNrObjsRequiredByFilter(filters[filter])==1 && trigTools::getMinNrObjsRequiredByFilter(filters[filter+1])==1 && filters[filter+1].find("Mass")!=filters[filter+1].npos)
132  || (filter<filters.size()-1 && trigTools::getMinNrObjsRequiredByFilter(filters[filter])==1 && trigTools::getMinNrObjsRequiredByFilter(filters[filter+1])==1 && filters[filter+1].find("FEM")!=filters[filter+1].npos)
133  || (filter<filters.size()-1 && trigTools::getMinNrObjsRequiredByFilter(filters[filter])==1 && trigTools::getMinNrObjsRequiredByFilter(filters[filter+1])==1 && filters[filter+1].find("PFMT")!=filters[filter+1].npos)
134  || filter==filters.size()-1 ){
135  activeFilters.push_back(filters[filter]); //saves all modules with saveTags=true
136  //std::cout<<" Module Name: "<<filters[filter]<<" filter#: "<<int(filter)<<"/"<<filters.size()<<" ncandcut: "<<trigTools::getMinNrObjsRequiredByFilter(filters[filter])<<std::endl;
137  if(pathName.find("Photon")!=pathName.npos || pathName.find("Activity")!=pathName.npos || pathName.find("Physics")!=pathName.npos ){
138  activePhoFilters.push_back(filters[filter]);//saves all "Photon" paths into photon set
139  int posPho = pathName.find("Pho")+1;
140  if( pathName.find("Pho",posPho)!=pathName.npos || pathName.find("SC",posPho)!=pathName.npos ){
141  //This saves all "x_Photon_x_Photon_x" and "x_Photon_x_SC_x" path filters into 2leg photon set
142  activePho2LegFilters.push_back(filters[filter]);
143  //std::cout<<"Pho2LegPath: "<<pathName<<std::endl;
144  }
145  }
146  if(pathName.find("Ele")!=pathName.npos || pathName.find("Activity")!=pathName.npos || pathName.find("Physics")!=pathName.npos ){
147  activeEleFilters.push_back(filters[filter]);//saves all "Ele" paths into electron set
148  int posEle = pathName.find("Ele")+1;
149  if( pathName.find("Ele",posEle)!=pathName.npos || pathName.find("SC",posEle)!=pathName.npos ){
150  if(
151  (filter<filters.size()-1 && trigTools::getMinNrObjsRequiredByFilter(filters[filter])==1 && trigTools::getMinNrObjsRequiredByFilter(filters[filter+1])==2)
152  || (filter<filters.size()-1 && trigTools::getMinNrObjsRequiredByFilter(filters[filter])==1 && trigTools::getMinNrObjsRequiredByFilter(filters[filter+1])==1 && filters[filter+1].find("Mass")!=filters[filter+1].npos)
153  || (filter<filters.size()-1 && trigTools::getMinNrObjsRequiredByFilter(filters[filter])==1 && trigTools::getMinNrObjsRequiredByFilter(filters[filter+1])==1 && filters[filter+1].find("SC")!=filters[filter+1].npos)
154  || (filter<filters.size()-1 && trigTools::getMinNrObjsRequiredByFilter(filters[filter])==1 && trigTools::getMinNrObjsRequiredByFilter(filters[filter+1])==1 && filters[filter+1].find("FEM")!=filters[filter+1].npos)
155  ){
156  //This saves all "x_Ele_x_Ele_x" and "x_Ele_x_SC_x" path filters into 2leg electron set
157  activeEle2LegFilters.push_back(filters[filter]+"::"+filters[filter+1]);
158  //std::cout<<"Ele2LegPath: "<<pathName<<std::endl;
159  }
160  }
161  }
162  }
163  }
164  }
165  //std::cout<<filters[filters.size()-2]<<std::endl;
166  //}else activeFilters.push_back(filters.back());
167  }
168  }
169  }//end hlt path check
170  }//end path loop over
171  /*for(size_t i=0;i<activeEle2LegFilters.size();i++){
172  std::cout<<"Leg1: "<<activeEle2LegFilters[i].substr(0,activeEle2LegFilters[i].find("::"))<<std::endl;
173  std::cout<<"Leg2: "<<activeEle2LegFilters[i].substr(activeEle2LegFilters[i].find("::")+2)<<std::endl<<std::endl;
174  }*/
175  std::sort(activeFilters.begin(),activeFilters.end());
176  std::sort(activeEleFilters.begin(),activeEleFilters.end());
177  std::sort(activeEle2LegFilters.begin(),activeEle2LegFilters.end());
178  std::sort(activePhoFilters.begin(),activePhoFilters.end());
179 }
HLTConfigProvider::size
unsigned int size() const
number of trigger paths in trigger table
Definition: HLTConfigProvider.h:88
HLTConfigProvider::triggerName
const std::string & triggerName(unsigned int triggerIndex) const
Definition: HLTConfigProvider.h:103
egHLT::trigTools::getMinNrObjsRequiredByFilter
int getMinNrObjsRequiredByFilter(const std::string &filterName)
Definition: EgHLTTrigTools.cc:39
alcazmumu_cfi.filter
tuple filter
Definition: alcazmumu_cfi.py:34
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:7
filters
std::vector< TPRegexp > filters
Definition: eve_filter.cc:25
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
HLTConfigProvider::saveTagsModules
const std::vector< std::string > & saveTagsModules(unsigned int trigger) const
Definition: HLTConfigProvider.h:119
EgammaValidation_cff.pathName
tuple pathName
Definition: EgammaValidation_cff.py:269
float egHLT::trigTools::getEtThresFromName ( const std::string &  trigName)

Definition at line 283 of file EgHLTTrigTools.cc.

Referenced by egHLT::MonElemFuncs::addTightLooseTrigHist().

284 {
285  size_t etStrPos = trigName.find("Et");
286  while(etStrPos!=std::string::npos && trigName.find_first_of("1234567890",etStrPos)!=etStrPos+2){
287  etStrPos = trigName.find("Et",etStrPos+1);
288  }
289  if(etStrPos!=std::string::npos && trigName.find_first_of("1234567890",etStrPos)==etStrPos+2){
290  size_t endOfEtValStr = trigName.find_first_not_of("1234567890",etStrPos+2);
291 
292  std::istringstream etValStr(trigName.substr(etStrPos+2,endOfEtValStr-etStrPos-2));
293  float etVal;
294  etValStr>> etVal;
295  return etVal;
296 
297  }
298  return 0;
299 
300 }
TrigCodes::TrigBitSet egHLT::trigTools::getFiltersPassed ( const std::vector< std::pair< std::string, int > > &  filters,
const trigger::TriggerEvent *  trigEvt,
const std::string &  hltTag 
)

Definition at line 10 of file EgHLTTrigTools.cc.

References alcaDQMUpload::encode(), trigger::TriggerEvent::filterIndex(), trigger::TriggerEvent::filterKeys(), filters, first, egHLT::TrigCodes::getCode(), and edm::second().

Referenced by egHLT::OffHelper::setTrigInfo().

11 {
12  TrigCodes::TrigBitSet evtTrigs;
13  for(size_t filterNrInVec=0;filterNrInVec<filters.size();filterNrInVec++){
14  size_t filterNrInEvt = trigEvt->filterIndex(edm::InputTag(filters[filterNrInVec].first,"",hltTag).encode());
15  const TrigCodes::TrigBitSet filterCode = TrigCodes::getCode(filters[filterNrInVec].first.c_str());
16  if(filterNrInEvt<trigEvt->sizeFilters()){ //filter found in event, however this only means that something passed the previous filter
17  const trigger::Keys& trigKeys = trigEvt->filterKeys(filterNrInEvt);
18  if(static_cast<int>(trigKeys.size())>=filters[filterNrInVec].second){
19  evtTrigs |=filterCode; //filter was passed
20  }
21  }//end check if filter is present
22  }//end loop over all filters
23 
24 
25  return evtTrigs;
26 
27 }
trigger::TriggerEvent::filterKeys
const Keys & filterKeys(trigger::size_type index) const
Definition: TriggerEvent.h:111
trigger::TriggerEvent::filterIndex
trigger::size_type filterIndex(const edm::InputTag &filterTag) const
find index of filter in data-member vector from filter tag
Definition: TriggerEvent.h:123
filters
std::vector< TPRegexp > filters
Definition: eve_filter.cc:25
edm::second
U second(std::pair< T, U > const &p)
Definition: ParameterSet.cc:250
first
bool first
Definition: L1TdeRCT.cc:79
patRefSel_triggerSelection_cff.hltTag
tuple hltTag
Definition: patRefSel_triggerSelection_cff.py:4
trigger::Keys
std::vector< size_type > Keys
Definition: TriggerTypeDefs.h:20
egHLT::TrigCodes::TrigBitSet
std::bitset< maxNrBits_ > TrigBitSet
Definition: EgHLTTrigCodes.h:25
std::string egHLT::trigTools::getL1SeedFilterOfPath ( const HLTConfigProvider &  hltConfig,
const std::string &  path 
)

Definition at line 269 of file EgHLTTrigTools.cc.

References HLTConfigProvider::moduleLabels(), edm::moduleName(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by egHLT::OffHelper::setupTriggers().

270 {
271  const std::vector<std::string>& modules = hltConfig.moduleLabels(path);
272 
273  for(size_t moduleNr=0;moduleNr<modules.size();moduleNr++){
274  const std::string& moduleName=modules[moduleNr];
275  if(moduleName.find("hltL1s")==0) return moduleName; //found l1 seed module
276  }
277  std::string dummy;
278  return dummy;
279 
280 }
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
edm::moduleName
std::string moduleName(Provenance const &provenance)
Definition: Provenance.cc:27
HLTConfigProvider::moduleLabels
const std::vector< std::string > & moduleLabels(unsigned int trigger) const
label(s) of module(s) on a trigger path
Definition: HLTConfigProvider.h:113
int egHLT::trigTools::getMinNrObjsRequiredByFilter ( const std::string &  filterName)

Definition at line 39 of file EgHLTTrigTools.cc.

References edm::pset::Registry::begin(), edm::pset::Registry::end(), reco_skim_cfg_mod::filterName, edm::pset::Registry::instance(), and NULL.

Referenced by getActiveFilters(), and egHLT::OffHelper::setupTriggers().

40 {
41 
42  //will return out of for loop once its found it to save time
43  const edm::pset::Registry* psetRegistry = edm::pset::Registry::instance();
44  if(psetRegistry==NULL) return -1;
45  for(edm::pset::Registry::const_iterator psetIt=psetRegistry->begin();psetIt!=psetRegistry->end();++psetIt){ //loop over every pset for every module ever run
46  const std::map<std::string,edm::Entry>& mapOfPara = psetIt->second.tbl(); //contains the parameter name and value for all the parameters of the pset
47  const std::map<std::string,edm::Entry>::const_iterator itToModLabel = mapOfPara.find("@module_label");
48  if(itToModLabel!=mapOfPara.end()){
49  if(itToModLabel->second.getString()==filterName){ //moduleName is the filter name, we have found filter, we will now return something
50  std::map<std::string,edm::Entry>::const_iterator itToCandCut = mapOfPara.find("ncandcut");
51  if(itToCandCut!=mapOfPara.end() && itToCandCut->second.typeCode()=='I') return itToCandCut->second.getInt32();
52  else{ //checks if nZcandcut exists and is int32, if not return -1
53  itToCandCut = mapOfPara.find("nZcandcut");
54  if(itToCandCut!=mapOfPara.end() && itToCandCut->second.typeCode()=='I') return itToCandCut->second.getInt32();
55  else{ //checks if MinN exists and is int32, if not return -1
56  itToCandCut = mapOfPara.find("MinN");
57  if(itToCandCut!=mapOfPara.end() && itToCandCut->second.typeCode()=='I') return itToCandCut->second.getInt32();
58  else{ //checks if minN exists and is int32, if not return -1
59  itToCandCut = mapOfPara.find("minN");
60  if(itToCandCut!=mapOfPara.end() && itToCandCut->second.typeCode()=='I') return itToCandCut->second.getInt32();
61  else return -1;
62  }
63  }
64  }
65  }
66  }
67  }
68  return -1;
69 }
edm::pset::Registry::const_iterator
map_type::const_iterator const_iterator
Definition: Registry.h:63
edm::pset::Registry::end
const_iterator end() const
Definition: Registry.h:69
reco_skim_cfg_mod.filterName
tuple filterName
Definition: reco_skim_cfg_mod.py:160
NULL
#define NULL
Definition: scimark2.h:8
edm::pset::Registry::begin
const_iterator begin() const
Definition: Registry.h:65
edm::pset::Registry::instance
static Registry * instance()
Definition: Registry.cc:14
float egHLT::trigTools::getSecondEtThresFromName ( const std::string &  trigName)

Definition at line 304 of file EgHLTTrigTools.cc.

305 {//std::cout<<"What the heck is this trigName?:"<<trigName<<std::endl;
306  bool isEle=false,isPhoton=false,isEG=false,isEle2=false,isPhoton2=false,isEG2=false,isSC2=false;
307  size_t etStrPos = trigName.npos;
308  if( trigName.find("Ele")<trigName.find("Photon") && trigName.find("Ele")<trigName.find("EG") ){
309  etStrPos=trigName.find("Ele");isEle=true; }
310  else if( trigName.find("EG")<trigName.find("Photon") && trigName.find("EG")<trigName.find("Ele") ){
311  etStrPos=trigName.find("EG");isEG=true; }
312  else if( trigName.find("Photon")<trigName.find("Ele") && trigName.find("Photon")<trigName.find("EG") ){
313  etStrPos=trigName.find("Photon");isPhoton=true; }
314  //size_t etStrPos = trigName.find("Et");
315  //std::cout<<"Got Original Et spot; etStrPos="<<etStrPos<<std::endl;
316  /*while(etStrPos!=std::string::npos && trigName.find_first_of("1234567890",etStrPos)!=etStrPos+2){
317  etStrPos = trigName.find("Et",etStrPos+1);//std::cout<<"Got first Et spot; etStrPos="<<etStrPos<<std::endl;
318  }*/
319  if(etStrPos!=trigName.npos
320  && ( trigName.find("Ele",etStrPos+1)!=trigName.npos || trigName.find("EG",etStrPos+1)!=trigName.npos
321  || trigName.find("Photon",etStrPos+1)!=trigName.npos || trigName.find("SC",etStrPos+1)!=trigName.npos)){
322  if(trigName.find("Ele",etStrPos+1)<trigName.find("Photon",etStrPos+1) && trigName.find("Ele",etStrPos+1)<trigName.find("EG",etStrPos+1) && trigName.find("Ele",etStrPos+1)<trigName.find("SC",etStrPos+1) ){etStrPos=trigName.find("Ele",etStrPos+1);isEle2=true;}
323  else if(trigName.find("EG",etStrPos+1)<trigName.find("Photon",etStrPos+1) && trigName.find("EG",etStrPos+1)<trigName.find("Ele",etStrPos+1) && trigName.find("EG",etStrPos+1)<trigName.find("SC",etStrPos+1)){etStrPos=trigName.find("EG",etStrPos+1);isEG2=true;}
324  else if(trigName.find("Photon",etStrPos+1)<trigName.find("EG",etStrPos+1) && trigName.find("Photon",etStrPos+1)<trigName.find("Ele",etStrPos+1) && trigName.find("Photon",etStrPos+1)<trigName.find("SC",etStrPos+1)){etStrPos=trigName.find("Photon",etStrPos+1);isPhoton2=true;}
325  else if(trigName.find("SC",etStrPos+1)<trigName.find("EG",etStrPos+1) && trigName.find("SC",etStrPos+1)<trigName.find("Ele",etStrPos+1) && trigName.find("SC",etStrPos+1)<trigName.find("Photon",etStrPos+1)){etStrPos=trigName.find("Photon",etStrPos+1);isSC2=true;}
326  //std::cout<<"Got second Et spot; etStrPos="<<etStrPos<<std::endl;//}//get second instance. if it dne, keep first
327 
328  if(isEle2){
329  if(etStrPos!=trigName.npos && trigName.find_first_of("1234567890",etStrPos)==etStrPos+3){//std::cout<<"In if"<<std::endl;
330  size_t endOfEtValStr = trigName.find_first_not_of("1234567890",etStrPos+3);
331 
332  std::istringstream etValStr(trigName.substr(etStrPos+3,endOfEtValStr-etStrPos-3));
333  float etVal;
334  etValStr>> etVal;//std::cout<<"TrigName= "<<trigName<<" etVal= "<<etVal<<std::endl;
335  return etVal;
336  }
337  }
338  if(isEG2 || isSC2){
339  if(etStrPos!=trigName.npos && trigName.find_first_of("1234567890",etStrPos)==etStrPos+2){//std::cout<<"In if"<<std::endl;
340  size_t endOfEtValStr = trigName.find_first_not_of("1234567890",etStrPos+2);
341 
342  std::istringstream etValStr(trigName.substr(etStrPos+2,endOfEtValStr-etStrPos-2));
343  float etVal;
344  etValStr>> etVal;//std::cout<<"TrigName= "<<trigName<<" etVal= "<<etVal<<std::endl;
345  return etVal;
346  }
347  }
348 
349  if(isPhoton2){
350  if(etStrPos!=trigName.npos && trigName.find_first_of("1234567890",etStrPos)==etStrPos+6){//std::cout<<"In if"<<std::endl;
351  size_t endOfEtValStr = trigName.find_first_not_of("1234567890",etStrPos+6);
352 
353  std::istringstream etValStr(trigName.substr(etStrPos+6,endOfEtValStr-etStrPos-6));
354  float etVal;
355  etValStr>> etVal;//std::cout<<"TrigName= "<<trigName<<" etVal= "<<etVal<<std::endl;
356  return etVal;
357  }
358  }
359  }
360  else if(etStrPos!=trigName.npos){
361  if(isEle){
362  if(etStrPos!=trigName.npos && trigName.find_first_of("1234567890",etStrPos)==etStrPos+3){//std::cout<<"In if"<<std::endl;
363  size_t endOfEtValStr = trigName.find_first_not_of("1234567890",etStrPos+3);
364 
365  std::istringstream etValStr(trigName.substr(etStrPos+3,endOfEtValStr-etStrPos-3));
366  float etVal;
367  etValStr>> etVal;//std::cout<<"TrigName= "<<trigName<<" etVal= "<<etVal<<std::endl;
368  return etVal;
369  }
370  }
371  if(isEG){
372  if(etStrPos!=trigName.npos && trigName.find_first_of("1234567890",etStrPos)==etStrPos+2){//std::cout<<"In if"<<std::endl;
373  size_t endOfEtValStr = trigName.find_first_not_of("1234567890",etStrPos+2);
374 
375  std::istringstream etValStr(trigName.substr(etStrPos+2,endOfEtValStr-etStrPos-2));
376  float etVal;
377  etValStr>> etVal;//std::cout<<"TrigName= "<<trigName<<" etVal= "<<etVal<<std::endl;
378  return etVal;
379  }
380  }
381 
382  if(isPhoton){
383  if(etStrPos!=trigName.npos && trigName.find_first_of("1234567890",etStrPos)==etStrPos+6){//std::cout<<"In if"<<std::endl;
384  size_t endOfEtValStr = trigName.find_first_not_of("1234567890",etStrPos+6);
385 
386  std::istringstream etValStr(trigName.substr(etStrPos+6,endOfEtValStr-etStrPos-6));
387  float etVal;
388  etValStr>> etVal;//std::cout<<"TrigName= "<<trigName<<" etVal= "<<etVal<<std::endl;
389  return etVal;
390  }
391  }
392  }
393 
394  return 0;
395 
396 }
template<class T >
void egHLT::trigTools::setFiltersObjPasses ( std::vector< T > &  objs,
const std::vector< std::string > &  filters,
const std::vector< std::pair< std::string, std::string > > &  l1PreAndSeedFilters,
const TrigCodes::TrigBitSet &  evtTrigBits,
const trigger::TriggerEvent *  trigEvt,
const std::string &  hltTag 
)

Definition at line 47 of file EgHLTTrigTools.h.

References reco::deltaR(), alcaDQMUpload::encode(), eta(), trigger::TriggerEvent::filterIndex(), trigger::TriggerEvent::filterKeys(), first, egHLT::TrigCodes::getCode(), trigger::TriggerEvent::getObjects(), phi, and edm::second().

Referenced by egHLT::OffHelper::setTrigInfo().

51  {
52  std::vector<TrigCodes::TrigBitSet> partTrigBits(particles.size());
53  const double maxDeltaR=0.1;
54  for(size_t filterNrInVec=0;filterNrInVec<filters.size();filterNrInVec++){
55  size_t filterNrInEvt = trigEvt->filterIndex(edm::InputTag(filters[filterNrInVec],"",hltTag).encode());
56  const TrigCodes::TrigBitSet filterCode = TrigCodes::getCode(filters[filterNrInVec].c_str());
57 
58  if(filterNrInEvt<trigEvt->sizeFilters()){ //filter found in event, something passes it
59  const trigger::Keys& trigKeys = trigEvt->filterKeys(filterNrInEvt); //trigger::Keys is actually a vector<uint16_t> holding the position of trigger objects in the trigger collection passing the filter
60  const trigger::TriggerObjectCollection & trigObjColl(trigEvt->getObjects());
61  for(size_t partNr=0;partNr<particles.size();partNr++){
62  for(trigger::Keys::const_iterator keyIt=trigKeys.begin();keyIt!=trigKeys.end();++keyIt){
63  float trigObjEta = trigObjColl[*keyIt].eta();
64  float trigObjPhi = trigObjColl[*keyIt].phi();
65  if (reco::deltaR(particles[partNr].eta(),particles[partNr].phi(),trigObjEta,trigObjPhi) < maxDeltaR){
66  partTrigBits[partNr] |= filterCode;
67  }//end dR<maxDeltaR trig obj match test
68  }//end loop over all objects passing filter
69  }//end loop over particles
70  }//end check if filter is present
71  }//end loop over all filters
72 
73  //okay the first element is the key, the second is the filter that exists in trigger event
74  for(size_t l1FilterNrInVec=0;l1FilterNrInVec<l1PreAndSeedFilters.size();l1FilterNrInVec++){
75  const TrigCodes::TrigBitSet filterCode = TrigCodes::getCode(l1PreAndSeedFilters[l1FilterNrInVec].first.c_str());
76  if((filterCode&evtTrigBits)==filterCode){ //check that filter has fired in the event
77 
78  size_t filterNrInEvt = trigEvt->filterIndex(edm::InputTag(l1PreAndSeedFilters[l1FilterNrInVec].second,"",hltTag).encode());
79 
80  if(filterNrInEvt<trigEvt->sizeFilters()){ //filter found in event, something passes it
81  const trigger::Keys& trigKeys = trigEvt->filterKeys(filterNrInEvt); //trigger::Keys is actually a vector<uint16_t> holding the position of trigger objects in the trigger collection passing the filter
82  const trigger::TriggerObjectCollection & trigObjColl(trigEvt->getObjects());
83  for(size_t partNr=0;partNr<particles.size();partNr++){
84  for(trigger::Keys::const_iterator keyIt=trigKeys.begin();keyIt!=trigKeys.end();++keyIt){
85  float trigObjEta = trigObjColl[*keyIt].eta();
86  float trigObjPhi = trigObjColl[*keyIt].phi();
87  if (reco::deltaR(particles[partNr].eta(),particles[partNr].phi(),trigObjEta,trigObjPhi) < maxDeltaR){
88  partTrigBits[partNr] |= filterCode;
89  }//end dR<maxDeltaR trig obj match test
90  }//end loop over all objects passing filter
91  }//end loop over particles
92  }//end check if filter is present
93  }//end check if path has fired in the event
94  }//end loop over all filters
95 
96  for(size_t partNr=0;partNr<particles.size();partNr++) particles[partNr].setTrigBits(partTrigBits[partNr]);
97 
98  }
trigger::TriggerEvent::filterKeys
const Keys & filterKeys(trigger::size_type index) const
Definition: TriggerEvent.h:111
trigger::TriggerEvent::filterIndex
trigger::size_type filterIndex(const edm::InputTag &filterTag) const
find index of filter in data-member vector from filter tag
Definition: TriggerEvent.h:123
eta
T eta() const
Definition: Basic3DVectorLD.h:179
filters
std::vector< TPRegexp > filters
Definition: eve_filter.cc:25
edm::second
U second(std::pair< T, U > const &p)
Definition: ParameterSet.cc:250
trigger::TriggerEvent::getObjects
const TriggerObjectCollection & getObjects() const
Definition: TriggerEvent.h:98
reco::deltaR
auto deltaR(const T1 &t1, const T2 &t2) -> decltype(t1.eta())
Definition: deltaR.h:30
first
bool first
Definition: L1TdeRCT.cc:79
trigger::TriggerObjectCollection
std::vector< TriggerObject > TriggerObjectCollection
collection of trigger physics objects (e.g., all isolated muons)
Definition: TriggerObject.h:81
patRefSel_triggerSelection_cff.hltTag
tuple hltTag
Definition: patRefSel_triggerSelection_cff.py:4
trigger::Keys
std::vector< size_type > Keys
Definition: TriggerTypeDefs.h:20
phi
Definition: DDAxes.h:10
void egHLT::trigTools::translateFiltersToPathNames ( const HLTConfigProvider &  hltConfig,
const std::vector< std::string > &  filters,
std::vector< std::string > &  paths 
)

Definition at line 234 of file EgHLTTrigTools.cc.

References HLTConfigProvider::moduleLabels(), EgammaValidation_cff::pathName, HLTConfigProvider::size(), python.multivaluedict::sort(), AlCaHLTBitMon_QueryRunRegistry::string, and HLTConfigProvider::triggerName().

Referenced by EgHLTOfflineSummaryClient::beginRun(), and egHLT::OffHelper::setupTriggers().

235 {
236 
237  paths.clear();
238  std::vector<std::pair<std::string,std::string> > filtersAndPaths;
239 
240  for(size_t pathNr=0;pathNr<hltConfig.size();pathNr++){
241  const std::string& pathName = hltConfig.triggerName(pathNr);
242  if(pathName.find("HLT_")==0){ //hlt path as they all start with HLT_XXXX
243 
244  std::string lastFilter;
245  const std::vector<std::string>& pathFilters = hltConfig.moduleLabels(pathNr);
246  if(!pathFilters.empty()){
247  if(pathFilters.back()=="hltBoolEnd" && pathFilters.size()>=2){
248  //2nd to last element is the last filter, useally the case as last is hltBool except for ES bits
249  filtersAndPaths.push_back(std::make_pair(pathFilters[pathFilters.size()-2],pathName));
250  }else filtersAndPaths.push_back(std::make_pair(pathFilters.back(),pathName));
251  }
252  }//end hlt path check
253  }//end path loop over
254 
255  std::sort(filtersAndPaths.begin(),filtersAndPaths.end(),StringPairCompare());
256 
257  for(size_t filterNr=0;filterNr<filters.size();filterNr++){
258  typedef std::vector<std::pair<std::string,std::string> >::const_iterator VecIt;
259  std::pair<VecIt,VecIt> searchResult = std::equal_range(filtersAndPaths.begin(),filtersAndPaths.end(),filters[filterNr],StringPairCompare());
260  if(searchResult.first!=searchResult.second) paths.push_back(searchResult.first->second);
261  else paths.push_back(filters[filterNr]);//if cant find the path, just write the filter
262  //---Morse-----
263  //std::cout<<filtersAndPaths[filterNr].first<<" "<<filtersAndPaths[filterNr].second<<std::endl;
264  //-------------
265  }
266 
267 }
HLTConfigProvider::size
unsigned int size() const
number of trigger paths in trigger table
Definition: HLTConfigProvider.h:88
HLTConfigProvider::triggerName
const std::string & triggerName(unsigned int triggerIndex) const
Definition: HLTConfigProvider.h:103
EgammaValidation_cff.paths
tuple paths
Definition: EgammaValidation_cff.py:18
filters
std::vector< TPRegexp > filters
Definition: eve_filter.cc:25
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
EgammaValidation_cff.pathName
tuple pathName
Definition: EgammaValidation_cff.py:269
HLTConfigProvider::moduleLabels
const std::vector< std::string > & moduleLabels(unsigned int trigger) const
label(s) of module(s) on a trigger path
Definition: HLTConfigProvider.h:113
Generated for CMSSW Reference Manual by doxygen 1.8.5