#include <HLTPMMassFilter.h>

Inheritance diagram for HLTPMMassFilter:

Public Member Functions
virtual bool	hltFilter (edm::Event &, const edm::EventSetup &, trigger::TriggerFilterObjectWithRefs &filterproduct) const override

	HLTPMMassFilter (const edm::ParameterSet &)

	~HLTPMMassFilter ()

Public Member Functions inherited from HLTFilter
	HLTFilter (const edm::ParameterSet &config)

int	module (edm::Event const &) const

const std::string *	moduleLabel () const

int	path (edm::Event const &) const

const std::string *	pathName (edm::Event const &) const

std::pair< int, int >	pmid (edm::Event const &) const

bool	saveTags () const

virtual	~HLTFilter ()

Public Member Functions inherited from edm::global::EDFilter<>
	EDFilter ()=default

Public Member Functions inherited from edm::global::EDFilterBase
	EDFilterBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDFilterBase ()

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< ConsumesInfo >	consumesInfo () 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 (const std::string &iProcessName, std::vector< const char * > &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 ()

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

Static Public Member Functions inherited from HLTFilter
static void	makeHLTFilterDescription (edm::ParameterSetDescription &desc)

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
TLorentzVector	approxMomAtVtx (const MagneticField *magField, const GlobalPoint &xvert, const reco::SuperClusterRef sc, int charge) const

Private Attributes
edm::InputTag	beamSpot_

edm::EDGetTokenT< reco::BeamSpot >	beamSpotToken_

edm::InputTag	candTag_

edm::EDGetTokenT < trigger::TriggerFilterObjectWithRefs >	candToken_

bool	isElectron1_

bool	isElectron2_

edm::InputTag	L1IsoCollTag_

edm::InputTag	L1NonIsoCollTag_

double	lowerMassCut_

int	nZcandcut_

bool	relaxed_

bool	reqOppCharge_

double	upperMassCut_

Additional Inherited Members
Public Types inherited from edm::global::EDFilterBase
typedef EDFilterBase	ModuleType

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

Original Author: Jeremy Werner Institution: Princeton University, USA Contact: Jerem.nosp@m.y.We.nosp@m.rner@.nosp@m.cern.nosp@m..ch Date: February 21, 2007

Definition at line 55 of file HLTPMMassFilter.h.

Constructor & Destructor Documentation

HLTPMMassFilter::HLTPMMassFilter ( const edm::ParameterSet & iConfig )

explicit

Definition at line 16 of file HLTPMMassFilter.cc.

References beamSpot_, beamSpotToken_, candTag_, candToken_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), isElectron1_, isElectron2_, L1IsoCollTag_, L1NonIsoCollTag_, lowerMassCut_, nZcandcut_, relaxed_, reqOppCharge_, and upperMassCut_.

                                                                : HLTFilter(iConfig)
 {
   candTag_            = iConfig.getParameter< edm::InputTag > ("candTag");
   beamSpot_           = iConfig.getParameter< edm::InputTag > ("beamSpot");
   lowerMassCut_       = iConfig.getParameter<double> ("lowerMassCut");
   upperMassCut_       = iConfig.getParameter<double> ("upperMassCut");
   nZcandcut_          = iConfig.getParameter<int> ("nZcandcut");
   reqOppCharge_       = iConfig.getUntrackedParameter<bool> ("reqOppCharge",false);
   isElectron1_ = iConfig.getUntrackedParameter<bool> ("isElectron1",true) ;
   isElectron2_ = iConfig.getUntrackedParameter<bool> ("isElectron2",true) ;
   relaxed_ = iConfig.getUntrackedParameter<bool> ("relaxed",true) ;
   L1IsoCollTag_= iConfig.getParameter< edm::InputTag > ("L1IsoCand");
   L1NonIsoCollTag_= iConfig.getParameter< edm::InputTag > ("L1NonIsoCand");
   candToken_ = consumes<trigger::TriggerFilterObjectWithRefs>(candTag_);
   beamSpotToken_ = consumes<reco::BeamSpot>(beamSpot_);
 }

HLTPMMassFilter::~HLTPMMassFilter ( )

Definition at line 33 of file HLTPMMassFilter.cc.

33 {}

Member Function Documentation

TLorentzVector HLTPMMassFilter::approxMomAtVtx	(	const MagneticField *	magField,
		const GlobalPoint &	xvert,
		const reco::SuperClusterRef	sc,
		int	charge
	)		const

private

Definition at line 186 of file HLTPMMassFilter.cc.

References relval_parameters_module::energy, FTSFromVertexToPointFactory::get(), FreeTrajectoryState::momentum(), mathSSE::sqrt(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by hltFilter().

 {
     GlobalPoint xsc(sc->position().x(),
             sc->position().y(),
             sc->position().z());
     float energy = sc->energy();
     FreeTrajectoryState theFTS = FTSFromVertexToPointFactory::get(*magField, xsc, xvert, energy, charge);
     float theApproxMomMod = theFTS.momentum().x()*theFTS.momentum().x() +
                             theFTS.momentum().y()*theFTS.momentum().y() +
                             theFTS.momentum().z()*theFTS.momentum().z();
     TLorentzVector theApproxMom(theFTS.momentum().x(),
                 theFTS.momentum().y(),
                 theFTS.momentum().z(),
                                 sqrt(theApproxMomMod + 2.61121E-7));
     return theApproxMom ;
 }

void HLTPMMassFilter::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 36 of file HLTPMMassFilter.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addUntracked(), HLT_25ns14e33_v1_cff::InputTag, and HLTFilter::makeHLTFilterDescription().

                                                                             {
   edm::ParameterSetDescription desc;
   makeHLTFilterDescription(desc);
   desc.add<edm::InputTag>("candTag",edm::InputTag("hltL1NonIsoDoublePhotonEt5UpsHcalIsolFilter"));
   desc.add<edm::InputTag>("beamSpot",edm::InputTag("hltOfflineBeamSpot"));
   desc.add<double>("lowerMassCut",8.0);
   desc.add<double>("upperMassCut",11.0);
   desc.add<int>("nZcandcut",1);
   desc.addUntracked<bool>("reqOppCharge",true);
   desc.addUntracked<bool>("isElectron1",false);
   desc.addUntracked<bool>("isElectron2",false);
   desc.addUntracked<bool>("relaxed",true);
   desc.add<edm::InputTag>("L1IsoCand",edm::InputTag("hltL1IsoRecoEcalCandidate"));
   desc.add<edm::InputTag>("L1NonIsoCand",edm::InputTag("hltL1IsoRecoEcalCandidate"));
   descriptions.add("hltPMMassFilter",desc);
 }

bool HLTPMMassFilter::hltFilter	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		trigger::TriggerFilterObjectWithRefs &	filterproduct
	)		const

overridevirtual

Implements HLTFilter.

Definition at line 55 of file HLTPMMassFilter.cc.

References accept(), trigger::TriggerFilterObjectWithRefs::addCollectionTag(), trigger::TriggerRefsCollections::addObject(), approxMomAtVtx(), beamSpotToken_, candToken_, RecoTauCleanerPlugins::charge, HI_PhotonSkim_cff::electrons, edm::EventSetup::get(), edm::Event::getByToken(), i, cuy::ii, isElectron1_, isElectron2_, findQualityFiles::jj, L1IsoCollTag_, L1NonIsoCollTag_, lowerMassCut_, gen::n, nZcandcut_, dt_dqm_sourceclient_common_cff::reco, relaxed_, reqOppCharge_, HLTFilter::saveTags(), trigger::TriggerCluster, trigger::TriggerElectron, trigger::TriggerPhoton, and upperMassCut_.

 {
   using namespace std;
   using namespace edm;
   using namespace reco;
   using namespace trigger;
 
   // The filter object
   if (saveTags()) {
     filterproduct.addCollectionTag(L1IsoCollTag_);
     if (relaxed_) filterproduct.addCollectionTag(L1NonIsoCollTag_);
   }
 
   edm::ESHandle<MagneticField> theMagField;
   iSetup.get<IdealMagneticFieldRecord>().get(theMagField);
 
   edm::Handle<trigger::TriggerFilterObjectWithRefs> PrevFilterOutput;
   iEvent.getByToken (candToken_,PrevFilterOutput);
 
   // beam spot
   edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
   iEvent.getByToken(beamSpotToken_,recoBeamSpotHandle);
   // gets its position
   const GlobalPoint vertexPos(recoBeamSpotHandle->position().x(),
                   recoBeamSpotHandle->position().y(),
                   recoBeamSpotHandle->position().z());
 
   int n = 0;
 
   // REMOVED USAGE OF STATIC ARRAYS
   // double px[66];
   // double py[66];
   // double pz[66];
   // double energy[66];
   std::vector<TLorentzVector> pEleCh1;
   std::vector<TLorentzVector> pEleCh2;
   std::vector<double> charge;
   std::vector<double> etaOrig;
 
   if (isElectron1_ && isElectron2_) {
 
     Ref< ElectronCollection > refele;
 
     vector< Ref< ElectronCollection > > electrons;
     PrevFilterOutput->getObjects(TriggerElectron, electrons);
 
     for (unsigned int i=0; i<electrons.size(); i++) {
 
       refele = electrons[i];
 
       TLorentzVector pThisEle(refele->px(), refele->py(),
                   refele->pz(), refele->energy() );
       pEleCh1.push_back( pThisEle );
       charge.push_back( refele->charge() );
     }
 
     for(unsigned int jj=0;jj<electrons.size();jj++){
 
       TLorentzVector p1Ele = pEleCh1.at(jj);
       for(unsigned int ii=jj+1;ii<electrons.size();ii++){
     
     TLorentzVector p2Ele = pEleCh1.at(ii);
     
     if(fabs(p1Ele.E() - p2Ele.E()) < 0.00001) continue;
     if(reqOppCharge_ && charge[jj]*charge[ii] > 0) continue;
     
     TLorentzVector pTot = p1Ele + p2Ele;
     double mass = pTot.M();
     
     if(mass>=lowerMassCut_ && mass<=upperMassCut_){
       n++;
       refele = electrons[ii];
       filterproduct.addObject(TriggerElectron, refele);
       refele = electrons[jj];
       filterproduct.addObject(TriggerElectron, refele);
     }
       }
     }
 
   } else {
 
     Ref< RecoEcalCandidateCollection > refsc;
 
     vector< Ref< RecoEcalCandidateCollection > > scs;
     PrevFilterOutput->getObjects(TriggerCluster, scs);
     if(scs.empty()) PrevFilterOutput->getObjects(TriggerPhoton, scs);  //we dont know if its type trigger cluster or trigger photon
 
     for (unsigned int i=0; i<scs.size(); i++) {
 
       refsc = scs[i];
       const reco::SuperClusterRef sc = refsc->superCluster();
       TLorentzVector pscPos = approxMomAtVtx(theMagField.product(), vertexPos, sc, 1);
       pEleCh1.push_back( pscPos );
 
       TLorentzVector pscEle = approxMomAtVtx(theMagField.product(), vertexPos, sc, -1);
       pEleCh2.push_back( pscEle );
       etaOrig.push_back( sc->eta() );
 
     }
 
     for(unsigned int jj=0;jj<scs.size();jj++){
 
       TLorentzVector p1Ele = pEleCh1.at(jj);
       for(unsigned int ii=0;ii<scs.size();ii++){
     
     TLorentzVector p2Ele = pEleCh2.at(ii);
     
     if(fabs(p1Ele.E() - p2Ele.E()) < 0.00001) continue;
     
     TLorentzVector pTot = p1Ele + p2Ele;
     double mass = pTot.M();
     
     if(mass>= lowerMassCut_ && mass<=upperMassCut_){
       n++;
       refsc = scs[ii];
       filterproduct.addObject(TriggerCluster, refsc);
       refsc = scs[jj];
       filterproduct.addObject(TriggerCluster, refsc);
     }
       }
     }
   }
 
 
   // filter decision
   bool accept(n>=nZcandcut_);
   // if (accept) std::cout << "n size = " << n << std::endl;
 
   return accept;
 }