Inheritance diagram for IsoTrackCalibration:

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

	~IsoTrackCalibration () override

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::SharedResources >
	EDAnalyzer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

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

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

void	beginJob () override

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

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

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

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

Private Attributes
double	a_charIsoR_

double	a_coneR_

double	a_mipR_

double	constTrackPt_

edm::Service< TFileService >	fs_

TH1I *	h_iEta

TH1F *	h_jetp

TH1F *	h_jetpt [4]

TH1F *	h_nTrk

TH1F *	h_nVtx

TH1F *	h_Rechit_E

TProfile *	h_RecHit_iEta

TProfile *	h_RecHit_num

TH1I *	h_tketa0 [5]

TH1I *	h_tketa1 [5]

TH1I *	h_tketa2 [5]

TH1I *	h_tketa3 [5]

TH1I *	h_tketa4 [5]

TH1I *	h_tketa5 [5]

TH1I *	h_tketav1 [5][6]

TH1I *	h_tketav2 [5][6]

HLTConfigProvider	hltConfig_

std::vector< std::string >	HLTNames_

bool	isAOD_

bool	isMC_

bool	isQCD_

double	maxEcalEnr_

double	maxNeighborTrackEnr_

int	nRun_

spr::trackSelectionParameters	selectionParameters_

double	slopeTrackPt_

std::vector< unsigned int > *	t_DetIds

std::vector< unsigned int > *	t_DetIds1

std::vector< unsigned int > *	t_DetIds3

double	t_eHcal

double	t_eHcal10

double	t_eHcal30

double	t_eMipDR

int	t_Event

double	t_EventWeight

std::vector< double > *	t_HitEnergies

std::vector< double > *	t_HitEnergies1

std::vector< double > *	t_HitEnergies3

double	t_hmaxNearP

int	t_ieta

int	t_nTrk

int	t_nVtx

double	t_p

double	t_phi

double	t_pt

bool	t_qltyMissFlag

bool	t_qltyPVFlag

int	t_Run

bool	t_selectTk

std::string	theTrackQuality_

edm::InputTag	theTriggerResultsLabel_

edm::EDGetTokenT< reco::BeamSpot >	tok_bs_

edm::EDGetTokenT< EcalRecHitCollection >	tok_EB_

edm::EDGetTokenT< EcalRecHitCollection >	tok_EE_

edm::EDGetTokenT< GenEventInfoProduct >	tok_ew_

edm::EDGetTokenT< reco::TrackCollection >	tok_genTrack_

edm::EDGetTokenT< HBHERecHitCollection >	tok_hbhe_

edm::EDGetTokenT< reco::VertexCollection >	tok_recVtx_

edm::EDGetTokenT< trigger::TriggerEvent >	tok_trigEvt_

edm::EDGetTokenT< edm::TriggerResults >	tok_trigRes_

TTree *	tree

edm::InputTag	triggerEvent_

std::vector< std::string >	trigNames_

int	verbosity_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 74 of file IsoTrackCalibration.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 132 of file IsoTrackCalibration.cc.

                                                                        : 
   nRun_(0) {
 
   usesResource("TFileService");
   //now do whatever initialization is needed
   const double isolationRadius(28.9);
   verbosity_                          = iConfig.getUntrackedParameter<int>("Verbosity",0);
   trigNames_                          = iConfig.getUntrackedParameter<std::vector<std::string> >("Triggers");
   theTrackQuality_                    = iConfig.getUntrackedParameter<std::string>("TrackQuality","highPurity");
   reco::TrackBase::TrackQuality trackQuality_=reco::TrackBase::qualityByName(theTrackQuality_);
   constTrackPt_                       = iConfig.getUntrackedParameter<double>("ConstTrackPt", 10.0);
   slopeTrackPt_                       = iConfig.getUntrackedParameter<double>("SlopeTrackPt", 0.16);
   selectionParameters_.minPt          = constTrackPt_;
   selectionParameters_.minQuality     = trackQuality_;
   selectionParameters_.maxDxyPV       = iConfig.getUntrackedParameter<double>("MaxDxyPV", 0.02);
   selectionParameters_.maxDzPV        = iConfig.getUntrackedParameter<double>("MaxDzPV",  0.02);
   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);
   a_coneR_                            = iConfig.getUntrackedParameter<double>("ConeRadius",34.98);
   a_charIsoR_                         = a_coneR_ + isolationRadius;
   a_mipR_                             = iConfig.getUntrackedParameter<double>("ConeRadiusMIP",14.0);
   maxEcalEnr_                         = iConfig.getUntrackedParameter<double>("MaxEcalEnergyInCone",2.5);
   maxNeighborTrackEnr_                = iConfig.getUntrackedParameter<double>("MaxNeighborTrackEnergy",40.0);
   isMC_                               = iConfig.getUntrackedParameter<bool>("IsMC", false);
   isQCD_                              = iConfig.getUntrackedParameter<bool>("IsQCD", false);
   isAOD_                              = iConfig.getUntrackedParameter<bool>("IsAOD", true);
   triggerEvent_                       = edm::InputTag("hltTriggerSummaryAOD","","HLT");
   theTriggerResultsLabel_             = edm::InputTag("TriggerResults","","HLT");
 
   // define tokens for access
   tok_trigEvt_  = consumes<trigger::TriggerEvent>(triggerEvent_);
   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_ew_       = consumes<GenEventInfoProduct>(edm::InputTag("generator")); 
  
   if (isAOD_) {
     tok_EB_     = consumes<EcalRecHitCollection>(edm::InputTag("reducedEcalRecHitsEB"));
     tok_EE_     = consumes<EcalRecHitCollection>(edm::InputTag("reducedEcalRecHitsEE"));
     tok_hbhe_   = consumes<HBHERecHitCollection>(edm::InputTag("reducedHcalRecHits", "hbhereco"));
   } else {
     tok_EB_     = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit","EcalRecHitsEB"));
     tok_EE_     = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit","EcalRecHitsEE"));
     tok_hbhe_   = consumes<HBHERecHitCollection>(edm::InputTag("hbhereco"));
   }
 
 #ifdef DebugLog
   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_mipR "          << a_mipR_
           <<"\t isMC "            << isMC_
           <<"\t isQCD "           << isQCD_
           <<"\t isAOD "           << isAOD_
           << std::endl;
     std::cout << trigNames_.size() << " triggers to be studied:";
     for (unsigned int k=0; k<trigNames_.size(); ++k)
       std::cout << " " << trigNames_[k];
     std::cout << std::endl;
   }
 #endif
 }

IsoTrackCalibration::~IsoTrackCalibration ( )

override

Definition at line 212 of file IsoTrackCalibration.cc.

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

Member Function Documentation

void IsoTrackCalibration::analyze	(	edm::Event const &	,
		edm::EventSetup const &
	)

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 217 of file IsoTrackCalibration.cc.

References a_charIsoR_, a_coneR_, a_mipR_, funct::abs(), edm::HLTGlobalStatus::accept(), ecalTB2006H4_GenSimDigiReco_cfg::bField, edm::EventBase::bunchCrossing(), spr::chargeIsolationCone(), constTrackPt_, gather_cfg::cout, spr::eCone_ecal(), spr::eCone_hcal(), reco::TrackBase::eta(), edm::EventID::event(), spr::find(), particlelevel_cff::genEventInfo, edm::EventSetup::get(), edm::Event::getByToken(), spr::goodTrack(), h_nTrk, h_nVtx, photonIsolationHIProducer_cfi::hbhe, ValidationMatrix::hlt, mps_fire::i, edm::EventBase::id(), HcalDetId::ieta(), isMC_, isQCD_, edm::HandleBase::isValid(), gen::k, edm::InputTag::label(), edm::EventBase::luminosityBlock(), spr::trackSelectionParameters::maxDxyPV, spr::trackSelectionParameters::maxDzPV, maxEcalEnr_, spr::trackSelectionParameters::maxInMiss, maxNeighborTrackEnr_, spr::trackSelectionParameters::maxOutMiss, spr::trackSelectionParameters::minPt, nTracks(), reco::TrackBase::p(), reco::TrackBase::phi(), reco::BeamSpot::position(), edm::ESHandle< T >::product(), spr::propagateCALO(), reco::TrackBase::pt(), edm::EventID::run(), selectionParameters_, edm::HLTGlobalStatus::size(), slopeTrackPt_, t_DetIds, t_DetIds1, t_DetIds3, t_eHcal, t_eHcal10, t_eHcal30, t_eMipDR, t_Event, t_EventWeight, t_HitEnergies, t_HitEnergies1, t_HitEnergies3, t_hmaxNearP, t_ieta, t_nTrk, t_nVtx, t_p, t_phi, t_pt, t_qltyMissFlag, t_qltyPVFlag, t_Run, t_selectTk, theTrackQuality_, tok_bs_, tok_EB_, tok_EE_, tok_ew_, tok_genTrack_, tok_hbhe_, tok_recVtx_, tok_trigEvt_, tok_trigRes_, triggerEvent_, edm::TriggerNames::triggerNames(), edm::Event::triggerNames(), triggerResults, trigNames_, verbosity_, GenEventInfoProduct::weight(), x, y, and z.

                                                {
 
   t_Run   = iEvent.id().run();
   t_Event = iEvent.id().event();
 #ifdef DebugLog
   if (verbosity_%10 > 0) 
     std::cout << "Run " << t_Run << " Event " << t_Event << " Luminosity " 
           << iEvent.luminosityBlock() << " Bunch " << iEvent.bunchCrossing()
           << " starts ==========" << std::endl;
 #endif
   //Get magnetic field and ECAL channel status
   edm::ESHandle<MagneticField> bFieldH;
   iSetup.get<IdealMagneticFieldRecord>().get(bFieldH);
   const MagneticField *bField = bFieldH.product();
   
   edm::ESHandle<EcalSeverityLevelAlgo> sevlv;
   iSetup.get<EcalSeverityLevelAlgoRcd>().get(sevlv);
 
   // get handles to calogeometry and calotopology
   edm::ESHandle<CaloGeometry> pG;
   iSetup.get<CaloGeometryRecord>().get(pG);
   const CaloGeometry* geo = pG.product();
   
   //Get track collection
   edm::Handle<reco::TrackCollection> trkCollection;
   iEvent.getByToken(tok_genTrack_, trkCollection);
   reco::TrackCollection::const_iterator trkItr;
  
   //event weight for FLAT sample and PU information
   t_EventWeight = 1.0;
   edm::Handle<GenEventInfoProduct> genEventInfo;
   iEvent.getByToken(tok_ew_, genEventInfo);
   if (genEventInfo.isValid()) t_EventWeight = genEventInfo->weight();  
 
   //Define the best vertex and the beamspot
   edm::Handle<reco::VertexCollection> recVtxs;
   iEvent.getByToken(tok_recVtx_, recVtxs);  
   edm::Handle<reco::BeamSpot> beamSpotH;
   iEvent.getByToken(tok_bs_, beamSpotH);
   math::XYZPoint leadPV(0,0,0);
 
   t_nVtx = recVtxs->size();
   h_nVtx->Fill(t_nVtx);
 
   if (!recVtxs->empty() && !((*recVtxs)[0].isFake())) {
     leadPV = math::XYZPoint( (*recVtxs)[0].x(),(*recVtxs)[0].y(), (*recVtxs)[0].z() );
   } else if (beamSpotH.isValid()) {
     leadPV = beamSpotH->position();
   }
 #ifdef DebugLog
   if ((verbosity_/100)%10>0) {
     std::cout << "Primary Vertex " << leadPV;
     if (beamSpotH.isValid()) std::cout << " Beam Spot " 
                        << beamSpotH->position();
     std::cout << std::endl;
   }
 #endif  
 
   // RecHits
   edm::Handle<EcalRecHitCollection> barrelRecHitsHandle;
   edm::Handle<EcalRecHitCollection> endcapRecHitsHandle;
   iEvent.getByToken(tok_EB_, barrelRecHitsHandle);
   iEvent.getByToken(tok_EE_, endcapRecHitsHandle);
   edm::Handle<HBHERecHitCollection> hbhe;
   iEvent.getByToken(tok_hbhe_, hbhe);
 
   //Trigger
   bool triggerOK = false;
   if (isMC_ && !isQCD_) {
     triggerOK = true; // ignore HLT for single pion MC
   } else {
     trigger::TriggerEvent triggerEvent;
     edm::Handle<trigger::TriggerEvent> triggerEventHandle;
     iEvent.getByToken(tok_trigEvt_, triggerEventHandle);
     if (!triggerEventHandle.isValid()) {
 #ifdef DebugLog
       std::cout << "Error! Can't get the product "<< triggerEvent_.label() 
         << std::endl;
 #endif
     } else {
       edm::Handle<edm::TriggerResults> triggerResults;
       iEvent.getByToken(tok_trigRes_, triggerResults);
       if (triggerResults.isValid()) {
     const edm::TriggerNames & triggerNames = iEvent.triggerNames(*triggerResults);
     const std::vector<std::string> & triggerNames_ = triggerNames.triggerNames();
     for (unsigned int iHLT=0; iHLT<triggerResults->size(); iHLT++) {
       int hlt    = triggerResults->accept(iHLT);
       if (hlt > 0) {
         for (unsigned int i=0; i<trigNames_.size(); ++i) {
           if (triggerNames_[iHLT].find(trigNames_[i])!=std::string::npos) {
         triggerOK = true;
 #ifdef DebugLog
         if (verbosity_%10 > 0)
           std::cout << "This is the trigger we are looking for "
                 << triggerNames_[iHLT] << " Flag " << hlt 
                 << std::endl;
 #endif
           }
         }
       }
     }
       }
     }
   }
 
   if (triggerOK) {
     //Propagate tracks to calorimeter surface)
     std::vector<spr::propagatedTrackDirection> trkCaloDirections;
     spr::propagateCALO(trkCollection, geo, bField, theTrackQuality_,
                trkCaloDirections, ((verbosity_/100)%10>2));
     //Loop over tracks
     std::vector<spr::propagatedTrackDirection>::const_iterator trkDetItr;
     unsigned int nTracks(0), nselTracks(0);
     t_nTrk = trkCaloDirections.size();
 
     for (trkDetItr = trkCaloDirections.begin(),nTracks=0; 
      trkDetItr != trkCaloDirections.end(); trkDetItr++, nTracks++) {
 
       const reco::Track* pTrack = &(*(trkDetItr->trkItr));
 #ifdef DebugLog
       if (verbosity_%10> 0) 
     std::cout << "This track : " << nTracks << " (pt/eta/phi/p) :" 
           << pTrack->pt() << "/" << pTrack->eta() << "/" 
           << pTrack->phi() << "/" << pTrack->p() << std::endl;
 #endif
 
       t_ieta = 0;
       if (trkDetItr->okHCAL) {
     HcalDetId detId = (HcalDetId)(trkDetItr->detIdHCAL);
     t_ieta = detId.ieta();
       }
       // ---------- eta-dependent restriction on Pt ----------------------------
       selectionParameters_.minPt = constTrackPt_ - std::abs((double)t_ieta)*slopeTrackPt_;
       // -----------------------------------------------------------------------
         
       //Selection of good track
       t_selectTk = spr::goodTrack(pTrack,leadPV,selectionParameters_,((verbosity_/100)%10>2));
       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>2));
       oneCutParameters           = selectionParameters_;
       oneCutParameters.maxDxyPV  = 10;
       oneCutParameters.maxDzPV   = 100;
       t_qltyMissFlag = spr::goodTrack(pTrack,leadPV,oneCutParameters,((verbosity_/100)%10>2));
       oneCutParameters           = selectionParameters_;
       oneCutParameters.maxInMiss = 2;
       oneCutParameters.maxOutMiss= 2;
       t_qltyPVFlag   = spr::goodTrack(pTrack,leadPV,oneCutParameters,((verbosity_/100)%10>2));
         
 #ifdef DebugLog
       if (verbosity_%10 > 0) 
     std::cout << "qltyFlag|okECAL|okHCAL : " << qltyFlag << "|" 
           << trkDetItr->okECAL << "/" << trkDetItr->okHCAL 
           << std::endl;
 #endif
       if (qltyFlag && trkDetItr->okECAL && trkDetItr->okHCAL) {
     nselTracks++;
     int nRH_eMipDR(0), nNearTRKs(0);
     //------ ecal energy around track -------------------------------
     t_eMipDR = spr::eCone_ecal(geo, barrelRecHitsHandle, 
                    endcapRecHitsHandle,
                    trkDetItr->pointHCAL,
                    trkDetItr->pointECAL,
                    a_mipR_, trkDetItr->directionECAL, 
                    nRH_eMipDR);
     //---- isolation criteria ----------------------------------------------
     t_hmaxNearP = spr::chargeIsolationCone(nTracks,trkCaloDirections,
                            a_charIsoR_, nNearTRKs, 
                            ((verbosity_/100)%10>2));
 
     if ( t_eMipDR < maxEcalEnr_ && t_hmaxNearP < maxNeighborTrackEnr_ ) {
       //------------ HCAL --------------------------------------------------
       //------ initialize arrays of DetID and hit energies -----------------
       t_DetIds->clear();
       t_DetIds1->clear();
       t_DetIds3->clear();
       t_HitEnergies->clear();
       t_HitEnergies1->clear();
       t_HitEnergies3->clear();
       int nRecHits(-999), nRecHits1(-999), nRecHits3(-999);       
       std::vector<DetId> ids, ids1, ids3;
       //------ hcal energy in the main cone -------------------------------
       t_eHcal = spr::eCone_hcal(geo, hbhe, trkDetItr->pointHCAL, 
                     trkDetItr->pointECAL, a_coneR_, 
                     trkDetItr->directionHCAL,nRecHits, 
                     ids, *t_HitEnergies);
       t_DetIds->reserve(ids.size());
       for (unsigned int k = 0; k < ids.size(); ++k) {
         t_DetIds->push_back(ids[k].rawId());
       }
       //----- hcal energy in the extended cone 1 (a_coneR+10) --------------
       t_eHcal10 = spr::eCone_hcal(geo, hbhe, trkDetItr->pointHCAL,
                       trkDetItr->pointECAL, a_coneR_+10,
                       trkDetItr->directionHCAL,nRecHits1,
                       ids1, *t_HitEnergies1);
       t_DetIds1->reserve(ids1.size());
       for (unsigned int k = 0; k < ids1.size(); ++k) {
         t_DetIds1->push_back(ids1[k].rawId());
       }
       //----- hcal energy in the extended cone 3 (a_coneR+30) --------------
       t_eHcal30 = spr::eCone_hcal(geo, hbhe, trkDetItr->pointHCAL,
                       trkDetItr->pointECAL, a_coneR_+30,
                       trkDetItr->directionHCAL,nRecHits3,
                       ids3, *t_HitEnergies3);
       t_DetIds3->reserve(ids3.size());
       for (unsigned int k = 0; k < ids3.size(); ++k) {
         t_DetIds3->push_back(ids3[k].rawId());
       }
 
       t_p = pTrack->p();
       t_pt = pTrack->pt();
       t_phi = pTrack->phi();
 #ifdef DebugLog
       if (verbosity_%10 > 0) {
         std::cout << "This track : " << nTracks << " (pt/eta/phi/p) :" 
               << pTrack->pt() << "/" << pTrack->eta() << "/" 
               << pTrack->phi() << "/" << t_p << std::endl;
         std::cout << "e_MIP " << t_eMipDR << " Chg Isolation "
               << t_hmaxNearP << " eHcal" << t_eHcal << " ieta " 
               << t_ieta << " Quality " << t_qltyMissFlag
               << ":" << t_qltyPVFlag << ":" << t_selectTk 
               << std::endl;
         for (unsigned int lll=0;lll<t_DetIds->size();lll++) {
           std::cout << "det id is = " << t_DetIds->at(lll)  << "  " 
             << " hit enery is  = "  << t_HitEnergies->at(lll) 
             << std::endl;         
         }
       }
 #endif
       tree->Fill();
     } // end of conditions on t_eMipDR and t_hmaxNearP
       } // end of loose check of track quality
     } // end of loop over tracks
 
     h_nTrk->Fill(nTracks);
   } // end of triggerOK
 }

void IsoTrackCalibration::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 460 of file IsoTrackCalibration.cc.

References fs_, h_nTrk, h_nVtx, TFileService::make(), t_DetIds, t_DetIds1, t_DetIds3, t_eHcal, t_eHcal10, t_eHcal30, t_eMipDR, t_Event, t_EventWeight, t_HitEnergies, t_HitEnergies1, t_HitEnergies3, t_hmaxNearP, t_ieta, t_nTrk, t_nVtx, t_p, t_phi, t_pt, t_qltyMissFlag, t_qltyPVFlag, t_Run, and t_selectTk.

                                    {
 
   h_nVtx = fs_->make<TH1F>("h_nVtx","h_nVtx",100,0,100);
   h_nTrk = fs_->make<TH1F>("h_nTrk","h_nTrk",100,0,2000);
 
   tree = fs_->make<TTree>("CalibTree", "CalibTree");
    
   tree->Branch("t_Run",         &t_Run,         "t_Run/I");
   tree->Branch("t_Event",       &t_Event,       "t_Event/I");
   tree->Branch("t_nVtx",        &t_nVtx,        "t_nVtx/I");
   tree->Branch("t_nTrk",        &t_nTrk,        "t_nTrk/I");
   tree->Branch("t_EventWeight", &t_EventWeight, "t_EventWeight/D");
   tree->Branch("t_p",           &t_p,           "t_p/D");
   tree->Branch("t_pt",          &t_pt,          "t_pt/D");
   tree->Branch("t_ieta",        &t_ieta,        "t_ieta/I");
   tree->Branch("t_phi",         &t_phi,         "t_phi/D");
   tree->Branch("t_eMipDR",      &t_eMipDR,      "t_eMipDR/D");
   tree->Branch("t_eHcal",       &t_eHcal,       "t_eHcal/D");
   tree->Branch("t_eHcal10",     &t_eHcal10,     "t_eHcal10/D");
   tree->Branch("t_eHcal30",     &t_eHcal30,     "t_eHcal30/D");
   tree->Branch("t_hmaxNearP",   &t_hmaxNearP,   "t_hmaxNearP/D");
   tree->Branch("t_selectTk",    &t_selectTk,    "t_selectTk/O");
   tree->Branch("t_qltyMissFlag",&t_qltyMissFlag,"t_qltyMissFlag/O");
   tree->Branch("t_qltyPVFlag",  &t_qltyPVFlag,  "t_qltyPVFlag/O)");
 
   t_DetIds       = new std::vector<unsigned int>();
   t_DetIds1      = new std::vector<unsigned int>();
   t_DetIds3      = new std::vector<unsigned int>();
   t_HitEnergies  = new std::vector<double>();
   t_HitEnergies1 = new std::vector<double>();
   t_HitEnergies3 = new std::vector<double>();
 
   tree->Branch("t_DetIds",       "std::vector<unsigned int>", &t_DetIds);
 //tree->Branch("t_DetIds1",      "std::vector<unsigned int>", &t_DetIds1);
 //tree->Branch("t_DetIds3",      "std::vector<unsigned int>", &t_DetIds3);
   tree->Branch("t_HitEnergies",  "std::vector<double>",       &t_HitEnergies);
 //tree->Branch("t_HitEnergies1", "std::vector<double>",       &t_HitEnergies1);
 //tree->Branch("t_HitEnergies3", "std::vector<double>",       &t_HitEnergies3);
 }

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

privatevirtual

Definition at line 535 of file IsoTrackCalibration.cc.

535 {}

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

overrideprivate

Definition at line 501 of file IsoTrackCalibration.cc.

References RemoveAddSevLevel::flag, hltConfig_, HLTConfigProvider::init(), gen::n, nRun_, edm::RunBase::run(), HLTConfigProvider::size(), HLTConfigProvider::triggerIndex(), and trigNames_.

                                                                                   {
 
   bool changed_(true);
   bool flag =  hltConfig_.init(iRun,iSetup,"HLT",changed_);
   edm::LogInfo("HcalIsoTrack") << "Run[" << nRun_ << "] " << iRun.run() 
                                << " process HLT init flag "
                                << flag << " change flag " << changed_;
 
   // check if trigger names in (new) config                       
   if (changed_) {
     changed_ = false;
     edm::LogInfo("HcalIsoTrack") << "New trigger menu found !!!";
     const unsigned int n(hltConfig_.size());
     for (unsigned itrig=0; itrig<trigNames_.size(); itrig++) {
       unsigned int triggerindx = hltConfig_.triggerIndex(trigNames_[itrig]);
       if (triggerindx >= n) {
         edm::LogWarning("HcalIsoTrack") << trigNames_[itrig] << " " 
                                         << triggerindx << " does not exist in "
                                         << "the current menu";
       } else {
         edm::LogInfo("HcalIsoTrack") << trigNames_[itrig] << " " 
                                      << triggerindx << " exists";
       }
     }
   }
 }

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

privatevirtual

Definition at line 537 of file IsoTrackCalibration.cc.

537 {}

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

overrideprivate

Definition at line 529 of file IsoTrackCalibration.cc.

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

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

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

static

Definition at line 539 of file IsoTrackCalibration.cc.

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

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

Member Data Documentation

double IsoTrackCalibration::a_charIsoR_

private

Definition at line 96 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

double IsoTrackCalibration::a_coneR_

private

Definition at line 96 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

double IsoTrackCalibration::a_mipR_

private

Definition at line 96 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

double IsoTrackCalibration::constTrackPt_

private

Definition at line 98 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

edm::Service<TFileService> IsoTrackCalibration::fs_

private

Definition at line 90 of file IsoTrackCalibration.cc.

Referenced by beginJob().

TH1I* IsoTrackCalibration::h_iEta

private

Definition at line 125 of file IsoTrackCalibration.cc.

TH1F * IsoTrackCalibration::h_jetp

private

Definition at line 127 of file IsoTrackCalibration.cc.

TH1F* IsoTrackCalibration::h_jetpt[4]

private

Definition at line 128 of file IsoTrackCalibration.cc.

TH1F* IsoTrackCalibration::h_nTrk

private

Definition at line 123 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

TH1F * IsoTrackCalibration::h_nVtx

private

Definition at line 123 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

TH1F* IsoTrackCalibration::h_Rechit_E

private

Definition at line 127 of file IsoTrackCalibration.cc.

TProfile* IsoTrackCalibration::h_RecHit_iEta

private

Definition at line 124 of file IsoTrackCalibration.cc.

TProfile * IsoTrackCalibration::h_RecHit_num

private

Definition at line 124 of file IsoTrackCalibration.cc.

TH1I * IsoTrackCalibration::h_tketa0[5]

private

Definition at line 125 of file IsoTrackCalibration.cc.

TH1I * IsoTrackCalibration::h_tketa1[5]

private

Definition at line 125 of file IsoTrackCalibration.cc.

TH1I * IsoTrackCalibration::h_tketa2[5]

private

Definition at line 125 of file IsoTrackCalibration.cc.

TH1I* IsoTrackCalibration::h_tketa3[5]

private

Definition at line 126 of file IsoTrackCalibration.cc.

TH1I * IsoTrackCalibration::h_tketa4[5]

private

Definition at line 126 of file IsoTrackCalibration.cc.

TH1I * IsoTrackCalibration::h_tketa5[5]

private

Definition at line 126 of file IsoTrackCalibration.cc.

TH1I* IsoTrackCalibration::h_tketav1[5][6]

private

Definition at line 129 of file IsoTrackCalibration.cc.

TH1I * IsoTrackCalibration::h_tketav2[5][6]

private

Definition at line 129 of file IsoTrackCalibration.cc.

HLTConfigProvider IsoTrackCalibration::hltConfig_

private

Definition at line 91 of file IsoTrackCalibration.cc.

Referenced by beginRun().

std::vector<std::string> IsoTrackCalibration::HLTNames_

private

Definition at line 92 of file IsoTrackCalibration.cc.

bool IsoTrackCalibration::isAOD_

private

Definition at line 97 of file IsoTrackCalibration.cc.

Referenced by IsoTrackCalibration().

bool IsoTrackCalibration::isMC_

private

Definition at line 97 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

bool IsoTrackCalibration::isQCD_

private

Definition at line 97 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

double IsoTrackCalibration::maxEcalEnr_

private

Definition at line 99 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

double IsoTrackCalibration::maxNeighborTrackEnr_

private

Definition at line 100 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

int IsoTrackCalibration::nRun_

private

Definition at line 101 of file IsoTrackCalibration.cc.

Referenced by beginRun(), and endRun().

spr::trackSelectionParameters IsoTrackCalibration::selectionParameters_

private

Definition at line 94 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

double IsoTrackCalibration::slopeTrackPt_

private

Definition at line 98 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

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

private

Definition at line 120 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

std::vector<unsigned int> * IsoTrackCalibration::t_DetIds1

private

Definition at line 120 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

std::vector<unsigned int> * IsoTrackCalibration::t_DetIds3

private

Definition at line 120 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalibration::t_eHcal

private

Definition at line 117 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalibration::t_eHcal10

private

Definition at line 117 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalibration::t_eHcal30

private

Definition at line 117 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalibration::t_eMipDR

private

Definition at line 118 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

int IsoTrackCalibration::t_Event

private

Definition at line 115 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalibration::t_EventWeight

private

Definition at line 116 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

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

private

Definition at line 121 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

std::vector<double> * IsoTrackCalibration::t_HitEnergies1

private

Definition at line 121 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

std::vector<double> * IsoTrackCalibration::t_HitEnergies3

private

Definition at line 121 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalibration::t_hmaxNearP

private

Definition at line 118 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

int IsoTrackCalibration::t_ieta

private

Definition at line 115 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

int IsoTrackCalibration::t_nTrk

private

Definition at line 115 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

int IsoTrackCalibration::t_nVtx

private

Definition at line 115 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalibration::t_p

private

Definition at line 116 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalibration::t_phi

private

Definition at line 116 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalibration::t_pt

private

Definition at line 116 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

bool IsoTrackCalibration::t_qltyMissFlag

private

Definition at line 119 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

bool IsoTrackCalibration::t_qltyPVFlag

private

Definition at line 119 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

int IsoTrackCalibration::t_Run

private

Definition at line 115 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

bool IsoTrackCalibration::t_selectTk

private

Definition at line 119 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

std::string IsoTrackCalibration::theTrackQuality_

private

Definition at line 95 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

edm::InputTag IsoTrackCalibration::theTriggerResultsLabel_

private

Definition at line 102 of file IsoTrackCalibration.cc.

Referenced by IsoTrackCalibration().

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

private

Definition at line 108 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

edm::EDGetTokenT<EcalRecHitCollection> IsoTrackCalibration::tok_EB_

private

Definition at line 109 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

edm::EDGetTokenT<EcalRecHitCollection> IsoTrackCalibration::tok_EE_

private

Definition at line 110 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

edm::EDGetTokenT<GenEventInfoProduct> IsoTrackCalibration::tok_ew_

private

Definition at line 112 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

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

private

Definition at line 106 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

edm::EDGetTokenT<HBHERecHitCollection> IsoTrackCalibration::tok_hbhe_

private

Definition at line 111 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

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

private

Definition at line 107 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

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

private

Definition at line 103 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

edm::EDGetTokenT<edm::TriggerResults> IsoTrackCalibration::tok_trigRes_

private

Definition at line 104 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

TTree* IsoTrackCalibration::tree

private

Definition at line 114 of file IsoTrackCalibration.cc.

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

edm::InputTag IsoTrackCalibration::triggerEvent_

private

Definition at line 102 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

std::vector<std::string> IsoTrackCalibration::trigNames_

private

Definition at line 92 of file IsoTrackCalibration.cc.

Referenced by analyze(), beginRun(), and IsoTrackCalibration().

int IsoTrackCalibration::verbosity_

private

Definition at line 93 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

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