EopElecTreeWriter Class Reference

Inheritance diagram for EopElecTreeWriter:

Public Member Functions
	EopElecTreeWriter (const edm::ParameterSet &)
	~EopElecTreeWriter () override
Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::SharedResources >
	EDAnalyzer ()=default
	EDAnalyzer (const EDAnalyzer &)=delete
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
const EDAnalyzer &	operator= (const EDAnalyzer &)=delete
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
bool	wantsInputProcessBlocks () const final
bool	wantsProcessBlocks () 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
ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
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::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &)
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 (const edm::Event &, const edm::EventSetup &) override
void	beginRun (const edm::Run &, const edm::EventSetup &) override
void	endJob () override
void	endRun (const edm::Run &, const edm::EventSetup &) override

Private Attributes
const edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeomToken_
bool	debugTriggerSelection_
edm::Service< TFileService >	fs_
TH1D *	h_counter1
TH1D *	h_counter2
TH1D *	h_cut_OneSCmatch
TH1D *	h_cut_Ptmin
TH1D *	h_distToClosestSC
TH1D *	h_distToClosestSCgsf
TH1D *	h_EcalEnergy
TH1D *	h_Eop_InnerNegative
TH1D *	h_Eop_InnerPositive
TH1D *	h_fBREM
TH1D *	h_Momentum
TH1D *	h_nEvents
TH1D *	h_nEvents2Elec
TH1D *	h_nEventsHLTelectron
TH1D *	h_nEventsHLTFilter
TH1D *	h_nEventsHLTrejected
TH1D *	h_nEventsTriggered
TH1D *	h_nEventsWithVertex
TH1D *	h_nHLTelectrons
TH1D *	h_nTracks
TH1D *	h_nTracksFiltered
TH1D *	h_nTrkRejectedPerEvt
TH1D *	h_nTrkSelectedPerEvt
TH1D *	HcalEnergy
TH2D *	HcalVSEcal
HLTConfigProvider	hltConfig_
const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	magFieldToken_
MultiTrajectoryStateTransform *	mtsTransform_
const edm::EDGetTokenT< reco::SuperClusterCollection >	theBarrelSupClusCollectionToken_
const edm::EDGetTokenT< EcalRecHitCollection >	theEcalRecHitCollectionToken_
const edm::EDGetTokenT< reco::SuperClusterCollection >	theEndCapSupClusCollectionToken_
std::string	theFilter_
const edm::EDGetTokenT< reco::GsfElectronCoreCollection >	theGsfElectronCoreCollectionToken_
const edm::EDGetTokenT< reco::GsfTrackCollection >	theGsfTrackCollectionToken_
const edm::EDGetTokenT< HBHERecHitCollection >	theHBHERecHitCollectionToken_
std::string	theTrigger_
const edm::EDGetTokenT< trigger::TriggerEvent >	theTriggerEventToken_
const edm::EDGetTokenT< edm::TriggerResults >	theTriggerResultsToken_
const edm::EDGetTokenT< reco::VertexCollection >	theVertexCollectionToken_
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	tkGeomToken_
TTree *	tree_
EopElecVariables *	treeMemPtr_
std::vector< std::string >	triggerNames_

Static Private Attributes
static constexpr float	k_etaBarrel = 1.55
static constexpr float	k_etaEndcap = 1.44

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)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
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)
void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 101 of file EopElecTreeWriter.cc.

Constructor & Destructor Documentation

EopElecTreeWriter()

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

explicit

Definition at line 222 of file EopElecTreeWriter.cc.

References fs_, h_counter1, h_counter2, h_cut_OneSCmatch, h_cut_Ptmin, h_distToClosestSC, h_distToClosestSCgsf, h_EcalEnergy, h_Eop_InnerNegative, h_Eop_InnerPositive, h_fBREM, h_Momentum, h_nEvents, h_nEvents2Elec, h_nEventsHLTelectron, h_nEventsHLTFilter, h_nEventsHLTrejected, h_nEventsTriggered, h_nEventsWithVertex, h_nHLTelectrons, h_nTracks, h_nTracksFiltered, h_nTrkRejectedPerEvt, h_nTrkSelectedPerEvt, HcalEnergy, HcalVSEcal, TFileService::kSharedResource, TFileService::make(), tree_, and treeMemPtr_.

    : magFieldToken_(esConsumes()),
      tkGeomToken_(esConsumes()),
      caloGeomToken_(esConsumes()),
      theVertexCollectionToken_(consumes<reco::VertexCollection>(edm::InputTag("offlinePrimaryVertices"))),
      theHBHERecHitCollectionToken_(consumes<HBHERecHitCollection>(edm::InputTag("hbhereco", ""))),
      theEcalRecHitCollectionToken_(consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit", "EcalRecHitsEB"))),
      theBarrelSupClusCollectionToken_(
          consumes<reco::SuperClusterCollection>(edm::InputTag("hybridSuperClusters", ""))),
      theEndCapSupClusCollectionToken_(consumes<reco::SuperClusterCollection>(
          edm::InputTag("multi5x5SuperClusters", "multi5x5EndcapSuperClusters"))),
      theTriggerResultsToken_(consumes<edm::TriggerResults>(edm::InputTag("TriggerResults", "", "HLT"))),
      theTriggerEventToken_(consumes<trigger::TriggerEvent>(edm::InputTag("hltTriggerSummaryAOD"))),
      theGsfTrackCollectionToken_(consumes<reco::GsfTrackCollection>(iConfig.getParameter<edm::InputTag>("src"))),
      theGsfElectronCoreCollectionToken_(
          consumes<reco::GsfElectronCoreCollection>(edm::InputTag("gedGsfElectronCores"))),
      theTrigger_(iConfig.getParameter<std::string>("triggerPath")),
      theFilter_(iConfig.getParameter<std::string>("hltFilter")),
      debugTriggerSelection_(iConfig.getParameter<bool>("debugTriggerSelection")) {
  usesResource(TFileService::kSharedResource);

  // TTree creation
  tree_ = fs_->make<TTree>("EopTree", "EopTree");
  treeMemPtr_ = new EopElecVariables;
  tree_->Branch("EopElecVariables", &treeMemPtr_);  // address of pointer!

  // Control histograms declaration
  h_distToClosestSC = fs_->make<TH1D>("distToClosestSC", "distToClosestSC", 100, 0, 0.1);
  h_distToClosestSCgsf = fs_->make<TH1D>("distToClosestSCgsf", "distToClosestSCgsf", 100, 0, 0.1);
  h_EcalEnergy = fs_->make<TH1D>("EcalEnergy", "EcalEnergy", 100, 0, 200);
  h_Momentum = fs_->make<TH1D>("Momentum", "Momentum", 100, 0, 200);
  HcalEnergy = fs_->make<TH1D>("HcalEnergy", "HcalEnergy", 100, 0, 40);
  h_fBREM = fs_->make<TH1D>("fBREM", "fBREM", 100, -0.2, 1);
  h_Eop_InnerNegative = fs_->make<TH1D>("Eop_InnerNegative", "Eop_InnerNegative", 100, 0, 3);
  h_Eop_InnerPositive = fs_->make<TH1D>("Eop_InnerPositive", "Eop_InnerPositive", 100, 0, 3);
  HcalVSEcal = fs_->make<TH2D>("HcalVSEcal", "HcalVSEcal", 100, 0, 160, 100, 0, 10);

  h_nEvents = fs_->make<TH1D>("nEvents", "nEvents", 1, 0, 1);
  h_nEventsWithVertex = fs_->make<TH1D>("nEventsWithVertex", "nEventsWithVertex", 1, 0, 1);
  h_nEventsTriggered = fs_->make<TH1D>("nEventsTriggered", "nEventsTriggered", 1, 0, 1);
  h_nEventsHLTFilter = fs_->make<TH1D>("nEventsHLTFilter", "nEventsHLTFilter", 1, 0, 1);
  h_nEventsHLTelectron = fs_->make<TH1D>("nEventsHLTelectron", "nEventsHLTelectron", 1, 0, 1);
  h_nEventsHLTrejected = fs_->make<TH1D>("nEventsHLTrejected", "nEventsHLTrejected", 1, 0, 1);
  h_nEvents2Elec = fs_->make<TH1D>("nEvents2Elec", "nEvents2Elec", 1, 0, 1);

  h_nHLTelectrons = fs_->make<TH1D>("nHLTelectrons", "nHLTelectrons", 20, 0, 20);
  h_nTrkRejectedPerEvt = fs_->make<TH1D>("nTrkRejectedPerEvt", "nTrkRejectedPerEvt", 20, 0, 20);
  h_nTrkSelectedPerEvt = fs_->make<TH1D>("nTrkSelectedPerEvt", "nTrkSelectedPerEvt", 20, 0, 20);

  h_nTracks = fs_->make<TH1D>("nTracks", "nTracks", 1, 0, 1);
  h_nTracksFiltered = fs_->make<TH1D>("nTracksFiltered", "nTracksFiltered", 1, 0, 1);
  h_cut_Ptmin = fs_->make<TH1D>("cut_Ptmin", "cut_Ptmin", 1, 0, 1);
  h_cut_OneSCmatch = fs_->make<TH1D>("cut_OneSCmatch", "cut_OneSCmatch", 1, 0, 1);

  h_counter1 = fs_->make<TH1D>("counter1", "counter1", 1, 0, 1);
  h_counter2 = fs_->make<TH1D>("counter2", "counter2", 1, 0, 1);
}

~EopElecTreeWriter()

EopElecTreeWriter::~EopElecTreeWriter ( )

override

Definition at line 280 of file EopElecTreeWriter.cc.

References h_distToClosestSC, h_distToClosestSCgsf, h_EcalEnergy, h_Eop_InnerNegative, h_Eop_InnerPositive, h_Momentum, HcalEnergy, and HcalVSEcal.

                                      {
  // control histograms
  h_distToClosestSC->SetXTitle("distance from track to closest SuperCluster in eta-phi plan (weighted matching)");
  h_distToClosestSC->SetYTitle("# Tracks");

  h_distToClosestSCgsf->SetXTitle("distance from track to closest SuperCluster in eta-phi plan (gsfElectronCore)");
  h_distToClosestSCgsf->SetYTitle("# Tracks");

  h_EcalEnergy->SetXTitle("Ecal energy deposit (GeV)");
  h_EcalEnergy->SetYTitle("# tracks");

  HcalEnergy->SetXTitle("Hcal energy deposit (GeV)");
  HcalEnergy->SetYTitle("# tracks");

  h_Momentum->SetXTitle("Momentum magnitude (GeV)");
  h_Momentum->SetYTitle("# tracks");

  h_Eop_InnerNegative->SetXTitle("E/p");
  h_Eop_InnerNegative->SetYTitle("# tracks");

  h_Eop_InnerPositive->SetXTitle("E/p");
  h_Eop_InnerPositive->SetYTitle("# tracks");

  HcalVSEcal->SetXTitle("Ecal energy (GeV)");
  HcalVSEcal->SetYTitle("Hcal energy (GeV)");
}

Member Function Documentation

analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 311 of file EopElecTreeWriter.cc.

References funct::abs(), edm::HLTGlobalStatus::accept(), SuperClusterShapeAlgo::Calculate_Covariances(), caloGeomToken_, EopElecVariables::charge, EopElecVariables::chi2, debugTriggerSelection_, reco::deltaR(), HGC3DClusterGenMatchSelector_cfi::dR, EopElecVariables::dRto1stSC, EopElecVariables::dRto2ndSC, DetId::Ecal, EcalBarrel, eopUtils::ecalEta(), pwdgSkimBPark_cfi::electrons, edm::HLTGlobalStatus::error(), PVValHelper::eta, EopElecVariables::eta, PV3DBase< T, PVType, FrameType >::eta(), EopElecVariables::etaIn, EopElecVariables::etaOut, SuperClusterShapeAlgo::etaWidth(), electrons_cff::etaWidth, EopElecVariables::EvtNumber, EopElecVariables::fbrem, electrons_cff::fbrem, filters, EopTriggerType::fired, edm::EventSetup::getData(), CaloGeometry::getPosition(), CaloGeometry::getSubdetectorGeometry(), h_counter1, h_counter2, h_cut_OneSCmatch, h_cut_Ptmin, h_distToClosestSC, h_distToClosestSCgsf, h_EcalEnergy, h_Eop_InnerNegative, h_Eop_InnerPositive, h_fBREM, h_Momentum, h_nEvents, h_nEvents2Elec, h_nEventsHLTelectron, h_nEventsHLTFilter, h_nEventsHLTrejected, h_nEventsTriggered, h_nEventsWithVertex, h_nHLTelectrons, h_nTracks, h_nTracksFiltered, h_nTrkRejectedPerEvt, h_nTrkSelectedPerEvt, HcalEnergy, EopElecVariables::HcalEnergyIn01, EopElecVariables::HcalEnergyIn02, EopElecVariables::HcalEnergyIn03, EopElecVariables::HcalEnergyIn04, EopElecVariables::HcalEnergyIn05, HcalVSEcal, hltConfig_, mps_fire::i, iEvent, EopTriggerType::index, EopElecVariables::innerOk, MultiTrajectoryStateTransform::innerStateOnSurface(), PixelPluginsPhase0_cfi::isBarrel, EopElecVariables::isEcalDriven, electrons_cff::isEcalDriven, GeomDetEnumerators::isEndcap(), edm::Ref< C, T, F >::isNull(), EopElecVariables::isTrackerDriven, TrajectoryStateOnSurface::isValid(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, dqmiolumiharvest::j, k_etaBarrel, k_etaEndcap, PV3DBase< T, PVType, FrameType >::mag(), magFieldToken_, EopElecVariables::MaxPtIn01, EopElecVariables::MaxPtIn02, EopElecVariables::MaxPtIn03, EopElecVariables::MaxPtIn04, EopElecVariables::MaxPtIn05, HLTConfigProvider::moduleLabels(), multiTrajectoryStateMode::momentumFromModePPhiEta(), EopElecVariables::nHits, EopElecVariables::nLostHits, EopElecVariables::normalizedChi2, EopElecVariables::NoTrackIn0015, EopElecVariables::NoTrackIn0020, EopElecVariables::NoTrackIn0025, EopElecVariables::NoTrackIn0030, EopElecVariables::NoTrackIn0035, EopElecVariables::NoTrackIn0040, EopElecVariables::outerRadius, MultiTrajectoryStateTransform::outerStateOnSurface(), EopElecVariables::p, EopElecVariables::p_rejected_track, phi, PV3DBase< T, PVType, FrameType >::phi(), EopElecVariables::phi, EopElecVariables::phiIn, EopElecVariables::phiOut, SuperClusterShapeAlgo::phiWidth(), electrons_cff::phiWidth, EopElecVariables::pIn, reco::Vertex::position(), EopElecVariables::pOut, EopTriggerType::prescale, HLTConfigProvider::prescaleSize(), HLTConfigProvider::prescaleValue(), DiDispStaMuonMonitor_cfi::pt, EopElecVariables::pt, EopElecVariables::ptError, EopElecVariables::px, EopElecVariables::px_rejected_track, EopElecVariables::py, EopElecVariables::py_rejected_track, EopElecVariables::pz, EopElecVariables::pz_rejected_track, rho, EopElecVariables::RunNumber, EopElecVariables::SC_algoID, EopElecVariables::SC_energy, EopElecVariables::SC_eta, EopElecVariables::SC_etaWidth, EopElecVariables::SC_isBarrel, EopElecVariables::SC_isEndcap, EopElecVariables::SC_nBasicClus, EopElecVariables::SC_phi, EopElecVariables::SC_phiWidth, fileCollector::seed, AlCaHLTBitMon_QueryRunRegistry::string, EopElecVariables::SumPtIn01, EopElecVariables::SumPtIn02, EopElecVariables::SumPtIn03, EopElecVariables::SumPtIn04, EopElecVariables::SumPtIn05, theBarrelSupClusCollectionToken_, theEcalRecHitCollectionToken_, theEndCapSupClusCollectionToken_, theFilter_, theGsfElectronCoreCollectionToken_, theGsfTrackCollectionToken_, theHBHERecHitCollectionToken_, EopElecVariables::theta, theTrigger_, theTriggerEventToken_, theTriggerResultsToken_, theVertexCollectionToken_, tkGeomToken_, HLT_2024v12_cff::track, DiMuonV_cfg::tracks, tree_, treeMemPtr_, PDWG_DiPhoton_SD_cff::triggerEvent, HLTConfigProvider::triggerIndex(), PDWG_DiPhoton_SD_cff::triggerName, triggerNames_, ewkElecDQM_cfi::TrigTag, bphysicsOniaDQM_cfi::vertex, edm::HLTGlobalStatus::wasrun(), and reco::Vertex::z().

                                                                                   {
  using namespace edm;
  h_nEvents->Fill(0.5);

  Double_t EnergyHcalIn01;
  Double_t EnergyHcalIn02;
  Double_t EnergyHcalIn03;
  Double_t EnergyHcalIn04;
  Double_t EnergyHcalIn05;
  Double_t etaWidth;
  Double_t phiWidth;
  Int_t algo_ID;
  Double_t EnergyEcal;
  Double_t fbrem;
  Double_t pin;
  Double_t etaIn;
  Double_t phiIn;
  Double_t pout;
  Double_t etaOut;
  Double_t phiOut;
  Double_t MaxPtIn01;
  Double_t SumPtIn01;
  Bool_t NoTrackIn0015;
  Double_t MaxPtIn02;
  Double_t SumPtIn02;
  Bool_t NoTrackIn0020;
  Double_t MaxPtIn03;
  Double_t SumPtIn03;
  Bool_t NoTrackIn0025;
  Double_t MaxPtIn04;
  Double_t SumPtIn04;
  Bool_t NoTrackIn0030;
  Double_t MaxPtIn05;
  Double_t SumPtIn05;
  Bool_t NoTrackIn0035;
  Bool_t NoTrackIn0040;
  Double_t dRSC_first;
  Double_t dRSC_second;
  Double_t etaSC;
  Double_t phiSC;
  Int_t nbSC;        //to count the nb of SuperCluster matching with a given track
  Int_t nBasicClus;  //to count the nb of basic cluster in a given superCluster
  Bool_t isEcalDriven;
  Bool_t isTrackerDriven;
  Bool_t isBarrel;
  Bool_t isEndcap;
  Bool_t TrigTag;

  //---------------   for GsfTrack propagation through tracker  ---------------

  const MagneticField* magField_ = &iSetup.getData(magFieldToken_);
  const TrackerGeometry* trackerGeom_ = &iSetup.getData(tkGeomToken_);
  MultiTrajectoryStateTransform mtsTransform(trackerGeom_, magField_);

  //---------------    Super Cluster    -----------------

  // getting primary vertex (necessary to convert eta track to eta detector
  const reco::VertexCollection& vertex = iEvent.get(theVertexCollectionToken_);

  if (vertex.empty()) {
    edm::LogError("EopElecTreeWriter") << "Error: no primary vertex found!";
    return;
  }
  const reco::Vertex& vert = vertex.front();
  h_nEventsWithVertex->Fill(0.5);

  // getting calorimeter geometry
  const CaloGeometry* geo = &iSetup.getData(caloGeomToken_);
  const CaloSubdetectorGeometry* subGeo = geo->getSubdetectorGeometry(DetId::Ecal, EcalBarrel);
  if (subGeo == nullptr)
    edm::LogError("EopElecTreeWriter") << "ERROR: unable to find SubDetector geometry!!!";

  // getting Hcal rechits
  const HBHERecHitCollection* HcalHits = &iEvent.get(theHBHERecHitCollectionToken_);

  // getting Ecal rechits
  const EcalRecHitCollection* rhc = &iEvent.get(theEcalRecHitCollectionToken_);
  if (rhc == nullptr)
    edm::LogError("EopElecTreeWriter") << "ERROR: unable to find the EcalRecHit collection !!!";

  // getting SuperCluster
  const reco::SuperClusterCollection* BarrelSupClusCollection = &iEvent.get(theBarrelSupClusCollectionToken_);
  const reco::SuperClusterCollection* EndcapSupClusCollection = &iEvent.get(theEndCapSupClusCollectionToken_);

  // necessary to re-calculate phi and eta width of SuperClusters
  SuperClusterShapeAlgo SCShape(rhc, subGeo);

  //---------------    Trigger   -----------------
  TrigTag = false;
  const edm::TriggerResults* trigRes = &iEvent.get(theTriggerResultsToken_);

  // trigger event
  const trigger::TriggerEvent* triggerEvent = &iEvent.get(theTriggerEventToken_);

  // our trigger table
  std::map<std::string, EopTriggerType> HLTpaths;
  for (const auto& triggerName : triggerNames_) {
    if (triggerName.find(theTrigger_) != 0)
      continue;
    EopTriggerType myTrigger;

    const unsigned int prescaleSize = hltConfig_.prescaleSize();
    for (unsigned int ps = 0; ps < prescaleSize; ps++) {
      auto const prescaleValue = hltConfig_.prescaleValue<double>(ps, triggerName);
      if (prescaleValue != 1) {
        myTrigger.prescale = prescaleValue;
      }
    }

    myTrigger.index = hltConfig_.triggerIndex(triggerName);
    if (myTrigger.index == -1)
      continue;
    myTrigger.fired =
        trigRes->wasrun(myTrigger.index) && trigRes->accept(myTrigger.index) && !trigRes->error(myTrigger.index);
    HLTpaths[triggerName] = myTrigger;
  }

  // First cut : trigger cut
  std::string firstFiredPath = "";
  for (const auto& it : HLTpaths) {
    if (it.second.fired) {
      TrigTag = true;
      firstFiredPath = it.first;
      break;
    }
  }
  if (!TrigTag)
    return;
  h_nEventsTriggered->Fill(0.5);

  // Displaying filters label from the first fired trigger
  // Useful for finding the good filter label
  std::vector<std::string> filters = hltConfig_.moduleLabels(firstFiredPath);

  if (debugTriggerSelection_) {
    edm::LogInfo("EopElecTreeWriter") << "filters : ";
    for (unsigned int i = 0; i < filters.size(); i++) {
      edm::LogInfo("EopElecTreeWriter") << filters[i] << " ";
    }
  }

  // Getting HLT electrons
  edm::InputTag testTag(theFilter_, "", "HLT");
  int testindex = triggerEvent->filterIndex(testTag);

  if (testindex >= triggerEvent->sizeFilters())
    return;
  h_nEventsHLTFilter->Fill(0.5);

  const trigger::Keys& KEYS_el(triggerEvent->filterKeys(testindex));

  std::vector<const trigger::TriggerObject*> HLTelectrons;
  for (unsigned int i = 0; i != KEYS_el.size(); ++i) {
    const trigger::TriggerObject* triggerObject_el = &(triggerEvent->getObjects().at(KEYS_el[i]));
    HLTelectrons.push_back(triggerObject_el);
  }

  h_nHLTelectrons->Fill(HLTelectrons.size());

  if (HLTelectrons.empty())
    return;
  h_nEventsHLTelectron->Fill(0.5);

  // finding the HLT electron with highest pt and saving the corresponding index
  unsigned int HighPtIndex = 0;
  double maxPtHLT = -5.;
  for (unsigned int j = 0; j < HLTelectrons.size(); j++) {
    if (HLTelectrons[j]->pt() > maxPtHLT) {
      maxPtHLT = HLTelectrons[j]->pt();
      HighPtIndex = j;
    }
  }

  //-----------------   Tracks   -------------------

  // getting GsfTrack
  const reco::GsfTrackCollection& tracks = iEvent.get(theGsfTrackCollectionToken_);

  // filtering track
  int nRejected = 0;
  int nSelected = 0;
  std::vector<const reco::GsfTrack*> filterTracks;
  for (const auto& track : tracks) {
    h_nTracks->Fill(0.5);
    double deltar =
        reco::deltaR(track.eta(), track.phi(), HLTelectrons[HighPtIndex]->eta(), HLTelectrons[HighPtIndex]->phi());
    // remove the triggered electron with highest pt
    if (deltar < 0.025) {
      treeMemPtr_->px_rejected_track = track.px();
      treeMemPtr_->py_rejected_track = track.py();
      treeMemPtr_->pz_rejected_track = track.pz();
      treeMemPtr_->p_rejected_track = track.p();
      nRejected++;
      continue;
    }
    filterTracks.push_back(&track);  // we use all the others
    nSelected++;
    h_nTracksFiltered->Fill(0.5);
  }
  h_nTrkRejectedPerEvt->Fill(nRejected);
  h_nTrkSelectedPerEvt->Fill(nSelected);

  if (nRejected == 0)
    return;
  h_nEventsHLTrejected->Fill(0.5);

  if (filterTracks.empty())
    return;
  h_nEvents2Elec->Fill(0.5);

  //-------- test:Matching SC/track using gsfElectonCore collection --------

  const reco::GsfElectronCoreCollection* electrons = &iEvent.get(theGsfElectronCoreCollectionToken_);
  for (const auto& elec : *electrons) {
    double etaGSF = eopUtils::ecalEta((elec.gsfTrack())->eta(), vert.z(), (vert.position()).rho());
    if ((elec.gsfTrack())->pt() < 10.)
      continue;

    double DELTAR = 0;
    DELTAR = reco::deltaR((elec.superCluster())->eta(), (elec.superCluster())->phi(), etaGSF, (elec.gsfTrack())->phi());

    if (DELTAR < 0.1)
      h_distToClosestSCgsf->Fill(DELTAR);
  }

  //------------------------------------------------------------
  //--------------    Loop on tracks   -------------------------

  for (const auto& track : filterTracks) {
    // initializing variables
    isEcalDriven = false;
    isTrackerDriven = false;
    isBarrel = false;
    isEndcap = false;
    etaWidth = 0;
    phiWidth = 0;
    etaSC = 0;
    phiSC = 0;
    fbrem = 0;
    pin = 0;
    etaIn = 0;
    phiIn = 0;
    pout = 0;
    etaOut = 0;
    phiOut = 0;
    algo_ID = 0;
    EnergyEcal = 0;
    dRSC_first = 999;
    dRSC_second = 9999;
    nbSC = 0;
    nBasicClus = 0;

    // First cut on momentum magnitude
    h_Momentum->Fill(track->p());
    if (track->pt() < 10.)
      continue;
    h_cut_Ptmin->Fill(0.5);

    // calculating track parameters at innermost and outermost for Gsf tracks
    TrajectoryStateOnSurface inTSOS = mtsTransform.innerStateOnSurface(*track);
    TrajectoryStateOnSurface outTSOS = mtsTransform.outerStateOnSurface(*track);
    if (inTSOS.isValid() && outTSOS.isValid()) {
      GlobalVector inMom;
      //multiTrajectoryStateMode::momentumFromModeCartesian(inTSOS, inMom_);
      multiTrajectoryStateMode::momentumFromModePPhiEta(inTSOS, inMom);
      pin = inMom.mag();
      etaIn = inMom.eta();
      phiIn = inMom.phi();
      GlobalVector outMom;
      //multiTrajectoryStateMode::momentumFromModeCartesian(outTSOS, outMom_);
      multiTrajectoryStateMode::momentumFromModePPhiEta(outTSOS, outMom);
      pout = outMom.mag();
      etaOut = outMom.eta();
      phiOut = outMom.phi();
      fbrem = (pin - pout) / pin;
      h_fBREM->Fill(fbrem);
    }

    // Matching track with Hcal rec hits
    EnergyHcalIn01 = 0;
    EnergyHcalIn02 = 0;
    EnergyHcalIn03 = 0;
    EnergyHcalIn04 = 0;
    EnergyHcalIn05 = 0;

    //for (std::vector<HBHERecHit>::const_iterator hcal = (*HcalHits).begin(); hcal != (*HcalHits).end(); hcal++) {
    for (const auto& hcal : *HcalHits) {
      GlobalPoint posH = geo->getPosition(hcal.detid());
      double phihit = posH.phi();
      double etahit = posH.eta();
      double dR = reco::deltaR(etahit, phihit, etaOut, phiOut);

      // saving Hcal energy deposit measured for different eta-phi radius
      if (dR < 0.1)
        EnergyHcalIn01 += hcal.energy();
      if (dR < 0.2)
        EnergyHcalIn02 += hcal.energy();
      if (dR < 0.3)
        EnergyHcalIn03 += hcal.energy();
      if (dR < 0.4)
        EnergyHcalIn04 += hcal.energy();
      if (dR < 0.5)
        EnergyHcalIn05 += hcal.energy();
    }

    HcalEnergy->Fill(EnergyHcalIn02);

    //Isolation against charged particles
    MaxPtIn01 = 0.;
    SumPtIn01 = 0.;
    NoTrackIn0015 = true;
    MaxPtIn02 = 0.;
    SumPtIn02 = 0.;
    NoTrackIn0020 = true;
    MaxPtIn03 = 0.;
    SumPtIn03 = 0.;
    NoTrackIn0025 = true;
    MaxPtIn04 = 0.;
    SumPtIn04 = 0.;
    NoTrackIn0030 = true;
    MaxPtIn05 = 0.;
    SumPtIn05 = 0.;
    NoTrackIn0035 = true;
    NoTrackIn0040 = true;

    for (const auto& track1 : filterTracks) {
      if (track == track1)
        continue;

      double etaIn1 = 0.;
      double phiIn1 = 0.;

      TrajectoryStateOnSurface inTSOS1 = mtsTransform.innerStateOnSurface(*track1);
      if (inTSOS1.isValid()) {
        GlobalVector inMom1;
        multiTrajectoryStateMode::momentumFromModePPhiEta(inTSOS1, inMom1);
        etaIn1 = inMom1.eta();
        phiIn1 = inMom1.phi();
      }

      if (etaIn1 == 0 && phiIn1 == 0)
        continue;

      double dR = reco::deltaR(etaIn1, phiIn1, etaIn, phiIn);

      // different radius of inner isolation cone
      if (dR < 0.015)
        NoTrackIn0015 = false;
      if (dR < 0.020)
        NoTrackIn0020 = false;
      if (dR < 0.025)
        NoTrackIn0025 = false;
      if (dR < 0.030)
        NoTrackIn0030 = false;
      if (dR < 0.035)
        NoTrackIn0035 = false;
      if (dR < 0.040)
        NoTrackIn0040 = false;

      //calculate maximum Pt and sum Pt inside cones of different radius
      if (dR < 0.1) {
        SumPtIn01 += track1->pt();
        if (track1->pt() > MaxPtIn01) {
          MaxPtIn01 = track1->pt();
        }
      }

      if (dR < 0.2) {
        SumPtIn02 += track1->pt();
        if (track1->pt() > MaxPtIn02) {
          MaxPtIn02 = track1->pt();
        }
      }
      if (dR < 0.3) {
        SumPtIn03 += track1->pt();
        if (track1->pt() > MaxPtIn03) {
          MaxPtIn03 = track1->pt();
        }
      }
      if (dR < 0.4) {
        SumPtIn04 += track1->pt();
        if (track1->pt() > MaxPtIn04) {
          MaxPtIn04 = track1->pt();
        }
      }
      if (dR < 0.5) {
        SumPtIn05 += track1->pt();
        if (track1->pt() > MaxPtIn05) {
          MaxPtIn05 = track1->pt();
        }
      }
    }

    // Track-SuperCluster matching

    double dRSC;
    double etaAtEcal = 0;  // to convert eta from track to detector frame
    etaAtEcal = eopUtils::ecalEta(etaIn, vert.z(), (vert.position()).rho());

    //Barrel
    if (std::abs(track->eta()) < k_etaBarrel) {
      for (const auto& SC : *BarrelSupClusCollection) {
        dRSC = reco::deltaR(SC.eta(), SC.phi(), etaAtEcal, phiIn);

        if (dRSC < dRSC_first) {
          dRSC_first = dRSC;
        }  // distance in eta-phi plan to closest SC
        else if (dRSC < dRSC_second) {
          dRSC_second = dRSC;
        }  // to next closest SC

        if (dRSC < 0.09) {
          //Calculate phiWidth & etaWidth for associated SuperClusters
          SCShape.Calculate_Covariances(SC);
          phiWidth = SCShape.phiWidth();
          etaWidth = SCShape.etaWidth();

          etaSC = SC.eta();
          phiSC = SC.phi();

          algo_ID = SC.algoID();
          EnergyEcal = SC.energy();
          nBasicClus = SC.clustersSize();
          nbSC++;
          isBarrel = true;
        }
      }
    }

    // Endcap
    if (std::abs(track->eta()) > k_etaEndcap) {
      for (const auto& SC : *EndcapSupClusCollection) {
        dRSC = reco::deltaR(SC.eta(), SC.phi(), etaAtEcal, phiIn);

        if (dRSC < dRSC_first) {
          dRSC_first = dRSC;
        }  // distance in eta-phi plan to closest SC
        else if (dRSC < dRSC_second) {
          dRSC_second = dRSC;
        }  // to next closest SC

        if (dRSC < 0.09) {
          //Calculate phiWidth & etaWidth for associated SuperClusters
          SCShape.Calculate_Covariances(SC);
          phiWidth = SCShape.phiWidth();
          etaWidth = SCShape.etaWidth();

          etaSC = SC.eta();
          phiSC = SC.phi();

          algo_ID = SC.algoID();
          EnergyEcal = SC.energy();
          nBasicClus = SC.clustersSize();
          nbSC++;
          isEndcap = true;
        }
      }
    }

    if (dRSC_first < 0.1)
      h_distToClosestSC->Fill(dRSC_first);
    if (nbSC == 1)
      h_counter1->Fill(0.5);
    if (nbSC == 0)
      h_counter2->Fill(0.5);

    if (nbSC > 1 || nbSC == 0)
      continue;
    h_cut_OneSCmatch->Fill(0.5);

    if (isBarrel && isEndcap) {
      edm::LogError("EopElecTreeWriter") << "Error: Super Cluster double matching!";
      return;
    }

    h_EcalEnergy->Fill(EnergyEcal);
    HcalVSEcal->Fill(EnergyEcal, EnergyHcalIn02);

    // E over p plots
    if (track->charge() < 0)
      h_Eop_InnerNegative->Fill(EnergyEcal / pin);
    if (track->charge() > 0)
      h_Eop_InnerPositive->Fill(EnergyEcal / pin);

    //Check if track-SuperCluster matching is Ecal or Tracker driven
    edm::RefToBase<TrajectorySeed> seed = track->extra()->seedRef();
    if (seed.isNull()) {
      edm::LogError("GsfElectronCore") << "The GsfTrack has no seed ?!";
    } else {
      reco::ElectronSeedRef elseed = seed.castTo<reco::ElectronSeedRef>();
      if (elseed.isNull()) {
        edm::LogError("GsfElectronCore") << "The GsfTrack seed is not an ElectronSeed ?!";
      } else {
        isEcalDriven = elseed->isEcalDriven();
        isTrackerDriven = elseed->isTrackerDriven();
      }
    }

    treeMemPtr_->charge = track->charge();
    treeMemPtr_->nHits = track->numberOfValidHits();
    treeMemPtr_->nLostHits = track->numberOfLostHits();
    treeMemPtr_->innerOk = track->innerOk();
    treeMemPtr_->outerRadius = track->outerRadius();
    treeMemPtr_->chi2 = track->chi2();
    treeMemPtr_->normalizedChi2 = track->normalizedChi2();
    treeMemPtr_->px = track->px();
    treeMemPtr_->py = track->py();
    treeMemPtr_->pz = track->pz();
    treeMemPtr_->p = track->p();
    treeMemPtr_->pIn = pin;
    treeMemPtr_->etaIn = etaIn;
    treeMemPtr_->phiIn = phiIn;
    treeMemPtr_->pOut = pout;
    treeMemPtr_->etaOut = etaOut;
    treeMemPtr_->phiOut = phiOut;
    treeMemPtr_->pt = track->pt();
    treeMemPtr_->ptError = track->ptError();
    treeMemPtr_->theta = track->theta();
    treeMemPtr_->eta = track->eta();
    treeMemPtr_->phi = track->phi();
    treeMemPtr_->fbrem = fbrem;
    treeMemPtr_->MaxPtIn01 = MaxPtIn01;
    treeMemPtr_->SumPtIn01 = SumPtIn01;
    treeMemPtr_->NoTrackIn0015 = NoTrackIn0015;
    treeMemPtr_->MaxPtIn02 = MaxPtIn02;
    treeMemPtr_->SumPtIn02 = SumPtIn02;
    treeMemPtr_->NoTrackIn0020 = NoTrackIn0020;
    treeMemPtr_->MaxPtIn03 = MaxPtIn03;
    treeMemPtr_->SumPtIn03 = SumPtIn03;
    treeMemPtr_->NoTrackIn0025 = NoTrackIn0025;
    treeMemPtr_->MaxPtIn04 = MaxPtIn04;
    treeMemPtr_->SumPtIn04 = SumPtIn04;
    treeMemPtr_->NoTrackIn0030 = NoTrackIn0030;
    treeMemPtr_->MaxPtIn05 = MaxPtIn05;
    treeMemPtr_->SumPtIn05 = SumPtIn05;
    treeMemPtr_->NoTrackIn0035 = NoTrackIn0035;
    treeMemPtr_->NoTrackIn0040 = NoTrackIn0040;
    treeMemPtr_->SC_algoID = algo_ID;
    treeMemPtr_->SC_energy = EnergyEcal;
    treeMemPtr_->SC_nBasicClus = nBasicClus;
    treeMemPtr_->SC_etaWidth = etaWidth;
    treeMemPtr_->SC_phiWidth = phiWidth;
    treeMemPtr_->SC_eta = etaSC;
    treeMemPtr_->SC_phi = phiSC;
    treeMemPtr_->SC_isBarrel = isBarrel;
    treeMemPtr_->SC_isEndcap = isEndcap;
    treeMemPtr_->dRto1stSC = dRSC_first;
    treeMemPtr_->dRto2ndSC = dRSC_second;
    treeMemPtr_->HcalEnergyIn01 = EnergyHcalIn01;
    treeMemPtr_->HcalEnergyIn02 = EnergyHcalIn02;
    treeMemPtr_->HcalEnergyIn03 = EnergyHcalIn03;
    treeMemPtr_->HcalEnergyIn04 = EnergyHcalIn04;
    treeMemPtr_->HcalEnergyIn05 = EnergyHcalIn05;
    treeMemPtr_->isEcalDriven = isEcalDriven;
    treeMemPtr_->isTrackerDriven = isTrackerDriven;
    treeMemPtr_->RunNumber = iEvent.id().run();
    treeMemPtr_->EvtNumber = iEvent.id().event();

    tree_->Fill();

  }  // loop on tracks
}  // analyze

beginRun()

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

overrideprivate

Definition at line 881 of file EopElecTreeWriter.cc.

References debugTriggerSelection_, hltConfig_, mps_fire::i, HLTConfigProvider::init(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, HLTConfigProvider::triggerNames(), and triggerNames_.

                                                                                {
  bool changed = true;

  // Load trigger configuration at each begin of run
  // Fill the trigger names
  if (hltConfig_.init(iRun, iSetup, "HLT", changed)) {
    if (changed) {
      triggerNames_ = hltConfig_.triggerNames();
    }
  }

  // Displaying the trigger names
  if (debugTriggerSelection_) {
    unsigned int i = 0;
    for (const auto& it : triggerNames_) {
      edm::LogInfo("EopElecTreeWriter") << "HLTpath: " << (i++) << " = " << it;
    }
  }
}

endJob()

void EopElecTreeWriter::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 875 of file EopElecTreeWriter.cc.

References treeMemPtr_.

Referenced by o2olib.O2ORunMgr::executeJob().

                               {
  delete treeMemPtr_;
  treeMemPtr_ = nullptr;
}

endRun()

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

inlineoverrideprivate

Definition at line 111 of file EopElecTreeWriter.cc.

{};

fillDescriptions()

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

static

Definition at line 902 of file EopElecTreeWriter.cc.

References edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, ProducerED_cfi::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

{
  edm::ParameterSetDescription desc;
  desc.setComment("Generate tree for Tracker Alignment E/p validation");
  desc.add<edm::InputTag>("src", edm::InputTag("electronGsfTracks"));
  desc.add<std::string>("triggerPath", "HLT_Ele");
  desc.add<std::string>("hltFilter", "hltDiEle27L1DoubleEGWPTightHcalIsoFilter");
  desc.add<bool>("debugTriggerSelection", false);
  descriptions.addWithDefaultLabel(desc);
}

Member Data Documentation

caloGeomToken_

const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> EopElecTreeWriter::caloGeomToken_

private

Definition at line 119 of file EopElecTreeWriter.cc.

Referenced by analyze().

debugTriggerSelection_

bool EopElecTreeWriter::debugTriggerSelection_

private

Definition at line 139 of file EopElecTreeWriter.cc.

Referenced by analyze(), and beginRun().

fs_

edm::Service<TFileService> EopElecTreeWriter::fs_

private

Definition at line 140 of file EopElecTreeWriter.cc.

Referenced by EopElecTreeWriter().

h_counter1

TH1D* EopElecTreeWriter::h_counter1

private

Definition at line 167 of file EopElecTreeWriter.cc.

Referenced by analyze(), and EopElecTreeWriter().

h_counter2

TH1D* EopElecTreeWriter::h_counter2

private

Definition at line 168 of file EopElecTreeWriter.cc.

Referenced by analyze(), and EopElecTreeWriter().

h_cut_OneSCmatch

TH1D* EopElecTreeWriter::h_cut_OneSCmatch

private

Definition at line 165 of file EopElecTreeWriter.cc.

Referenced by analyze(), and EopElecTreeWriter().

h_cut_Ptmin

TH1D* EopElecTreeWriter::h_cut_Ptmin

private

Definition at line 164 of file EopElecTreeWriter.cc.

Referenced by analyze(), and EopElecTreeWriter().

h_distToClosestSC

TH1D* EopElecTreeWriter::h_distToClosestSC

private

Definition at line 171 of file EopElecTreeWriter.cc.

Referenced by analyze(), EopElecTreeWriter(), and ~EopElecTreeWriter().

h_distToClosestSCgsf

TH1D* EopElecTreeWriter::h_distToClosestSCgsf

private

Definition at line 170 of file EopElecTreeWriter.cc.

Referenced by analyze(), EopElecTreeWriter(), and ~EopElecTreeWriter().

h_EcalEnergy

TH1D* EopElecTreeWriter::h_EcalEnergy

private

Definition at line 172 of file EopElecTreeWriter.cc.

Referenced by analyze(), EopElecTreeWriter(), and ~EopElecTreeWriter().

h_Eop_InnerNegative

TH1D* EopElecTreeWriter::h_Eop_InnerNegative

private

Definition at line 176 of file EopElecTreeWriter.cc.

Referenced by analyze(), EopElecTreeWriter(), and ~EopElecTreeWriter().

h_Eop_InnerPositive

TH1D* EopElecTreeWriter::h_Eop_InnerPositive

private

Definition at line 177 of file EopElecTreeWriter.cc.

Referenced by analyze(), EopElecTreeWriter(), and ~EopElecTreeWriter().

h_fBREM

TH1D* EopElecTreeWriter::h_fBREM

private

Definition at line 175 of file EopElecTreeWriter.cc.

Referenced by analyze(), and EopElecTreeWriter().

h_Momentum

TH1D* EopElecTreeWriter::h_Momentum

private

Definition at line 173 of file EopElecTreeWriter.cc.

Referenced by analyze(), EopElecTreeWriter(), and ~EopElecTreeWriter().

h_nEvents

TH1D* EopElecTreeWriter::h_nEvents

private

Definition at line 150 of file EopElecTreeWriter.cc.

Referenced by analyze(), and EopElecTreeWriter().

h_nEvents2Elec

TH1D* EopElecTreeWriter::h_nEvents2Elec

private

Definition at line 156 of file EopElecTreeWriter.cc.

Referenced by analyze(), and EopElecTreeWriter().

h_nEventsHLTelectron

TH1D* EopElecTreeWriter::h_nEventsHLTelectron

private

Definition at line 154 of file EopElecTreeWriter.cc.

Referenced by analyze(), and EopElecTreeWriter().

h_nEventsHLTFilter

TH1D* EopElecTreeWriter::h_nEventsHLTFilter

private

Definition at line 153 of file EopElecTreeWriter.cc.

Referenced by analyze(), and EopElecTreeWriter().

h_nEventsHLTrejected

TH1D* EopElecTreeWriter::h_nEventsHLTrejected

private

Definition at line 155 of file EopElecTreeWriter.cc.

Referenced by analyze(), and EopElecTreeWriter().

h_nEventsTriggered

TH1D* EopElecTreeWriter::h_nEventsTriggered

private

Definition at line 152 of file EopElecTreeWriter.cc.

Referenced by analyze(), and EopElecTreeWriter().

h_nEventsWithVertex

TH1D* EopElecTreeWriter::h_nEventsWithVertex

private

Definition at line 151 of file EopElecTreeWriter.cc.

Referenced by analyze(), and EopElecTreeWriter().

h_nHLTelectrons

TH1D* EopElecTreeWriter::h_nHLTelectrons

private

Definition at line 157 of file EopElecTreeWriter.cc.

Referenced by analyze(), and EopElecTreeWriter().

h_nTracks

TH1D* EopElecTreeWriter::h_nTracks

private

Definition at line 162 of file EopElecTreeWriter.cc.

Referenced by analyze(), and EopElecTreeWriter().

h_nTracksFiltered

TH1D* EopElecTreeWriter::h_nTracksFiltered

private

Definition at line 163 of file EopElecTreeWriter.cc.

Referenced by analyze(), and EopElecTreeWriter().

h_nTrkRejectedPerEvt

TH1D* EopElecTreeWriter::h_nTrkRejectedPerEvt

private

Definition at line 158 of file EopElecTreeWriter.cc.

Referenced by analyze(), and EopElecTreeWriter().

h_nTrkSelectedPerEvt

TH1D* EopElecTreeWriter::h_nTrkSelectedPerEvt

private

Definition at line 159 of file EopElecTreeWriter.cc.

Referenced by analyze(), and EopElecTreeWriter().

HcalEnergy

TH1D* EopElecTreeWriter::HcalEnergy

private

Definition at line 174 of file EopElecTreeWriter.cc.

Referenced by analyze(), EopElecTreeWriter(), and ~EopElecTreeWriter().

HcalVSEcal

TH2D* EopElecTreeWriter::HcalVSEcal

private

Definition at line 179 of file EopElecTreeWriter.cc.

Referenced by analyze(), EopElecTreeWriter(), and ~EopElecTreeWriter().

hltConfig_

HLTConfigProvider EopElecTreeWriter::hltConfig_

private

Definition at line 143 of file EopElecTreeWriter.cc.

Referenced by analyze(), and beginRun().

k_etaBarrel

constexpr float EopElecTreeWriter::k_etaBarrel = 1.55

staticprivate

Definition at line 133 of file EopElecTreeWriter.cc.

Referenced by analyze().

k_etaEndcap

constexpr float EopElecTreeWriter::k_etaEndcap = 1.44

staticprivate

Definition at line 134 of file EopElecTreeWriter.cc.

Referenced by analyze().

magFieldToken_

const edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> EopElecTreeWriter::magFieldToken_

private

Definition at line 117 of file EopElecTreeWriter.cc.

Referenced by analyze().

mtsTransform_

MultiTrajectoryStateTransform* EopElecTreeWriter::mtsTransform_

private

Definition at line 145 of file EopElecTreeWriter.cc.

theBarrelSupClusCollectionToken_

const edm::EDGetTokenT<reco::SuperClusterCollection> EopElecTreeWriter::theBarrelSupClusCollectionToken_

private

Definition at line 125 of file EopElecTreeWriter.cc.

Referenced by analyze().

theEcalRecHitCollectionToken_

const edm::EDGetTokenT<EcalRecHitCollection> EopElecTreeWriter::theEcalRecHitCollectionToken_

private

Definition at line 124 of file EopElecTreeWriter.cc.

Referenced by analyze().

theEndCapSupClusCollectionToken_

const edm::EDGetTokenT<reco::SuperClusterCollection> EopElecTreeWriter::theEndCapSupClusCollectionToken_

private

Definition at line 126 of file EopElecTreeWriter.cc.

Referenced by analyze().

theFilter_

std::string EopElecTreeWriter::theFilter_

private

Definition at line 138 of file EopElecTreeWriter.cc.

Referenced by analyze().

theGsfElectronCoreCollectionToken_

const edm::EDGetTokenT<reco::GsfElectronCoreCollection> EopElecTreeWriter::theGsfElectronCoreCollectionToken_

private

Definition at line 130 of file EopElecTreeWriter.cc.

Referenced by analyze().

theGsfTrackCollectionToken_

const edm::EDGetTokenT<reco::GsfTrackCollection> EopElecTreeWriter::theGsfTrackCollectionToken_

private

Definition at line 129 of file EopElecTreeWriter.cc.

Referenced by analyze().

theHBHERecHitCollectionToken_

const edm::EDGetTokenT<HBHERecHitCollection> EopElecTreeWriter::theHBHERecHitCollectionToken_

private

Definition at line 123 of file EopElecTreeWriter.cc.

Referenced by analyze().

theTrigger_

std::string EopElecTreeWriter::theTrigger_

private

Definition at line 137 of file EopElecTreeWriter.cc.

Referenced by analyze().

theTriggerEventToken_

const edm::EDGetTokenT<trigger::TriggerEvent> EopElecTreeWriter::theTriggerEventToken_

private

Definition at line 128 of file EopElecTreeWriter.cc.

Referenced by analyze().

theTriggerResultsToken_

const edm::EDGetTokenT<edm::TriggerResults> EopElecTreeWriter::theTriggerResultsToken_

private

Definition at line 127 of file EopElecTreeWriter.cc.

Referenced by analyze().

theVertexCollectionToken_

const edm::EDGetTokenT<reco::VertexCollection> EopElecTreeWriter::theVertexCollectionToken_

private

Definition at line 122 of file EopElecTreeWriter.cc.

Referenced by analyze().

tkGeomToken_

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> EopElecTreeWriter::tkGeomToken_

private

Definition at line 118 of file EopElecTreeWriter.cc.

Referenced by analyze().

tree_

TTree* EopElecTreeWriter::tree_

private

Definition at line 141 of file EopElecTreeWriter.cc.

Referenced by analyze(), and EopElecTreeWriter().

treeMemPtr_

EopElecVariables* EopElecTreeWriter::treeMemPtr_

private

Definition at line 142 of file EopElecTreeWriter.cc.

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

triggerNames_

std::vector<std::string> EopElecTreeWriter::triggerNames_

private

Definition at line 144 of file EopElecTreeWriter.cc.

Referenced by analyze(), and beginRun().