|
|
|
Producer for particle flow reconstructed particles (PFCandidates) More...
#include <PFProducer.h>
Public Member Functions | |
| virtual void | beginJob () |
| virtual void | beginRun (edm::Run &, const edm::EventSetup &) |
| PFProducer (const edm::ParameterSet &) | |
| virtual void | produce (edm::Event &, const edm::EventSetup &) |
| ~PFProducer () | |
Private Attributes | |
| std::string | electronExtraOutputCol_ |
| std::string | electronOutputCol_ |
| edm::InputTag | inputTagBlocks_ |
| std::vector< edm::InputTag > | inputTagCleanedHF_ |
| edm::InputTag | inputTagEgammaElectrons_ |
| edm::InputTag | inputTagMuons_ |
| std::auto_ptr< PFAlgo > | pfAlgo_ |
| particle flow algorithm | |
| bool | postMuonCleaning_ |
| bool | useCalibrationsFromDB_ |
| bool | useEGammaElectrons_ |
| bool | usePFElectrons_ |
| bool | usePFPhotons_ |
| bool | useVerticesForNeutral_ |
| bool | verbose_ |
| verbose ? | |
| edm::InputTag | vertices_ |
Producer for particle flow reconstructed particles (PFCandidates)
This producer makes use of PFAlgo, the particle flow algorithm.
Definition at line 29 of file PFProducer.h.
| PFProducer::PFProducer | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 27 of file PFProducer.cc.
References edm::ParameterSet::getParameter(), and edm::ParameterSet::getUntrackedParameter().
{
//--ab: get calibration factors for HF:
bool calibHF_use;
std::vector<double> calibHF_eta_step;
std::vector<double> calibHF_a_EMonly;
std::vector<double> calibHF_b_HADonly;
std::vector<double> calibHF_a_EMHAD;
std::vector<double> calibHF_b_EMHAD;
calibHF_use = iConfig.getParameter<bool>("calibHF_use");
calibHF_eta_step = iConfig.getParameter<std::vector<double> >("calibHF_eta_step");
calibHF_a_EMonly = iConfig.getParameter<std::vector<double> >("calibHF_a_EMonly");
calibHF_b_HADonly = iConfig.getParameter<std::vector<double> >("calibHF_b_HADonly");
calibHF_a_EMHAD = iConfig.getParameter<std::vector<double> >("calibHF_a_EMHAD");
calibHF_b_EMHAD = iConfig.getParameter<std::vector<double> >("calibHF_b_EMHAD");
boost::shared_ptr<PFEnergyCalibrationHF>
thepfEnergyCalibrationHF ( new PFEnergyCalibrationHF(calibHF_use,calibHF_eta_step,calibHF_a_EMonly,calibHF_b_HADonly,calibHF_a_EMHAD,calibHF_b_EMHAD) ) ;
//-----------------
inputTagBlocks_
= iConfig.getParameter<InputTag>("blocks");
//Post cleaning of the muons
inputTagMuons_
= iConfig.getParameter<InputTag>("muons");
postMuonCleaning_
= iConfig.getParameter<bool>("postMuonCleaning");
usePFElectrons_
= iConfig.getParameter<bool>("usePFElectrons");
usePFPhotons_
= iConfig.getParameter<bool>("usePFPhotons");
useEGammaElectrons_
= iConfig.getParameter<bool>("useEGammaElectrons");
if( useEGammaElectrons_) {
inputTagEgammaElectrons_ = iConfig.getParameter<edm::InputTag>("egammaElectrons");
}
electronOutputCol_
= iConfig.getParameter<std::string>("pf_electron_output_col");
bool usePFSCEleCalib;
std::vector<double> calibPFSCEle_Fbrem_barrel;
std::vector<double> calibPFSCEle_Fbrem_endcap;
std::vector<double> calibPFSCEle_barrel;
std::vector<double> calibPFSCEle_endcap;
usePFSCEleCalib = iConfig.getParameter<bool>("usePFSCEleCalib");
calibPFSCEle_Fbrem_barrel = iConfig.getParameter<std::vector<double> >("calibPFSCEle_Fbrem_barrel");
calibPFSCEle_Fbrem_endcap = iConfig.getParameter<std::vector<double> >("calibPFSCEle_Fbrem_endcap");
calibPFSCEle_barrel = iConfig.getParameter<std::vector<double> >("calibPFSCEle_barrel");
calibPFSCEle_endcap = iConfig.getParameter<std::vector<double> >("calibPFSCEle_endcap");
boost::shared_ptr<PFSCEnergyCalibration>
thePFSCEnergyCalibration ( new PFSCEnergyCalibration(calibPFSCEle_Fbrem_barrel,calibPFSCEle_Fbrem_endcap,
calibPFSCEle_barrel,calibPFSCEle_endcap ));
bool useEGammaSupercluster = iConfig.getParameter<bool>("useEGammaSupercluster");
double sumEtEcalIsoForEgammaSC_barrel = iConfig.getParameter<double>("sumEtEcalIsoForEgammaSC_barrel");
double sumEtEcalIsoForEgammaSC_endcap = iConfig.getParameter<double>("sumEtEcalIsoForEgammaSC_endcap");
double coneEcalIsoForEgammaSC = iConfig.getParameter<double>("coneEcalIsoForEgammaSC");
double sumPtTrackIsoForEgammaSC_barrel = iConfig.getParameter<double>("sumPtTrackIsoForEgammaSC_barrel");
double sumPtTrackIsoForEgammaSC_endcap = iConfig.getParameter<double>("sumPtTrackIsoForEgammaSC_endcap");
double coneTrackIsoForEgammaSC = iConfig.getParameter<double>("coneTrackIsoForEgammaSC");
unsigned int nTrackIsoForEgammaSC = iConfig.getParameter<unsigned int>("nTrackIsoForEgammaSC");
// register products
produces<reco::PFCandidateCollection>();
produces<reco::PFCandidateCollection>("CleanedHF");
produces<reco::PFCandidateCollection>("CleanedCosmicsMuons");
produces<reco::PFCandidateCollection>("CleanedTrackerAndGlobalMuons");
produces<reco::PFCandidateCollection>("CleanedFakeMuons");
produces<reco::PFCandidateCollection>("CleanedPunchThroughMuons");
produces<reco::PFCandidateCollection>("CleanedPunchThroughNeutralHadrons");
produces<reco::PFCandidateCollection>("AddedMuonsAndHadrons");
if (usePFElectrons_) {
produces<reco::PFCandidateCollection>(electronOutputCol_);
produces<reco::PFCandidateElectronExtraCollection>(electronExtraOutputCol_);
}
double nSigmaECAL
= iConfig.getParameter<double>("pf_nsigma_ECAL");
double nSigmaHCAL
= iConfig.getParameter<double>("pf_nsigma_HCAL");
//PFElectrons Configuration
double mvaEleCut
= iConfig.getParameter<double>("pf_electron_mvaCut");
string mvaWeightFileEleID
= iConfig.getParameter<string>("pf_electronID_mvaWeightFile");
bool applyCrackCorrectionsForElectrons
= iConfig.getParameter<bool>("pf_electronID_crackCorrection");
string path_mvaWeightFileEleID;
if(usePFElectrons_)
{
path_mvaWeightFileEleID = edm::FileInPath ( mvaWeightFileEleID.c_str() ).fullPath();
}
//PFPhoton Configuration
string path_mvaWeightFileConvID;
string mvaWeightFileConvID;
double mvaConvCut=-99.;
if(usePFPhotons_)
{
mvaWeightFileConvID =iConfig.getParameter<string>("pf_convID_mvaWeightFile");
mvaConvCut = iConfig.getParameter<double>("pf_conv_mvaCut");
path_mvaWeightFileConvID = edm::FileInPath ( mvaWeightFileConvID.c_str() ).fullPath();
}
//Secondary tracks and displaced vertices parameters
bool rejectTracks_Bad
= iConfig.getParameter<bool>("rejectTracks_Bad");
bool rejectTracks_Step45
= iConfig.getParameter<bool>("rejectTracks_Step45");
bool usePFNuclearInteractions
= iConfig.getParameter<bool>("usePFNuclearInteractions");
bool usePFConversions
= iConfig.getParameter<bool>("usePFConversions");
bool usePFDecays
= iConfig.getParameter<bool>("usePFDecays");
double dptRel_DispVtx
= iConfig.getParameter<double>("dptRel_DispVtx");
edm::ParameterSet iCfgCandConnector
= iConfig.getParameter<edm::ParameterSet>("iCfgCandConnector");
// fToRead = iConfig.getUntrackedParameter<vector<string> >("toRead");
useCalibrationsFromDB_
= iConfig.getParameter<bool>("useCalibrationsFromDB");
boost::shared_ptr<PFEnergyCalibration>
calibration( new PFEnergyCalibration() );
int algoType
= iConfig.getParameter<unsigned>("algoType");
switch(algoType) {
case 0:
pfAlgo_.reset( new PFAlgo);
break;
default:
assert(0);
}
pfAlgo_->setParameters( nSigmaECAL,
nSigmaHCAL,
calibration,
thepfEnergyCalibrationHF);
//PFElectrons: call the method setpfeleparameters
pfAlgo_->setPFEleParameters(mvaEleCut,
path_mvaWeightFileEleID,
usePFElectrons_,
thePFSCEnergyCalibration,
calibration,
sumEtEcalIsoForEgammaSC_barrel,
sumEtEcalIsoForEgammaSC_endcap,
coneEcalIsoForEgammaSC,
sumPtTrackIsoForEgammaSC_barrel,
sumPtTrackIsoForEgammaSC_endcap,
nTrackIsoForEgammaSC,
coneTrackIsoForEgammaSC,
applyCrackCorrectionsForElectrons,
usePFSCEleCalib,
useEGammaElectrons_,
useEGammaSupercluster);
// pfAlgo_->setPFConversionParameters(usePFConversions);
// PFPhotons:
pfAlgo_->setPFPhotonParameters(usePFPhotons_,
path_mvaWeightFileConvID,
mvaConvCut,
calibration);
//Secondary tracks and displaced vertices parameters
pfAlgo_->setDisplacedVerticesParameters(rejectTracks_Bad,
rejectTracks_Step45,
usePFNuclearInteractions,
usePFConversions,
usePFDecays,
dptRel_DispVtx);
if (usePFNuclearInteractions)
pfAlgo_->setCandConnectorParameters( iCfgCandConnector );
// Muon parameters
std::vector<double> muonHCAL
= iConfig.getParameter<std::vector<double> >("muon_HCAL");
std::vector<double> muonECAL
= iConfig.getParameter<std::vector<double> >("muon_ECAL");
assert ( muonHCAL.size() == 2 && muonECAL.size() == 2 );
// Fake track parameters
double nSigmaTRACK
= iConfig.getParameter<double>("nsigma_TRACK");
double ptError
= iConfig.getParameter<double>("pt_Error");
std::vector<double> factors45
= iConfig.getParameter<std::vector<double> >("factors_45");
assert ( factors45.size() == 2 );
bool usePFMuonMomAssign
= iConfig.getParameter<bool>("usePFMuonMomAssign");
// Set muon and fake track parameters
pfAlgo_->setPFMuonAndFakeParameters(muonHCAL,
muonECAL,
nSigmaTRACK,
ptError,
factors45,
usePFMuonMomAssign);
//Post cleaning of the HF
bool postHFCleaning
= iConfig.getParameter<bool>("postHFCleaning");
double minHFCleaningPt
= iConfig.getParameter<double>("minHFCleaningPt");
double minSignificance
= iConfig.getParameter<double>("minSignificance");
double maxSignificance
= iConfig.getParameter<double>("maxSignificance");
double minSignificanceReduction
= iConfig.getParameter<double>("minSignificanceReduction");
double maxDeltaPhiPt
= iConfig.getParameter<double>("maxDeltaPhiPt");
double minDeltaMet
= iConfig.getParameter<double>("minDeltaMet");
// Set post HF cleaning muon parameters
pfAlgo_->setPostHFCleaningParameters(postHFCleaning,
minHFCleaningPt,
minSignificance,
maxSignificance,
minSignificanceReduction,
maxDeltaPhiPt,
minDeltaMet);
// Input tags for HF cleaned rechits
inputTagCleanedHF_
= iConfig.getParameter< std::vector<edm::InputTag> >("cleanedHF");
//MIKE: Vertex Parameters
vertices_ = iConfig.getParameter<edm::InputTag>("vertexCollection");
useVerticesForNeutral_ = iConfig.getParameter<bool>("useVerticesForNeutral");
verbose_ =
iConfig.getUntrackedParameter<bool>("verbose",false);
bool debug_ =
iConfig.getUntrackedParameter<bool>("debug",false);
pfAlgo_->setDebug( debug_ );
}
| PFProducer::~PFProducer | ( | ) |
Definition at line 310 of file PFProducer.cc.
{}
| void PFProducer::beginJob | ( | void | ) | [virtual] |
| void PFProducer::beginRun | ( | edm::Run & | run, |
| const edm::EventSetup & | es | ||
| ) | [virtual] |
Reimplemented from edm::EDProducer.
Definition at line 317 of file PFProducer.cc.
References edm::EventSetup::get(), and edm::ESHandle< T >::product().
{
/*
static map<string, PerformanceResult::ResultType> functType;
functType["PFfa_BARREL"] = PerformanceResult::PFfa_BARREL;
functType["PFfa_ENDCAP"] = PerformanceResult::PFfa_ENDCAP;
functType["PFfb_BARREL"] = PerformanceResult::PFfb_BARREL;
functType["PFfb_ENDCAP"] = PerformanceResult::PFfb_ENDCAP;
functType["PFfc_BARREL"] = PerformanceResult::PFfc_BARREL;
functType["PFfc_ENDCAP"] = PerformanceResult::PFfc_ENDCAP;
functType["PFfaEta_BARREL"] = PerformanceResult::PFfaEta_BARREL;
functType["PFfaEta_ENDCAP"] = PerformanceResult::PFfaEta_ENDCAP;
functType["PFfbEta_BARREL"] = PerformanceResult::PFfbEta_BARREL;
functType["PFfbEta_ENDCAP"] = PerformanceResult::PFfbEta_ENDCAP;
*/
if ( useCalibrationsFromDB_ ) {
// Read the PFCalibration functions from the global tags.
edm::ESHandle<PerformancePayload> perfH;
es.get<PFCalibrationRcd>().get(perfH);
const PerformancePayloadFromTFormula *pfCalibrations = static_cast< const PerformancePayloadFromTFormula *>(perfH.product());
pfAlgo_->thePFEnergyCalibration()->setCalibrationFunctions(pfCalibrations);
}
/*
for(vector<string>::const_iterator name = fToRead.begin(); name != fToRead.end(); ++name) {
cout << "Function: " << *name << endl;
PerformanceResult::ResultType fType = functType[*name];
pfCalibrations->printFormula(fType);
// evaluate it @ 10 GeV
float energy = 10.;
BinningPointByMap point;
point.insert(BinningVariables::JetEt, energy);
if(pfCalibrations->isInPayload(fType, point)) {
float value = pfCalibrations->getResult(fType, point);
cout << " Energy before:: " << energy << " after: " << value << endl;
} else cout << "outside limits!" << endl;
}
*/
}
| void PFProducer::produce | ( | edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [virtual] |
Implements edm::EDProducer.
Definition at line 372 of file PFProducer.cc.
References runregparse::blocks, gather_cfg::cout, edm::EventID::event(), Exception, newFWLiteAna::found, edm::Event::getByLabel(), edm::EventBase::id(), edm::HandleBase::isValid(), LogDebug, patZpeak::muons, edm::Event::put(), and edm::EventID::run().
{
LogDebug("PFProducer")<<"START event: "
<<iEvent.id().event()
<<" in run "<<iEvent.id().run()<<endl;
// Get The vertices from the event
// and assign dynamic vertex parameters
edm::Handle<reco::VertexCollection> vertices;
bool gotVertices = iEvent.getByLabel(vertices_,vertices);
if(!gotVertices) {
ostringstream err;
err<<"Cannot find vertices for this event.Continuing Without them ";
LogError("PFProducer")<<err.str()<<endl;
}
//Assign the PFAlgo Parameters
pfAlgo_->setPFVertexParameters(useVerticesForNeutral_,*vertices);
// get the collection of blocks
Handle< reco::PFBlockCollection > blocks;
LogDebug("PFProducer")<<"getting blocks"<<endl;
bool found = iEvent.getByLabel( inputTagBlocks_, blocks );
if(!found ) {
ostringstream err;
err<<"cannot find blocks: "<<inputTagBlocks_;
LogError("PFProducer")<<err.str()<<endl;
throw cms::Exception( "MissingProduct", err.str());
}
// get the collection of muons
Handle< reco::MuonCollection > muons;
if ( postMuonCleaning_ ) {
LogDebug("PFProducer")<<"getting muons"<<endl;
found = iEvent.getByLabel( inputTagMuons_, muons );
if(!found) {
ostringstream err;
err<<"cannot find muons: "<<inputTagMuons_;
LogError("PFProducer")<<err.str()<<endl;
throw cms::Exception( "MissingProduct", err.str());
}
}
if (useEGammaElectrons_) {
Handle < reco::GsfElectronCollection > egelectrons;
LogDebug("PFProducer")<<" Reading e/gamma electrons activated "<<endl;
found = iEvent.getByLabel( inputTagEgammaElectrons_, egelectrons );
if(!found) {
ostringstream err;
err<<"cannot find electrons: "<<inputTagEgammaElectrons_;
LogError("PFProducer")<<err.str()<<endl;
throw cms::Exception( "MissingProduct", err.str());
}
pfAlgo_->setEGElectronCollection(*egelectrons);
}
LogDebug("PFProducer")<<"particle flow is starting"<<endl;
assert( blocks.isValid() );
pfAlgo_->reconstructParticles( blocks );
if(verbose_) {
ostringstream str;
str<<(*pfAlgo_)<<endl;
// cout << (*pfAlgo_) << endl;
LogInfo("PFProducer") <<str.str()<<endl;
}
if ( postMuonCleaning_ )
pfAlgo_->postMuonCleaning( muons, *vertices );
// Florian 5/01/2011
// Save the PFElectron Extra Collection First as to be able to create valid References
if(usePFElectrons_) {
auto_ptr< reco::PFCandidateElectronExtraCollection >
pOutputElectronCandidateExtraCollection( pfAlgo_->transferElectronExtra() );
const edm::OrphanHandle<reco::PFCandidateElectronExtraCollection > electronExtraProd=
iEvent.put(pOutputElectronCandidateExtraCollection,electronExtraOutputCol_);
pfAlgo_->setElectronExtraRef(electronExtraProd);
}
// Save cosmic cleaned muon candidates
auto_ptr< reco::PFCandidateCollection >
pCosmicsMuonCleanedCandidateCollection( pfAlgo_->transferCosmicsMuonCleanedCandidates() );
// Save tracker/global cleaned muon candidates
auto_ptr< reco::PFCandidateCollection >
pTrackerAndGlobalCleanedMuonCandidateCollection( pfAlgo_->transferCleanedTrackerAndGlobalMuonCandidates() );
// Save fake cleaned muon candidates
auto_ptr< reco::PFCandidateCollection >
pFakeCleanedMuonCandidateCollection( pfAlgo_->transferFakeMuonCleanedCandidates() );
// Save punch-through cleaned muon candidates
auto_ptr< reco::PFCandidateCollection >
pPunchThroughMuonCleanedCandidateCollection( pfAlgo_->transferPunchThroughMuonCleanedCandidates() );
// Save punch-through cleaned neutral hadron candidates
auto_ptr< reco::PFCandidateCollection >
pPunchThroughHadronCleanedCandidateCollection( pfAlgo_->transferPunchThroughHadronCleanedCandidates() );
// Save added muon candidates
auto_ptr< reco::PFCandidateCollection >
pAddedMuonCandidateCollection( pfAlgo_->transferAddedMuonCandidates() );
// Check HF overcleaning
reco::PFRecHitCollection hfCopy;
for ( unsigned ihf=0; ihf<inputTagCleanedHF_.size(); ++ihf ) {
Handle< reco::PFRecHitCollection > hfCleaned;
bool foundHF = iEvent.getByLabel( inputTagCleanedHF_[ihf], hfCleaned );
if (!foundHF) continue;
for ( unsigned jhf=0; jhf<(*hfCleaned).size(); ++jhf ) {
hfCopy.push_back( (*hfCleaned)[jhf] );
}
}
pfAlgo_->checkCleaning( hfCopy );
// Save recovered HF candidates
auto_ptr< reco::PFCandidateCollection >
pCleanedCandidateCollection( pfAlgo_->transferCleanedCandidates() );
// Save the final PFCandidate collection
auto_ptr< reco::PFCandidateCollection >
pOutputCandidateCollection( pfAlgo_->transferCandidates() );
LogDebug("PFProducer")<<"particle flow: putting products in the event"<<endl;
if ( verbose_ ) std::cout <<"particle flow: putting products in the event. Here the full list"<<endl;
int nC=0;
for( reco::PFCandidateCollection::const_iterator itCand = (*pOutputCandidateCollection).begin(); itCand != (*pOutputCandidateCollection).end(); itCand++) {
nC++;
if (verbose_ ) std::cout << nC << ")" << (*itCand).particleId() << std::endl;
}
// Write in the event
iEvent.put(pOutputCandidateCollection);
iEvent.put(pCleanedCandidateCollection,"CleanedHF");
if ( postMuonCleaning_ ) {
iEvent.put(pCosmicsMuonCleanedCandidateCollection,"CleanedCosmicsMuons");
iEvent.put(pTrackerAndGlobalCleanedMuonCandidateCollection,"CleanedTrackerAndGlobalMuons");
iEvent.put(pFakeCleanedMuonCandidateCollection,"CleanedFakeMuons");
iEvent.put(pPunchThroughMuonCleanedCandidateCollection,"CleanedPunchThroughMuons");
iEvent.put(pPunchThroughHadronCleanedCandidateCollection,"CleanedPunchThroughNeutralHadrons");
iEvent.put(pAddedMuonCandidateCollection,"AddedMuonsAndHadrons");
}
if(usePFElectrons_)
{
auto_ptr< reco::PFCandidateCollection >
pOutputElectronCandidateCollection( pfAlgo_->transferElectronCandidates() );
iEvent.put(pOutputElectronCandidateCollection,electronOutputCol_);
}
}
std::string PFProducer::electronExtraOutputCol_ [private] |
Definition at line 46 of file PFProducer.h.
std::string PFProducer::electronOutputCol_ [private] |
Definition at line 45 of file PFProducer.h.
edm::InputTag PFProducer::inputTagBlocks_ [private] |
Definition at line 40 of file PFProducer.h.
std::vector<edm::InputTag> PFProducer::inputTagCleanedHF_ [private] |
Definition at line 44 of file PFProducer.h.
Definition at line 43 of file PFProducer.h.
edm::InputTag PFProducer::inputTagMuons_ [private] |
Definition at line 41 of file PFProducer.h.
std::auto_ptr<PFAlgo> PFProducer::pfAlgo_ [private] |
particle flow algorithm
Definition at line 73 of file PFProducer.h.
bool PFProducer::postMuonCleaning_ [private] |
Definition at line 52 of file PFProducer.h.
bool PFProducer::useCalibrationsFromDB_ [private] |
Definition at line 67 of file PFProducer.h.
bool PFProducer::useEGammaElectrons_ [private] |
Definition at line 61 of file PFProducer.h.
bool PFProducer::usePFElectrons_ [private] |
Definition at line 55 of file PFProducer.h.
bool PFProducer::usePFPhotons_ [private] |
Definition at line 58 of file PFProducer.h.
bool PFProducer::useVerticesForNeutral_ [private] |
Definition at line 64 of file PFProducer.h.
bool PFProducer::verbose_ [private] |
verbose ?
Definition at line 49 of file PFProducer.h.
edm::InputTag PFProducer::vertices_ [private] |
Definition at line 42 of file PFProducer.h.
1.7.3