DiJetVarProducer Class Reference

#include <DiJetVarProducer.h>

Inheritance diagram for DiJetVarProducer:

Public Member Functions
	DiJetVarProducer (const edm::ParameterSet &)
virtual void	produce (edm::Event &, const edm::EventSetup &)
	~DiJetVarProducer ()
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 ()

Private Attributes
edm::InputTag	inputJetTag_
double	wideJetDeltaR_

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
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
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

Definition at line 19 of file DiJetVarProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 24 of file DiJetVarProducer.cc.

References edm::InputTag::encode(), inputJetTag_, LogDebug, and wideJetDeltaR_.

                                                                 :
  inputJetTag_    (iConfig.getParameter<edm::InputTag>("inputJetTag")),
  wideJetDeltaR_   (iConfig.getParameter<double>("wideJetDeltaR"))
{
  //register your products
  //produces<std::vector<double> >("dijetvariables");
  produces<std::vector<math::PtEtaPhiMLorentzVector> >("widejets");
  
  LogDebug("") << "Input Jet Tag: " << inputJetTag_.encode() << " ";
  LogDebug("") << "Radius Parameter Wide Jet: " << wideJetDeltaR_ << ".";
}

DiJetVarProducer::~DiJetVarProducer

(

)

Definition at line 36 of file DiJetVarProducer.cc.

{
}

Member Function Documentation

void DiJetVarProducer::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

virtual

Implements edm::EDProducer.

Definition at line 42 of file DiJetVarProducer.cc.

References edm::Event::getByLabel(), inputJetTag_, edm::Event::put(), dt_dqm_sourceclient_common_cff::reco, and wideJetDeltaR_.

Referenced by JSONExport.JsonExport::export(), HTMLExport.HTMLExport::export(), and HTMLExport.HTMLExportStatic::export().

{
   using namespace std;
   using namespace edm;
   using namespace reco;

   // ## The output collections        
   //std::auto_ptr<std::vector<double> > dijetvariables(new std::vector<double>); 
   std::auto_ptr<std::vector<math::PtEtaPhiMLorentzVector> > widejets(new std::vector<math::PtEtaPhiMLorentzVector>); 

   // ## Get jet collection
   edm::Handle<reco::CaloJetCollection> calojets_handle;
   iEvent.getByLabel(inputJetTag_,calojets_handle);
   //cout << "size: " << calojets_handle->size() << endl;              
   
   // ## Wide Jet Algorithm
   // At least two jets
   if( calojets_handle->size() >=2 )
     {
       TLorentzVector wj1_tmp;
       TLorentzVector wj2_tmp;
       TLorentzVector wj1;
       TLorentzVector wj2;
       TLorentzVector wdijet;
       
       //        // Loop over all the input jets
       //        for(reco::CaloJetCollection::const_iterator it = calojets_handle->begin(); it != calojets_handle->end(); ++it)
       //            {
       //      cout << "jet: " << it->pt() << " " << it->eta() << " " << it->phi() << endl;
       //            }
       
       // Find two leading jets
       TLorentzVector jet1, jet2;

       reco::CaloJetCollection::const_iterator j1 = calojets_handle->begin();
       reco::CaloJetCollection::const_iterator j2 = j1; ++j2;       

       jet1.SetPtEtaPhiM(j1->pt(),j1->eta(),j1->phi(),j1->mass());
       jet2.SetPtEtaPhiM(j2->pt(),j2->eta(),j2->phi(),j2->mass());       

       //cout << "j1: " << jet1.Pt() << " " << jet1.Eta() << " " << jet1.Phi() << endl;
       //cout << "j2: " << jet2.Pt() << " " << jet2.Eta() << " " << jet2.Phi() << endl;

       // Create wide jets (radiation recovery algorithm)
       for(reco::CaloJetCollection::const_iterator it = calojets_handle->begin(); it != calojets_handle->end(); ++it)
     {
       TLorentzVector currentJet;
       currentJet.SetPtEtaPhiM(it->pt(),it->eta(),it->phi(),it->mass());

       double DeltaR1 = currentJet.DeltaR(jet1);
       double DeltaR2 = currentJet.DeltaR(jet2);

       if(DeltaR1 < DeltaR2 && DeltaR1 < wideJetDeltaR_) {
         wj1_tmp += currentJet;
       } else if(DeltaR2 < wideJetDeltaR_) {
         wj2_tmp += currentJet;
       }
     }       

       // Re-order the wide jets in pT
       if( wj1_tmp.Pt() > wj2_tmp.Pt() )    
     { 
       wj1 = wj1_tmp;
       wj2 = wj2_tmp;
     }
       else
     {
       wj1 = wj2_tmp;
       wj2 = wj1_tmp;
     }

       // Create dijet system
       wdijet = wj1 + wj2;

       //        cout << "j1 wide: " << wj1.Pt() << " " << wj1.Eta() << " " << wj1.Phi() << " " << wj1.M() << endl;
       //        cout << "j2 wide: " << wj2.Pt() << " " << wj2.Eta() << " " << wj2.Phi() << " " << wj2.M() << endl;
       //        cout << "MJJWide: " << wdijet.M() << endl;
       //        cout << "DeltaEtaJJWide: " << fabs(wj1.Eta()-wj2.Eta()) << endl;
       //        cout << "DeltaPhiJJWide: " << fabs(wj1.DeltaPhi(wj2)) << endl;
 
       //        // Put variables in the container
       //        dijetvariables->push_back( wdijet.M() );                 //0 = MJJWide
       //        dijetvariables->push_back( fabs(wj1.Eta()-wj2.Eta()) );  //1 = DeltaEtaJJWide       
       //        dijetvariables->push_back( fabs(wj1.DeltaPhi(wj2)) );    //2 = DeltaPhiJJWide             

       // Put widejets in the container 
       math::PtEtaPhiMLorentzVector wj1math(wj1.Pt(), wj1.Eta(), wj1.Phi(), wj1.M());
       math::PtEtaPhiMLorentzVector wj2math(wj2.Pt(), wj2.Eta(), wj2.Phi(), wj2.M());
       widejets->push_back( wj1math );
       widejets->push_back( wj2math );
     }
   //    else
   //      {
   //        // Put variables in the container
   //        dijetvariables->push_back( -1 );                //0 = MJJWide
   //        dijetvariables->push_back( -1 );                //1 = DeltaEtaJJWide       
   //        dijetvariables->push_back( -1 );                //2 = DeltaPhiJJWide       
   //      }

   // ## Put objects in the Event
   //iEvent.put(dijetvariables, "dijetvariables");
   iEvent.put(widejets, "widejets");

}

Member Data Documentation

edm::InputTag DiJetVarProducer::inputJetTag_

private

Definition at line 27 of file DiJetVarProducer.h.

Referenced by DiJetVarProducer(), and produce().

double DiJetVarProducer::wideJetDeltaR_

private

Definition at line 28 of file DiJetVarProducer.h.

Referenced by DiJetVarProducer(), and produce().