Inheritance diagram for AlCaIsoTracksProducer:

Public Member Functions
	AlCaIsoTracksProducer (const edm::ParameterSet &)

virtual void	produce (edm::Event &, const edm::EventSetup &)

	~AlCaIsoTracksProducer ()

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
std::vector< ConsumesInfo >	consumesInfo () const

	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

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
virtual void	beginJob ()

virtual void	beginRun (edm::Run const &, edm::EventSetup const &)

virtual void	endJob ()

virtual void	endRun (edm::Run const &, edm::EventSetup const &)

reco::HcalIsolatedTrackCandidateCollection *	select (edm::Handle< edm::TriggerResults > &triggerResults, const std::vector< std::string > &triggerNames_, edm::Handle< reco::TrackCollection > &trkCollection, math::XYZPoint &leadPV, edm::Handle< EcalRecHitCollection > &barrelRecHitsHandle, edm::Handle< EcalRecHitCollection > &endcapRecHitsHandle, edm::Handle< HBHERecHitCollection > &hbhe)

void	setPtEtaPhi (std::vector< edm::Ref< l1extra::L1JetParticleCollection > > &objref, double &ptL1, double &etaL1, double &phiL1)

Private Attributes
double	a_charIsoR

double	a_coneR

double	a_mipR

const MagneticField *	bField

double	eEcalMax_

double	eIsolation_

double	etaL1

const CaloGeometry *	geo

HLTConfigProvider	hltConfig_

std::vector< std::string >	HLTNames

std::string	l1Filter

std::string	l2Filter

std::string	l3Filter

edm::InputTag	labelBS_

edm::InputTag	labelEB_

edm::InputTag	labelEE_

edm::InputTag	labelGenTrack_

edm::InputTag	labelHBHE_

edm::InputTag	labelHltGT_

edm::InputTag	labelRecVtx_

edm::InputTag	labelTriggerEvent_

edm::InputTag	labelTriggerResults_

int	nAll

int	nGood

int	nRun

double	phiL1

std::string	processName

double	ptL1

double	pTrackMin_

spr::trackSelectionParameters	selectionParameters

std::string	theTrackQuality

edm::EDGetTokenT< reco::BeamSpot >	tok_bs_

edm::EDGetTokenT < EcalRecHitCollection >	tok_EB_

edm::EDGetTokenT < EcalRecHitCollection >	tok_EE_

edm::EDGetTokenT < reco::TrackCollection >	tok_genTrack_

edm::EDGetTokenT < HBHERecHitCollection >	tok_hbhe_

edm::EDGetTokenT < trigger::TriggerFilterObjectWithRefs >	tok_hltGT_

edm::EDGetTokenT < reco::VertexCollection >	tok_recVtx_

edm::EDGetTokenT < trigger::TriggerEvent >	tok_trigEvt_

edm::EDGetTokenT < edm::TriggerResults >	tok_trigRes_

std::vector< int >	trigKount

std::vector< std::string >	trigNames

std::vector< int >	trigPass

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

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 79 of file AlCaIsoTracksProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 124 of file AlCaIsoTracksProducer.cc.

                                                                            :
   nRun(0), nAll(0), nGood(0), ptL1(0), etaL1(0), phiL1(0) {
   //Get the run parameters
   const double isolationRadius(28.9);
   trigNames                           = iConfig.getParameter<std::vector<std::string> >("Triggers");
   theTrackQuality                     = iConfig.getParameter<std::string>("TrackQuality");
   processName                         = iConfig.getParameter<std::string>("ProcessName");
   l1Filter                            = iConfig.getParameter<std::string>("L1Filter");
   l2Filter                            = iConfig.getParameter<std::string>("L2Filter");
   l3Filter                            = iConfig.getParameter<std::string>("L3Filter");
   selectionParameters.minPt           = iConfig.getParameter<double>("MinTrackPt");
   selectionParameters.minQuality      = reco::TrackBase::qualityByName(theTrackQuality);
   selectionParameters.maxDxyPV        = iConfig.getParameter<double>("MaxDxyPV");
   selectionParameters.maxDzPV         = iConfig.getParameter<double>("MaxDzPV");
   selectionParameters.maxChi2         = iConfig.getParameter<double>("MaxChi2");
   selectionParameters.maxDpOverP      = iConfig.getParameter<double>("MaxDpOverP");
   selectionParameters.minOuterHit     = iConfig.getParameter<int>("MinOuterHit");
   selectionParameters.minLayerCrossed = iConfig.getParameter<int>("MinLayerCrossed");
   selectionParameters.maxInMiss       = iConfig.getParameter<int>("MaxInMiss");
   selectionParameters.maxOutMiss      = iConfig.getParameter<int>("MaxOutMiss");
   a_coneR                             = iConfig.getParameter<double>("ConeRadius");
   a_charIsoR                          = a_coneR + isolationRadius;
   a_mipR                              = iConfig.getParameter<double>("ConeRadiusMIP");
   pTrackMin_                          = iConfig.getParameter<double>("MinimumTrackP");
   eEcalMax_                           = iConfig.getParameter<double>("MaximumEcalEnergy");
   eIsolation_                         = iConfig.getParameter<double>("IsolationEnergy");
   labelGenTrack_                      = iConfig.getParameter<edm::InputTag>("TrackLabel");
   labelRecVtx_                        = iConfig.getParameter<edm::InputTag>("VertexLabel");
   labelBS_                            = iConfig.getParameter<edm::InputTag>("BeamSpotLabel");
   labelEB_                            = iConfig.getParameter<edm::InputTag>("EBRecHitLabel");
   labelEE_                            = iConfig.getParameter<edm::InputTag>("EERecHitLabel");
   labelHBHE_                          = iConfig.getParameter<edm::InputTag>("HBHERecHitLabel");
   labelHltGT_                         = iConfig.getParameter<edm::InputTag>("L1GTSeedLabel");
   labelTriggerEvent_                  = iConfig.getParameter<edm::InputTag>("TriggerEventLabel");
   labelTriggerResults_                = iConfig.getParameter<edm::InputTag>("TriggerResultLabel");
 
   // define tokens for access
   tok_hltGT_    = consumes<trigger::TriggerFilterObjectWithRefs>(labelHltGT_);
   tok_trigEvt_  = consumes<trigger::TriggerEvent>(labelTriggerEvent_);
   tok_trigRes_  = consumes<edm::TriggerResults>(labelTriggerResults_);
   tok_genTrack_ = consumes<reco::TrackCollection>(labelGenTrack_);
   tok_recVtx_   = consumes<reco::VertexCollection>(labelRecVtx_);
   tok_bs_       = consumes<reco::BeamSpot>(labelBS_);
   tok_EB_       = consumes<EcalRecHitCollection>(labelEB_);
   tok_EE_       = consumes<EcalRecHitCollection>(labelEE_);
   tok_hbhe_     = consumes<HBHERecHitCollection>(labelHBHE_);
 
   edm::LogInfo("HcalIsoTrack") <<"Parameters read from config file \n" 
                    <<"\t minPt "           << selectionParameters.minPt   
                    <<"\t theTrackQuality " << theTrackQuality
                    <<"\t minQuality "      << selectionParameters.minQuality
                    <<"\t maxDxyPV "        << selectionParameters.maxDxyPV          
                    <<"\t maxDzPV "         << selectionParameters.maxDzPV          
                    <<"\t maxChi2 "         << selectionParameters.maxChi2          
                    <<"\t maxDpOverP "      << selectionParameters.maxDpOverP
                    <<"\t minOuterHit "     << selectionParameters.minOuterHit
                    <<"\t minLayerCrossed " << selectionParameters.minLayerCrossed
                    <<"\t maxInMiss "       << selectionParameters.maxInMiss
                    <<"\t maxOutMiss "      << selectionParameters.maxOutMiss
                    <<"\t a_coneR "         << a_coneR          
                    <<"\t a_charIsoR "      << a_charIsoR          
                    <<"\t a_mipR "          << a_mipR
                    <<"\t pTrackMin "       << pTrackMin_
                    <<"\t eEcalMax "        << eEcalMax_
                    <<"\t eIsolation "      << eIsolation_;
   edm::LogInfo("HcalIsoTrack") << "Process " << processName << " L1Filter:" 
                    << l1Filter << " L2Filter:" << l2Filter 
                    << " L3Filter:" << l3Filter;
   for (unsigned int k=0; k<trigNames.size(); ++k)
     edm::LogInfo("HcalIsoTrack") << "Trigger[" << k << "] " << trigNames[k];
 
   std::vector<int> dummy(trigNames.size(),0);
   trigKount = trigPass = dummy;
 
   //create also IsolatedPixelTrackCandidateCollection which contains isolation info and reference to primary track
   produces<reco::HcalIsolatedTrackCandidateCollection>("HcalIsolatedTrackCollection");
   produces<reco::TrackCollection>(labelGenTrack_.encode());
   produces<reco::VertexCollection>(labelRecVtx_.encode());
   produces<EcalRecHitCollection>(labelEB_.encode());
   produces<EcalRecHitCollection>(labelEE_.encode());
   produces<HBHERecHitCollection>(labelHBHE_.encode());
 }

AlCaIsoTracksProducer::~AlCaIsoTracksProducer ( )

Definition at line 208 of file AlCaIsoTracksProducer.cc.

208 { }

Member Function Documentation

void AlCaIsoTracksProducer::beginJob ( void )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 340 of file AlCaIsoTracksProducer.cc.

340 { }

void AlCaIsoTracksProducer::beginRun	(	edm::Run const &	iRun,
		edm::EventSetup const &	iSetup
	)

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 351 of file AlCaIsoTracksProducer.cc.

References bField, geo, edm::EventSetup::get(), hltConfig_, HLTConfigProvider::init(), nRun, processName, edm::ESHandle< class >::product(), and edm::RunBase::run().

                                                                                     {
   bool changed(false);
   edm::LogInfo("HcalIsoTrack") << "Run[" << nRun << "] " << iRun.run() 
                    << " hltconfig.init " << hltConfig_.init(iRun,iSetup,processName,changed);
 
   edm::ESHandle<MagneticField> bFieldH;
   iSetup.get<IdealMagneticFieldRecord>().get(bFieldH);
   bField = bFieldH.product();
   edm::ESHandle<CaloGeometry> pG;
   iSetup.get<CaloGeometryRecord>().get(pG);
   geo    = pG.product();
 }

void AlCaIsoTracksProducer::endJob ( void )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 342 of file AlCaIsoTracksProducer.cc.

References relval_steps::k, nAll, nGood, nRun, trigKount, trigNames, and trigPass.

                                    {
   edm::LogInfo("HcalIsoTrack") << "Finds " << nGood << " good tracks in " 
                    << nAll << " events from " << nRun << " runs";
   for (unsigned int k=0; k<trigNames.size(); ++k)
     edm::LogInfo("HcalIsoTrack") << "Trigger[" << k << "]: " << trigNames[k] 
                  << " Events " << trigKount[k] << " Passed " 
                  << trigPass[k];
 }

void AlCaIsoTracksProducer::endRun	(	edm::Run const &	iRun,
		edm::EventSetup const &
	)

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 364 of file AlCaIsoTracksProducer.cc.

References nRun, and edm::RunBase::run().

                                                                            {
   nRun++;
   edm::LogInfo("HcalIsoTrack") << "endRun[" << nRun << "] " << iRun.run();
 }

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

virtual

Implements edm::EDProducer.

Definition at line 211 of file AlCaIsoTracksProducer.cc.

References edm::EventBase::bunchCrossing(), edm::InputTag::encode(), etaL1, edm::EventID::event(), edm::Event::getByToken(), edm::EventBase::id(), edm::HandleBase::isValid(), edm::InputTag::label(), labelEB_, labelEE_, labelGenTrack_, labelHBHE_, labelRecVtx_, labelTriggerEvent_, labelTriggerResults_, LogDebug, edm::EventBase::luminosityBlock(), nAll, phiL1, edm::Handle< T >::product(), ptL1, edm::Event::put(), edm::EventID::run(), select(), setPtEtaPhi(), tok_bs_, tok_EB_, tok_EE_, tok_genTrack_, tok_hbhe_, tok_hltGT_, tok_recVtx_, tok_trigEvt_, tok_trigRes_, trigger::TriggerL1CenJet, trigger::TriggerL1ForJet, trigger::TriggerL1TauJet, edm::TriggerNames::triggerNames(), edm::Event::triggerNames(), triggerResults, x, detailsBasic3DVector::y, and detailsBasic3DVector::z.

Referenced by JSONExport.JsonExport::export(), HTMLExport.HTMLExport::export(), and HTMLExport.HTMLExportStatic::export().

                                                                                  {
 
   nAll++;
   LogDebug("HcalIsoTrack") << "Run " << iEvent.id().run() << " Event " 
                << iEvent.id().event() << " Luminosity " 
                << iEvent.luminosityBlock() << " Bunch " 
                << iEvent.bunchCrossing();
   
   //Step1: Get all the relevant containers
   trigger::TriggerEvent triggerEvent;
   edm::Handle<trigger::TriggerEvent> triggerEventHandle;
   iEvent.getByToken(tok_trigEvt_, triggerEventHandle);
   if (!triggerEventHandle.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelTriggerEvent_;
     return;
   }
   edm::Handle<edm::TriggerResults> triggerResults;
   iEvent.getByToken(tok_trigRes_, triggerResults);
   if (!triggerResults.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelTriggerResults_;
     return;
   }
 
   edm::Handle<reco::TrackCollection> trkCollection;
   iEvent.getByToken(tok_genTrack_, trkCollection);
   if (!trkCollection.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelGenTrack_;
     return;
   }
   reco::TrackCollection::const_iterator trkItr;
   edm::Handle<reco::VertexCollection> recVtxs;
   iEvent.getByToken(tok_recVtx_, recVtxs);  
   if (!trkCollection.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelGenTrack_;
     return;
   }
 
   edm::Handle<reco::BeamSpot> beamSpotH;
   iEvent.getByToken(tok_bs_, beamSpotH);
   math::XYZPoint leadPV(0,0,0);
   if (recVtxs->size()>0 && !((*recVtxs)[0].isFake())) {
     leadPV = math::XYZPoint((*recVtxs)[0].x(),(*recVtxs)[0].y(),
                 (*recVtxs)[0].z());
   } else if (beamSpotH.isValid()) {
     leadPV = beamSpotH->position();
   }
   LogDebug("HcalIsoTrack") << "Primary Vertex " << leadPV;
 
   edm::Handle<EcalRecHitCollection> barrelRecHitsHandle;
   iEvent.getByToken(tok_EB_, barrelRecHitsHandle);
   if (!barrelRecHitsHandle.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelEB_;
     return;
   }
   edm::Handle<EcalRecHitCollection> endcapRecHitsHandle;
   iEvent.getByToken(tok_EE_, endcapRecHitsHandle);
   if (!endcapRecHitsHandle.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelEE_;
     return;
   }
   edm::Handle<HBHERecHitCollection> hbhe;
   iEvent.getByToken(tok_hbhe_, hbhe);
   if (!hbhe.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelHBHE_;
     return;
   }
 
   //Get L1 trigger object
   ptL1 = etaL1 = phiL1 = 0;
   edm::Handle<trigger::TriggerFilterObjectWithRefs> l1trigobj;
   iEvent.getByToken(tok_hltGT_, l1trigobj);
 
   std::vector< edm::Ref<l1extra::L1JetParticleCollection> > l1tauobjref;
   l1trigobj->getObjects(trigger::TriggerL1TauJet, l1tauobjref);
   setPtEtaPhi(l1tauobjref,ptL1,etaL1,phiL1);
 
   std::vector< edm::Ref<l1extra::L1JetParticleCollection> > l1jetobjref;
   l1trigobj->getObjects(trigger::TriggerL1CenJet, l1jetobjref);
   setPtEtaPhi(l1jetobjref,ptL1,etaL1,phiL1);
 
   std::vector< edm::Ref<l1extra::L1JetParticleCollection> > l1forjetobjref;
   l1trigobj->getObjects(trigger::TriggerL1ForJet, l1forjetobjref);
   setPtEtaPhi(l1forjetobjref,ptL1,etaL1,phiL1);
 
   //For valid HLT record
   if (!triggerEventHandle.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Error! Can't get the product "
                     << labelTriggerEvent_.label() ;
   } else {
     trigger::TriggerEvent triggerEvent = *(triggerEventHandle.product());
     if (triggerResults.isValid()) {
       const edm::TriggerNames & triggerNames = iEvent.triggerNames(*triggerResults);
       const std::vector<std::string> & triggerNames_ = triggerNames.triggerNames();
       reco::HcalIsolatedTrackCandidateCollection* isotk = select(triggerResults,triggerNames_,trkCollection,leadPV, barrelRecHitsHandle, endcapRecHitsHandle, hbhe);
 
       if (isotk->size() > 0) {
     std::auto_ptr<reco::HcalIsolatedTrackCandidateCollection> outputHcalIsoTrackColl(isotk);
     std::auto_ptr<reco::TrackCollection>  outputTColl(new reco::TrackCollection);
     std::auto_ptr<reco::VertexCollection> outputVColl(new reco::VertexCollection);
     std::auto_ptr<EBRecHitCollection>     outputEBColl(new EBRecHitCollection);
     std::auto_ptr<EERecHitCollection>     outputEEColl(new EERecHitCollection);
     std::auto_ptr<HBHERecHitCollection>   outputHBHEColl(new HBHERecHitCollection);
     
     for (reco::TrackCollection::const_iterator track=trkCollection->begin(); track!=trkCollection->end(); ++track)
       outputTColl->push_back(*track);
 
     for (reco::VertexCollection::const_iterator vtx=recVtxs->begin(); vtx!=recVtxs->end(); ++vtx)
       outputVColl->push_back(*vtx);
       
     for (edm::SortedCollection<EcalRecHit>::const_iterator ehit=barrelRecHitsHandle->begin(); ehit!=barrelRecHitsHandle->end(); ++ehit)
       outputEBColl->push_back(*ehit);
     
     for (edm::SortedCollection<EcalRecHit>::const_iterator ehit=endcapRecHitsHandle->begin(); ehit!=endcapRecHitsHandle->end(); ++ehit)
       outputEEColl->push_back(*ehit);
     
     for (std::vector<HBHERecHit>::const_iterator hhit=hbhe->begin(); hhit!=hbhe->end(); ++hhit)
       outputHBHEColl->push_back(*hhit);
 
     iEvent.put(outputHcalIsoTrackColl, "HcalIsolatedTrackCollection");
     iEvent.put(outputTColl,    labelGenTrack_.encode());
     iEvent.put(outputVColl,    labelRecVtx_.encode());
     iEvent.put(outputEBColl,   labelEB_.encode());
     iEvent.put(outputEEColl,   labelEE_.encode());
     iEvent.put(outputHBHEColl, labelHBHE_.encode());
       }
     }
   }
 }

reco::HcalIsolatedTrackCandidateCollection * AlCaIsoTracksProducer::select	(	edm::Handle< edm::TriggerResults > &	triggerResults,
		const std::vector< std::string > &	triggerNames_,
		edm::Handle< reco::TrackCollection > &	trkCollection,
		math::XYZPoint &	leadPV,
		edm::Handle< EcalRecHitCollection > &	barrelRecHitsHandle,
		edm::Handle< EcalRecHitCollection > &	endcapRecHitsHandle,
		edm::Handle< HBHERecHitCollection > &	hbhe
	)

private

Definition at line 370 of file AlCaIsoTracksProducer.cc.

References a_charIsoR, a_mipR, bField, spr::chargeIsolationCone(), spr::eCone_ecal(), eEcalMax_, eIsolation_, reco::TrackBase::eta(), etaL1, spr::find(), geo, spr::goodTrack(), runregparse::hlt, i, HcalDetId::ieta(), HcalDetId::iphi(), LogDebug, convertSQLiteXML::ok, reco::TrackBase::p(), phi, reco::TrackBase::phi(), phiL1, spr::propagateCALO(), reco::TrackBase::pt(), ptL1, pTrackMin_, reco::TrackBase::px(), reco::TrackBase::py(), reco::TrackBase::pz(), selectionParameters, reco::HcalIsolatedTrackCandidate::SetEnergyEcal(), reco::HcalIsolatedTrackCandidate::SetEtaPhiEcal(), reco::HcalIsolatedTrackCandidate::SetEtaPhiHcal(), reco::HcalIsolatedTrackCandidate::setL1(), reco::HcalIsolatedTrackCandidate::SetMaxP(), reco::HcalIsolatedTrackCandidate::setTrack(), theTrackQuality, HLT_50ns_5e33_v1_cff::trackCollection, trigKount, trigNames, and trigPass.

Referenced by Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::find(), Vispa.Views.LineDecayView.LineDecayContainer::mousePressEvent(), Vispa.Gui.PortConnection.PointToPointConnection::mousePressEvent(), Vispa.Gui.VispaWidget.VispaWidget::mousePressEvent(), produce(), and Vispa.Views.AbstractView.AbstractView::restoreSelection().

                                                          {
   
   reco::HcalIsolatedTrackCandidateCollection* trackCollection=new reco::HcalIsolatedTrackCandidateCollection;
   bool ok(false);
 
   // Find a good HLT trigger
   for (unsigned int iHLT=0; iHLT<triggerResults->size(); iHLT++) {
     int hlt    = triggerResults->accept(iHLT);
     for (unsigned int i=0; i<trigNames.size(); ++i) {
       if (triggerNames_[iHLT].find(trigNames[i].c_str())!=std::string::npos) {
     trigKount[i]++;
     if (hlt > 0) {
       ok = true;
       trigPass[i]++;
     }
     LogDebug("HcalIsoTrack") <<"This is the trigger we are looking for "
                  << triggerNames_[iHLT] << " Flag " << hlt 
                  << ":" << ok;
       }
     }
   }
   if (ok) {
     //Propagate tracks to calorimeter surface)
     std::vector<spr::propagatedTrackDirection> trkCaloDirections;
     spr::propagateCALO(trkCollection, geo, bField, theTrackQuality,
                trkCaloDirections, false);
 
     std::vector<spr::propagatedTrackDirection>::const_iterator trkDetItr;
     unsigned int nTracks(0), nselTracks(0);
     for (trkDetItr = trkCaloDirections.begin(),nTracks=0; 
      trkDetItr != trkCaloDirections.end(); trkDetItr++,nTracks++) {
       const reco::Track* pTrack = &(*(trkDetItr->trkItr));
       math::XYZTLorentzVector v4(pTrack->px(), pTrack->py(), 
                  pTrack->pz(), pTrack->p());
       LogDebug("HcalIsoTrack") << "This track : " << nTracks 
                    << " (pt|eta|phi|p) :" << pTrack->pt() 
                    << "|" << pTrack->eta() << "|" 
                    << pTrack->phi() << "|" << pTrack->p();
         
       //Selection of good track
       bool qltyFlag  = spr::goodTrack(pTrack,leadPV,selectionParameters,false);
       LogDebug("HcalIsoTrack") << "qltyFlag|okECAL|okHCAL : " << qltyFlag
                    << "|" << trkDetItr->okECAL << "|" 
                    << trkDetItr->okHCAL;
       if (qltyFlag && trkDetItr->okECAL && trkDetItr->okHCAL) {
     double t_p        = pTrack->p();
     nselTracks++;
     int    nRH_eMipDR(0), nNearTRKs(0);
     double eMipDR = spr::eCone_ecal(geo, barrelRecHitsHandle, 
                     endcapRecHitsHandle,
                     trkDetItr->pointHCAL,
                     trkDetItr->pointECAL, a_mipR,
                     trkDetItr->directionECAL, 
                     nRH_eMipDR);
     double hmaxNearP = spr::chargeIsolationCone(nTracks,
                             trkCaloDirections,
                             a_charIsoR, 
                             nNearTRKs, false);
     LogDebug("HcalIsoTrack") << "This track : " << nTracks 
                  << " (pt|eta|phi|p) :"  << pTrack->pt() 
                  << "|" << pTrack->eta() << "|" 
                  << pTrack->phi() << "|" << t_p
                  << "e_MIP " << eMipDR 
                  << " Chg Isolation " << hmaxNearP;
     if (t_p>pTrackMin_ && eMipDR<eEcalMax_ && hmaxNearP<eIsolation_) {
       reco::HcalIsolatedTrackCandidate newCandidate(v4);
       newCandidate.SetMaxP(hmaxNearP);
       newCandidate.SetEnergyEcal(eMipDR);
       newCandidate.setL1(ptL1,etaL1,phiL1);
       newCandidate.SetEtaPhiEcal((trkDetItr->pointECAL).eta(),
                      (trkDetItr->pointECAL).phi());
       HcalDetId detId = HcalDetId(trkDetItr->detIdHCAL);
       newCandidate.SetEtaPhiHcal((trkDetItr->pointHCAL).eta(),
                      (trkDetItr->pointHCAL).phi(),
                      detId.ieta(), detId.iphi());
       int indx(0);
       reco::TrackCollection::const_iterator trkItr1;
       for (trkItr1=trkCollection->begin(); trkItr1 != trkCollection->end(); ++trkItr1,++indx) {
         const reco::Track* pTrack1 = &(*trkItr1);
         if (pTrack1 == pTrack) {
           reco::TrackRef tRef = reco::TrackRef(trkCollection,indx);
           newCandidate.setTrack(tRef);
           break;
         }
       }
       trackCollection->push_back(newCandidate);
     }
       }
     }
   }
   return trackCollection;
 }

void AlCaIsoTracksProducer::setPtEtaPhi	(	std::vector< edm::Ref< l1extra::L1JetParticleCollection > > &	objref,
		double &	ptL1,
		double &	etaL1,
		double &	phiL1
	)

private

Definition at line 468 of file AlCaIsoTracksProducer.cc.

References AlCaHLTBitMon_ParallelJobs::p, EnergyCorrector::pt, and ptL1.

Referenced by produce().

                                                                                                                                                {
 
   for (unsigned int p=0; p<objref.size(); p++) {
     if (objref[p]->pt()>ptL1) {
       ptL1  = objref[p]->pt(); 
       phiL1 = objref[p]->phi();
       etaL1 = objref[p]->eta();
     }
   }
 }