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, const TrigCodes &trigCodes)
 
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 TrigCodes &trigCodes, 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 103 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().

108 {
109  std::vector<TrigCodes::TrigBitSet> partTrigBits(1);
110  const double maxDeltaR=0.1;
111  for(size_t filterNrInVec=0;filterNrInVec<filters.size();filterNrInVec++){
112  size_t filterNrInEvt = trigEvt->filterIndex(edm::InputTag(filters[filterNrInVec],"",hltTag).encode());
113  //const TrigCodes::TrigBitSet filterCode = trigCodes.getCode(filters[filterNrInVec].c_str());
114  if(filterNrInEvt<trigEvt->sizeFilters()){ //filter found in event, something passes it
115  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
116  const trigger::TriggerObjectCollection & trigObjColl(trigEvt->getObjects());
117  for(trigger::Keys::const_iterator keyIt=trigKeys.begin();keyIt!=trigKeys.end();++keyIt){
118  float trigObjEta = trigObjColl[*keyIt].eta();
119  float trigObjPhi = trigObjColl[*keyIt].phi();
120  float trigObjE = trigObjColl[*keyIt].energy();
121  if (reco::deltaR(particle.superCluster()->eta(),particle.superCluster()->phi(),trigObjEta,trigObjPhi) < maxDeltaR){
122  hltData.HLTeta=trigObjEta;
123  hltData.HLTphi=trigObjPhi;
124  hltData.HLTenergy=trigObjE;
125  }//end dR<maxDeltaR trig obj match test
126  }//end loop over all objects passing filter`
127  }//end check if filter is present
128  }//end check if path has fired in the event
129 }//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
reco::deltaR
double deltaR(const T1 &t1, const T2 &t2)
Definition: deltaR.h:48
filters
std::vector< TPRegexp > filters
Definition: eve_filter.cc:25
trigger::TriggerEvent::getObjects
const TriggerObjectCollection & getObjects() const
Definition: TriggerEvent.h:98
trigger::TriggerObjectCollection
std::vector< TriggerObject > TriggerObjectCollection
collection of trigger physics objects (e.g., all isolated muons)
Definition: TriggerObject.h:81
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 206 of file EgHLTTrigTools.cc.

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

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

207 {
208  //tempory vector to store the filtered results
209  std::vector<std::string> filteredNames;
210 
211  for(size_t inputFilterNr=0;inputFilterNr<namesToFilter.size();inputFilterNr++){
212  std::vector<std::string> names;
213  boost::split(names,namesToFilter[inputFilterNr],boost::is_any_of(std::string(":")));
214  if(names.size()!=2) continue; //format incorrect, reject it
215  if(std::binary_search(activeFilters.begin(),activeFilters.end(),names[0]) &&
216  std::binary_search(activeFilters.begin(),activeFilters.end(),names[1])){ //both filters are valid
217  filteredNames.push_back(namesToFilter[inputFilterNr]);
218  }
219  }
220 
221  namesToFilter.swap(filteredNames);
222 }
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 189 of file EgHLTTrigTools.cc.

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

190 {
191  //tempory vector to store the filtered results
192  std::vector<std::string> filteredNames;
193  /*
194  for(size_t inputFilterNr=0;inputFilterNr<namesToFilter.size();inputFilterNr++){
195  if(std::binary_search(activeFilters.begin(),activeFilters.end(),namesToFilter[inputFilterNr])){
196  filteredNames.push_back(namesToFilter[inputFilterNr]);
197  }//std::cout<<filteredNames[inputFilterNr]<<std::endl;
198  }
199  */
200  //namesToFilter.swap(activeFilters);
201  filteredNames = activeFilters;
202  namesToFilter.swap(filteredNames);
203 }
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 105 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().

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

Definition at line 287 of file EgHLTTrigTools.cc.

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

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

Definition at line 10 of file EgHLTTrigTools.cc.

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

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

15 {
16  TrigCodes::TrigBitSet evtTrigs;
17  for(size_t filterNrInVec=0;filterNrInVec<filters.size();filterNrInVec++){
18  size_t filterNrInEvt = trigEvt->filterIndex(edm::InputTag(filters[filterNrInVec].first,"",hltTag).encode());
19  const TrigCodes::TrigBitSet filterCode = trigCodes.getCode(filters[filterNrInVec].first.c_str());
20  if(filterNrInEvt<trigEvt->sizeFilters()){ //filter found in event, however this only means that something passed the previous filter
21  const trigger::Keys& trigKeys = trigEvt->filterKeys(filterNrInEvt);
22  if(static_cast<int>(trigKeys.size())>=filters[filterNrInVec].second){
23  evtTrigs |=filterCode; //filter was passed
24  }
25  }//end check if filter is present
26  }//end loop over all filters
27 
28 
29  return evtTrigs;
30 
31 }
egHLT::TrigCodes::getCode
TrigBitSet getCode(const char *descript) const
Definition: EgHLTTrigCodes.cc:44
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:248
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 273 of file EgHLTTrigTools.cc.

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

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

274 {
275  const std::vector<std::string>& modules = hltConfig.moduleLabels(path);
276 
277  for(size_t moduleNr=0;moduleNr<modules.size();moduleNr++){
278  const std::string& moduleName=modules[moduleNr];
279  if(moduleName.find("hltL1s")==0) return moduleName; //found l1 seed module
280  }
281  std::string dummy;
282  return dummy;
283 
284 }
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
edm::moduleName
std::string moduleName(Provenance const &provenance)
Definition: Provenance.cc:27
cmsHarvester.path
tuple path
else: Piece not in the list, fine.
Definition: cmsHarvester.py:1591
HLTConfigProvider::moduleLabels
const std::vector< std::string > & moduleLabels(unsigned int trigger) const
label(s) of module(s) on a trigger path
Definition: HLTConfigProvider.h:105
int egHLT::trigTools::getMinNrObjsRequiredByFilter ( const std::string &  filterName)

Definition at line 43 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().

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

Definition at line 308 of file EgHLTTrigTools.cc.

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

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

52  {
53  std::vector<TrigCodes::TrigBitSet> partTrigBits(particles.size());
54  const double maxDeltaR=0.1;
55  for(size_t filterNrInVec=0;filterNrInVec<filters.size();filterNrInVec++){
56  size_t filterNrInEvt = trigEvt->filterIndex(edm::InputTag(filters[filterNrInVec],"",hltTag).encode());
57  const TrigCodes::TrigBitSet filterCode = trigCodes.getCode(filters[filterNrInVec].c_str());
58 
59  if(filterNrInEvt<trigEvt->sizeFilters()){ //filter found in event, something passes it
60  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
61  const trigger::TriggerObjectCollection & trigObjColl(trigEvt->getObjects());
62  for(size_t partNr=0;partNr<particles.size();partNr++){
63  for(trigger::Keys::const_iterator keyIt=trigKeys.begin();keyIt!=trigKeys.end();++keyIt){
64  float trigObjEta = trigObjColl[*keyIt].eta();
65  float trigObjPhi = trigObjColl[*keyIt].phi();
66  if (reco::deltaR(particles[partNr].eta(),particles[partNr].phi(),trigObjEta,trigObjPhi) < maxDeltaR){
67  partTrigBits[partNr] |= filterCode;
68  }//end dR<maxDeltaR trig obj match test
69  }//end loop over all objects passing filter
70  }//end loop over particles
71  }//end check if filter is present
72  }//end loop over all filters
73 
74  //okay the first element is the key, the second is the filter that exists in trigger event
75  for(size_t l1FilterNrInVec=0;l1FilterNrInVec<l1PreAndSeedFilters.size();l1FilterNrInVec++){
76  const TrigCodes::TrigBitSet filterCode = trigCodes.getCode(l1PreAndSeedFilters[l1FilterNrInVec].first.c_str());
77  if((filterCode&evtTrigBits)==filterCode){ //check that filter has fired in the event
78 
79  size_t filterNrInEvt = trigEvt->filterIndex(edm::InputTag(l1PreAndSeedFilters[l1FilterNrInVec].second,"",hltTag).encode());
80 
81  if(filterNrInEvt<trigEvt->sizeFilters()){ //filter found in event, something passes it
82  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
83  const trigger::TriggerObjectCollection & trigObjColl(trigEvt->getObjects());
84  for(size_t partNr=0;partNr<particles.size();partNr++){
85  for(trigger::Keys::const_iterator keyIt=trigKeys.begin();keyIt!=trigKeys.end();++keyIt){
86  float trigObjEta = trigObjColl[*keyIt].eta();
87  float trigObjPhi = trigObjColl[*keyIt].phi();
88  if (reco::deltaR(particles[partNr].eta(),particles[partNr].phi(),trigObjEta,trigObjPhi) < maxDeltaR){
89  partTrigBits[partNr] |= filterCode;
90  }//end dR<maxDeltaR trig obj match test
91  }//end loop over all objects passing filter
92  }//end loop over particles
93  }//end check if filter is present
94  }//end check if path has fired in the event
95  }//end loop over all filters
96 
97  for(size_t partNr=0;partNr<particles.size();partNr++) particles[partNr].setTrigBits(partTrigBits[partNr]);
98 
99  }
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
reco::deltaR
double deltaR(const T1 &t1, const T2 &t2)
Definition: deltaR.h:48
eta
T eta() const
Definition: Basic3DVectorLD.h:177
filters
std::vector< TPRegexp > filters
Definition: eve_filter.cc:25
edm::second
U second(std::pair< T, U > const &p)
Definition: ParameterSet.cc:248
trigger::TriggerEvent::getObjects
const TriggerObjectCollection & getObjects() const
Definition: TriggerEvent.h:98
trigger::TriggerObjectCollection
std::vector< TriggerObject > TriggerObjectCollection
collection of trigger physics objects (e.g., all isolated muons)
Definition: TriggerObject.h:81
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 238 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().

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