test
CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Classes | Public Member Functions | Private Member Functions | Private Attributes
pat::PATTriggerProducer Class Reference

Produces the full or stand-alone PAT trigger information collections. More...

#include "PhysicsTools/PatAlgos/plugins/PATTriggerProducer.h"

Inheritance diagram for pat::PATTriggerProducer:
edm::stream::EDProducer<> edm::stream::EDProducerBase edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Classes

class  ModuleLabelToPathAndFlags
 

Public Member Functions

 PATTriggerProducer (const edm::ParameterSet &iConfig)
 
 ~PATTriggerProducer ()
 
- Public Member Functions inherited from edm::stream::EDProducer<>
 EDProducer ()=default
 
- Public Member Functions inherited from edm::stream::EDProducerBase
 EDProducerBase ()
 
ModuleDescription const & moduleDescription () const
 
virtual ~EDProducerBase ()
 
- Public Member Functions inherited from edm::ProducerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 ProducerBase ()
 
void registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
 
std::function< void(BranchDescription
const &)> 
registrationCallback () const
 used by the fwk to register list of products More...
 
virtual ~ProducerBase ()
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfoconsumesInfo () const
 
 EDConsumerBase ()
 
ProductHolderIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
std::vector
< ProductHolderIndexAndSkipBit >
const & 
itemsToGetFromEvent () const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
bool registeredToConsume (ProductHolderIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
 
virtual ~EDConsumerBase ()
 

Private Member Functions

virtual void beginLuminosityBlock (const edm::LuminosityBlock &iLuminosityBlock, const edm::EventSetup &iSetup) override
 
virtual void beginRun (const edm::Run &iRun, const edm::EventSetup &iSetup) override
 
virtual void produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override
 

Private Attributes

bool addL1Algos_
 
bool addPathModuleLabels_
 
bool autoProcessName_
 
bool autoProcessNameL1ExtraCenJet_
 
bool autoProcessNameL1ExtraETM_
 
bool autoProcessNameL1ExtraForJet_
 
bool autoProcessNameL1ExtraHTM_
 
bool autoProcessNameL1ExtraIsoEG_
 
bool autoProcessNameL1ExtraMu_
 
bool autoProcessNameL1ExtraNoIsoEG_
 
bool autoProcessNameL1ExtraTauJet_
 
std::vector< std::string > exludeCollections_
 
bool firstInRun_
 
bool hltConfigInit_
 
std::string hltPrescaleLabel_
 
HLTPrescaleProvider hltPrescaleProvider_
 
edm::GetterOfProducts
< trigger::HLTPrescaleTable
hltPrescaleTableEventGetter_
 
trigger::HLTPrescaleTable hltPrescaleTableLumi_
 
edm::GetterOfProducts
< trigger::HLTPrescaleTable
hltPrescaleTableLumiGetter_
 
trigger::HLTPrescaleTable hltPrescaleTableRun_
 
edm::GetterOfProducts
< trigger::HLTPrescaleTable
hltPrescaleTableRunGetter_
 
edm::GetterOfProducts
< l1extra::L1JetParticleCollection
l1ExtraCenJetGetter_
 
edm::GetterOfProducts
< l1extra::L1EtMissParticleCollection
l1ExtraETMGetter_
 
edm::GetterOfProducts
< l1extra::L1JetParticleCollection
l1ExtraForJetGetter_
 
edm::GetterOfProducts
< l1extra::L1EtMissParticleCollection
l1ExtraHTMGetter_
 
edm::GetterOfProducts
< l1extra::L1EmParticleCollection
l1ExtraIsoEGGetter_
 
edm::GetterOfProducts
< l1extra::L1MuonParticleCollection
l1ExtraMuGetter_
 
edm::GetterOfProducts
< l1extra::L1EmParticleCollection
l1ExtraNoIsoEGGetter_
 
edm::GetterOfProducts
< l1extra::L1JetParticleCollection
l1ExtraTauJetGetter_
 
edm::EDGetTokenT
< L1GlobalTriggerObjectMaps
l1GlobalTriggerObjectMapsToken_
 
edm::ParameterSetl1PSet_
 
std::string labelHltPrescaleTable_
 
bool mainBxOnly_
 
ModuleLabelToPathAndFlags moduleLabelToPathAndFlags_
 
std::string nameProcess_
 
bool onlyStandAlone_
 
bool packPathNames_
 
bool packPrescales_
 
bool saveL1Refs_
 
edm::InputTag tagL1ExtraCenJet_
 
edm::InputTag tagL1ExtraETM_
 
edm::InputTag tagL1ExtraForJet_
 
edm::InputTag tagL1ExtraHTM_
 
edm::InputTag tagL1ExtraIsoEG_
 
edm::InputTag tagL1ExtraMu_
 
edm::InputTag tagL1ExtraNoIsoEG_
 
edm::InputTag tagL1ExtraTauJet_
 
edm::InputTag tagL1GlobalTriggerObjectMaps_
 
edm::InputTag tagTriggerEvent_
 
edm::InputTag tagTriggerResults_
 
edm::GetterOfProducts
< trigger::TriggerEvent
triggerEventGetter_
 
edm::GetterOfProducts
< edm::TriggerResults
triggerResultsGetter_
 

Additional Inherited Members

- Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T...> CacheTypes
 
typedef CacheTypes::GlobalCache GlobalCache
 
typedef AbilityChecker< T...> HasAbility
 
typedef
CacheTypes::LuminosityBlockCache 
LuminosityBlockCache
 
typedef
LuminosityBlockContextT
< LuminosityBlockCache,
RunCache, GlobalCache
LuminosityBlockContext
 
typedef
CacheTypes::LuminosityBlockSummaryCache 
LuminosityBlockSummaryCache
 
typedef CacheTypes::RunCache RunCache
 
typedef RunContextT< RunCache,
GlobalCache
RunContext
 
typedef CacheTypes::RunSummaryCache RunSummaryCache
 
- Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase ModuleType
 
- Public Types inherited from edm::ProducerBase
typedef
ProductRegistryHelper::TypeLabelList 
TypeLabelList
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::stream::EDProducerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
- 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 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

Produces the full or stand-alone PAT trigger information collections.

This producer extracts the trigger information from

For me information, s. https://twiki.cern.ch/twiki/bin/view/CMS/SWGuidePATTrigger

Author
Volker Adler
Version
Id:
PATTriggerProducer.h,v 1.20 2012/09/11 22:45:29 vadler Exp

Definition at line 58 of file PATTriggerProducer.h.

Constructor & Destructor Documentation

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

Definition at line 46 of file PATTriggerProducer.cc.

References addL1Algos_, addPathModuleLabels_, autoProcessName_, autoProcessNameL1ExtraCenJet_, autoProcessNameL1ExtraETM_, autoProcessNameL1ExtraForJet_, autoProcessNameL1ExtraHTM_, autoProcessNameL1ExtraIsoEG_, autoProcessNameL1ExtraMu_, autoProcessNameL1ExtraNoIsoEG_, autoProcessNameL1ExtraTauJet_, edm::ProducerBase::callWhenNewProductsRegistered(), edm::ParameterSet::exists(), exludeCollections_, edm::ParameterSet::getParameter(), hltPrescaleLabel_, hltPrescaleTableEventGetter_, hltPrescaleTableLumiGetter_, hltPrescaleTableRunGetter_, edm::InLumi, HLT_25ns10e33_v2_cff::InputTag, edm::InRun, edm::InputTag::instance(), l1ExtraCenJetGetter_, l1ExtraETMGetter_, l1ExtraForJetGetter_, l1ExtraHTMGetter_, l1ExtraIsoEGGetter_, l1ExtraMuGetter_, l1ExtraNoIsoEGGetter_, l1ExtraTauJetGetter_, l1GlobalTriggerObjectMapsToken_, edm::InputTag::label(), labelHltPrescaleTable_, mainBxOnly_, nameProcess_, onlyStandAlone_, packPrescales_, edm::InputTag::process(), saveL1Refs_, AlCaHLTBitMon_QueryRunRegistry::string, tagL1ExtraCenJet_, tagL1ExtraETM_, tagL1ExtraForJet_, tagL1ExtraHTM_, tagL1ExtraIsoEG_, tagL1ExtraMu_, tagL1ExtraNoIsoEG_, tagL1ExtraTauJet_, tagL1GlobalTriggerObjectMaps_, tagTriggerEvent_, tagTriggerResults_, triggerEventGetter_, and triggerResultsGetter_.

46  :
47  nameProcess_( iConfig.getParameter< std::string >( "processName" ) ),
49  onlyStandAlone_( iConfig.getParameter< bool >( "onlyStandAlone" ) ),
50  firstInRun_( true ),
51  // L1 configuration parameters
52  addL1Algos_( false ),
53  tagL1GlobalTriggerObjectMaps_( "l1L1GtObjectMap" ),
54  tagL1ExtraMu_(),
70  mainBxOnly_( true ),
71  saveL1Refs_( false ),
72  hltPrescaleProvider_(iConfig, consumesCollector(), *this),
73  hltConfigInit_( false ),
74  // HLT configuration parameters
75  tagTriggerResults_( "TriggerResults" ),
76  tagTriggerEvent_( "hltTriggerSummaryAOD" ),
81  addPathModuleLabels_( false ),
82  packPathNames_( iConfig.existsAs<bool>("packTriggerPathNames") ? iConfig.getParameter<bool>("packTriggerPathNames") : false ),
83  packPrescales_( iConfig.existsAs<bool>("packTriggerPrescales") ? iConfig.getParameter<bool>("packTriggerPrescales") : true )
84 {
85 
86  // L1 configuration parameters
87  if ( iConfig.exists( "addL1Algos" ) ) addL1Algos_ = iConfig.getParameter< bool >( "addL1Algos" );
88  if ( iConfig.exists( "l1GlobalTriggerObjectMaps" ) ) tagL1GlobalTriggerObjectMaps_ = iConfig.getParameter< InputTag >( "l1GlobalTriggerObjectMaps" );
89  l1GlobalTriggerObjectMapsToken_ = mayConsume< L1GlobalTriggerObjectMaps >( tagL1GlobalTriggerObjectMaps_ );
90  if ( iConfig.exists( "l1ExtraMu" ) ) {
91  tagL1ExtraMu_ = iConfig.getParameter< InputTag >( "l1ExtraMu" );
92  if ( tagL1ExtraMu_.process() == "*" ) {
95  }
97  }
98  if ( iConfig.exists( "l1ExtraNoIsoEG" ) ) {
99  tagL1ExtraNoIsoEG_ = iConfig.getParameter< InputTag >( "l1ExtraNoIsoEG" );
100  if ( tagL1ExtraNoIsoEG_.process() == "*" ) {
103  }
105  }
106  if ( iConfig.exists( "l1ExtraIsoEG" ) ) {
107  tagL1ExtraIsoEG_ = iConfig.getParameter< InputTag >( "l1ExtraIsoEG" );
108  if ( tagL1ExtraIsoEG_.process() == "*" ) {
111  }
113  }
114  if ( iConfig.exists( "l1ExtraCenJet" ) ) {
115  tagL1ExtraCenJet_ = iConfig.getParameter< InputTag >( "l1ExtraCenJet" );
116  if ( tagL1ExtraCenJet_.process() == "*" ) {
119  }
121  }
122  if ( iConfig.exists( "l1ExtraForJet" ) ) {
123  tagL1ExtraForJet_ = iConfig.getParameter< InputTag >( "l1ExtraForJet" );
124  if ( tagL1ExtraForJet_.process() == "*" ) {
127  }
129  }
130  if ( iConfig.exists( "l1ExtraTauJet" ) ) {
131  tagL1ExtraTauJet_ = iConfig.getParameter< InputTag >( "l1ExtraTauJet" );
132  if ( tagL1ExtraTauJet_.process() == "*" ) {
135  }
137  }
138  if ( iConfig.exists( "l1ExtraETM" ) ) {
139  tagL1ExtraETM_ = iConfig.getParameter< InputTag >( "l1ExtraETM" );
140  if ( tagL1ExtraETM_.process() == "*" ) {
143  }
145  }
146  if ( iConfig.exists( "l1ExtraHTM" ) ) {
147  tagL1ExtraHTM_ = iConfig.getParameter< InputTag >( "l1ExtraHTM" );
148  if ( tagL1ExtraHTM_.process() == "*" ) {
151  }
153  }
154  if ( iConfig.exists( "mainBxOnly" ) ) mainBxOnly_ = iConfig.getParameter< bool >( "mainBxOnly" );
155  if ( iConfig.exists( "saveL1Refs" ) ) saveL1Refs_ = iConfig.getParameter< bool >( "saveL1Refs" );
156 
157  // HLT configuration parameters
158  if ( iConfig.exists( "triggerResults" ) ) tagTriggerResults_ = iConfig.getParameter< InputTag >( "triggerResults" );
160  if ( iConfig.exists( "triggerEvent" ) ) tagTriggerEvent_ = iConfig.getParameter< InputTag >( "triggerEvent" );
162  if ( iConfig.exists( "hltPrescaleLabel" ) ) hltPrescaleLabel_ = iConfig.getParameter< std::string >( "hltPrescaleLabel" );
163  if ( iConfig.exists( "hltPrescaleTable" ) ) {
164  labelHltPrescaleTable_ = iConfig.getParameter< std::string >( "hltPrescaleTable" );
168  }
169  if ( iConfig.exists( "addPathModuleLabels" ) ) addPathModuleLabels_ = iConfig.getParameter< bool >( "addPathModuleLabels" );
170  exludeCollections_.clear();
171  if ( iConfig.exists( "exludeCollections" ) ) exludeCollections_ = iConfig.getParameter< std::vector< std::string > >( "exludeCollections" );
172 
173  callWhenNewProductsRegistered( [ this, &iConfig ]( BranchDescription const& bd ) {
174  if ( iConfig.exists( "l1ExtraMu" ) ) l1ExtraMuGetter_( bd );
175  if ( iConfig.exists( "l1ExtraNoIsoEG" ) ) l1ExtraNoIsoEGGetter_( bd );
176  if ( iConfig.exists( "l1ExtraIsoEG" ) ) l1ExtraIsoEGGetter_( bd );
177  if ( iConfig.exists( "l1ExtraCenJet" ) ) l1ExtraCenJetGetter_( bd );
178  if ( iConfig.exists( "l1ExtraForJet" ) ) l1ExtraForJetGetter_( bd );
179  if ( iConfig.exists( "l1ExtraTauJet" ) ) l1ExtraTauJetGetter_( bd );
180  if ( iConfig.exists( "l1ExtraETM" ) ) l1ExtraETMGetter_( bd );
181  if ( iConfig.exists( "l1ExtraHTM" ) ) l1ExtraHTMGetter_( bd );
182  triggerResultsGetter_( bd );
183  triggerEventGetter_( bd );
184  if ( iConfig.exists( "hltPrescaleTable" ) ) {
188  }
189  } );
190 
191  if ( ! onlyStandAlone_ ) {
192  produces< TriggerAlgorithmCollection >();
193  produces< TriggerConditionCollection >();
194  produces< TriggerPathCollection >();
195  produces< TriggerFilterCollection >();
196  produces< TriggerObjectCollection >();
197  }
198  if (packPrescales_) {
199  produces< PackedTriggerPrescales >();
200  produces< PackedTriggerPrescales >("l1max");
201  produces< PackedTriggerPrescales >("l1min");
202  }
203  produces< TriggerObjectStandAloneCollection >();
204 
205 }
T getParameter(std::string const &) const
void callWhenNewProductsRegistered(std::function< void(BranchDescription const &)> const &func)
Definition: ProducerBase.h:50
edm::GetterOfProducts< l1extra::L1MuonParticleCollection > l1ExtraMuGetter_
edm::GetterOfProducts< trigger::HLTPrescaleTable > hltPrescaleTableEventGetter_
bool existsAs(std::string const &parameterName, bool trackiness=true) const
checks if a parameter exists as a given type
Definition: ParameterSet.h:186
edm::GetterOfProducts< trigger::TriggerEvent > triggerEventGetter_
bool exists(std::string const &parameterName) const
checks if a parameter exists
edm::GetterOfProducts< trigger::HLTPrescaleTable > hltPrescaleTableLumiGetter_
edm::GetterOfProducts< l1extra::L1JetParticleCollection > l1ExtraForJetGetter_
ConsumesCollector consumesCollector()
Use a ConsumesCollector to gather consumes information from helper functions.
edm::GetterOfProducts< l1extra::L1JetParticleCollection > l1ExtraTauJetGetter_
trigger::HLTPrescaleTable hltPrescaleTableLumi_
edm::GetterOfProducts< l1extra::L1EtMissParticleCollection > l1ExtraHTMGetter_
edm::EDGetTokenT< L1GlobalTriggerObjectMaps > l1GlobalTriggerObjectMapsToken_
edm::GetterOfProducts< l1extra::L1EmParticleCollection > l1ExtraNoIsoEGGetter_
edm::GetterOfProducts< l1extra::L1EtMissParticleCollection > l1ExtraETMGetter_
edm::GetterOfProducts< trigger::HLTPrescaleTable > hltPrescaleTableRunGetter_
edm::InputTag tagL1GlobalTriggerObjectMaps_
edm::GetterOfProducts< l1extra::L1EmParticleCollection > l1ExtraIsoEGGetter_
edm::GetterOfProducts< l1extra::L1JetParticleCollection > l1ExtraCenJetGetter_
std::string const & label() const
Definition: InputTag.h:36
std::string const & process() const
Definition: InputTag.h:40
edm::GetterOfProducts< edm::TriggerResults > triggerResultsGetter_
trigger::HLTPrescaleTable hltPrescaleTableRun_
std::vector< std::string > exludeCollections_
std::string const & instance() const
Definition: InputTag.h:37
edm::InputTag tagL1ExtraNoIsoEG_
HLTPrescaleProvider hltPrescaleProvider_
pat::PATTriggerProducer::~PATTriggerProducer ( )
inline

Definition at line 63 of file PATTriggerProducer.h.

63 {};

Member Function Documentation

void PATTriggerProducer::beginLuminosityBlock ( const edm::LuminosityBlock iLuminosityBlock,
const edm::EventSetup iSetup 
)
overrideprivatevirtual

Reimplemented from edm::stream::EDProducerBase.

Definition at line 291 of file PATTriggerProducer.cc.

References edm::LuminosityBlock::getByLabel(), hltConfigInit_, hltPrescaleTableLumi_, hltPrescaleTableRun_, HLT_25ns10e33_v2_cff::InputTag, edm::HandleBase::isValid(), labelHltPrescaleTable_, trigger::HLTPrescaleTable::labels(), nameProcess_, trigger::HLTPrescaleTable::set(), and trigger::HLTPrescaleTable::table().

292 {
293 
294  // Terminate, if auto process name determination failed
295  if ( nameProcess_ == "*" ) return;
296 
297  // Extract pre-scales
298  if ( hltConfigInit_ ) {
299  // Start from run
301  // Try lumi product, if configured and available
302  if ( ! labelHltPrescaleTable_.empty() ) {
303  Handle< trigger::HLTPrescaleTable > handleHltPrescaleTable;
304  iLuminosityBlock.getByLabel( InputTag( labelHltPrescaleTable_, "Lumi", nameProcess_ ), handleHltPrescaleTable );
305  if ( handleHltPrescaleTable.isValid() ) {
306  hltPrescaleTableLumi_ = trigger::HLTPrescaleTable( handleHltPrescaleTable->set(), handleHltPrescaleTable->labels(), handleHltPrescaleTable->table() );
307  }
308  }
309  }
310 
311 }
unsigned int set() const
low-level const accessors for data members
The single EDProduct containing the HLT Prescale Table.
bool getByLabel(std::string const &label, Handle< PROD > &result) const
const std::map< std::string, std::vector< unsigned int > > & table() const
const std::vector< std::string > & labels() const
trigger::HLTPrescaleTable hltPrescaleTableLumi_
bool isValid() const
Definition: HandleBase.h:75
trigger::HLTPrescaleTable hltPrescaleTableRun_
void PATTriggerProducer::beginRun ( const edm::Run iRun,
const edm::EventSetup iSetup 
)
overrideprivatevirtual

Reimplemented from edm::stream::EDProducerBase.

Definition at line 208 of file PATTriggerProducer.cc.

References autoProcessName_, autoProcessNameL1ExtraCenJet_, autoProcessNameL1ExtraETM_, autoProcessNameL1ExtraForJet_, autoProcessNameL1ExtraHTM_, autoProcessNameL1ExtraIsoEG_, autoProcessNameL1ExtraMu_, autoProcessNameL1ExtraNoIsoEG_, autoProcessNameL1ExtraTauJet_, edm::ParameterSet::exists(), firstInRun_, edm::Run::getByLabel(), HltComparatorCreateWorkflow::hltConfig, hltConfigInit_, HLTPrescaleProvider::hltConfigProvider(), hltPrescaleProvider_, hltPrescaleTableRun_, HLTPrescaleProvider::init(), pat::PATTriggerProducer::ModuleLabelToPathAndFlags::init(), HLT_25ns10e33_v2_cff::InputTag, edm::InputTag::instance(), instance, edm::HandleBase::isValid(), l1PSet_, edm::InputTag::label(), labelHltPrescaleTable_, LogDebug, moduleLabelToPathAndFlags_, nameProcess_, edm::InputTag::process(), edm::processConfiguration, edm::Run::processHistory(), tagL1ExtraCenJet_, tagL1ExtraETM_, tagL1ExtraForJet_, tagL1ExtraHTM_, tagL1ExtraIsoEG_, tagL1ExtraMu_, tagL1ExtraNoIsoEG_, tagL1ExtraTauJet_, tagTriggerEvent_, and tagTriggerResults_.

209 {
210 
211  // Initialize
212  firstInRun_ = true;
213  l1PSet_ = 0;
214  hltConfigInit_ = false;
215 
216  // Initialize process name
217  if ( autoProcessName_ ) {
218  // reset
219  nameProcess_ = "*";
220  // determine process name from last run TriggerSummaryProducerAOD module in process history of input
221  const ProcessHistory & processHistory( iRun.processHistory() );
223  ParameterSet processPSet;
224  // unbroken loop, which relies on time ordering (accepts the last found entry)
225  for ( ProcessHistory::const_iterator iHist = processHistory.begin(); iHist != processHistory.end(); ++iHist ) {
226  if ( processHistory.getConfigurationForProcess( iHist->processName(), processConfiguration ) &&
227  pset::Registry::instance()->getMapped( processConfiguration.parameterSetID(), processPSet ) &&
228  processPSet.exists( tagTriggerEvent_.label() )
229  ) {
230  nameProcess_ = iHist->processName();
231  LogDebug( "autoProcessName" ) << "HLT process name '" << nameProcess_ << "' discovered";
232  }
233  }
234  // terminate, if nothing is found
235  if ( nameProcess_ == "*" ) {
236  LogError( "autoProcessName" ) << "trigger::TriggerEvent product with label '" << tagTriggerEvent_.label() << "' not produced according to process history of input data\n"
237  << "No trigger information produced";
238  return;
239  }
240  LogInfo( "autoProcessName" ) << "HLT process name' " << nameProcess_ << "' used for PAT trigger information";
241  }
242  // adapt configuration of used input tags
243  if ( tagTriggerResults_.process().empty() || tagTriggerResults_.process() == "*" ) {
245  } else if ( tagTriggerEvent_.process() != nameProcess_ ) {
246  LogWarning( "inputTags" ) << "TriggerResults process name '" << tagTriggerResults_.process() << "' differs from HLT process name '" << nameProcess_ << "'";
247  }
248  if ( tagTriggerEvent_.process().empty() || tagTriggerEvent_.process() == "*" ) {
250  } else if ( tagTriggerEvent_.process() != nameProcess_ ) {
251  LogWarning( "inputTags" ) << "TriggerEvent process name '" << tagTriggerEvent_.process() << "' differs from HLT process name '" << nameProcess_ << "'";
252  }
261 
262  // Initialize HLTConfigProvider
264  bool changed( true );
265  if ( ! hltPrescaleProvider_.init( iRun, iSetup, nameProcess_, changed ) ) {
266  LogError( "hltConfig" ) << "HLT config extraction error with process name '" << nameProcess_ << "'";
267  } else if ( hltConfig.size() <= 0 ) {
268  LogError( "hltConfig" ) << "HLT config size error";
269  } else hltConfigInit_ = true;
270 
271  // Update mapping from filter names to path names
272  if (hltConfigInit_ && changed) moduleLabelToPathAndFlags_.init( hltConfig );
273 
274  // Extract pre-scales
275  if ( hltConfigInit_ ) {
276  // Start empty
278  // Try run product, if configured
279  if ( ! labelHltPrescaleTable_.empty() ) {
280  Handle< trigger::HLTPrescaleTable > handleHltPrescaleTable;
281  iRun.getByLabel( InputTag( labelHltPrescaleTable_, "Run", nameProcess_ ), handleHltPrescaleTable );
282  if ( handleHltPrescaleTable.isValid() ) {
283  hltPrescaleTableRun_ = trigger::HLTPrescaleTable( handleHltPrescaleTable->set(), handleHltPrescaleTable->labels(), handleHltPrescaleTable->table() );
284  }
285  }
286  }
287 
288 }
#define LogDebug(id)
collection_type::const_iterator const_iterator
bool getByLabel(std::string const &label, Handle< PROD > &result) const
Definition: Run.h:219
static PFTauRenderPlugin instance
The single EDProduct containing the HLT Prescale Table.
bool exists(std::string const &parameterName) const
checks if a parameter exists
processConfiguration
Definition: Schedule.cc:374
ProcessHistory const & processHistory() const
Definition: Run.cc:99
bool isValid() const
Definition: HandleBase.h:75
ModuleLabelToPathAndFlags moduleLabelToPathAndFlags_
HLTConfigProvider const & hltConfigProvider() const
std::string const & label() const
Definition: InputTag.h:36
std::string const & process() const
Definition: InputTag.h:40
trigger::HLTPrescaleTable hltPrescaleTableRun_
std::string const & instance() const
Definition: InputTag.h:37
edm::InputTag tagL1ExtraNoIsoEG_
bool init(const edm::Run &iRun, const edm::EventSetup &iSetup, const std::string &processName, bool &changed)
HLTPrescaleProvider hltPrescaleProvider_
edm::ParameterSet * l1PSet_
void PATTriggerProducer::produce ( edm::Event iEvent,
const edm::EventSetup iSetup 
)
overrideprivatevirtual

Implements edm::stream::EDProducerBase.

Definition at line 314 of file PATTriggerProducer.cc.

References pat::TriggerAlgorithm::addConditionKey(), pat::TriggerObjectStandAlone::addFilterLabel(), addL1Algos_, pat::TriggerFilter::addObjectKey(), addPathModuleLabels_, pat::TriggerObjectStandAlone::addPathName(), pat::TriggerFilter::addTriggerObjectType(), pat::TriggerObject::addTriggerObjectType(), L1GtUtils::AlgorithmTrigger, assert(), python.rootplot.argparse::category, CenJet, pat::TriggerAlgorithm::decision(), relativeConstraints::empty, edm::InputTag::encode(), relativeConstraints::error, trigger::TriggerObject::eta(), ETM, edm::ParameterSet::exists(), exludeCollections_, firstInRun_, ForJet, edm::EventSetup::get(), edm::Event::getByLabel(), edm::Event::getByToken(), L1GlobalTriggerObjectMaps::ConditionsInAlgorithm::getConditionResult(), L1GlobalTriggerObjectMaps::CombinationsInCondition::getObjectIndex(), edm::ParameterSet::getParameter(), pat::TriggerObjectStandAlone::hasCollection(), HltComparatorCreateWorkflow::hltConfig, hltConfigInit_, HLTPrescaleProvider::hltConfigProvider(), HLTConfigProvider::hltL1GTSeeds(), hltPrescaleLabel_, hltPrescaleProvider_, hltPrescaleTableLumi_, HTM, i, trigger::TriggerObject::id(), HLT_25ns10e33_v2_cff::InputTag, instance, IsoEG, edm::EventBase::isRealData(), edm::HandleBase::isValid(), relval_steps::key, relativeConstraints::keys, L1GtUtils::l1AlgoTechTrigBitNumber(), l1GlobalTriggerObjectMapsToken_, HLTPrescaleProvider::l1GtUtils(), l1PSet_, L1GtUtils::l1Results(), edm::InputTag::label(), diffTwoXMLs::label, labelHltPrescaleTable_, trigger::HLTPrescaleTable::labels(), LogDebug, mainBxOnly_, trigger::TriggerObject::mass(), bookConverter::max, HLTConfigProvider::moduleIndex(), HLTConfigProvider::moduleLabel(), moduleLabelToPathAndFlags_, HLTConfigProvider::moduleType(), Mu, gen::n, pat::TriggerAlgorithm::name(), nameProcess_, cscdqm::h::names, L1GlobalTriggerObjectMaps::CombinationsInCondition::nCombinations(), L1GlobalTriggerObjectMaps::ConditionsInAlgorithm::nConditions(), L1GlobalTriggerObjectMaps::CombinationsInCondition::nObjectsPerCombination(), NoIsoEG, L1GlobalTriggerReadoutSetup::NumberPhysTriggers, L1GlobalTriggerReadoutSetup::NumberTechnicalTriggers, getGTfromDQMFile::obj, onlyStandAlone_, AlCaHLTBitMon_ParallelJobs::p, packPathNames_, packPrescales_, mergeAndRegister::paths, trigger::TriggerObject::phi(), HLTConfigProvider::prescaleLabels(), HLTPrescaleProvider::prescaleSet(), HLTConfigProvider::prescaleSize(), HLTConfigProvider::prescaleTable(), HLTConfigProvider::prescaleValue(), HLTPrescaleProvider::prescaleValuesInDetail(), trigger::TriggerObject::pt(), edm::Event::put(), saveL1Refs_, HLTConfigProvider::saveTags(), HLTConfigProvider::saveTagsModules(), trigger::HLTPrescaleTable::set(), pat::TriggerCondition::setCategory(), pat::TriggerObject::setCollection(), pat::TriggerAlgorithm::setGtlResult(), pat::TriggerAlgorithm::setLogicalExpression(), pat::TriggerFilter::setSaveTags(), pat::TriggerFilter::setStatus(), pat::TriggerFilter::setType(), HLTConfigProvider::size(), edm::TriggerNames::size(), AlCaHLTBitMon_QueryRunRegistry::string, trigger::HLTPrescaleTable::table(), tagL1ExtraCenJet_, tagL1ExtraETM_, tagL1ExtraForJet_, tagL1ExtraHTM_, tagL1ExtraIsoEG_, tagL1ExtraMu_, tagL1ExtraNoIsoEG_, tagL1ExtraTauJet_, tagL1GlobalTriggerObjectMaps_, tagTriggerEvent_, tagTriggerResults_, TauJet, L1GtUtils::TechnicalTrigger, JetPDSkim_cfg::triggerConditions, trigger::TriggerL1CenJet, trigger::TriggerL1ETM, trigger::TriggerL1ForJet, trigger::TriggerL1HTM, trigger::TriggerL1IsoEG, trigger::TriggerL1Mu, trigger::TriggerL1NoIsoEG, trigger::TriggerL1TauJet, edm::TriggerNames::triggerName(), HLTConfigProvider::triggerNames(), edm::Event::triggerNames(), and triggerPaths.

315 {
316 
317  // Terminate, if auto process name determination failed
318  if ( nameProcess_ == "*" ) return;
319 
320  std::auto_ptr< TriggerObjectCollection > triggerObjects( new TriggerObjectCollection() );
321  std::auto_ptr< TriggerObjectStandAloneCollection > triggerObjectsStandAlone( new TriggerObjectStandAloneCollection() );
322  std::auto_ptr< PackedTriggerPrescales > packedPrescales, packedPrescalesL1min, packedPrescalesL1max;
323 
324  // HLT
326 
327  // Get and check HLT event data
328  Handle< trigger::TriggerEvent > handleTriggerEvent;
329  iEvent.getByLabel( tagTriggerEvent_, handleTriggerEvent );
330  Handle< TriggerResults > handleTriggerResults;
331  iEvent.getByLabel( tagTriggerResults_, handleTriggerResults );
332  bool goodHlt( hltConfigInit_ );
333  if ( goodHlt ) {
334  if( ! handleTriggerResults.isValid() ) {
335  LogError( "triggerResultsValid" ) << "TriggerResults product with InputTag '" << tagTriggerResults_.encode() << "' not in event\n"
336  << "No HLT information produced";
337  goodHlt = false;
338  } else if ( ! handleTriggerEvent.isValid() ) {
339  LogError( "triggerEventValid" ) << "trigger::TriggerEvent product with InputTag '" << tagTriggerEvent_.encode() << "' not in event\n"
340  << "No HLT information produced";
341  goodHlt = false;
342  }
343  }
344 
345  // Produce HLT paths and determine status of modules
346 
347  if ( goodHlt ) {
348 
349  // Extract pre-scales
350  // Start from lumi
352  // Try event product, if configured and available
353  if ( ! labelHltPrescaleTable_.empty() ) {
354  Handle< trigger::HLTPrescaleTable > handleHltPrescaleTable;
355  iEvent.getByLabel( InputTag( labelHltPrescaleTable_, "Event", nameProcess_ ), handleHltPrescaleTable );
356  if ( handleHltPrescaleTable.isValid() ) {
357  hltPrescaleTable = trigger::HLTPrescaleTable( handleHltPrescaleTable->set(), handleHltPrescaleTable->labels(), handleHltPrescaleTable->table() );
358  }
359  }
360  // Try event setup, if no product
361  if ( hltPrescaleTable.size() == 0 ) {
362  if ( ! labelHltPrescaleTable_.empty() ) {
363  LogWarning( "hltPrescaleInputTag" ) << "HLTPrescaleTable product with label '" << labelHltPrescaleTable_ << "' not found in process" << nameProcess_ << "\n"
364  << "Using default from event setup";
365  }
366  if ( hltConfig.prescaleSize() > 0 ) {
367  if ( hltPrescaleProvider_.prescaleSet( iEvent, iSetup ) != -1 ) {
368  hltPrescaleTable = trigger::HLTPrescaleTable( hltPrescaleProvider_.prescaleSet( iEvent, iSetup ), hltConfig.prescaleLabels(), hltConfig.prescaleTable() );
369  LogDebug( "hltPrescaleTable" ) << "HLT prescale table found in event setup";
370  } else {
371  LogWarning( "hltPrescaleSet" ) << "HLTPrescaleTable from event setup has error";
372  }
373  }
374  }
375  unsigned set( hltPrescaleTable.set() );
376  if ( hltPrescaleTable.size() > 0 ) {
377  if ( hltPrescaleLabel_.size() > 0 ) {
378  bool foundPrescaleLabel( false );
379  for ( unsigned iLabel = 0; iLabel < hltPrescaleTable.labels().size(); ++iLabel ) {
380  if ( hltPrescaleTable.labels().at( iLabel ) == hltPrescaleLabel_ ) {
381  set = iLabel;
382  foundPrescaleLabel = true;
383  break;
384  }
385  }
386  if ( ! foundPrescaleLabel ) {
387  LogWarning( "hltPrescaleLabel" ) << "HLT prescale label '" << hltPrescaleLabel_ << "' not in prescale table\n"
388  << "Using default";
389  }
390  }
391  } else if ( iEvent.isRealData() ) {
392  if ( ( labelHltPrescaleTable_.empty() && firstInRun_ ) || ! labelHltPrescaleTable_.empty() ) {
393  LogError( "hltPrescaleTable" ) << "No HLT prescale table found\n"
394  << "Using default empty table with all prescales 1";
395  }
396  }
397 
398  const unsigned sizePaths( hltConfig.size() );
399  const unsigned sizeFilters( handleTriggerEvent->sizeFilters() );
400  const unsigned sizeObjects( handleTriggerEvent->sizeObjects() );
401 
402  std::map< std::string, int > moduleStates;
403 
404  if ( ! onlyStandAlone_ ) {
405  std::auto_ptr< TriggerPathCollection > triggerPaths( new TriggerPathCollection() );
406  triggerPaths->reserve( sizePaths );
407  const std::vector<std::string> & pathNames = hltConfig.triggerNames();
408  for ( size_t indexPath = 0; indexPath < sizePaths; ++indexPath ) {
409  const std::string & namePath = pathNames.at( indexPath );
410  unsigned indexLastFilterPathModules( handleTriggerResults->index( indexPath ) + 1 );
411  unsigned indexLastFilterFilters( sizeFilters );
412  while ( indexLastFilterPathModules > 0 ) {
413  --indexLastFilterPathModules;
414  const std::string & labelLastFilterPathModules( hltConfig.moduleLabel( indexPath, indexLastFilterPathModules ) );
415  indexLastFilterFilters = handleTriggerEvent->filterIndex( InputTag( labelLastFilterPathModules, "", nameProcess_ ) );
416  if ( indexLastFilterFilters < sizeFilters ) {
417  if ( hltConfig.moduleType( labelLastFilterPathModules ) == "HLTBool" ) continue;
418  break;
419  }
420  }
421  TriggerPath triggerPath( namePath, indexPath, hltConfig.prescaleValue( set, namePath ), handleTriggerResults->wasrun( indexPath ), handleTriggerResults->accept( indexPath ), handleTriggerResults->error( indexPath ), indexLastFilterPathModules, hltConfig.saveTagsModules( namePath ).size() );
422  // add module names to path and states' map
423  const unsigned sizeModulesPath( hltConfig.size( indexPath ) );
424  assert( indexLastFilterPathModules < sizeModulesPath );
425  std::map< unsigned, std::string > indicesModules;
426  for ( size_t iM = 0; iM < sizeModulesPath; ++iM ) {
427  const std::string nameModule( hltConfig.moduleLabel( indexPath, iM ) );
428  if ( addPathModuleLabels_ ) {
429  triggerPath.addModule( nameModule );
430  }
431  const unsigned indexFilter( handleTriggerEvent->filterIndex( InputTag( nameModule, "", nameProcess_ ) ) );
432  if ( indexFilter < sizeFilters ) {
433  triggerPath.addFilterIndex( indexFilter );
434  }
435  const unsigned slotModule( hltConfig.moduleIndex( indexPath, nameModule ) );
436  indicesModules.insert( std::pair< unsigned, std::string >( slotModule, nameModule ) );
437  }
438  // add L1 seeds
439  const L1SeedCollection l1Seeds( hltConfig.hltL1GTSeeds( namePath ) );
440  for ( L1SeedCollection::const_iterator iSeed = l1Seeds.begin(); iSeed != l1Seeds.end(); ++iSeed ) {
441  triggerPath.addL1Seed( *iSeed );
442  }
443  // store path
444  triggerPaths->push_back( triggerPath );
445  // cache module states to be used for the filters
446  for ( std::map< unsigned, std::string >::const_iterator iM = indicesModules.begin(); iM != indicesModules.end(); ++iM ) {
447  if ( iM->first < indexLastFilterPathModules ) {
448  moduleStates[ iM->second ] = 1;
449  } else if ( iM->first == indexLastFilterPathModules ) {
450  moduleStates[ iM->second ] = handleTriggerResults->accept( indexPath );
451  } else if ( moduleStates.find( iM->second ) == moduleStates.end() ) {
452  moduleStates[ iM->second ] = -1;
453  }
454  }
455  }
456  // Put HLT paths to event
457  iEvent.put( triggerPaths );
458  }
459 
460  // Store used trigger objects and their types for HLT filters
461  // (only active filter(s) available from trigger::TriggerEvent)
462 
463  std::multimap< trigger::size_type, int > objectTypes;
464  std::multimap< trigger::size_type, std::string > filterLabels;
465 
466  for ( size_t iF = 0; iF < sizeFilters; ++iF ) {
467  const std::string nameFilter( handleTriggerEvent->filterLabel( iF ) );
468  const trigger::Keys & keys = handleTriggerEvent->filterKeys( iF );
469  const trigger::Vids & types = handleTriggerEvent->filterIds( iF );
470  assert( types.size() == keys.size() );
471  for ( size_t iK = 0; iK < keys.size(); ++iK ) {
472  filterLabels.insert( std::pair< trigger::size_type, std::string >( keys[ iK ], nameFilter ) );
473  objectTypes.insert( std::pair< trigger::size_type, int >( keys[ iK ], types[ iK ] ) );
474  }
475  }
476 
477  // HLT objects
478 
479  triggerObjects->reserve( onlyStandAlone_ ? 0 : sizeObjects );
480  triggerObjectsStandAlone->reserve( sizeObjects );
481 
482  const trigger::Keys & collectionKeys( handleTriggerEvent->collectionKeys() );
483  std::map< trigger::size_type, trigger::size_type > newObjectKeys;
484  for ( size_t iO = 0, iC = 0, nC = handleTriggerEvent->sizeCollections(); iO < sizeObjects && iC < nC; ++iO ) {
485  const trigger::TriggerObject tobj = handleTriggerEvent->getObjects().at( iO );
486  TriggerObject triggerObject( reco::Particle::PolarLorentzVector(tobj.pt(), tobj.eta(), tobj.phi(), tobj.mass()), tobj.id() );
487  // set collection
488  while ( iO >= collectionKeys[ iC ] ) ++iC; // relies on well ordering of trigger objects with respect to the collections
489  triggerObject.setCollection( handleTriggerEvent->collectionTagEncoded( iC ) );
490  // set filter ID
491  typedef std::multimap< trigger::size_type, int >::const_iterator it_type;
492  for (std::pair<it_type,it_type> trange = objectTypes.equal_range(iO);
493  trange.first != trange.second; ++trange.first) {
494  triggerObject.addTriggerObjectType( trange.first->second );
495  }
496 
497  // stand-alone trigger object
498  TriggerObjectStandAlone triggerObjectStandAlone( triggerObject );
499  // check for excluded collections
500  bool excluded( false );
501  for ( size_t iE = 0; iE < exludeCollections_.size(); ++iE ) {
502  if ( triggerObjectStandAlone.hasCollection( exludeCollections_.at( iE ) ) ) {
503  if ( ! onlyStandAlone_ ) newObjectKeys[ iO ] = trigger::size_type( sizeObjects );
504  excluded = true;
505  break;
506  }
507  }
508  if ( excluded ) continue;
509  typedef std::multimap< trigger::size_type, std::string >::const_iterator it_fl;
510  for (std::pair<it_fl,it_fl> frange = filterLabels.equal_range(iO); frange.first != frange.second; ++frange.first) {
511  triggerObjectStandAlone.addFilterLabel( frange.first->second );
512  const std::vector<ModuleLabelToPathAndFlags::PathAndFlags> & paths = moduleLabelToPathAndFlags_[frange.first->second];
513  for (std::vector<ModuleLabelToPathAndFlags::PathAndFlags>::const_iterator iP = paths.begin(); iP != paths.end(); ++iP) {
514  bool pathFired = handleTriggerResults->wasrun( iP->pathIndex ) && handleTriggerResults->accept( iP->pathIndex );
515  triggerObjectStandAlone.addPathName( iP->pathName, pathFired && iP->lastFilter, pathFired && iP->l3Filter );
516  }
517  }
518 
519  triggerObjectsStandAlone->push_back( triggerObjectStandAlone );
520  if ( ! onlyStandAlone_ ) {
521  triggerObjects->push_back( triggerObject );
522  newObjectKeys[ iO ] = trigger::size_type( triggerObjects->size() - 1 );
523  }
524  }
525 
526  // Re-iterate HLT filters and finally produce them in order to account for optionally skipped objects
527  if ( ! onlyStandAlone_ ) {
528  std::auto_ptr< TriggerFilterCollection > triggerFilters( new TriggerFilterCollection() );
529  triggerFilters->reserve( sizeFilters );
530  for ( size_t iF = 0; iF < sizeFilters; ++iF ) {
531  const std::string nameFilter( handleTriggerEvent->filterTag( iF ).label() );
532  const trigger::Keys & keys = handleTriggerEvent->filterKeys( iF ); // not cached
533  const trigger::Vids & types = handleTriggerEvent->filterIds( iF ); // not cached
534  TriggerFilter triggerFilter( nameFilter );
535  // set filter type
536  const std::string typeFilter( hltConfig.moduleType( nameFilter ) );
537  triggerFilter.setType( typeFilter );
538  triggerFilter.setSaveTags( hltConfig.saveTags( nameFilter ) );
539  // set keys and trigger object types of used objects
540  for ( size_t iK = 0; iK < keys.size(); ++iK ) { // identical to types.size()
541  // check, if current object is excluded
542  if ( newObjectKeys.find( keys.at( iK ) ) != newObjectKeys.end() ) {
543  if ( newObjectKeys[ keys.at( iK ) ] == sizeObjects ) continue;
544  triggerFilter.addObjectKey( newObjectKeys[ keys.at( iK ) ] );
545  triggerFilter.addTriggerObjectType( types.at( iK ) );
546  } else {
547  LogWarning( "triggerObjectKey" ) << "TriggerFilter '" << nameFilter << "' requests non-existing TriggerObject key " << keys.at( iK ) << "\n"
548  << "Skipping object assignment";
549  }
550  }
551  // set status from path info
552  std::map< std::string, int >::iterator iS( moduleStates.find( nameFilter ) );
553  if ( iS != moduleStates.end() ) {
554  if ( ! triggerFilter.setStatus( iS->second ) ) {
555  triggerFilter.setStatus( -1 ); // FIXME different code for "unvalid status determined" needed?
556  }
557  } else {
558  triggerFilter.setStatus( -1 ); // FIXME different code for "unknown" needed?
559  }
560  // store filter
561  triggerFilters->push_back( triggerFilter );
562  }
563  // put HLT filters to event
564  iEvent.put( triggerFilters );
565  }
566 
567  if (packPrescales_) {
568  packedPrescales.reset(new PackedTriggerPrescales(handleTriggerResults));
569  packedPrescalesL1min.reset(new PackedTriggerPrescales(handleTriggerResults));
570  packedPrescalesL1max.reset(new PackedTriggerPrescales(handleTriggerResults));
571  const edm::TriggerNames & names = iEvent.triggerNames(*handleTriggerResults);
572  //std::cout << "Run " << iEvent.id().run() << ", LS " << iEvent.id().luminosityBlock() << ": pset " << set << std::endl;
573  for (unsigned int i = 0, n = names.size(); i < n; ++i) {
574  auto pvdet = hltPrescaleProvider_.prescaleValuesInDetail( iEvent, iSetup, names.triggerName(i) );
575  //int hltprescale = hltConfig_.prescaleValue(set, names.triggerName(i));
576  if (pvdet.first.empty()) {
577  packedPrescalesL1max->addPrescaledTrigger(i, 1);
578  packedPrescalesL1min->addPrescaledTrigger(i, 1);
579  } else {
580  int pmin = -1, pmax = -1;
581  for (const auto & p : pvdet.first) {
582  pmax = std::max(pmax, p.second);
583  if (p.second > 0 && (pmin == -1 || pmin > p.second)) pmin = p.second;
584  }
585  packedPrescalesL1max->addPrescaledTrigger(i, pmax);
586  packedPrescalesL1min->addPrescaledTrigger(i, pmin);
587  //std::cout << "\tTrigger " << names.triggerName(i) << ", L1 ps " << pmin << "-" << pmax << ", HLT ps " << hltprescale << std::endl;
588  }
589  packedPrescales->addPrescaledTrigger(i, pvdet.second);
590  //assert( hltprescale == pvdet.second );
591  }
592  iEvent.put( packedPrescales );
593  iEvent.put( packedPrescalesL1max, "l1max" );
594  iEvent.put( packedPrescalesL1min, "l1min" );
595  }
596 
597  } // if ( goodHlt )
598 
599  // L1 objects
600  // (needs to be done after HLT objects, since their x-links with filters rely on their collection keys)
601 
602  // map for assignments of objects to conditions
603  std::map< L1GtObject, std::vector< unsigned > > l1ObjectTypeMap;
604  if ( ! tagL1ExtraMu_.label().empty() ) {
606  iEvent.getByLabel( tagL1ExtraMu_, handleL1ExtraMu );
607  if ( handleL1ExtraMu.isValid() ) {
608  std::vector< unsigned > muKeys;
609  for ( size_t l1Mu = 0; l1Mu < handleL1ExtraMu->size(); ++l1Mu ) {
610  if ( mainBxOnly_ && handleL1ExtraMu->at( l1Mu ).bx() != 0 ) continue;
611  TriggerObject triggerObject;
612  if ( saveL1Refs_ ) {
613  const reco::CandidateBaseRef leafCandRef( l1extra::L1MuonParticleRef( handleL1ExtraMu, l1Mu ) );
614  triggerObject = TriggerObject( leafCandRef );
615  } else {
616  const reco::LeafCandidate leafCandidate( *( handleL1ExtraMu->at( l1Mu ).reco::LeafCandidate::clone() ) );
617  triggerObject = TriggerObject( leafCandidate );
618  }
619  triggerObject.setCollection( tagL1ExtraMu_ );
621  if ( ! onlyStandAlone_ ) triggerObjects->push_back( triggerObject );
622  TriggerObjectStandAlone triggerObjectStandAlone( triggerObject );
623  triggerObjectsStandAlone->push_back( triggerObjectStandAlone );
624  if ( handleL1ExtraMu->at( l1Mu ).bx() == 0 ) muKeys.push_back( triggerObjectsStandAlone->size() - 1 );
625  }
626  l1ObjectTypeMap.insert( std::make_pair( Mu, muKeys ) );
627  } else LogError( "l1ExtraValid" ) << "l1extra::L1MuonParticleCollection product with InputTag '" << tagL1ExtraMu_.encode() << "' not in event";
628  }
629  if ( ! tagL1ExtraNoIsoEG_.label().empty() ) {
630  Handle< l1extra::L1EmParticleCollection > handleL1ExtraNoIsoEG;
631  iEvent.getByLabel( tagL1ExtraNoIsoEG_, handleL1ExtraNoIsoEG );
632  if ( handleL1ExtraNoIsoEG.isValid() ) {
633  std::vector< unsigned > noIsoEGKeys;
634  for ( size_t l1NoIsoEG = 0; l1NoIsoEG < handleL1ExtraNoIsoEG->size(); ++l1NoIsoEG ) {
635  if ( mainBxOnly_ && handleL1ExtraNoIsoEG->at( l1NoIsoEG ).bx() != 0 ) continue;
636  TriggerObject triggerObject;
637  if ( saveL1Refs_ ) {
638  const reco::CandidateBaseRef leafCandRef( l1extra::L1EmParticleRef( handleL1ExtraNoIsoEG, l1NoIsoEG ) );
639  triggerObject = TriggerObject( leafCandRef );
640  } else {
641  const reco::LeafCandidate leafCandidate( *( handleL1ExtraNoIsoEG->at( l1NoIsoEG ).reco::LeafCandidate::clone() ) );
642  triggerObject = TriggerObject( leafCandidate );
643  }
644  triggerObject.setCollection( tagL1ExtraNoIsoEG_ );
646  if ( ! onlyStandAlone_ ) triggerObjects->push_back( triggerObject );
647  TriggerObjectStandAlone triggerObjectStandAlone( triggerObject );
648  triggerObjectsStandAlone->push_back( triggerObjectStandAlone );
649  if ( handleL1ExtraNoIsoEG->at( l1NoIsoEG ).bx() == 0 ) noIsoEGKeys.push_back( triggerObjectsStandAlone->size() - 1 );
650  }
651  l1ObjectTypeMap.insert( std::make_pair( NoIsoEG, noIsoEGKeys ) );
652  } else LogError( "l1ExtraValid" ) << "l1extra::L1EmParticleCollection product with InputTag '" << tagL1ExtraNoIsoEG_.encode() << "' not in event";
653  }
654  if ( ! tagL1ExtraIsoEG_.label().empty() ) {
656  iEvent.getByLabel( tagL1ExtraIsoEG_, handleL1ExtraIsoEG );
657  if ( handleL1ExtraIsoEG.isValid() ) {
658  std::vector< unsigned > isoEGKeys;
659  for ( size_t l1IsoEG = 0; l1IsoEG < handleL1ExtraIsoEG->size(); ++l1IsoEG ) {
660  if ( mainBxOnly_ && handleL1ExtraIsoEG->at( l1IsoEG ).bx() != 0 ) continue;
661  TriggerObject triggerObject;
662  if ( saveL1Refs_ ) {
663  const reco::CandidateBaseRef leafCandRef( l1extra::L1EmParticleRef( handleL1ExtraIsoEG, l1IsoEG ) );
664  triggerObject = TriggerObject( leafCandRef );
665  } else {
666  const reco::LeafCandidate leafCandidate( *( handleL1ExtraIsoEG->at( l1IsoEG ).reco::LeafCandidate::clone() ) );
667  triggerObject = TriggerObject( leafCandidate );
668  }
669  triggerObject.setCollection( tagL1ExtraIsoEG_ );
671  if ( ! onlyStandAlone_ ) triggerObjects->push_back( triggerObject );
672  TriggerObjectStandAlone triggerObjectStandAlone( triggerObject );
673  triggerObjectsStandAlone->push_back( triggerObjectStandAlone );
674  if ( handleL1ExtraIsoEG->at( l1IsoEG ).bx() == 0 ) isoEGKeys.push_back( triggerObjectsStandAlone->size() - 1 );
675  }
676  l1ObjectTypeMap.insert( std::make_pair( IsoEG, isoEGKeys ) );
677  } else LogError( "l1ExtraValid" ) << "l1extra::L1EmParticleCollection product with InputTag '" << tagL1ExtraIsoEG_.encode() << "' not in event";
678  }
679  if ( ! tagL1ExtraCenJet_.label().empty() ) {
680  Handle< l1extra::L1JetParticleCollection > handleL1ExtraCenJet;
681  iEvent.getByLabel( tagL1ExtraCenJet_, handleL1ExtraCenJet );
682  if ( handleL1ExtraCenJet.isValid() ) {
683  std::vector< unsigned > cenJetKeys;
684  for ( size_t l1CenJet = 0; l1CenJet < handleL1ExtraCenJet->size(); ++l1CenJet ) {
685  if ( mainBxOnly_ && handleL1ExtraCenJet->at( l1CenJet ).bx() != 0 ) continue;
686  TriggerObject triggerObject;
687  if ( saveL1Refs_ ) {
688  const reco::CandidateBaseRef leafCandRef( l1extra::L1JetParticleRef( handleL1ExtraCenJet, l1CenJet ) );
689  triggerObject = TriggerObject( leafCandRef );
690  } else {
691  const reco::LeafCandidate leafCandidate( *( handleL1ExtraCenJet->at( l1CenJet ).reco::LeafCandidate::clone() ) );
692  triggerObject = TriggerObject( leafCandidate );
693  }
694  triggerObject.setCollection( tagL1ExtraCenJet_ );
696  if ( ! onlyStandAlone_ ) triggerObjects->push_back( triggerObject );
697  TriggerObjectStandAlone triggerObjectStandAlone( triggerObject );
698  triggerObjectsStandAlone->push_back( triggerObjectStandAlone );
699  if ( handleL1ExtraCenJet->at( l1CenJet ).bx() == 0 ) cenJetKeys.push_back( triggerObjectsStandAlone->size() - 1 );
700  }
701  l1ObjectTypeMap.insert( std::make_pair( CenJet, cenJetKeys ) );
702  } else LogError( "l1ExtraValid" ) << "l1extra::L1JetParticleCollection product with InputTag '" << tagL1ExtraCenJet_.encode() << "' not in event";
703  }
704  if ( ! tagL1ExtraForJet_.label().empty() ) {
705  Handle< l1extra::L1JetParticleCollection > handleL1ExtraForJet;
706  iEvent.getByLabel( tagL1ExtraForJet_, handleL1ExtraForJet );
707  if ( handleL1ExtraForJet.isValid() ) {
708  std::vector< unsigned > forJetKeys;
709  for ( size_t l1ForJet = 0; l1ForJet < handleL1ExtraForJet->size(); ++l1ForJet ) {
710  if ( mainBxOnly_ && handleL1ExtraForJet->at( l1ForJet ).bx() != 0 ) continue;
711  TriggerObject triggerObject;
712  if ( saveL1Refs_ ) {
713  const reco::CandidateBaseRef leafCandRef( l1extra::L1JetParticleRef( handleL1ExtraForJet, l1ForJet ) );
714  triggerObject = TriggerObject( leafCandRef );
715  } else {
716  const reco::LeafCandidate leafCandidate( *( handleL1ExtraForJet->at( l1ForJet ).reco::LeafCandidate::clone() ) );
717  triggerObject = TriggerObject( leafCandidate );
718  }
719  triggerObject.setCollection( tagL1ExtraForJet_ );
721  if ( ! onlyStandAlone_ ) triggerObjects->push_back( triggerObject );
722  TriggerObjectStandAlone triggerObjectStandAlone( triggerObject );
723  triggerObjectsStandAlone->push_back( triggerObjectStandAlone );
724  if ( handleL1ExtraForJet->at( l1ForJet ).bx() == 0 ) forJetKeys.push_back( triggerObjectsStandAlone->size() - 1 );
725  }
726  l1ObjectTypeMap.insert( std::make_pair( ForJet, forJetKeys ) );
727  } else LogError( "l1ExtraValid" ) << "l1extra::L1JetParticleCollection product with InputTag '" << tagL1ExtraForJet_.encode() << "' not in event";
728  }
729  if ( ! tagL1ExtraTauJet_.label().empty() ) {
730  Handle< l1extra::L1JetParticleCollection > handleL1ExtraTauJet;
731  iEvent.getByLabel( tagL1ExtraTauJet_, handleL1ExtraTauJet );
732  if ( handleL1ExtraTauJet.isValid() ) {
733  std::vector< unsigned > tauJetKeys;
734  for ( size_t l1TauJet = 0; l1TauJet < handleL1ExtraTauJet->size(); ++l1TauJet ) {
735  if ( mainBxOnly_ && handleL1ExtraTauJet->at( l1TauJet ).bx() != 0 ) continue;
736  TriggerObject triggerObject;
737  if ( saveL1Refs_ ) {
738  const reco::CandidateBaseRef leafCandRef( l1extra::L1JetParticleRef( handleL1ExtraTauJet, l1TauJet ) );
739  triggerObject = TriggerObject( leafCandRef );
740  } else {
741  const reco::LeafCandidate leafCandidate( *( handleL1ExtraTauJet->at( l1TauJet ).reco::LeafCandidate::clone() ) );
742  triggerObject = TriggerObject( leafCandidate );
743  }
744  triggerObject.setCollection( tagL1ExtraTauJet_ );
746  if ( ! onlyStandAlone_ ) triggerObjects->push_back( triggerObject );
747  TriggerObjectStandAlone triggerObjectStandAlone( triggerObject );
748  triggerObjectsStandAlone->push_back( triggerObjectStandAlone );
749  if ( handleL1ExtraTauJet->at( l1TauJet ).bx() == 0 ) tauJetKeys.push_back( triggerObjectsStandAlone->size() - 1 );
750  }
751  l1ObjectTypeMap.insert( std::make_pair( TauJet, tauJetKeys ) );
752  } else LogError( "l1ExtraValid" ) << "l1extra::L1JetParticleCollection product with InputTag '" << tagL1ExtraTauJet_.encode() << "' not in event";
753  }
754  if ( ! tagL1ExtraETM_ .label().empty()) {
756  iEvent.getByLabel( tagL1ExtraETM_, handleL1ExtraETM );
757  if ( handleL1ExtraETM.isValid() ) {
758  std::vector< unsigned > etmKeys;
759  for ( size_t l1ETM = 0; l1ETM < handleL1ExtraETM->size(); ++l1ETM ) {
760  if ( mainBxOnly_ && handleL1ExtraETM->at( l1ETM ).bx() != 0 ) continue;
761  TriggerObject triggerObject;
762  if ( saveL1Refs_ ) {
763  const reco::CandidateBaseRef leafCandRef( l1extra::L1EtMissParticleRef( handleL1ExtraETM, l1ETM ) );
764  triggerObject = TriggerObject( leafCandRef );
765  } else {
766  const reco::LeafCandidate leafCandidate( *( handleL1ExtraETM->at( l1ETM ).reco::LeafCandidate::clone() ) );
767  triggerObject = TriggerObject( leafCandidate );
768  }
769  triggerObject.setCollection( tagL1ExtraETM_ );
771  if ( ! onlyStandAlone_ ) triggerObjects->push_back( triggerObject );
772  TriggerObjectStandAlone triggerObjectStandAlone( triggerObject );
773  triggerObjectsStandAlone->push_back( triggerObjectStandAlone );
774  if ( handleL1ExtraETM->at( l1ETM ).bx() == 0 ) etmKeys.push_back( triggerObjectsStandAlone->size() - 1 );
775  }
776  l1ObjectTypeMap.insert( std::make_pair( ETM, etmKeys ) );
777  } else LogError( "l1ExtraValid" ) << "l1extra::L1EtMissParticleCollection product with InputTag '" << tagL1ExtraETM_.encode() << "' not in event";
778  }
779  if ( ! tagL1ExtraHTM_.label().empty() ) {
781  iEvent.getByLabel( tagL1ExtraHTM_, handleL1ExtraHTM );
782  if ( handleL1ExtraHTM.isValid() ) {
783  std::vector< unsigned > htmKeys;
784  for ( size_t l1HTM = 0; l1HTM < handleL1ExtraHTM->size(); ++l1HTM ) {
785  if ( mainBxOnly_ && handleL1ExtraHTM->at( l1HTM ).bx() != 0 ) continue;
786  TriggerObject triggerObject;
787  if ( saveL1Refs_ ) {
788  const reco::CandidateBaseRef leafCandRef( l1extra::L1EtMissParticleRef( handleL1ExtraHTM, l1HTM ) );
789  triggerObject = TriggerObject( leafCandRef );
790  } else {
791  const reco::LeafCandidate leafCandidate( *( handleL1ExtraHTM->at( l1HTM ).reco::LeafCandidate::clone() ) );
792  triggerObject = TriggerObject( leafCandidate );
793  }
794  triggerObject.setCollection( tagL1ExtraHTM_ );
796  if ( ! onlyStandAlone_ ) triggerObjects->push_back( triggerObject );
797  TriggerObjectStandAlone triggerObjectStandAlone( triggerObject );
798  triggerObjectsStandAlone->push_back( triggerObjectStandAlone );
799  if ( handleL1ExtraHTM->at( l1HTM ).bx() == 0 ) htmKeys.push_back( triggerObjectsStandAlone->size() - 1 );
800  }
801  l1ObjectTypeMap.insert( std::make_pair( HTM, htmKeys ) );
802  } else LogError( "l1ExtraValid" ) << "l1extra::L1EtMissParticleCollection product with InputTag '" << tagL1ExtraHTM_.encode() << "' not in event";
803  }
804 
805  // Put trigger objects to event
806  if ( ! onlyStandAlone_ ) iEvent.put( triggerObjects );
807 
808  // L1 algorithms
809  if ( ! onlyStandAlone_ ) {
810  std::auto_ptr< TriggerAlgorithmCollection > triggerAlgos( new TriggerAlgorithmCollection() );
811  std::auto_ptr< TriggerConditionCollection > triggerConditions( new TriggerConditionCollection() );
812  if ( addL1Algos_ ) {
813  // create trigger object types transalation map (yes, it's ugly!)
814  std::map< L1GtObject, trigger::TriggerObjectType > mapObjectTypes;
815  mapObjectTypes.insert( std::make_pair( Mu , trigger::TriggerL1Mu ) );
816  mapObjectTypes.insert( std::make_pair( NoIsoEG, trigger::TriggerL1NoIsoEG ) );
817  mapObjectTypes.insert( std::make_pair( IsoEG , trigger::TriggerL1IsoEG ) );
818  mapObjectTypes.insert( std::make_pair( CenJet , trigger::TriggerL1CenJet ) );
819  mapObjectTypes.insert( std::make_pair( ForJet , trigger::TriggerL1ForJet ) );
820  mapObjectTypes.insert( std::make_pair( TauJet , trigger::TriggerL1TauJet ) );
821  mapObjectTypes.insert( std::make_pair( ETM , trigger::TriggerL1ETM ) );
822  mapObjectTypes.insert( std::make_pair( HTM , trigger::TriggerL1HTM ) );
823  // get and cache L1 menu
824  L1GtUtils const& l1GtUtils = hltPrescaleProvider_.l1GtUtils();
825  ESHandle< L1GtTriggerMenu > handleL1GtTriggerMenu;
826  iSetup.get< L1GtTriggerMenuRcd >().get( handleL1GtTriggerMenu );
827  auto const & l1GtAlgorithms = handleL1GtTriggerMenu->gtAlgorithmMap();
828  auto const & l1GtTechTriggers = handleL1GtTriggerMenu->gtTechnicalTriggerMap();
829  auto const & l1GtConditionsVector = handleL1GtTriggerMenu->gtConditionMap();
830  // cache conditions in one single condition map
831  ConditionMap l1GtConditions;
832  for ( size_t iCv = 0; iCv < l1GtConditionsVector.size(); ++iCv ) {
833  l1GtConditions.insert( l1GtConditionsVector.at( iCv ).begin(), l1GtConditionsVector.at( iCv ).end() );
834  }
835  triggerAlgos->reserve( l1GtAlgorithms.size() + l1GtTechTriggers.size() );
836  Handle< L1GlobalTriggerObjectMaps > handleL1GlobalTriggerObjectMaps;
837  iEvent.getByToken( l1GlobalTriggerObjectMapsToken_, handleL1GlobalTriggerObjectMaps );
838  if( ! handleL1GlobalTriggerObjectMaps.isValid() ) {
839  LogError( "l1ObjectMap" ) << "L1GlobalTriggerObjectMaps product with InputTag '" << tagL1GlobalTriggerObjectMaps_.encode() << "' not in event\n"
840  << "No L1 objects and GTL results available for physics algorithms";
841  }
842  handleL1GlobalTriggerObjectMaps->consistencyCheck();
843  if ( firstInRun_ ) {
844  l1PSet_ = ( ParameterSet* )( pset::Registry::instance()->getMapped(handleL1GlobalTriggerObjectMaps->namesParameterSetID()) );
845  if (l1PSet_ == 0) {
846  LogError( "l1ObjectMap" ) << "ParameterSet registry not available\n"
847  << "Skipping conditions for all L1 physics algorithm names in this run";
848  }
849  } else {
850  if (l1PSet_ == 0) {
851  LogInfo( "l1ObjectMap" ) << "ParameterSet registry not available\n"
852  << "Skipping conditions for all L1 physics algorithm names in this event";
853  }
854  }
855  // physics algorithms
856  for ( CItAlgo iAlgo = l1GtAlgorithms.begin(); iAlgo != l1GtAlgorithms.end(); ++iAlgo ) {
857  const std::string & algoName( iAlgo->second.algoName() );
858  if ( ! ( iAlgo->second.algoBitNumber() < int( L1GlobalTriggerReadoutSetup::NumberPhysTriggers ) ) ) {
859  LogError( "l1Algo" ) << "L1 physics algorithm '" << algoName << "' has bit number " << iAlgo->second.algoBitNumber() << " >= " << L1GlobalTriggerReadoutSetup::NumberPhysTriggers << "\n"
860  << "Skipping";
861  continue;
862  }
864  int bit;
865  if ( ! l1GtUtils.l1AlgoTechTrigBitNumber( algoName, category, bit ) ) {
866  LogError( "l1Algo" ) << "L1 physics algorithm '" << algoName << "' not found in the L1 menu\n"
867  << "Skipping";
868  continue;
869  }
870  if ( category != L1GtUtils::AlgorithmTrigger ) {
871  LogError( "l1Algo" ) << "L1 physics algorithm '" << algoName << "' does not have category 'AlgorithmTrigger' from 'L1GtUtils'\n"
872  << "Skipping";
873  continue;
874  }
875  bool decisionBeforeMask;
876  bool decisionAfterMask;
877  int prescale;
878  int mask;
879  int error( l1GtUtils.l1Results( iEvent, algoName, decisionBeforeMask, decisionAfterMask, prescale, mask ) );
880  if ( error ) {
881  LogError( "l1Algo" ) << "L1 physics algorithm '" << algoName << "' decision has error code " << error << " from 'L1GtUtils'\n"
882  << "Skipping";
883  continue;
884  }
885  TriggerAlgorithm triggerAlgo( algoName, iAlgo->second.algoAlias(), category == L1GtUtils::TechnicalTrigger, (unsigned)bit, (unsigned)prescale, (bool)mask, decisionBeforeMask, decisionAfterMask );
886  triggerAlgo.setLogicalExpression( iAlgo->second.algoLogicalExpression() );
887  // GTL result and used conditions in physics algorithm
888  if( ! handleL1GlobalTriggerObjectMaps.isValid() ) {
889  triggerAlgos->push_back( triggerAlgo );
890  continue; // LogWarning already earlier (before loop)
891  }
892  if ( ! handleL1GlobalTriggerObjectMaps->algorithmExists(bit)) {
893  LogError( "l1ObjectMap" ) << "L1 physics algorithm '" << algoName << "' is missing in L1GlobalTriggerObjectMaps\n"
894  << "Skipping conditions and GTL result";
895  triggerAlgos->push_back( triggerAlgo );
896  continue;
897  }
898  bool algorithmResult = handleL1GlobalTriggerObjectMaps->algorithmResult(bit);
899 // if ( ( algorithmResult != decisionBeforeMask ) && ( decisionBeforeMask == true || prescale == 1 ) ) {
900  if ( ( algorithmResult != decisionBeforeMask ) && ( decisionBeforeMask == true ) ) { // FIXME: understand the difference for un-prescaled algos 118, 119, 123
901  LogInfo( "l1ObjectMap" ) << "L1 physics algorithm '" << algoName << "' with different decisions in\n"
902  << "L1GlobalTriggerObjectMaps (GTL result) : " << algorithmResult << "\n"
903  << "L1GlobalTriggerReadoutRecord (decision before mask): " << decisionBeforeMask;
904  }
905  triggerAlgo.setGtlResult( algorithmResult );
906  // conditions in algorithm
907  L1GlobalTriggerObjectMaps::ConditionsInAlgorithm conditions = handleL1GlobalTriggerObjectMaps->getConditionsInAlgorithm(bit);
908  if (l1PSet_ == 0) {
909  triggerAlgos->push_back( triggerAlgo );
910  continue;
911  }
912  if (!l1PSet_->exists(algoName)) {
913  if ( firstInRun_ ) {
914  LogError( "l1ObjectMap" ) << "L1 physics algorithm name '" << algoName << "' not available in ParameterSet registry\n"
915  << "Skipping conditions for this algorithm in this run";
916  } else {
917  LogInfo( "l1ObjectMap" ) << "L1 physics algorithm name '" << algoName << "' not available in ParameterSet registry\n"
918  << "Skipping conditions for this algorithm in this event";
919  }
920  triggerAlgos->push_back( triggerAlgo );
921  continue;
922  }
923  std::vector<std::string> conditionNames( l1PSet_->getParameter<std::vector<std::string> >(algoName) );
924 
925  for ( unsigned iT = 0; iT < conditionNames.size(); ++iT ) {
926  size_t key( triggerConditions->size() );
927  for ( size_t iC = 0; iC < triggerConditions->size(); ++iC ) {
928  if ( conditionNames.at(iT) == triggerConditions->at( iC ).name() ) {
929  key = iC;
930  break;
931  }
932  }
933  if ( key == triggerConditions->size() ) {
934  if ( iT >= conditions.nConditions() ) {
935  LogError( "l1CondMap" ) << "More condition names from ParameterSet registry than the " << conditions.nConditions() << " conditions in L1GlobalTriggerObjectMaps\n"
936  << "Skipping condition " << conditionNames.at(iT) << " in algorithm " << algoName;
937  break;
938  }
939  TriggerCondition triggerCond( conditionNames[iT], conditions.getConditionResult(iT) );
940  if ( l1GtConditions.find( triggerCond.name() ) != l1GtConditions.end() ) {
941  triggerCond.setCategory( l1GtConditions[ triggerCond.name() ]->condCategory() );
942  triggerCond.setType( l1GtConditions[ triggerCond.name() ]->condType() );
943  const std::vector< L1GtObject > l1ObjectTypes( l1GtConditions[ triggerCond.name() ]->objectType() );
944  for ( size_t iType = 0 ; iType < l1ObjectTypes.size(); ++iType ) {
945  triggerCond.addTriggerObjectType( mapObjectTypes[ l1ObjectTypes.at( iType ) ] );
946  }
947  // objects in condition
948  L1GlobalTriggerObjectMaps::CombinationsInCondition combinations = handleL1GlobalTriggerObjectMaps->getCombinationsInCondition(bit, iT);
949  for ( size_t iVV = 0; iVV < combinations.nCombinations(); ++iVV ) {
950  for ( size_t iV = 0; iV < combinations.nObjectsPerCombination(); ++iV ) {
951 
952  unsigned objectIndex = combinations.getObjectIndex(iVV, iV);
953  if ( iV >= l1ObjectTypes.size() ) {
954  LogError( "l1CondMap" ) << "Index " << iV << " in combinations vector overshoots size " << l1ObjectTypes.size() << " of types vector in conditions map\n"
955  << "Skipping object key in condition " << triggerCond.name();
956  } else if ( l1ObjectTypeMap.find( l1ObjectTypes.at( iV ) ) != l1ObjectTypeMap.end() ) {
957  if ( objectIndex >= l1ObjectTypeMap[ l1ObjectTypes.at( iV ) ].size() ) {
958  LogError( "l1CondMap" ) << "Index " << objectIndex << " in combination overshoots number " << l1ObjectTypeMap[ l1ObjectTypes.at( iV ) ].size() << "of according trigger objects\n"
959  << "Skipping object key in condition " << triggerCond.name();
960  }
961  const unsigned objectKey( l1ObjectTypeMap[ l1ObjectTypes.at( iV ) ].at( objectIndex ) );
962  triggerCond.addObjectKey( objectKey );
963  // add current condition and algorithm also to the according stand-alone trigger object
964  triggerObjectsStandAlone->at( objectKey ).addAlgorithmName( triggerAlgo.name(), ( triggerAlgo.decision() && triggerCond.wasAccept() ) );
965  triggerObjectsStandAlone->at( objectKey ).addConditionName( triggerCond.name() );
966  }
967  }
968  }
969  } else {
970  LogWarning( "l1CondMap" ) << "L1 conditions '" << triggerCond.name() << "' not found in the L1 menu\n"
971  << "Remains incomplete";
972  }
973  triggerConditions->push_back( triggerCond );
974  }
975  triggerAlgo.addConditionKey( key );
976  }
977  triggerAlgos->push_back( triggerAlgo );
978  }
979  // technical triggers
980  for ( CItAlgo iAlgo = l1GtTechTriggers.begin(); iAlgo != l1GtTechTriggers.end(); ++iAlgo ) {
981  const std::string & algoName( iAlgo->second.algoName() );
982  if ( ! ( iAlgo->second.algoBitNumber() < int( L1GlobalTriggerReadoutSetup::NumberTechnicalTriggers ) ) ) {
983  LogError( "l1Algo" ) << "L1 technical trigger '" << algoName << "' has bit number " << iAlgo->second.algoBitNumber() << " >= " << L1GlobalTriggerReadoutSetup::NumberTechnicalTriggers << "\n"
984  << "Skipping";
985  continue;
986  }
988  int bit;
989  if ( ! l1GtUtils.l1AlgoTechTrigBitNumber( algoName, category, bit ) ) {
990  LogError( "l1Algo" ) << "L1 technical trigger '" << algoName << "' not found in the L1 menu\n"
991  << "Skipping";
992  continue;
993  }
994  if ( category != L1GtUtils::TechnicalTrigger ) {
995  LogError( "l1Algo" ) << "L1 technical trigger '" << algoName << "' does not have category 'TechnicalTrigger' from 'L1GtUtils'\n"
996  << "Skipping";
997  continue;
998  }
999  bool decisionBeforeMask;
1000  bool decisionAfterMask;
1001  int prescale;
1002  int mask;
1003  int error( l1GtUtils.l1Results( iEvent, algoName, decisionBeforeMask, decisionAfterMask, prescale, mask ) );
1004  if ( error ) {
1005  LogError( "l1Algo" ) << "L1 technical trigger '" << algoName << "' decision has error code " << error << " from 'L1GtUtils'\n"
1006  << "Skipping";
1007  continue;
1008  }
1009  TriggerAlgorithm triggerAlgo( algoName, iAlgo->second.algoAlias(), category == L1GtUtils::TechnicalTrigger, (unsigned)bit, (unsigned)prescale, (bool)mask, decisionBeforeMask, decisionAfterMask );
1010  triggerAlgo.setLogicalExpression( iAlgo->second.algoLogicalExpression() );
1011  triggerAlgos->push_back( triggerAlgo );
1012  }
1013  }
1014 
1015  // Put L1 algorithms and conditions to event
1016  iEvent.put( triggerAlgos );
1017  iEvent.put( triggerConditions );
1018  }
1019 
1020 
1021  if (packPathNames_) {
1022  const edm::TriggerNames & names = iEvent.triggerNames(*handleTriggerResults);
1023  for (pat::TriggerObjectStandAlone &obj : *triggerObjectsStandAlone) {
1024  obj.packPathNames(names);
1025  }
1026  }
1027  // Put (finally) stand-alone trigger objects to event
1028  iEvent.put( triggerObjectsStandAlone );
1029 
1030  firstInRun_ = false;
1031 
1032 }
#define LogDebug(id)
unsigned int set() const
low-level const accessors for data members
unsigned int size() const
number of trigger paths in trigger table
T getParameter(std::string const &) const
unsigned char getObjectIndex(unsigned combination, unsigned object) const
int i
Definition: DBlmapReader.cc:9
bool saveTags(const std::string &module) const
Is module an L3 filter (ie, tracked saveTags=true)
virtual edm::TriggerNames const & triggerNames(edm::TriggerResults const &triggerResults) const
Definition: Event.cc:215
std::pair< std::vector< std::pair< std::string, int > >, int > prescaleValuesInDetail(const edm::Event &iEvent, const edm::EventSetup &iSetup, const std::string &trigger)
Analysis-level L1 trigger algorithm class.
const std::string moduleType(const std::string &module) const
C++ class name of module.
int id() const
getters
Definition: TriggerObject.h:55
Definition: L1GtObject.h:39
static const HistoName names[]
static PFTauRenderPlugin instance
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:462
enum start value shifted to 81 so as to avoid clashes with PDG codes
float phi() const
Definition: TriggerObject.h:58
std::map< std::string, L1GtCondition * > ConditionMap
map containing the conditions
Definition: L1GtObject.h:36
tuple triggerConditions
const bool l1AlgoTechTrigBitNumber(const std::string &nameAlgoTechTrig, TriggerCategory &trigCategory, int &bitNumber) const
Definition: L1GtUtils.cc:376
static const unsigned int NumberTechnicalTriggers
const std::vector< std::vector< std::pair< bool, std::string > > > & hltL1GTSeeds() const
assert(m_qm.get())
const std::vector< std::string > & triggerNames() const
names of trigger paths
The single EDProduct containing the HLT Prescale Table.
bool exists(std::string const &parameterName) const
checks if a parameter exists
Strings::size_type size() const
Definition: TriggerNames.cc:39
float eta() const
Definition: TriggerObject.h:57
std::vector< TriggerObjectStandAlone > TriggerObjectStandAloneCollection
Collection of TriggerObjectStandAlone.
const std::map< std::string, std::vector< unsigned int > > & prescaleTable() const
Analysis-level trigger object class.
Definition: TriggerObject.h:47
uint16_t size_type
bool isRealData() const
Definition: EventBase.h:63
const std::map< std::string, std::vector< unsigned int > > & table() const
std::string encode() const
Definition: InputTag.cc:164
const std::string & moduleLabel(unsigned int trigger, unsigned int module) const
void put(edm::Event &evt, double value, const char *instanceName)
AlgorithmMap::const_iterator CItAlgo
iterators through map containing the algorithms
const std::vector< std::string > & saveTagsModules(unsigned int trigger) const
Single trigger physics object (e.g., an isolated muon)
Definition: TriggerObject.h:22
U second(std::pair< T, U > const &p)
std::vector< TriggerAlgorithm > TriggerAlgorithmCollection
Collection of TriggerAlgorithm.
Analysis-level HLTrigger filter class.
Definition: TriggerFilter.h:35
unsigned int prescaleValue(unsigned int set, const std::string &trigger) const
HLT prescale value in specific prescale set for a specific trigger path.
unsigned int moduleIndex(unsigned int trigger, const std::string &module) const
slot position of module on trigger path (0 to size-1)
void setCollection(const std::string &collName)
Methods.
Definition: TriggerObject.h:84
Analysis-level L1 trigger condition class.
Analysis-level HLTrigger path class.
Definition: TriggerPath.h:39
OrphanHandle< PROD > put(std::auto_ptr< PROD > product)
Put a new product.
Definition: Event.h:121
const std::vector< std::string > & labels() const
trigger::HLTPrescaleTable hltPrescaleTableLumi_
const int l1Results(const edm::Event &iEvent, const std::string &nameAlgoTechTrig, bool &decisionBeforeMask, bool &decisionAfterMask, int &prescaleFactor, int &triggerMask) const
Definition: L1GtUtils.cc:713
static const unsigned int NumberPhysTriggers
bool isValid() const
Definition: HandleBase.h:75
string key
FastSim: produces sample of signal events, overlayed with premixed minbias events.
edm::EDGetTokenT< L1GlobalTriggerObjectMaps > l1GlobalTriggerObjectMapsToken_
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:418
int prescaleSet(const edm::Event &iEvent, const edm::EventSetup &iSetup)
edm::InputTag tagL1GlobalTriggerObjectMaps_
std::vector< TriggerFilter > TriggerFilterCollection
Collection of TriggerFilter.
ModuleLabelToPathAndFlags moduleLabelToPathAndFlags_
L1GtUtils const & l1GtUtils() const
void setCategory(L1GtConditionCategory category)
Set the condition category.
Definition: L1GtObject.h:30
std::string const & triggerName(unsigned int index) const
Definition: TriggerNames.cc:27
std::vector< size_type > Keys
const T & get() const
Definition: EventSetup.h:56
const std::vector< std::string > & prescaleLabels() const
low-level data member access
HLTConfigProvider const & hltConfigProvider() const
std::string const & label() const
Definition: InputTag.h:36
if(dp >Float(M_PI)) dp-
std::vector< TriggerObject > TriggerObjectCollection
Collection of TriggerObject.
math::PtEtaPhiMLorentzVector PolarLorentzVector
Lorentz vector.
Definition: Particle.h:23
unsigned int prescaleSize() const
std::vector< TriggerCondition > TriggerConditionCollection
Collection of TriggerCondition.
std::vector< std::string > exludeCollections_
std::vector< L1Seed > L1SeedCollection
Collection of L1Seed.
Definition: TriggerPath.h:37
float mass() const
Definition: TriggerObject.h:59
void addTriggerObjectType(trigger::TriggerObjectType triggerObjectType)
Add a new trigger object type identifier.
Definition: TriggerObject.h:87
static std::string const triggerPaths
Definition: EdmProvDump.cc:42
std::vector< TriggerPath > TriggerPathCollection
Collection of TriggerPath.
Definition: TriggerPath.h:153
edm::InputTag tagL1ExtraNoIsoEG_
std::vector< int > Vids
HLTPrescaleProvider hltPrescaleProvider_
Analysis-level trigger object class (stand-alone)
edm::ParameterSet * l1PSet_

Member Data Documentation

bool pat::PATTriggerProducer::addL1Algos_
private

Definition at line 77 of file PATTriggerProducer.h.

Referenced by PATTriggerProducer(), and produce().

bool pat::PATTriggerProducer::addPathModuleLabels_
private

Definition at line 120 of file PATTriggerProducer.h.

Referenced by PATTriggerProducer(), and produce().

bool pat::PATTriggerProducer::autoProcessName_
private

Definition at line 72 of file PATTriggerProducer.h.

Referenced by beginRun(), and PATTriggerProducer().

bool pat::PATTriggerProducer::autoProcessNameL1ExtraCenJet_
private

Definition at line 99 of file PATTriggerProducer.h.

Referenced by beginRun(), and PATTriggerProducer().

bool pat::PATTriggerProducer::autoProcessNameL1ExtraETM_
private

Definition at line 102 of file PATTriggerProducer.h.

Referenced by beginRun(), and PATTriggerProducer().

bool pat::PATTriggerProducer::autoProcessNameL1ExtraForJet_
private

Definition at line 100 of file PATTriggerProducer.h.

Referenced by beginRun(), and PATTriggerProducer().

bool pat::PATTriggerProducer::autoProcessNameL1ExtraHTM_
private

Definition at line 103 of file PATTriggerProducer.h.

Referenced by beginRun(), and PATTriggerProducer().

bool pat::PATTriggerProducer::autoProcessNameL1ExtraIsoEG_
private

Definition at line 98 of file PATTriggerProducer.h.

Referenced by beginRun(), and PATTriggerProducer().

bool pat::PATTriggerProducer::autoProcessNameL1ExtraMu_
private

Definition at line 96 of file PATTriggerProducer.h.

Referenced by beginRun(), and PATTriggerProducer().

bool pat::PATTriggerProducer::autoProcessNameL1ExtraNoIsoEG_
private

Definition at line 97 of file PATTriggerProducer.h.

Referenced by beginRun(), and PATTriggerProducer().

bool pat::PATTriggerProducer::autoProcessNameL1ExtraTauJet_
private

Definition at line 101 of file PATTriggerProducer.h.

Referenced by beginRun(), and PATTriggerProducer().

std::vector< std::string > pat::PATTriggerProducer::exludeCollections_
private

Definition at line 121 of file PATTriggerProducer.h.

Referenced by PATTriggerProducer(), and produce().

bool pat::PATTriggerProducer::firstInRun_
private

Definition at line 74 of file PATTriggerProducer.h.

Referenced by beginRun(), and produce().

bool pat::PATTriggerProducer::hltConfigInit_
private

Definition at line 108 of file PATTriggerProducer.h.

Referenced by beginLuminosityBlock(), beginRun(), and produce().

std::string pat::PATTriggerProducer::hltPrescaleLabel_
private

Definition at line 113 of file PATTriggerProducer.h.

Referenced by PATTriggerProducer(), and produce().

HLTPrescaleProvider pat::PATTriggerProducer::hltPrescaleProvider_
private

Definition at line 107 of file PATTriggerProducer.h.

Referenced by beginRun(), and produce().

edm::GetterOfProducts< trigger::HLTPrescaleTable > pat::PATTriggerProducer::hltPrescaleTableEventGetter_
private

Definition at line 117 of file PATTriggerProducer.h.

Referenced by PATTriggerProducer().

trigger::HLTPrescaleTable pat::PATTriggerProducer::hltPrescaleTableLumi_
private

Definition at line 119 of file PATTriggerProducer.h.

Referenced by beginLuminosityBlock(), and produce().

edm::GetterOfProducts< trigger::HLTPrescaleTable > pat::PATTriggerProducer::hltPrescaleTableLumiGetter_
private

Definition at line 116 of file PATTriggerProducer.h.

Referenced by PATTriggerProducer().

trigger::HLTPrescaleTable pat::PATTriggerProducer::hltPrescaleTableRun_
private

Definition at line 118 of file PATTriggerProducer.h.

Referenced by beginLuminosityBlock(), and beginRun().

edm::GetterOfProducts< trigger::HLTPrescaleTable > pat::PATTriggerProducer::hltPrescaleTableRunGetter_
private

Definition at line 115 of file PATTriggerProducer.h.

Referenced by PATTriggerProducer().

edm::GetterOfProducts< l1extra::L1JetParticleCollection > pat::PATTriggerProducer::l1ExtraCenJetGetter_
private

Definition at line 87 of file PATTriggerProducer.h.

Referenced by PATTriggerProducer().

edm::GetterOfProducts< l1extra::L1EtMissParticleCollection > pat::PATTriggerProducer::l1ExtraETMGetter_
private

Definition at line 93 of file PATTriggerProducer.h.

Referenced by PATTriggerProducer().

edm::GetterOfProducts< l1extra::L1JetParticleCollection > pat::PATTriggerProducer::l1ExtraForJetGetter_
private

Definition at line 89 of file PATTriggerProducer.h.

Referenced by PATTriggerProducer().

edm::GetterOfProducts< l1extra::L1EtMissParticleCollection > pat::PATTriggerProducer::l1ExtraHTMGetter_
private

Definition at line 95 of file PATTriggerProducer.h.

Referenced by PATTriggerProducer().

edm::GetterOfProducts< l1extra::L1EmParticleCollection > pat::PATTriggerProducer::l1ExtraIsoEGGetter_
private

Definition at line 85 of file PATTriggerProducer.h.

Referenced by PATTriggerProducer().

edm::GetterOfProducts< l1extra::L1MuonParticleCollection > pat::PATTriggerProducer::l1ExtraMuGetter_
private

Definition at line 81 of file PATTriggerProducer.h.

Referenced by PATTriggerProducer().

edm::GetterOfProducts< l1extra::L1EmParticleCollection > pat::PATTriggerProducer::l1ExtraNoIsoEGGetter_
private

Definition at line 83 of file PATTriggerProducer.h.

Referenced by PATTriggerProducer().

edm::GetterOfProducts< l1extra::L1JetParticleCollection > pat::PATTriggerProducer::l1ExtraTauJetGetter_
private

Definition at line 91 of file PATTriggerProducer.h.

Referenced by PATTriggerProducer().

edm::EDGetTokenT< L1GlobalTriggerObjectMaps > pat::PATTriggerProducer::l1GlobalTriggerObjectMapsToken_
private

Definition at line 79 of file PATTriggerProducer.h.

Referenced by PATTriggerProducer(), and produce().

edm::ParameterSet* pat::PATTriggerProducer::l1PSet_
private

Definition at line 76 of file PATTriggerProducer.h.

Referenced by beginRun(), and produce().

std::string pat::PATTriggerProducer::labelHltPrescaleTable_
private

Definition at line 114 of file PATTriggerProducer.h.

Referenced by beginLuminosityBlock(), beginRun(), PATTriggerProducer(), and produce().

bool pat::PATTriggerProducer::mainBxOnly_
private

Definition at line 104 of file PATTriggerProducer.h.

Referenced by PATTriggerProducer(), and produce().

ModuleLabelToPathAndFlags pat::PATTriggerProducer::moduleLabelToPathAndFlags_
private

Definition at line 148 of file PATTriggerProducer.h.

Referenced by beginRun(), and produce().

std::string pat::PATTriggerProducer::nameProcess_
private

Definition at line 71 of file PATTriggerProducer.h.

Referenced by beginLuminosityBlock(), beginRun(), PATTriggerProducer(), and produce().

bool pat::PATTriggerProducer::onlyStandAlone_
private

Definition at line 73 of file PATTriggerProducer.h.

Referenced by PATTriggerProducer(), and produce().

bool pat::PATTriggerProducer::packPathNames_
private

Definition at line 122 of file PATTriggerProducer.h.

Referenced by produce().

bool pat::PATTriggerProducer::packPrescales_
private

Definition at line 123 of file PATTriggerProducer.h.

Referenced by PATTriggerProducer(), and produce().

bool pat::PATTriggerProducer::saveL1Refs_
private

Definition at line 105 of file PATTriggerProducer.h.

Referenced by PATTriggerProducer(), and produce().

edm::InputTag pat::PATTriggerProducer::tagL1ExtraCenJet_
private

Definition at line 86 of file PATTriggerProducer.h.

Referenced by beginRun(), PATTriggerProducer(), and produce().

edm::InputTag pat::PATTriggerProducer::tagL1ExtraETM_
private

Definition at line 92 of file PATTriggerProducer.h.

Referenced by beginRun(), PATTriggerProducer(), and produce().

edm::InputTag pat::PATTriggerProducer::tagL1ExtraForJet_
private

Definition at line 88 of file PATTriggerProducer.h.

Referenced by beginRun(), PATTriggerProducer(), and produce().

edm::InputTag pat::PATTriggerProducer::tagL1ExtraHTM_
private

Definition at line 94 of file PATTriggerProducer.h.

Referenced by beginRun(), PATTriggerProducer(), and produce().

edm::InputTag pat::PATTriggerProducer::tagL1ExtraIsoEG_
private

Definition at line 84 of file PATTriggerProducer.h.

Referenced by beginRun(), PATTriggerProducer(), and produce().

edm::InputTag pat::PATTriggerProducer::tagL1ExtraMu_
private

Definition at line 80 of file PATTriggerProducer.h.

Referenced by beginRun(), PATTriggerProducer(), and produce().

edm::InputTag pat::PATTriggerProducer::tagL1ExtraNoIsoEG_
private

Definition at line 82 of file PATTriggerProducer.h.

Referenced by beginRun(), PATTriggerProducer(), and produce().

edm::InputTag pat::PATTriggerProducer::tagL1ExtraTauJet_
private

Definition at line 90 of file PATTriggerProducer.h.

Referenced by beginRun(), PATTriggerProducer(), and produce().

edm::InputTag pat::PATTriggerProducer::tagL1GlobalTriggerObjectMaps_
private

Definition at line 78 of file PATTriggerProducer.h.

Referenced by PATTriggerProducer(), and produce().

edm::InputTag pat::PATTriggerProducer::tagTriggerEvent_
private

Definition at line 111 of file PATTriggerProducer.h.

Referenced by beginRun(), PATTriggerProducer(), and produce().

edm::InputTag pat::PATTriggerProducer::tagTriggerResults_
private

Definition at line 109 of file PATTriggerProducer.h.

Referenced by beginRun(), PATTriggerProducer(), and produce().

edm::GetterOfProducts< trigger::TriggerEvent > pat::PATTriggerProducer::triggerEventGetter_
private

Definition at line 112 of file PATTriggerProducer.h.

Referenced by PATTriggerProducer().

edm::GetterOfProducts< edm::TriggerResults > pat::PATTriggerProducer::triggerResultsGetter_
private

Definition at line 110 of file PATTriggerProducer.h.

Referenced by PATTriggerProducer().