GEDPhotonCoreProducer Class Reference

#include <GEDPhotonCoreProducer.h>

Inheritance diagram for GEDPhotonCoreProducer:

Public Member Functions
	GEDPhotonCoreProducer (const edm::ParameterSet &ps)
virtual void	produce (edm::Event &evt, const edm::EventSetup &es)
	~GEDPhotonCoreProducer ()
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
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
	EDConsumerBase ()
ProductHolderIndex	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndex > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndex > &) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Member Functions
void	createSingleLegConversions (reco::CaloClusterPtr, std::vector< reco::TrackRef >, std::vector< float >, reco::ConversionCollection &oneLegConversions)

Private Attributes
edm::ParameterSet	conf_
std::string	GEDPhotonCoreCollection_
double	minSCEt_
std::string	PFConversionCollection_
edm::InputTag	pfEgammaCandidates_
std::string	pixelSeedProducer_
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
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

Id:

GEDPhotonCoreProducer.h,v 1.1 2013/05/07 12:34:14 nancy Exp

Date:

2013/05/07 12:34:14

Revision:

1.1

Author

Nancy Marinelli, U. of Notre Dame, US

Id:

GEDPhotonCoreProducer.cc,v 1.1 2013/05/07 12:33:37 nancy Exp

Date:

2013/05/07 12:33:37

Revision:

1.1

Author

Nancy Marinelli, U. of Notre Dame, US

Definition at line 25 of file GEDPhotonCoreProducer.h.

Constructor & Destructor Documentation

GEDPhotonCoreProducer::GEDPhotonCoreProducer

(

const edm::ParameterSet &

ps

)

Definition at line 39 of file GEDPhotonCoreProducer.cc.

References conf_, GEDPhotonCoreCollection_, edm::ParameterSet::getParameter(), PFConversionCollection_, pfEgammaCandidates_, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                          : 
  conf_(config)

{

  // use onfiguration file to setup input/output collection names
  pfEgammaCandidates_      = conf_.getParameter<edm::InputTag>("pfEgammaCandidates");
  GEDPhotonCoreCollection_ = conf_.getParameter<std::string>("gedPhotonCoreCollection");

  // Register the product
  produces< reco::PhotonCoreCollection >(GEDPhotonCoreCollection_);
  produces<reco::ConversionCollection>(PFConversionCollection_);

}

GEDPhotonCoreProducer::~GEDPhotonCoreProducer

(

)

Definition at line 54 of file GEDPhotonCoreProducer.cc.

{}

Member Function Documentation

void GEDPhotonCoreProducer::createSingleLegConversions ( reco::CaloClusterPtr  sc, std::vector< reco::TrackRef >  conv, std::vector< float >  mva, reco::ConversionCollection &  oneLegConversions  )

private

Definition at line 159 of file GEDPhotonCoreProducer.cc.

References error, reco::Conversion::pflow, edm::PtrVector< T >::push_back(), reco::Conversion::setOneLegMVA(), X, and Gflash::Z.

Referenced by produce().

                                                                                                                                                                                {
  // this method translates the single track into the Conversion Data Format

  math::Error<3>::type error;
  for (unsigned int itk=0; itk<conv.size(); itk++){
    const reco::Vertex  * convVtx = new reco::Vertex(conv[itk]->innerPosition(), error);
    std::vector<reco::TrackRef> OneLegConvVector;
    OneLegConvVector.push_back(conv[itk]);
    std::vector< float > OneLegMvaVector;
    OneLegMvaVector.push_back(mva[itk]);
    std::vector<reco::CaloClusterPtr> dummymatchingBC;
    reco::CaloClusterPtrVector scPtrVec;
    scPtrVec.push_back(sc);


    std::vector<math::XYZPointF>trackPositionAtEcalVec;
    std::vector<math::XYZPointF>innPointVec;
    std::vector<math::XYZVectorF>trackPinVec;
    std::vector<math::XYZVectorF>trackPoutVec;
    math::XYZPointF trackPositionAtEcal(conv[itk]->outerPosition().X(), conv[itk]->outerPosition().Y(), conv[itk]->outerPosition().Z());
    trackPositionAtEcalVec.push_back(trackPositionAtEcal);

    math::XYZPointF innPoint(conv[itk]->innerPosition().X(), conv[itk]->innerPosition().Y(), conv[itk]->innerPosition().Z());
    innPointVec.push_back(innPoint);

    math::XYZVectorF trackPin(conv[itk]->innerMomentum().X(), conv[itk]->innerMomentum().Y(), conv[itk]->innerMomentum().Z());
    trackPinVec.push_back(trackPin);

    math::XYZVectorF trackPout(conv[itk]->outerMomentum().X(), conv[itk]->outerMomentum().Y(), conv[itk]->outerMomentum().Z());
    trackPoutVec.push_back( trackPout );

    float DCA = conv[itk]->d0() ;
    reco::Conversion singleLegConvCandidate(scPtrVec, 
                    OneLegConvVector,
                    trackPositionAtEcalVec,
                    *convVtx,
                    dummymatchingBC,
                    DCA,
                    innPointVec,
                    trackPinVec,
                    trackPoutVec,
                    mva[itk],             
                    reco::Conversion::pflow);
    singleLegConvCandidate.setOneLegMVA(OneLegMvaVector); 
    oneLegConversions.push_back(singleLegConvCandidate);

  } 



}

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

virtual

Implements edm::EDProducer.

Definition at line 59 of file GEDPhotonCoreProducer.cc.

References reco::PhotonCore::addConversion(), createSingleLegConversions(), reco::PFCandidate::gamma, GEDPhotonCoreCollection_, edm::Event::getByLabel(), edm::Ref< C, T, F >::id(), edm::Ref< C, T, F >::key(), edm::InputTag::label(), PFConversionCollection_, pfEgammaCandidates_, edm::Event::put(), edm::refToPtr(), reco::PhotonCore::setPFlowPhoton(), reco::PhotonCore::setPflowSuperCluster(), reco::PhotonCore::setStandardPhoton(), reco::PhotonCore::setSuperCluster(), and edm::RefVector< C, T, F >::size().

                                                                                          {


  using namespace edm;
  //  nEvt_++;

  reco::PhotonCoreCollection outputPhotonCoreCollection;
  std::auto_ptr< reco::PhotonCoreCollection > outputPhotonCoreCollection_p(new reco::PhotonCoreCollection);

  reco::ConversionCollection outputOneLegConversionCollection;
  std::auto_ptr<reco::ConversionCollection> SingleLeg_p(new reco::ConversionCollection(outputOneLegConversionCollection));  

  // Get the  PF refined cluster  collection
  Handle<reco::PFCandidateCollection> pfCandidateHandle;
  theEvent.getByLabel(pfEgammaCandidates_,pfCandidateHandle);
  if (!pfCandidateHandle.isValid()) {
    edm::LogError("GEDPhotonCoreProducer") << "Error! Can't get the product "<<pfEgammaCandidates_.label();
  }

 
  //  std::cout <<  "  GEDPhotonCoreProducer::produce input PFcandidate size " <<   pfCandidateHandle->size() << std::endl;


  // Loop over PF candidates and get only photons

  for(unsigned int lCand=0; lCand < pfCandidateHandle->size(); lCand++) {
    reco::PFCandidateRef candRef (reco::PFCandidateRef(pfCandidateHandle,lCand));
    if(candRef->particleId()!=reco::PFCandidate::gamma) continue;

    // Retrieve stuff from the pfPhoton
    reco::PFCandidateEGammaExtraRef pfPhoRef =  candRef->egammaExtraRef();
    reco::SuperClusterRef  refinedSC= pfPhoRef->superClusterRef();
    reco::SuperClusterRef  boxSC= pfPhoRef->superClusterBoxRef();
    const reco::ConversionRefVector & doubleLegConv = pfPhoRef->conversionRef();
    reco::CaloClusterPtr refinedSCPtr= edm::refToPtr(refinedSC);

    //    std::cout << "newCandidate  doubleLegConv="<<doubleLegConv.size()<< std::endl;
    //std::cout << "newCandidate  singleLegConv="<<  pfPhoRef->singleLegConvTrackRef().size()<< std::endl;

    reco::PhotonCore newCandidate;
    newCandidate.setPFlowPhoton(true);
    newCandidate.setStandardPhoton(false);
    newCandidate.setSuperCluster(refinedSC);
    newCandidate.setPflowSuperCluster(boxSC);
    // fill conversion infos
    

    for(unsigned int lConv=0; lConv < doubleLegConv.size(); lConv++) {
      newCandidate.addConversion(doubleLegConv[lConv]);
    } 

    //    std::cout << "newCandidate pf refined SC energy="<< newCandidate.superCluster()->energy()<<std::endl;
    //std::cout << "newCandidate pf SC energy="<< newCandidate.pfSuperCluster()->energy()<<std::endl;
    //std::cout << "newCandidate  nconv2leg="<<newCandidate.conversions().size()<< std::endl;



    //
    createSingleLegConversions(refinedSCPtr,  pfPhoRef->singleLegConvTrackRef(), pfPhoRef->singleLegConvMva(), outputOneLegConversionCollection); 
    outputPhotonCoreCollection.push_back(newCandidate);
  }


  const edm::OrphanHandle<reco::ConversionCollection> singleLegConvOrhpHandle = theEvent.put(SingleLeg_p,PFConversionCollection_);

  //  std::cout <<  "  GEDPhotonCoreProducer::produce orphanHandle to single legs " <<  singleLegConvOrhpHandle->size() << std::endl;
  //std::cout <<  "  GEDPhotonCoreProducer::produce photon size " <<  outputPhotonCoreCollection.size() << std::endl;

  

  int ipho=0;
  for (reco::PhotonCoreCollection::iterator gamIter=outputPhotonCoreCollection.begin(); gamIter != outputPhotonCoreCollection.end(); ++gamIter){

    for( unsigned int icp = 0;  icp <  singleLegConvOrhpHandle->size(); icp++) {
      const reco::ConversionRef cpRef(reco::ConversionRef(singleLegConvOrhpHandle,icp));
      if ( !cpRef->caloCluster().size()) continue; 
      if (!( gamIter->superCluster().id() == cpRef->caloCluster()[0].id() && gamIter->superCluster().key() == cpRef->caloCluster()[0].key() )) continue; 
      gamIter->addOneLegConversion(cpRef);     
    }
    // debug
    //    std::cout << "PhotonCoreCollection i="<<ipho<<" pf refined SC energy="<<gamIter->superCluster()->energy()<<std::endl;
    //std::cout << "PhotonCoreCollection i="<<ipho<<" pf SC energy="<<gamIter->pfSuperCluster()->energy()<<std::endl;
    //std::cout << "PhotonCoreCollection i="<<ipho<<" nconv2leg="<<gamIter->conversions().size()<<" nconv1leg="<<gamIter->conversionsOneLeg().size()<<std::endl;
    ipho++;
  }



  // put the product in the event
  //  edm::LogInfo("GEDPhotonCoreProducer") << " Put in the event " << iSC << " Photon Candidates \n";
  outputPhotonCoreCollection_p->assign(outputPhotonCoreCollection.begin(),outputPhotonCoreCollection.end());
  theEvent.put( outputPhotonCoreCollection_p, GEDPhotonCoreCollection_);
  


  
} 

Member Data Documentation

edm::ParameterSet GEDPhotonCoreProducer::conf_

private

Definition at line 45 of file GEDPhotonCoreProducer.h.

Referenced by GEDPhotonCoreProducer().

std::string GEDPhotonCoreProducer::GEDPhotonCoreCollection_

private

Definition at line 38 of file GEDPhotonCoreProducer.h.

Referenced by GEDPhotonCoreProducer(), and produce().

double GEDPhotonCoreProducer::minSCEt_

private

Definition at line 42 of file GEDPhotonCoreProducer.h.

std::string GEDPhotonCoreProducer::PFConversionCollection_

private

Definition at line 39 of file GEDPhotonCoreProducer.h.

Referenced by GEDPhotonCoreProducer(), and produce().

edm::InputTag GEDPhotonCoreProducer::pfEgammaCandidates_

private

Definition at line 40 of file GEDPhotonCoreProducer.h.

Referenced by GEDPhotonCoreProducer(), and produce().

std::string GEDPhotonCoreProducer::pixelSeedProducer_

private

Definition at line 44 of file GEDPhotonCoreProducer.h.

bool GEDPhotonCoreProducer::validConversions_

private

Definition at line 43 of file GEDPhotonCoreProducer.h.

bool GEDPhotonCoreProducer::validPixelSeeds_

private

Definition at line 46 of file GEDPhotonCoreProducer.h.