#include <AlCaGammaJetProducer.h>

Inheritance diagram for AlCaGammaJetProducer:

Public Member Functions
	AlCaGammaJetProducer (const edm::ParameterSet &)

virtual void	beginJob ()

virtual void	produce (edm::Event &, const edm::EventSetup &)

	~AlCaGammaJetProducer ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	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 ()

Private Attributes
bool	allowMissingInputs_

std::string	correctedIslandBarrelSuperClusterCollection_

std::string	correctedIslandBarrelSuperClusterProducer_

std::string	correctedIslandEndcapSuperClusterCollection_

std::string	correctedIslandEndcapSuperClusterProducer_

std::vector< edm::InputTag >	ecalLabels_

const CaloGeometry *	geo

std::vector< edm::InputTag >	mInputCalo

edm::EDGetTokenT < reco::SuperClusterCollection >	tok_EBSC_

edm::EDGetTokenT < reco::SuperClusterCollection >	tok_EESC_

edm::EDGetTokenT < HBHERecHitCollection >	tok_hbhe_

edm::EDGetTokenT < HFRecHitCollection >	tok_hf_

edm::EDGetTokenT < HORecHitCollection >	tok_ho_

edm::EDGetTokenT < reco::TrackCollection >	tok_inputTrack_

std::vector< edm::EDGetTokenT < reco::CaloJetCollection > >	toks_calo_

std::vector< edm::EDGetTokenT < EcalRecHitCollection > >	toks_ecal_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

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::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

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

Definition at line 41 of file AlCaGammaJetProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 13 of file AlCaGammaJetProducer.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), i, HLT_25ns14e33_v1_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

 {
    // Take input 
    
    tok_hbhe_= consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("hbheInput"));
    tok_ho_ = consumes<HORecHitCollection>(iConfig.getParameter<edm::InputTag>("hoInput"));
    tok_hf_ = consumes<HFRecHitCollection>(iConfig.getParameter<edm::InputTag>("hfInput"));
 
    mInputCalo = iConfig.getParameter<std::vector<edm::InputTag> >("srcCalo");
    ecalLabels_=iConfig.getParameter<std::vector<edm::InputTag> >("ecalInputs");
 
    unsigned nLabels = mInputCalo.size();
    for ( unsigned i=0; i != nLabels; i++ ) 
       toks_calo_.push_back( consumes<reco::CaloJetCollection>( mInputCalo[i] ) );
 
    nLabels = ecalLabels_.size();
    for ( unsigned i=0; i != nLabels; i++ )
       toks_ecal_.push_back( consumes<EcalRecHitCollection>( ecalLabels_[i] ) );
 
 
    tok_inputTrack_ =  consumes<reco::TrackCollection>(iConfig.getUntrackedParameter<std::string>("inputTrackLabel","generalTracks"));
    correctedIslandBarrelSuperClusterCollection_ = iConfig.getParameter<std::string>("correctedIslandBarrelSuperClusterCollection");
    correctedIslandBarrelSuperClusterProducer_   = iConfig.getParameter<std::string>("correctedIslandBarrelSuperClusterProducer");
     tok_EBSC_ = consumes<reco::SuperClusterCollection>(
         edm::InputTag(correctedIslandBarrelSuperClusterProducer_,
     correctedIslandBarrelSuperClusterCollection_)); 
 
    correctedIslandEndcapSuperClusterCollection_ = iConfig.getParameter<std::string>("correctedIslandEndcapSuperClusterCollection");
    correctedIslandEndcapSuperClusterProducer_   = iConfig.getParameter<std::string>("correctedIslandEndcapSuperClusterProducer");  
    tok_EESC_ = consumes<reco::SuperClusterCollection>(
     edm::InputTag(correctedIslandEndcapSuperClusterProducer_,
     correctedIslandEndcapSuperClusterCollection_));
 
    allowMissingInputs_=iConfig.getUntrackedParameter<bool>("AllowMissingInputs",true);
    
     
    //register your products
    produces<reco::TrackCollection>("GammaJetTracksCollection");
    produces<EcalRecHitCollection>("GammaJetEcalRecHitCollection");
    produces<HBHERecHitCollection>("GammaJetHBHERecHitCollection");
    produces<HORecHitCollection>("GammaJetHORecHitCollection");
    produces<HFRecHitCollection>("GammaJetHFRecHitCollection");
    produces<CaloJetCollection>("GammaJetJetBackToBackCollection");
    produces<reco::SuperClusterCollection>("GammaJetGammaBackToBackCollection");
 
     
    
 }

AlCaGammaJetProducer::~AlCaGammaJetProducer ( )

Definition at line 65 of file AlCaGammaJetProducer.cc.

 {
  
 }

Member Function Documentation

void AlCaGammaJetProducer::beginJob ( void )

virtual

Reimplemented from edm::EDProducer.

Definition at line 61 of file AlCaGammaJetProducer.cc.

62 {

63 }

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

virtual

Implements edm::EDProducer.

Definition at line 74 of file AlCaGammaJetProducer.cc.

References edm::SortedCollection< T, SORT >::begin(), correctedIslandBarrelSuperClusters_cfi::correctedIslandBarrelSuperClusters, correctedIslandEndcapSuperClusters_cfi::correctedIslandEndcapSuperClusters, gather_cfg::cout, edm::SortedCollection< T, SORT >::end(), PV3DBase< T, PVType, FrameType >::eta(), edm::EventSetup::get(), edm::Event::getByToken(), i, edm::HandleBase::isValid(), metsig::jet, fwrapper::jets, PV3DBase< T, PVType, FrameType >::phi(), edm::Event::put(), query::result, mathSSE::sqrt(), and HLT_50ns_5e33_v1_cff::trackCollection.

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

 {
    using namespace edm;
    using namespace std;
 
    edm::ESHandle<CaloGeometry> pG;
    iSetup.get<CaloGeometryRecord>().get(pG);
    geo = pG.product();
 
 // Produced collections 
   
   std::auto_ptr<reco::SuperClusterCollection> result (new reco::SuperClusterCollection); //Corrected jets
   std::auto_ptr<CaloJetCollection> resultjet (new CaloJetCollection); //Corrected jets
   std::auto_ptr<EcalRecHitCollection> miniEcalRecHitCollection(new EcalRecHitCollection);
   std::auto_ptr<HBHERecHitCollection> miniHBHERecHitCollection(new HBHERecHitCollection);
   std::auto_ptr<HORecHitCollection> miniHORecHitCollection(new HORecHitCollection);
   std::auto_ptr<HFRecHitCollection> miniHFRecHitCollection(new HFRecHitCollection);
   std::auto_ptr<reco::TrackCollection> outputTColl(new reco::TrackCollection);
    
    
 // Get  Corrected supercluster collection
   int nclusb = 0;
   int ncluse = 0;
   reco::SuperClusterCollection::const_iterator maxclusbarrel;
   reco::SuperClusterCollection::const_iterator maxclusendcap;
 
   double vetmax = -100.;
   Handle<reco::SuperClusterCollection> pCorrectedIslandBarrelSuperClusters;
   iEvent.getByToken(tok_EBSC_,
             pCorrectedIslandBarrelSuperClusters);  
   if (!pCorrectedIslandBarrelSuperClusters.isValid()) {
     // can't find it!
     if (!allowMissingInputs_) {
       *pCorrectedIslandBarrelSuperClusters;  // will throw the proper exception
     }
   } else {
     const reco::SuperClusterCollection* correctedIslandBarrelSuperClusters = pCorrectedIslandBarrelSuperClusters.product();
     
     // loop over the super clusters and find the highest
     maxclusbarrel = correctedIslandBarrelSuperClusters->begin();
     for(reco::SuperClusterCollection::const_iterator aClus = correctedIslandBarrelSuperClusters->begin();
     aClus != correctedIslandBarrelSuperClusters->end(); aClus++) {
       double vet = aClus->energy()/cosh(aClus->eta());
           cout<<" Barrel supercluster " << vet <<" energy "<<aClus->energy()<<" eta "<<aClus->eta()<<endl;
       if(vet>20.) {
     if(vet > vetmax)
       {
         vetmax = vet;
         maxclusbarrel = aClus;
         nclusb = 1;
       }
       }
     }
     
   }
 
   Handle<reco::SuperClusterCollection> pCorrectedIslandEndcapSuperClusters;
   iEvent.getByToken(tok_EESC_,
             pCorrectedIslandEndcapSuperClusters);  
   if (!pCorrectedIslandEndcapSuperClusters.isValid()) {
     // can't find it!
     if (!allowMissingInputs_) {
       *pCorrectedIslandEndcapSuperClusters;  // will throw the proper exception
     }
   } else {
     const reco::SuperClusterCollection* correctedIslandEndcapSuperClusters = pCorrectedIslandEndcapSuperClusters.product();
     
     // loop over the super clusters and find the highest
     maxclusendcap = correctedIslandEndcapSuperClusters->end();
     double vetmaxe = vetmax;
     for(reco::SuperClusterCollection::const_iterator aClus = correctedIslandEndcapSuperClusters->begin();
     aClus != correctedIslandEndcapSuperClusters->end(); aClus++) {
       double vet = aClus->energy()/cosh(aClus->eta());
       //   cout<<" Endcap supercluster " << vet <<" energy "<<aClus->energy()<<" eta "<<aClus->eta()<<endl;
       if(vet>20.) {
     if(vet > vetmaxe)
       {
         vetmaxe = vet;
         maxclusendcap = aClus;
         ncluse = 1;
       }
       }
     }
   }
    
   cout<<" Number of gammas "<<nclusb<<" "<<ncluse<<endl;  
   
   if( nclusb == 0 && ncluse == 0 ) {
    
   iEvent.put( outputTColl, "GammaJetTracksCollection");
   iEvent.put( miniEcalRecHitCollection, "GammaJetEcalRecHitCollection");
   iEvent.put( miniHBHERecHitCollection, "GammaJetHBHERecHitCollection");
   iEvent.put( miniHORecHitCollection, "GammaJetHORecHitCollection");
   iEvent.put( miniHFRecHitCollection, "GammaJetHFRecHitCollection");
   iEvent.put( result, "GammaJetGammaBackToBackCollection");
   iEvent.put( resultjet, "GammaJetJetBackToBackCollection");
   
   return;
   }
 
   double phigamma = -100.;
   double etagamma = -100.;
   if(ncluse == 1)
   {
     result->push_back(*maxclusendcap);
     phigamma = (*maxclusendcap).phi();
     etagamma = (*maxclusendcap).eta();
     
   } else
   {
     result->push_back(*maxclusbarrel);
     phigamma = (*maxclusbarrel).phi();
     etagamma = (*maxclusbarrel).eta();
   }
   
 //  cout<<" Size of egamma clusters "<<result->size()<<endl;
    
 //  
 // Jet Collection
 // Find jet in the angle ~ +- 170 degrees
 //
     
     double phijet =  -100.;
     double etajet = -100.;
     double phijet0 =  -100.;
     double etajet0 = -100.;
     
     std::vector<edm::EDGetTokenT<reco::CaloJetCollection> >::const_iterator ic;
     for (ic=toks_calo_.begin(); ic!=toks_calo_.end(); ic++) {
 
       edm::Handle<reco::CaloJetCollection> jets;
       iEvent.getByToken(*ic, jets);
       if (!jets.isValid()) {
     // can't find it!
     if (!allowMissingInputs_) {
       *jets;  // will throw the proper exception
     }
       } else {
        reco::CaloJetCollection::const_iterator jet = jets->begin ();
 
  //      cout<<" Size of jets "<<jets->size()<<endl;
 
        if( jets->size() == 0 ) continue;
 
          int iejet = 0;   
          int numjet = 0;
  
          for (; jet != jets->end (); jet++)
          {
            phijet0 = (*jet).phi();
            etajet0 = (*jet).eta();
 
 // Only 3 jets are kept
            numjet++;
            if(numjet > 3) break;
 
            cout<<" phi,eta "<< phigamma<<" "<< etagamma<<" "<<phijet0<<" "<<etajet0<<endl;
 
 // Find jet back to gamma
        
        double dphi = fabs(phigamma-phijet0); 
        if(dphi > 4.*atan(1.)) dphi = 8.*atan(1.) - dphi;
            double deta = fabs(etagamma-etajet0);
            double dr = sqrt(dphi*dphi+deta*deta);
            if(dr < 0.5 ) continue;
            resultjet->push_back ((*jet));
 
        dphi = dphi*180./(4.*atan(1.));
        if( fabs(dphi-180) < 30. )
        {
    //           cout<<" Jet is found "<<endl;
               iejet = 1;
               phijet = (*jet).phi();
               etajet = (*jet).eta();
        } // dphi
          } //jet collection
          if(iejet == 0) resultjet->clear(); 
       }
     } // Jet collection
      
      if( resultjet->size() == 0 ) {
 
       iEvent.put( outputTColl, "GammaJetTracksCollection");
       iEvent.put( miniEcalRecHitCollection, "GammaJetEcalRecHitCollection");
       iEvent.put( miniHBHERecHitCollection, "GammaJetHBHERecHitCollection");
       iEvent.put( miniHORecHitCollection, "GammaJetHORecHitCollection");
       iEvent.put( miniHFRecHitCollection, "GammaJetHFRecHitCollection");
       iEvent.put( result, "GammaJetGammaBackToBackCollection");
       iEvent.put( resultjet, "GammaJetJetBackToBackCollection");
 
      return;
      }
     // cout<<" Accepted event "<<resultjet->size()<<" PHI gamma "<<phigamma<<std::endl;
 //
 // Add Ecal, Hcal RecHits around Egamma caluster
 //
 
 // Load EcalRecHits
 
     
     std::vector<edm::EDGetTokenT<EcalRecHitCollection> >::const_iterator i;
     for (i=toks_ecal_.begin(); i!=toks_ecal_.end(); i++) {
       edm::Handle<EcalRecHitCollection> ec;
       iEvent.getByToken(*i,ec);
       if (!ec.isValid()) {
     // can't find it!
     if (!allowMissingInputs_) { 
       cout<<" No ECAL input "<<endl;
       *ec;  // will throw the proper exception
     }
       } else {
     for(EcalRecHitCollection::const_iterator recHit = (*ec).begin();
         recHit != (*ec).end(); ++recHit)
       {
 // EcalBarrel = 1, EcalEndcap = 2
         GlobalPoint pos = geo->getPosition(recHit->detid());
         double phihit = pos.phi();
         double etahit = pos.eta();
         
         double dphi = fabs(phigamma - phihit); 
         if(dphi > 4.*atan(1.)) dphi = 8.*atan(1.) - dphi;
         double deta = fabs(etagamma - etahit); 
         double dr = sqrt(dphi*dphi + deta*deta);
         if(dr<1.)  miniEcalRecHitCollection->push_back(*recHit);
         
         dphi = fabs(phijet - phihit); 
         if(dphi > 4.*atan(1.)) dphi = 8.*atan(1.) - dphi;
         deta = fabs(etajet - etahit); 
         dr = sqrt(dphi*dphi + deta*deta);
         if(dr<1.4)  miniEcalRecHitCollection->push_back(*recHit);
         
       }
       }
     }
 
 //   cout<<" Ecal is done "<<endl; 
 
       edm::Handle<HBHERecHitCollection> hbhe;
       iEvent.getByToken(tok_hbhe_,hbhe);
       if (!hbhe.isValid()) {
     // can't find it!
     if (!allowMissingInputs_) {
       *hbhe;  // will throw the proper exception
     }
       } else {
     const HBHERecHitCollection Hithbhe = *(hbhe.product());
     for(HBHERecHitCollection::const_iterator hbheItr=Hithbhe.begin(); hbheItr!=Hithbhe.end(); hbheItr++)
       {
         GlobalPoint pos = geo->getPosition(hbheItr->detid());
         double phihit = pos.phi();
         double etahit = pos.eta();
         
         double dphi = fabs(phigamma - phihit); 
         if(dphi > 4.*atan(1.)) dphi = 8.*atan(1.) - dphi;
         double deta = fabs(etagamma - etahit); 
         double dr = sqrt(dphi*dphi + deta*deta);
         
         
         if(dr<1.)  miniHBHERecHitCollection->push_back(*hbheItr);
         dphi = fabs(phijet - phihit); 
         if(dphi > 4.*atan(1.)) dphi = 8.*atan(1.) - dphi;
         deta = fabs(etajet - etahit); 
         dr = sqrt(dphi*dphi + deta*deta);
         if(dr<1.4)  miniHBHERecHitCollection->push_back(*hbheItr);
       }
       }
 
 //   cout<<" HBHE is done "<<endl; 
     
       edm::Handle<HORecHitCollection> ho;
       iEvent.getByToken(tok_ho_,ho);
       if (!ho.isValid()) {
     // can't find it!
     if (!allowMissingInputs_) {
       *ho;  // will throw the proper exception
     }
       } else {
     const HORecHitCollection Hitho = *(ho.product());
     for(HORecHitCollection::const_iterator hoItr=Hitho.begin(); hoItr!=Hitho.end(); hoItr++)
       {
         GlobalPoint pos = geo->getPosition(hoItr->detid());
         double phihit = pos.phi();
         double etahit = pos.eta();
         
         double dphi = fabs(phigamma - phihit); 
         if(dphi > 4.*atan(1.)) dphi = 8.*atan(1.) - dphi;
         double deta = fabs(etagamma - etahit); 
         double dr = sqrt(dphi*dphi + deta*deta);
         
         if(dr<1.)  miniHORecHitCollection->push_back(*hoItr);
         dphi = fabs(phijet - phihit); 
         if(dphi > 4.*atan(1.)) dphi = 8.*atan(1.) - dphi;
         deta = fabs(etajet - etahit); 
         dr = sqrt(dphi*dphi + deta*deta);
         if(dr<1.4)  miniHORecHitCollection->push_back(*hoItr);
         
       }
       }
 
  //  cout<<" HO is done "<<endl; 
    
       edm::Handle<HFRecHitCollection> hf;
       iEvent.getByToken(tok_hf_,hf);
       if (!hf.isValid()) {
     // can't find it!
     if (!allowMissingInputs_) {
       *hf;  // will throw the proper exception
     }
       } else {
     const HFRecHitCollection Hithf = *(hf.product());
     //  cout<<" Size of HF collection "<<Hithf.size()<<endl;
     for(HFRecHitCollection::const_iterator hfItr=Hithf.begin(); hfItr!=Hithf.end(); hfItr++)
       {
         GlobalPoint pos = geo->getPosition(hfItr->detid());
         
         double phihit = pos.phi();
         double etahit = pos.eta();
         
         double dphi = fabs(phigamma - phihit); 
         if(dphi > 4.*atan(1.)) dphi = 8.*atan(1.) - dphi;
         double deta = fabs(etagamma - etahit); 
         double dr = sqrt(dphi*dphi + deta*deta);
         
         if(dr<1.)  miniHFRecHitCollection->push_back(*hfItr);
         dphi = fabs(phijet - phihit); 
         if(dphi > 4.*atan(1.)) dphi = 8.*atan(1.) - dphi;
         deta = fabs(etajet - etahit); 
         dr = sqrt(dphi*dphi + deta*deta);
         
         if( dr < 1.4 )  miniHFRecHitCollection->push_back(*hfItr);
       }
       }
 
  //  cout<<" Size of mini HF collection "<<miniHFRecHitCollection->size()<<endl;
      
 
 // Track Collection   
       edm::Handle<reco::TrackCollection> trackCollection;
       iEvent.getByToken(tok_inputTrack_,trackCollection);
       if (!trackCollection.isValid()) {
     // can't find it!
     if (!allowMissingInputs_) {
       *trackCollection;  // will throw the proper exception
     }
       } else {
     const reco::TrackCollection tC = *(trackCollection.product());
   
     //Create empty output collections
 
 
     for (reco::TrackCollection::const_iterator track=tC.begin(); track!=tC.end(); track++)
       {
         
         double deta = track->momentum().eta() - etagamma; 
         
         double dphi = fabs(track->momentum().phi() - phigamma);
         
         if (dphi > atan(1.)*4.) dphi = 8.*atan(1.) - dphi;
         double ddir1 = sqrt(deta*deta+dphi*dphi);
         
         
         
         deta = track->momentum().eta() - etajet;
         dphi = fabs(track->momentum().phi() - phijet);
         if (dphi > atan(1.)*4.) dphi = 8.*atan(1.) - dphi;
         double ddir2 = sqrt(deta*deta+dphi*dphi);
         
         if( ddir1 < 1.4  || ddir2 < 1.4)      
           {
         outputTColl->push_back(*track);
           } 
       }
       }
    
   //Put selected information in the event
   
   iEvent.put( outputTColl, "GammaJetTracksCollection");
   iEvent.put( miniEcalRecHitCollection, "GammaJetEcalRecHitCollection");
   iEvent.put( miniHBHERecHitCollection, "GammaJetHBHERecHitCollection");
   iEvent.put( miniHORecHitCollection, "GammaJetHORecHitCollection");
   iEvent.put( miniHFRecHitCollection, "GammaJetHFRecHitCollection");
   iEvent.put( result, "GammaJetGammaBackToBackCollection");
   iEvent.put( resultjet, "GammaJetJetBackToBackCollection");
 }

Member Data Documentation

bool AlCaGammaJetProducer::allowMissingInputs_

private

Definition at line 68 of file AlCaGammaJetProducer.h.

std::string AlCaGammaJetProducer::correctedIslandBarrelSuperClusterCollection_

private

Definition at line 56 of file AlCaGammaJetProducer.h.

std::string AlCaGammaJetProducer::correctedIslandBarrelSuperClusterProducer_

private

Definition at line 57 of file AlCaGammaJetProducer.h.

std::string AlCaGammaJetProducer::correctedIslandEndcapSuperClusterCollection_

private

Definition at line 58 of file AlCaGammaJetProducer.h.

std::string AlCaGammaJetProducer::correctedIslandEndcapSuperClusterProducer_

private

Definition at line 59 of file AlCaGammaJetProducer.h.

std::vector<edm::InputTag> AlCaGammaJetProducer::ecalLabels_

private

Definition at line 50 of file AlCaGammaJetProducer.h.

const CaloGeometry* AlCaGammaJetProducer::geo

private

Definition at line 71 of file AlCaGammaJetProducer.h.

std::vector<edm::InputTag> AlCaGammaJetProducer::mInputCalo

private

Definition at line 51 of file AlCaGammaJetProducer.h.

edm::EDGetTokenT<reco::SuperClusterCollection> AlCaGammaJetProducer::tok_EBSC_

private

Definition at line 65 of file AlCaGammaJetProducer.h.

edm::EDGetTokenT<reco::SuperClusterCollection> AlCaGammaJetProducer::tok_EESC_

private

Definition at line 66 of file AlCaGammaJetProducer.h.

edm::EDGetTokenT<HBHERecHitCollection> AlCaGammaJetProducer::tok_hbhe_

private

Definition at line 61 of file AlCaGammaJetProducer.h.

edm::EDGetTokenT<HFRecHitCollection> AlCaGammaJetProducer::tok_hf_

private

Definition at line 63 of file AlCaGammaJetProducer.h.

edm::EDGetTokenT<HORecHitCollection> AlCaGammaJetProducer::tok_ho_

private

Definition at line 62 of file AlCaGammaJetProducer.h.

edm::EDGetTokenT<reco::TrackCollection> AlCaGammaJetProducer::tok_inputTrack_

private

Definition at line 64 of file AlCaGammaJetProducer.h.

std::vector<edm::EDGetTokenT<reco::CaloJetCollection> > AlCaGammaJetProducer::toks_calo_

private

Definition at line 54 of file AlCaGammaJetProducer.h.

std::vector<edm::EDGetTokenT<EcalRecHitCollection> > AlCaGammaJetProducer::toks_ecal_

private

Definition at line 53 of file AlCaGammaJetProducer.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation