|
|
|
Abstract. More...
#include <PFElecTkProducer.h>
Abstract.
PFElecTkProducer reads the merged GsfTracks collection built with the TrackerDriven and EcalDriven seeds and transform them in PFGsfRecTracks.
Definition at line 35 of file PFElecTkProducer.h.
| PFElecTkProducer::PFElecTkProducer | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Constructor.
Definition at line 50 of file PFElecTkProducer.cc.
References applyAngularGsfClean_, applyGsfClean_, applySel_, debugGsfCleaning_, detaCutGsfClean_, detaGsfSC_, dphiCutGsfClean_, dphiGsfSC_, edm::ParameterSet::getParameter(), gsfTrackLabel_, maxPtConvReco_, modemomentum_, mvaConvBremFinderID_, path_mvaWeightFileConvBrem_, pfConv_, pfEcalClusters_, pfNuclear_, pfTrackLabel_, pfV0_, primVtxLabel_, SCEne_, trajinev_, useConvBremFinder_, useConversions_, useFifthStepForEcalDriven_, useFifthStepForTrackDriven_, useNuclear_, and useV0_.
: conf_(iConfig), pfTransformer_(0), convBremFinder_(0) { LogInfo("PFElecTkProducer")<<"PFElecTkProducer started"; gsfTrackLabel_ = iConfig.getParameter<InputTag> ("GsfTrackModuleLabel"); pfTrackLabel_ = iConfig.getParameter<InputTag> ("PFRecTrackLabel"); primVtxLabel_ = iConfig.getParameter<InputTag> ("PrimaryVertexLabel"); pfEcalClusters_ = iConfig.getParameter<InputTag> ("PFEcalClusters"); pfNuclear_ = iConfig.getParameter<InputTag> ("PFNuclear"); pfConv_ = iConfig.getParameter<InputTag> ("PFConversions"); pfV0_ = iConfig.getParameter<InputTag> ("PFV0"); useNuclear_ = iConfig.getParameter<bool>("useNuclear"); useConversions_ = iConfig.getParameter<bool>("useConversions"); useV0_ = iConfig.getParameter<bool>("useV0"); debugGsfCleaning_ = iConfig.getParameter<bool>("debugGsfCleaning"); produces<GsfPFRecTrackCollection>(); produces<GsfPFRecTrackCollection>( "Secondary" ).setBranchAlias( "secondary" ); trajinev_ = iConfig.getParameter<bool>("TrajInEvents"); modemomentum_ = iConfig.getParameter<bool>("ModeMomentum"); applySel_ = iConfig.getParameter<bool>("applyEGSelection"); applyGsfClean_ = iConfig.getParameter<bool>("applyGsfTrackCleaning"); applyAngularGsfClean_ = iConfig.getParameter<bool>("applyAlsoGsfAngularCleaning"); detaCutGsfClean_ = iConfig.getParameter<double>("maxDEtaGsfAngularCleaning"); dphiCutGsfClean_ = iConfig.getParameter<double>("maxDPhiBremTangGsfAngularCleaning"); useFifthStepForTrackDriven_ = iConfig.getParameter<bool>("useFifthStepForTrackerDrivenGsf"); useFifthStepForEcalDriven_ = iConfig.getParameter<bool>("useFifthStepForEcalDrivenGsf"); maxPtConvReco_ = iConfig.getParameter<double>("MaxConvBremRecoPT"); detaGsfSC_ = iConfig.getParameter<double>("MinDEtaGsfSC"); dphiGsfSC_ = iConfig.getParameter<double>("MinDPhiGsfSC"); SCEne_ = iConfig.getParameter<double>("MinSCEnergy"); // set parameter for convBremFinder useConvBremFinder_ = iConfig.getParameter<bool>("useConvBremFinder"); mvaConvBremFinderID_ = iConfig.getParameter<double>("pf_convBremFinderID_mvaCut"); string mvaWeightFileConvBrem = iConfig.getParameter<string>("pf_convBremFinderID_mvaWeightFile"); if(useConvBremFinder_) path_mvaWeightFileConvBrem_ = edm::FileInPath ( mvaWeightFileConvBrem.c_str() ).fullPath(); }
| PFElecTkProducer::~PFElecTkProducer | ( | ) |
Destructor.
Definition at line 116 of file PFElecTkProducer.cc.
References convBremFinder_, and pfTransformer_.
{
delete pfTransformer_;
delete convBremFinder_;
}
| bool PFElecTkProducer::applySelection | ( | reco::GsfTrack | gsftk | ) | [private] |
Definition at line 553 of file PFElecTkProducer.cc.
References detaGsfSC_, dphiGsfSC_, reco::GsfTrack::etaMode(), reco::Track::extra(), edm::Ref< C, T, F >::isNonnull(), edm::Ref< C, T, F >::isNull(), reco::GsfTrack::phiMode(), Pi, SCEne_, reco::Track::seedRef(), and TwoPi.
Referenced by produce().
{
if (&(*gsftk.seedRef())==0) return false;
ElectronSeedRef ElSeedRef=gsftk.extra()->seedRef().castTo<ElectronSeedRef>();
bool passCut = false;
if (ElSeedRef->ctfTrack().isNull()){
if(ElSeedRef->caloCluster().isNull()) return passCut;
SuperClusterRef scRef = ElSeedRef->caloCluster().castTo<SuperClusterRef>();
//do this just to know if exist a SC?
if(scRef.isNonnull()) {
float caloEne = scRef->energy();
float feta = fabs(scRef->eta()-gsftk.etaMode());
float fphi = fabs(scRef->phi()-gsftk.phiMode());
if (fphi>TMath::Pi()) fphi-= TMath::TwoPi();
if(caloEne > SCEne_ && feta < detaGsfSC_ && fabs(fphi) < dphiGsfSC_)
passCut = true;
}
}
else {
// get all the gsf found by tracker driven
passCut = true;
}
return passCut;
}
| void PFElecTkProducer::beginRun | ( | edm::Run & | run, |
| const edm::EventSetup & | iSetup | ||
| ) | [private, virtual] |
Reimplemented from edm::EDProducer.
Definition at line 1175 of file PFElecTkProducer.cc.
References convBremFinder_, edm::EventSetup::get(), mtsTransform_, mvaConvBremFinderID_, path_mvaWeightFileConvBrem_, pfTransformer_, edm::ESHandle< T >::product(), patCandidatesForDimuonsSequences_cff::tracker, and useConvBremFinder_.
{
ESHandle<MagneticField> magneticField;
iSetup.get<IdealMagneticFieldRecord>().get(magneticField);
ESHandle<TrackerGeometry> tracker;
iSetup.get<TrackerDigiGeometryRecord>().get(tracker);
mtsTransform_ = MultiTrajectoryStateTransform(tracker.product(),magneticField.product());
pfTransformer_= new PFTrackTransformer(math::XYZVector(magneticField->inTesla(GlobalPoint(0,0,0))));
edm::ESHandle<TransientTrackBuilder> builder;
iSetup.get<TransientTrackRecord>().get("TransientTrackBuilder", builder);
TransientTrackBuilder thebuilder = *(builder.product());
if(useConvBremFinder_) {
FILE * fileConvBremID = fopen(path_mvaWeightFileConvBrem_.c_str(), "r");
if (fileConvBremID) {
fclose(fileConvBremID);
}
else {
string err = "PFElecTkProducer: cannot open weight file '";
err += path_mvaWeightFileConvBrem_;
err += "'";
throw invalid_argument( err );
}
}
convBremFinder_ = new ConvBremPFTrackFinder(thebuilder,mvaConvBremFinderID_,path_mvaWeightFileConvBrem_);
}
| void PFElecTkProducer::createGsfPFRecTrackRef | ( | const edm::OrphanHandle< reco::GsfPFRecTrackCollection > & | gsfPfHandle, |
| std::vector< reco::GsfPFRecTrack > & | gsfPFRecTrackPrimary, | ||
| const std::map< unsigned int, std::vector< reco::GsfPFRecTrack > > & | MapPrimSec | ||
| ) | [private] |
Definition at line 413 of file PFElecTkProducer.cc.
References Association::map.
Referenced by produce().
{
unsigned int cgsf=0;
unsigned int csecgsf=0;
for (std::map<unsigned int, std::vector<reco::GsfPFRecTrack> >::const_iterator igsf = MapPrimSec.begin();
igsf != MapPrimSec.end(); igsf++,cgsf++) {
vector<reco::GsfPFRecTrack> SecGsfPF = igsf->second;
for (unsigned int iSecGsf=0; iSecGsf < SecGsfPF.size(); iSecGsf++) {
edm::Ref<reco::GsfPFRecTrackCollection> refgprt(gsfPfHandle,csecgsf);
gsfPFRecTrackPrimary[cgsf].addConvBremGsfPFRecTrackRef(refgprt);
++csecgsf;
}
}
return;
}
| void PFElecTkProducer::endRun | ( | void | ) | [private, virtual] |
Definition at line 1214 of file PFElecTkProducer.cc.
References pfTransformer_.
{
delete pfTransformer_;
}
| int PFElecTkProducer::FindPfRef | ( | const reco::PFRecTrackCollection & | PfRTkColl, |
| reco::GsfTrack | gsftk, | ||
| bool | otherColl | ||
| ) | [private] |
Definition at line 432 of file PFElecTkProducer.cc.
References TrackingRecHit::all, reco::TrackBase::eta(), reco::Track::extra(), edm::Ref< C, T, F >::isNull(), reco::TrackBase::phi(), Pi, reco::Track::seedRef(), mathSSE::sqrt(), and TwoPi.
Referenced by produce().
{
if (&(*gsftk.seedRef())==0) return -1;
ElectronSeedRef ElSeedRef=gsftk.extra()->seedRef().castTo<ElectronSeedRef>();
//CASE 1 ELECTRONSEED DOES NOT HAVE A REF TO THE CKFTRACK
if (ElSeedRef->ctfTrack().isNull()){
reco::PFRecTrackCollection::const_iterator pft=PfRTkColl.begin();
reco::PFRecTrackCollection::const_iterator pftend=PfRTkColl.end();
unsigned int i_pf=0;
int ibest=-1;
unsigned int ish_max=0;
float dr_min=1000;
//SEARCH THE PFRECTRACK THAT SHARES HITS WITH THE ELECTRON SEED
// Here the cpu time can be improved.
for(;pft!=pftend;++pft){
unsigned int ish=0;
float dph= fabs(pft->trackRef()->phi()-gsftk.phi());
if (dph>TMath::Pi()) dph-= TMath::TwoPi();
float det=fabs(pft->trackRef()->eta()-gsftk.eta());
float dr =sqrt(dph*dph+det*det);
trackingRecHit_iterator hhit=
pft->trackRef()->recHitsBegin();
trackingRecHit_iterator hhit_end=
pft->trackRef()->recHitsEnd();
for(;hhit!=hhit_end;++hhit){
if (!(*hhit)->isValid()) continue;
TrajectorySeed::const_iterator hit=
gsftk.seedRef()->recHits().first;
TrajectorySeed::const_iterator hit_end=
gsftk.seedRef()->recHits().second;
for(;hit!=hit_end;++hit){
if (!(hit->isValid())) continue;
if((*hhit)->sharesInput(&*(hit),TrackingRecHit::all)) ish++;
// if((hit->geographicalId()==(*hhit)->geographicalId())&&
// (((*hhit)->localPosition()-hit->localPosition()).mag()<0.01)) ish++;
}
}
if ((ish>ish_max)||
((ish==ish_max)&&(dr<dr_min))){
ish_max=ish;
dr_min=dr;
ibest=i_pf;
}
i_pf++;
}
if (ibest<0) return -1;
if((ish_max==0) || (dr_min>0.05))return -1;
if(otherColl && (ish_max==0)) return -1;
return ibest;
}
else{
//ELECTRON SEED HAS A REFERENCE
reco::PFRecTrackCollection::const_iterator pft=PfRTkColl.begin();
reco::PFRecTrackCollection::const_iterator pftend=PfRTkColl.end();
unsigned int i_pf=0;
for(;pft!=pftend;++pft){
//REF COMPARISON
if (pft->trackRef()==ElSeedRef->ctfTrack()){
return i_pf;
}
i_pf++;
}
}
return -1;
}
| bool PFElecTkProducer::isFifthStep | ( | reco::PFRecTrackRef | pfKfTrack | ) | [private] |
Definition at line 514 of file PFElecTkProducer.cc.
Referenced by produce().
{
bool isFithStep = false;
TrackRef kfref = pfKfTrack->trackRef();
unsigned int Algo = 0;
switch (kfref->algo()) {
case TrackBase::undefAlgorithm:
case TrackBase::ctf:
case TrackBase::iter0:
case TrackBase::iter1:
case TrackBase::iter2:
Algo = 0;
break;
case TrackBase::iter3:
Algo = 1;
break;
case TrackBase::iter4:
Algo = 2;
break;
case TrackBase::iter5:
Algo = 3;
break;
case TrackBase::iter6:
Algo = 4;
break;
default:
Algo = 5;
break;
}
if ( Algo >= 4 ) {
isFithStep = true;
}
return isFithStep;
}
| bool PFElecTkProducer::isInnerMost | ( | const reco::GsfTrackRef & | nGsfTrack, |
| const reco::GsfTrackRef & | iGsfTrack, | ||
| bool & | sameLayer | ||
| ) | [private] |
Definition at line 1109 of file PFElecTkProducer.cc.
References reco::HitPattern::getHitPattern(), reco::HitPattern::getLayer(), and reco::HitPattern::getSubStructure().
{
// copied by the class RecoEgamma/EgammaElectronAlgos/src/EgAmbiguityTools.cc
// obsolete but the code is kept: now using lost hits method
reco::HitPattern gsfHitPattern1 = nGsfTrack->hitPattern();
reco::HitPattern gsfHitPattern2 = iGsfTrack->hitPattern();
// retrieve first valid hit
int gsfHitCounter1 = 0 ;
trackingRecHit_iterator elHitsIt1 ;
for
( elHitsIt1 = nGsfTrack->recHitsBegin() ;
elHitsIt1 != nGsfTrack->recHitsEnd() ;
elHitsIt1++, gsfHitCounter1++ )
{ if (((**elHitsIt1).isValid())) break ; }
int gsfHitCounter2 = 0 ;
trackingRecHit_iterator elHitsIt2 ;
for
( elHitsIt2 = iGsfTrack->recHitsBegin() ;
elHitsIt2 != iGsfTrack->recHitsEnd() ;
elHitsIt2++, gsfHitCounter2++ )
{ if (((**elHitsIt2).isValid())) break ; }
uint32_t gsfHit1 = gsfHitPattern1.getHitPattern(gsfHitCounter1) ;
uint32_t gsfHit2 = gsfHitPattern2.getHitPattern(gsfHitCounter2) ;
if (gsfHitPattern1.getSubStructure(gsfHit1)!=gsfHitPattern2.getSubStructure(gsfHit2))
{
return (gsfHitPattern2.getSubStructure(gsfHit2)<gsfHitPattern1.getSubStructure(gsfHit1));
}
else if (gsfHitPattern1.getLayer(gsfHit1)!=gsfHitPattern2.getLayer(gsfHit2))
{
return (gsfHitPattern2.getLayer(gsfHit2)<gsfHitPattern1.getLayer(gsfHit1));
}
else
{
sameLayer = true;
return false;
}
}
| bool PFElecTkProducer::isInnerMostWithLostHits | ( | const reco::GsfTrackRef & | nGsfTrack, |
| const reco::GsfTrackRef & | iGsfTrack, | ||
| bool & | sameLayer | ||
| ) | [private] |
Definition at line 1154 of file PFElecTkProducer.cc.
Referenced by resolveGsfTracks().
{
// define closest using the lost hits on the expectedhitsineer
unsigned int nLostHits = nGsfTrack->trackerExpectedHitsInner().numberOfLostHits();
unsigned int iLostHits = iGsfTrack->trackerExpectedHitsInner().numberOfLostHits();
if (nLostHits!=iLostHits) {
return (nLostHits > iLostHits);
}
else {
sameLayer = true;
return false;
}
}
| bool PFElecTkProducer::isSameEgSC | ( | const reco::ElectronSeedRef & | nSeedRef, |
| const reco::ElectronSeedRef & | iSeedRef, | ||
| bool & | bothGsfEcalDriven, | ||
| float & | SCEnergy | ||
| ) | [private] |
Definition at line 892 of file PFElecTkProducer.cc.
References edm::Ref< C, T, F >::isNonnull().
Referenced by resolveGsfTracks().
{
bool isSameSC = false;
if(nSeedRef->caloCluster().isNonnull() && iSeedRef->caloCluster().isNonnull()) {
SuperClusterRef nscRef = nSeedRef->caloCluster().castTo<SuperClusterRef>();
SuperClusterRef iscRef = iSeedRef->caloCluster().castTo<SuperClusterRef>();
if(nscRef.isNonnull() && iscRef.isNonnull()) {
bothGsfEcalDriven = true;
if(nscRef == iscRef) {
isSameSC = true;
// retrieve the supercluster energy
SCEnergy = nscRef->energy();
}
}
}
return isSameSC;
}
| bool PFElecTkProducer::isSharingEcalEnergyWithEgSC | ( | const reco::GsfPFRecTrack & | nGsfPFRecTrack, |
| const reco::GsfPFRecTrack & | iGsfPFRecTrack, | ||
| const reco::ElectronSeedRef & | nSeedRef, | ||
| const reco::ElectronSeedRef & | iSeedRef, | ||
| const reco::PFClusterCollection & | theEClus, | ||
| bool & | bothGsfTrackerDriven, | ||
| bool & | nEcalDriven, | ||
| bool & | iEcalDriven, | ||
| float & | nEnergy, | ||
| float & | iEnergy | ||
| ) | [private] |
Definition at line 915 of file PFElecTkProducer.cc.
References reco::PFCluster::calculatePositionREP(), reco::PFTrajectoryPoint::ECALShowerMax, reco::PFCluster::energy(), reco::PFTrack::extrapolatedPoint(), reco::GsfPFRecTrack::gsfTrackRef(), edm::Ref< C, T, F >::isNonnull(), reco::PFTrajectoryPoint::isValid(), muon::overlap(), reco::GsfPFRecTrack::PFRecBrem(), Pi, reco::CaloCluster::position(), LinkByRecHit::testTrackAndClusterByRecHit(), and TwoPi.
Referenced by resolveGsfTracks().
{
bool isSharingEnergy = false;
//which is EcalDriven?
bool oneEcalDriven = true;
SuperClusterRef scRef;
GsfPFRecTrack gsfPfTrack;
if(nSeedRef->caloCluster().isNonnull()) {
scRef = nSeedRef->caloCluster().castTo<SuperClusterRef>();
nEnergy = scRef->energy();
nEcalDriven = true;
gsfPfTrack = iGsfPFRecTrack;
}
else if(iSeedRef->caloCluster().isNonnull()){
scRef = iSeedRef->caloCluster().castTo<SuperClusterRef>();
iEnergy = scRef->energy();
iEcalDriven = true;
gsfPfTrack = nGsfPFRecTrack;
}
else{
oneEcalDriven = false;
}
if(oneEcalDriven) {
//run a basic reconstruction for the particle flow
vector<PFCluster> vecPFClusters;
vecPFClusters.clear();
for (PFClusterCollection::const_iterator clus = theEClus.begin();
clus != theEClus.end();
clus++ ) {
PFCluster clust = *clus;
clust.calculatePositionREP();
float deta = fabs(scRef->position().eta() - clust.position().eta());
float dphi = fabs(scRef->position().phi() - clust.position().phi());
if (dphi>TMath::Pi()) dphi-= TMath::TwoPi();
// Angle preselection between the supercluster and pfclusters
// this is needed just to save some cpu-time for this is hard-coded
if(deta < 0.5 && fabs(dphi) < 1.0) {
bool foundLink = false;
double distGsf = gsfPfTrack.extrapolatedPoint( reco::PFTrajectoryPoint::ECALShowerMax ).isValid() ?
LinkByRecHit::testTrackAndClusterByRecHit(gsfPfTrack , clust ) : -1.;
// check if it touch the GsfTrack
if(distGsf > 0.) {
if(nEcalDriven)
iEnergy += clust.energy();
else
nEnergy += clust.energy();
vecPFClusters.push_back(clust);
foundLink = true;
}
// check if it touch the Brem-tangents
if(foundLink == false) {
vector<PFBrem> primPFBrem = gsfPfTrack.PFRecBrem();
for(unsigned ipbrem = 0; ipbrem < primPFBrem.size(); ipbrem++) {
if(primPFBrem[ipbrem].indTrajPoint() == 99) continue;
const reco::PFRecTrack& pfBremTrack = primPFBrem[ipbrem];
double dist = pfBremTrack.extrapolatedPoint( reco::PFTrajectoryPoint::ECALShowerMax ).isValid() ?
LinkByRecHit::testTrackAndClusterByRecHit(pfBremTrack , clust, true ) : -1.;
if(dist > 0.) {
if(nEcalDriven)
iEnergy += clust.energy();
else
nEnergy += clust.energy();
vecPFClusters.push_back(clust);
foundLink = true;
}
}
}
} // END if anble preselection
} // PFClusters Loop
if(vecPFClusters.size() > 0 ) {
for(unsigned int pf = 0; pf < vecPFClusters.size(); pf++) {
bool isCommon = ClusterClusterMapping::overlap(vecPFClusters[pf],*scRef);
if(isCommon) {
isSharingEnergy = true;
}
break;
}
}
}
else {
// both tracks are trackerDriven, try ECAL energy matching also in this case.
bothGsfTrackerDriven = true;
vector<PFCluster> nPFCluster;
vector<PFCluster> iPFCluster;
nPFCluster.clear();
iPFCluster.clear();
for (PFClusterCollection::const_iterator clus = theEClus.begin();
clus != theEClus.end();
clus++ ) {
PFCluster clust = *clus;
clust.calculatePositionREP();
float ndeta = fabs(nGsfPFRecTrack.gsfTrackRef()->eta() - clust.position().eta());
float ndphi = fabs(nGsfPFRecTrack.gsfTrackRef()->phi() - clust.position().phi());
if (ndphi>TMath::Pi()) ndphi-= TMath::TwoPi();
// Apply loose preselection with the track
// just to save cpu time, for this hard-coded
if(ndeta < 0.5 && fabs(ndphi) < 1.0) {
bool foundNLink = false;
double distGsf = nGsfPFRecTrack.extrapolatedPoint( reco::PFTrajectoryPoint::ECALShowerMax ).isValid() ?
LinkByRecHit::testTrackAndClusterByRecHit(nGsfPFRecTrack , clust ) : -1.;
if(distGsf > 0.) {
nPFCluster.push_back(clust);
nEnergy += clust.energy();
foundNLink = true;
}
if(foundNLink == false) {
const vector<PFBrem>& primPFBrem = nGsfPFRecTrack.PFRecBrem();
for(unsigned ipbrem = 0; ipbrem < primPFBrem.size(); ipbrem++) {
if(primPFBrem[ipbrem].indTrajPoint() == 99) continue;
const reco::PFRecTrack& pfBremTrack = primPFBrem[ipbrem];
if(foundNLink == false) {
double dist = pfBremTrack.extrapolatedPoint( reco::PFTrajectoryPoint::ECALShowerMax ).isValid() ?
LinkByRecHit::testTrackAndClusterByRecHit(pfBremTrack , clust, true ) : -1.;
if(dist > 0.) {
nPFCluster.push_back(clust);
nEnergy += clust.energy();
foundNLink = true;
}
}
}
}
}
float ideta = fabs(iGsfPFRecTrack.gsfTrackRef()->eta() - clust.position().eta());
float idphi = fabs(iGsfPFRecTrack.gsfTrackRef()->phi() - clust.position().phi());
if (idphi>TMath::Pi()) idphi-= TMath::TwoPi();
// Apply loose preselection with the track
// just to save cpu time, for this hard-coded
if(ideta < 0.5 && fabs(idphi) < 1.0) {
bool foundILink = false;
double dist = iGsfPFRecTrack.extrapolatedPoint( reco::PFTrajectoryPoint::ECALShowerMax ).isValid() ?
LinkByRecHit::testTrackAndClusterByRecHit(iGsfPFRecTrack , clust ) : -1.;
if(dist > 0.) {
iPFCluster.push_back(clust);
iEnergy += clust.energy();
foundILink = true;
}
if(foundILink == false) {
vector<PFBrem> primPFBrem = iGsfPFRecTrack.PFRecBrem();
for(unsigned ipbrem = 0; ipbrem < primPFBrem.size(); ipbrem++) {
if(primPFBrem[ipbrem].indTrajPoint() == 99) continue;
const reco::PFRecTrack& pfBremTrack = primPFBrem[ipbrem];
if(foundILink == false) {
double dist = LinkByRecHit::testTrackAndClusterByRecHit(pfBremTrack , clust, true );
if(dist > 0.) {
iPFCluster.push_back(clust);
iEnergy += clust.energy();
foundILink = true;
}
}
}
}
}
}
if(nPFCluster.size() > 0 && iPFCluster.size() > 0) {
for(unsigned int npf = 0; npf < nPFCluster.size(); npf++) {
for(unsigned int ipf = 0; ipf < iPFCluster.size(); ipf++) {
bool isCommon = ClusterClusterMapping::overlap(nPFCluster[npf],iPFCluster[ipf]);
if(isCommon) {
isSharingEnergy = true;
break;
}
}
if(isSharingEnergy)
break;
}
}
}
return isSharingEnergy;
}
| float PFElecTkProducer::minTangDist | ( | const reco::GsfPFRecTrack & | primGsf, |
| const reco::GsfPFRecTrack & | secGsf | ||
| ) | [private] |
Definition at line 848 of file PFElecTkProducer.cc.
References reco::PFTrajectoryPoint::ECALEntrance, reco::PFTrajectoryPoint::isValid(), reco::GsfPFRecTrack::PFRecBrem(), Pi, reco::PFTrajectoryPoint::positionREP(), and TwoPi.
Referenced by resolveGsfTracks().
{
float minDphi = 1000.;
std::vector<reco::PFBrem> primPFBrem = primGsf.PFRecBrem();
std::vector<reco::PFBrem> secPFBrem = secGsf.PFRecBrem();
unsigned int cbrem = 0;
for (unsigned isbrem = 0; isbrem < secPFBrem.size(); isbrem++) {
if(secPFBrem[isbrem].indTrajPoint() == 99) continue;
const reco::PFTrajectoryPoint& atSecECAL
= secPFBrem[isbrem].extrapolatedPoint( reco::PFTrajectoryPoint::ECALEntrance );
if( ! atSecECAL.isValid() ) continue;
float secPhi = atSecECAL.positionREP().Phi();
unsigned int sbrem = 0;
for(unsigned ipbrem = 0; ipbrem < primPFBrem.size(); ipbrem++) {
if(primPFBrem[ipbrem].indTrajPoint() == 99) continue;
const reco::PFTrajectoryPoint& atPrimECAL
= primPFBrem[ipbrem].extrapolatedPoint( reco::PFTrajectoryPoint::ECALEntrance );
if( ! atPrimECAL.isValid() ) continue;
sbrem++;
if(sbrem <= 3) {
float primPhi = atPrimECAL.positionREP().Phi();
float dphi = fabs(primPhi - secPhi);
if (dphi>TMath::Pi()) dphi-= TMath::TwoPi();
if(fabs(dphi) < minDphi) {
minDphi = fabs(dphi);
}
}
}
cbrem++;
if(cbrem == 3)
break;
}
return minDphi;
}
| void PFElecTkProducer::produce | ( | edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [private, virtual] |
Produce the PFRecTrack collection.
Implements edm::EDProducer.
Definition at line 130 of file PFElecTkProducer.cc.
References PFTrackTransformer::addPointsAndBrems(), applyGsfClean_, applySel_, applySelection(), DeDxDiscriminatorTools::charge(), convBremFinder_, createGsfPFRecTrackRef(), edm::EventID::event(), FindPfRef(), newFWLiteAna::found, ConvBremPFTrackFinder::foundConvBremPFRecTrack(), edm::Event::getByLabel(), ConvBremPFTrackFinder::getConvBremPFRecTracks(), reco::PFRecTrack::GSF, gsfTrackLabel_, edm::EventBase::id(), isFifthStep(), edm::Ref< C, T, F >::isNonnull(), edm::Ref< C, T, F >::isNull(), LogDebug, modemomentum_, mtsTransform_, pfConv_, pfConversions_cfi::pfConversions, pfEcalClusters_, pfNuclear_, pftrack_, pfTrackLabel_, pfTransformer_, pfV0_cfi::pfV0, pfV0_, primVtxLabel_, edm::Handle< T >::product(), edm::Event::put(), resolveGsfTracks(), edm::EventID::run(), secpftrack_, trajinev_, useConversions_, useFifthStepForEcalDriven_, useFifthStepForTrackDriven_, useNuclear_, and useV0_.
{
LogDebug("PFElecTkProducer")<<"START event: "<<iEvent.id().event()
<<" in run "<<iEvent.id().run();
//create the empty collections
auto_ptr< GsfPFRecTrackCollection >
gsfPFRecTrackCollection(new GsfPFRecTrackCollection);
auto_ptr< GsfPFRecTrackCollection >
gsfPFRecTrackCollectionSecondary(new GsfPFRecTrackCollection);
//read collections of tracks
Handle<GsfTrackCollection> gsftrackscoll;
iEvent.getByLabel(gsfTrackLabel_,gsftrackscoll);
//read collections of trajectories
Handle<vector<Trajectory> > TrajectoryCollection;
//read pfrectrack collection
Handle<PFRecTrackCollection> thePfRecTrackCollection;
iEvent.getByLabel(pfTrackLabel_,thePfRecTrackCollection);
const PFRecTrackCollection& PfRTkColl = *(thePfRecTrackCollection.product());
// PFClusters
Handle<PFClusterCollection> theECPfClustCollection;
iEvent.getByLabel(pfEcalClusters_,theECPfClustCollection);
const PFClusterCollection& theEcalClusters = *(theECPfClustCollection.product());
//Primary Vertexes
Handle<reco::VertexCollection> thePrimaryVertexColl;
iEvent.getByLabel(primVtxLabel_,thePrimaryVertexColl);
// Displaced Vertex
Handle< reco::PFDisplacedTrackerVertexCollection > pfNuclears;
if( useNuclear_ ) {
bool found = iEvent.getByLabel(pfNuclear_, pfNuclears);
if(!found )
LogError("PFElecTkProducer")<<" cannot get PFNuclear : "
<< pfNuclear_
<< " please set useNuclear=False in RecoParticleFlow/PFTracking/python/pfTrackElec_cfi.py" << endl;
}
// Conversions
Handle< reco::PFConversionCollection > pfConversions;
if( useConversions_ ) {
bool found = iEvent.getByLabel(pfConv_,pfConversions);
if(!found )
LogError("PFElecTkProducer")<<" cannot get PFConversions : "
<< pfConv_
<< " please set useConversions=False in RecoParticleFlow/PFTracking/python/pfTrackElec_cfi.py" << endl;
}
// V0
Handle< reco::PFV0Collection > pfV0;
if( useV0_ ) {
bool found = iEvent.getByLabel(pfV0_, pfV0);
if(!found )
LogError("PFElecTkProducer")<<" cannot get PFV0 : "
<< pfV0_
<< " please set useV0=False RecoParticleFlow/PFTracking/python/pfTrackElec_cfi.py" << endl;
}
GsfTrackCollection gsftracks = *(gsftrackscoll.product());
vector<Trajectory> tjvec(0);
if (trajinev_){
bool foundTraj = iEvent.getByLabel(gsfTrackLabel_,TrajectoryCollection);
if(!foundTraj)
LogError("PFElecTkProducer")
<<" cannot get Trajectories of : "
<< gsfTrackLabel_
<< " please set TrajInEvents = False in RecoParticleFlow/PFTracking/python/pfTrackElec_cfi.py" << endl;
tjvec= *(TrajectoryCollection.product());
}
vector<reco::GsfPFRecTrack> selGsfPFRecTracks;
vector<reco::GsfPFRecTrack> primaryGsfPFRecTracks;
std::map<unsigned int, std::vector<reco::GsfPFRecTrack> > GsfPFMap;
for (unsigned int igsf=0; igsf<gsftracks.size();igsf++) {
GsfTrackRef trackRef(gsftrackscoll, igsf);
int kf_ind=FindPfRef(PfRTkColl,gsftracks[igsf],false);
if (kf_ind>=0) {
PFRecTrackRef kf_ref(thePfRecTrackCollection,
kf_ind);
// remove fifth step tracks
if( useFifthStepForEcalDriven_ == false
|| useFifthStepForTrackDriven_ == false) {
bool isFifthStepTrack = isFifthStep(kf_ref);
bool isEcalDriven = true;
bool isTrackerDriven = true;
if (&(*trackRef->seedRef())==0) {
isEcalDriven = false;
isTrackerDriven = false;
}
else {
ElectronSeedRef SeedRef= trackRef->extra()->seedRef().castTo<ElectronSeedRef>();
if(SeedRef->caloCluster().isNull())
isEcalDriven = false;
if(SeedRef->ctfTrack().isNull())
isTrackerDriven = false;
}
//note: the same track could be both ecalDriven and trackerDriven
if(isFifthStepTrack &&
isEcalDriven &&
isTrackerDriven == false &&
useFifthStepForEcalDriven_ == false) {
continue;
}
if(isFifthStepTrack &&
isTrackerDriven &&
isEcalDriven == false &&
useFifthStepForTrackDriven_ == false) {
continue;
}
if(isFifthStepTrack &&
isTrackerDriven &&
isEcalDriven &&
useFifthStepForTrackDriven_ == false &&
useFifthStepForEcalDriven_ == false) {
continue;
}
}
pftrack_=GsfPFRecTrack( gsftracks[igsf].charge(),
reco::PFRecTrack::GSF,
igsf, trackRef,
kf_ref);
} else {
PFRecTrackRef dummyRef;
pftrack_=GsfPFRecTrack( gsftracks[igsf].charge(),
reco::PFRecTrack::GSF,
igsf, trackRef,
dummyRef);
}
bool validgsfbrem = false;
if(trajinev_) {
validgsfbrem = pfTransformer_->addPointsAndBrems(pftrack_,
gsftracks[igsf],
tjvec[igsf],
modemomentum_);
} else {
validgsfbrem = pfTransformer_->addPointsAndBrems(pftrack_,
gsftracks[igsf],
mtsTransform_);
}
bool passSel = true;
if(applySel_)
passSel = applySelection(gsftracks[igsf]);
if(validgsfbrem && passSel)
selGsfPFRecTracks.push_back(pftrack_);
}
unsigned int count_primary = 0;
if(selGsfPFRecTracks.size() > 0) {
for(unsigned int ipfgsf=0; ipfgsf<selGsfPFRecTracks.size();ipfgsf++) {
vector<unsigned int> secondaries(0);
secondaries.clear();
bool keepGsf = true;
if(applyGsfClean_) {
keepGsf = resolveGsfTracks(selGsfPFRecTracks,ipfgsf,secondaries,theEcalClusters);
}
//is primary?
if(keepGsf == true) {
// Find kf tracks from converted brem photons
if(convBremFinder_->foundConvBremPFRecTrack(thePfRecTrackCollection,thePrimaryVertexColl,
pfNuclears,pfConversions,pfV0,
useNuclear_,useConversions_,useV0_,
theEcalClusters,selGsfPFRecTracks[ipfgsf])) {
const vector<PFRecTrackRef>& convBremPFRecTracks(convBremFinder_->getConvBremPFRecTracks());
for(unsigned int ii = 0; ii<convBremPFRecTracks.size(); ii++) {
selGsfPFRecTracks[ipfgsf].addConvBremPFRecTrackRef(convBremPFRecTracks[ii]);
}
}
// save primaries gsf tracks
// gsfPFRecTrackCollection->push_back(selGsfPFRecTracks[ipfgsf]);
primaryGsfPFRecTracks.push_back(selGsfPFRecTracks[ipfgsf]);
// NOTE:: THE TRACKID IS USED TO LINK THE PRIMARY GSF TRACK. THIS NEEDS
// TO BE CHANGED AS SOON AS IT IS POSSIBLE TO CHANGE DATAFORMATS
// A MODIFICATION HERE IMPLIES A MODIFICATION IN PFBLOCKALGO.CC/H
unsigned int primGsfIndex = selGsfPFRecTracks[ipfgsf].trackId();
vector<reco::GsfPFRecTrack> trueGsfPFRecTracks;
if(secondaries.size() > 0) {
// loop on secondaries gsf tracks (from converted brems)
for(unsigned int isecpfgsf=0; isecpfgsf<secondaries.size();isecpfgsf++) {
PFRecTrackRef refsecKF = selGsfPFRecTracks[(secondaries[isecpfgsf])].kfPFRecTrackRef();
unsigned int secGsfIndex = selGsfPFRecTracks[(secondaries[isecpfgsf])].trackId();
GsfTrackRef secGsfRef = selGsfPFRecTracks[(secondaries[isecpfgsf])].gsfTrackRef();
if(refsecKF.isNonnull()) {
// NOTE::IT SAVED THE TRACKID OF THE PRIMARY!!! THIS IS USED IN PFBLOCKALGO.CC/H
secpftrack_= GsfPFRecTrack( gsftracks[secGsfIndex].charge(),
reco::PFRecTrack::GSF,
primGsfIndex, secGsfRef,
refsecKF);
}
else{
PFRecTrackRef dummyRef;
// NOTE::IT SAVED THE TRACKID OF THE PRIMARY!!! THIS IS USED IN PFBLOCKALGO.CC/H
secpftrack_= GsfPFRecTrack( gsftracks[secGsfIndex].charge(),
reco::PFRecTrack::GSF,
primGsfIndex, secGsfRef,
dummyRef);
}
bool validgsfbrem = false;
if(trajinev_) {
validgsfbrem = pfTransformer_->addPointsAndBrems(secpftrack_,
gsftracks[secGsfIndex],
tjvec[secGsfIndex],
modemomentum_);
} else {
validgsfbrem = pfTransformer_->addPointsAndBrems(secpftrack_,
gsftracks[secGsfIndex],
mtsTransform_);
}
if(validgsfbrem) {
gsfPFRecTrackCollectionSecondary->push_back(secpftrack_);
trueGsfPFRecTracks.push_back(secpftrack_);
}
}
}
GsfPFMap.insert(pair<unsigned int,std::vector<reco::GsfPFRecTrack> >(count_primary,trueGsfPFRecTracks));
trueGsfPFRecTracks.clear();
count_primary++;
}
}
}
const edm::OrphanHandle<GsfPFRecTrackCollection> gsfPfRefProd =
iEvent.put(gsfPFRecTrackCollectionSecondary,"Secondary");
//now the secondary GsfPFRecTracks are in the event, the Ref can be created
createGsfPFRecTrackRef(gsfPfRefProd,primaryGsfPFRecTracks,GsfPFMap);
for(unsigned int iGSF = 0; iGSF<primaryGsfPFRecTracks.size();iGSF++){
gsfPFRecTrackCollection->push_back(primaryGsfPFRecTracks[iGSF]);
}
iEvent.put(gsfPFRecTrackCollection);
selGsfPFRecTracks.clear();
GsfPFMap.clear();
primaryGsfPFRecTracks.clear();
}
| bool PFElecTkProducer::resolveGsfTracks | ( | const std::vector< reco::GsfPFRecTrack > & | GsfPFVec, |
| unsigned int | ngsf, | ||
| std::vector< unsigned int > & | secondaries, | ||
| const reco::PFClusterCollection & | theEClus | ||
| ) | [private] |
Definition at line 578 of file PFElecTkProducer.cc.
References applyAngularGsfClean_, gather_cfg::cout, debugGsfCleaning_, detaCutGsfClean_, dphiCutGsfClean_, MultiTrajectoryStateTransform::innerStateOnSurface(), isInnerMostWithLostHits(), isSameEgSC(), isSharingEcalEnergyWithEgSC(), TrajectoryStateOnSurface::isValid(), PV3DBase< T, PVType, FrameType >::mag(), minTangDist(), MultiTrajectoryStateMode::momentumFromModeCartesian(), mtsMode_, mtsTransform_, Pi, and TwoPi.
Referenced by produce().
{
bool debugCleaning = debugGsfCleaning_;
bool n_keepGsf = true;
reco::GsfTrackRef nGsfTrack = GsfPFVec[ngsf].gsfTrackRef();
if (&(*nGsfTrack->seedRef())==0) return false;
ElectronSeedRef nElSeedRef=nGsfTrack->extra()->seedRef().castTo<ElectronSeedRef>();
TrajectoryStateOnSurface inTSOS = mtsTransform_.innerStateOnSurface((*nGsfTrack));
GlobalVector ninnMom;
float nPin = nGsfTrack->pMode();
if(inTSOS.isValid()){
mtsMode_->momentumFromModeCartesian(inTSOS,ninnMom);
nPin = ninnMom.mag();
}
float neta = nGsfTrack->innerMomentum().eta();
float nphi = nGsfTrack->innerMomentum().phi();
if(debugCleaning)
cout << " PFElecTkProducer:: considering track " << nGsfTrack->pt()
<< " eta,phi " << nGsfTrack->eta() << ", " << nGsfTrack->phi() << endl;
for (unsigned int igsf=0; igsf< GsfPFVec.size();igsf++) {
if(igsf != ngsf ) {
reco::GsfTrackRef iGsfTrack = GsfPFVec[igsf].gsfTrackRef();
if(debugCleaning)
cout << " PFElecTkProducer:: and comparing with track " << iGsfTrack->pt()
<< " eta,phi " << iGsfTrack->eta() << ", " << iGsfTrack->phi() << endl;
float ieta = iGsfTrack->innerMomentum().eta();
float iphi = iGsfTrack->innerMomentum().phi();
float feta = fabs(neta - ieta);
float fphi = fabs(nphi - iphi);
if (fphi>TMath::Pi()) fphi-= TMath::TwoPi();
// apply a superloose preselection only to avoid un-useful cpu time: hard-coded for this reason
if(feta < 0.5 && fabs(fphi) < 1.0) {
if(debugCleaning)
cout << " Entering angular superloose preselection " << endl;
TrajectoryStateOnSurface i_inTSOS = mtsTransform_.innerStateOnSurface((*iGsfTrack));
GlobalVector i_innMom;
float iPin = iGsfTrack->pMode();
if(i_inTSOS.isValid()){
mtsMode_->momentumFromModeCartesian(i_inTSOS,i_innMom);
iPin = i_innMom.mag();
}
if (&(*iGsfTrack->seedRef())==0) continue;
ElectronSeedRef iElSeedRef=iGsfTrack->extra()->seedRef().castTo<ElectronSeedRef>();
float SCEnergy = -1.;
// Check if two tracks match the same SC
bool areBothGsfEcalDriven = false;;
bool isSameSC = isSameEgSC(nElSeedRef,iElSeedRef,areBothGsfEcalDriven,SCEnergy);
// CASE1 both GsfTracks ecalDriven and match the same SC
if(areBothGsfEcalDriven ) {
if(isSameSC) {
float nEP = SCEnergy/nPin;
float iEP = SCEnergy/iPin;
if(debugCleaning)
cout << " Entering SAME supercluster case "
<< " nEP " << nEP
<< " iEP " << iEP << endl;
// if same SC take the closest or if same
// distance the best E/p
// Innermost using LostHits technology
bool isSameLayer = false;
bool iGsfInnermostWithLostHits =
isInnerMostWithLostHits(nGsfTrack,iGsfTrack,isSameLayer);
if(debugCleaning)
cout << " iGsf is InnerMostWithLostHits " << iGsfInnermostWithLostHits
<< " isSameLayer " << isSameLayer << endl;
if (iGsfInnermostWithLostHits) {
n_keepGsf = false;
return n_keepGsf;
}
else if(isSameLayer){
if(fabs(iEP-1) < fabs(nEP-1)) {
n_keepGsf = false;
return n_keepGsf;
}
else{
secondaries.push_back(igsf);
}
}
else {
// save secondaries gsf track (put selection)
secondaries.push_back(igsf);
}
} // end same SC case
}
else {
// enter in the condition where at least one track is trackerDriven
float minBremDphi = minTangDist(GsfPFVec[ngsf],GsfPFVec[igsf]);
float nETot = 0.;
float iETot = 0.;
bool isBothGsfTrackerDriven = false;
bool nEcalDriven = false;
bool iEcalDriven = false;
bool isSameScEgPf = isSharingEcalEnergyWithEgSC(GsfPFVec[ngsf],
GsfPFVec[igsf],
nElSeedRef,
iElSeedRef,
theEClus,
isBothGsfTrackerDriven,
nEcalDriven,
iEcalDriven,
nETot,
iETot);
// check if the first hit of iGsfTrack < nGsfTrack
bool isSameLayer = false;
bool iGsfInnermostWithLostHits =
isInnerMostWithLostHits(nGsfTrack,iGsfTrack,isSameLayer);
if(isSameScEgPf) {
// CASE 2 : One Gsf has reference to a SC and the other one not or both not
if(debugCleaning) {
cout << " Sharing ECAL energy passed "
<< " nEtot " << nETot
<< " iEtot " << iETot << endl;
if(isBothGsfTrackerDriven)
cout << " Both Track are trackerDriven " << endl;
}
// Innermost using LostHits technology
if (iGsfInnermostWithLostHits) {
n_keepGsf = false;
return n_keepGsf;
}
else if(isSameLayer){
// Thirt Case: One Gsf has reference to a SC and the other one not or both not
// gsf tracks starts from the same layer
// check number of sharing modules (at least 50%)
// check number of sharing hits (at least 2)
// check charge flip inner/outer
// they share energy
if(isBothGsfTrackerDriven == false) {
// if at least one Gsf track is EcalDriven choose that one.
if(iEcalDriven) {
n_keepGsf = false;
return n_keepGsf;
}
else {
secondaries.push_back(igsf);
}
}
else {
// if both tracks are tracker driven choose that one with the best E/p
// with ETot = max(En,Ei)
float ETot = -1;
if(nETot != iETot) {
if(nETot > iETot)
ETot = nETot;
else
ETot = iETot;
}
else {
ETot = nETot;
}
float nEP = ETot/nPin;
float iEP = ETot/iPin;
if(debugCleaning)
cout << " nETot " << nETot
<< " iETot " << iETot
<< " ETot " << ETot << endl
<< " nPin " << nPin
<< " iPin " << iPin
<< " nEP " << nEP
<< " iEP " << iEP << endl;
if(fabs(iEP-1) < fabs(nEP-1)) {
n_keepGsf = false;
return n_keepGsf;
}
else{
secondaries.push_back(igsf);
}
}
}
else {
secondaries.push_back(igsf);
}
}
else if(feta < detaCutGsfClean_ && minBremDphi < dphiCutGsfClean_) {
// very close tracks
bool secPushedBack = false;
if(nEcalDriven == false && nETot == 0.) {
n_keepGsf = false;
return n_keepGsf;
}
else if(iEcalDriven == false && iETot == 0.) {
secondaries.push_back(igsf);
secPushedBack = true;
}
if(debugCleaning)
cout << " Close Tracks "
<< " feta " << feta << " fabs(fphi) " << fabs(fphi)
<< " minBremDphi " << minBremDphi
<< " nETot " << nETot
<< " iETot " << iETot
<< " nLostHits " << nGsfTrack->trackerExpectedHitsInner().numberOfLostHits()
<< " iLostHits " << iGsfTrack->trackerExpectedHitsInner().numberOfLostHits() << endl;
// apply selection only if one track has lost hits
if(applyAngularGsfClean_) {
if (iGsfInnermostWithLostHits) {
n_keepGsf = false;
return n_keepGsf;
}
else if(isSameLayer == false) {
if(secPushedBack == false)
secondaries.push_back(igsf);
}
}
}
else if(feta < 0.1 && minBremDphi < 0.2){
// failed all the conditions, discard only tracker driven tracks
// with no PFClusters linked.
if(debugCleaning)
cout << " Close Tracks and failed all the conditions "
<< " feta " << feta << " fabs(fphi) " << fabs(fphi)
<< " minBremDphi " << minBremDphi
<< " nETot " << nETot
<< " iETot " << iETot
<< " nLostHits " << nGsfTrack->trackerExpectedHitsInner().numberOfLostHits()
<< " iLostHits " << iGsfTrack->trackerExpectedHitsInner().numberOfLostHits() << endl;
if(nEcalDriven == false && nETot == 0.) {
n_keepGsf = false;
return n_keepGsf;
}
// Here I do not push back the secondary because considered fakes...
}
}
}
}
}
return n_keepGsf;
}
bool PFElecTkProducer::applyAngularGsfClean_ [private] |
Definition at line 109 of file PFElecTkProducer.h.
Referenced by PFElecTkProducer(), and resolveGsfTracks().
bool PFElecTkProducer::applyGsfClean_ [private] |
Definition at line 124 of file PFElecTkProducer.h.
Referenced by PFElecTkProducer(), and produce().
bool PFElecTkProducer::applySel_ [private] |
Definition at line 123 of file PFElecTkProducer.h.
Referenced by PFElecTkProducer(), and produce().
edm::ParameterSet PFElecTkProducer::conf_ [private] |
Definition at line 98 of file PFElecTkProducer.h.
Definition at line 117 of file PFElecTkProducer.h.
Referenced by beginRun(), produce(), and ~PFElecTkProducer().
bool PFElecTkProducer::debugGsfCleaning_ [private] |
Definition at line 129 of file PFElecTkProducer.h.
Referenced by PFElecTkProducer(), and resolveGsfTracks().
double PFElecTkProducer::detaCutGsfClean_ [private] |
Definition at line 110 of file PFElecTkProducer.h.
Referenced by PFElecTkProducer(), and resolveGsfTracks().
double PFElecTkProducer::detaGsfSC_ [private] |
Definition at line 131 of file PFElecTkProducer.h.
Referenced by applySelection(), and PFElecTkProducer().
double PFElecTkProducer::dphiCutGsfClean_ [private] |
Definition at line 111 of file PFElecTkProducer.h.
Referenced by PFElecTkProducer(), and resolveGsfTracks().
double PFElecTkProducer::dphiGsfSC_ [private] |
Definition at line 132 of file PFElecTkProducer.h.
Referenced by applySelection(), and PFElecTkProducer().
Definition at line 99 of file PFElecTkProducer.h.
Referenced by PFElecTkProducer(), and produce().
double PFElecTkProducer::maxPtConvReco_ [private] |
Definition at line 133 of file PFElecTkProducer.h.
Referenced by PFElecTkProducer().
bool PFElecTkProducer::modemomentum_ [private] |
Definition at line 122 of file PFElecTkProducer.h.
Referenced by PFElecTkProducer(), and produce().
const MultiTrajectoryStateMode* PFElecTkProducer::mtsMode_ [private] |
Definition at line 115 of file PFElecTkProducer.h.
Referenced by resolveGsfTracks().
Definition at line 116 of file PFElecTkProducer.h.
Referenced by beginRun(), produce(), and resolveGsfTracks().
double PFElecTkProducer::mvaConvBremFinderID_ [private] |
Definition at line 137 of file PFElecTkProducer.h.
Referenced by beginRun(), and PFElecTkProducer().
std::string PFElecTkProducer::path_mvaWeightFileConvBrem_ [private] |
Definition at line 138 of file PFElecTkProducer.h.
Referenced by beginRun(), and PFElecTkProducer().
edm::InputTag PFElecTkProducer::pfConv_ [private] |
Definition at line 104 of file PFElecTkProducer.h.
Referenced by PFElecTkProducer(), and produce().
Definition at line 102 of file PFElecTkProducer.h.
Referenced by PFElecTkProducer(), and produce().
edm::InputTag PFElecTkProducer::pfNuclear_ [private] |
Definition at line 103 of file PFElecTkProducer.h.
Referenced by PFElecTkProducer(), and produce().
Definition at line 96 of file PFElecTkProducer.h.
Referenced by produce().
edm::InputTag PFElecTkProducer::pfTrackLabel_ [private] |
Definition at line 100 of file PFElecTkProducer.h.
Referenced by PFElecTkProducer(), and produce().
Definition at line 114 of file PFElecTkProducer.h.
Referenced by beginRun(), endRun(), produce(), and ~PFElecTkProducer().
edm::InputTag PFElecTkProducer::pfV0_ [private] |
Definition at line 105 of file PFElecTkProducer.h.
Referenced by PFElecTkProducer(), and produce().
edm::InputTag PFElecTkProducer::primVtxLabel_ [private] |
Definition at line 101 of file PFElecTkProducer.h.
Referenced by PFElecTkProducer(), and produce().
double PFElecTkProducer::SCEne_ [private] |
Definition at line 130 of file PFElecTkProducer.h.
Referenced by applySelection(), and PFElecTkProducer().
Definition at line 97 of file PFElecTkProducer.h.
Referenced by produce().
bool PFElecTkProducer::trajinev_ [private] |
Trajectory of GSfTracks in the event?
Definition at line 121 of file PFElecTkProducer.h.
Referenced by PFElecTkProducer(), and produce().
bool PFElecTkProducer::useConvBremFinder_ [private] |
Conv Brem Finder.
Definition at line 136 of file PFElecTkProducer.h.
Referenced by beginRun(), and PFElecTkProducer().
bool PFElecTkProducer::useConversions_ [private] |
Definition at line 107 of file PFElecTkProducer.h.
Referenced by PFElecTkProducer(), and produce().
bool PFElecTkProducer::useFifthStep_ [private] |
Definition at line 125 of file PFElecTkProducer.h.
bool PFElecTkProducer::useFifthStepForEcalDriven_ [private] |
Definition at line 126 of file PFElecTkProducer.h.
Referenced by PFElecTkProducer(), and produce().
bool PFElecTkProducer::useFifthStepForTrackDriven_ [private] |
Definition at line 127 of file PFElecTkProducer.h.
Referenced by PFElecTkProducer(), and produce().
bool PFElecTkProducer::useNuclear_ [private] |
Definition at line 106 of file PFElecTkProducer.h.
Referenced by PFElecTkProducer(), and produce().
bool PFElecTkProducer::useV0_ [private] |
Definition at line 108 of file PFElecTkProducer.h.
Referenced by PFElecTkProducer(), and produce().
1.7.3