|
|
|
#include <PhotonCoreProducer.h>
Definition at line 34 of file PhotonCoreProducer.h.
| PhotonCoreProducer::PhotonCoreProducer | ( | const edm::ParameterSet & | ps | ) |
Definition at line 20 of file PhotonCoreProducer.cc.
References conf_, conversionProducer_, edm::ParameterSet::getParameter(), minSCEt_, PhotonCoreCollection_, pixelSeedProducer_, risolveAmbiguity_, scHybridBarrelProducer_, and scIslandEndcapProducer_.
: conf_(config) { // use onfiguration file to setup input/output collection names scHybridBarrelProducer_ = conf_.getParameter<edm::InputTag>("scHybridBarrelProducer"); scIslandEndcapProducer_ = conf_.getParameter<edm::InputTag>("scIslandEndcapProducer"); conversionProducer_ = conf_.getParameter<edm::InputTag>("conversionProducer"); PhotonCoreCollection_ = conf_.getParameter<std::string>("photonCoreCollection"); pixelSeedProducer_ = conf_.getParameter<std::string>("pixelSeedProducer"); minSCEt_ = conf_.getParameter<double>("minSCEt"); risolveAmbiguity_ = conf_.getParameter<bool>("risolveConversionAmbiguity"); // Register the product produces< reco::PhotonCoreCollection >(PhotonCoreCollection_); }
| PhotonCoreProducer::~PhotonCoreProducer | ( | ) |
Definition at line 39 of file PhotonCoreProducer.cc.
{}
| void PhotonCoreProducer::beginRun | ( | edm::Run & | r, |
| edm::EventSetup const & | es | ||
| ) | [virtual] |
| void PhotonCoreProducer::endRun | ( | edm::Run & | r, |
| edm::EventSetup const & | theEventSetup | ||
| ) | [virtual] |
| void PhotonCoreProducer::fillPhotonCollection | ( | edm::Event & | evt, |
| edm::EventSetup const & | es, | ||
| const edm::Handle< reco::SuperClusterCollection > & | scHandle, | ||
| const edm::Handle< reco::ConversionCollection > & | conversionHandle, | ||
| const edm::Handle< reco::ElectronSeedCollection > & | pixelSeeds, | ||
| reco::PhotonCoreCollection & | outputCollection, | ||
| int & | iSC | ||
| ) | [private] |
Definition at line 131 of file PhotonCoreProducer.cc.
References reco::PhotonCore::addConversion(), reco::PhotonCore::addElectronPixelSeed(), edm::Ref< C, T, F >::id(), edm::Ref< C, T, F >::isNonnull(), edm::Ref< C, T, F >::key(), minSCEt_, risolveAmbiguity_, solveAmbiguity(), validConversions_, and validPixelSeeds_.
Referenced by produce().
{
reco::ElectronSeedCollection::const_iterator pixelSeedItr;
for(unsigned int lSC=0; lSC < scHandle->size(); lSC++) {
// get SuperClusterRef
reco::SuperClusterRef scRef(reco::SuperClusterRef(scHandle, lSC));
iSC++;
//const reco::SuperCluster* pClus=&(*scRef);
// SC energy preselection
if (scRef->energy()/cosh(scRef->eta()) <= minSCEt_) continue;
reco::PhotonCore newCandidate(scRef);
if ( validConversions_) {
if ( risolveAmbiguity_ ) {
reco::ConversionRef bestRef=solveAmbiguity( conversionHandle , scRef);
if (bestRef.isNonnull() ) newCandidate.addConversion(bestRef);
} else {
for( unsigned int icp = 0; icp < conversionHandle->size(); icp++) {
reco::ConversionRef cpRef(reco::ConversionRef(conversionHandle,icp));
if ( !cpRef->caloCluster().size()) continue;
if (!( scRef.id() == cpRef->caloCluster()[0].id() && scRef.key() == cpRef->caloCluster()[0].key() )) continue;
if ( !cpRef->isConverted() ) continue;
newCandidate.addConversion(cpRef);
}
} // solve or not the ambiguity of many conversion candidates
}
if ( validPixelSeeds_) {
for( unsigned int icp = 0; icp < pixelSeedHandle->size(); icp++) {
reco::ElectronSeedRef cpRef(reco::ElectronSeedRef(pixelSeedHandle,icp));
if ( ! cpRef->isEcalDriven() ) continue;
if (!( scRef.id() == cpRef->caloCluster().id() && scRef.key() == cpRef->caloCluster().key() )) continue;
newCandidate.addElectronPixelSeed(cpRef);
}
}
outputPhotonCoreCollection.push_back(newCandidate);
}
}
| void PhotonCoreProducer::produce | ( | edm::Event & | evt, |
| const edm::EventSetup & | es | ||
| ) | [virtual] |
Implements edm::EDProducer.
Definition at line 54 of file PhotonCoreProducer.cc.
References printConversionInfo::conversionHandle, conversionProducer_, fillPhotonCollection(), edm::Event::getByLabel(), edm::InputTag::label(), PhotonCoreCollection_, pixelSeedProducer_, edm::Event::put(), scHybridBarrelProducer_, scIslandEndcapProducer_, validConversions_, and validPixelSeeds_.
{
using namespace edm;
// nEvt_++;
reco::PhotonCoreCollection outputPhotonCoreCollection;
std::auto_ptr< reco::PhotonCoreCollection > outputPhotonCoreCollection_p(new reco::PhotonCoreCollection);
// Get the Barrel Super Cluster collection
bool validBarrelSCHandle=true;
Handle<reco::SuperClusterCollection> scBarrelHandle;
theEvent.getByLabel(scHybridBarrelProducer_,scBarrelHandle);
if (!scBarrelHandle.isValid()) {
edm::LogError("PhotonCoreProducer") << "Error! Can't get the product "<<scHybridBarrelProducer_.label();
validBarrelSCHandle=false;
}
// Get the Endcap Super Cluster collection
bool validEndcapSCHandle=true;
Handle<reco::SuperClusterCollection> scEndcapHandle;
theEvent.getByLabel(scIslandEndcapProducer_,scEndcapHandle);
if (!scEndcapHandle.isValid()) {
edm::LogError("PhotonCoreProducer") << "Error! Can't get the product "<<scIslandEndcapProducer_.label();
validEndcapSCHandle=false;
}
validConversions_=true;
edm::Handle<reco::ConversionCollection> conversionHandle;
theEvent.getByLabel(conversionProducer_, conversionHandle);
if (!conversionHandle.isValid()) {
//edm::LogError("PhotonCoreProducer") << "Error! Can't get the product "<< conversionProducer_.label() << "\n" ;
validConversions_=false;
}
// Get ElectronPixelSeeds
validPixelSeeds_=true;
Handle<reco::ElectronSeedCollection> pixelSeedHandle;
reco::ElectronSeedCollection pixelSeeds;
theEvent.getByLabel(pixelSeedProducer_, pixelSeedHandle);
if (!pixelSeedHandle.isValid()) {
validPixelSeeds_=false;
}
// if ( validPixelSeeds_) pixelSeeds = *(pixelSeedHandle.product());
int iSC=0; // index in photon collection
// Loop over barrel and endcap SC collections and fill the photon collection
if ( validBarrelSCHandle) fillPhotonCollection(theEvent,
theEventSetup,
scBarrelHandle,
conversionHandle,
pixelSeedHandle,
outputPhotonCoreCollection,
iSC);
if ( validEndcapSCHandle) fillPhotonCollection(theEvent,
theEventSetup,
scEndcapHandle,
conversionHandle,
pixelSeedHandle,
outputPhotonCoreCollection,
iSC);
// put the product in the event
edm::LogInfo("PhotonCoreProducer") << " Put in the event " << iSC << " Photon Candidates \n";
outputPhotonCoreCollection_p->assign(outputPhotonCoreCollection.begin(),outputPhotonCoreCollection.end());
theEvent.put( outputPhotonCoreCollection_p, PhotonCoreCollection_);
}
| reco::ConversionRef PhotonCoreProducer::solveAmbiguity | ( | const edm::Handle< reco::ConversionCollection > & | conversionHandle, |
| reco::SuperClusterRef & | sc | ||
| ) | [private] |
only one track reconstructed. Pick the one with best E/P
Definition at line 193 of file PhotonCoreProducer.cc.
References edm::Ref< C, T, F >::id(), edm::Ref< C, T, F >::key(), AlCaHLTBitMon_ParallelJobs::p, mathSSE::sqrt(), and testEve_cfg::tracks.
Referenced by fillPhotonCollection().
{
std::multimap<reco::ConversionRef, double > convMap;
for ( unsigned int icp=0; icp< conversionHandle->size(); icp++) {
reco::ConversionRef cpRef(reco::ConversionRef(conversionHandle,icp));
if (!( scRef.id() == cpRef->caloCluster()[0].id() && scRef.key() == cpRef->caloCluster()[0].key() )) continue;
if ( !cpRef->isConverted() ) continue;
double like = cpRef->MVAout();
convMap.insert ( std::make_pair(cpRef,like) );
}
std::multimap<reco::ConversionRef, double >::iterator iMap;
double max_lh = -1.;
reco::ConversionRef bestRef;
// std::cout << " Pick up the best conv " << std::endl;
for (iMap=convMap.begin(); iMap!=convMap.end(); iMap++) {
double like = iMap->second;
if (like > max_lh) {
max_lh = like;
bestRef=iMap->first;
}
}
//std::cout << " Best conv like " << max_lh << std::endl;
float ep=0;
if ( max_lh <0 ) {
// std::cout << " Candidates with only one track " << std::endl;
float epMin=999;
for (iMap=convMap.begin(); iMap!=convMap.end(); iMap++) {
reco::ConversionRef convRef=iMap->first;
// std::vector<reco::TrackRef> tracks = convRef->tracks();
const std::vector<edm::RefToBase<reco::Track> > tracks = convRef->tracks();
float px=tracks[0]->innerMomentum().x();
float py=tracks[0]->innerMomentum().y();
float pz=tracks[0]->innerMomentum().z();
float p=sqrt(px*px+py*py+pz*pz);
ep=fabs(1.-convRef->caloCluster()[0]->energy()/p);
// std::cout << " 1-E/P = " << ep << std::endl;
if ( ep<epMin) {
epMin=ep;
bestRef=iMap->first;
}
}
// std::cout << " Best conv 1-E/P " << ep << std::endl;
}
return bestRef;
}
edm::ParameterSet PhotonCoreProducer::conf_ [private] |
Definition at line 70 of file PhotonCoreProducer.h.
Referenced by PhotonCoreProducer().
Definition at line 65 of file PhotonCoreProducer.h.
Referenced by PhotonCoreProducer(), and produce().
double PhotonCoreProducer::minSCEt_ [private] |
Definition at line 67 of file PhotonCoreProducer.h.
Referenced by fillPhotonCollection(), and PhotonCoreProducer().
std::string PhotonCoreProducer::PhotonCoreCollection_ [private] |
Definition at line 59 of file PhotonCoreProducer.h.
Referenced by PhotonCoreProducer(), and produce().
std::string PhotonCoreProducer::pixelSeedProducer_ [private] |
Definition at line 69 of file PhotonCoreProducer.h.
Referenced by PhotonCoreProducer(), and produce().
bool PhotonCoreProducer::risolveAmbiguity_ [private] |
Definition at line 72 of file PhotonCoreProducer.h.
Referenced by fillPhotonCollection(), and PhotonCoreProducer().
Definition at line 62 of file PhotonCoreProducer.h.
Definition at line 60 of file PhotonCoreProducer.h.
Referenced by PhotonCoreProducer(), and produce().
Definition at line 63 of file PhotonCoreProducer.h.
Definition at line 61 of file PhotonCoreProducer.h.
Referenced by PhotonCoreProducer(), and produce().
bool PhotonCoreProducer::validConversions_ [private] |
Definition at line 68 of file PhotonCoreProducer.h.
Referenced by fillPhotonCollection(), and produce().
bool PhotonCoreProducer::validPixelSeeds_ [private] |
Definition at line 71 of file PhotonCoreProducer.h.
Referenced by fillPhotonCollection(), and produce().
1.7.3