#include <SUSY_HLT_Razor.h>

Inheritance diagram for SUSY_HLT_Razor:

Public Member Functions
	SUSY_HLT_Razor (const edm::ParameterSet &ps)

	~SUSY_HLT_Razor () override

Public Member Functions inherited from DQMEDAnalyzer
void	beginRun (edm::Run const &, edm::EventSetup const &) final

void	beginStream (edm::StreamID id) final

	DQMEDAnalyzer ()

void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

uint32_t	streamId () const

Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default

Public Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static double	CalcMR (TLorentzVector ja, TLorentzVector jb)

static double	CalcR (double MR, TLorentzVector ja, TLorentzVector jb, edm::Handle< edm::View< reco::MET > > met, const std::vector< math::XYZTLorentzVector > &muons)

static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr< dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)

static std::shared_ptr< dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)

static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const , dqmDetails::NoCache )

static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const , dqmDetails::NoCache )

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions
void	analyze (edm::Event const &e, edm::EventSetup const &eSetup) override

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &eSetup) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

void	dqmBeginRun (edm::Run const &, edm::EventSetup const &) override

void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &eSetup) override

void	endRun (edm::Run const &run, edm::EventSetup const &eSetup) override

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)

Private Member Functions
void	bookHistos (DQMStore::IBooker &)

Private Attributes
edm::InputTag	caloFilter_

MonitorElement *	h_calo_mr_vs_mr

MonitorElement *	h_calo_rsq_vs_rsq

MonitorElement *	h_ht

MonitorElement *	h_htMet

MonitorElement *	h_met

MonitorElement *	h_mr

MonitorElement *	h_mr_tight

MonitorElement *	h_mrRsq

MonitorElement *	h_online_mr_vs_mr

MonitorElement *	h_online_rsq_vs_rsq

MonitorElement *	h_rsq

MonitorElement *	h_rsq_loose

MonitorElement *	h_rsq_tight

edm::EDGetTokenT< std::vector< math::XYZTLorentzVector > >	theHemispheres_

edm::EDGetTokenT< edm::View< reco::Jet > >	theJetCollection_

edm::EDGetTokenT< edm::View< reco::MET > >	theMETCollection_

edm::EDGetTokenT< trigger::TriggerEvent >	theTrigSummary_

edm::InputTag	triggerFilter_

std::string	triggerPath_

edm::EDGetTokenT< edm::TriggerResults >	triggerResults_

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
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::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Detailed Description

Definition at line 37 of file SUSY_HLT_Razor.h.

Constructor & Destructor Documentation

SUSY_HLT_Razor::SUSY_HLT_Razor ( const edm::ParameterSet & ps )

Definition at line 9 of file SUSY_HLT_Razor.cc.

References caloFilter_, edm::ParameterSet::getParameter(), AlCaHLTBitMon_QueryRunRegistry::string, theHemispheres_, theJetCollection_, theMETCollection_, theTrigSummary_, triggerFilter_, triggerPath_, and triggerResults_.

 {
   edm::LogInfo("SUSY_HLT_Razor") << "Constructor SUSY_HLT_Razor::SUSY_HLT_Razor " << std::endl;
   // Get parameters from configuration file
   theTrigSummary_ = consumes<trigger::TriggerEvent>(ps.getParameter<edm::InputTag>("trigSummary"));
   theMETCollection_ = consumes<edm::View<reco::MET> >(ps.getParameter<edm::InputTag>("METCollection"));
   theHemispheres_ = consumes<std::vector<math::XYZTLorentzVector> >(ps.getParameter<edm::InputTag>("hemispheres"));
   triggerResults_ = consumes<edm::TriggerResults>(ps.getParameter<edm::InputTag>("TriggerResults"));
   triggerPath_ = ps.getParameter<std::string>("TriggerPath");
   triggerFilter_ = ps.getParameter<edm::InputTag>("TriggerFilter");
   caloFilter_ = ps.getParameter<edm::InputTag>("CaloFilter");
   theJetCollection_ = consumes<edm::View<reco::Jet> >(ps.getParameter<edm::InputTag>("jetCollection"));
 }

SUSY_HLT_Razor::~SUSY_HLT_Razor ( )

override

Definition at line 23 of file SUSY_HLT_Razor.cc.

 {
    edm::LogInfo("SUSY_HLT_Razor") << "Destructor SUSY_HLT_Razor::~SUSY_HLT_Razor " << std::endl;
 }

Member Function Documentation

void SUSY_HLT_Razor::analyze	(	edm::Event const &	e,
		edm::EventSetup const &	eSetup
	)

overrideprotectedvirtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 59 of file SUSY_HLT_Razor.cc.

References funct::abs(), edm::HLTGlobalStatus::accept(), CalcMR(), CalcR(), caloFilter_, MonitorElement::Fill(), trigger::TriggerEvent::filterIndex(), trigger::TriggerEvent::filterKeys(), edm::Event::getByToken(), trigger::TriggerEvent::getObjects(), h_calo_mr_vs_mr, h_calo_rsq_vs_rsq, h_ht, h_htMet, h_met, h_mr, h_mr_tight, h_mrRsq, h_online_mr_vs_mr, h_online_rsq_vs_rsq, h_rsq, h_rsq_loose, h_rsq_tight, mps_fire::i, trigger::TriggerObject::id(), edm::HandleBase::isValid(), relativeConstraints::keys, EnergyCorrector::pt, trigger::TriggerObject::px(), trigger::TriggerObject::py(), dttmaxenums::R, edm::TriggerNames::size(), trigger::TriggerEvent::sizeFilters(), theHemispheres_, theJetCollection_, theMETCollection_, theTrigSummary_, triggerFilter_, edm::TriggerNames::triggerName(), edm::Event::triggerNames(), TriggerAnalyzer::triggerObjects, triggerPath_, triggerResults_, trigNames, and edm::HLTGlobalStatus::wasrun().

                                                                           {
     
   edm::LogInfo("SUSY_HLT_Razor") << "SUSY_HLT_Razor::analyze" << std::endl;
 
   using namespace std;
   using namespace edm;
   using namespace reco;
 
   // get hold of collection of objects
   Handle< vector<math::XYZTLorentzVector> > hemispheres;
   e.getByToken (theHemispheres_,hemispheres);
   // get hold of the MET Collection
   edm::Handle<edm::View<reco::MET> > inputMet;
   e.getByToken(theMETCollection_,inputMet);
   if ( !inputMet.isValid() ){
     edm::LogError ("SUSY_HLT_Razor") << "invalid collection: MET" << "\n";
    return;
   }
   //  edm::Handle<reco::PFJetCollection> jetCollection;
   edm::Handle<edm::View<reco::Jet> > jetCollection;
   e.getByToken (theJetCollection_,jetCollection);
   if ( !jetCollection.isValid() ){
     edm::LogError ("SUSY_HLT_Razor") << "invalid collection: jets" << "\n";
     return;
   }
 
   //check what is in the menu
   edm::Handle<edm::TriggerResults> hltresults;
   e.getByToken(triggerResults_,hltresults);
   if(!hltresults.isValid()){
     edm::LogError ("SUSY_HLT_Razor") << "invalid collection: TriggerResults" << "\n";
     return;
   }
   
   //-------------------------------
   //--- Trigger
   //-------------------------------
   edm::Handle<trigger::TriggerEvent> triggerSummary;
   e.getByToken(theTrigSummary_, triggerSummary);
   if(!triggerSummary.isValid()) {
     edm::LogError ("SUSY_HLT_Razor") << "invalid collection: TriggerSummary" << "\n";
     return;
   }
   //get online objects
   //HLTriggerOffline/Egamma/python/TriggerTypeDefs.py contains the trigger object IDs
   double onlineMR = 0, onlineRsq = 0;
   double caloMR = 0, caloRsq = 0; //online razor variables computed using calo quantities
   size_t filterIndex = triggerSummary->filterIndex( triggerFilter_ );
   size_t caloFilterIndex = triggerSummary->filterIndex( caloFilter_ );
   //search for online MR and Rsq objects
   trigger::TriggerObjectCollection triggerObjects = triggerSummary->getObjects();
   if( !(filterIndex >= triggerSummary->sizeFilters()) ){
       const trigger::Keys& keys = triggerSummary->filterKeys( filterIndex );
       for( size_t j = 0; j < keys.size(); ++j ){
           trigger::TriggerObject foundObject = triggerObjects[keys[j]];
           if(foundObject.id() == 0){ //the MET object containing MR and Rsq will show up with ID = 0
               onlineMR = foundObject.px(); //razor variables stored in dummy reco::MET objects
               onlineRsq = foundObject.py();
           }
       }
   }
 
   //search for calo MR and Rsq objects
   if( !(caloFilterIndex >= triggerSummary->sizeFilters()) ){
       const trigger::Keys& keys = triggerSummary->filterKeys( caloFilterIndex );
       for( size_t j = 0; j < keys.size(); ++j ){
           trigger::TriggerObject foundObject = triggerObjects[keys[j]];
           if(foundObject.id() == 0){ 
               caloMR = foundObject.px(); //razor variables stored in dummy reco::MET objects
               caloRsq = foundObject.py();
           }
       }
   }
 
   bool hasFired = false;
   const edm::TriggerNames& trigNames = e.triggerNames(*hltresults);
   unsigned int numTriggers = trigNames.size();
   for( unsigned int hltIndex=0; hltIndex<numTriggers; ++hltIndex ){
       if (trigNames.triggerName(hltIndex).find(triggerPath_) != std::string::npos && hltresults->wasrun(hltIndex) && hltresults->accept(hltIndex)) 
     {
       hasFired = true;
     }      
   }
 
   float HT = 0.0;
   
   for (unsigned int i=0; i<jetCollection->size(); i++) {    
     if(std::abs(jetCollection->at(i).eta()) < 3.0 && jetCollection->at(i).pt() >= 40.0){
       HT += jetCollection->at(i).pt();
     }
   }
   //for (reco::PFJetCollection::const_iterator i_jet = jetCollection->begin(); i_jet != jetCollection->end(); ++i_jet){
   //    if (i_pfjet->pt() < 40) continue;
   //    if (fabs(i_pfjet->eta()) > 3.0) continue;
   //    pfHT += i_pfjet->pt();
   //}
   float MET = (inputMet->front()).pt();
 
   //this part is adapted from HLTRFilter.cc 
 
   // check that the input collections are available
   if (not hemispheres.isValid()){
       //This is happening many times (which is normal) and it's noisy for the output, we removed the error message
       //edm::LogError("SUSY_HLT_Razor") << "Hemisphere object is invalid!" << "\n";
       return;
   }
 
   if(hasFired){
 
       if(hemispheres->empty()){  // the Hemisphere Maker will produce an empty collection of hemispheres if the number of jets in the
           edm::LogError("SUSY_HLT_Razor") << "Cannot calculate M_R and R^2 because there are too many jets! (trigger passed automatically without forming the hemispheres)" << endl;
           return;
       }
 
       //***********************************
       //Calculate R 
 
       if(!hemispheres->empty() && hemispheres->size() != 2 && hemispheres->size() != 5 && hemispheres->size() != 10){
           edm::LogError("SUSY_HLT_Razor") << "Invalid hemisphere collection!  hemispheres->size() = " << hemispheres->size() << endl;
           return;
       }
 
       TLorentzVector ja(hemispheres->at(0).x(),hemispheres->at(0).y(),hemispheres->at(0).z(),hemispheres->at(0).t());
       TLorentzVector jb(hemispheres->at(1).x(),hemispheres->at(1).y(),hemispheres->at(1).z(),hemispheres->at(1).t());
 
       //dummy vector (this trigger does not care about muons)
       std::vector<math::XYZTLorentzVector> muonVec;
 
       double MR = CalcMR(ja,jb);
       double R  = CalcR(MR,ja,jb,inputMet,muonVec);
       double Rsq = R*R;    
 
       if(Rsq > 0.15) h_mr->Fill(MR);
       if(MR > 300) h_rsq->Fill(Rsq);
       h_mrRsq->Fill(MR, Rsq);
 
       h_rsq_loose->Fill(Rsq);
 
       if(Rsq > 0.25){
           h_mr_tight->Fill(MR);
       }
       if(MR > 400){
           h_rsq_tight->Fill(Rsq);
       }
 
       h_ht->Fill(HT);
       h_met->Fill(MET);
       h_htMet->Fill(HT, MET);
 
       h_online_mr_vs_mr->Fill(MR, onlineMR);
       h_online_rsq_vs_rsq->Fill(Rsq, onlineRsq);
 
       h_calo_mr_vs_mr->Fill(MR, caloMR);
       h_calo_rsq_vs_rsq->Fill(Rsq, caloRsq);
 
   } 
 
 }

void SUSY_HLT_Razor::beginLuminosityBlock	(	edm::LuminosityBlock const &	lumi,
		edm::EventSetup const &	eSetup
	)

overrideprotectedvirtual

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 40 of file SUSY_HLT_Razor.cc.

 {
    edm::LogInfo("SUSY_HLT_Razor") << "SUSY_HLT_Razor::beginLuminosityBlock" << std::endl;
 }

void SUSY_HLT_Razor::bookHistograms	(	DQMStore::IBooker &	ibooker_,
		edm::Run const &	,
		edm::EventSetup const &
	)

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 33 of file SUSY_HLT_Razor.cc.

References bookHistos().

 {
   edm::LogInfo("SUSY_HLT_Razor") << "SUSY_HLT_Razor::bookHistograms" << std::endl;
   //book at beginRun
   bookHistos(ibooker_);
 }

void SUSY_HLT_Razor::bookHistos ( DQMStore::IBooker & ibooker_ )

private

Definition at line 229 of file SUSY_HLT_Razor.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::cd(), h_calo_mr_vs_mr, h_calo_rsq_vs_rsq, h_ht, h_htMet, h_met, h_mr, h_mr_tight, h_mrRsq, h_online_mr_vs_mr, h_online_rsq_vs_rsq, h_rsq, h_rsq_loose, h_rsq_tight, DQMStore::IBooker::setCurrentFolder(), and triggerPath_.

Referenced by bookHistograms().

 {
   ibooker_.cd();
   ibooker_.setCurrentFolder("HLT/SUSYBSM/" + triggerPath_);
 
   h_mr = ibooker_.book1D("mr", "M_{R} (R^{2} > 0.15) ; GeV", 100, 0.0, 4000);
   h_rsq = ibooker_.book1D("rsq", "R^{2} (M_{R} > 300)", 100, 0.0, 1.5);
   h_mrRsq = ibooker_.book2D("mrRsq", "R^{2} vs M_{R}; GeV; ", 100, 0.0, 4000.0, 100, 0.0, 1.5);
 
   h_mr_tight = ibooker_.book1D("mr_tight", "M_{R} (R^{2} > 0.25) ; GeV", 100, 0.0, 4000);
   h_rsq_tight = ibooker_.book1D("rsq_tight", "R^{2} (M_{R} > 400) ; ", 100, 0.0, 1.5);
 
   h_rsq_loose = ibooker_.book1D("rsq_loose", "R^{2} (M_{R} > 0) ; ", 100, 0.0, 1.5);
 
   h_ht = ibooker_.book1D("ht", "HT; GeV; ", 100, 0.0, 4000.0);
   h_met = ibooker_.book1D("met", "MET; GeV; ", 100, 0.0, 1000);
   h_htMet = ibooker_.book2D("htMet", "MET vs HT; GeV; ", 100, 0.0, 4000.0, 100, 0.0, 1000);
 
   h_online_mr_vs_mr = ibooker_.book2D("online_mr_vs_mr", "Online M_{R} vs  Offline M_{R} (events passing trigger); Offline M_{R} (GeV); Online M_{R} (GeV); ", 100, 0.0, 4000.0, 100, 0.0, 4000.0);
   h_calo_mr_vs_mr = ibooker_.book2D("calo_mr_vs_mr", "Calo M_{R} vs  Offline M_{R} (events passing trigger); Offline M_{R} (GeV); Calo M_{R} (GeV); ", 100, 0.0, 4000.0, 100, 0.0, 4000.0);
   h_online_rsq_vs_rsq = ibooker_.book2D("online_rsq_vs_rsq", "Online R^{2} vs Offline R^{2} (events passing trigger); Offline R^{2}; Online R^{2}; ", 100, 0.0, 1.5, 100, 0.0, 1.5);
   h_calo_rsq_vs_rsq = ibooker_.book2D("calo_rsq_vs_rsq", "Calo R^{2} vs Offline R^{2} (events passing trigger); Offline R^{2}; Calo R^{2}; ", 100, 0.0, 1.5, 100, 0.0, 1.5);
 
   ibooker_.cd();
 }

double SUSY_HLT_Razor::CalcMR	(	TLorentzVector	ja,
		TLorentzVector	jb
	)

static

Definition at line 257 of file SUSY_HLT_Razor.cc.

References patCaloMETCorrections_cff::A, TtFullHadDaughter::B, mathSSE::sqrt(), and groupFilesInBlocks::temp.

Referenced by analyze().

                                                           {
   if(ja.Pt()<=0.1) return -1;
 
   ja.SetPtEtaPhiM(ja.Pt(),ja.Eta(),ja.Phi(),0.0);
   jb.SetPtEtaPhiM(jb.Pt(),jb.Eta(),jb.Phi(),0.0);
   
   if(ja.Pt() > jb.Pt()){
     TLorentzVector temp = ja;
     ja = jb;
     jb = temp;
   }
   
   double A = ja.P();
   double B = jb.P();
   double az = ja.Pz();
   double bz = jb.Pz();
   TVector3 jaT, jbT;
   jaT.SetXYZ(ja.Px(),ja.Py(),0.0);
   jbT.SetXYZ(jb.Px(),jb.Py(),0.0);
   double ATBT = (jaT+jbT).Mag2();
   
   double MR = sqrt((A+B)*(A+B)-(az+bz)*(az+bz)-
            (jbT.Dot(jbT)-jaT.Dot(jaT))*(jbT.Dot(jbT)-jaT.Dot(jaT))/(jaT+jbT).Mag2());
   
   double mybeta = (jbT.Dot(jbT)-jaT.Dot(jaT))/
     sqrt(ATBT*((A+B)*(A+B)-(az+bz)*(az+bz)));
   
   double mygamma = 1./sqrt(1.-mybeta*mybeta);
   
   //use gamma times MRstar
   return MR*mygamma;  
 }

double SUSY_HLT_Razor::CalcR	(	double	MR,
		TLorentzVector	ja,
		TLorentzVector	jb,
		edm::Handle< edm::View< reco::MET > >	met,
		const std::vector< math::XYZTLorentzVector > &	muons
	)

static

Definition at line 291 of file SUSY_HLT_Razor.cc.

References DEFINE_FWK_MODULE, objects.autophobj::float, RazorAnalyzer::met, phi, EnergyCorrector::pt, mathSSE::sqrt(), and tmp.

Referenced by analyze().

                                                                                                                                                             {
   //now we can calculate MTR
   TVector3 met;
   met.SetPtEtaPhi((inputMet->front()).pt(),0.0,(inputMet->front()).phi());
   
   std::vector<math::XYZTLorentzVector>::const_iterator muonIt;
   for(muonIt = muons.begin(); muonIt!=muons.end(); muonIt++){
     TVector3 tmp;
     tmp.SetPtEtaPhi(muonIt->pt(),0,muonIt->phi());
     met-=tmp;
   }
 
   double MTR = sqrt(0.5*(met.Mag()*(ja.Pt()+jb.Pt()) - met.Dot(ja.Vect()+jb.Vect())));
   
   //filter events
   return float(MTR)/float(MR); //R
   
 }

void SUSY_HLT_Razor::dqmBeginRun	(	edm::Run const &	,
		edm::EventSetup const &
	)

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 28 of file SUSY_HLT_Razor.cc.

 {
   edm::LogInfo("SUSY_HLT_Razor") << "SUSY_HLT_Razor::beginRun" << std::endl;
 }

void SUSY_HLT_Razor::endLuminosityBlock	(	edm::LuminosityBlock const &	lumi,
		edm::EventSetup const &	eSetup
	)

overrideprotectedvirtual

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 218 of file SUSY_HLT_Razor.cc.

 {
   edm::LogInfo("SUSY_HLT_Razor") << "SUSY_HLT_Razor::endLuminosityBlock" << std::endl;
 }

void SUSY_HLT_Razor::endRun	(	edm::Run const &	run,
		edm::EventSetup const &	eSetup
	)

overrideprotectedvirtual

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 224 of file SUSY_HLT_Razor.cc.

 {
   edm::LogInfo("SUSY_HLT_Razor") << "SUSY_HLT_Razor::endRun" << std::endl;
 }

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

static

Definition at line 46 of file SUSY_HLT_Razor.cc.

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

                                                                                 {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("jetCollection",edm::InputTag("ak4PFJetsCHS"));
   desc.add<edm::InputTag>("CaloFilter",edm::InputTag("hltRsqMR200Rsq0p01MR100Calo","","HLT"))->setComment("Calo HLT filter module used to save razor variable objects");
   desc.add<edm::InputTag>("METCollection",edm::InputTag("pfMet"));
   desc.add<edm::InputTag>("hemispheres",edm::InputTag("hemispheres"))->setComment("hemisphere jets used to compute razor variables");
   desc.add<std::string>("TriggerPath","HLT_RsqMR300_Rsq0p09_MR200_v")->setComment("trigger path name");
   desc.add<edm::InputTag>("TriggerFilter",edm::InputTag("hltRsqMR300Rsq0p09MR200","","HLT"))->setComment("PF HLT filter module used to save razor variable objects");
   desc.add<edm::InputTag>("TriggerResults",edm::InputTag("TriggerResults","","HLT"));
   desc.add<edm::InputTag>("trigSummary",edm::InputTag("hltTriggerSummaryAOD"));
   descriptions.add("SUSY_HLT_Razor_Main",desc);
 }