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 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
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > 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
bool	registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices 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 ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
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 ( 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, HLT_FULL_cff::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 ( )

override

Definition at line 98 of file EnergyScaleAnalyzer.cc.

References rootFile_.

{
  delete rootFile_;
}

Member Function Documentation

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, fillTree(), fixedMatrixSuperClustersWithPreshower_token, MCTruth::genEvent, GenParticle::GenParticle, edm::Event::getByToken(), edm::Event::getManyByType(), hepMCLabel_, unifiedSCCollection_cfi::hybridSuperClusters, hybridSuperClusters_token, mps_fire::i, sistrip::SpyUtilities::isValid(), cmsLHEtoEOSManager::l, edm::EDConsumerBase::labelsForToken(), HLT_FULL_cff::max_eta, 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, AlCaHLTBitMon_ParallelJobs::p, EnergyScaleAnalyzer::tree_structure_::parID, funct::sin(), AlCaHLTBitMon_QueryRunRegistry::string, tree_, 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);

  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
}

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);
}

void EnergyScaleAnalyzer::endJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 468 of file EnergyScaleAnalyzer.cc.

References rootFile_.

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

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 338 of file EnergyScaleAnalyzer.cc.

References kinem::delta_R(), 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, relval_parameters_module::energy, energyMax_, etaMax_, etaMaxVtx_, eTMax_, eTMaxVtx_, log, mytree_, phiMax_, phiMaxVtx_, rClust_vtx_, funct::sin(), mathSSE::sqrt(), funct::tan(), thetaMax_, thetaMaxVtx_, 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

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

private

Definition at line 64 of file EnergyScaleAnalyzer.h.

Referenced by analyze(), and EnergyScaleAnalyzer().

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

private

Definition at line 65 of file EnergyScaleAnalyzer.h.

Referenced by analyze(), and EnergyScaleAnalyzer().

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

private

Definition at line 63 of file EnergyScaleAnalyzer.h.

Referenced by analyze(), and EnergyScaleAnalyzer().

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

private

Definition at line 62 of file EnergyScaleAnalyzer.h.

Referenced by analyze(), and EnergyScaleAnalyzer().

float EnergyScaleAnalyzer::energyMax_

private

Definition at line 137 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

float EnergyScaleAnalyzer::etaMax_

private

Definition at line 140 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

float EnergyScaleAnalyzer::etaMaxVtx_

private

Definition at line 141 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

float EnergyScaleAnalyzer::eTMax_

private

Definition at line 138 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

float EnergyScaleAnalyzer::eTMaxVtx_

private

Definition at line 139 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

int EnergyScaleAnalyzer::evtN

private

Definition at line 148 of file EnergyScaleAnalyzer.h.

Referenced by EnergyScaleAnalyzer().

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

private

Definition at line 66 of file EnergyScaleAnalyzer.h.

Referenced by analyze(), and EnergyScaleAnalyzer().

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

private

Definition at line 60 of file EnergyScaleAnalyzer.h.

Referenced by analyze(), and EnergyScaleAnalyzer().

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

private

Definition at line 61 of file EnergyScaleAnalyzer.h.

Referenced by analyze(), and EnergyScaleAnalyzer().

TTree* EnergyScaleAnalyzer::mytree_

private

Definition at line 72 of file EnergyScaleAnalyzer.h.

Referenced by beginJob(), and fillTree().

std::string EnergyScaleAnalyzer::outputFile_

private

Definition at line 58 of file EnergyScaleAnalyzer.h.

Referenced by EnergyScaleAnalyzer().

float EnergyScaleAnalyzer::phiMax_

private

Definition at line 142 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

float EnergyScaleAnalyzer::phiMaxVtx_

private

Definition at line 143 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

float EnergyScaleAnalyzer::rClust_vtx_

private

Definition at line 135 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

TFile* EnergyScaleAnalyzer::rootFile_

private

Definition at line 69 of file EnergyScaleAnalyzer.h.

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

float EnergyScaleAnalyzer::thetaMax_

private

Definition at line 144 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

float EnergyScaleAnalyzer::thetaMaxVtx_

private

Definition at line 145 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

tree_structure_ EnergyScaleAnalyzer::tree_

private

Definition at line 120 of file EnergyScaleAnalyzer.h.

Referenced by analyze(), and beginJob().

float EnergyScaleAnalyzer::xClust_vtx_

private

Definition at line 131 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

float EnergyScaleAnalyzer::xClust_zero_

private

Definition at line 127 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

float EnergyScaleAnalyzer::xVtx_

private

Definition at line 123 of file EnergyScaleAnalyzer.h.

Referenced by analyze().

float EnergyScaleAnalyzer::yClust_vtx_

private

Definition at line 132 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

float EnergyScaleAnalyzer::yClust_zero_

private

Definition at line 128 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

float EnergyScaleAnalyzer::yVtx_

private

Definition at line 124 of file EnergyScaleAnalyzer.h.

Referenced by analyze().

float EnergyScaleAnalyzer::zClust_vtx_

private

Definition at line 133 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

float EnergyScaleAnalyzer::zClust_zero_

private

Definition at line 129 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

float EnergyScaleAnalyzer::zVtx_

private

Definition at line 125 of file EnergyScaleAnalyzer.h.

Referenced by analyze().