PFTauTransverseImpactParameters Class Reference

Inheritance diagram for PFTauTransverseImpactParameters:

Public Types
enum	Alg { useInputPV =0, useFont }
enum	CMSSWPerigee {   aCurv =0, aTheta, aPhi, aTip,   aLip }
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Member Functions
	PFTauTransverseImpactParameters (const edm::ParameterSet &iConfig)
virtual void	produce (edm::Event &, const edm::EventSetup &)
	~PFTauTransverseImpactParameters ()
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
	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
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Attributes
edm::EDGetTokenT < edm::AssociationVector < PFTauRefProd, std::vector < reco::VertexRef > > >	PFTauPVAToken_
edm::EDGetTokenT < edm::AssociationVector < PFTauRefProd, std::vector < std::vector< reco::VertexRef > > > >	PFTauSVAToken_
edm::EDGetTokenT< std::vector < reco::PFTau > >	PFTauToken_
bool	useFullCalculation_

Additional Inherited Members
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 53 of file PFTauTransverseImpactParameters.cc.

Member Enumeration Documentation

enum PFTauTransverseImpactParameters::Alg

Enumerator
useInputPV
useFont

Definition at line 55 of file PFTauTransverseImpactParameters.cc.

{useInputPV=0, useFont};

enum PFTauTransverseImpactParameters::CMSSWPerigee

Enumerator
aCurv
aTheta
aPhi
aTip
aLip

Definition at line 56 of file PFTauTransverseImpactParameters.cc.

{aCurv=0,aTheta,aPhi,aTip,aLip};

Constructor & Destructor Documentation

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

explicit

Definition at line 67 of file PFTauTransverseImpactParameters.cc.

                                                                                              :
  PFTauToken_(consumes<std::vector<reco::PFTau> >(iConfig.getParameter<edm::InputTag>("PFTauTag"))),
  PFTauPVAToken_(consumes<edm::AssociationVector<PFTauRefProd, std::vector<reco::VertexRef> > >(iConfig.getParameter<edm::InputTag>("PFTauPVATag"))),
  PFTauSVAToken_(consumes<edm::AssociationVector<PFTauRefProd,std::vector<std::vector<reco::VertexRef> > > >(iConfig.getParameter<edm::InputTag>("PFTauSVATag"))),
  useFullCalculation_(iConfig.getParameter<bool>("useFullCalculation"))
{
  produces<edm::AssociationVector<PFTauRefProd, std::vector<reco::PFTauTransverseImpactParameterRef> > >(); 
  produces<PFTauTransverseImpactParameterCollection>("PFTauTIP");
}

PFTauTransverseImpactParameters::~PFTauTransverseImpactParameters

(

)

Definition at line 77 of file PFTauTransverseImpactParameters.cc.

                                                                 {

}

Member Function Documentation

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

virtual

Implements edm::EDProducer.

Definition at line 81 of file PFTauTransverseImpactParameters.cc.

References aTip, PerigeeTrajectoryError::covarianceMatrix(), reco::Vertex::dimension, edm::EventSetup::get(), edm::Event::getByToken(), edm::Event::getRefBeforePut(), i, edm::Ref< C, T, F >::isNonnull(), edm::HandleBase::isValid(), j, edm::Ref< C, T, F >::key(), TrajectoryStateClosestToPoint::perigeeError(), TrajectoryStateClosestToPoint::perigeeParameters(), PFTauPVAToken_, PFTauSVAToken_, PFTauToken_, TrajectoryStateClosestToPoint::position(), edm::Event::put(), MetAnalyzer::pv(), reco::TransientTrack::trajectoryStateClosestToPoint(), useFullCalculation_, findQualityFiles::v, PerigeeTrajectoryParameters::vector(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                         {
  // Obtain Collections
  edm::ESHandle<TransientTrackBuilder> transTrackBuilder;
  iSetup.get<TransientTrackRecord>().get("TransientTrackBuilder",transTrackBuilder);
  
  edm::Handle<std::vector<reco::PFTau> > Tau;
  iEvent.getByToken(PFTauToken_,Tau);

  edm::Handle<edm::AssociationVector<PFTauRefProd, std::vector<reco::VertexRef> > > PFTauPVA;
  iEvent.getByToken(PFTauPVAToken_,PFTauPVA);

  edm::Handle<edm::AssociationVector<PFTauRefProd,std::vector<std::vector<reco::VertexRef> > > > PFTauSVA;
  iEvent.getByToken(PFTauSVAToken_,PFTauSVA);

  // Set Association Map
  auto_ptr<edm::AssociationVector<PFTauRefProd, std::vector<reco::PFTauTransverseImpactParameterRef> > > AVPFTauTIP(new edm::AssociationVector<PFTauRefProd, std::vector<reco::PFTauTransverseImpactParameterRef> >(PFTauRefProd(Tau)));
  std::auto_ptr<PFTauTransverseImpactParameterCollection>  TIPCollection_out= std::auto_ptr<PFTauTransverseImpactParameterCollection>(new PFTauTransverseImpactParameterCollection());
  reco::PFTauTransverseImpactParameterRefProd TIPRefProd_out = iEvent.getRefBeforePut<reco::PFTauTransverseImpactParameterCollection>("PFTauTIP");


  // For each Tau Run Algorithim
  if(Tau.isValid()) {
    for(reco::PFTauCollection::size_type iPFTau = 0; iPFTau < Tau->size(); iPFTau++) {
      reco::PFTauRef RefPFTau(Tau, iPFTau);
      const reco::VertexRef PV=PFTauPVA->value(RefPFTau.key());
      const std::vector<reco::VertexRef> SV=PFTauSVA->value(RefPFTau.key());
      double dxy(-999), dxy_err(-999);
      reco::Vertex::Point poca(0,0,0);
      if(RefPFTau->leadPFChargedHadrCand().isNonnull()){
    if(RefPFTau->leadPFChargedHadrCand()->trackRef().isNonnull()){
      if(useFullCalculation_){
        reco::TransientTrack transTrk=transTrackBuilder->build(RefPFTau->leadPFChargedHadrCand()->trackRef());
        GlobalPoint pv(PV->position().x(),PV->position().y(),PV->position().z());
        dxy=-transTrk.trajectoryStateClosestToPoint(pv).perigeeParameters().vector()(aTip);
        dxy_err=transTrk.trajectoryStateClosestToPoint(pv).perigeeError().covarianceMatrix()(aTip,aTip);
        GlobalPoint pos=transTrk.trajectoryStateClosestToPoint(pv).position();
        poca=reco::Vertex::Point(pos.x(),pos.y(),pos.z());
      }
      else{
        dxy_err=RefPFTau->leadPFChargedHadrCand()->trackRef()->d0Error();
        dxy=RefPFTau->leadPFChargedHadrCand()->trackRef()->dxy(PV->position());
      }
    }
      }
      if(SV.size()>0){
    reco::Vertex::CovarianceMatrix cov;
    reco::Vertex::Point v(SV.at(0)->x()-PV->x(),SV.at(0)->y()-PV->y(),SV.at(0)->z()-PV->z());
    for(int i=0;i<reco::Vertex::dimension;i++){
      for(int j=0;j<reco::Vertex::dimension;j++){
        cov(i,j)=SV.at(0)->covariance(i,j)+PV->covariance(i,j);
      }
    }
    GlobalVector direction(RefPFTau->px(),RefPFTau->py(),RefPFTau->pz());
    double vSig = SecondaryVertex::computeDist3d(*PV,*SV.at(0),direction,true).significance();
    reco::PFTauTransverseImpactParameter TIPV(poca,dxy,dxy_err,PV,v,vSig,SV.at(0));
    reco::PFTauTransverseImpactParameterRef TIPVRef=reco::PFTauTransverseImpactParameterRef(TIPRefProd_out,TIPCollection_out->size());
        TIPCollection_out->push_back(TIPV);
        AVPFTauTIP->setValue(iPFTau,TIPVRef);
      }
      else{ 
    reco::PFTauTransverseImpactParameter TIPV(poca,dxy,dxy_err,PV);
    reco::PFTauTransverseImpactParameterRef TIPVRef=reco::PFTauTransverseImpactParameterRef(TIPRefProd_out,TIPCollection_out->size());
    TIPCollection_out->push_back(TIPV);
    AVPFTauTIP->setValue(iPFTau,TIPVRef);
      }
    }
  }
  iEvent.put(TIPCollection_out,"PFTauTIP");
  iEvent.put(AVPFTauTIP);
}

Member Data Documentation

edm::EDGetTokenT<edm::AssociationVector<PFTauRefProd, std::vector<reco::VertexRef> > > PFTauTransverseImpactParameters::PFTauPVAToken_

private

Definition at line 62 of file PFTauTransverseImpactParameters.cc.

Referenced by produce().

edm::EDGetTokenT<edm::AssociationVector<PFTauRefProd,std::vector<std::vector<reco::VertexRef> > > > PFTauTransverseImpactParameters::PFTauSVAToken_

private

Definition at line 63 of file PFTauTransverseImpactParameters.cc.

Referenced by produce().

edm::EDGetTokenT<std::vector<reco::PFTau> > PFTauTransverseImpactParameters::PFTauToken_

private

Definition at line 61 of file PFTauTransverseImpactParameters.cc.

Referenced by produce().

bool PFTauTransverseImpactParameters::useFullCalculation_

private

Definition at line 64 of file PFTauTransverseImpactParameters.cc.

Referenced by produce().