EnergyScaleAnalyzer Class Reference

#include <EnergyScaleAnalyzer.h>

Inheritance diagram for EnergyScaleAnalyzer:

Classes
struct	tree_structure_

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	beginJob () override
void	endJob () override
	EnergyScaleAnalyzer (const edm::ParameterSet &)
	~EnergyScaleAnalyzer () override
Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default
	EDAnalyzer (const EDAnalyzer &)=delete
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
const EDAnalyzer &	operator= (const EDAnalyzer &)=delete
bool	wantsGlobalLuminosityBlocks () const noexcept final
bool	wantsGlobalRuns () const noexcept final
bool	wantsInputProcessBlocks () const noexcept final
bool	wantsProcessBlocks () const noexcept 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 noexcept
bool	wantsStreamRuns () const noexcept
	~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)

Private Member Functions
void	fillTree (const reco::SuperClusterCollection *scColl, const reco::SuperClusterCollection *corrSCColl, HepMC::GenParticle *mc, tree_structure_ &tree_, float xV, float yV, float zV, int scType)

Private Attributes
edm::EDGetTokenT< reco::SuperClusterCollection >	correctedDynamicHybridSuperClusters_token
edm::EDGetTokenT< reco::SuperClusterCollection >	correctedFixedMatrixSuperClustersWithPreshower_token
edm::EDGetTokenT< reco::SuperClusterCollection >	correctedHybridSuperClusters_token
edm::EDGetTokenT< reco::SuperClusterCollection >	dynamicHybridSuperClusters_token
float	energyMax_
float	etaMax_
float	etaMaxVtx_
float	eTMax_
float	eTMaxVtx_
int	evtN
edm::EDGetTokenT< reco::SuperClusterCollection >	fixedMatrixSuperClustersWithPreshower_token
edm::EDGetTokenT< edm::HepMCProduct >	hepMCLabel_
edm::EDGetTokenT< reco::SuperClusterCollection >	hybridSuperClusters_token
TTree *	mytree_
std::string	outputFile_
float	phiMax_
float	phiMaxVtx_
float	rClust_vtx_
TFile *	rootFile_
float	thetaMax_
float	thetaMaxVtx_
tree_structure_	tree_
float	xClust_vtx_
float	xClust_zero_
float	xVtx_
float	yClust_vtx_
float	yClust_zero_
float	yVtx_
float	zClust_vtx_
float	zClust_zero_
float	zVtx_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
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

Validation/EcalClusters/src/EnergyScaleAnalyzer.cc

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 48 of file EnergyScaleAnalyzer.h.

Constructor & Destructor Documentation

EnergyScaleAnalyzer()

EnergyScaleAnalyzer::EnergyScaleAnalyzer ( const edm::ParameterSet &  ps )

explicit

Definition at line 70 of file EnergyScaleAnalyzer.cc.

References correctedDynamicHybridSuperClusters_token, correctedFixedMatrixSuperClustersWithPreshower_token, correctedHybridSuperClusters_token, dynamicHybridSuperClusters_token, evtN, fixedMatrixSuperClustersWithPreshower_token, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hepMCLabel_, hybridSuperClusters_token, ProducerED_cfi::InputTag, outputFile_, rootFile_, and AlCaHLTBitMon_QueryRunRegistry::string.

{
  hepMCLabel_ = consumes<edm::HepMCProduct>(ps.getParameter<std::string>("hepMCLabel"));
  hybridSuperClusters_token = consumes<reco::SuperClusterCollection>(
      ps.getUntrackedParameter<edm::InputTag>("hybridSuperClusters", edm::InputTag("hybridSuperClusters")));
  dynamicHybridSuperClusters_token = consumes<reco::SuperClusterCollection>(ps.getUntrackedParameter<edm::InputTag>(
      "dynamicHybridSuperClusters", edm::InputTag("dynamicHybridSuperClusters")));
  correctedHybridSuperClusters_token = consumes<reco::SuperClusterCollection>(ps.getUntrackedParameter<edm::InputTag>(
      "correctedHybridSuperClusters", edm::InputTag("correctedHybridSuperClusters")));
  correctedDynamicHybridSuperClusters_token =
      consumes<reco::SuperClusterCollection>(ps.getUntrackedParameter<edm::InputTag>(
          "correctedDynamicHybridSuperClusters", edm::InputTag("correctedDynamicHybridSuperClusters")));
  correctedFixedMatrixSuperClustersWithPreshower_token = consumes<reco::SuperClusterCollection>(
      ps.getUntrackedParameter<edm::InputTag>("correctedFixedMatrixSuperClustersWithPreshower",
                                              edm::InputTag("correctedFixedMatrixSuperClustersWithPreshower")));
  fixedMatrixSuperClustersWithPreshower_token =
      consumes<reco::SuperClusterCollection>(ps.getUntrackedParameter<edm::InputTag>(
          "fixedMatrixSuperClustersWithPreshower", edm::InputTag("fixedMatrixSuperClustersWithPreshower")));

  outputFile_ = ps.getParameter<std::string>("outputFile");
  rootFile_ = TFile::Open(outputFile_.c_str(),
                          "RECREATE");  // open output file to store histograms

  evtN = 0;
}

~EnergyScaleAnalyzer()

EnergyScaleAnalyzer::~EnergyScaleAnalyzer ( )

override

Definition at line 98 of file EnergyScaleAnalyzer.cc.

References rootFile_.

{
  delete rootFile_;
}

Member Function Documentation

analyze()

void EnergyScaleAnalyzer::analyze ( const edm::Event &  evt, const edm::EventSetup &  es  )

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 131 of file EnergyScaleAnalyzer.cc.

References funct::abs(), correctedDynamicHybridSuperClusters_token, correctedFixedMatrixSuperClustersWithPreshower_token, correctedHybridSuperClusters_token, kinem::delta_R(), dynamicHybridSuperClusters_cfi::dynamicHybridSuperClusters, dynamicHybridSuperClusters_token, Exception, fillTree(), fixedMatrixSuperClustersWithPreshower_token, genWeightsTable_cfi::genEvent, GenParticle::GenParticle, edm::Event::getByToken(), MuonClassifier_cff::hepMC, hepMCLabel_, unifiedSCCollection_cfi::hybridSuperClusters, hybridSuperClusters_token, mps_fire::i, sistrip::SpyUtilities::isValid(), MainPageGenerator::l, edm::EDConsumerBase::labelsForToken(), HLT_2024v14_cff::max_eta, CaloTowersParam_cfi::mc, EnergyScaleAnalyzer::tree_structure_::mc_e, EnergyScaleAnalyzer::tree_structure_::mc_et, EnergyScaleAnalyzer::tree_structure_::mc_eta, EnergyScaleAnalyzer::tree_structure_::mc_npar, EnergyScaleAnalyzer::tree_structure_::mc_phi, EnergyScaleAnalyzer::tree_structure_::mc_sep, EnergyScaleAnalyzer::tree_structure_::mc_theta, HerwigMaxPtPartonFilter_cfi::moduleLabel, AlCaHLTBitMon_ParallelJobs::p, EnergyScaleAnalyzer::tree_structure_::parID, edm::Handle< T >::product(), crab_script::provenance, funct::sin(), AlCaHLTBitMon_QueryRunRegistry::string, tree_, L1BJetProducer_cff::vtx, xVtx_, yVtx_, and zVtx_.

                                                                              {
  using namespace edm;  // needed for all fwk related classes
  using namespace std;

  //  std::cout << "Proceccing event # " << ++evtN << std::endl;

  // Get containers for MC truth, SC etc.
  // ===================================================
  // =======================================================================================
  // =======================================================================================
  Handle<HepMCProduct> hepMC;
  evt.getByToken(hepMCLabel_, hepMC);

  Labels l;
  labelsForToken(hepMCLabel_, l);

  const HepMC::GenEvent *genEvent = hepMC->GetEvent();
  if (!(hepMC.isValid())) {
    LogInfo("EnergyScaleAnalyzer") << "Could not get MC Product!";
    return;
  }

  //=======================For Vertex correction
  std::vector<Handle<HepMCProduct>> evtHandles;

  //evt.getManyByType(evtHandles);
  throw cms::Exception("UnsupportedFunction") << "EnergyScaleAnalyzer::analyze: "
                                              << "getManyByType has not been supported by the Framework since 2015. "
                                              << "This module has been broken since then. Maybe it should be deleted. "
                                              << "Another possibility is to upgrade to use GetterOfProducts instead.";

  for (unsigned int i = 0; i < evtHandles.size(); ++i) {
    if (evtHandles[i].isValid()) {
      const HepMC::GenEvent *evt = evtHandles[i]->GetEvent();

      // take only 1st vertex for now - it's been tested only of PGuns...
      //
      HepMC::GenEvent::vertex_const_iterator vtx = evt->vertices_begin();
      if (evtHandles[i].provenance()->moduleLabel() == std::string(l.module)) {
        // Corrdinates of Vertex w.r.o. the point (0,0,0)
        xVtx_ = 0.1 * (*vtx)->position().x();
        yVtx_ = 0.1 * (*vtx)->position().y();
        zVtx_ = 0.1 * (*vtx)->position().z();
      }
    }
  }
  //==============================================================================
  // Get handle to SC collections

  Handle<reco::SuperClusterCollection> hybridSuperClusters;
  try {
    evt.getByToken(hybridSuperClusters_token, hybridSuperClusters);
  } catch (cms::Exception &ex) {
    edm::LogError("EnergyScaleAnalyzer") << "Can't get collection with producer hybridSuperClusters.";
  }

  Handle<reco::SuperClusterCollection> dynamicHybridSuperClusters;
  try {
    evt.getByToken(dynamicHybridSuperClusters_token, dynamicHybridSuperClusters);
  } catch (cms::Exception &ex) {
    edm::LogError("EnergyScaleAnalyzer") << "Can't get collection with producer dynamicHybridSuperClusters.";
  }

  Handle<reco::SuperClusterCollection> fixedMatrixSuperClustersWithPS;
  try {
    evt.getByToken(fixedMatrixSuperClustersWithPreshower_token, fixedMatrixSuperClustersWithPS);
  } catch (cms::Exception &ex) {
    edm::LogError("EnergyScaleAnalyzer") << "Can't get collection with producer "
                                            "fixedMatrixSuperClustersWithPreshower.";
  }

  // Corrected collections
  Handle<reco::SuperClusterCollection> correctedHybridSC;
  try {
    evt.getByToken(correctedHybridSuperClusters_token, correctedHybridSC);
  } catch (cms::Exception &ex) {
    edm::LogError("EnergyScaleAnalyzer") << "Can't get collection with producer correctedHybridSuperClusters.";
  }

  Handle<reco::SuperClusterCollection> correctedDynamicHybridSC;
  try {
    evt.getByToken(correctedDynamicHybridSuperClusters_token, correctedDynamicHybridSC);
  } catch (cms::Exception &ex) {
    edm::LogError("EnergyScaleAnalyzer") << "Can't get collection with producer "
                                            "correctedDynamicHybridSuperClusters.";
  }

  Handle<reco::SuperClusterCollection> correctedFixedMatrixSCWithPS;
  try {
    evt.getByToken(correctedFixedMatrixSuperClustersWithPreshower_token, correctedFixedMatrixSCWithPS);
  } catch (cms::Exception &ex) {
    edm::LogError("EnergyScaleAnalyzer") << "Can't get collection with producer "
                                            "correctedFixedMatrixSuperClustersWithPreshower.";
  }

  const reco::SuperClusterCollection *dSC = dynamicHybridSuperClusters.product();
  const reco::SuperClusterCollection *hSC = hybridSuperClusters.product();
  const reco::SuperClusterCollection *fmSC = fixedMatrixSuperClustersWithPS.product();
  const reco::SuperClusterCollection *chSC = correctedHybridSC.product();
  const reco::SuperClusterCollection *cdSC = correctedDynamicHybridSC.product();
  const reco::SuperClusterCollection *cfmSC = correctedFixedMatrixSCWithPS.product();

  // -----------------------  Print outs for debugging
  /*
  std::cout << "MC truth" << std::endl;
  int counterI = 0;
  for(  HepMC::GenEvent::particle_const_iterator p =
  genEvent->particles_begin(); p != genEvent->particles_end(); ++p ) { if (
  fabs((*p)->momentum().eta()) < 1.5 ) { std::cout << ++counterI << " " <<
  (*p)->momentum().e() << " "
                << (*p)->momentum().phi() << " " << (*p)->momentum().eta() <<
  std::endl;
    }
  }

  std::cout << "Standard clusters" << std::endl;
  counterI = 0;
  for(reco::SuperClusterCollection::const_iterator em = hSC->begin();
      em != hSC->end(); ++em )
    std::cout << ++counterI << " " << em->energy() << " " <<
  em->position().phi() << " " << em->position().eta() << std::endl;

  std::cout << "Dynamic clusters" << std::endl;
  counterI = 0;
  for(reco::SuperClusterCollection::const_iterator em = dSC->begin();
      em != dSC->end(); ++em )
    std::cout << ++counterI << " " << em->energy() << " " <<
  em->position().phi() << " " << em->position().eta() << std::endl;

  std::cout << "FixedMatrix clusters with PS" << std::endl;
  counterI = 0;
  for(reco::SuperClusterCollection::const_iterator em = fmSC->begin();
      em != fmSC->end(); ++em )
    std::cout << ++counterI << " " << em->energy() << " " <<
  em->position().phi() << " " << em->position().eta() << std::endl;
  */
  // -----------------------------
  //=====================================================================
  // All containers are loaded, perform the analysis
  //====================================================================

  // --------------------------- Store MC particles
  HepMC::GenEvent::particle_const_iterator p = genEvent->particles_begin();

  // Search for MC electrons or photons that satisfy the criteria
  float min_eT = 5;
  float max_eta = 2.5;

  std::vector<HepMC::GenParticle *> mcParticles;
  // int counter = 0;
  for (HepMC::GenEvent::particle_const_iterator p = genEvent->particles_begin(); p != genEvent->particles_end(); ++p) {
    // LogInfo("EnergyScaleAnalyzer") << "Particle " << ++counter
    //<< " PDG ID = " << (*p)->pdg_id() << " pT = " << (*p)->momentum().perp();
    // require photon or electron
    if ((*p)->pdg_id() != 22 && abs((*p)->pdg_id()) != 11)
      continue;

    // require selection criteria
    bool satisfySelectionCriteria = (*p)->momentum().perp() > min_eT && fabs((*p)->momentum().eta()) < max_eta;

    if (!satisfySelectionCriteria)
      continue;

    // EM MC particle is found, save it in the vector
    mcParticles.push_back(*p);
  }
  // separation in dR between 2 first MC particles
  // should not be used for MC samples with > 2 em objects generated!
  if (mcParticles.size() == 2) {
    HepMC::GenParticle *mc1 = mcParticles[0];
    HepMC::GenParticle *mc2 = mcParticles[1];
    tree_.mc_sep =
        kinem::delta_R(mc1->momentum().eta(), mc1->momentum().phi(), mc2->momentum().eta(), mc2->momentum().phi());
  } else
    tree_.mc_sep = -100;

  // now loop over MC particles, find the match with SC and do everything we
  // need then save info in the tree for every MC particle
  for (std::vector<HepMC::GenParticle *>::const_iterator p = mcParticles.begin(); p != mcParticles.end(); ++p) {
    HepMC::GenParticle *mc = *p;

    // Fill MC information
    tree_.mc_npar = mcParticles.size();
    tree_.parID = mc->pdg_id();
    tree_.mc_e = mc->momentum().e();
    tree_.mc_et = mc->momentum().e() * sin(mc->momentum().theta());
    tree_.mc_phi = mc->momentum().phi();
    tree_.mc_eta = mc->momentum().eta();
    tree_.mc_theta = mc->momentum().theta();

    // Call function to fill tree
    // scType coprreponds:
    // HybridSuperCluster                     -- 1
    // DynamicHybridSuperCluster              -- 2
    // FixedMatrixSuperClustersWithPreshower  -- 3

    fillTree(hSC, chSC, mc, tree_, xVtx_, yVtx_, zVtx_, 1);
    //    std::cout << " TYPE " << 1 << " : " << tree_.em_e << " : " <<
    //    tree_.em_phi << " : " << tree_.em_eta << std::endl;

    fillTree(dSC, cdSC, mc, tree_, xVtx_, yVtx_, zVtx_, 2);
    //    std::cout << " TYPE " << 2 << " : " << tree_.em_e << " : " <<
    //    tree_.em_phi << " : " << tree_.em_eta << std::endl;

    fillTree(fmSC, cfmSC, mc, tree_, xVtx_, yVtx_, zVtx_, 3);
    //    std::cout << " TYPE " << 3 << " : " << tree_.em_e << " : " <<
    //    tree_.em_phi << " : " << tree_.em_eta << std::endl;

    //   mytree_->Fill();
  }  // loop over particles
}

beginJob()

void EnergyScaleAnalyzer::beginJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 105 of file EnergyScaleAnalyzer.cc.

References EnergyScaleAnalyzer::tree_structure_::mc_npar, mytree_, and tree_.

                                   {
  //========================================================================

  mytree_ = new TTree("energyScale", "");
  TString treeVariables =
      "mc_npar/I:parID:mc_sep/"
      "F:mc_e:mc_et:mc_phi:mc_eta:mc_theta:";  // MC information
  treeVariables += "em_dR/F:";                 // MC <-> EM matching information
  treeVariables +=
      "em_isInCrack/I:em_scType:em_e/F:em_et:em_phi:em_eta:em_theta:em_nCell/"
      "I:em_nBC:";                                       // EM SC info
  treeVariables += "em_pet/F:em_pe:em_peta:em_ptheta:";  // EM SC physics (eta corrected
                                                         // information)

  treeVariables += "emCorr_e/F:emCorr_et:emCorr_eta:emCorr_phi:emCorr_theta:";  // CMSSW
                                                                                // standard
                                                                                // corrections
  treeVariables += "emCorr_pet/F:emCorr_peta:emCorr_ptheta:";                   // CMSSW standard physics

  treeVariables += "em_pw/F:em_ew:em_br";  // EM widths pw -- phiWidth, ew --
                                           // etaWidth, ratios of pw/ew

  mytree_->Branch("energyScale", &(tree_.mc_npar), treeVariables);
}

endJob()

void EnergyScaleAnalyzer::endJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 473 of file EnergyScaleAnalyzer.cc.

References rootFile_.

                                 {
  //========================================================================
  // Fill ROOT tree
  rootFile_->Write();
}

fillTree()

void EnergyScaleAnalyzer::fillTree ( const reco::SuperClusterCollection *  scColl, const reco::SuperClusterCollection *  corrSCColl, HepMC::GenParticle *  mc, tree_structure_ &  tree_, float  xV, float  yV, float  zV, int  scType  )

private

fill tree with kinematic variables of corrected Super Cluster

Definition at line 343 of file EnergyScaleAnalyzer.cc.

References kinem::delta_R(), HGC3DClusterGenMatchSelector_cfi::dR, EnergyScaleAnalyzer::tree_structure_::em_br, EnergyScaleAnalyzer::tree_structure_::em_dR, EnergyScaleAnalyzer::tree_structure_::em_e, EnergyScaleAnalyzer::tree_structure_::em_et, EnergyScaleAnalyzer::tree_structure_::em_eta, EnergyScaleAnalyzer::tree_structure_::em_ew, EnergyScaleAnalyzer::tree_structure_::em_isInCrack, EnergyScaleAnalyzer::tree_structure_::em_nBC, EnergyScaleAnalyzer::tree_structure_::em_nCell, EnergyScaleAnalyzer::tree_structure_::em_pe, EnergyScaleAnalyzer::tree_structure_::em_pet, EnergyScaleAnalyzer::tree_structure_::em_peta, EnergyScaleAnalyzer::tree_structure_::em_phi, EnergyScaleAnalyzer::tree_structure_::em_ptheta, EnergyScaleAnalyzer::tree_structure_::em_pw, EnergyScaleAnalyzer::tree_structure_::em_scType, EnergyScaleAnalyzer::tree_structure_::em_theta, EnergyScaleAnalyzer::tree_structure_::emCorr_e, EnergyScaleAnalyzer::tree_structure_::emCorr_et, EnergyScaleAnalyzer::tree_structure_::emCorr_eta, EnergyScaleAnalyzer::tree_structure_::emCorr_pet, EnergyScaleAnalyzer::tree_structure_::emCorr_peta, EnergyScaleAnalyzer::tree_structure_::emCorr_phi, EnergyScaleAnalyzer::tree_structure_::emCorr_ptheta, EnergyScaleAnalyzer::tree_structure_::emCorr_theta, hcalRecHitTable_cff::energy, energyMax_, etaMax_, etaMaxVtx_, eTMax_, eTMaxVtx_, dqm-mbProfile::log, CaloTowersParam_cfi::mc, mytree_, phiMax_, phiMaxVtx_, rClust_vtx_, funct::sin(), mathSSE::sqrt(), funct::tan(), thetaMax_, thetaMaxVtx_, tree_, xClust_vtx_, xClust_zero_, yClust_vtx_, yClust_zero_, zClust_vtx_, and zClust_zero_.

Referenced by analyze().

                                               {
  // -----------------------------  SuperClusters before energy correction
  reco::SuperClusterCollection::const_iterator em = scColl->end();
  float energyMax = -100.0;  // dummy energy of the matched SC
  for (reco::SuperClusterCollection::const_iterator aClus = scColl->begin(); aClus != scColl->end(); ++aClus) {
    // check the matching
    float dR =
        kinem::delta_R(mc->momentum().eta(), mc->momentum().phi(), aClus->position().eta(), aClus->position().phi());
    if (dR < 0.4) {  // a rather loose matching cut
      float energy = aClus->energy();
      if (energy < energyMax)
        continue;
      energyMax = energy;
      em = aClus;
    }
  }

  if (em == scColl->end()) {
    //    std::cout << "No matching SC with type " << scType << " was found for
    //    MC particle! "  << std::endl; std::cout << "Going to next type of SC.
    //    " << std::endl;
    return;
  }
  //  ------------

  tree_.em_scType = scType;

  tree_.em_isInCrack = 0;
  double emAbsEta = fabs(em->position().eta());
  // copied from
  // RecoEgama/EgammaElectronAlgos/src/EgammaElectronClassification.cc
  if (emAbsEta < 0.018 || (emAbsEta > 0.423 && emAbsEta < 0.461) || (emAbsEta > 0.770 && emAbsEta < 0.806) ||
      (emAbsEta > 1.127 && emAbsEta < 1.163) || (emAbsEta > 1.460 && emAbsEta < 1.558))
    tree_.em_isInCrack = 1;

  tree_.em_dR = kinem::delta_R(mc->momentum().eta(), mc->momentum().phi(), em->position().eta(), em->position().phi());
  tree_.em_e = em->energy();
  tree_.em_et = em->energy() * sin(em->position().theta());
  tree_.em_phi = em->position().phi();
  tree_.em_eta = em->position().eta();
  tree_.em_theta = em->position().theta();
  tree_.em_nCell = em->size();
  tree_.em_nBC = em->clustersSize();

  // Get physics e, et etc:
  // Coordinates of EM object with respect of the point (0,0,0)
  xClust_zero_ = em->position().x();
  yClust_zero_ = em->position().y();
  zClust_zero_ = em->position().z();
  // Coordinates of EM object w.r.o. the Vertex position
  xClust_vtx_ = xClust_zero_ - xV;
  yClust_vtx_ = yClust_zero_ - yV;
  zClust_vtx_ = zClust_zero_ - zV;

  energyMax_ = em->energy();
  thetaMax_ = em->position().theta();
  etaMax_ = em->position().eta();
  phiMax_ = em->position().phi();
  eTMax_ = energyMax_ * sin(thetaMax_);
  if (phiMax_ < 0)
    phiMax_ += 2 * 3.14159;

  rClust_vtx_ = sqrt(xClust_vtx_ * xClust_vtx_ + yClust_vtx_ * yClust_vtx_ + zClust_vtx_ * zClust_vtx_);
  thetaMaxVtx_ = acos(zClust_vtx_ / rClust_vtx_);
  etaMaxVtx_ = -log(tan(thetaMaxVtx_ / 2));
  eTMaxVtx_ = energyMax_ * sin(thetaMaxVtx_);
  phiMaxVtx_ = atan2(yClust_vtx_, xClust_vtx_);
  if (phiMaxVtx_ < 0)
    phiMaxVtx_ += 2 * 3.14159;
  //=============================
  // parametres of EM object after vertex correction
  tree_.em_pet = eTMaxVtx_;
  tree_.em_pe = tree_.em_pet / sin(thetaMaxVtx_);
  tree_.em_peta = etaMaxVtx_;
  tree_.em_ptheta = thetaMaxVtx_;

  //-------------------------------   Get SC after energy correction
  em = corrSCColl->end();
  energyMax = -100.0;  // dummy energy of the matched SC
  for (reco::SuperClusterCollection::const_iterator aClus = corrSCColl->begin(); aClus != corrSCColl->end(); ++aClus) {
    // check the matching
    float dR =
        kinem::delta_R(mc->momentum().eta(), mc->momentum().phi(), aClus->position().eta(), aClus->position().phi());
    if (dR < 0.4) {
      float energy = aClus->energy();
      if (energy < energyMax)
        continue;
      energyMax = energy;
      em = aClus;
    }
  }

  if (em == corrSCColl->end()) {
    //    std::cout << "No matching corrected SC with type " << scType << " was
    //    found for MC particle! "  << std::endl; std::cout << "Going to next
    //    type of SC. " << std::endl;
    return;
  }
  //  ------------

  tree_.emCorr_e = em->energy();
  tree_.emCorr_et = em->energy() * sin(em->position().theta());
  tree_.emCorr_phi = em->position().phi();
  tree_.emCorr_eta = em->position().eta();
  tree_.emCorr_theta = em->position().theta();

  // =========== Eta and Theta wrt Vertex does not change after energy
  // corrections are applied
  // =========== So, no need to calculate them again

  tree_.emCorr_peta = tree_.em_peta;
  tree_.emCorr_ptheta = tree_.em_ptheta;
  tree_.emCorr_pet = tree_.emCorr_e / cosh(tree_.emCorr_peta);

  tree_.em_pw = em->phiWidth();
  tree_.em_ew = em->etaWidth();
  tree_.em_br = tree_.em_pw / tree_.em_ew;

  mytree_->Fill();
}

Member Data Documentation

correctedDynamicHybridSuperClusters_token

edm::EDGetTokenT<reco::SuperClusterCollection> EnergyScaleAnalyzer::correctedDynamicHybridSuperClusters_token

private

Definition at line 64 of file EnergyScaleAnalyzer.h.

Referenced by analyze(), and EnergyScaleAnalyzer().

correctedFixedMatrixSuperClustersWithPreshower_token

edm::EDGetTokenT<reco::SuperClusterCollection> EnergyScaleAnalyzer::correctedFixedMatrixSuperClustersWithPreshower_token

private

Definition at line 65 of file EnergyScaleAnalyzer.h.

Referenced by analyze(), and EnergyScaleAnalyzer().

correctedHybridSuperClusters_token

edm::EDGetTokenT<reco::SuperClusterCollection> EnergyScaleAnalyzer::correctedHybridSuperClusters_token

private

Definition at line 63 of file EnergyScaleAnalyzer.h.

Referenced by analyze(), and EnergyScaleAnalyzer().

dynamicHybridSuperClusters_token

edm::EDGetTokenT<reco::SuperClusterCollection> EnergyScaleAnalyzer::dynamicHybridSuperClusters_token

private

Definition at line 62 of file EnergyScaleAnalyzer.h.

Referenced by analyze(), and EnergyScaleAnalyzer().

energyMax_

float EnergyScaleAnalyzer::energyMax_

private

Definition at line 137 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

etaMax_

float EnergyScaleAnalyzer::etaMax_

private

Definition at line 140 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

etaMaxVtx_

float EnergyScaleAnalyzer::etaMaxVtx_

private

Definition at line 141 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

eTMax_

float EnergyScaleAnalyzer::eTMax_

private

Definition at line 138 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

eTMaxVtx_

float EnergyScaleAnalyzer::eTMaxVtx_

private

Definition at line 139 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

evtN

int EnergyScaleAnalyzer::evtN

private

Definition at line 148 of file EnergyScaleAnalyzer.h.

Referenced by EnergyScaleAnalyzer().

fixedMatrixSuperClustersWithPreshower_token

edm::EDGetTokenT<reco::SuperClusterCollection> EnergyScaleAnalyzer::fixedMatrixSuperClustersWithPreshower_token

private

Definition at line 66 of file EnergyScaleAnalyzer.h.

Referenced by analyze(), and EnergyScaleAnalyzer().

hepMCLabel_

edm::EDGetTokenT<edm::HepMCProduct> EnergyScaleAnalyzer::hepMCLabel_

private

Definition at line 60 of file EnergyScaleAnalyzer.h.

Referenced by analyze(), and EnergyScaleAnalyzer().

hybridSuperClusters_token

edm::EDGetTokenT<reco::SuperClusterCollection> EnergyScaleAnalyzer::hybridSuperClusters_token

private

Definition at line 61 of file EnergyScaleAnalyzer.h.

Referenced by analyze(), and EnergyScaleAnalyzer().

mytree_

TTree* EnergyScaleAnalyzer::mytree_

private

Definition at line 72 of file EnergyScaleAnalyzer.h.

Referenced by beginJob(), and fillTree().

outputFile_

std::string EnergyScaleAnalyzer::outputFile_

private

Definition at line 58 of file EnergyScaleAnalyzer.h.

Referenced by EnergyScaleAnalyzer().

phiMax_

float EnergyScaleAnalyzer::phiMax_

private

Definition at line 142 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

phiMaxVtx_

float EnergyScaleAnalyzer::phiMaxVtx_

private

Definition at line 143 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

rClust_vtx_

float EnergyScaleAnalyzer::rClust_vtx_

private

Definition at line 135 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

rootFile_

TFile* EnergyScaleAnalyzer::rootFile_

private

Definition at line 69 of file EnergyScaleAnalyzer.h.

Referenced by endJob(), EnergyScaleAnalyzer(), and ~EnergyScaleAnalyzer().

thetaMax_

float EnergyScaleAnalyzer::thetaMax_

private

Definition at line 144 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

thetaMaxVtx_

float EnergyScaleAnalyzer::thetaMaxVtx_

private

Definition at line 145 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

tree_

tree_structure_ EnergyScaleAnalyzer::tree_

private

Definition at line 120 of file EnergyScaleAnalyzer.h.

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

xClust_vtx_

float EnergyScaleAnalyzer::xClust_vtx_

private

Definition at line 131 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

xClust_zero_

float EnergyScaleAnalyzer::xClust_zero_

private

Definition at line 127 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

xVtx_

float EnergyScaleAnalyzer::xVtx_

private

Definition at line 123 of file EnergyScaleAnalyzer.h.

Referenced by analyze().

yClust_vtx_

float EnergyScaleAnalyzer::yClust_vtx_

private

Definition at line 132 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

yClust_zero_

float EnergyScaleAnalyzer::yClust_zero_

private

Definition at line 128 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

yVtx_

float EnergyScaleAnalyzer::yVtx_

private

Definition at line 124 of file EnergyScaleAnalyzer.h.

Referenced by analyze().

zClust_vtx_

float EnergyScaleAnalyzer::zClust_vtx_

private

Definition at line 133 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

zClust_zero_

float EnergyScaleAnalyzer::zClust_zero_

private

Definition at line 129 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

zVtx_

float EnergyScaleAnalyzer::zVtx_

private

Definition at line 125 of file EnergyScaleAnalyzer.h.

Referenced by analyze().