AlCaIsoTracksProducer Class Reference

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_
std::string	labelIsoTk_
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 125 of file AlCaIsoTracksProducer.cc.

References a_charIsoR, a_coneR, a_mipR, eEcalMax_, eIsolation_, edm::ParameterSet::getParameter(), edm::InputTag::instance(), relval_steps::k, l1Filter, l2Filter, l3Filter, edm::InputTag::label(), labelBS_, labelEB_, labelEE_, labelGenTrack_, labelHBHE_, labelHltGT_, labelIsoTk_, labelRecVtx_, labelTriggerEvent_, labelTriggerResults_, spr::trackSelectionParameters::maxChi2, spr::trackSelectionParameters::maxDpOverP, spr::trackSelectionParameters::maxDxyPV, spr::trackSelectionParameters::maxDzPV, spr::trackSelectionParameters::maxInMiss, spr::trackSelectionParameters::maxOutMiss, spr::trackSelectionParameters::minLayerCrossed, spr::trackSelectionParameters::minOuterHit, spr::trackSelectionParameters::minPt, spr::trackSelectionParameters::minQuality, processName, pTrackMin_, reco::TrackBase::qualityByName(), selectionParameters, AlCaHLTBitMon_QueryRunRegistry::string, theTrackQuality, tok_bs_, tok_EB_, tok_EE_, tok_genTrack_, tok_hbhe_, tok_hltGT_, tok_recVtx_, tok_trigEvt_, tok_trigRes_, trigKount, trigNames, and trigPass.

                                                                           :
  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");
  labelIsoTk_                         = iConfig.getParameter<std::string>("IsoTrackLabel");

  // 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>(labelIsoTk_);
  produces<reco::VertexCollection>(labelRecVtx_.label());
  produces<EcalRecHitCollection>(labelEB_.instance());
  produces<EcalRecHitCollection>(labelEE_.instance());
  produces<HBHERecHitCollection>(labelHBHE_.label());

  edm::LogInfo("HcalIsoTrack") << " Expected to produce the collections:\n"
                   << "reco::HcalIsolatedTrackCandidateCollection "
                   << " with label HcalIsolatedTrackCollection\n"
                   << "reco::VertexCollection with label " << labelRecVtx_.label() << "\n"
                   << "EcalRecHitCollection with label EcalRecHitsEB\n"
                   << "EcalRecHitCollection with label EcalRecHitsEE\n"
                   << "HBHERecHitCollection with label " << labelHBHE_.label();
}

AlCaIsoTracksProducer::~AlCaIsoTracksProducer

(

)

Definition at line 217 of file AlCaIsoTracksProducer.cc.

{ }

Member Function Documentation

void AlCaIsoTracksProducer::beginJob ( void  )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 350 of file AlCaIsoTracksProducer.cc.

{ }

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

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 361 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 352 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 374 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 220 of file AlCaIsoTracksProducer.cc.

References edm::EventBase::bunchCrossing(), etaL1, edm::EventID::event(), edm::Event::getByToken(), edm::EventBase::id(), edm::InputTag::instance(), edm::HandleBase::isValid(), edm::InputTag::label(), labelEB_, labelEE_, labelGenTrack_, labelHBHE_, labelIsoTk_, labelRecVtx_, labelTriggerEvent_, labelTriggerResults_, 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++;
#ifdef DebugLog
  edm::LogInfo("HcalIsoTrack") << "Run " << iEvent.id().run() << " Event " 
                   << iEvent.id().event() << " Luminosity " 
                   << iEvent.luminosityBlock() << " Bunch " 
                   << iEvent.bunchCrossing();
#endif
  //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();
  }
#ifdef DebugLog
  edm::LogInfo("HcalIsoTrack") << "Primary Vertex " << leadPV;
#endif

  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);
#ifdef DebugLog
      edm::LogInfo("HcalIsoTrack") << "AlCaIsoTracksProducer::select returns "
                   << isotk->size() << " isolated tracks";
#endif
      std::auto_ptr<reco::HcalIsolatedTrackCandidateCollection> outputHcalIsoTrackColl(isotk);
      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);
    
      if (isotk->size() > 0) {
    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, labelIsoTk_);
      iEvent.put(outputVColl,            labelRecVtx_.label());
      iEvent.put(outputEBColl,           labelEB_.instance());
      iEvent.put(outputEEColl,           labelEE_.instance());
      iEvent.put(outputHBHEColl,         labelHBHE_.label());
    }
  }
}

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 380 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(), 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_25ns14e33_v3_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]++;
    }
    edm::LogInfo("HcalIsoTrack") << "The trigger we are looking for "
                     << triggerNames_[iHLT] << " Flag " 
                     << hlt << ":" << 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());
#ifdef DebugLog
    edm::LogInfo("HcalIsoTrack") << "This track : " << nTracks 
                 << " (pt|eta|phi|p) :" << pTrack->pt() 
                 << "|" << pTrack->eta() << "|" 
                 << pTrack->phi() << "|" << pTrack->p();
#endif      
    //Selection of good track
    bool qltyFlag  = spr::goodTrack(pTrack,leadPV,selectionParameters,false);
#ifdef DebugLog
    edm::LogInfo("HcalIsoTrack") << "qltyFlag|okECAL|okHCAL : " << qltyFlag
                 << "|" << trkDetItr->okECAL << "|"  
                 << trkDetItr->okHCAL;
#endif
    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);
#ifdef DebugLog
      edm::LogInfo("HcalIsoTrack") << "This track : " << nTracks 
                   << " (pt|eta|phi|p) :"  << pTrack->pt() 
                   << "|" << pTrack->eta() << "|" 
                   << pTrack->phi() << "|" << t_p 
                   << " e_MIP " << eMipDR 
                   << " Chg Isolation " << hmaxNearP;
#endif
      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 482 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();
    }
  }
}

Member Data Documentation

double AlCaIsoTracksProducer::a_charIsoR

private

Definition at line 102 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), and select().

double AlCaIsoTracksProducer::a_coneR

private

Definition at line 102 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer().

double AlCaIsoTracksProducer::a_mipR

private

Definition at line 102 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), and select().

const MagneticField* AlCaIsoTracksProducer::bField

private

Definition at line 109 of file AlCaIsoTracksProducer.cc.

Referenced by beginRun(), and select().

double AlCaIsoTracksProducer::eEcalMax_

private

Definition at line 103 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), and select().

double AlCaIsoTracksProducer::eIsolation_

private

Definition at line 103 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), and select().

double AlCaIsoTracksProducer::etaL1

private

Definition at line 111 of file AlCaIsoTracksProducer.cc.

Referenced by produce(), and select().

const CaloGeometry* AlCaIsoTracksProducer::geo

private

Definition at line 110 of file AlCaIsoTracksProducer.cc.

Referenced by beginRun(), and select().

HLTConfigProvider AlCaIsoTracksProducer::hltConfig_

private

Definition at line 96 of file AlCaIsoTracksProducer.cc.

Referenced by beginRun().

std::vector<std::string> AlCaIsoTracksProducer::HLTNames

private

Definition at line 97 of file AlCaIsoTracksProducer.cc.

std::string AlCaIsoTracksProducer::l1Filter

private

Definition at line 101 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer().

std::string AlCaIsoTracksProducer::l2Filter

private

Definition at line 101 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer().

std::string AlCaIsoTracksProducer::l3Filter

private

Definition at line 101 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer().

edm::InputTag AlCaIsoTracksProducer::labelBS_

private

Definition at line 107 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer().

edm::InputTag AlCaIsoTracksProducer::labelEB_

private

Definition at line 107 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), and produce().

edm::InputTag AlCaIsoTracksProducer::labelEE_

private

Definition at line 107 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), and produce().

edm::InputTag AlCaIsoTracksProducer::labelGenTrack_

private

Definition at line 106 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), and produce().

edm::InputTag AlCaIsoTracksProducer::labelHBHE_

private

Definition at line 107 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), and produce().

edm::InputTag AlCaIsoTracksProducer::labelHltGT_

private

Definition at line 106 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer().

std::string AlCaIsoTracksProducer::labelIsoTk_

private

Definition at line 108 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), and produce().

edm::InputTag AlCaIsoTracksProducer::labelRecVtx_

private

Definition at line 106 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), and produce().

edm::InputTag AlCaIsoTracksProducer::labelTriggerEvent_

private

Definition at line 105 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), and produce().

edm::InputTag AlCaIsoTracksProducer::labelTriggerResults_

private

Definition at line 105 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), and produce().

int AlCaIsoTracksProducer::nAll

private

Definition at line 104 of file AlCaIsoTracksProducer.cc.

Referenced by endJob(), and produce().

int AlCaIsoTracksProducer::nGood

private

Definition at line 104 of file AlCaIsoTracksProducer.cc.

Referenced by endJob().

int AlCaIsoTracksProducer::nRun

private

Definition at line 104 of file AlCaIsoTracksProducer.cc.

Referenced by beginRun(), endJob(), and endRun().

double AlCaIsoTracksProducer::phiL1

private

Definition at line 111 of file AlCaIsoTracksProducer.cc.

Referenced by produce(), and select().

std::string AlCaIsoTracksProducer::processName

private

Definition at line 100 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), beginRun(), core.TriggerBitAnalyzer.TriggerBitAnalyzer::declareHandles(), and core.TriggerBitFilter.TriggerBitFilter::declareHandles().

double AlCaIsoTracksProducer::ptL1

private

Definition at line 111 of file AlCaIsoTracksProducer.cc.

Referenced by produce(), select(), and setPtEtaPhi().

double AlCaIsoTracksProducer::pTrackMin_

private

Definition at line 103 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), and select().

spr::trackSelectionParameters AlCaIsoTracksProducer::selectionParameters

private

Definition at line 99 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), and select().

std::string AlCaIsoTracksProducer::theTrackQuality

private

Definition at line 100 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), and select().

edm::EDGetTokenT<reco::BeamSpot> AlCaIsoTracksProducer::tok_bs_

private

Definition at line 118 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), and produce().

edm::EDGetTokenT<EcalRecHitCollection> AlCaIsoTracksProducer::tok_EB_

private

Definition at line 119 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), and produce().

edm::EDGetTokenT<EcalRecHitCollection> AlCaIsoTracksProducer::tok_EE_

private

Definition at line 120 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), and produce().

edm::EDGetTokenT<reco::TrackCollection> AlCaIsoTracksProducer::tok_genTrack_

private

Definition at line 116 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), and produce().

edm::EDGetTokenT<HBHERecHitCollection> AlCaIsoTracksProducer::tok_hbhe_

private

Definition at line 121 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), and produce().

edm::EDGetTokenT<trigger::TriggerFilterObjectWithRefs> AlCaIsoTracksProducer::tok_hltGT_

private

Definition at line 113 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), and produce().

edm::EDGetTokenT<reco::VertexCollection> AlCaIsoTracksProducer::tok_recVtx_

private

Definition at line 117 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), and produce().

edm::EDGetTokenT<trigger::TriggerEvent> AlCaIsoTracksProducer::tok_trigEvt_

private

Definition at line 114 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), and produce().

edm::EDGetTokenT<edm::TriggerResults> AlCaIsoTracksProducer::tok_trigRes_

private

Definition at line 115 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), and produce().

std::vector<int> AlCaIsoTracksProducer::trigKount

private

Definition at line 98 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), endJob(), and select().

std::vector<std::string> AlCaIsoTracksProducer::trigNames

private

Definition at line 97 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), endJob(), and select().

std::vector<int> AlCaIsoTracksProducer::trigPass

private

Definition at line 98 of file AlCaIsoTracksProducer.cc.

Referenced by AlCaIsoTracksProducer(), endJob(), and select().