#include <HLTJetCollForElePlusJets.h>

Inheritance diagram for HLTJetCollForElePlusJets< T >:

Public Member Functions
	HLTJetCollForElePlusJets (const edm::ParameterSet &)

	~HLTJetCollForElePlusJets ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	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
	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 (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
virtual void	produce (edm::Event &, const edm::EventSetup &)

Private Attributes
edm::InputTag	hltElectronTag

edm::EDGetTokenT < trigger::TriggerFilterObjectWithRefs >	m_theElectronToken

edm::EDGetTokenT< std::vector < T > >	m_theJetToken

double	maxAbsJetEta_

double	minDeltaEta_

double	minDeltaR_

double	minJetPt_

unsigned int	minNJets_

double	minSoftJetPt_

edm::InputTag	sourceJetTag

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

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

template<typename T>
class HLTJetCollForElePlusJets< T >

This class is an EDProducer implementing an HLT trigger for electron and jet objects, cutting on variables relating to the jet 4-momentum representation. The producer checks for overlaps between electrons and jets and if a combination of one electron + jets cleaned against this electrons satisfy the cuts. These jets are then added to a cleaned jet collection which is put into the event.

Author: Lukasz Kreczko

Definition at line 41 of file HLTJetCollForElePlusJets.h.

Constructor & Destructor Documentation

template<typename T >

HLTJetCollForElePlusJets< T >::HLTJetCollForElePlusJets ( const edm::ParameterSet & iConfig )

explicit

Definition at line 16 of file HLTJetCollForElePlusJets.cc.

References HLTJetCollForElePlusJets< T >::hltElectronTag, HLTJetCollForElePlusJets< T >::m_theElectronToken, HLTJetCollForElePlusJets< T >::m_theJetToken, and HLTJetCollForElePlusJets< T >::sourceJetTag.

                                                                                    :
   hltElectronTag(iConfig.getParameter< edm::InputTag > ("HltElectronTag")),
   sourceJetTag(iConfig.getParameter< edm::InputTag > ("SourceJetTag")),
   minJetPt_(iConfig.getParameter<double> ("MinJetPt")),
   maxAbsJetEta_(iConfig.getParameter<double> ("MaxAbsJetEta")),
   minNJets_(iConfig.getParameter<unsigned int> ("MinNJets")),
   minDeltaR_(iConfig.getParameter< double > ("minDeltaR")),
   //Only for VBF
   minSoftJetPt_(iConfig.getParameter< double > ("MinSoftJetPt")),
   minDeltaEta_(iConfig.getParameter< double > ("MinDeltaEta"))
 {
   typedef std::vector<T> TCollection;
   m_theElectronToken = consumes<trigger::TriggerFilterObjectWithRefs>(hltElectronTag);
   m_theJetToken = consumes<TCollection>(sourceJetTag);
   produces<TCollection>();
 }

template<typename T >

HLTJetCollForElePlusJets< T >::~HLTJetCollForElePlusJets ( )

Definition at line 35 of file HLTJetCollForElePlusJets.cc.

 {
    // do anything here that needs to be done at desctruction time
    // (e.g. close files, deallocate resources etc.)
 
 }

Member Function Documentation

template<typename T >

void HLTJetCollForElePlusJets< T >::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 43 of file HLTJetCollForElePlusJets.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), mergeVDriftHistosByStation::name, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                              {
     edm::ParameterSetDescription desc;
     desc.add<edm::InputTag> ("HltElectronTag", edm::InputTag("triggerFilterObjectWithRefs"));
     desc.add<edm::InputTag> ("SourceJetTag", edm::InputTag("jetCollection"));
     desc.add<double> ("MinJetPt", 30.);
     desc.add<double> ("MaxAbsJetEta", 2.6);
     desc.add<unsigned int> ("MinNJets", 1);
     desc.add<double> ("minDeltaR", 0.5);
     //Only for VBF
     desc.add<double> ("MinSoftJetPt", 25.);
     desc.add<double> ("MinDeltaEta", -1.);    
     descriptions.add(std::string("hlt")+std::string(typeid(HLTJetCollForElePlusJets<T>).name()), desc);
 }

template<typename T >

void HLTJetCollForElePlusJets< T >::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Implements edm::EDProducer.

Definition at line 65 of file HLTJetCollForElePlusJets.cc.

References funct::abs(), edm::Event::getByToken(), i, j, gen::k, position, edm::Event::put(), python.multivaluedict::sort(), trigger::TriggerCluster, trigger::TriggerElectron, x, detailsBasic3DVector::y, and detailsBasic3DVector::z.

 {
   using namespace edm;
   using namespace std;
 
   typedef vector<T> TCollection;
   typedef Ref<TCollection> TRef;
   typedef edm::RefVector<TCollection> TRefVector;
   typedef std::vector<edm::RefVector<std::vector<T>,T,edm::refhelper::FindUsingAdvance<std::vector<T>,T> > > TCollectionVector;
 
   edm::Handle<trigger::TriggerFilterObjectWithRefs> PrevFilterOutput;
   iEvent.getByToken(m_theElectronToken,PrevFilterOutput);
  
   //its easier on the if statement flow if I try everything at once, shouldnt add to timing
   std::vector<edm::Ref<reco::RecoEcalCandidateCollection> > clusCands;
   PrevFilterOutput->getObjects(trigger::TriggerCluster,clusCands);
 
   std::vector<edm::Ref<reco::ElectronCollection> > eleCands;
   PrevFilterOutput->getObjects(trigger::TriggerElectron,eleCands);
   
   //prepare the collection of 3-D vector for electron momenta
   std::vector<TVector3> ElePs;
 
   if(!clusCands.empty()){ //try trigger cluster
     for(size_t candNr=0;candNr<clusCands.size();candNr++){
       TVector3 positionVector(
                   clusCands[candNr]->superCluster()->position().x(),
                   clusCands[candNr]->superCluster()->position().y(),
                   clusCands[candNr]->superCluster()->position().z());
       ElePs.push_back(positionVector);
     }
   }else if(!eleCands.empty()){ // try trigger electrons
     for(size_t candNr=0;candNr<eleCands.size();candNr++){
       TVector3 positionVector(
                   eleCands[candNr]->superCluster()->position().x(),
                   eleCands[candNr]->superCluster()->position().y(),
                   eleCands[candNr]->superCluster()->position().z());
       ElePs.push_back(positionVector);
     }
   }
   
   edm::Handle<TCollection> theJetCollectionHandle;
   iEvent.getByToken(m_theJetToken, theJetCollectionHandle);
  
   const TCollection & theJetCollection = *theJetCollectionHandle;
   
   std::auto_ptr< TCollection >  theFilteredJetCollection(new TCollection);
   
   std::auto_ptr < TCollectionVector > allSelections(new TCollectionVector());
   
   bool foundSolution(false);
 
   for (unsigned int i = 0; i < ElePs.size(); i++) {
     
     bool VBFJetPair = false;
     std::vector<int> store_jet;
     TRefVector refVector;
     
     for (unsigned int j = 0; j < theJetCollection.size(); j++) {
       TVector3 JetP(theJetCollection[j].px(), theJetCollection[j].py(),
                     theJetCollection[j].pz());
       double DR = ElePs[i].DeltaR(JetP);
       
       if (JetP.Pt() > minJetPt_ && std::abs(JetP.Eta()) < maxAbsJetEta_ && DR > minDeltaR_) {
     store_jet.push_back(j);
     // The VBF part of the filter
     if ( minDeltaEta_ > 0 ) {
       for ( unsigned int k = j+1; k < theJetCollection.size(); k++ ) {
         TVector3 SoftJetP(theJetCollection[k].px(), theJetCollection[k].py(),
                   theJetCollection[k].pz());
         double softDR = ElePs[i].DeltaR(SoftJetP);
         
         if (SoftJetP.Pt() > minSoftJetPt_ && std::abs(SoftJetP.Eta()) < maxAbsJetEta_ && softDR > minDeltaR_)
           if ( std::abs(SoftJetP.Eta() - JetP.Eta()) > minDeltaEta_ ) {
         store_jet.push_back(k);
         VBFJetPair = true;
           }
       }
     }
       }
       
     }
     
     // Now remove duplicates from the jet collection to store
     std::sort( store_jet.begin(), store_jet.end() );
     store_jet.erase( unique( store_jet.begin(), store_jet.end() ), store_jet.end() );
     
     // Now save the cleaned jets
     for ( unsigned int ijet = 0; ijet < store_jet.size(); ijet++ )
       {
     //store all selections
     refVector.push_back(TRef(theJetCollectionHandle, store_jet.at(ijet)));
     //store first selection which matches the criteria
     if(!foundSolution)
       theFilteredJetCollection->push_back(theJetCollection[store_jet.at(ijet)]);
       }
     //store all selections
     allSelections->push_back(refVector);
     
     if (theFilteredJetCollection->size() >= minNJets_ && minDeltaEta_ < 0)
       foundSolution = true;
     else if (VBFJetPair && minDeltaEta_ > 0)
       foundSolution = true;
     else if (!foundSolution)
       theFilteredJetCollection->clear();
     
     
   }
   
   iEvent.put(theFilteredJetCollection);
   
   return;
   
 }