#include <PhotonCoreProducer.h>

Inheritance diagram for PhotonCoreProducer:

Public Member Functions
virtual void	beginRun (edm::Run &r, edm::EventSetup const &es)

virtual void	endRun (edm::Run &, edm::EventSetup const &)

	PhotonCoreProducer (const edm::ParameterSet &ps)

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

	~PhotonCoreProducer ()

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 Member Functions
void	fillPhotonCollection (edm::Event &evt, edm::EventSetup const &es, const edm::Handle< reco::SuperClusterCollection > &scHandle, const edm::Handle< reco::ConversionCollection > &conversionHandle, const edm::Handle< reco::ElectronSeedCollection > &pixelSeeds, reco::PhotonCoreCollection &outputCollection, int &iSC)

reco::ConversionRef	solveAmbiguity (const edm::Handle< reco::ConversionCollection > &conversionHandle, reco::SuperClusterRef &sc)

Private Attributes
edm::ParameterSet	conf_

edm::InputTag	conversionProducer_

double	minSCEt_

std::string	PhotonCoreCollection_

std::string	pixelSeedProducer_

bool	risolveAmbiguity_

edm::InputTag	scHybridBarrelCollection_

edm::InputTag	scHybridBarrelProducer_

edm::InputTag	scIslandEndcapCollection_

edm::InputTag	scIslandEndcapProducer_

bool	validConversions_

bool	validPixelSeeds_

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

Id:: PhotonCoreProducer.h,v 1.3 2010/02/22 19:30:49 nancy Exp

Date:: 2010/02/22 19:30:49

Revision:: 1.3

Author: Nancy Marinelli, U. of Notre Dame, US

Definition at line 34 of file PhotonCoreProducer.h.

Constructor & Destructor Documentation

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

Definition at line 20 of file PhotonCoreProducer.cc.

References conf_, conversionProducer_, edm::ParameterSet::getParameter(), minSCEt_, PhotonCoreCollection_, pixelSeedProducer_, risolveAmbiguity_, scHybridBarrelProducer_, and scIslandEndcapProducer_.

                                                                     : 
   conf_(config)
 
 {
 
   // use onfiguration file to setup input/output collection names
   scHybridBarrelProducer_       = conf_.getParameter<edm::InputTag>("scHybridBarrelProducer");
   scIslandEndcapProducer_       = conf_.getParameter<edm::InputTag>("scIslandEndcapProducer");
   conversionProducer_ = conf_.getParameter<edm::InputTag>("conversionProducer");
   PhotonCoreCollection_ = conf_.getParameter<std::string>("photonCoreCollection");
   pixelSeedProducer_   = conf_.getParameter<std::string>("pixelSeedProducer");
   minSCEt_        = conf_.getParameter<double>("minSCEt");
   risolveAmbiguity_ = conf_.getParameter<bool>("risolveConversionAmbiguity");
 
   // Register the product
   produces< reco::PhotonCoreCollection >(PhotonCoreCollection_);
 
 }

PhotonCoreProducer::~PhotonCoreProducer ( )

Definition at line 39 of file PhotonCoreProducer.cc.

39 {}

Member Function Documentation

void PhotonCoreProducer::beginRun	(	edm::Run &	r,
		edm::EventSetup const &	es
	)

virtual

Reimplemented from edm::EDProducer.

Definition at line 43 of file PhotonCoreProducer.cc.

                                                                                   {
 
 }

void PhotonCoreProducer::endRun	(	edm::Run &	r,
		edm::EventSetup const &	theEventSetup
	)

virtual

Reimplemented from edm::EDProducer.

Definition at line 47 of file PhotonCoreProducer.cc.

                                                                                 {
 
 }

void PhotonCoreProducer::fillPhotonCollection	(	edm::Event &	evt,
		edm::EventSetup const &	es,
		const edm::Handle< reco::SuperClusterCollection > &	scHandle,
		const edm::Handle< reco::ConversionCollection > &	conversionHandle,
		const edm::Handle< reco::ElectronSeedCollection > &	pixelSeeds,
		reco::PhotonCoreCollection &	outputCollection,
		int &	iSC
	)

private

Definition at line 131 of file PhotonCoreProducer.cc.

References reco::PhotonCore::addConversion(), reco::PhotonCore::addElectronPixelSeed(), edm::Ref< C, T, F >::id(), edm::Ref< C, T, F >::isNonnull(), edm::Ref< C, T, F >::key(), minSCEt_, risolveAmbiguity_, solveAmbiguity(), validConversions_, and validPixelSeeds_.

Referenced by produce().

                                                                                            {
   
 
   
   reco::ElectronSeedCollection::const_iterator pixelSeedItr;
   for(unsigned int lSC=0; lSC < scHandle->size(); lSC++) {
     
     // get SuperClusterRef
     reco::SuperClusterRef scRef(reco::SuperClusterRef(scHandle, lSC));
     iSC++;
     //const reco::SuperCluster* pClus=&(*scRef);
     
     // SC energy preselection
     if (scRef->energy()/cosh(scRef->eta()) <= minSCEt_) continue;
     
     reco::PhotonCore newCandidate(scRef);
     if ( validConversions_) {    
 
       if ( risolveAmbiguity_ ) {
     reco::ConversionRef bestRef=solveAmbiguity( conversionHandle , scRef);  
     if (bestRef.isNonnull() ) newCandidate.addConversion(bestRef);     
     
       } else {
     
 
     for( unsigned int icp = 0;  icp < conversionHandle->size(); icp++) {
       reco::ConversionRef cpRef(reco::ConversionRef(conversionHandle,icp));
           if ( !cpRef->caloCluster().size()) continue; 
       if (!( scRef.id() == cpRef->caloCluster()[0].id() && scRef.key() == cpRef->caloCluster()[0].key() )) continue; 
       if ( !cpRef->isConverted() ) continue;  
       newCandidate.addConversion(cpRef);     
 
     }     
     
       } // solve or not the ambiguity of many conversion candidates 
     
     }  
 
     if ( validPixelSeeds_) {
       for( unsigned int icp = 0;  icp < pixelSeedHandle->size(); icp++) {
     reco::ElectronSeedRef cpRef(reco::ElectronSeedRef(pixelSeedHandle,icp));
     if ( ! cpRef->isEcalDriven() ) continue;
     if (!( scRef.id() == cpRef->caloCluster().id() && scRef.key() == cpRef->caloCluster().key() )) continue; 
     newCandidate.addElectronPixelSeed(cpRef);     
       } 
     }
     
     outputPhotonCoreCollection.push_back(newCandidate);
    
   }    
   
 
 }

void PhotonCoreProducer::produce	(	edm::Event &	evt,
		const edm::EventSetup &	es
	)

virtual

Implements edm::EDProducer.

Definition at line 54 of file PhotonCoreProducer.cc.

References printConversionInfo::conversionHandle, conversionProducer_, fillPhotonCollection(), edm::Event::getByLabel(), edm::InputTag::label(), PhotonCoreCollection_, pixelSeedProducer_, edm::Event::put(), scHybridBarrelProducer_, scIslandEndcapProducer_, validConversions_, and validPixelSeeds_.

                                                                                        {
 
 
   using namespace edm;
   //  nEvt_++;
 
   reco::PhotonCoreCollection outputPhotonCoreCollection;
   std::auto_ptr< reco::PhotonCoreCollection > outputPhotonCoreCollection_p(new reco::PhotonCoreCollection);
 
   // Get the  Barrel Super Cluster collection
   bool validBarrelSCHandle=true;
   Handle<reco::SuperClusterCollection> scBarrelHandle;
   theEvent.getByLabel(scHybridBarrelProducer_,scBarrelHandle);
   if (!scBarrelHandle.isValid()) {
     edm::LogError("PhotonCoreProducer") << "Error! Can't get the product "<<scHybridBarrelProducer_.label();
     validBarrelSCHandle=false;
   }
 
 
  // Get the  Endcap Super Cluster collection
   bool validEndcapSCHandle=true;
   Handle<reco::SuperClusterCollection> scEndcapHandle;
   theEvent.getByLabel(scIslandEndcapProducer_,scEndcapHandle);
   if (!scEndcapHandle.isValid()) {
     edm::LogError("PhotonCoreProducer") << "Error! Can't get the product "<<scIslandEndcapProducer_.label();
     validEndcapSCHandle=false;
   }
 
 
   validConversions_=true;
   edm::Handle<reco::ConversionCollection> conversionHandle; 
   theEvent.getByLabel(conversionProducer_, conversionHandle);
   if (!conversionHandle.isValid()) {
     //edm::LogError("PhotonCoreProducer") << "Error! Can't get the product "<< conversionProducer_.label() << "\n" ;
     validConversions_=false;
   }
  
 
 
 
   // Get ElectronPixelSeeds
   validPixelSeeds_=true;
   Handle<reco::ElectronSeedCollection> pixelSeedHandle;
   reco::ElectronSeedCollection pixelSeeds;
   theEvent.getByLabel(pixelSeedProducer_, pixelSeedHandle);
   if (!pixelSeedHandle.isValid()) {
     validPixelSeeds_=false;
   }
   //  if ( validPixelSeeds_) pixelSeeds = *(pixelSeedHandle.product());
 
 
 
   int iSC=0; // index in photon collection
   // Loop over barrel and endcap SC collections and fill the  photon collection
   if ( validBarrelSCHandle) fillPhotonCollection(theEvent,
                          theEventSetup,
                          scBarrelHandle,
                          conversionHandle,
                          pixelSeedHandle,
                          outputPhotonCoreCollection,
                          iSC);
   if ( validEndcapSCHandle) fillPhotonCollection(theEvent,
                          theEventSetup,
                          scEndcapHandle,
                          conversionHandle,
                          pixelSeedHandle,
                          outputPhotonCoreCollection,
                          iSC);
 
   // put the product in the event
   edm::LogInfo("PhotonCoreProducer") << " Put in the event " << iSC << " Photon Candidates \n";
   outputPhotonCoreCollection_p->assign(outputPhotonCoreCollection.begin(),outputPhotonCoreCollection.end());
   theEvent.put( outputPhotonCoreCollection_p, PhotonCoreCollection_);
 
 }

reco::ConversionRef PhotonCoreProducer::solveAmbiguity	(	const edm::Handle< reco::ConversionCollection > &	conversionHandle,
		reco::SuperClusterRef &	sc
	)

private

only one track reconstructed. Pick the one with best E/P

Definition at line 193 of file PhotonCoreProducer.cc.

References edm::Ref< C, T, F >::id(), edm::Ref< C, T, F >::key(), AlCaHLTBitMon_ParallelJobs::p, mathSSE::sqrt(), and testEve_cfg::tracks.

Referenced by fillPhotonCollection().

                                                                                                                                               {
 
 
   std::multimap<reco::ConversionRef, double >   convMap;
   for ( unsigned int icp=0; icp< conversionHandle->size(); icp++) {
 
     reco::ConversionRef cpRef(reco::ConversionRef(conversionHandle,icp));
 
     if (!( scRef.id() == cpRef->caloCluster()[0].id() && scRef.key() == cpRef->caloCluster()[0].key() )) continue;    
     if ( !cpRef->isConverted() ) continue;  
     double like = cpRef->MVAout();
     convMap.insert ( std::make_pair(cpRef,like) );
 
   }          
   
   
   
   std::multimap<reco::ConversionRef, double >::iterator  iMap; 
   double max_lh = -1.;
   reco::ConversionRef bestRef;
   //  std::cout << " Pick up the best conv " << std::endl;
   for (iMap=convMap.begin();  iMap!=convMap.end(); iMap++) {
     double like = iMap->second;
     if (like > max_lh) { 
       max_lh = like;
       bestRef=iMap->first;
     }
   }            
   
   //std::cout << " Best conv like " << max_lh << std::endl;    
   
   float ep=0;
   if ( max_lh <0 ) {
     //    std::cout << " Candidates with only one track " << std::endl;
     float epMin=999; 
     
     for (iMap=convMap.begin();  iMap!=convMap.end(); iMap++) {
       reco::ConversionRef convRef=iMap->first;
       // std::vector<reco::TrackRef> tracks = convRef->tracks();    
       const std::vector<edm::RefToBase<reco::Track> > tracks = convRef->tracks();   
         float px=tracks[0]->innerMomentum().x();
         float py=tracks[0]->innerMomentum().y();
         float pz=tracks[0]->innerMomentum().z();
         float p=sqrt(px*px+py*py+pz*pz);
         ep=fabs(1.-convRef->caloCluster()[0]->energy()/p);
         //    std::cout << " 1-E/P = " << ep << std::endl;
         if ( ep<epMin) {
           epMin=ep;
           bestRef=iMap->first;
         }
     }
     //  std::cout << " Best conv 1-E/P " << ep << std::endl;    
             
   }
   
 
   return bestRef;
   
   
 } 