#include <AlCaDiJetsProducer.h>

Inheritance diagram for cms::AlCaDiJetsProducer:

Public Member Functions
	AlCaDiJetsProducer (const edm::ParameterSet &)
virtual void	beginJob ()
virtual void	produce (edm::Event &, const edm::EventSetup &)
	~AlCaDiJetsProducer ()
Private Attributes
bool	allowMissingInputs_
std::vector< edm::InputTag >	ecalLabels_
edm::InputTag	hbheInput_
edm::InputTag	hfInput_
edm::InputTag	hoInput_
edm::InputTag	jetsInput_

Detailed Description

Definition at line 38 of file AlCaDiJetsProducer.h.

Constructor & Destructor Documentation

cms::AlCaDiJetsProducer::AlCaDiJetsProducer ( const edm::ParameterSet & iConfig ) [explicit]

Definition at line 17 of file AlCaDiJetsProducer.cc.

References edm::ParameterSet::getParameter().

{
   jetsInput_ = iConfig.getParameter<edm::InputTag>("jetsInput");
   ecalLabels_=iConfig.getParameter<std::vector<edm::InputTag> >("ecalInputs");
   hbheInput_ = iConfig.getParameter<edm::InputTag>("hbheInput");
   hoInput_ = iConfig.getParameter<edm::InputTag>("hoInput");
   hfInput_ = iConfig.getParameter<edm::InputTag>("hfInput"); 
   allowMissingInputs_ = true;

//register your products
   produces<CaloJetCollection>("DiJetsBackToBackCollection");
   produces<EcalRecHitCollection>("DiJetsEcalRecHitCollection");
   produces<HBHERecHitCollection>("DiJetsHBHERecHitCollection");
   produces<HORecHitCollection>("DiJetsHORecHitCollection");
   produces<HFRecHitCollection>("DiJetsHFRecHitCollection");

}

cms::AlCaDiJetsProducer::~AlCaDiJetsProducer ( )

Definition at line 38 of file AlCaDiJetsProducer.cc.

{
 

}

Member Function Documentation

void cms::AlCaDiJetsProducer::beginJob ( void ) [virtual]

Reimplemented from edm::EDProducer.

Definition at line 34 of file AlCaDiJetsProducer.cc.

{
}

void cms::AlCaDiJetsProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)		`[virtual]`

Implements edm::EDProducer.

Definition at line 47 of file AlCaDiJetsProducer.cc.

References eta(), PV3DBase< T, PVType, FrameType >::eta(), edm::EventSetup::get(), edm::Event::getByLabel(), CaloGeometry::getPosition(), i, fwrapper::jets, reco::LeafCandidate::phi(), PV3DBase< T, PVType, FrameType >::phi(), phi, pi, pos, edm::ESHandle< T >::product(), edm::Event::put(), query::result, and mathSSE::sqrt().

{

  double pi = 4.*atan(1.);

   std::auto_ptr<CaloJetCollection> result (new CaloJetCollection); //Corrected jets
   std::auto_ptr<EcalRecHitCollection> miniDiJetsEcalRecHitCollection(new EcalRecHitCollection);
   std::auto_ptr<HBHERecHitCollection> miniDiJetsHBHERecHitCollection(new HBHERecHitCollection);
   std::auto_ptr<HORecHitCollection> miniDiJetsHORecHitCollection(new HORecHitCollection);
   std::auto_ptr<HFRecHitCollection> miniDiJetsHFRecHitCollection(new HFRecHitCollection);

   edm::ESHandle<CaloGeometry> pG;
   iSetup.get<CaloGeometryRecord>().get(pG);
   const CaloGeometry* geo = pG.product();

  // Jets Collections 

   vector<CaloJet> jetv; 

   CaloJet fJet1, fJet2, fJet3;
   edm::Handle<CaloJetCollection> jets;                
   iEvent.getByLabel(jetsInput_, jets);    
   int iflag_select = 0; 
   if(jets->size()>1){
    fJet1 = (*jets)[0];
    fJet2 = (*jets)[1];
    double dphi = fabs(fJet1.phi() - fJet2.phi());  
    if(dphi > pi){dphi = 2*pi - dphi;}
    double degreedphi = dphi*180./pi; 
    if(fabs(degreedphi-180)<30.){iflag_select = 1;}
   }
   if(iflag_select == 1){
     result->push_back(fJet1);
     result->push_back(fJet2);
     jetv.push_back(fJet1); 
     jetv.push_back(fJet2); 
     if(jets->size()>2){
     fJet3 = (*jets)[2]; 
     result->push_back(fJet3);
     jetv.push_back(fJet3);
     }     
   } else {
     iEvent.put( result, "DiJetsBackToBackCollection");
     iEvent.put( miniDiJetsEcalRecHitCollection,"DiJetsEcalRecHitCollection");
     iEvent.put( miniDiJetsHBHERecHitCollection, "DiJetsHBHERecHitCollection");
     iEvent.put( miniDiJetsHORecHitCollection, "DiJetsHORecHitCollection");
     iEvent.put( miniDiJetsHFRecHitCollection, "DiJetsHFRecHitCollection");
     return;
   }  
  
  // Ecal Collections 

   std::vector<edm::InputTag>::const_iterator i;
   for (i=ecalLabels_.begin(); i!=ecalLabels_.end(); i++) {
   edm::Handle<EcalRecHitCollection> ec;
   iEvent.getByLabel(*i,ec);
   for(EcalRecHitCollection::const_iterator ecItr = (*ec).begin();
                                                 ecItr != (*ec).end(); ++ecItr)
        {
// EcalBarrel = 1, EcalEndcap = 2
          GlobalPoint pos = geo->getPosition(ecItr->detid());
          double phihit = pos.phi();
          double etahit = pos.eta();
          int iflag_select = 0;  
          for(unsigned int i=0; i<jetv.size(); i++){
            double deta = fabs(etahit - jetv[i].eta());
            double dphi = fabs(phihit - jetv[i].phi());
            if(dphi > pi) dphi = 2*pi - dphi;
            double dr = sqrt(deta*deta+dphi*dphi);
            if(dr < 1.4){iflag_select = 1;}
          }  
          if(iflag_select==1){miniDiJetsEcalRecHitCollection->push_back(*ecItr);}
       }

   }

  // HB & HE Collections 
   
   edm::Handle<HBHERecHitCollection> hbhe;
   iEvent.getByLabel(hbheInput_,hbhe);
   for(HBHERecHitCollection::const_iterator hbheItr=hbhe->begin(); 
                                                  hbheItr!=hbhe->end(); hbheItr++)
        {
          GlobalPoint pos = geo->getPosition(hbheItr->detid());
          double phihit = pos.phi();
          double etahit = pos.eta();
          int iflag_select = 0;  
          for(unsigned int i=0; i<jetv.size(); i++){
            double deta = fabs(etahit - jetv[i].eta());
            double dphi = fabs(phihit - jetv[i].phi());
            if(dphi > pi) dphi = 2*pi - dphi;
            double dr = sqrt(deta*deta+dphi*dphi);
            if(dr < 1.4){iflag_select = 1;}
          }  
          if(iflag_select==1){miniDiJetsHBHERecHitCollection->push_back(*hbheItr);}
        }
  

// HO Collections 


   edm::Handle<HORecHitCollection> ho;
   iEvent.getByLabel(hoInput_,ho);
   for(HORecHitCollection::const_iterator hoItr=ho->begin(); 
                                                hoItr!=ho->end(); hoItr++)
     {
          GlobalPoint pos = geo->getPosition(hoItr->detid());
          double phihit = pos.phi();
          double etahit = pos.eta();
          int iflag_select = 0;  
          for(unsigned int i=0; i<jetv.size(); i++){
            double deta = fabs(etahit - jetv[i].eta());
            double dphi = fabs(phihit - jetv[i].phi());
            if(dphi > pi) dphi = 2*pi - dphi;
            double dr = sqrt(deta*deta+dphi*dphi);
            if(dr < 1.4){iflag_select = 1;}
          }  
          if(iflag_select==1){miniDiJetsHORecHitCollection->push_back(*hoItr);}
        }
  
  // HF Collection

 
   edm::Handle<HFRecHitCollection> hf;
   iEvent.getByLabel(hfInput_,hf);
   for(HFRecHitCollection::const_iterator hfItr=hf->begin(); 
                                                hfItr!=hf->end(); hfItr++)
       {
          GlobalPoint pos = geo->getPosition(hfItr->detid());
          double phihit = pos.phi();
          double etahit = pos.eta();
          int iflag_select = 0;  
          for(unsigned int i=0; i<jetv.size(); i++){
            double deta = fabs(etahit - jetv[i].eta());
            double dphi = fabs(phihit - jetv[i].phi());
            if(dphi > pi) dphi = 2*pi - dphi;
            double dr = sqrt(deta*deta+dphi*dphi);
            if(dr < 1.4){iflag_select = 1;}
          }  
          if(iflag_select==1){miniDiJetsHFRecHitCollection->push_back(*hfItr);}
       }
 

  //Put selected information in the event

   iEvent.put( result, "DiJetsBackToBackCollection");
   iEvent.put( miniDiJetsEcalRecHitCollection,"DiJetsEcalRecHitCollection");
   iEvent.put( miniDiJetsHBHERecHitCollection, "DiJetsHBHERecHitCollection");
   iEvent.put( miniDiJetsHORecHitCollection, "DiJetsHORecHitCollection");
   iEvent.put( miniDiJetsHFRecHitCollection, "DiJetsHFRecHitCollection");
}