CMS 3D CMS Logo
List of all members | Public Types | Public Member Functions | Private Attributes
PFTauTransverseImpactParameters Class Reference
Inheritance diagram for PFTauTransverseImpactParameters:
edm::stream::EDProducer<>

Public Types

enum  CMSSWPerigee {
  aCurv =0, aTheta, aPhi, aTip,
  aLip
}
 
- Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... > CacheTypes
 
typedef CacheTypes::GlobalCache GlobalCache
 
typedef AbilityChecker< T... > HasAbility
 
typedef CacheTypes::LuminosityBlockCache LuminosityBlockCache
 
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache > LuminosityBlockContext
 
typedef CacheTypes::LuminosityBlockSummaryCache LuminosityBlockSummaryCache
 
typedef CacheTypes::RunCache RunCache
 
typedef RunContextT< RunCache, GlobalCache > RunContext
 
typedef CacheTypes::RunSummaryCache RunSummaryCache
 

Public Member Functions

 PFTauTransverseImpactParameters (const edm::ParameterSet &iConfig)
 
void produce (edm::Event &, const edm::EventSetup &) override
 
 ~PFTauTransverseImpactParameters () override
 
- Public Member Functions inherited from edm::stream::EDProducer<>
 EDProducer ()=default
 
bool hasAbilityToProduceInLumis () const final
 
bool hasAbilityToProduceInRuns () const final
 

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_
 

Detailed Description

Definition at line 51 of file PFTauTransverseImpactParameters.cc.

Member Enumeration Documentation

enum PFTauTransverseImpactParameters::CMSSWPerigee
Enumerator
aCurv 
aTheta 
aPhi 
aTip 
aLip 

Definition at line 53 of file PFTauTransverseImpactParameters.cc.

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

Constructor & Destructor Documentation

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

Definition at line 64 of file PFTauTransverseImpactParameters.cc.

64  :
65  PFTauToken_(consumes<std::vector<reco::PFTau> >(iConfig.getParameter<edm::InputTag>("PFTauTag"))),
66  PFTauPVAToken_(consumes<edm::AssociationVector<PFTauRefProd, std::vector<reco::VertexRef> > >(iConfig.getParameter<edm::InputTag>("PFTauPVATag"))),
67  PFTauSVAToken_(consumes<edm::AssociationVector<PFTauRefProd,std::vector<std::vector<reco::VertexRef> > > >(iConfig.getParameter<edm::InputTag>("PFTauSVATag"))),
68  useFullCalculation_(iConfig.getParameter<bool>("useFullCalculation"))
69 {
70  produces<edm::AssociationVector<PFTauRefProd, std::vector<reco::PFTauTransverseImpactParameterRef> > >();
71  produces<PFTauTransverseImpactParameterCollection>("PFTauTIP");
72 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
PFTauTransverseImpactParameters::PFTauToken_
edm::EDGetTokenT< std::vector< reco::PFTau > > PFTauToken_
Definition: PFTauTransverseImpactParameters.cc:58
PFTauTransverseImpactParameters::PFTauSVAToken_
edm::EDGetTokenT< edm::AssociationVector< PFTauRefProd, std::vector< std::vector< reco::VertexRef > > > > PFTauSVAToken_
Definition: PFTauTransverseImpactParameters.cc:60
PFTauTransverseImpactParameters::PFTauPVAToken_
edm::EDGetTokenT< edm::AssociationVector< PFTauRefProd, std::vector< reco::VertexRef > > > PFTauPVAToken_
Definition: PFTauTransverseImpactParameters.cc:59
PFTauTransverseImpactParameters::~PFTauTransverseImpactParameters ( )
override

Definition at line 74 of file PFTauTransverseImpactParameters.cc.

74  {
75 
76 }

Member Function Documentation

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

Definition at line 78 of file PFTauTransverseImpactParameters.cc.

References TransientTrackBuilder::build(), RecoVertex::convertPos(), reco::TrackBase::d0Error(), DEFINE_FWK_MODULE, reco::Vertex::dimension, PVValHelper::dxy, reco::TrackBase::dxy(), reco::TrackBase::dz(), reco::TrackBase::dzError(), reco::TransientTrack::field(), edm::EventSetup::get(), edm::Ref< C, T, F >::get(), edm::Event::getByToken(), edm::Event::getRefBeforePut(), mps_fire::i, reco::TransientTrack::impactPointState(), electrons_cff::ip3d, edm::Ref< C, T, F >::isNonnull(), edm::HandleBase::isValid(), edm::Ref< C, T, F >::key(), eostools::move(), PFTauPVAToken_, PFTauSVAToken_, PFTauToken_, edm::Event::put(), nanoDQM_cfi::PV, IPTools::signedImpactParameter3D(), IPTools::signedTransverseImpactParameter(), met_cff::significance, nanoDQM_cfi::SV, HiIsolationCommonParameters_cff::track, useFullCalculation_, and findQualityFiles::v.

78  {
79  // Obtain Collections
80  edm::ESHandle<TransientTrackBuilder> transTrackBuilder;
81  iSetup.get<TransientTrackRecord>().get("TransientTrackBuilder",transTrackBuilder);
82 
83  edm::Handle<std::vector<reco::PFTau> > Tau;
84  iEvent.getByToken(PFTauToken_,Tau);
85 
86  edm::Handle<edm::AssociationVector<PFTauRefProd, std::vector<reco::VertexRef> > > PFTauPVA;
87  iEvent.getByToken(PFTauPVAToken_,PFTauPVA);
88 
89  edm::Handle<edm::AssociationVector<PFTauRefProd,std::vector<std::vector<reco::VertexRef> > > > PFTauSVA;
90  iEvent.getByToken(PFTauSVAToken_,PFTauSVA);
91 
92  // Set Association Map
93  auto AVPFTauTIP = std::make_unique< edm::AssociationVector<PFTauRefProd, std::vector<reco::PFTauTransverseImpactParameterRef>>>(PFTauRefProd(Tau));
94  auto TIPCollection_out = std::make_unique<PFTauTransverseImpactParameterCollection>();
95  reco::PFTauTransverseImpactParameterRefProd TIPRefProd_out = iEvent.getRefBeforePut<reco::PFTauTransverseImpactParameterCollection>("PFTauTIP");
96 
97  // For each Tau Run Algorithim
98  if(Tau.isValid()) {
99  for(reco::PFTauCollection::size_type iPFTau = 0; iPFTau < Tau->size(); iPFTau++) {
100  reco::PFTauRef RefPFTau(Tau, iPFTau);
101  const reco::VertexRef PV=PFTauPVA->value(RefPFTau.key());
102  const std::vector<reco::VertexRef> SV=PFTauSVA->value(RefPFTau.key());
103  double dxy(-999), dxy_err(-999);
104  reco::Vertex::Point poca(0,0,0);
105  double ip3d(-999), ip3d_err(-999);
106  reco::Vertex::Point ip3d_poca(0,0,0);
107  if(RefPFTau->leadPFChargedHadrCand().isNonnull()){
108  const reco::Track* track = nullptr;
109  if(RefPFTau->leadPFChargedHadrCand()->trackRef().isNonnull())
110  track = RefPFTau->leadPFChargedHadrCand()->trackRef().get();
111  else if(RefPFTau->leadPFChargedHadrCand()->gsfTrackRef().isNonnull())
112  track = RefPFTau->leadPFChargedHadrCand()->gsfTrackRef().get();
113  if(track != nullptr){
114  if(useFullCalculation_){
115  reco::TransientTrack transTrk=transTrackBuilder->build(*track);
116  GlobalVector direction(RefPFTau->p4().px(), RefPFTau->p4().py(), RefPFTau->p4().pz()); //To compute sign of IP
117  std::pair<bool,Measurement1D> signed_IP2D = IPTools::signedTransverseImpactParameter(transTrk, direction, (*PV));
118  dxy=signed_IP2D.second.value();
119  dxy_err=signed_IP2D.second.error();
120  std::pair<bool,Measurement1D> signed_IP3D = IPTools::signedImpactParameter3D(transTrk, direction, (*PV));
121  ip3d=signed_IP3D.second.value();
122  ip3d_err=signed_IP3D.second.error();
123  TransverseImpactPointExtrapolator extrapolator(transTrk.field());
124  GlobalPoint pos = extrapolator.extrapolate(transTrk.impactPointState(), RecoVertex::convertPos(PV->position())).globalPosition();
125  poca=reco::Vertex::Point(pos.x(),pos.y(),pos.z());
126  AnalyticalImpactPointExtrapolator extrapolator3D(transTrk.field());
127  GlobalPoint pos3d = extrapolator3D.extrapolate(transTrk.impactPointState(),RecoVertex::convertPos(PV->position())).globalPosition();
128  ip3d_poca=reco::Vertex::Point(pos3d.x(),pos3d.y(),pos3d.z());
129  }
130  else{
131  dxy_err=track->d0Error();
132  dxy=track->dxy(PV->position());
133  ip3d_err=track->dzError(); //store dz, ip3d not available
134  ip3d=track->dz(PV->position()); //store dz, ip3d not available
135  }
136  }
137  }
138  if(!SV.empty()){
139  reco::Vertex::CovarianceMatrix cov;
140  reco::Vertex::Point v(SV.at(0)->x()-PV->x(),SV.at(0)->y()-PV->y(),SV.at(0)->z()-PV->z());
141  for(int i=0;i<reco::Vertex::dimension;i++){
142  for(int j=0;j<reco::Vertex::dimension;j++){
143  cov(i,j)=SV.at(0)->covariance(i,j)+PV->covariance(i,j);
144  }
145  }
146  GlobalVector direction(RefPFTau->px(),RefPFTau->py(),RefPFTau->pz());
147  double vSig = SecondaryVertex::computeDist3d(*PV,*SV.at(0),direction,true).significance();
148  reco::PFTauTransverseImpactParameter TIPV(poca,dxy,dxy_err,ip3d_poca,ip3d,ip3d_err,PV,v,vSig,SV.at(0));
149  reco::PFTauTransverseImpactParameterRef TIPVRef=reco::PFTauTransverseImpactParameterRef(TIPRefProd_out,TIPCollection_out->size());
150  TIPCollection_out->push_back(TIPV);
151  AVPFTauTIP->setValue(iPFTau,TIPVRef);
152  }
153  else{
154  reco::PFTauTransverseImpactParameter TIPV(poca,dxy,dxy_err,ip3d_poca,ip3d,ip3d_err,PV);
155  reco::PFTauTransverseImpactParameterRef TIPVRef=reco::PFTauTransverseImpactParameterRef(TIPRefProd_out,TIPCollection_out->size());
156  TIPCollection_out->push_back(TIPV);
157  AVPFTauTIP->setValue(iPFTau,TIPVRef);
158  }
159  }
160  }
161  iEvent.put(std::move(TIPCollection_out),"PFTauTIP");
162  iEvent.put(std::move(AVPFTauTIP));
163 }
RecoVertex::convertPos
reco::Vertex::Point convertPos(const GlobalPoint &p)
Definition: ConvertToFromReco.h:7
reco::TrackBase::d0Error
double d0Error() const
error on d0
Definition: TrackBase.h:802
edm::Event::put
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:137
mps_fire.i
i
Definition: mps_fire.py:277
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:579
nanoDQM_cfi.PV
PV
Definition: nanoDQM_cfi.py:319
IPTools::signedTransverseImpactParameter
std::pair< bool, Measurement1D > signedTransverseImpactParameter(const reco::TransientTrack &track, const GlobalVector &direction, const reco::Vertex &vertex)
Definition: IPTools.cc:50
TransientTrackBuilder::build
reco::TransientTrack build(const reco::Track *p) const
Definition: TransientTrackBuilder.cc:20
IPTools::signedImpactParameter3D
std::pair< bool, Measurement1D > signedImpactParameter3D(const reco::TransientTrack &track, const GlobalVector &direction, const reco::Vertex &vertex)
Definition: IPTools.cc:71
reco::TransientTrack::field
const MagneticField * field() const
Definition: TransientTrack.h:113
reco::PFTauRefProd
edm::RefProd< PFTauCollection > PFTauRefProd
references to PFTau collection
Definition: PFTauFwd.h:15
electrons_cff.ip3d
ip3d
Definition: electrons_cff.py:285
trigger::size_type
uint16_t size_type
Definition: TriggerTypeDefs.h:19
reco::Vertex::CovarianceMatrix
math::Error< dimension >::type CovarianceMatrix
covariance error matrix (3x3)
Definition: Vertex.h:45
PFTauTransverseImpactParameters::PFTauToken_
edm::EDGetTokenT< std::vector< reco::PFTau > > PFTauToken_
Definition: PFTauTransverseImpactParameters.cc:58
PFTauTransverseImpactParameters::PFTauSVAToken_
edm::EDGetTokenT< edm::AssociationVector< PFTauRefProd, std::vector< std::vector< reco::VertexRef > > > > PFTauSVAToken_
Definition: PFTauTransverseImpactParameters.cc:60
nanoDQM_cfi.SV
SV
Definition: nanoDQM_cfi.py:390
PFTauTransverseImpactParameters::PFTauPVAToken_
edm::EDGetTokenT< edm::AssociationVector< PFTauRefProd, std::vector< reco::VertexRef > > > PFTauPVAToken_
Definition: PFTauTransverseImpactParameters.cc:59
reco::Vertex::Point
math::XYZPoint Point
point in the space
Definition: Vertex.h:39
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:74
edm::Event::getRefBeforePut
RefProd< PROD > getRefBeforePut()
Definition: Event.h:167
reco::TrackBase::dz
double dz() const
dz parameter (= dsz/cos(lambda)). This is the track z0 w.r.t (0,0,0) only if the refPoint is close to...
Definition: TrackBase.h:609
reco::TrackBase::dzError
double dzError() const
error on dz
Definition: TrackBase.h:814
HiIsolationCommonParameters_cff.track
track
Definition: HiIsolationCommonParameters_cff.py:9
reco::PFTauTransverseImpactParameterRef
edm::Ref< PFTauTransverseImpactParameterCollection > PFTauTransverseImpactParameterRef
presistent reference to a PFTauTransverseImpactParameter
Definition: PFTauTransverseImpactParameterFwd.h:13
met_cff.significance
significance
Definition: met_cff.py:30
edm::EventSetup::get
T get() const
Definition: EventSetup.h:62
reco::TransientTrack::impactPointState
TrajectoryStateOnSurface impactPointState() const
Definition: TransientTrack.h:97
reco::PFTauTransverseImpactParameterCollection
std::vector< reco::PFTauTransverseImpactParameter > PFTauTransverseImpactParameterCollection
collection of PFTauTransverseImpactParameter objects
Definition: PFTauTransverseImpactParameterFwd.h:9
reco::TrackBase::dxy
double dxy() const
dxy parameter. (This is the transverse impact parameter w.r.t. to (0,0,0) ONLY if refPoint is close t...
Definition: TrackBase.h:591
eostools.move
def move(src, dest)
Definition: eostools.py:511

Member Data Documentation

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

Definition at line 59 of file PFTauTransverseImpactParameters.cc.

Referenced by produce().

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

Definition at line 60 of file PFTauTransverseImpactParameters.cc.

Referenced by produce().

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

Definition at line 58 of file PFTauTransverseImpactParameters.cc.

Referenced by produce().

bool PFTauTransverseImpactParameters::useFullCalculation_
private

Definition at line 61 of file PFTauTransverseImpactParameters.cc.

Referenced by produce().

Generated for CMSSW Reference Manual by doxygen 1.8.11