Inheritance diagram for HcalIsoTrkAnalyzer:

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

	~HcalIsoTrkAnalyzer ()

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

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 ()

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)

virtual void	beginJob ()

virtual void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

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

double	dR (math::XYZTLorentzVector &, math::XYZTLorentzVector &)

virtual void	endJob ()

virtual void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

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

Private Attributes
double	a_charIsoR

double	a_coneR

double	a_mipR

bool	changed

double	eEcalMax_

double	eIsolation_

edm::Service< TFileService >	fs

HLTConfigProvider	hltConfig_

std::vector< std::string >	HLTNames

std::string	l1Filter

std::string	l2Filter

std::string	l3Filter

edm::InputTag	labelEB_

edm::InputTag	labelEE_

edm::InputTag	labelGenTrack_

edm::InputTag	labelHBHE_

edm::InputTag	labelRecVtx_

int	nAll

int	nGood

int	nRun

std::vector< double >	pbin

std::string	processName

double	pTrackMin_

spr::trackSelectionParameters	selectionParameters

std::vector< unsigned int > *	t_DetIds

double	t_eHcal

double	t_eMipDR

int	t_Event

double	t_EventWeight

std::vector< double > *	t_HitEnergies

double	t_hmaxNearP

int	t_ieta

double	t_l1eta

double	t_l1phi

double	t_l1pt

double	t_l3eta

double	t_l3phi

double	t_l3pt

double	t_mindR1

double	t_mindR2

double	t_p

bool	t_qltyFlag

bool	t_qltyMissFlag

bool	t_qltyPVFlag

int	t_Run

bool	t_selectTk

std::vector< bool > *	t_trgbits

std::string	theTrackQuality

edm::InputTag	theTriggerResultsLabel

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 < reco::VertexCollection >	tok_recVtx_

edm::EDGetTokenT < trigger::TriggerEvent >	tok_trigEvt

edm::EDGetTokenT < edm::TriggerResults >	tok_trigRes

TTree *	tree

std::vector< std::string >	trgnames

edm::InputTag	triggerEvent_

std::vector< int >	trigKount

std::vector< std::string >	trigNames

std::vector< int >	trigPass

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

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 63 of file HcalIsoTrkAnalyzer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 117 of file HcalIsoTrkAnalyzer.cc.

                                                                      : 
   changed(false), nRun(0), nAll(0), nGood(0) {
   //now do whatever initialization is needed
   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");
   reco::TrackBase::TrackQuality trackQuality_=reco::TrackBase::qualityByName(theTrackQuality);
   selectionParameters.minPt           = iConfig.getParameter<double>("MinTrackPt");
   selectionParameters.minQuality      = trackQuality_;
   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 + 28.9;
   a_mipR                              = iConfig.getParameter<double>("ConeRadiusMIP");
   pTrackMin_                          = iConfig.getParameter<double>("MinimumTrackP");
   eEcalMax_                           = iConfig.getParameter<double>("MaximumEcalEnergy");
   eIsolation_                         = iConfig.getParameter<double>("IsolationEnergy");
   triggerEvent_                       = iConfig.getParameter<edm::InputTag>("TriggerEventLabel");
   theTriggerResultsLabel              = iConfig.getParameter<edm::InputTag>("TriggerResultLabel");
   labelGenTrack_                      = iConfig.getParameter<edm::InputTag>("TrackLabel");
   labelRecVtx_                        = iConfig.getParameter<edm::InputTag>("VertexLabel");
   labelEB_                            = iConfig.getParameter<edm::InputTag>("EBRecHitLabel");
   labelEE_                            = iConfig.getParameter<edm::InputTag>("EERecHitLabel");
   labelHBHE_                          = iConfig.getParameter<edm::InputTag>("HBHERecHitLabel");
 
   // define tokens for access
   tok_trigEvt   = consumes<trigger::TriggerEvent>(triggerEvent_);
   tok_trigRes   = consumes<edm::TriggerResults>(theTriggerResultsLabel);
   tok_genTrack_ = consumes<reco::TrackCollection>(labelGenTrack_);
   tok_recVtx_   = consumes<reco::VertexCollection>(labelRecVtx_);
   tok_bs_       = consumes<reco::BeamSpot>(iConfig.getParameter<edm::InputTag>("BeamSpotLabel"));
  
   tok_EB_       = consumes<EcalRecHitCollection>(labelEB_);
   tok_EE_       = consumes<EcalRecHitCollection>(labelEE_);
   tok_hbhe_     = consumes<HBHERecHitCollection>(labelHBHE_);
 
   std::vector<int> dummy(trigNames.size(),0);
   trigKount = trigPass = dummy;
   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;
   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];
 }

HcalIsoTrkAnalyzer::~HcalIsoTrkAnalyzer ( )

Definition at line 186 of file HcalIsoTrkAnalyzer.cc.

                                         {
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
 }

Member Function Documentation

void HcalIsoTrkAnalyzer::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 191 of file HcalIsoTrkAnalyzer.cc.

References a_charIsoR, a_coneR, a_mipR, ecalTB2006H4_GenSimDigiReco_cfg::bField, edm::EventBase::bunchCrossing(), changed, spr::chargeIsolationCone(), dR(), spr::eCone_ecal(), spr::eCone_hcal(), eEcalMax_, eIsolation_, trigger::TriggerObject::energy(), trigger::TriggerObject::eta(), reco::TrackBase::eta(), edm::EventID::event(), trigger::TriggerEvent::filterKeys(), trigger::TriggerEvent::filterTag(), spr::find(), edm::EventSetup::get(), edm::Event::getByToken(), trigger::TriggerEvent::getObjects(), spr::goodTrack(), runregparse::hlt, hltConfig_, i, trigger::TriggerObject::id(), edm::EventBase::id(), HcalDetId::ieta(), edm::HandleBase::isValid(), relval_steps::k, l1Filter, l2Filter, l3Filter, edm::InputTag::label(), diffTwoXMLs::label, labelEB_, labelEE_, labelGenTrack_, labelHBHE_, LogDebug, edm::EventBase::luminosityBlock(), trigger::TriggerObject::mass(), spr::trackSelectionParameters::maxDxyPV, spr::trackSelectionParameters::maxDzPV, spr::trackSelectionParameters::maxInMiss, spr::trackSelectionParameters::maxOutMiss, HLTConfigProvider::moduleLabels(), edmConvertToStreamModule::modules, gen::n, nAll, nGood, convertSQLiteXML::ok, reco::TrackBase::p(), trigger::TriggerObject::phi(), reco::TrackBase::phi(), edm::Handle< T >::product(), edm::ESHandle< class >::product(), spr::propagateCALO(), trigger::TriggerObject::pt(), reco::TrackBase::pt(), pTrackMin_, trigger::TriggerObject::px(), reco::TrackBase::px(), trigger::TriggerObject::py(), reco::TrackBase::py(), trigger::TriggerObject::pz(), reco::TrackBase::pz(), edm::EventID::run(), selectionParameters, HLTConfigProvider::size(), findQualityFiles::size, trigger::TriggerEvent::sizeFilters(), AlCaHLTBitMon_QueryRunRegistry::string, t_DetIds, t_eHcal, t_eMipDR, t_Event, t_EventWeight, t_HitEnergies, t_hmaxNearP, t_ieta, t_l1eta, t_l1phi, t_l1pt, t_l3eta, t_l3phi, t_l3pt, t_mindR1, t_mindR2, t_p, t_qltyFlag, t_qltyMissFlag, t_qltyPVFlag, t_Run, t_selectTk, t_trgbits, theTrackQuality, tok_bs_, tok_EB_, tok_EE_, tok_genTrack_, tok_hbhe_, tok_recVtx_, tok_trigEvt, tok_trigRes, tree, trgnames, triggerEvent_, HLTConfigProvider::triggerIndex(), edm::TriggerNames::triggerNames(), edm::Event::triggerNames(), triggerResults, trigKount, trigNames, trigPass, x, detailsBasic3DVector::y, and detailsBasic3DVector::z.

                                                                                     {
 
   t_Run   = iEvent.id().run();
   t_Event = iEvent.id().event();
   nAll++;
   LogDebug("HcalIsoTrack") << "Run " << t_Run << " Event " << t_Event 
                << " Luminosity " << iEvent.luminosityBlock() 
                << " Bunch " << iEvent.bunchCrossing();
 
   //Get magnetic field and ECAL channel status
   edm::ESHandle<MagneticField> bFieldH;
   iSetup.get<IdealMagneticFieldRecord>().get(bFieldH);
   const MagneticField *bField = bFieldH.product();
 
   // get handles to calogeometry and calotopology
   edm::ESHandle<CaloGeometry> pG;
   iSetup.get<CaloGeometryRecord>().get(pG);
   const CaloGeometry* geo = pG.product();
   
   bool okC(true);
   //Get track collection
   edm::Handle<reco::TrackCollection> trkCollection;
   iEvent.getByToken(tok_genTrack_, trkCollection);
   reco::TrackCollection::const_iterator trkItr;
   if (!trkCollection.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelGenTrack_;
     okC = false;
   }
  
   //event weight for FLAT sample
   t_EventWeight = 1.0;
 
   //Define the best vertex and the beamspot
   edm::Handle<reco::VertexCollection> recVtxs;
   iEvent.getByToken(tok_recVtx_, recVtxs);  
   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;
   if (beamSpotH.isValid()) LogDebug("HcalIsoTrack") << "Beam Spot " 
                             << beamSpotH->position();
   
   // RecHits
   edm::Handle<EcalRecHitCollection> barrelRecHitsHandle;
   iEvent.getByToken(tok_EB_, barrelRecHitsHandle);
   if (!barrelRecHitsHandle.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelEB_;
     okC = false;
   }
   edm::Handle<EcalRecHitCollection> endcapRecHitsHandle;
   iEvent.getByToken(tok_EE_, endcapRecHitsHandle);
   if (!endcapRecHitsHandle.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelEE_;
     okC = false;
   }
   edm::Handle<HBHERecHitCollection> hbhe;
   iEvent.getByToken(tok_hbhe_, hbhe);
   if (!hbhe.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelHBHE_;
     okC = false;
   }
       
   //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;
 
   //Trigger
   trgnames.clear();
   t_trgbits->clear();
   for (unsigned int i=0; i<trigNames.size(); ++i) t_trgbits->push_back(false);
 
   trigger::TriggerEvent triggerEvent;
   edm::Handle<trigger::TriggerEvent> triggerEventHandle;
   iEvent.getByToken(tok_trigEvt, triggerEventHandle);
   if (!triggerEventHandle.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Error! Can't get the product "
                     << triggerEvent_.label() ;
   } else if (okC) {
     triggerEvent = *(triggerEventHandle.product());
     
     const trigger::TriggerObjectCollection& TOC(triggerEvent.getObjects());
     edm::Handle<edm::TriggerResults> triggerResults;
     iEvent.getByToken(tok_trigRes, triggerResults);
     if (triggerResults.isValid()) {
       std::vector<std::string> modules;
       const edm::TriggerNames & triggerNames = iEvent.triggerNames(*triggerResults);
       const std::vector<std::string> & triggerNames_ = triggerNames.triggerNames();
       for (unsigned int iHLT=0; iHLT<triggerResults->size(); iHLT++) {
     bool ok(false);
     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) {
         t_trgbits->at(i) = (hlt>0);
         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) {
       unsigned int triggerindx = hltConfig_.triggerIndex(triggerNames_[iHLT]);
       const std::vector<std::string>& moduleLabels(hltConfig_.moduleLabels(triggerindx));
       std::vector<math::XYZTLorentzVector> vec[3];
       //loop over all trigger filters in event (i.e. filters passed)
       for (unsigned int ifilter=0; ifilter<triggerEvent.sizeFilters();
            ++ifilter) {  
         std::vector<int> Keys;
         std::string label = triggerEvent.filterTag(ifilter).label();
         //loop over keys to objects passing this filter
         for (unsigned int imodule=0; imodule<moduleLabels.size(); 
          imodule++) {
           if (label.find(moduleLabels[imodule]) != std::string::npos) {
         LogDebug("HcalIsoTrack") << "FilterName " << label;
         for (unsigned int ifiltrKey=0; ifiltrKey<triggerEvent.filterKeys(ifilter).size(); ++ifiltrKey) {
           Keys.push_back(triggerEvent.filterKeys(ifilter)[ifiltrKey]);
           const trigger::TriggerObject& TO(TOC[Keys[ifiltrKey]]);
           math::XYZTLorentzVector v4(TO.px(), TO.py(), TO.pz(), TO.energy());
           if (label.find(l2Filter) != std::string::npos) {
             vec[1].push_back(v4);
           } else if (label.find(l3Filter) != std::string::npos) {
             vec[2].push_back(v4);
           } else if (label.find(l1Filter) != std::string::npos ||
                  l1Filter == "") {
             vec[0].push_back(v4);
           }
           LogDebug("HcalIsoTrack") << "key " << ifiltrKey << " : pt " 
                        << TO.pt() << " eta " << TO.eta() 
                        << " phi " << TO.phi() << " mass " 
                        << TO.mass() << " Id " << TO.id();
         }
         LogDebug("HcalIsoTrack") << "sizes " << vec[0].size() << ":" 
                      << vec[1].size() <<":" <<vec[2].size();
           }
         }
       }
       double dr;
       math::XYZTLorentzVector mindRvec1;
       double mindR1(999);
       for (int lvl=1; lvl<3; lvl++) {
         for (unsigned int i=0; i<vec[lvl].size(); i++) {
           dr   = dR(vec[0][0],vec[lvl][i]);
           LogDebug("HcalIsoTrack") << "lvl " <<lvl << " i " << i << " dR " 
                        << dr;
           if (dr<mindR1) {
         mindR1    = dr;
         mindRvec1 = vec[lvl][i];
           }
         }
       }
       
       t_l1pt  = vec[0][0].pt();
       t_l1eta = vec[0][0].eta();
       t_l1phi = vec[0][0].phi();
           if (vec[2].size()>0) {
         t_l3pt  = vec[2][0].pt();
         t_l3eta = vec[2][0].eta();
         t_l3phi = vec[2][0].phi();
       }
       //Loop over tracks
       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();
         math::XYZTLorentzVector mindRvec2;
         t_mindR2 = 999;
 
         for (unsigned int k=0; k<vec[2].size(); ++k) {
           dr   = dR(vec[2][k],v4); //changed 1 to 2
           if (dr<t_mindR2) {
         t_mindR2  = dr;
         mindRvec2 = vec[2][k];
           }
         }
         t_mindR1 = dR(vec[0][0],v4);
         LogDebug("HcalIsoTrack") << "Closest L3 object at mindr :" 
                      << t_mindR2 << " is " << mindRvec2 
                      << " and from L1 " << t_mindR1;
         
         //Selection of good track
         t_selectTk = spr::goodTrack(pTrack,leadPV,selectionParameters,false);
         spr::trackSelectionParameters oneCutParameters = selectionParameters;
         oneCutParameters.maxDxyPV  = 10;
         oneCutParameters.maxDzPV   = 100;
         oneCutParameters.maxInMiss = 2;
         oneCutParameters.maxOutMiss= 2;
         bool qltyFlag  = spr::goodTrack(pTrack,leadPV,oneCutParameters,false);
         oneCutParameters           = selectionParameters;
         oneCutParameters.maxDxyPV  = 10;
         oneCutParameters.maxDzPV   = 100;
         t_qltyMissFlag = spr::goodTrack(pTrack,leadPV,oneCutParameters,false);
         oneCutParameters           = selectionParameters;
         oneCutParameters.maxInMiss = 2;
         oneCutParameters.maxOutMiss= 2;
         t_qltyPVFlag   = spr::goodTrack(pTrack,leadPV,oneCutParameters,false);
         t_ieta = 0;
         if (trkDetItr->okHCAL) {
           HcalDetId detId = (HcalDetId)(trkDetItr->detIdHCAL);
           t_ieta = detId.ieta();
         }
         LogDebug("HcalIsoTrack") << "qltyFlag|okECAL|okHCAL : " << qltyFlag
                      << "|" << trkDetItr->okECAL << "|" 
                      << trkDetItr->okHCAL;
         t_qltyFlag = (qltyFlag && trkDetItr->okECAL && trkDetItr->okHCAL);
         t_p        = pTrack->p();
         if (t_qltyFlag) {
               nselTracks++;
           int nRH_eMipDR(0), nNearTRKs(0), nRecHits(-999);
           t_eMipDR = spr::eCone_ecal(geo, barrelRecHitsHandle, 
                      endcapRecHitsHandle,
                      trkDetItr->pointHCAL,
                      trkDetItr->pointECAL,
                      a_mipR, trkDetItr->directionECAL, 
                      nRH_eMipDR);
           t_DetIds->clear(); t_HitEnergies->clear();
           std::vector<DetId> ids;
           t_eHcal = spr::eCone_hcal(geo, hbhe, trkDetItr->pointHCAL, 
                     trkDetItr->pointECAL, a_coneR, 
                     trkDetItr->directionHCAL,nRecHits, 
                     ids, *t_HitEnergies);
           for (unsigned int k=0; k<ids.size(); ++k) {
         t_DetIds->push_back(ids[k].rawId());
           }
           t_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;
           LogDebug("HcalIsoTrack") << "e_MIP " << t_eMipDR 
                        << " Chg Isolation " << t_hmaxNearP 
                        << " eHcal" << t_eHcal << " ieta " 
                        << t_ieta << " Quality " 
                        << t_qltyMissFlag << ":" 
                        << t_qltyPVFlag << ":" << t_selectTk;
           for (unsigned int lll=0;lll<t_DetIds->size();lll++) {
         LogDebug("HcalIsoTrack") << "det id is = " <<t_DetIds->at(lll)
                      << "   hit enery is  = "  
                      << t_HitEnergies->at(lll) ;
           }
           if (t_p>pTrackMin_ && t_eMipDR<eEcalMax_ && 
           t_hmaxNearP<eIsolation_) {
         tree->Fill();
         nGood++;
         for (unsigned int k=0; k<t_trgbits->size(); k++) 
           LogDebug("HcalIsoTrack") << "trigger bit is  = " 
                        << t_trgbits->at(k);
           }
         }
       }
     }
       }
       // check if trigger names in (new) config                       
       if (changed) {
     changed = false;
     LogDebug("HcalIsoTrack") << "New trigger menu found !!!";
     const unsigned int n(hltConfig_.size());
     for (unsigned itrig=0; itrig<triggerNames_.size(); itrig++) {
       unsigned int triggerindx = hltConfig_.triggerIndex(triggerNames_[itrig]);
       if (triggerindx >= n)
         LogDebug("HcalIsoTrack") << triggerNames_[itrig] << " " 
                      << triggerindx << " does not exist in "
                      << "the current menu";
       else
         LogDebug("HcalIsoTrack") << triggerNames_[itrig] << " " 
                      << triggerindx << " exists";
     }
       }
     }
   }
 }

void HcalIsoTrkAnalyzer::beginJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 484 of file HcalIsoTrkAnalyzer.cc.

References fs, TFileService::make(), t_DetIds, t_eHcal, t_eMipDR, t_Event, t_EventWeight, t_HitEnergies, t_hmaxNearP, t_ieta, t_l1eta, t_l1phi, t_l1pt, t_l3eta, t_l3phi, t_l3pt, t_mindR1, t_mindR2, t_p, t_qltyFlag, t_qltyMissFlag, t_qltyPVFlag, t_Run, t_selectTk, t_trgbits, and tree.

                                   {
 
   tree = fs->make<TTree>("CalibTree", "CalibTree");
      
   tree->Branch("t_Run",         &t_Run,         "t_Run/I");
   tree->Branch("t_Event",       &t_Event,       "t_Event/I");
   tree->Branch("t_ieta",        &t_ieta,        "t_ieta/I");
   tree->Branch("t_EventWeight", &t_EventWeight, "t_EventWeight/D");
   tree->Branch("t_l1pt",        &t_l1pt,        "t_l1pt/D");
   tree->Branch("t_l1eta",       &t_l1eta,       "t_l1eta/D");
   tree->Branch("t_l1phi",       &t_l1phi,       "t_l1phi/D"); 
   tree->Branch("t_l3pt",        &t_l3pt,        "t_l3pt/D");
   tree->Branch("t_l3eta",       &t_l3eta,       "t_l3eta/D"); 
   tree->Branch("t_l3phi",       &t_l3phi,       "t_l3phi/D");
   tree->Branch("t_p",           &t_p,           "t_p/D");
   tree->Branch("t_mindR1",      &t_mindR1,      "t_mindR1/D");
   tree->Branch("t_mindR2",      &t_mindR2,      "t_mindR2/D");
   tree->Branch("t_eMipDR",      &t_eMipDR,      "t_eMipDR/D");
   tree->Branch("t_eHcal",       &t_eHcal,       "t_eHcal/D");
   tree->Branch("t_hmaxNearP",   &t_hmaxNearP,   "t_hmaxNearP/D");
   tree->Branch("t_selectTk",    &t_selectTk,    "t_selectTk/O");
   tree->Branch("t_qltyFlag",    &t_qltyFlag,    "t_qltyFlag/O");
   tree->Branch("t_qltyMissFlag",&t_qltyMissFlag,"t_qltyMissFlag/O");
   tree->Branch("t_qltyPVFlag",  &t_qltyPVFlag,  "t_qltyPVFlag/O)");
 
   t_DetIds      = new std::vector<unsigned int>();
   t_HitEnergies = new std::vector<double>();
   t_trgbits     = new std::vector<bool>();
   tree->Branch("t_DetIds",      "std::vector<unsigned int>", &t_DetIds);
   tree->Branch("t_HitEnergies", "std::vector<double>",       &t_HitEnergies);
   tree->Branch("t_trgbits","std::vector<bool>", &t_trgbits); 
 }

void HcalIsoTrkAnalyzer::beginLuminosityBlock	(	edm::LuminosityBlock const &	,
		edm::EventSetup const &
	)

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 540 of file HcalIsoTrkAnalyzer.cc.

540 {}

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

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 528 of file HcalIsoTrkAnalyzer.cc.

References changed, hltConfig_, HLTConfigProvider::init(), nRun, processName, and edm::RunBase::run().

                                                                                  {
   edm::LogInfo("HcalIsoTrack") << "Run[" << nRun << "] " << iRun.run() 
                    << " hltconfig.init " << hltConfig_.init(iRun,iSetup,processName,changed);
 }

double HcalIsoTrkAnalyzer::dR	(	math::XYZTLorentzVector &	vec1,
		math::XYZTLorentzVector &	vec2
	)

private

Definition at line 552 of file HcalIsoTrkAnalyzer.cc.

References reco::deltaR().

Referenced by analyze().

                                                                                       {
   return reco::deltaR(vec1.eta(),vec1.phi(),vec2.eta(),vec2.phi());
 }

void HcalIsoTrkAnalyzer::endJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 518 of file HcalIsoTrkAnalyzer.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 HcalIsoTrkAnalyzer::endLuminosityBlock	(	edm::LuminosityBlock const &	,
		edm::EventSetup const &
	)

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 542 of file HcalIsoTrkAnalyzer.cc.

542 {}

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

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 534 of file HcalIsoTrkAnalyzer.cc.

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

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

void HcalIsoTrkAnalyzer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 544 of file HcalIsoTrkAnalyzer.cc.

References edm::ConfigurationDescriptions::addDefault(), and edm::ParameterSetDescription::setUnknown().

                                                                                     {
   //The following says we do not know what parameters are allowed so do no validation
   // Please change this to state exactly what you do use, even if it is no parameters
   edm::ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }

Member Data Documentation

double HcalIsoTrkAnalyzer::a_charIsoR

private

Definition at line 90 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

double HcalIsoTrkAnalyzer::a_coneR

private

Definition at line 90 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

double HcalIsoTrkAnalyzer::a_mipR

private

Definition at line 90 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

bool HcalIsoTrkAnalyzer::changed

private

Definition at line 82 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginRun().

double HcalIsoTrkAnalyzer::eEcalMax_

private

Definition at line 91 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

double HcalIsoTrkAnalyzer::eIsolation_

private

Definition at line 91 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::Service<TFileService> HcalIsoTrkAnalyzer::fs

private

Definition at line 83 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob().

HLTConfigProvider HcalIsoTrkAnalyzer::hltConfig_

private

Definition at line 84 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginRun().

std::vector<std::string> HcalIsoTrkAnalyzer::HLTNames

private

Definition at line 85 of file HcalIsoTrkAnalyzer.cc.

std::string HcalIsoTrkAnalyzer::l1Filter

private

Definition at line 89 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

std::string HcalIsoTrkAnalyzer::l2Filter

private

Definition at line 89 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

std::string HcalIsoTrkAnalyzer::l3Filter

private

Definition at line 89 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::InputTag HcalIsoTrkAnalyzer::labelEB_

private

Definition at line 94 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::InputTag HcalIsoTrkAnalyzer::labelEE_

private

Definition at line 94 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::InputTag HcalIsoTrkAnalyzer::labelGenTrack_

private

Definition at line 94 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::InputTag HcalIsoTrkAnalyzer::labelHBHE_

private

Definition at line 93 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::InputTag HcalIsoTrkAnalyzer::labelRecVtx_

private

Definition at line 94 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

int HcalIsoTrkAnalyzer::nAll

private

Definition at line 92 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and endJob().

int HcalIsoTrkAnalyzer::nGood

private

Definition at line 92 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and endJob().

int HcalIsoTrkAnalyzer::nRun

private

Definition at line 92 of file HcalIsoTrkAnalyzer.cc.

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

std::vector<double> HcalIsoTrkAnalyzer::pbin

private

Definition at line 112 of file HcalIsoTrkAnalyzer.cc.

std::string HcalIsoTrkAnalyzer::processName

private

Definition at line 88 of file HcalIsoTrkAnalyzer.cc.

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

double HcalIsoTrkAnalyzer::pTrackMin_

private

Definition at line 91 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

spr::trackSelectionParameters HcalIsoTrkAnalyzer::selectionParameters

private

Definition at line 87 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

std::vector<unsigned int>* HcalIsoTrkAnalyzer::t_DetIds

private

Definition at line 111 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_eHcal

private

Definition at line 109 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_eMipDR

private

Definition at line 109 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

int HcalIsoTrkAnalyzer::t_Event

private

Definition at line 106 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_EventWeight

private

Definition at line 107 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

std::vector<double>* HcalIsoTrkAnalyzer::t_HitEnergies

private

Definition at line 112 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_hmaxNearP

private

Definition at line 109 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

int HcalIsoTrkAnalyzer::t_ieta

private

Definition at line 106 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_l1eta

private

Definition at line 107 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_l1phi

private

Definition at line 107 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_l1pt

private

Definition at line 107 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_l3eta

private

Definition at line 108 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_l3phi

private

Definition at line 108 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_l3pt

private

Definition at line 108 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_mindR1

private

Definition at line 108 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_mindR2

private

Definition at line 109 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_p

private

Definition at line 108 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

bool HcalIsoTrkAnalyzer::t_qltyFlag

private

Definition at line 110 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

bool HcalIsoTrkAnalyzer::t_qltyMissFlag

private

Definition at line 110 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

bool HcalIsoTrkAnalyzer::t_qltyPVFlag

private

Definition at line 110 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

int HcalIsoTrkAnalyzer::t_Run

private

Definition at line 106 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

bool HcalIsoTrkAnalyzer::t_selectTk

private

Definition at line 110 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

std::vector<bool>* HcalIsoTrkAnalyzer::t_trgbits

private

Definition at line 113 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

std::string HcalIsoTrkAnalyzer::theTrackQuality

private

Definition at line 88 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::InputTag HcalIsoTrkAnalyzer::theTriggerResultsLabel

private

Definition at line 93 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

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

private

Definition at line 100 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::EDGetTokenT<EcalRecHitCollection> HcalIsoTrkAnalyzer::tok_EB_

private

Definition at line 101 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::EDGetTokenT<EcalRecHitCollection> HcalIsoTrkAnalyzer::tok_EE_

private

Definition at line 102 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

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

private

Definition at line 98 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::EDGetTokenT<HBHERecHitCollection> HcalIsoTrkAnalyzer::tok_hbhe_

private

Definition at line 103 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

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

private

Definition at line 99 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::EDGetTokenT<trigger::TriggerEvent> HcalIsoTrkAnalyzer::tok_trigEvt

private

Definition at line 95 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::EDGetTokenT<edm::TriggerResults> HcalIsoTrkAnalyzer::tok_trigRes

private

Definition at line 96 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

TTree* HcalIsoTrkAnalyzer::tree

private

Definition at line 105 of file HcalIsoTrkAnalyzer.cc.

Referenced by eventstfile.Events::__iter__(), svgfig.SVG::__str__(), analyze(), beginJob(), MetTreeProducer.MetTreeProducer::declareVariables(), core.AutoFillTreeProducer.AutoFillTreeProducer::declareVariables(), core.AutoFillTreeProducer.AutoFillTreeProducer::fillTree(), MetTreeProducer.MetTreeProducer::process(), python.cmstools.EventTree::SetAlias(), and eventstfile.Events::to().

std::vector<std::string> HcalIsoTrkAnalyzer::trgnames

private

Definition at line 114 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze().

edm::InputTag HcalIsoTrkAnalyzer::triggerEvent_

private

Definition at line 93 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

std::vector<int> HcalIsoTrkAnalyzer::trigKount

private

Definition at line 86 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), endJob(), and HcalIsoTrkAnalyzer().

std::vector<std::string> HcalIsoTrkAnalyzer::trigNames

private

Definition at line 85 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), endJob(), and HcalIsoTrkAnalyzer().

std::vector<int> HcalIsoTrkAnalyzer::trigPass

private

Definition at line 86 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), endJob(), and HcalIsoTrkAnalyzer().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation