#include <IsoTrig/IsoTrig/src/IsoTrig.cc>

Inheritance diagram for IsoTrig:

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

	~IsoTrig ()

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

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

void	chgIsolation (double &etaTriggered, double &phiTriggered, edm::Handle< reco::TrackCollection > &trkCollection, const edm::Event &theEvent)

void	clearChgIsolnTreeVectors ()

void	clearMipCutTreeVectors ()

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

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

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

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

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

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

std::pair< double, double >	etaPhiTrigger ()

void	fillCuts (int, double, double, double, math::XYZTLorentzVector &, int, bool)

void	fillDifferences (int, math::XYZTLorentzVector &, math::XYZTLorentzVector &, bool)

void	fillEnergy (int, int, double, double, math::XYZTLorentzVector &)

void	fillHist (int, math::XYZTLorentzVector &)

double	getDistInCM (double eta1, double phi1, double eta2, double phi2)

std::pair< double, double >	GetEtaPhiAtEcal (double etaIP, double phiIP, double pT, int charge, double vtxZ)

void	pushChgIsolnTreeVecs (math::XYZTLorentzVector &Pixcand, math::XYZTLorentzVector &Trkcand, std::vector< double > &PixMaxP, double &TrkMaxP, bool &selTk)

void	pushMipCutTreeVecs (math::XYZTLorentzVector &NFcand, math::XYZTLorentzVector &Trkcand, double &EmipNFcand, double &EmipTrkcand, double &mindR, double &mindP1, std::vector< bool > &Flags, double hCone)

void	studyIsolation (edm::Handle< reco::TrackCollection > &, std::vector< reco::TrackCollection::const_iterator > &)

void	studyMipCut (edm::Handle< reco::TrackCollection > &trkCollection, edm::Handle< reco::IsolatedPixelTrackCandidateCollection > &L2cands)

void	studyTiming (const edm::Event &theEvent)

void	studyTrigger (edm::Handle< reco::TrackCollection > &, std::vector< reco::TrackCollection::const_iterator > &)

void	StudyTrkEbyP (edm::Handle< reco::TrackCollection > &trkCollection)

Private Attributes
double	a_charIsoR

double	a_coneR

double	a_mipR

double	a_neutIsoR

double	a_neutR1

double	a_neutR2

edm::Handle< EcalRecHitCollection >	barrelRecHitsHandle

edm::Handle< reco::BeamSpot >	beamSpotH

edm::ESHandle< MagneticField >	bFieldH

double	bfVal

bool	changed

TTree *	ChgIsolnTree

double	cutCharge

double	cutMip

double	cutNeutral

bool	doChgIsolTree

bool	doL2L3

bool	doMipCutTree

bool	doStudyIsol

bool	doTiming

bool	doTrkResTree

double	dr_L1

edm::Handle< EcalRecHitCollection >	endcapRecHitsHandle

edm::Service< TFileService >	fs

TH1I *	g_Accepts

TH1I *	g_Pre

TH1I *	g_PreHLT

TH1I *	g_PreL1

TH1D *	h_dEta [9]

TH1D *	h_dEtaL1 [2]

TH1D *	h_dinvPt [9]

TH1D *	h_dP [9]

TH1D *	h_dPhi [9]

TH1D *	h_dPhiL1 [2]

TH1D *	h_dPt [9]

TH1D *	h_dRL1 [2]

TH1D *	h_eCalo [5][6][48]

TH1D *	h_eHcal [5][6][48]

TH1D *	h_eMaxNearP [2]

TH1D *	h_eMip [2]

TH1D *	h_eNeutIso [2]

TH1D *	h_EnIn

TH1D *	h_EnOut

TH1D *	h_eta [20]

TH1D *	h_etaCalibTracks [5][2][2]

TH1D *	h_etaMipTracks [5][2][2]

TH1I *	h_Filters

TH1I *	h_HLT

TH1D *	h_L1ObjEnergy

TH1D *	h_mindR [9]

TH2D *	h_MipEnMatch

TH2D *	h_MipEnNoMatch

TH1I *	h_nHLT

TH1I *	h_nL3Objs

TH1D *	h_p [20]

TH1D *	h_phi [20]

TH1I *	h_Pre

TH1I *	h_PreHLT

TH1D *	h_PreHLTwt

TH1I *	h_PreL1

TH1D *	h_PreL1wt

TH1D *	h_pt [20]

edm::Handle< HBHERecHitCollection >	hbhe

HLTConfigProvider	hltConfig_

edm::InputTag	L1candTag_

edm::InputTag	L2candTag_

math::XYZPoint	leadPV

int	maxRunNo

double	minPTrackValue_

int	minRunNo

TTree *	MipCutTree

edm::ESHandle< CaloGeometry >	pG

edm::InputTag	PixcandTag_

std::vector< double >	pixelIsolationConeSizeAtEC_

std::vector< reco::TrackRef >	pixelTrackRefsHB

std::vector< reco::TrackRef >	pixelTrackRefsHE

std::vector< edm::InputTag >	pixelTracksSources_

double	pLimits [6]

double	prelimCone_

std::string	processName

double	rEB_

edm::Handle < reco::VertexCollection >	recVtxs

spr::trackSelectionParameters	selectionParameters

std::vector< double > *	t_ECone

std::vector< double > *	t_NFcandEmip

std::vector< double > *	t_NFcandEta

std::vector< double > *	t_NFcandP

std::vector< double > *	t_NFcandPhi

std::vector< double > *	t_NFcandPt

std::vector< double > *	t_NFTrkcandEmip

std::vector< double > *	t_NFTrkcandEta

std::vector< double > *	t_NFTrkcandP

std::vector< double > *	t_NFTrkcandPhi

std::vector< double > *	t_NFTrkcandPt

std::vector< bool > *	t_NFTrkChgIsoFlag

std::vector< double > *	t_NFTrkMinDP1

std::vector< double > *	t_NFTrkMinDR

std::vector< bool > *	t_NFTrkMipFlag

std::vector< bool > *	t_NFTrkMissFlag

std::vector< bool > *	t_NFTrkNeuIsoFlag

std::vector< bool > *	t_NFTrkPropFlag

std::vector< bool > *	t_NFTrkPVFlag

std::vector< bool > *	t_NFTrkqltyFlag

std::vector< bool > *	t_NFTrkselTkFlag

std::vector< int > *	t_nPixCand

std::vector< int > *	t_nPixSeed

std::vector< double > *	t_PixcandEta

std::vector< std::vector < double > > *	t_PixcandMaxP

std::vector< double > *	t_PixcandP

std::vector< double > *	t_PixcandPhi

std::vector< double > *	t_PixcandPt

std::vector< double > *	t_PixTrkcandEta

std::vector< double > *	t_PixTrkcandMaxP

std::vector< double > *	t_PixTrkcandP

std::vector< double > *	t_PixTrkcandPhi

std::vector< double > *	t_PixTrkcandPt

std::vector< bool > *	t_PixTrkcandselTk

std::vector< double > *	t_timeL2Prod

std::vector< double > *	t_TrkhCone

std::vector< bool > *	t_TrkMissFlag

std::vector< bool > *	t_TrkNuIsolFlag

std::vector< double > *	t_TrkP

std::vector< bool > *	t_TrkPVFlag

std::vector< bool > *	t_TrkqltyFlag

std::vector< bool > *	t_TrkselTkFlag

double	tauUnbiasCone_

std::string	theTrackQuality

TTree *	TimingTree

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_hlt_

edm::EDGetTokenT < trigger::TriggerFilterObjectWithRefs >	tok_l1cand_

edm::EDGetTokenT < reco::IsolatedPixelTrackCandidateCollection >	tok_l2cand_

edm::EDGetTokenT< LumiDetails >	tok_lumi

edm::EDGetTokenT < reco::IsolatedPixelTrackCandidateCollection >	tok_pixtk_

std::vector< edm::EDGetTokenT < reco::TrackCollection > >	tok_pixtks_

edm::EDGetTokenT < reco::VertexCollection >	tok_recVtx_

edm::EDGetTokenT < trigger::TriggerEvent >	tok_trigEvt

edm::EDGetTokenT < edm::TriggerResults >	tok_trigRes

edm::EDGetTokenT < reco::VertexCollection >	tok_verthb_

edm::EDGetTokenT < reco::VertexCollection >	tok_verthe_

std::map< unsigned int, unsigned int >	TrigList

std::vector< std::string >	trigNames

std::map< unsigned int, const std::pair< int, int > >	TrigPreList

TTree *	TrkResTree

std::vector < math::XYZTLorentzVector >	vec [3]

int	verbosity

double	vtxCutIsol_

double	vtxCutSeed_

double	zEE_

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

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 65 of file IsoTrig.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 29 of file IsoTrig.cc.

                                                : changed(false) {
   //now do whatever initialization is neededOA
   trigNames                           = iConfig.getUntrackedParameter<std::vector<std::string> >("Triggers");
   PixcandTag_                          = iConfig.getParameter<edm::InputTag> ("PixcandTag");
   L1candTag_                          = iConfig.getParameter<edm::InputTag> ("L1candTag");
   L2candTag_                          = iConfig.getParameter<edm::InputTag> ("L2candTag");
   doL2L3                              = iConfig.getUntrackedParameter<bool>("DoL2L3",true);
   doTiming                            = iConfig.getUntrackedParameter<bool>("DoTimingTree",true);
   doMipCutTree                        = iConfig.getUntrackedParameter<bool>("DoMipCutTree",true);
   doTrkResTree                        = iConfig.getUntrackedParameter<bool>("DoTrkResTree",true);
   doChgIsolTree                        = iConfig.getUntrackedParameter<bool>("DoChgIsolTree",true);
   doStudyIsol                         = iConfig.getUntrackedParameter<bool>("DoStudyIsol",true);
   verbosity                           = iConfig.getUntrackedParameter<int>("Verbosity",0);
   theTrackQuality                     = iConfig.getUntrackedParameter<std::string>("TrackQuality","highPurity");
   processName                         = iConfig.getUntrackedParameter<std::string>("ProcessName","HLT");
   reco::TrackBase::TrackQuality trackQuality_=reco::TrackBase::qualityByName(theTrackQuality);
   selectionParameters.minPt           = iConfig.getUntrackedParameter<double>("MinTrackPt", 10.0);
   selectionParameters.minQuality      = trackQuality_;
   selectionParameters.maxDxyPV        = iConfig.getUntrackedParameter<double>("MaxDxyPV", 0.2);
   selectionParameters.maxDzPV         = iConfig.getUntrackedParameter<double>("MaxDzPV",  5.0);
   selectionParameters.maxChi2         = iConfig.getUntrackedParameter<double>("MaxChi2",  5.0);
   selectionParameters.maxDpOverP      = iConfig.getUntrackedParameter<double>("MaxDpOverP",  0.1);
   selectionParameters.minOuterHit     = iConfig.getUntrackedParameter<int>("MinOuterHit", 4);
   selectionParameters.minLayerCrossed = iConfig.getUntrackedParameter<int>("MinLayerCrossed", 8);
   selectionParameters.maxInMiss       = iConfig.getUntrackedParameter<int>("MaxInMiss", 0);
   selectionParameters.maxOutMiss      = iConfig.getUntrackedParameter<int>("MaxOutMiss", 0);
   dr_L1                               = iConfig.getUntrackedParameter<double>("IsolationL1",1.0);
   a_coneR                             = iConfig.getUntrackedParameter<double>("ConeRadius",34.98);
   a_charIsoR                          = a_coneR + 28.9;
   a_neutIsoR                          = a_charIsoR*0.726;
   a_mipR                              = iConfig.getUntrackedParameter<double>("ConeRadiusMIP",14.0);
   a_neutR1                            = iConfig.getUntrackedParameter<double>("ConeRadiusNeut1",21.0);
   a_neutR2                            = iConfig.getUntrackedParameter<double>("ConeRadiusNeut2",29.0);
   cutMip                              = iConfig.getUntrackedParameter<double>("MIPCut", 1.0);
   cutCharge                           = iConfig.getUntrackedParameter<double>("ChargeIsolation",  2.0);
   cutNeutral                          = iConfig.getUntrackedParameter<double>("NeutralIsolation",  2.0);
   minRunNo                            = iConfig.getUntrackedParameter<int>("minRun");
   maxRunNo                            = iConfig.getUntrackedParameter<int>("maxRun");
   pixelTracksSources_                 = iConfig.getUntrackedParameter<std::vector<edm::InputTag> >("PixelTracksSources");
   pixelIsolationConeSizeAtEC_         = iConfig.getUntrackedParameter<std::vector<double> >("PixelIsolationConeSizeAtEC");
   minPTrackValue_                     = iConfig.getUntrackedParameter<double>("MinPTrackValue");
   vtxCutSeed_                         = iConfig.getUntrackedParameter<double>("VertexCutSeed");
   vtxCutIsol_                         = iConfig.getUntrackedParameter<double>("VertexCutIsol");
   tauUnbiasCone_                      = iConfig.getUntrackedParameter<double>("TauUnbiasCone");
   prelimCone_                         = iConfig.getUntrackedParameter<double>("PrelimCone");
   // define tokens for access
   tok_lumi      = consumes<LumiDetails, edm::InLumi>(edm::InputTag("lumiProducer"));
   edm::InputTag triggerEvent_ ("hltTriggerSummaryAOD","",processName); 
   tok_trigEvt   = consumes<trigger::TriggerEvent>(triggerEvent_);
   edm::InputTag theTriggerResultsLabel ("TriggerResults","",processName); 
   tok_trigRes   = consumes<edm::TriggerResults>(theTriggerResultsLabel);
   tok_genTrack_ = consumes<reco::TrackCollection>(edm::InputTag("generalTracks"));
   tok_recVtx_   = consumes<reco::VertexCollection>(edm::InputTag("offlinePrimaryVertices"));
   tok_bs_       = consumes<reco::BeamSpot>(edm::InputTag("offlineBeamSpot"));
   tok_EB_     = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit","EcalRecHitsEB"));
   tok_EE_     = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit","EcalRecHitsEE"));
   tok_hbhe_   = consumes<HBHERecHitCollection>(edm::InputTag("hbhereco"));
   tok_pixtk_    = consumes<reco::IsolatedPixelTrackCandidateCollection>(PixcandTag_);
   tok_l1cand_   = consumes<trigger::TriggerFilterObjectWithRefs>(L1candTag_);
   tok_l2cand_   = consumes<reco::IsolatedPixelTrackCandidateCollection>(L2candTag_);
   if(doTiming){
     tok_verthb_ = consumes<reco::VertexCollection>(edm::InputTag("hltHITPixelVerticesHB"));
     tok_verthe_ = consumes<reco::VertexCollection>(edm::InputTag("hltHITPixelVerticesHB"));
     tok_hlt_ = consumes<trigger::TriggerFilterObjectWithRefs>(edm::InputTag("hltL1sL1SingleJet68"));
   }
   if(doChgIsolTree) {
     for (unsigned int k=0; k<pixelTracksSources_.size(); ++k) {
       //      edm::InputTag  pix (pixelTracksSources_[k],"",processName);
       //      tok_pixtks_.push_back(consumes<reco::TrackCollection>(pix));
       tok_pixtks_.push_back(consumes<reco::TrackCollection>(pixelTracksSources_[k]));
     }
   }
   if (verbosity>=0) {
     std::cout <<"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_neutIsoR "      << a_neutIsoR          
           <<"\t a_mipR "          << a_mipR 
           <<"\t a_neutR "         << a_neutR1 << ":" << a_neutR2
               <<"\t cuts (MIP "       << cutMip << " : Charge " << cutCharge
           <<" : Neutral "         << cutNeutral << ")"
           << std::endl;
     std::cout <<"Charge Isolation parameters:"
           <<"\t minPTrackValue " << minPTrackValue_
           <<"\t vtxCutSeed "     << vtxCutSeed_
           <<"\t vtxCutIsol "     << vtxCutIsol_
           <<"\t tauUnbiasCone "  << tauUnbiasCone_
           <<"\t prelimCone "     << prelimCone_
           <<"\t pixelIsolationConeSizeAtEC";
     for (unsigned int k=0; k<pixelIsolationConeSizeAtEC_.size(); ++k)
       std::cout << " " << pixelIsolationConeSizeAtEC_[k];
     std::cout << std::endl;
   }
   double pl[] = {20,30,40,60,80,120};
   for (int i=0; i<6; ++i) pLimits[i] = pl[i];
   rEB_ = 123.8;
   zEE_ = 317.0;
 }

IsoTrig::~IsoTrig ( )

Definition at line 139 of file IsoTrig.cc.

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

Member Function Documentation

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

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 153 of file IsoTrig.cc.

                                                                          {
   if (verbosity%10 > 0) std::cout << "Event starts====================================" << std::endl;
   int RunNo = iEvent.id().run();
   int EvtNo = iEvent.id().event();
   int Lumi  = iEvent.luminosityBlock();
   int Bunch = iEvent.bunchCrossing();
 
   iSetup.get<IdealMagneticFieldRecord>().get(bFieldH);
   iSetup.get<CaloGeometryRecord>().get(pG);
   const MagneticField *bField = bFieldH.product();   
   GlobalVector BField=bField->inTesla(GlobalPoint(0,0,0));
   bfVal = BField.mag();
 
   trigger::TriggerEvent triggerEvent;
   edm::Handle<trigger::TriggerEvent> triggerEventHandle;
   iEvent.getByToken(tok_trigEvt, triggerEventHandle);
   if (!triggerEventHandle.isValid()) {
     if (verbosity%10 > 0)   std::cout << "Error! Can't get the product hltTriggerSummaryAOD"<< std::endl;
   } else {
     triggerEvent = *(triggerEventHandle.product());
   }
   const trigger::TriggerObjectCollection& TOC(triggerEvent.getObjects());
   edm::Handle<edm::TriggerResults> triggerResults;
   iEvent.getByToken(tok_trigRes, triggerResults);  
 
   edm::Handle<reco::TrackCollection> trkCollection;
   iEvent.getByToken(tok_genTrack_, trkCollection);
 
   iEvent.getByToken(tok_EB_, barrelRecHitsHandle);
   iEvent.getByToken(tok_EE_, endcapRecHitsHandle);
 
   iEvent.getByToken(tok_hbhe_, hbhe);
 
   iEvent.getByToken(tok_recVtx_, recVtxs);  
   iEvent.getByToken(tok_bs_, beamSpotH);
   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();
   }
   if ((verbosity/100)%10>0) {
     std::cout << "Primary Vertex " << leadPV;
     if (beamSpotH.isValid()) std::cout << " Beam Spot " 
                        << beamSpotH->position();
     std::cout << std::endl;
   }
   pixelTrackRefsHE.clear(); pixelTrackRefsHB.clear();
   for (unsigned int iPix=0; iPix<pixelTracksSources_.size(); iPix++) {
     edm::Handle<reco::TrackCollection> iPixCol;
     iEvent.getByToken(tok_pixtks_[iPix],iPixCol); 
     if(iPixCol.isValid()){
       for (reco::TrackCollection::const_iterator pit=iPixCol->begin(); pit!=iPixCol->end(); pit++) {
     if(iPix==0) 
       pixelTrackRefsHB.push_back(reco::TrackRef(iPixCol,pit-iPixCol->begin()));
     pixelTrackRefsHE.push_back(reco::TrackRef(iPixCol,pit-iPixCol->begin()));
       }
     }
   }
   if(doTiming)  
     studyTiming(iEvent);
   for (unsigned int ifilter=0; ifilter<triggerEvent.sizeFilters(); 
        ++ifilter) {  
      std::string FilterNames[7] = {"hltL1sL1SingleJet68", "hltIsolPixelTrackL2FilterHE", "ecalIsolPixelTrackFilterHE", "hltIsolPixelTrackL3FilterHE",
                   "hltIsolPixelTrackL2FilterHB", "ecalIsolPixelTrackFilterHB", "hltIsolPixelTrackL3FilterHB"};
     std::string label = triggerEvent.filterTag(ifilter).label();
     for(int i=0; i<7; i++) {
       if(label==FilterNames[i]) h_Filters->Fill(i);
     }
   }
   edm::InputTag lumiProducer("LumiProducer", "", "RECO");
   edm::Handle<LumiDetails> Lumid;
   iEvent.getLuminosityBlock().getByToken(tok_lumi, Lumid);
   float mybxlumi=-1;
   if (Lumid.isValid()) 
     mybxlumi=Lumid->lumiValue(LumiDetails::kOCC1,iEvent.bunchCrossing())*6.37;
   if (verbosity%10 > 0)
     std::cout << "RunNo " << RunNo << " EvtNo " << EvtNo << " Lumi " << Lumi 
           << " Bunch " << Bunch << " mybxlumi " << mybxlumi << std::endl;
 
   if (!triggerResults.isValid()) {
     std::cout << "Error! Can't get the product triggerResults" << std::endl;
     //      boost::shared_ptr<cms::Exception> const & error = triggerResults.whyFailed();
     //      edm::LogWarning(error->category()) << error->what();
   } else {
     std::vector<std::string> modules;
     h_nHLT->Fill(triggerResults->size());
     const edm::TriggerNames & triggerNames = iEvent.triggerNames(*triggerResults);
     
     const std::vector<std::string> & triggerNames_ = triggerNames.triggerNames();
     //      std::cout << "number of HLTs " << triggerNames_.size() << std::endl;
     int hlt(-1), preL1(-1), preHLT(-1), prescale(-1);
     for (unsigned int i=0; i<triggerResults->size(); i++) {
       unsigned int triggerindx = hltConfig_.triggerIndex(triggerNames_[i]);
       const std::vector<std::string>& moduleLabels(hltConfig_.moduleLabels(triggerindx));
       
       for (unsigned int in=0; in<trigNames.size(); ++in) {
     //    if (triggerNames_[i].find(trigNames[in].c_str())!=std::string::npos || triggerNames_[i]==" ") {
     if (triggerNames_[i].find(trigNames[in].c_str())!=std::string::npos) {
       if (verbosity%10 > 0) std::cout << "trigger that i want " << triggerNames_[i] << " accept " << triggerResults->accept(i) << std::endl;
       hlt    = triggerResults->accept(i);
       h_HLT      -> Fill(hlt);
       //        if (hlt>0 || triggerNames_[i]==" ") {
       if (hlt>0) {
         edm::Handle<reco::IsolatedPixelTrackCandidateCollection> Pixcands;
         iEvent.getByToken(tok_pixtk_,Pixcands); 
         edm::Handle<trigger::TriggerFilterObjectWithRefs> L1cands;
         iEvent.getByToken(tok_l1cand_, L1cands); 
         
         const std::pair<int,int> prescales(hltConfig_.prescaleValues(iEvent,iSetup,triggerNames_[i]));
         preL1  = prescales.first;
         preHLT = prescales.second;
         prescale = preL1*preHLT;
         if (verbosity%10 > 0)
           std::cout << triggerNames_[i] << " accept " << hlt << " preL1 " 
             << preL1 << " preHLT " << preHLT << std::endl;
         
         for (int iv=0; iv<3; ++iv) vec[iv].clear();
         if (TrigList.find(RunNo) != TrigList.end() ) {
           TrigList[RunNo] += 1;
           } else {
           TrigList.insert(std::pair<unsigned int, unsigned int>(RunNo,1));
           TrigPreList.insert(std::pair<unsigned int, std::pair<int, int>>(RunNo,prescales));
         }
         //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) {
           if (verbosity%10 > 0) std::cout << "FILTERNAME " << label << std::endl;
           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 {
               vec[0].push_back(v4);
               h_L1ObjEnergy->Fill(TO.energy());
             }
             if (verbosity%10 > 0)
               std::cout << "key " << ifiltrKey << " : pt " << TO.pt() << " eta " << TO.eta() << " phi " << TO.phi() << " mass " << TO.mass() << " Id " << TO.id() << std::endl;
           }
         }
           }
         }
         std::vector<reco::TrackCollection::const_iterator> goodTks;
         if (doL2L3) {
           h_nL3Objs  -> Fill(vec[2].size());
           studyTrigger(trkCollection, goodTks);
         } else {
           if (trkCollection.isValid()) {
         reco::TrackCollection::const_iterator trkItr;
         for (trkItr=trkCollection->begin(); 
              trkItr!=trkCollection->end(); trkItr++) 
             goodTks.push_back(trkItr);
           }
         }
         // Now study isolation etc
         if(doStudyIsol) studyIsolation(trkCollection, goodTks);
         
         if(doTrkResTree) StudyTrkEbyP(trkCollection);
         
         std::pair<double,double> etaphi = etaPhiTrigger();
         edm::Handle<reco::IsolatedPixelTrackCandidateCollection> L2cands;
         iEvent.getByToken(tok_l2cand_,L2cands); 
         if (!L2cands.isValid()) {
           if (verbosity%10 > 0) std::cout << " trigCand is not valid " << std::endl;
         } else {
           if(doMipCutTree) studyMipCut(trkCollection, L2cands);
         }
         if (pixelTracksSources_.size()>0)
           if(doChgIsolTree && pixelTrackRefsHE.size()>0) chgIsolation(etaphi.first, etaphi.second, trkCollection, iEvent);
       }
       break;
     }
       }
     }
     h_PreL1    -> Fill(preL1);
     h_PreHLT   -> Fill(preHLT);
     h_Pre      -> Fill(prescale);
     h_PreL1wt  -> Fill(preL1, mybxlumi);
     h_PreHLTwt -> Fill(preHLT, mybxlumi);
     
     // check if trigger names in (new) config                       
     //      std::cout << "changed " <<changed << std::endl;
     if (changed) {
       changed = false;
       if ((verbosity/10)%10 > 1) {
     std::cout<<"New trigger menu found !!!" << std::endl;
     const unsigned int n(hltConfig_.size());
     for (unsigned itrig=0; itrig<triggerNames_.size(); itrig++) {
       unsigned int triggerindx = hltConfig_.triggerIndex(triggerNames_[itrig]);
       std::cout << triggerNames_[itrig] << " " << triggerindx << " ";
       if (triggerindx >= n)
         std::cout << "does not exist in the current menu" << std::endl;
       else
         std::cout << "exists" << std::endl;
     }
       }
     }
   }  
 }

void IsoTrig::beginJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 448 of file IsoTrig.cc.

References ChgIsolnTree, Reconstruction_Data_cff::cuts, doChgIsolTree, doL2L3, doMipCutTree, doTiming, doTrkResTree, fs, h_dEta, h_dEtaL1, h_dinvPt, h_dP, h_dPhi, h_dPhiL1, h_dPt, h_dRL1, h_eCalo, h_eHcal, h_eMaxNearP, h_eMip, h_eNeutIso, h_EnIn, h_EnOut, h_eta, h_etaCalibTracks, h_etaMipTracks, h_Filters, h_HLT, h_L1ObjEnergy, h_mindR, h_MipEnMatch, h_MipEnNoMatch, h_nHLT, h_nL3Objs, h_p, h_phi, h_Pre, h_PreHLT, h_PreHLTwt, h_PreL1, h_PreL1wt, h_pt, i, TFileService::make(), MipCutTree, filesave_online::pairs, pLimits, AlCaHLTBitMon_QueryRunRegistry::string, t_ECone, t_NFcandEmip, t_NFcandEta, t_NFcandP, t_NFcandPhi, t_NFcandPt, t_NFTrkcandEmip, t_NFTrkcandEta, t_NFTrkcandP, t_NFTrkcandPhi, t_NFTrkcandPt, t_NFTrkChgIsoFlag, t_NFTrkMinDP1, t_NFTrkMinDR, t_NFTrkMipFlag, t_NFTrkMissFlag, t_NFTrkNeuIsoFlag, t_NFTrkPropFlag, t_NFTrkPVFlag, t_NFTrkqltyFlag, t_NFTrkselTkFlag, t_nPixCand, t_nPixSeed, t_PixcandEta, t_PixcandMaxP, t_PixcandP, t_PixcandPhi, t_PixcandPt, t_PixTrkcandEta, t_PixTrkcandMaxP, t_PixTrkcandP, t_PixTrkcandPhi, t_PixTrkcandPt, t_PixTrkcandselTk, t_timeL2Prod, t_TrkhCone, t_TrkMissFlag, t_TrkNuIsolFlag, t_TrkP, t_TrkPVFlag, t_TrkqltyFlag, t_TrkselTkFlag, TimingTree, and TrkResTree.

                        {
   char hname[100], htit[100];
   std::string levels[20] = {"L1", "L2", "L3", 
                 "Reco", "RecoMatch", "RecoNoMatch", 
                 "L2Match", "L2NoMatch", "L3Match", "L3NoMatch", 
                 "HLTTrk", "HLTGoodTrk", "HLTIsoTrk", "HLTMip", "HLTSelect",
                 "nonHLTTrk", "nonHLTGoodTrk", "nonHLTIsoTrk", "nonHLTMip", "nonHLTSelect"};
   if(doTiming) {
     TimingTree = fs->make<TTree>("TimingTree", "TimingTree");
     t_timeL2Prod = new std::vector<double>();
     t_nPixCand   = new std::vector<int>();
     t_nPixSeed   = new std::vector<int>();
 
     TimingTree->Branch("t_timeL2Prod", "std::vector<double>", &t_timeL2Prod);
     TimingTree->Branch("t_nPixCand", "std::vector<int>", &t_nPixCand);
     TimingTree->Branch("t_nPixSeed", "std::vector<int>", &t_nPixSeed);
   }
   if(doTrkResTree) {
     TrkResTree = fs->make<TTree>("TrkRestree", "TrkRestree");
     t_TrkhCone     = new std::vector<double>();
     t_TrkP         = new std::vector<double>();
     t_TrkselTkFlag = new std::vector<bool>();
     t_TrkqltyFlag  = new std::vector<bool>();
     t_TrkMissFlag  = new std::vector<bool>();
     t_TrkPVFlag    = new std::vector<bool>();
     t_TrkNuIsolFlag= new std::vector<bool>();
     
     TrkResTree->Branch("t_TrkhCone",     "std::vector<double>", &t_TrkhCone);
     TrkResTree->Branch("t_TrkP",         "std::vector<double>", &t_TrkP);
     TrkResTree->Branch("t_TrkselTkFlag", "std::vector<bool>",   &t_TrkselTkFlag);
     TrkResTree->Branch("t_TrkqltyFlag",  "std::vector<bool>",   &t_TrkqltyFlag);
     TrkResTree->Branch("t_TrkMissFlag",  "std::vector<bool>",   &t_TrkMissFlag);
     TrkResTree->Branch("t_TrkPVFlag",    "std::vector<bool>",   &t_TrkPVFlag);
     TrkResTree->Branch("t_TrkNuIsolFlag","std::vector<bool>",   &t_TrkNuIsolFlag);
   }
   if(doChgIsolTree) {
     ChgIsolnTree = fs->make<TTree>("ChgIsolnTree", "ChgIsolntree");
     t_PixcandP        = new std::vector<double>();
     t_PixcandPt       = new std::vector<double>();
     t_PixcandEta      = new std::vector<double>();
     t_PixcandPhi      = new std::vector<double>();
     t_PixcandMaxP     = new std::vector<std::vector<double> >();
     t_PixTrkcandP     = new std::vector<double>();
     t_PixTrkcandPt    = new std::vector<double>();
     t_PixTrkcandEta   = new std::vector<double>();
     t_PixTrkcandPhi   = new std::vector<double>();
     t_PixTrkcandMaxP  = new std::vector<double>();
     t_PixTrkcandselTk = new std::vector<bool>();
     
     ChgIsolnTree->Branch("t_PixcandP",       "std::vector<double>", &t_PixcandP);
     ChgIsolnTree->Branch("t_PixcandPt",      "std::vector<double>", &t_PixcandPt);
     ChgIsolnTree->Branch("t_PixcandEta",     "std::vector<double>", &t_PixcandEta);
     ChgIsolnTree->Branch("t_PixcandPhi",     "std::vector<double>", &t_PixcandPhi);
     ChgIsolnTree->Branch("t_PixcandMaxP",    "std::vector<std::vector<double> >", &t_PixcandMaxP);
     ChgIsolnTree->Branch("t_PixTrkcandP",    "std::vector<double>", &t_PixTrkcandP);
     ChgIsolnTree->Branch("t_PixTrkcandPt",   "std::vector<double>", &t_PixTrkcandPt  );
     ChgIsolnTree->Branch("t_PixTrkcandEta",  "std::vector<double>", &t_PixTrkcandEta );
     ChgIsolnTree->Branch("t_PixTrkcandPhi",  "std::vector<double>", &t_PixTrkcandPhi );
     ChgIsolnTree->Branch("t_PixTrkcandMaxP", "std::vector<double>", &t_PixTrkcandMaxP);
     ChgIsolnTree->Branch("t_PixTrkcandselTk","std::vector<bool>", &t_PixTrkcandselTk);
   }
   if(doMipCutTree) {
     MipCutTree = fs->make<TTree>("MipCutTree", "MipCuttree");
     t_NFcandP        = new std::vector<double>();
     t_NFcandPt       = new std::vector<double>();
     t_NFcandEta      = new std::vector<double>();
     t_NFcandPhi      = new std::vector<double>();
     t_NFcandEmip     = new std::vector<double>();
     t_NFTrkcandP     = new std::vector<double>();
     t_NFTrkcandPt    = new std::vector<double>();
     t_NFTrkcandEta   = new std::vector<double>();
     t_NFTrkcandPhi   = new std::vector<double>();
     t_NFTrkcandEmip  = new std::vector<double>();
     t_NFTrkMinDR     = new std::vector<double>();
     t_NFTrkMinDP1    = new std::vector<double>();
     t_NFTrkselTkFlag = new std::vector<bool>();
     t_NFTrkqltyFlag  = new std::vector<bool>();
     t_NFTrkMissFlag  = new std::vector<bool>();
     t_NFTrkPVFlag    = new std::vector<bool>();
     t_NFTrkPropFlag  = new std::vector<bool>();
     t_NFTrkChgIsoFlag= new std::vector<bool>();
     t_NFTrkNeuIsoFlag= new std::vector<bool>();
     t_NFTrkMipFlag   = new std::vector<bool>();
     t_ECone          = new std::vector<double>();
     
     MipCutTree->Branch("t_NFcandP",        "std::vector<double>", &t_NFcandP);
     MipCutTree->Branch("t_NFcandPt",       "std::vector<double>", &t_NFcandPt);
     MipCutTree->Branch("t_NFcandEta",      "std::vector<double>", &t_NFcandEta);
     MipCutTree->Branch("t_NFcandPhi",      "std::vector<double>", &t_NFcandPhi);
     MipCutTree->Branch("t_NFcandEmip",     "std::vector<double>", &t_NFcandEmip);
     MipCutTree->Branch("t_NFTrkcandP",     "std::vector<double>", &t_NFTrkcandP);
     MipCutTree->Branch("t_NFTrkcandPt",    "std::vector<double>", &t_NFTrkcandPt  );
     MipCutTree->Branch("t_NFTrkcandEta",   "std::vector<double>", &t_NFTrkcandEta );
     MipCutTree->Branch("t_NFTrkcandPhi",   "std::vector<double>", &t_NFTrkcandPhi );
     MipCutTree->Branch("t_NFTrkcandEmip",  "std::vector<double>", &t_NFTrkcandEmip);
     MipCutTree->Branch("t_NFTrkMinDR",     "std::vector<double>", &t_NFTrkMinDR);
     MipCutTree->Branch("t_NFTrkMinDP1",    "std::vector<double>", &t_NFTrkMinDP1);
     MipCutTree->Branch("t_NFTrkselTkFlag", "std::vector<bool>",   &t_NFTrkselTkFlag);
     MipCutTree->Branch("t_NFTrkqltyFlag",  "std::vector<bool>",   &t_NFTrkqltyFlag);
     MipCutTree->Branch("t_NFTrkMissFlag",  "std::vector<bool>",   &t_NFTrkMissFlag);
     MipCutTree->Branch("t_NFTrkPVFlag",    "std::vector<bool>",   &t_NFTrkPVFlag);
     MipCutTree->Branch("t_NFTrkPropFlag",  "std::vector<bool>",   &t_NFTrkPropFlag);
     MipCutTree->Branch("t_NFTrkChgIsoFlag","std::vector<bool>",   &t_NFTrkChgIsoFlag);
     MipCutTree->Branch("t_NFTrkNeuIsoFlag","std::vector<bool>",   &t_NFTrkNeuIsoFlag);
     MipCutTree->Branch("t_NFTrkMipFlag",   "std::vector<bool>",   &t_NFTrkMipFlag);
     MipCutTree->Branch("t_ECone",          "std::vector<double>", &t_ECone);
   }
   h_Filters     = fs->make<TH1I>("h_Filters", "Filter Accepts", 10, 0, 10);
   std::string FilterNames[7] = {"hltL1sL1SingleJet68", "hltIsolPixelTrackL2FilterHE", "ecalIsolPixelTrackFilterHE", "hltIsolPixelTrackL3FilterHE",
                 "hltIsolPixelTrackL2FilterHB", "ecalIsolPixelTrackFilterHB", "hltIsolPixelTrackL3FilterHB"};
   for(int i=0; i<7; i++) {
     h_Filters->GetXaxis()->SetBinLabel(i+1, FilterNames[i].c_str());
   }
 
   h_nHLT        = fs->make<TH1I>("h_nHLT" , "Size of trigger Names", 1000, 1, 1000);
   h_HLT         = fs->make<TH1I>("h_HLT"  ,  "HLT accept", 3, -1, 2);
   h_PreL1       = fs->make<TH1I>("h_PreL1",  "L1 Prescale", 500, 0, 500);
   h_PreHLT      = fs->make<TH1I>("h_PreHLT", "HLT Prescale", 50, 0, 50);
   h_Pre         = fs->make<TH1I>("h_Pre",    "Prescale", 3000, 0, 3000);
 
   h_PreL1wt     = fs->make<TH1D>("h_PreL1wt", "Weighted L1 Prescale", 500, 0, 500);
   h_PreHLTwt    = fs->make<TH1D>("h_PreHLTwt", "Weighted HLT Prescale", 500, 0, 100);
   h_L1ObjEnergy = fs->make<TH1D>("h_L1ObjEnergy", "Energy of L1Object", 500, 0.0, 500.0);
 
   h_EnIn = fs->make<TH1D>("h_EnInEcal", "EnergyIn Ecal", 200, 0.0, 20.0);
   h_EnOut = fs->make<TH1D>("h_EnOutEcal", "EnergyOut Ecal", 200, 0.0, 20.0);
   h_MipEnMatch = fs->make<TH2D>("h_MipEnMatch", "MipEn: HLT level vs Reco Level (Matched)", 200, 0.0, 20.0, 200, 0.0, 20.0);
   h_MipEnNoMatch = fs->make<TH2D>("h_MipEnNoMatch", "MipEn: HLT level vs Reco Level (No Match Found)", 200, 0.0, 20.0, 200, 0.0, 20.0);
 
   if (doL2L3) {
     h_nL3Objs     = fs->make<TH1I>("h_nL3Objs", "Number of L3 objects", 10, 0, 10);
     
     std::string pairs[9] = {"L2L3", "L2L3Match", "L2L3NoMatch", "L3Reco", "L3RecoMatch", "L3RecoNoMatch", "NewFilterReco", "NewFilterRecoMatch", "NewFilterRecoNoMatch"};
     for (int ipair=0; ipair<9; ipair++) {
       sprintf(hname, "h_dEta%s", pairs[ipair].c_str());
       sprintf(htit, "#Delta#eta for %s", pairs[ipair].c_str());
       h_dEta[ipair]        = fs->make<TH1D>(hname, htit, 200, -10.0, 10.0);
       h_dEta[ipair]->GetXaxis()->SetTitle("d#eta");
       
       sprintf(hname, "h_dPhi%s", pairs[ipair].c_str());
       sprintf(htit, "#Delta#phi for %s", pairs[ipair].c_str());
       h_dPhi[ipair]        = fs->make<TH1D>(hname, htit, 140, -7.0, 7.0);
       h_dPhi[ipair]->GetXaxis()->SetTitle("d#phi");
 
       sprintf(hname, "h_dPt%s", pairs[ipair].c_str());
       sprintf(htit, "#Delta dp_{T} for %s objects", pairs[ipair].c_str());
       h_dPt[ipair]         = fs->make<TH1D>(hname, htit, 400, -200.0, 200.0);
       h_dPt[ipair]->GetXaxis()->SetTitle("dp_{T} (GeV)");
 
       sprintf(hname, "h_dP%s", pairs[ipair].c_str());
       sprintf(htit, "#Delta p for %s objects", pairs[ipair].c_str());
       h_dP[ipair]         = fs->make<TH1D>(hname, htit, 400, -200.0, 200.0);
       h_dP[ipair]->GetXaxis()->SetTitle("dP (GeV)");
 
       sprintf(hname, "h_dinvPt%s", pairs[ipair].c_str());
       sprintf(htit, "#Delta (1/p_{T}) for %s objects", pairs[ipair].c_str());
       h_dinvPt[ipair]      = fs->make<TH1D>(hname, htit, 500, -0.4, 0.1);
       h_dinvPt[ipair]->GetXaxis()->SetTitle("d(1/p_{T})");
       sprintf(hname, "h_mindR%s", pairs[ipair].c_str());
       sprintf(htit, "min(#Delta R) for %s objects", pairs[ipair].c_str());
       h_mindR[ipair]       = fs->make<TH1D>(hname, htit, 500, 0.0, 1.0);
       h_mindR[ipair]->GetXaxis()->SetTitle("dR");
     }
 
     for (int lvl=0; lvl<2; lvl++) {
       sprintf(hname, "h_dEtaL1%s", levels[lvl+1].c_str());
       sprintf(htit, "#Delta#eta for L1 and %s objects", levels[lvl+1].c_str());
       h_dEtaL1[lvl] = fs->make<TH1D>(hname, htit, 400, -10.0, 10.0);
 
       sprintf(hname, "h_dPhiL1%s", levels[lvl+1].c_str());
       sprintf(htit, "#Delta#phi for L1 and %s objects", levels[lvl+1].c_str());
       h_dPhiL1[lvl] = fs->make<TH1D>(hname, htit, 280, -7.0, 7.0);
 
       sprintf(hname, "h_dRL1%s", levels[lvl+1].c_str());
       sprintf(htit, "#Delta R for L1 and %s objects", levels[lvl+1].c_str());
       h_dRL1[lvl] = fs->make<TH1D>(hname, htit, 100, 0.0, 10.0);
     }
   }
 
   int levmin = (doL2L3 ? 0 : 10);
   for (int ilevel=levmin; ilevel<20; ilevel++) {
     sprintf(hname, "h_p%s", levels[ilevel].c_str());
     sprintf(htit, "p for %s objects", levels[ilevel].c_str());
     h_p[ilevel] = fs->make<TH1D>(hname, htit, 100, 0.0, 500.0);
     h_p[ilevel]->GetXaxis()->SetTitle("p (GeV)");
     
     sprintf(hname, "h_pt%s", levels[ilevel].c_str());
     sprintf(htit, "p_{T} for %s objects", levels[ilevel].c_str());
     h_pt[ilevel] = fs->make<TH1D>(hname, htit, 100, 0.0, 500.0);
     h_pt[ilevel]->GetXaxis()->SetTitle("p_{T} (GeV)");
     
     sprintf(hname, "h_eta%s", levels[ilevel].c_str());
     sprintf(htit, "#eta for %s objects", levels[ilevel].c_str());
     h_eta[ilevel] = fs->make<TH1D>(hname, htit, 100, -5.0, 5.0);
     h_eta[ilevel]->GetXaxis()->SetTitle("#eta");
     
     sprintf(hname, "h_phi%s", levels[ilevel].c_str());
     sprintf(htit, "#phi for %s objects", levels[ilevel].c_str());
     h_phi[ilevel] = fs->make<TH1D>(hname, htit, 70, -3.5, 3.50);
     h_phi[ilevel]->GetXaxis()->SetTitle("#phi");
   }
   
   std::string cuts[2]  = {"HLTMatched", "HLTNotMatched"};
   std::string cuts2[2] = {"All", "Away from L1"};
   for (int icut=0; icut<2; icut++) {
     sprintf(hname, "h_eMip%s", cuts[icut].c_str());
     sprintf(htit, "eMip for %s tracks", cuts[icut].c_str());
     h_eMip[icut]     =fs->make<TH1D>(hname, htit, 200, 0.0, 10.0);
     h_eMip[icut]->GetXaxis()->SetTitle("E_{Mip} (GeV)");
 
     sprintf(hname, "h_eMaxNearP%s", cuts[icut].c_str());
     sprintf(htit, "eMaxNearP for %s tracks", cuts[icut].c_str());
     h_eMaxNearP[icut]=fs->make<TH1D>(hname, htit, 240, -2.0, 10.0);
     h_eMaxNearP[icut]->GetXaxis()->SetTitle("E_{MaxNearP} (GeV)");
 
     sprintf(hname, "h_eNeutIso%s", cuts[icut].c_str());
     sprintf(htit, "eNeutIso for %s ", cuts[icut].c_str());
     h_eNeutIso[icut] =fs->make<TH1D>(hname, htit, 200, 0.0, 10.0);
     h_eNeutIso[icut]->GetXaxis()->SetTitle("E_{NeutIso} (GeV)");
 
     for (int kcut=0; kcut<2; ++kcut) {
       for (int lim=0; lim<5; ++lim) {
     sprintf(hname, "h_etaCalibTracks%sCut%dLim%d", cuts[icut].c_str(), kcut, lim);
     sprintf(htit, "#eta for %s isolated MIP tracks (%4.1f < p < %5.1f Gev/c %s)", cuts[icut].c_str(), pLimits[lim], pLimits[lim+1], cuts2[kcut].c_str());
     h_etaCalibTracks[lim][icut][kcut]=fs->make<TH1D>(hname, htit, 60, -30.0, 30.0);
     h_etaCalibTracks[lim][icut][kcut]->GetXaxis()->SetTitle("i#eta");
 
     sprintf(hname, "h_etaMipTracks%sCut%dLim%d", cuts[icut].c_str(), kcut, lim);
     sprintf(htit, "#eta for %s charge isolated MIP tracks (%4.1f < p < %5.1f Gev/c %s)", cuts[icut].c_str(), pLimits[lim], pLimits[lim+1], cuts2[kcut].c_str());
     h_etaMipTracks[lim][icut][kcut]=fs->make<TH1D>(hname, htit, 60, -30.0, 30.0);
     h_etaMipTracks[lim][icut][kcut]->GetXaxis()->SetTitle("i#eta");
       }
     }
   }
 
   std::string ecut1[3] = {"all","HLTMatched","HLTNotMatched"};
   std::string ecut2[2] = {"without","with"};
   int etac[48] = {-1,-2,-3,-4,-5,-6,-7,-8,-9,-10,-11,-12,-13,-14,-15,-16,-17,-18,-19,-20,-21,-22,-23,-24,
           1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24};
   for (int icut=0; icut<6; icut++) {
     //    int i1 = (icut>3 ? 1 : 0);
     int i1 = (icut>2 ? 1 : 0);
     int i2 = icut - i1*3;
     for (int kcut=0; kcut<48; kcut++) {
       for (int lim=0; lim<5; ++lim) {
     sprintf(hname, "h_eta%dEnHcal%s%s%d", etac[kcut], ecut1[i2].c_str(), ecut2[i1].c_str(), lim);
     sprintf(htit, "HCAL energy for #eta=%d for %s tracks (p=%4.1f:%5.1f Gev) %s neutral isolation", etac[kcut], ecut1[i2].c_str(), pLimits[lim], pLimits[lim+1], ecut2[i1].c_str());
     h_eHcal[lim][icut][kcut]=fs->make<TH1D>(hname, htit, 750, 0.0, 150.0);
     h_eHcal[lim][icut][kcut]->GetXaxis()->SetTitle("Energy (GeV)");
     sprintf(hname, "h_eta%dEnCalo%s%s%d", etac[kcut], ecut1[i2].c_str(), ecut2[i1].c_str(), lim);
     sprintf(htit, "Calorimter energy for #eta=%d for %s tracks (p=%4.1f:%5.1f Gev) %s neutral isolation", etac[kcut], ecut1[i2].c_str(), pLimits[lim], pLimits[lim+1], ecut2[i1].c_str());
     h_eCalo[lim][icut][kcut]=fs->make<TH1D>(hname, htit, 750, 0.0, 150.0);
     h_eCalo[lim][icut][kcut]->GetXaxis()->SetTitle("Energy (GeV)");
       }
     }
   }
 }

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

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 742 of file IsoTrig.cc.

742 {}

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

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 733 of file IsoTrig.cc.

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

                                                                       {
   std::cout  << "Run " << iRun.run() << " hltconfig.init " 
          << hltConfig_.init(iRun,iSetup,processName,changed) << std::endl;;
 }

void IsoTrig::chgIsolation	(	double &	etaTriggered,
		double &	phiTriggered,
		edm::Handle< reco::TrackCollection > &	trkCollection,
		const edm::Event &	theEvent
	)

private

Definition at line 1309 of file IsoTrig.cc.

References a_charIsoR, ecalTB2006H4_GenSimDigiReco_cfg::bField, bFieldH, spr::chargeIsolationCone(), ChgIsolnTree, clearChgIsolnTreeVectors(), gather_cfg::cout, deltaR(), dR(), eta(), getDistInCM(), GetEtaPhiAtEcal(), spr::goodTrack(), i, j, roll_playback::k, prof2calltree::l, leadPV, minPTrackValue_, AlCaHLTBitMon_ParallelJobs::p, reco::TrackBase::p(), pG, phi, pixelIsolationConeSizeAtEC_, pixelTrackRefsHE, position, prelimCone_, edm::ESHandle< class >::product(), spr::propagateCALO(), EnergyCorrector::pt, pushChgIsolnTreeVecs(), reco::TrackBase::px(), reco::TrackBase::py(), reco::TrackBase::pz(), dttmaxenums::R, recVtxs, edm::second(), selectionParameters, tauUnbiasCone_, theTrackQuality, verbosity, vtxCutIsol_, vtxCutSeed_, and detailsBasic3DVector::z.

Referenced by analyze().

                                          {
   clearChgIsolnTreeVectors();
   if (verbosity%10>0)  std::cout << "Inside chgIsolation() with eta/phi Triggered: " << etaTriggered << "/" << phiTriggered << std::endl;
   std::vector<double> maxP;
   
   std::vector<spr::propagatedTrackDirection>::const_iterator trkDetItr; 
   const MagneticField *bField = bFieldH.product();
   const CaloGeometry* geo = pG.product();
   std::vector<spr::propagatedTrackDirection> trkCaloDirections1;
   spr::propagateCALO(trkCollection, geo, bField, theTrackQuality, trkCaloDirections1, ((verbosity/100)%10>2));
   if (verbosity%10>0) std::cout << "Propagated TrkCollection" << std::endl;
   for (unsigned int k=0; k<pixelIsolationConeSizeAtEC_.size(); ++k)
     maxP.push_back(0);
   unsigned i = pixelTrackRefsHE.size();
   std::vector<std::pair<unsigned int, std::pair<double, double>>>  VecSeedsatEC;
   //loop to select isolated tracks
   for (unsigned iS=0; iS<pixelTrackRefsHE.size(); iS++) {
     if (pixelTrackRefsHE[iS]->p() > minPTrackValue_) {
       bool vtxMatch = false;
       //associate to vertex (in Z) 
       unsigned int ivSel = recVtxs->size();
       double minDZ = 100;
       for (unsigned int iv = 0; iv < recVtxs->size(); ++iv) {
     if (fabs(pixelTrackRefsHE[iS]->dz((*recVtxs)[iv].position()))<minDZ) {
       minDZ  = fabs(pixelTrackRefsHE[iS]->dz((*recVtxs)[iv].position()));
       ivSel  = iv;
     }
       }
       //cut on dYX:
       if (ivSel == recVtxs->size()) {
     vtxMatch = true;
       } else if (fabs(pixelTrackRefsHE[iS]->dxy((*recVtxs)[ivSel].position()))<vtxCutSeed_){
     vtxMatch = true;
       }
       //select tracks not matched to triggered L1 jet
       double R = deltaR(etaTriggered, phiTriggered, pixelTrackRefsHE[iS]->eta(), 
             pixelTrackRefsHE[iS]->phi());
       if (R > tauUnbiasCone_ && vtxMatch) {
     //propagate seed track to ECAL surface:
     std::pair<double,double> seedCooAtEC;
     // in case vertex is found:
     if (minDZ!=100) seedCooAtEC=GetEtaPhiAtEcal(pixelTrackRefsHE[iS]->eta(), pixelTrackRefsHE[iS]->phi(), pixelTrackRefsHE[iS]->pt(), pixelTrackRefsHE[iS]->charge(), (*recVtxs)[ivSel].z());
     //in case vertex is not found:
     else            seedCooAtEC=GetEtaPhiAtEcal(pixelTrackRefsHE[iS]->eta(), pixelTrackRefsHE[iS]->phi(), pixelTrackRefsHE[iS]->pt(), pixelTrackRefsHE[iS]->charge(), 0);
     VecSeedsatEC.push_back(std::make_pair(iS, seedCooAtEC));
       }
     }
   }
 //  std::cout << " index of the matched L2cand: " << i << " eta/phi/p: " << pixelTrackRefs[i]->eta()
 //      << "/" << pixelTrackRefs[i]->phi() << "/" << pixelTrackRefs[i]->p() << std::endl;;
 //  std::cout << " VecSeedsatEC.size() " << VecSeedsatEC.size() << " pixelTrackRefs.size() " << pixelTrackRefs.size() << std::endl;
   for (unsigned int l=0; l<VecSeedsatEC.size(); l++) {
     unsigned int iSeed =  VecSeedsatEC[l].first;
     //    std::cout << l << "th Pixcandidate iseed:" << iSeed << std::endl; 
     math::XYZTLorentzVector v1(pixelTrackRefsHE[iSeed]->px(),pixelTrackRefsHE[iSeed]->py(),
                    pixelTrackRefsHE[iSeed]->pz(),pixelTrackRefsHE[iSeed]->p());
 
     for (unsigned int j=0; j<VecSeedsatEC.size(); j++) {
       unsigned int iSurr = VecSeedsatEC[j].first;
       if (iSeed != iSurr) {
     //define preliminary cone around seed track impact point from which tracks will be extrapolated:
     //  edm::Ref<reco::IsolatedPixelTrackCandidateCollection> cand2ref =
     //    edm::Ref<reco::IsolatedPixelTrackCandidateCollection>(L2cands, iSurr);
     if (deltaR(pixelTrackRefsHE[iSeed]->eta(), pixelTrackRefsHE[iSeed]->phi(), pixelTrackRefsHE[iSurr]->eta(), 
            pixelTrackRefsHE[iSurr]->phi()) < prelimCone_) {
       unsigned int ivSel = recVtxs->size();
       double minDZ2=100;
       for (unsigned int iv = 0; iv < recVtxs->size(); ++iv) {
         if (fabs(pixelTrackRefsHE[iSurr]->dz((*recVtxs)[iv].position()))<minDZ2) {
           minDZ2  = fabs(pixelTrackRefsHE[iSurr]->dz((*recVtxs)[iv].position()));
           ivSel   = iv;
         }
       }
       //      std::cout << "minDZ2 " << minDZ2 << " dxy " << pixelTrackRefs[iSurr]->dxy((*recVtxs)[ivSel].position()) << " vtxCutIsol_ " << vtxCutIsol_ <<std::endl; 
       //cut ot dXY:
       if (minDZ2==100 || fabs(pixelTrackRefsHE[iSurr]->dxy((*recVtxs)[ivSel].position()))<vtxCutIsol_) {
         //calculate distance at ECAL surface and update isolation: 
         double dist = getDistInCM(VecSeedsatEC[i].second.first, VecSeedsatEC[i].second.second, VecSeedsatEC[j].second.first, VecSeedsatEC[j].second.second);
         //      std::cout << "pixelIsolationConeSizeAtEC_.size() " <<pixelIsolationConeSizeAtEC_.size() << std::endl;
         for (unsigned int k=0; k<pixelIsolationConeSizeAtEC_.size(); ++k) {
           if (dist<pixelIsolationConeSizeAtEC_[k]) {
         if (pixelTrackRefsHE[iSurr]->p() > maxP[k]) 
           maxP[k] = pixelTrackRefsHE[iSurr]->p();
           }
         }
       }
     }
       }
     }
   
     double conehmaxNearP = -1; bool selectTk=false;
     double mindR=999.9; int nTracks=0;
     math::XYZTLorentzVector mindRvec;
     for (trkDetItr = trkCaloDirections1.begin(); trkDetItr != trkCaloDirections1.end(); trkDetItr++, nTracks++){
       int nNearTRKs=0;
       const reco::Track* pTrack = &(*(trkDetItr->trkItr));
       math::XYZTLorentzVector v2(pTrack->px(), pTrack->py(),  
                  pTrack->pz(), pTrack->p());
       double dr   = dR(v1,v2);
       //      double dp1  = std::abs(1./v1.r() - 1./v2.r());
       //    std::cout <<"THis recotrack(eta/phi/pt) " << pTrack->eta() << "/" 
       //          << pTrack->phi() << "/" << pTrack->pt() << " has dr/dp= " 
       //          << dr << "/" << dp << "/" << dp1 << std::endl;
       if (dr<mindR) {
     selectTk = spr::goodTrack(pTrack,leadPV,selectionParameters,((verbosity/100)%10>1));
     conehmaxNearP = spr::chargeIsolationCone(nTracks, trkCaloDirections1, a_charIsoR, nNearTRKs, ((verbosity/100)%10>1));
     mindR = dr;
     mindRvec = v2;
       }
     }
     pushChgIsolnTreeVecs(v1, mindRvec, maxP, conehmaxNearP, selectTk);
     //    std::cout << " conehmaxNearP " << conehmaxNearP << std::endl;
     //    for (unsigned int k=0; k<pixelIsolationConeSizeAtEC_.size(); ++k) 
     //      std::cout << "k/maxP[k]: " << k << "/" << maxP[k] << std::endl;   
   }
   ChgIsolnTree->Fill();
 }

void IsoTrig::clearChgIsolnTreeVectors ( )

private

Definition at line 363 of file IsoTrig.cc.

References hitfit::clear(), i, t_PixcandEta, t_PixcandMaxP, t_PixcandP, t_PixcandPhi, t_PixcandPt, t_PixTrkcandEta, t_PixTrkcandMaxP, t_PixTrkcandP, t_PixTrkcandPhi, t_PixTrkcandPt, and t_PixTrkcandselTk.

Referenced by chgIsolation().

                                        { 
   t_PixcandP        ->clear();
   t_PixcandPt       ->clear();
   t_PixcandEta      ->clear();
   t_PixcandPhi      ->clear();
   for(unsigned int i=0; i< t_PixcandMaxP->size(); i++)
     t_PixcandMaxP[i].clear();
   t_PixcandMaxP     ->clear();
   t_PixTrkcandP     ->clear();
   t_PixTrkcandPt    ->clear();
   t_PixTrkcandEta   ->clear();
   t_PixTrkcandPhi   ->clear();
   t_PixTrkcandMaxP  ->clear();
   t_PixTrkcandselTk  ->clear();
 }

void IsoTrig::clearMipCutTreeVectors ( )

private

Definition at line 379 of file IsoTrig.cc.

References t_ECone, t_NFcandEmip, t_NFcandEta, t_NFcandP, t_NFcandPhi, t_NFcandPt, t_NFTrkcandEmip, t_NFTrkcandEta, t_NFTrkcandP, t_NFTrkcandPhi, t_NFTrkcandPt, t_NFTrkChgIsoFlag, t_NFTrkMinDP1, t_NFTrkMinDR, t_NFTrkMipFlag, t_NFTrkMissFlag, t_NFTrkNeuIsoFlag, t_NFTrkPropFlag, t_NFTrkPVFlag, t_NFTrkqltyFlag, and t_NFTrkselTkFlag.

Referenced by studyMipCut().

                                      {
   t_NFcandP        ->clear();
   t_NFcandPt       ->clear();
   t_NFcandEta      ->clear();
   t_NFcandPhi      ->clear();
   t_NFcandEmip     ->clear();
   t_NFTrkcandP     ->clear();
   t_NFTrkcandPt    ->clear();
   t_NFTrkcandEta   ->clear();
   t_NFTrkcandPhi   ->clear();
   t_NFTrkcandEmip  ->clear();
   t_NFTrkMinDR     ->clear();
   t_NFTrkMinDP1    ->clear();
   t_NFTrkselTkFlag ->clear();
   t_NFTrkqltyFlag  ->clear();
   t_NFTrkMissFlag  ->clear();
   t_NFTrkPVFlag    ->clear();
   t_NFTrkPropFlag  ->clear();
   t_NFTrkChgIsoFlag->clear();
   t_NFTrkNeuIsoFlag->clear();
   t_NFTrkMipFlag   ->clear();
   t_ECone          ->clear();
 }

double IsoTrig::dEta	(	math::XYZTLorentzVector &	vec1,
		math::XYZTLorentzVector &	vec2
	)

private

Definition at line 1485 of file IsoTrig.cc.

Referenced by dR(), fillDifferences(), and studyTrigger().

                                                                              {
   return (vec1.eta()-vec2.eta());
 }

double IsoTrig::dinvPt	(	math::XYZTLorentzVector &	vec1,
		math::XYZTLorentzVector &	vec2
	)

private

Definition at line 1515 of file IsoTrig.cc.

Referenced by fillDifferences().

                                                                                {
   return ((1/vec1.pt())-(1/vec2.pt()));
 }

double IsoTrig::dP	(	math::XYZTLorentzVector &	vec1,
		math::XYZTLorentzVector &	vec2
	)

private

Definition at line 1511 of file IsoTrig.cc.

References funct::abs().

Referenced by fillDifferences().

                                                                            {
   return (std::abs(vec1.r()-vec2.r()));
 }

double IsoTrig::dPhi	(	math::XYZTLorentzVector &	vec1,
		math::XYZTLorentzVector &	vec2
	)

private

Definition at line 1489 of file IsoTrig.cc.

References M_PI.

Referenced by dR(), fillDifferences(), and studyTrigger().

                                                                              {
 
   double phi1 = vec1.phi();
   if (phi1 < 0) phi1 += 2.0*M_PI;
   double phi2 = vec2.phi();
   if (phi2 < 0) phi2 += 2.0*M_PI;
   double dphi = phi1-phi2;
   if (dphi > M_PI)       dphi -= 2.*M_PI;
   else if (dphi < -M_PI) dphi += 2.*M_PI;
   return dphi;
 }

double IsoTrig::dPt	(	math::XYZTLorentzVector &	vec1,
		math::XYZTLorentzVector &	vec2
	)

private

Definition at line 1507 of file IsoTrig.cc.

Referenced by fillDifferences().

                                                                             {
   return (vec1.pt()-vec2.pt());
 }

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

private

Definition at line 1501 of file IsoTrig.cc.

References dEta(), dPhi(), and mathSSE::sqrt().

Referenced by chgIsolation(), fillDifferences(), studyIsolation(), studyMipCut(), and studyTrigger().

                                                                            {
   double deta = dEta(vec1,vec2);
   double dphi = dPhi(vec1,vec2);
   return std::sqrt(deta*deta + dphi*dphi);
 }

void IsoTrig::endJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 707 of file IsoTrig.cc.

References gather_cfg::cout, first, fs, g_Accepts, g_Pre, g_PreHLT, g_PreL1, i, TFileService::make(), maxRunNo, minRunNo, n, edm::second(), TrigList, trigNames, TrigPreList, and verbosity.

                      {
   unsigned int preL1, preHLT;
   std::map<unsigned int, unsigned int>::iterator itr;
   std::map<unsigned int, const std::pair<int, int>>::iterator itrPre;
   if (verbosity%10 > 0)  std::cout << trigNames.size() << "Triggers were run. RunNo vs HLT accepts for";
   for (unsigned int i=0; i<trigNames.size(); ++i) 
     std::cout << " " << trigNames[i];
   std::cout << std::endl;
   unsigned int n = maxRunNo - minRunNo +1;
   g_Pre = fs->make<TH1I>("h_PrevsRN", "PreScale Vs Run Number", n, minRunNo, maxRunNo);
   g_PreL1 = fs->make<TH1I>("h_PreL1vsRN", "L1 PreScale Vs Run Number", n, minRunNo, maxRunNo);
   g_PreHLT = fs->make<TH1I>("h_PreHLTvsRN", "HLT PreScale Vs Run Number", n, minRunNo, maxRunNo);
   g_Accepts = fs->make<TH1I>("h_HLTAcceptsvsRN", "HLT Accepts Vs Run Number", n, minRunNo, maxRunNo); 
 
   for (itr=TrigList.begin(), itrPre=TrigPreList.begin(); itr!=TrigList.end(); itr++, itrPre++) {
     preL1 = (itrPre->second).first;
     preHLT = (itrPre->second).second;
     std::cout << itr->first << " " << itr->second << " " <<  itrPre->first << " " << preL1 << " " << preHLT << std::endl;
     g_Accepts->Fill(itr->first, itr->second);
     g_PreL1->Fill(itr->first, preL1);
     g_PreHLT->Fill(itr->first, preHLT);
     g_Pre->Fill(itr->first, preL1*preHLT);
   }
 }

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

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 744 of file IsoTrig.cc.

744 {}

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

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 739 of file IsoTrig.cc.

739 {}

std::pair< double, double > IsoTrig::etaPhiTrigger ( )

private

Definition at line 1519 of file IsoTrig.cc.

References eta(), roll_playback::k, phi, EnergyCorrector::pt, and vec.

Referenced by analyze().

                                               {
   double eta(0), phi(0), ptmax(0);
   for (unsigned int k=0; k<vec[0].size(); ++k) {
     if (k == 0) {
       eta   = vec[0][k].eta();
       phi   = vec[0][k].phi();
       ptmax = vec[0][k].pt();
     } else if (vec[0][k].pt() > ptmax) {
       eta   = vec[0][k].eta();
       phi   = vec[0][k].phi();
       ptmax = vec[0][k].pt();
     }
   }
   return std::pair<double,double>(eta,phi);
 }

void IsoTrig::fillCuts	(	int	indx,
		double	eMipDR,
		double	conehmaxNearP,
		double	e_inCone,
		math::XYZTLorentzVector &	vec,
		int	ieta,
		bool	cut
	)

private

Definition at line 1453 of file IsoTrig.cc.

References cutCharge, cutMip, cutNeutral, h_eMaxNearP, h_eMip, h_eNeutIso, h_etaCalibTracks, h_etaMipTracks, and pLimits.

Referenced by studyIsolation().

                                                                                                                                      {
   h_eMip[indx]     ->Fill(eMipDR);
   h_eMaxNearP[indx]->Fill(conehmaxNearP);
   h_eNeutIso[indx] ->Fill(e_inCone);
   if ((conehmaxNearP < cutCharge) && (eMipDR < cutMip)) {
     for (int lim=0; lim<5; ++lim) {
       if ((vec.r()>pLimits[lim]) && (vec.r()<=pLimits[lim+1])) {
     h_etaMipTracks[lim][indx][0]->Fill((double)(ieta));
     if (cut) h_etaMipTracks[lim][indx][1]->Fill((double)(ieta));
     if (e_inCone < cutNeutral) {
       h_etaCalibTracks[lim][indx][0]->Fill((double)(ieta));
       if (cut) h_etaCalibTracks[lim][indx][1]->Fill((double)(ieta));
     }
       }
     }
   }
 }

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

static

Definition at line 145 of file IsoTrig.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);
 }

void IsoTrig::fillDifferences	(	int	indx,
		math::XYZTLorentzVector &	vec1,
		math::XYZTLorentzVector &	vec2,
		bool	debug
	)

private

Definition at line 1436 of file IsoTrig.cc.

References gather_cfg::cout, dEta(), dinvPt(), dP(), dPhi(), dPt(), dR(), h_dEta, h_dinvPt, h_dP, h_dPhi, h_dPt, and h_mindR.

Referenced by studyMipCut(), and studyTrigger().

                                                            {
   double dr     = dR(vec1,vec2);
   double deta   = dEta(vec1, vec2);
   double dphi   = dPhi(vec1, vec2);
   double dpt    = dPt(vec1, vec2);
   double dp     = dP(vec1, vec2);
   double dinvpt = dinvPt(vec1, vec2);
   h_dEta[indx]  ->Fill(deta);
   h_dPhi[indx]  ->Fill(dphi);
   h_dPt[indx]   ->Fill(dpt);
   h_dP[indx]    ->Fill(dp);
   h_dinvPt[indx]->Fill(dinvpt);
   h_mindR[indx] ->Fill(dr);
   if (debug) std::cout << "mindR for index " << indx << " is " << dr << " deta " << deta << " dphi " << dphi << " dpt " << dpt <<  " dinvpt " << dinvpt <<std::endl;
 }

void IsoTrig::fillEnergy	(	int	indx,
		int	ieta,
		double	hCone,
		double	eMipDR,
		math::XYZTLorentzVector &	vec
	)

private

Definition at line 1471 of file IsoTrig.cc.

References h_eCalo, h_eHcal, GetRecoTauVFromDQM_MC_cff::kk, and pLimits.

Referenced by studyIsolation().

                                                                                                     {
   int kk(-1);
   if      (ieta > 0 && ieta < 25)  kk = 23 + ieta;
   else if (ieta > -25 && ieta < 0) kk = -(ieta + 1);
   if (kk >= 0 && eMipDR > 0.01 && hCone > 1.0) {
     for (int lim=0; lim<5; ++lim) {
       if ((vec.r()>pLimits[lim]) && (vec.r()<=pLimits[lim+1])) {
     h_eHcal[lim][indx][kk]     ->Fill(hCone);
     h_eCalo[lim][indx][kk]     ->Fill(hCone+eMipDR);
       }
     }
   }
 }

void IsoTrig::fillHist	(	int	indx,
		math::XYZTLorentzVector &	vec
	)

private

Definition at line 1429 of file IsoTrig.cc.

References h_eta, h_p, h_phi, and h_pt.

Referenced by studyIsolation(), and studyTrigger().

                                                            {
   h_p[indx]->Fill(vec.r());
   h_pt[indx]->Fill(vec.pt());
   h_eta[indx]->Fill(vec.eta());
   h_phi[indx]->Fill(vec.phi());
 }

double IsoTrig::getDistInCM	(	double	eta1,
		double	phi1,
		double	eta2,
		double	phi2
	)

private

Definition at line 1594 of file IsoTrig.cc.

References angle(), funct::cos(), create_public_lumi_plots::exp, M_PI, rEB_, funct::sin(), funct::tan(), and zEE_.

Referenced by chgIsolation().

                                                                             {
   double Rec;
   double theta1=2*atan(exp(-eta1));
   double theta2=2*atan(exp(-eta2));
   if (fabs(eta1)<1.479) Rec=rEB_; //radius of ECAL barrel
   else if (fabs(eta1)>1.479&&fabs(eta1)<7.0) Rec=tan(theta1)*zEE_; //distance from IP to ECAL endcap
   else return 1000;
 
   //|vect| times tg of acos(scalar product)
   double angle=acos((sin(theta1)*sin(theta2)*(sin(phi1)*sin(phi2)+cos(phi1)*cos(phi2))+cos(theta1)*cos(theta2)));
   if (angle<0.5*M_PI) return fabs((Rec/sin(theta1))*tan(angle));
   else return 1000;
 }

std::pair< double, double > IsoTrig::GetEtaPhiAtEcal	(	double	etaIP,
		double	phiIP,
		double	pT,
		int	charge,
		double	vtxZ
	)

private

Definition at line 1535 of file IsoTrig.cc.

References alpha, bfVal, SiPixelRawToDigiRegional_cfi::deltaPhi, create_public_lumi_plots::exp, create_public_lumi_plots::log, M_PI, rEB_, funct::sin(), funct::tan(), theta(), detailsBasic3DVector::z, and zEE_.

Referenced by chgIsolation().

                                        {
 
   double deltaPhi=0;
   double etaEC = 100;
   double phiEC = 100;
 
   double Rcurv = 9999999;
   if (bfVal!=0) Rcurv=pT*33.3*100/(bfVal*10); //r(m)=pT(GeV)*33.3/B(kG)
 
   double ecDist = zEE_;
   double ecRad  = rEB_;  //radius of ECAL barrel (cm)
   double theta  = 2*atan(exp(-etaIP));
   double zNew   = 0;
   if (theta>0.5*M_PI) theta=M_PI-theta;
   if (fabs(etaIP)<1.479) {
     if ((0.5*ecRad/Rcurv)>1) {
       etaEC         = 10000;
       deltaPhi      = 0;
     } else {
       deltaPhi      =-charge*asin(0.5*ecRad/Rcurv);
       double alpha1 = 2*asin(0.5*ecRad/Rcurv);
       double z      = ecRad/tan(theta);
       if (etaIP>0) zNew = z*(Rcurv*alpha1)/ecRad+vtxZ; //new z-coordinate of track
       else         zNew =-z*(Rcurv*alpha1)/ecRad+vtxZ; //new z-coordinate of track
       double zAbs=fabs(zNew);
       if (zAbs<ecDist) {
         etaEC    = -log(tan(0.5*atan(ecRad/zAbs)));
         deltaPhi = -charge*asin(0.5*ecRad/Rcurv);
       }
       if (zAbs>ecDist) {
         zAbs           = (fabs(etaIP)/etaIP)*ecDist;
         double Zflight = fabs(zAbs-vtxZ);
         double alpha   = (Zflight*ecRad)/(z*Rcurv);
         double Rec     = 2*Rcurv*sin(alpha/2);
         deltaPhi       =-charge*alpha/2;
         etaEC          =-log(tan(0.5*atan(Rec/ecDist)));
       }
     }
   } else {
     zNew           = (fabs(etaIP)/etaIP)*ecDist;
     double Zflight = fabs(zNew-vtxZ);
     double Rvirt   = fabs(Zflight*tan(theta));
     double Rec     = 2*Rcurv*sin(Rvirt/(2*Rcurv));
     deltaPhi       =-(charge)*(Rvirt/(2*Rcurv));
     etaEC          =-log(tan(0.5*atan(Rec/ecDist)));
   }
 
   if (zNew<0) etaEC=-etaEC;
   phiEC            = phiIP+deltaPhi;
 
   if (phiEC<-M_PI) phiEC += 2*M_PI;
   if (phiEC>M_PI)  phiEC -= 2*M_PI;
 
   std::pair<double,double> retVal(etaEC,phiEC);
   return retVal;
 }

void IsoTrig::pushChgIsolnTreeVecs	(	math::XYZTLorentzVector &	Pixcand,
		math::XYZTLorentzVector &	Trkcand,
		std::vector< double > &	PixMaxP,
		double &	TrkMaxP,
		bool &	selTk
	)

private

Definition at line 403 of file IsoTrig.cc.

References t_PixcandEta, t_PixcandMaxP, t_PixcandP, t_PixcandPhi, t_PixcandPt, t_PixTrkcandEta, t_PixTrkcandMaxP, t_PixTrkcandP, t_PixTrkcandPhi, t_PixTrkcandPt, and t_PixTrkcandselTk.

Referenced by chgIsolation().

                                 {
   t_PixcandP         ->push_back(Pixcand.r());
   t_PixcandPt        ->push_back(Pixcand.pt());
   t_PixcandEta       ->push_back(Pixcand.eta());
   t_PixcandPhi       ->push_back(Pixcand.phi());
   t_PixcandMaxP      ->push_back(PixMaxP);
   t_PixTrkcandP      ->push_back(Trkcand.r());
   t_PixTrkcandPt     ->push_back(Trkcand.pt());
   t_PixTrkcandEta    ->push_back(Trkcand.eta());
   t_PixTrkcandPhi    ->push_back(Trkcand.phi());
   t_PixTrkcandMaxP   ->push_back(TrkMaxP); 
   t_PixTrkcandselTk  ->push_back(selTk); 
   
 }

void IsoTrig::pushMipCutTreeVecs	(	math::XYZTLorentzVector &	NFcand,
		math::XYZTLorentzVector &	Trkcand,
		double &	EmipNFcand,
		double &	EmipTrkcand,
		double &	mindR,
		double &	mindP1,
		std::vector< bool > &	Flags,
		double	hCone
	)

private

Definition at line 420 of file IsoTrig.cc.

References t_ECone, t_NFcandEmip, t_NFcandEta, t_NFcandP, t_NFcandPhi, t_NFcandPt, t_NFTrkcandEmip, t_NFTrkcandEta, t_NFTrkcandP, t_NFTrkcandPhi, t_NFTrkcandPt, t_NFTrkChgIsoFlag, t_NFTrkMinDP1, t_NFTrkMinDR, t_NFTrkMipFlag, t_NFTrkMissFlag, t_NFTrkNeuIsoFlag, t_NFTrkPropFlag, t_NFTrkPVFlag, t_NFTrkqltyFlag, and t_NFTrkselTkFlag.

Referenced by studyMipCut().

                                                        {
   t_NFcandP        ->push_back(NFcand.r());
   t_NFcandPt       ->push_back(NFcand.pt());
   t_NFcandEta      ->push_back(NFcand.eta());
   t_NFcandPhi      ->push_back(NFcand.phi());
   t_NFcandEmip     ->push_back(EmipNFcand);
   t_NFTrkcandP     ->push_back(Trkcand.r());
   t_NFTrkcandPt    ->push_back(Trkcand.pt());
   t_NFTrkcandEta   ->push_back(Trkcand.eta());
   t_NFTrkcandPhi   ->push_back(Trkcand.phi());
   t_NFTrkcandEmip  ->push_back(EmipTrkcand);
   t_NFTrkMinDR     ->push_back(mindR);
   t_NFTrkMinDP1    ->push_back(mindP1);
   t_NFTrkselTkFlag ->push_back(Flags[0]);
   t_NFTrkqltyFlag  ->push_back(Flags[1]);
   t_NFTrkMissFlag  ->push_back(Flags[2]);
   t_NFTrkPVFlag    ->push_back(Flags[3]);
   t_NFTrkPropFlag  ->push_back(Flags[4]);
   t_NFTrkChgIsoFlag->push_back(Flags[5]);
   t_NFTrkNeuIsoFlag->push_back(Flags[6]);
   t_NFTrkMipFlag   ->push_back(Flags[7]);
   t_ECone          ->push_back(hCone);
 }

void IsoTrig::studyIsolation	(	edm::Handle< reco::TrackCollection > &	trkCollection,
		std::vector< reco::TrackCollection::const_iterator > &	goodTks
	)

private

Definition at line 1195 of file IsoTrig.cc.

References a_charIsoR, a_coneR, a_mipR, a_neutR1, a_neutR2, barrelRecHitsHandle, ecalTB2006H4_GenSimDigiReco_cfg::bField, bFieldH, spr::chargeIsolationCone(), gather_cfg::cout, cutCharge, cutMip, cutNeutral, deltaR(), doL2L3, dR(), dr_L1, spr::eCone_ecal(), spr::eCone_hcal(), endcapRecHitsHandle, fillCuts(), fillEnergy(), fillHist(), spr::find(), spr::goodTrack(), hbhe, hit::id, HcalDetId::ieta(), edm::HandleBase::isValid(), roll_playback::k, leadPV, reco::TrackBase::p(), pG, edm::ESHandle< class >::product(), spr::propagateCALO(), reco::TrackBase::px(), reco::TrackBase::py(), reco::TrackBase::pz(), selectionParameters, theTrackQuality, vec, and verbosity.

Referenced by analyze().

                                                                         {
 
   if (trkCollection.isValid()) {
     // get handles to calogeometry and calotopology
     const CaloGeometry* geo = pG.product();
     const MagneticField *bField = bFieldH.product();
     std::vector<spr::propagatedTrackDirection> trkCaloDirections;
     spr::propagateCALO(trkCollection, geo, bField, theTrackQuality, trkCaloDirections, ((verbosity/100)%10>2));
 
     std::vector<spr::propagatedTrackDirection>::const_iterator trkDetItr;
   
     if ((verbosity/1000)%10 > 1) {
       std::cout << "n of barrelRecHitsHandle " << barrelRecHitsHandle->size() << std::endl;
       for (EcalRecHitCollection::const_iterator hit = barrelRecHitsHandle->begin(); hit != barrelRecHitsHandle->end(); ++hit) {
     std::cout << "hit : detid(ieta,iphi) " << (EBDetId)(hit->id()) << " time " << hit->time() << " energy " <<  hit->energy() << std::endl;
       }
       std::cout << "n of endcapRecHitsHandle " << endcapRecHitsHandle->size() << std::endl;
       for (EcalRecHitCollection::const_iterator hit = endcapRecHitsHandle->begin(); hit != endcapRecHitsHandle->end(); ++hit) {
     std::cout << "hit : detid(ieta,iphi) " << (EEDetId)(hit->id()) << " time " << hit->time() << " energy " <<  hit->energy() << std::endl;
       }
       std::cout << "n of hbhe " << hbhe->size() << std::endl;
       for (HBHERecHitCollection::const_iterator hit = hbhe->begin(); hit != hbhe->end(); ++hit) {
     std::cout << "hit : detid(ieta,iphi) " << hit->id() << " time " << hit->time() << " energy " <<  hit->energy() << std::endl;
       }
     }
     unsigned int nTracks=0, ngoodTk=0, nselTk=0;
     int          ieta=999;
     for (trkDetItr = trkCaloDirections.begin(); trkDetItr != trkCaloDirections.end(); trkDetItr++,nTracks++){
       bool l3Track  = (std::find(goodTks.begin(),goodTks.end(),trkDetItr->trkItr) != goodTks.end());
       const reco::Track* pTrack = &(*(trkDetItr->trkItr));
       math::XYZTLorentzVector v4(pTrack->px(), pTrack->py(), 
                  pTrack->pz(), pTrack->p());
       bool selectTk = spr::goodTrack(pTrack,leadPV,selectionParameters,((verbosity/100)%10>1));
       double eMipDR=9999., e_inCone=0, conehmaxNearP=0, mindR=999.9, hCone=0;
       if (trkDetItr->okHCAL) {
     HcalDetId detId = (HcalDetId)(trkDetItr->detIdHCAL);
     ieta = detId.ieta();
       }
       for (unsigned k=0; k<vec[0].size(); ++k) {
     double deltaR = dR(v4, vec[0][k]);
     if (deltaR<mindR) mindR = deltaR;
       }
       if ((verbosity/100)%10 > 1) std::cout << "Track ECAL " << trkDetItr->okECAL << " HCAL " << trkDetItr->okHCAL << " Flag " << selectTk << std::endl;
       if (selectTk && trkDetItr->okECAL && trkDetItr->okHCAL) {
     ngoodTk++;
     int nRH_eMipDR=0, nNearTRKs=0;
     double e1 =  spr::eCone_ecal(geo, barrelRecHitsHandle, endcapRecHitsHandle,
                      trkDetItr->pointHCAL, trkDetItr->pointECAL,
                      a_neutR1, trkDetItr->directionECAL, nRH_eMipDR);
     double e2 =  spr::eCone_ecal(geo, barrelRecHitsHandle, endcapRecHitsHandle,
                      trkDetItr->pointHCAL, trkDetItr->pointECAL,
                      a_neutR2, trkDetItr->directionECAL, nRH_eMipDR);
     eMipDR = spr::eCone_ecal(geo, barrelRecHitsHandle, endcapRecHitsHandle,
                  trkDetItr->pointHCAL, trkDetItr->pointECAL,
                  a_mipR, trkDetItr->directionECAL, nRH_eMipDR);
     conehmaxNearP = spr::chargeIsolationCone(nTracks, trkCaloDirections, a_charIsoR, nNearTRKs, ((verbosity/100)%10>1));
     e_inCone = e2 - e1;
     double             distFromHotCell=-99.0;
     int                nRecHitsCone=-99, ietaHotCell=-99, iphiHotCell=-99;
     GlobalPoint        gposHotCell(0.,0.,0.);
     std::vector<DetId> coneRecHitDetIds;
     hCone    = spr::eCone_hcal(geo, hbhe, trkDetItr->pointHCAL, 
                    trkDetItr->pointECAL,
                    a_coneR, trkDetItr->directionHCAL, 
                    nRecHitsCone, coneRecHitDetIds,
                    distFromHotCell, ietaHotCell, iphiHotCell,
                    gposHotCell, -1);
     if (eMipDR<1.0) nselTk++;
       }
       if (l3Track) {
     fillHist(10,v4);
     if (selectTk) {
       fillHist(11,v4);
       fillCuts(0, eMipDR, conehmaxNearP, e_inCone, v4, ieta, (mindR>dr_L1));
       if (conehmaxNearP < cutCharge) {
         fillHist(12,v4);
         if (eMipDR < cutMip) {
           fillHist(13,v4);
           fillEnergy(1, ieta, hCone, eMipDR, v4);
           fillEnergy(0, ieta, hCone, eMipDR, v4);
           if (e_inCone < cutNeutral) {
         fillHist(14, v4);
         fillEnergy(4, ieta, hCone, eMipDR, v4);
         fillEnergy(3, ieta, hCone, eMipDR, v4);
           }
         }
       }
     }
       } else if (doL2L3) {
     fillHist(15,v4);
     if (selectTk) {
       fillHist(16,v4);
       fillCuts(1, eMipDR, conehmaxNearP, e_inCone, v4, ieta, (mindR>dr_L1));
       if (conehmaxNearP < cutCharge) {
         fillHist(17,v4);
         if (eMipDR < cutMip) {
           fillHist(18,v4);
           fillEnergy(2, ieta, hCone, eMipDR, v4);
           fillEnergy(0, ieta, hCone, eMipDR, v4);
           if (e_inCone < cutNeutral) {
         fillHist(19, v4);
         fillEnergy(5, ieta, hCone, eMipDR, v4);
         fillEnergy(3, ieta, hCone, eMipDR, v4);
           }
         }
       }
     }
       }
     }
     //    std::cout << "Number of tracks selected offline " << nselTk << std::endl;      
   }
 }

void IsoTrig::studyMipCut	(	edm::Handle< reco::TrackCollection > &	trkCollection,
		edm::Handle< reco::IsolatedPixelTrackCandidateCollection > &	L2cands
	)

private

Definition at line 949 of file IsoTrig.cc.

References a_charIsoR, a_coneR, a_mipR, a_neutR1, a_neutR2, funct::abs(), barrelRecHitsHandle, ecalTB2006H4_GenSimDigiReco_cfg::bField, bFieldH, spr::chargeIsolationCone(), clearMipCutTreeVectors(), gather_cfg::cout, cutCharge, cutMip, cutNeutral, dR(), spr::eCone_ecal(), spr::eCone_hcal(), endcapRecHitsHandle, fillDifferences(), spr::goodTrack(), h_EnIn, h_EnOut, h_MipEnMatch, h_MipEnNoMatch, hbhe, i, edm::HandleBase::isValid(), leadPV, spr::trackSelectionParameters::maxDxyPV, spr::trackSelectionParameters::maxDzPV, spr::trackSelectionParameters::maxInMiss, spr::trackSelectionParameters::maxOutMiss, MipCutTree, reco::TrackBase::p(), pG, edm::ESHandle< class >::product(), spr::propagateCALO(), pushMipCutTreeVecs(), reco::TrackBase::px(), reco::TrackBase::py(), reco::TrackBase::pz(), selectionParameters, theTrackQuality, and verbosity.

Referenced by analyze().

                                                                              {
 
   clearMipCutTreeVectors();
   if (verbosity%10>0) std::cout << "inside studymipcut" << std::endl;
   if (!trkCollection.isValid()) {
     std::cout << "trkCollection.isValid is false" << std::endl;
   } else {
     std::vector<spr::propagatedTrackDirection>::const_iterator trkDetItr; 
     const MagneticField *bField = bFieldH.product();
     const CaloGeometry* geo = pG.product();
     std::vector<spr::propagatedTrackDirection> trkCaloDirections1;
     spr::propagateCALO(trkCollection, geo, bField, theTrackQuality, trkCaloDirections1, ((verbosity/100)%10>2));
     if (verbosity%10>0) std::cout << "Number of L2cands:" << L2cands->size() 
                      << " to be matched to something out of " 
                      << trkCaloDirections1.size() << " reco tracks" << std::endl;
     for (unsigned int i=0; i<L2cands->size(); i++) {    
       edm::Ref<reco::IsolatedPixelTrackCandidateCollection> candref =
     edm::Ref<reco::IsolatedPixelTrackCandidateCollection>(L2cands, i);  
       double enIn = candref->energyIn();
       h_EnIn->Fill(candref->energyIn());
       h_EnOut->Fill(candref->energyOut());
       math::XYZTLorentzVector v1(candref->track()->px(),candref->track()->py(),
                  candref->track()->pz(),candref->track()->p());
       if (verbosity%10>1) std::cout << "HLT Level Candidate eta/phi/pt/enIn:" 
                     << candref->track()->eta() << "/" << candref->track()->phi() 
                     << "/" << candref->track()->pt() << "/" << candref->energyIn() 
                     << std::endl;   
       math::XYZTLorentzVector mindRvec; 
       double mindR=999.9, mindP1=999.9, eMipDR=0.0;
       std::vector<bool> selFlags;
       unsigned int nTracks=0;
       double conehmaxNearP = 0, hCone=0;
       for (trkDetItr = trkCaloDirections1.begin(); trkDetItr != trkCaloDirections1.end(); trkDetItr++,nTracks++){
     const reco::Track* pTrack = &(*(trkDetItr->trkItr));
     math::XYZTLorentzVector v2(pTrack->px(), pTrack->py(),  
                    pTrack->pz(), pTrack->p());
     double dr   = dR(v1,v2);
     double dp1  = std::abs(1./v1.r() - 1./v2.r());
     //  std::cout <<"This recotrack(eta/phi/pt) " << pTrack->eta() << "/" 
     //        << pTrack->phi() << "/" << pTrack->pt() << " has dr/dp= " 
     //        << dr << "/" << dp << "/" << dp1 << std::endl;
     if (dr<mindR) {
       mindR = dr;
       mindP1= dp1;
       mindRvec=v2;
       int nRH_eMipDR=0, nNearTRKs=0;
 //    std::cout << "debug starts here ";
 //    std::cout << geo;
 //    std::cout << "_" << barrelRecHitsHandle;
 //    std::cout << "_" << endcapRecHitsHandle << "_" << std::endl;
       eMipDR = spr::eCone_ecal(geo, barrelRecHitsHandle, endcapRecHitsHandle,
                    trkDetItr->pointHCAL, trkDetItr->pointECAL,
                    a_mipR, trkDetItr->directionECAL, nRH_eMipDR);
       //      std::cout << "mip " << eMipDR << std::endl;
       bool selectTk = spr::goodTrack(pTrack,leadPV,selectionParameters,((verbosity/100)%10>1));
       spr::trackSelectionParameters oneCutParameters = selectionParameters;
       oneCutParameters.maxDxyPV  = 10;
       oneCutParameters.maxDzPV     = 100;
       oneCutParameters.maxInMiss = 2;
       oneCutParameters.maxOutMiss= 2;
       bool qltyFlag     =  spr::goodTrack(pTrack,leadPV,oneCutParameters,((verbosity/100)%10>1));
       oneCutParameters           = selectionParameters;
       oneCutParameters.maxDxyPV  = 10;
       oneCutParameters.maxDzPV     = 100;
       bool qltyMissFlag =  spr::goodTrack(pTrack,leadPV,oneCutParameters,((verbosity/100)%10>1));
       oneCutParameters           = selectionParameters;
       oneCutParameters.maxInMiss = 2;
       oneCutParameters.maxOutMiss= 2;
       bool qltyPVFlag   =  spr::goodTrack(pTrack,leadPV,oneCutParameters,((verbosity/100)%10>1));
 //    std::cout << "sel " << selectTk << std::endl;
 //    std::cout << "ntracks " << nTracks;
 //    std::cout << " a_charIsoR " << a_charIsoR;
 //    std::cout << " nNearTRKs " << nNearTRKs << std::endl; 
       conehmaxNearP = spr::chargeIsolationCone(nTracks, trkCaloDirections1, a_charIsoR, nNearTRKs, ((verbosity/100)%10>1));
       //      std::cout << "coneh " << conehmaxNearP << std::endl;
       //      std::cout << "ok " << trkDetItr->okECAL << " " << trkDetItr->okHCAL << std::endl;
       double e1 =  spr::eCone_ecal(geo, barrelRecHitsHandle, endcapRecHitsHandle,
                        trkDetItr->pointHCAL, trkDetItr->pointECAL,
                        a_neutR1, trkDetItr->directionECAL, nRH_eMipDR);
       //      std::cout << "e1 " << e1 << std::endl;
       double e2 =  spr::eCone_ecal(geo, barrelRecHitsHandle, endcapRecHitsHandle,
                        trkDetItr->pointHCAL, trkDetItr->pointECAL,
                        a_neutR2, trkDetItr->directionECAL, nRH_eMipDR);
       double e_inCone  = e2 - e1;
       bool chgIsolFlag = (conehmaxNearP < cutCharge);
       bool mipFlag     = (eMipDR < cutMip);
       bool neuIsolFlag = (e_inCone < cutNeutral);
       bool trkpropFlag = ((trkDetItr->okECAL) && (trkDetItr->okHCAL));
       selFlags.clear();
       selFlags.push_back(selectTk);        selFlags.push_back(qltyFlag);
       selFlags.push_back(qltyMissFlag);    selFlags.push_back(qltyPVFlag);
       selFlags.push_back(trkpropFlag);     selFlags.push_back(chgIsolFlag);
       selFlags.push_back(neuIsolFlag);     selFlags.push_back(mipFlag);
       //      std::cout << selectTk << " " << qltyFlag << " " << qltyMissFlag << " " << qltyPVFlag << " " << trkpropFlag << " " << chgIsolFlag << " " << neuIsolFlag << " " << mipFlag << std::endl;
     double             distFromHotCell=-99.0;
     int                nRecHitsCone=-99, ietaHotCell=-99, iphiHotCell=-99;
     GlobalPoint        gposHotCell(0.,0.,0.);
     std::vector<DetId> coneRecHitDetIds;
     hCone    = spr::eCone_hcal(geo, hbhe, trkDetItr->pointHCAL, 
                    trkDetItr->pointECAL,
                    a_coneR, trkDetItr->directionHCAL, 
                    nRecHitsCone, coneRecHitDetIds,
                    distFromHotCell, ietaHotCell, iphiHotCell,
                    gposHotCell, -1);
     }
       }
       pushMipCutTreeVecs(v1, mindRvec, enIn, eMipDR, mindR, mindP1, selFlags, hCone);
       fillDifferences(6, v1, mindRvec, (verbosity%10 >0));
       if (mindR>=0.05) {
     fillDifferences(8, v1, mindRvec, (verbosity%10 >0));
     // std::cout << "No Matching found between reco tracks and HLT Level Pixel Tracks: " << mindR << ":" << mindP1 << std::endl;
     h_MipEnNoMatch->Fill(candref->energyIn(), eMipDR);
       } else {
     fillDifferences(7, v1, mindRvec, (verbosity%10 >0));
     //  std::cout << "Reco Track seems to be macthing eta/phi/pt/eMipDR:" 
     //        << mindRvec.eta() << "/" << mindRvec.phi() << "/" 
     //        << mindRvec.pt() << "/" << eMipDR << " Match " << mindR
     //        << ":" << mindP1 << std::endl;   
     h_MipEnMatch->Fill(candref->energyIn(), eMipDR);
       }
     }
   }
   MipCutTree->Fill();  
 }

void IsoTrig::studyTiming ( const edm::Event & theEvent )

private

Definition at line 852 of file IsoTrig.cc.

References gather_cfg::cout, deltaR(), eta(), edm::Event::getByToken(), AlCaHLTBitMon_ParallelJobs::p, phi, pixelTrackRefsHB, pixelTrackRefsHE, EnergyCorrector::pt, FastTimerService::queryModuleTimeByLabel(), dttmaxenums::R, edm::Event::streamID(), t_nPixCand, t_nPixSeed, t_timeL2Prod, tauUnbiasCone_, TimingTree, tok_hlt_, tok_verthb_, tok_verthe_, trigger::TriggerL1CenJet, trigger::TriggerL1ForJet, trigger::TriggerL1TauJet, verbosity, and vtxCutSeed_.

Referenced by analyze().

                                                   {
   t_timeL2Prod->clear(); t_nPixCand->clear(); t_nPixSeed->clear();
 
   double tHB=0.0, tHE=0.0;
   int nCandHB=pixelTrackRefsHB.size(), nCandHE=pixelTrackRefsHE.size();
   int nSeedHB=0, nSeedHE=0;
 
   if(nCandHE>0) {
     edm::Handle<reco::VertexCollection> pVertHB, pVertHE;
     theEvent.getByToken(tok_verthb_,pVertHB);
     theEvent.getByToken(tok_verthe_,pVertHE);
     edm::Handle<trigger::TriggerFilterObjectWithRefs> l1trigobj;
     theEvent.getByToken(tok_hlt_, l1trigobj);
     
     std::vector< edm::Ref<l1extra::L1JetParticleCollection> > l1tauobjref;
     std::vector< edm::Ref<l1extra::L1JetParticleCollection> > l1jetobjref;
     std::vector< edm::Ref<l1extra::L1JetParticleCollection> > l1forjetobjref;
     
     l1trigobj->getObjects(trigger::TriggerL1TauJet, l1tauobjref);
     l1trigobj->getObjects(trigger::TriggerL1CenJet, l1jetobjref);
     l1trigobj->getObjects(trigger::TriggerL1ForJet, l1forjetobjref);
     
     double ptTriggered  = -10;
     double etaTriggered = -100;
     double phiTriggered = -100;
     for (unsigned int p=0; p<l1tauobjref.size(); p++) {
       if (l1tauobjref[p]->pt()>ptTriggered) {
     ptTriggered  = l1tauobjref[p]->pt(); 
     phiTriggered = l1tauobjref[p]->phi();
     etaTriggered = l1tauobjref[p]->eta();
       }
     }
     for (unsigned int p=0; p<l1jetobjref.size(); p++) {
       if (l1jetobjref[p]->pt()>ptTriggered) {
     ptTriggered  = l1jetobjref[p]->pt();
     phiTriggered = l1jetobjref[p]->phi();
     etaTriggered = l1jetobjref[p]->eta();
       }
     }
     for (unsigned int p=0; p<l1forjetobjref.size(); p++) {
       if (l1forjetobjref[p]->pt()>ptTriggered) {
     ptTriggered=l1forjetobjref[p]->pt();
     phiTriggered=l1forjetobjref[p]->phi();
     etaTriggered=l1forjetobjref[p]->eta();
       }
     }
     for(unsigned iS=0; iS<pixelTrackRefsHE.size(); iS++) {
       reco::VertexCollection::const_iterator vitSel;
       double minDZ = 100;
       bool vtxMatch;
       for (reco::VertexCollection::const_iterator vit=pVertHE->begin(); vit!=pVertHE->end(); vit++) {
     if (fabs(pixelTrackRefsHE[iS]->dz(vit->position()))<minDZ) {
       minDZ  = fabs(pixelTrackRefsHE[iS]->dz(vit->position()));
       vitSel = vit;
     }
       }
       //cut on dYX:
       if (minDZ!=100&&fabs(pixelTrackRefsHE[iS]->dxy(vitSel->position()))<vtxCutSeed_) vtxMatch=true;
       if (minDZ==100) vtxMatch=true;
       
       //select tracks not matched to triggered L1 jet
       double R=deltaR(etaTriggered, phiTriggered, pixelTrackRefsHE[iS]->eta(), pixelTrackRefsHE[iS]->phi());
       if (R>tauUnbiasCone_ && vtxMatch) nSeedHE++;
     }
     for(unsigned iS=0; iS<pixelTrackRefsHB.size(); iS++) {
       reco::VertexCollection::const_iterator vitSel;
       double minDZ = 100;
       bool vtxMatch;
       for (reco::VertexCollection::const_iterator vit=pVertHB->begin(); vit!=pVertHB->end(); vit++) {
     if (fabs(pixelTrackRefsHB[iS]->dz(vit->position()))<minDZ) {
       minDZ  = fabs(pixelTrackRefsHB[iS]->dz(vit->position()));
       vitSel = vit;
     }
       }
       //cut on dYX:
       if (minDZ!=100&&fabs(pixelTrackRefsHB[iS]->dxy(vitSel->position()))<101.0) vtxMatch=true;
       if (minDZ==100) vtxMatch=true;
       
       //select tracks not matched to triggered L1 jet
       double R=deltaR(etaTriggered, phiTriggered, pixelTrackRefsHB[iS]->eta(), pixelTrackRefsHB[iS]->phi());
       if (R>1.2 && vtxMatch) nSeedHB++;
     }
     
     edm::Service<FastTimerService> ft;
     tHE = ft->queryModuleTimeByLabel(theEvent.streamID(),"hltIsolPixelTrackProdHE") ;
     tHB = ft->queryModuleTimeByLabel(theEvent.streamID(),"hltIsolPixelTrackProdHB");
     
     if (verbosity%10>0) std::cout << "(HB/HE) time: " << tHB <<"/" << tHE << 
               " nCand: " << nCandHB << "/" << nCandHE << 
               "nSeed: " << nSeedHB << "/" << nSeedHE << std::endl;
   }
   t_timeL2Prod->push_back(tHB);  t_timeL2Prod->push_back(tHE);
   t_nPixSeed->push_back(nSeedHB); t_nPixSeed->push_back(nSeedHE);
   t_nPixCand->push_back(nCandHB); t_nPixCand->push_back(nCandHE);
 
   TimingTree->Fill();
 }

void IsoTrig::studyTrigger	(	edm::Handle< reco::TrackCollection > &	trkCollection,
		std::vector< reco::TrackCollection::const_iterator > &	goodTks
	)

private

Definition at line 1074 of file IsoTrig.cc.

References funct::abs(), gather_cfg::cout, deltaR(), dEta(), dPhi(), dR(), HcalObjRepresent::Fill(), fillDifferences(), fillHist(), h_dEtaL1, h_dPhiL1, h_dRL1, i, edm::HandleBase::isValid(), j, roll_playback::k, mathSSE::sqrt(), vec, and verbosity.

Referenced by analyze().

                                                                       {
 
   if (verbosity%10 > 0) std::cout << "Inside StudyTrigger" << std::endl;
   for (int j=0; j<3; j++) {
     for (unsigned int k=0; k<vec[j].size(); k++) {
       if (verbosity%10 > 0) std::cout << "vec[" << j << "][" << k << "] pt " << vec[j][k].pt() << " eta " << vec[j][k].eta() << " phi " << vec[j][k].phi() << std::endl;
       fillHist(j, vec[j][k]);
     }
   }
       
   double deta, dphi, dr;
   for (int lvl=1; lvl<3; lvl++) {
     for (unsigned int i=0; i<vec[lvl].size(); i++) {
       deta = dEta(vec[0][0],vec[lvl][i]);
       dphi = dPhi(vec[0][0],vec[lvl][i]);
       dr   = dR(vec[0][0],vec[lvl][i]);
       if (verbosity%10 > 1) std::cout << "lvl " <<lvl << " i " << i << " deta " << deta << " dphi " << dphi << " dR " << dr << std::endl;
       h_dEtaL1[lvl-1] -> Fill(deta);
       h_dPhiL1[lvl-1] -> Fill(dphi);
       h_dRL1[lvl-1]   -> Fill(std::sqrt(dr));
     }
   }
 
   math::XYZTLorentzVector mindRvec;
   double mindR;
   for (unsigned int k=0; k<vec[2].size(); ++k) {
     mindR=999.9;
     if (verbosity%10 > 1) std::cout << "L3obj: pt " << vec[2][k].pt() << " eta " << vec[2][k].eta() << " phi " << vec[2][k].phi() << std::endl;
     for (unsigned int j=0; j<vec[1].size(); j++) {
       dr   = dR(vec[2][k],vec[1][j]);
       if (dr<mindR) {
     mindR=dr;
     mindRvec=vec[1][j];
       }
     }
     fillDifferences(0, vec[2][k], mindRvec, (verbosity%10 >0));
     if (mindR < 0.03) {
       fillDifferences(1, vec[2][k], mindRvec, (verbosity%10 >0));
       fillHist(6, mindRvec);
       fillHist(8, vec[2][k]);
     } else {
       fillDifferences(2, vec[2][k], mindRvec, (verbosity%10 >0));
       fillHist(7, mindRvec);
       fillHist(9, vec[2][k]);
     }
                   
     mindR=999.9;
     if (verbosity%10 > 0) 
       std::cout << "Now Matching L3 track with reco: L3 Track (eta, phi) " 
         << vec[2][k].eta() << ":" << vec[2][k].phi() << " L2 Track "
         << mindRvec.eta() << ":" << mindRvec.phi() << " dR "
         << mindR << std::endl;
     reco::TrackCollection::const_iterator goodTk = trkCollection->end();
     if (trkCollection.isValid()) {
       double mindP(9999.9), mindP1(9999.9);
       reco::TrackCollection::const_iterator trkItr;
       for (trkItr=trkCollection->begin(); 
        trkItr!=trkCollection->end(); trkItr++) {
     math::XYZTLorentzVector v4(trkItr->px(), trkItr->py(), 
                    trkItr->pz(), trkItr->p());
     double deltaR = dR(v4, vec[2][k]);
     double dp     = std::abs(v4.r()/vec[2][k].r()-1.0);
     double dp1    = std::abs(1./v4.r()-1./vec[2][k].r());
     if (deltaR<mindR) {
       mindR    = deltaR;
       mindP    = dp;
       mindP1   = dp1;
       mindRvec = v4;
       goodTk   = trkItr;
     }
     if ((verbosity/10)%10>1 && deltaR<1.0) {
       std::cout << "reco track: pt " << v4.pt() << " eta " << v4.eta()
             << " phi " << v4.phi() << " DR " << deltaR 
             << std::endl;
     }
       }
       if (verbosity%10 > 0) 
     std::cout << "Now Matching at Reco level in step 1 DR: "  << mindR
           << ":" << mindP << ":" << mindP1 << " eta:phi " 
           << mindRvec.eta() << ":" << mindRvec.phi() << std::endl;
       if (mindR < 0.03 && mindP > 0.1) {
     for (trkItr=trkCollection->begin(); 
          trkItr!=trkCollection->end(); trkItr++) {
       math::XYZTLorentzVector v4(trkItr->px(), trkItr->py(), 
                      trkItr->pz(), trkItr->p());
       double deltaR = dR(v4, vec[2][k]);
       double dp     = std::abs(v4.r()/vec[2][k].r()-1.0);
       double dp1    = std::abs(1./v4.r()-1./vec[2][k].r());
       if (dp<mindP && deltaR<0.03) {
         mindR    = deltaR;
         mindP    = dp;
         mindP1   = dp1;
         mindRvec = v4;
         goodTk   = trkItr;
       }
     }
     if (verbosity%10 > 0) 
       std::cout << "Now Matching at Reco level in step 2 DR: "  << mindR
             << ":" << mindP << ":" << mindP1 << " eta:phi " 
             << mindRvec.eta() << ":" << mindRvec.phi() << std::endl;
 
       }
       fillDifferences(3, vec[2][k], mindRvec, (verbosity%10 >0));
       fillHist(3, mindRvec);
       if(mindR < 0.03) {
     fillDifferences(4, vec[2][k], mindRvec, (verbosity%10 >0));
     fillHist(4, mindRvec);
       } else {
     fillDifferences(5, vec[2][k], mindRvec, (verbosity%10 >0));
     fillHist(5, mindRvec);
       }
       if (goodTk != trkCollection->end()) goodTks.push_back(goodTk);
     }
   }
 }

void IsoTrig::StudyTrkEbyP ( edm::Handle< reco::TrackCollection > & trkCollection )

private

Definition at line 747 of file IsoTrig.cc.

References a_charIsoR, a_coneR, a_mipR, a_neutR1, a_neutR2, barrelRecHitsHandle, ecalTB2006H4_GenSimDigiReco_cfg::bField, bFieldH, spr::chargeIsolationCone(), gather_cfg::cout, cutCharge, cutMip, cutNeutral, spr::eCone_ecal(), spr::eCone_hcal(), endcapRecHitsHandle, spr::goodTrack(), hbhe, edm::HandleBase::isValid(), leadPV, spr::trackSelectionParameters::maxDxyPV, spr::trackSelectionParameters::maxDzPV, spr::trackSelectionParameters::maxInMiss, spr::trackSelectionParameters::maxOutMiss, reco::TrackBase::p(), pG, edm::ESHandle< class >::product(), spr::propagateCALO(), reco::TrackBase::px(), reco::TrackBase::py(), reco::TrackBase::pz(), selectionParameters, t_TrkhCone, t_TrkMissFlag, t_TrkNuIsolFlag, t_TrkP, t_TrkPVFlag, t_TrkqltyFlag, t_TrkselTkFlag, theTrackQuality, TrkResTree, and verbosity.

Referenced by analyze().

                                                                         {
 
   t_TrkselTkFlag->clear();
   t_TrkqltyFlag->clear();
   t_TrkMissFlag->clear();
   t_TrkPVFlag->clear();
   t_TrkNuIsolFlag->clear();
   t_TrkhCone->clear();
   t_TrkP->clear();
   
   if (!trkCollection.isValid()) {
     std::cout << "trkCollection.isValid is false" << std::endl;
   } else {
     std::vector<spr::propagatedTrackDirection>::const_iterator trkDetItr; 
     const MagneticField *bField = bFieldH.product();
     const CaloGeometry* geo = pG.product();
     std::vector<spr::propagatedTrackDirection> trkCaloDirections1;
     spr::propagateCALO(trkCollection, geo, bField, theTrackQuality, trkCaloDirections1, ((verbosity/100)%10>2));
     unsigned int nTracks=0;
     int nRH_eMipDR=0, nNearTRKs=0;
     std::vector<bool> selFlags;
     for (trkDetItr = trkCaloDirections1.begin(); trkDetItr != trkCaloDirections1.end(); trkDetItr++,nTracks++) {
       double conehmaxNearP = 0, hCone=0, eMipDR=0.0;
       const reco::Track* pTrack = &(*(trkDetItr->trkItr));
       if (verbosity%10>0) std::cout << "track no. " << nTracks << " p(): " << pTrack->p() << std::endl;
       if (pTrack->p() > 20) {
     math::XYZTLorentzVector v2(pTrack->px(), pTrack->py(),  
                    pTrack->pz(), pTrack->p());
     eMipDR = spr::eCone_ecal(geo, barrelRecHitsHandle, endcapRecHitsHandle,
                  trkDetItr->pointHCAL, trkDetItr->pointECAL,
                  a_mipR, trkDetItr->directionECAL, nRH_eMipDR);
     //    std::cout << "mip " << eMipDR << std::endl;
     bool selectTk = spr::goodTrack(pTrack,leadPV,selectionParameters,((verbosity/100)%10>1));
     spr::trackSelectionParameters oneCutParameters = selectionParameters;
     oneCutParameters.maxDxyPV  = 10;
     oneCutParameters.maxDzPV   = 100;
     oneCutParameters.maxInMiss = 2;
     oneCutParameters.maxOutMiss= 2;
     bool qltyFlag     =  spr::goodTrack(pTrack,leadPV,oneCutParameters,((verbosity/100)%10>1));
     oneCutParameters           = selectionParameters;
     oneCutParameters.maxDxyPV  = 10;
     oneCutParameters.maxDzPV   = 100;
     bool qltyMissFlag =  spr::goodTrack(pTrack,leadPV,oneCutParameters,((verbosity/100)%10>1));
     oneCutParameters           = selectionParameters;
     oneCutParameters.maxInMiss = 2;
     oneCutParameters.maxOutMiss= 2;
     bool qltyPVFlag   =  spr::goodTrack(pTrack,leadPV,oneCutParameters,((verbosity/100)%10>1));
     /*
     std::cout << "sel " << selectTk << std::endl;
     std::cout << "ntracks " << nTracks;
     std::cout << " a_charIsoR " << a_charIsoR;
     std::cout << " nNearTRKs " << nNearTRKs << std::endl; 
     */
     conehmaxNearP = spr::chargeIsolationCone(nTracks, trkCaloDirections1, a_charIsoR, nNearTRKs, ((verbosity/100)%10>1));
     /*
     std::cout << "coneh " << conehmaxNearP << std::endl;
     std::cout << "ok " << trkDetItr->okECAL << " " << trkDetItr->okHCAL << std::endl;
     */
     double e1 =  spr::eCone_ecal(geo, barrelRecHitsHandle, endcapRecHitsHandle,
                      trkDetItr->pointHCAL, trkDetItr->pointECAL,
                      a_neutR1, trkDetItr->directionECAL, nRH_eMipDR);
 //  std::cout << "e1 " << e1 << std::endl;
     double e2 =  spr::eCone_ecal(geo, barrelRecHitsHandle, endcapRecHitsHandle,
                      trkDetItr->pointHCAL, trkDetItr->pointECAL,
                      a_neutR2, trkDetItr->directionECAL, nRH_eMipDR);
     double e_inCone  = e2 - e1;
     bool chgIsolFlag = (conehmaxNearP < cutCharge);
     bool mipFlag     = (eMipDR < cutMip);
     bool neuIsolFlag = (e_inCone < cutNeutral);
     bool trkpropFlag = ((trkDetItr->okECAL) && (trkDetItr->okHCAL));
     selFlags.clear();
     selFlags.push_back(selectTk);        selFlags.push_back(qltyFlag);
     selFlags.push_back(qltyMissFlag);    selFlags.push_back(qltyPVFlag);
     if (verbosity%10>0) std::cout << "emip: " << eMipDR << "<" << cutMip << "(" << mipFlag << ")" 
                       << " ; ok: " << trkDetItr->okECAL << "/" << trkDetItr->okHCAL 
                       << " ; chgiso: " << conehmaxNearP << "<" << cutCharge << "(" << chgIsolFlag << ")" << std::endl;
     //  std::cout << selectTk << " " << qltyFlag << " " << qltyMissFlag << " " << qltyPVFlag << " " << trkpropFlag << " " << chgIsolFlag << " " << neuIsolFlag << " " << mipFlag << std::endl;
     
     if(chgIsolFlag && mipFlag && trkpropFlag) {
       double             distFromHotCell=-99.0;
       int                nRecHitsCone=-99, ietaHotCell=-99, iphiHotCell=-99;
       GlobalPoint        gposHotCell(0.,0.,0.);
       std::vector<DetId> coneRecHitDetIds;
       hCone    = spr::eCone_hcal(geo, hbhe, trkDetItr->pointHCAL, 
                      trkDetItr->pointECAL,
                      a_coneR, trkDetItr->directionHCAL, 
                      nRecHitsCone, coneRecHitDetIds,
                      distFromHotCell, ietaHotCell, iphiHotCell,
                      gposHotCell, -1);
       // push vectors into the Tree
       t_TrkselTkFlag ->push_back(selFlags[0]);
       t_TrkqltyFlag  ->push_back(selFlags[1]);
       t_TrkMissFlag  ->push_back(selFlags[2]);
       t_TrkPVFlag    ->push_back(selFlags[3]);
       t_TrkNuIsolFlag->push_back(neuIsolFlag);
       t_TrkhCone     ->push_back(hCone);
       t_TrkP         ->push_back(pTrack->p());
     }
       }
     }
     if (verbosity%10>0) std::cout << "Filling " << t_TrkP->size() << " tracks in TrkRestree out of " << nTracks << std::endl;
   }
   TrkResTree->Fill();
 }

Member Data Documentation

double IsoTrig::a_charIsoR

private

Definition at line 133 of file IsoTrig.h.

Referenced by chgIsolation(), IsoTrig(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

double IsoTrig::a_coneR

private

Definition at line 133 of file IsoTrig.h.

Referenced by IsoTrig(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

double IsoTrig::a_mipR

private

Definition at line 133 of file IsoTrig.h.

Referenced by IsoTrig(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

double IsoTrig::a_neutIsoR

private

Definition at line 133 of file IsoTrig.h.

Referenced by IsoTrig().

double IsoTrig::a_neutR1

private

Definition at line 134 of file IsoTrig.h.

Referenced by IsoTrig(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

double IsoTrig::a_neutR2

private

Definition at line 134 of file IsoTrig.h.

Referenced by IsoTrig(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

edm::Handle<EcalRecHitCollection> IsoTrig::barrelRecHitsHandle

private

Definition at line 159 of file IsoTrig.h.

Referenced by analyze(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

edm::Handle<reco::BeamSpot> IsoTrig::beamSpotH

private

Definition at line 161 of file IsoTrig.h.

Referenced by analyze().

edm::ESHandle<MagneticField> IsoTrig::bFieldH

private

Definition at line 156 of file IsoTrig.h.

Referenced by analyze(), chgIsolation(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

double IsoTrig::bfVal

private

Definition at line 127 of file IsoTrig.h.

Referenced by analyze(), and GetEtaPhiAtEcal().

bool IsoTrig::changed

private

Definition at line 167 of file IsoTrig.h.

Referenced by analyze(), and beginRun().

TTree * IsoTrig::ChgIsolnTree

private

Definition at line 169 of file IsoTrig.h.

Referenced by beginJob(), and chgIsolation().

double IsoTrig::cutCharge

private

Definition at line 134 of file IsoTrig.h.

Referenced by fillCuts(), IsoTrig(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

double IsoTrig::cutMip

private

Definition at line 134 of file IsoTrig.h.

Referenced by fillCuts(), IsoTrig(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

double IsoTrig::cutNeutral

private

Definition at line 134 of file IsoTrig.h.

Referenced by fillCuts(), IsoTrig(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

bool IsoTrig::doChgIsolTree

private

Definition at line 124 of file IsoTrig.h.

Referenced by analyze(), beginJob(), and IsoTrig().

bool IsoTrig::doL2L3

private

Definition at line 124 of file IsoTrig.h.

Referenced by analyze(), beginJob(), IsoTrig(), and studyIsolation().

bool IsoTrig::doMipCutTree

private

Definition at line 124 of file IsoTrig.h.

Referenced by analyze(), beginJob(), and IsoTrig().

bool IsoTrig::doStudyIsol

private

Definition at line 124 of file IsoTrig.h.

Referenced by analyze(), and IsoTrig().

bool IsoTrig::doTiming

private

Definition at line 124 of file IsoTrig.h.

Referenced by analyze(), beginJob(), and IsoTrig().

bool IsoTrig::doTrkResTree

private

Definition at line 124 of file IsoTrig.h.

Referenced by analyze(), beginJob(), and IsoTrig().

double IsoTrig::dr_L1

private

Definition at line 133 of file IsoTrig.h.

Referenced by IsoTrig(), and studyIsolation().

edm::Handle<EcalRecHitCollection> IsoTrig::endcapRecHitsHandle

private

Definition at line 160 of file IsoTrig.h.

Referenced by analyze(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

edm::Service<TFileService> IsoTrig::fs

private

Definition at line 168 of file IsoTrig.h.

Referenced by beginJob(), and endJob().

TH1I * IsoTrig::g_Accepts

private

Definition at line 228 of file IsoTrig.h.

Referenced by endJob().

TH1I* IsoTrig::g_Pre

private

Definition at line 228 of file IsoTrig.h.

Referenced by endJob().

TH1I * IsoTrig::g_PreHLT

private

Definition at line 228 of file IsoTrig.h.

Referenced by endJob().

TH1I * IsoTrig::g_PreL1

private

Definition at line 228 of file IsoTrig.h.

Referenced by endJob().

TH1D* IsoTrig::h_dEta[9]

private

Definition at line 223 of file IsoTrig.h.

Referenced by beginJob(), and fillDifferences().

TH1D* IsoTrig::h_dEtaL1[2]

private

Definition at line 222 of file IsoTrig.h.

Referenced by beginJob(), and studyTrigger().

TH1D* IsoTrig::h_dinvPt[9]

private

Definition at line 224 of file IsoTrig.h.

Referenced by beginJob(), and fillDifferences().

TH1D * IsoTrig::h_dP[9]

private

Definition at line 223 of file IsoTrig.h.

Referenced by beginJob(), and fillDifferences().

TH1D * IsoTrig::h_dPhi[9]

private

Definition at line 223 of file IsoTrig.h.

Referenced by beginJob(), and fillDifferences().

TH1D * IsoTrig::h_dPhiL1[2]

private

Definition at line 222 of file IsoTrig.h.

Referenced by beginJob(), and studyTrigger().

TH1D * IsoTrig::h_dPt[9]

private

Definition at line 223 of file IsoTrig.h.

Referenced by beginJob(), and fillDifferences().

TH1D * IsoTrig::h_dRL1[2]

private

Definition at line 222 of file IsoTrig.h.

Referenced by beginJob(), and studyTrigger().

TH1D * IsoTrig::h_eCalo[5][6][48]

private

Definition at line 227 of file IsoTrig.h.

Referenced by beginJob(), and fillEnergy().

TH1D* IsoTrig::h_eHcal[5][6][48]

private

Definition at line 227 of file IsoTrig.h.

Referenced by beginJob(), and fillEnergy().

TH1D* IsoTrig::h_eMaxNearP[2]

private

Definition at line 225 of file IsoTrig.h.

Referenced by beginJob(), and fillCuts().

TH1D * IsoTrig::h_eMip[2]

private

Definition at line 224 of file IsoTrig.h.

Referenced by beginJob(), and fillCuts().

TH1D * IsoTrig::h_eNeutIso[2]

private

Definition at line 225 of file IsoTrig.h.

Referenced by beginJob(), and fillCuts().

TH1D* IsoTrig::h_EnIn

private

Definition at line 216 of file IsoTrig.h.

Referenced by beginJob(), and studyMipCut().

TH1D * IsoTrig::h_EnOut

private

Definition at line 216 of file IsoTrig.h.

Referenced by beginJob(), and studyMipCut().

TH1D * IsoTrig::h_eta[20]

private

Definition at line 221 of file IsoTrig.h.

Referenced by beginJob(), and fillHist().

TH1D* IsoTrig::h_etaCalibTracks[5][2][2]

private

Definition at line 226 of file IsoTrig.h.

Referenced by beginJob(), and fillCuts().

TH1D * IsoTrig::h_etaMipTracks[5][2][2]

private

Definition at line 226 of file IsoTrig.h.

Referenced by beginJob(), and fillCuts().

TH1I * IsoTrig::h_Filters

private

Definition at line 218 of file IsoTrig.h.

Referenced by analyze(), and beginJob().

TH1I * IsoTrig::h_HLT

private

Definition at line 218 of file IsoTrig.h.

Referenced by analyze(), and beginJob().

TH1D * IsoTrig::h_L1ObjEnergy

private

Definition at line 220 of file IsoTrig.h.

Referenced by analyze(), and beginJob().

TH1D * IsoTrig::h_mindR[9]

private

Definition at line 224 of file IsoTrig.h.

Referenced by beginJob(), and fillDifferences().

TH2D* IsoTrig::h_MipEnMatch

private

Definition at line 217 of file IsoTrig.h.

Referenced by beginJob(), and studyMipCut().

TH2D * IsoTrig::h_MipEnNoMatch

private

Definition at line 217 of file IsoTrig.h.

Referenced by beginJob(), and studyMipCut().

TH1I* IsoTrig::h_nHLT

private

Definition at line 218 of file IsoTrig.h.

Referenced by analyze(), and beginJob().

TH1I * IsoTrig::h_nL3Objs

private

Definition at line 218 of file IsoTrig.h.

Referenced by analyze(), and beginJob().

TH1D* IsoTrig::h_p[20]

private

Definition at line 221 of file IsoTrig.h.

Referenced by beginJob(), and fillHist().

TH1D * IsoTrig::h_phi[20]

private

Definition at line 221 of file IsoTrig.h.

Referenced by beginJob(), and fillHist().

TH1I * IsoTrig::h_Pre

private

Definition at line 218 of file IsoTrig.h.

Referenced by analyze(), and beginJob().

TH1I * IsoTrig::h_PreHLT

private

Definition at line 218 of file IsoTrig.h.

Referenced by analyze(), and beginJob().

TH1D * IsoTrig::h_PreHLTwt

private

Definition at line 220 of file IsoTrig.h.

Referenced by analyze(), and beginJob().

TH1I * IsoTrig::h_PreL1

private

Definition at line 218 of file IsoTrig.h.

Referenced by analyze(), and beginJob().

TH1D* IsoTrig::h_PreL1wt

private

Definition at line 220 of file IsoTrig.h.

Referenced by analyze(), and beginJob().

TH1D * IsoTrig::h_pt[20]

private

Definition at line 221 of file IsoTrig.h.

Referenced by beginJob(), and fillHist().

edm::Handle<HBHERecHitCollection> IsoTrig::hbhe

private

Definition at line 158 of file IsoTrig.h.

Referenced by analyze(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

HLTConfigProvider IsoTrig::hltConfig_

private

Definition at line 120 of file IsoTrig.h.

Referenced by analyze(), and beginRun().

edm::InputTag IsoTrig::L1candTag_

private

Definition at line 122 of file IsoTrig.h.

Referenced by IsoTrig().

edm::InputTag IsoTrig::L2candTag_

private

Definition at line 122 of file IsoTrig.h.

Referenced by IsoTrig().

math::XYZPoint IsoTrig::leadPV

private

Definition at line 163 of file IsoTrig.h.

Referenced by analyze(), chgIsolation(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

int IsoTrig::maxRunNo

private

Definition at line 135 of file IsoTrig.h.

Referenced by endJob(), and IsoTrig().

double IsoTrig::minPTrackValue_

private

Definition at line 129 of file IsoTrig.h.

Referenced by chgIsolation(), and IsoTrig().

int IsoTrig::minRunNo

private

Definition at line 135 of file IsoTrig.h.

Referenced by endJob(), and IsoTrig().

TTree* IsoTrig::MipCutTree

private

Definition at line 169 of file IsoTrig.h.

Referenced by beginJob(), and studyMipCut().

edm::ESHandle<CaloGeometry> IsoTrig::pG

private

Definition at line 157 of file IsoTrig.h.

Referenced by analyze(), chgIsolation(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

edm::InputTag IsoTrig::PixcandTag_

private

Definition at line 122 of file IsoTrig.h.

Referenced by IsoTrig().

std::vector<double> IsoTrig::pixelIsolationConeSizeAtEC_

private

Definition at line 128 of file IsoTrig.h.

Referenced by chgIsolation(), and IsoTrig().

std::vector<reco::TrackRef> IsoTrig::pixelTrackRefsHB

private

Definition at line 155 of file IsoTrig.h.

Referenced by analyze(), and studyTiming().

std::vector<reco::TrackRef> IsoTrig::pixelTrackRefsHE

private

Definition at line 155 of file IsoTrig.h.

Referenced by analyze(), chgIsolation(), and studyTiming().

std::vector<edm::InputTag> IsoTrig::pixelTracksSources_

private

Definition at line 123 of file IsoTrig.h.

Referenced by analyze(), and IsoTrig().

double IsoTrig::pLimits[6]

private

Definition at line 136 of file IsoTrig.h.

Referenced by beginJob(), fillCuts(), fillEnergy(), and IsoTrig().

double IsoTrig::prelimCone_

private

Definition at line 130 of file IsoTrig.h.

Referenced by chgIsolation(), and IsoTrig().

std::string IsoTrig::processName

private

Definition at line 132 of file IsoTrig.h.

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

double IsoTrig::rEB_

private

Definition at line 127 of file IsoTrig.h.

Referenced by getDistInCM(), GetEtaPhiAtEcal(), and IsoTrig().

edm::Handle<reco::VertexCollection> IsoTrig::recVtxs

private

Definition at line 162 of file IsoTrig.h.

Referenced by analyze(), and chgIsolation().

spr::trackSelectionParameters IsoTrig::selectionParameters

private

Definition at line 131 of file IsoTrig.h.

Referenced by chgIsolation(), IsoTrig(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

std::vector<double>* IsoTrig::t_ECone

private

Definition at line 214 of file IsoTrig.h.

Referenced by beginJob(), clearMipCutTreeVectors(), and pushMipCutTreeVecs().

std::vector<double>* IsoTrig::t_NFcandEmip

private

Definition at line 198 of file IsoTrig.h.

Referenced by beginJob(), clearMipCutTreeVectors(), and pushMipCutTreeVecs().

std::vector<double>* IsoTrig::t_NFcandEta

private

Definition at line 196 of file IsoTrig.h.

Referenced by beginJob(), clearMipCutTreeVectors(), and pushMipCutTreeVecs().

std::vector<double>* IsoTrig::t_NFcandP

private

Definition at line 194 of file IsoTrig.h.

Referenced by beginJob(), clearMipCutTreeVectors(), and pushMipCutTreeVecs().

std::vector<double>* IsoTrig::t_NFcandPhi

private

Definition at line 197 of file IsoTrig.h.

Referenced by beginJob(), clearMipCutTreeVectors(), and pushMipCutTreeVecs().

std::vector<double>* IsoTrig::t_NFcandPt

private

Definition at line 195 of file IsoTrig.h.

Referenced by beginJob(), clearMipCutTreeVectors(), and pushMipCutTreeVecs().

std::vector<double>* IsoTrig::t_NFTrkcandEmip

private

Definition at line 203 of file IsoTrig.h.

Referenced by beginJob(), clearMipCutTreeVectors(), and pushMipCutTreeVecs().

std::vector<double>* IsoTrig::t_NFTrkcandEta

private

Definition at line 201 of file IsoTrig.h.

Referenced by beginJob(), clearMipCutTreeVectors(), and pushMipCutTreeVecs().

std::vector<double>* IsoTrig::t_NFTrkcandP

private

Definition at line 199 of file IsoTrig.h.

Referenced by beginJob(), clearMipCutTreeVectors(), and pushMipCutTreeVecs().

std::vector<double>* IsoTrig::t_NFTrkcandPhi

private

Definition at line 202 of file IsoTrig.h.

Referenced by beginJob(), clearMipCutTreeVectors(), and pushMipCutTreeVecs().

std::vector<double>* IsoTrig::t_NFTrkcandPt

private

Definition at line 200 of file IsoTrig.h.

Referenced by beginJob(), clearMipCutTreeVectors(), and pushMipCutTreeVecs().

std::vector<bool>* IsoTrig::t_NFTrkChgIsoFlag

private

Definition at line 211 of file IsoTrig.h.

Referenced by beginJob(), clearMipCutTreeVectors(), and pushMipCutTreeVecs().

std::vector<double>* IsoTrig::t_NFTrkMinDP1

private

Definition at line 205 of file IsoTrig.h.

Referenced by beginJob(), clearMipCutTreeVectors(), and pushMipCutTreeVecs().

std::vector<double>* IsoTrig::t_NFTrkMinDR

private

Definition at line 204 of file IsoTrig.h.

Referenced by beginJob(), clearMipCutTreeVectors(), and pushMipCutTreeVecs().

std::vector<bool>* IsoTrig::t_NFTrkMipFlag

private

Definition at line 213 of file IsoTrig.h.

Referenced by beginJob(), clearMipCutTreeVectors(), and pushMipCutTreeVecs().

std::vector<bool>* IsoTrig::t_NFTrkMissFlag

private

Definition at line 208 of file IsoTrig.h.

Referenced by beginJob(), clearMipCutTreeVectors(), and pushMipCutTreeVecs().

std::vector<bool>* IsoTrig::t_NFTrkNeuIsoFlag

private

Definition at line 212 of file IsoTrig.h.

Referenced by beginJob(), clearMipCutTreeVectors(), and pushMipCutTreeVecs().

std::vector<bool>* IsoTrig::t_NFTrkPropFlag

private

Definition at line 210 of file IsoTrig.h.

Referenced by beginJob(), clearMipCutTreeVectors(), and pushMipCutTreeVecs().

std::vector<bool>* IsoTrig::t_NFTrkPVFlag

private

Definition at line 209 of file IsoTrig.h.

Referenced by beginJob(), clearMipCutTreeVectors(), and pushMipCutTreeVecs().

std::vector<bool>* IsoTrig::t_NFTrkqltyFlag

private

Definition at line 207 of file IsoTrig.h.

Referenced by beginJob(), clearMipCutTreeVectors(), and pushMipCutTreeVecs().

std::vector<bool>* IsoTrig::t_NFTrkselTkFlag

private

Definition at line 206 of file IsoTrig.h.

Referenced by beginJob(), clearMipCutTreeVectors(), and pushMipCutTreeVecs().

std::vector<int>* IsoTrig::t_nPixCand

private

Definition at line 171 of file IsoTrig.h.

Referenced by beginJob(), and studyTiming().

std::vector<int>* IsoTrig::t_nPixSeed

private

Definition at line 172 of file IsoTrig.h.

Referenced by beginJob(), and studyTiming().

std::vector<double>* IsoTrig::t_PixcandEta

private

Definition at line 184 of file IsoTrig.h.

Referenced by beginJob(), clearChgIsolnTreeVectors(), and pushChgIsolnTreeVecs().

std::vector<std::vector<double> >* IsoTrig::t_PixcandMaxP

private

Definition at line 186 of file IsoTrig.h.

Referenced by beginJob(), clearChgIsolnTreeVectors(), and pushChgIsolnTreeVecs().

std::vector<double>* IsoTrig::t_PixcandP

private

Definition at line 182 of file IsoTrig.h.

Referenced by beginJob(), clearChgIsolnTreeVectors(), and pushChgIsolnTreeVecs().

std::vector<double>* IsoTrig::t_PixcandPhi

private

Definition at line 185 of file IsoTrig.h.

Referenced by beginJob(), clearChgIsolnTreeVectors(), and pushChgIsolnTreeVecs().

std::vector<double>* IsoTrig::t_PixcandPt

private

Definition at line 183 of file IsoTrig.h.

Referenced by beginJob(), clearChgIsolnTreeVectors(), and pushChgIsolnTreeVecs().

std::vector<double>* IsoTrig::t_PixTrkcandEta

private

Definition at line 189 of file IsoTrig.h.

Referenced by beginJob(), clearChgIsolnTreeVectors(), and pushChgIsolnTreeVecs().

std::vector<double>* IsoTrig::t_PixTrkcandMaxP

private

Definition at line 191 of file IsoTrig.h.

Referenced by beginJob(), clearChgIsolnTreeVectors(), and pushChgIsolnTreeVecs().

std::vector<double>* IsoTrig::t_PixTrkcandP

private

Definition at line 187 of file IsoTrig.h.

Referenced by beginJob(), clearChgIsolnTreeVectors(), and pushChgIsolnTreeVecs().

std::vector<double>* IsoTrig::t_PixTrkcandPhi

private

Definition at line 190 of file IsoTrig.h.

Referenced by beginJob(), clearChgIsolnTreeVectors(), and pushChgIsolnTreeVecs().

std::vector<double>* IsoTrig::t_PixTrkcandPt

private

Definition at line 188 of file IsoTrig.h.

Referenced by beginJob(), clearChgIsolnTreeVectors(), and pushChgIsolnTreeVecs().

std::vector<bool>* IsoTrig::t_PixTrkcandselTk

private

Definition at line 192 of file IsoTrig.h.

Referenced by beginJob(), clearChgIsolnTreeVectors(), and pushChgIsolnTreeVecs().

std::vector<double>* IsoTrig::t_timeL2Prod

private

Definition at line 170 of file IsoTrig.h.

Referenced by beginJob(), and studyTiming().

std::vector<double>* IsoTrig::t_TrkhCone

private

Definition at line 174 of file IsoTrig.h.

Referenced by beginJob(), and StudyTrkEbyP().

std::vector<bool>* IsoTrig::t_TrkMissFlag

private

Definition at line 178 of file IsoTrig.h.

Referenced by beginJob(), and StudyTrkEbyP().

std::vector<bool>* IsoTrig::t_TrkNuIsolFlag

private

Definition at line 180 of file IsoTrig.h.

Referenced by beginJob(), and StudyTrkEbyP().

std::vector<double>* IsoTrig::t_TrkP

private

Definition at line 175 of file IsoTrig.h.

Referenced by beginJob(), and StudyTrkEbyP().

std::vector<bool>* IsoTrig::t_TrkPVFlag

private

Definition at line 179 of file IsoTrig.h.

Referenced by beginJob(), and StudyTrkEbyP().

std::vector<bool>* IsoTrig::t_TrkqltyFlag

private

Definition at line 177 of file IsoTrig.h.

Referenced by beginJob(), and StudyTrkEbyP().

std::vector<bool>* IsoTrig::t_TrkselTkFlag

private

Definition at line 176 of file IsoTrig.h.

Referenced by beginJob(), and StudyTrkEbyP().

double IsoTrig::tauUnbiasCone_

private

Definition at line 130 of file IsoTrig.h.

Referenced by chgIsolation(), IsoTrig(), and studyTiming().

std::string IsoTrig::theTrackQuality

private

Definition at line 132 of file IsoTrig.h.

Referenced by chgIsolation(), IsoTrig(), studyIsolation(), studyMipCut(), and StudyTrkEbyP().

TTree * IsoTrig::TimingTree

private

Definition at line 169 of file IsoTrig.h.

Referenced by beginJob(), and studyTiming().

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

private

Definition at line 143 of file IsoTrig.h.

Referenced by analyze(), and IsoTrig().

edm::EDGetTokenT<EcalRecHitCollection> IsoTrig::tok_EB_

private

Definition at line 144 of file IsoTrig.h.

Referenced by analyze(), and IsoTrig().

edm::EDGetTokenT<EcalRecHitCollection> IsoTrig::tok_EE_

private

Definition at line 145 of file IsoTrig.h.

Referenced by analyze(), and IsoTrig().

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

private

Definition at line 141 of file IsoTrig.h.

Referenced by analyze(), and IsoTrig().

edm::EDGetTokenT<HBHERecHitCollection> IsoTrig::tok_hbhe_

private

Definition at line 146 of file IsoTrig.h.

Referenced by analyze(), and IsoTrig().

edm::EDGetTokenT<trigger::TriggerFilterObjectWithRefs> IsoTrig::tok_hlt_

private

Definition at line 147 of file IsoTrig.h.

Referenced by IsoTrig(), and studyTiming().

edm::EDGetTokenT<trigger::TriggerFilterObjectWithRefs> IsoTrig::tok_l1cand_

private

Definition at line 151 of file IsoTrig.h.

Referenced by analyze(), and IsoTrig().

edm::EDGetTokenT<reco::IsolatedPixelTrackCandidateCollection> IsoTrig::tok_l2cand_

private

Definition at line 152 of file IsoTrig.h.

Referenced by analyze(), and IsoTrig().

edm::EDGetTokenT<LumiDetails> IsoTrig::tok_lumi

private

Definition at line 137 of file IsoTrig.h.

Referenced by analyze(), and IsoTrig().

edm::EDGetTokenT<reco::IsolatedPixelTrackCandidateCollection> IsoTrig::tok_pixtk_

private

Definition at line 150 of file IsoTrig.h.

Referenced by analyze(), and IsoTrig().

std::vector<edm::EDGetTokenT<reco::TrackCollection> > IsoTrig::tok_pixtks_

private

Definition at line 153 of file IsoTrig.h.

Referenced by analyze(), and IsoTrig().

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

private

Definition at line 142 of file IsoTrig.h.

Referenced by analyze(), and IsoTrig().

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

private

Definition at line 138 of file IsoTrig.h.

Referenced by analyze(), and IsoTrig().

edm::EDGetTokenT<edm::TriggerResults> IsoTrig::tok_trigRes

private

Definition at line 139 of file IsoTrig.h.

Referenced by analyze(), and IsoTrig().

edm::EDGetTokenT<reco::VertexCollection> IsoTrig::tok_verthb_

private

Definition at line 148 of file IsoTrig.h.

Referenced by IsoTrig(), and studyTiming().

edm::EDGetTokenT<reco::VertexCollection> IsoTrig::tok_verthe_

private

Definition at line 148 of file IsoTrig.h.

Referenced by IsoTrig(), and studyTiming().

std::map<unsigned int, unsigned int> IsoTrig::TrigList

private

Definition at line 165 of file IsoTrig.h.

Referenced by analyze(), and endJob().

std::vector<std::string> IsoTrig::trigNames

private

Definition at line 121 of file IsoTrig.h.

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

std::map<unsigned int, const std::pair<int, int> > IsoTrig::TrigPreList

private

Definition at line 166 of file IsoTrig.h.

Referenced by analyze(), and endJob().

TTree * IsoTrig::TrkResTree

private

Definition at line 169 of file IsoTrig.h.

Referenced by beginJob(), and StudyTrkEbyP().

std::vector<math::XYZTLorentzVector> IsoTrig::vec[3]

private

Definition at line 229 of file IsoTrig.h.

Referenced by analyze(), etaPhiTrigger(), studyIsolation(), and studyTrigger().

int IsoTrig::verbosity

private

Definition at line 126 of file IsoTrig.h.

Referenced by analyze(), chgIsolation(), endJob(), IsoTrig(), studyIsolation(), studyMipCut(), studyTiming(), studyTrigger(), and StudyTrkEbyP().

double IsoTrig::vtxCutIsol_

private

Definition at line 129 of file IsoTrig.h.

Referenced by chgIsolation(), and IsoTrig().

double IsoTrig::vtxCutSeed_

private

Definition at line 129 of file IsoTrig.h.

Referenced by chgIsolation(), IsoTrig(), and studyTiming().

double IsoTrig::zEE_

private

Definition at line 127 of file IsoTrig.h.

Referenced by getDistInCM(), GetEtaPhiAtEcal(), and IsoTrig().

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