HLTJetCollForElePlusJets< T > Class Template Reference

#include <HLTJetCollForElePlusJets.h>

Inheritance diagram for HLTJetCollForElePlusJets< T >:

Public Member Functions
	HLTJetCollForElePlusJets (const edm::ParameterSet &)
	~HLTJetCollForElePlusJets ()
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
virtual	~EDProducer ()
Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
boost::function< void(const BranchDescription &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()

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
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
typedef WorkerT< EDProducer >	WorkerType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Protected Member Functions inherited from edm::EDProducer
CurrentProcessingContext const *	currentContext () const
Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

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.

Date:

2012/02/12 11:41:59

Revision:

1.2

Author

Lukasz Kreczko

Definition at line 40 of file HLTJetCollForElePlusJets.h.

Constructor & Destructor Documentation

template<typename T >

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

explicit

Definition at line 22 of file HLTJetCollForElePlusJets.cc.

                                                                                   :
  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;
  produces<TCollection>();
}

template<typename T >

HLTJetCollForElePlusJets< T >::~HLTJetCollForElePlusJets

(

)

Definition at line 39 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 47 of file HLTJetCollForElePlusJets.cc.

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

                                                                                             {
    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 69 of file HLTJetCollForElePlusJets.cc.

References abs, edm::Event::getByLabel(), i, j, gen::k, position, edm::Event::put(), python.multivaluedict::sort(), trigger::TriggerCluster, trigger::TriggerElectron, vdt::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.getByLabel(hltElectronTag,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.getByLabel(sourceJetTag, 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;
  
}

Member Data Documentation

template<typename T >

edm::InputTag HLTJetCollForElePlusJets< T >::hltElectronTag

private

Definition at line 49 of file HLTJetCollForElePlusJets.h.

template<typename T >

double HLTJetCollForElePlusJets< T >::maxAbsJetEta_

private

Definition at line 53 of file HLTJetCollForElePlusJets.h.

template<typename T >

double HLTJetCollForElePlusJets< T >::minDeltaEta_

private

Definition at line 59 of file HLTJetCollForElePlusJets.h.

template<typename T >

double HLTJetCollForElePlusJets< T >::minDeltaR_

private

Definition at line 56 of file HLTJetCollForElePlusJets.h.

template<typename T >

double HLTJetCollForElePlusJets< T >::minJetPt_

private

Definition at line 52 of file HLTJetCollForElePlusJets.h.

template<typename T >

unsigned int HLTJetCollForElePlusJets< T >::minNJets_

private

Definition at line 54 of file HLTJetCollForElePlusJets.h.

template<typename T >

double HLTJetCollForElePlusJets< T >::minSoftJetPt_

private

Definition at line 58 of file HLTJetCollForElePlusJets.h.

template<typename T >

edm::InputTag HLTJetCollForElePlusJets< T >::sourceJetTag

private

Definition at line 50 of file HLTJetCollForElePlusJets.h.