EcalSimHitStudy Class Reference

Inheritance diagram for EcalSimHitStudy:

Classes
struct	EcalHit

Public Member Functions
	EcalSimHitStudy (const edm::ParameterSet &ps)
	~EcalSimHitStudy () 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
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)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions
void	analyze (edm::Event const &, edm::EventSetup const &) override
void	analyzeHits (std::vector< PCaloHit > &, int)
void	beginJob () override
void	beginRun (edm::Run const &, edm::EventSetup const &) override
void	endRun (edm::Run const &, edm::EventSetup const &) override
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)

Private Attributes
TH1F *	edep_ [ndets_]
TH1F *	edepAll_ [ndets_]
TH1F *	edepEM_ [ndets_]
TH1F *	edepHad_ [ndets_]
TH1F *	eneInc_
TH1F *	etaInc_
TH1F *	etot_ [ndets_]
TH1F *	etotg_ [ndets_]
const std::string	g4Label_
const CaloGeometry *	geometry_
TH1F *	hit_ [ndets_]
const std::vector< std::string >	hitLab_
const double	maxEnergy_
TH1F *	phiInc_
TH2F *	poszn_ [ndets_]
TH2F *	poszp_ [ndets_]
TH1F *	ptInc_
TH1F *	r1by25_ [ndets_]
TH1F *	r1by9_ [ndets_]
TH1F *	r9by25_ [ndets_]
const int	selX_
TH1F *	sEtaEta_ [ndets_]
TH1F *	sEtaPhi_ [ndets_]
TH1F *	sPhiPhi_ [ndets_]
TH1F *	time_ [ndets_]
TH1F *	timeAll_ [ndets_]
const double	tmax_
const edm::EDGetTokenT< edm::HepMCProduct >	tok_evt_
const edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	tokGeom_
const std::vector< edm::EDGetTokenT< edm::PCaloHitContainer > >	toks_calo_
const double	w0_

Static Private Attributes
static const int	ndets_ = 2

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

Detailed Description

Definition at line 46 of file EcalSimHitStudy.cc.

Constructor & Destructor Documentation

EcalSimHitStudy()

EcalSimHitStudy::EcalSimHitStudy

(

const edm::ParameterSet &

ps

)

Definition at line 83 of file EcalSimHitStudy.cc.

References g4Label_, hitLab_, Skims_PA_cff::name, AlCaHLTBitMon_QueryRunRegistry::string, and edm::vector_transform().

    : tokGeom_(esConsumes<CaloGeometry, CaloGeometryRecord>()),
      g4Label_(ps.getUntrackedParameter<std::string>("ModuleLabel", "g4SimHits")),
      hitLab_(ps.getUntrackedParameter<std::vector<std::string>>("CaloCollection")),
      maxEnergy_(ps.getUntrackedParameter<double>("MaxEnergy", 200.0)),
      tmax_(ps.getUntrackedParameter<double>("TimeCut", 100.0)),
      w0_(ps.getUntrackedParameter<double>("W0", 4.7)),
      selX_(ps.getUntrackedParameter<int>("SelectX", -1)),
      tok_evt_(consumes<edm::HepMCProduct>(
          edm::InputTag(ps.getUntrackedParameter<std::string>("SourceLabel", "VtxSmeared")))),
      toks_calo_{edm::vector_transform(hitLab_, [this](const std::string& name) {
        return consumes<edm::PCaloHitContainer>(edm::InputTag{g4Label_, name});
      })} {
  usesResource(TFileService::kSharedResource);

  edm::LogVerbatim("HitStudy") << "Module Label: " << g4Label_ << "   Hits: " << hitLab_[0] << ", " << hitLab_[1]
                               << "   MaxEnergy: " << maxEnergy_ << "  Tmax: " << tmax_ << " Select " << selX_;
}

~EcalSimHitStudy()

EcalSimHitStudy::~EcalSimHitStudy ( )

inlineoverride

Definition at line 49 of file EcalSimHitStudy.cc.

{}

Member Function Documentation

analyze()

void EcalSimHitStudy::analyze ( edm::Event const &  e, edm::EventSetup const &  iS  )

overrideprotectedvirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 256 of file EcalSimHitStudy.cc.

References funct::abs(), analyzeHits(), MillePedeFileConverter_cfg::e, eneInc_, etaInc_, geometry_, edm::EventSetup::getData(), edm::HepMCProduct::GetEvent(), edm::HandleBase::isValid(), AlCaHLTBitMon_ParallelJobs::p, phiInc_, ptInc_, tok_evt_, tokGeom_, and toks_calo_.

                                                                        {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HitStudy") << "Run = " << e.id().run() << " Event = " << e.id().event();
#endif
  // get handles to calogeometry
  geometry_ = &iS.getData(tokGeom_);

  double eInc = 0, etaInc = 0, phiInc = 0;
  int type(-1);
  edm::Handle<edm::HepMCProduct> EvtHandle;
  e.getByToken(tok_evt_, EvtHandle);
  if (EvtHandle.isValid()) {
    const HepMC::GenEvent* myGenEvent = EvtHandle->GetEvent();

    HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin();
    if (p != myGenEvent->particles_end()) {
      eInc = (*p)->momentum().e();
      etaInc = (*p)->momentum().eta();
      phiInc = (*p)->momentum().phi();
    }
    double ptInc = eInc / std::cosh(etaInc);
    ptInc_->Fill(ptInc);
    eneInc_->Fill(eInc);
    etaInc_->Fill(etaInc);
    phiInc_->Fill(phiInc);

    if (std::abs(etaInc) < 1.46)
      type = 0;
    else if (std::abs(etaInc) > 1.49 && std::abs(etaInc) < 3.0)
      type = 1;
  }

  int typeMin = (type < 0) ? 0 : type;
  int typeMax = (type < 0) ? 1 : type;
  for (int type = typeMin; type <= typeMax; ++type) {
    edm::Handle<edm::PCaloHitContainer> hitsCalo;
    e.getByToken(toks_calo_[type], hitsCalo);
    bool getHits = (hitsCalo.isValid());
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HitStudy") << "EcalSimHitStudy: Input flags Hits " << getHits << " with " << hitsCalo->size()
                                 << " hits";
#endif
    if (getHits) {
      std::vector<PCaloHit> caloHits;
      caloHits.insert(caloHits.end(), hitsCalo->begin(), hitsCalo->end());
      if (!caloHits.empty())
        analyzeHits(caloHits, type);
    }
  }
}

analyzeHits()

void EcalSimHitStudy::analyzeHits ( std::vector< PCaloHit > &  hits, int  indx  )

protected

Definition at line 307 of file EcalSimHitStudy.cc.

References funct::abs(), TauDecayModes::dec, makePileupJSON::denom, HLT_2022v11_cff::dEta, HLT_2022v11_cff::dPhi, PVValHelper::dz, edep_, edepAll_, edepEM_, edepHad_, HCALHighEnergyHPDFilter_cfi::energy, etot_, etotg_, geometry_, CaloGeometry::getGeometry(), hit_, hfClusterShapes_cfi::hits, triggerObjects_cff::id, hit::id, LEDCalibrationChannels::iphi, dqmdumpme::k, dqm-mbProfile::log, M_PI, SiStripPI::max, position, poszn_, poszp_, r1by25_, r1by9_, r9by25_, singleTopDQM_cfi::select, selX_, sEtaEta_, sEtaPhi_, sPhiPhi_, protons_cff::time, time_, timeAll_, tmax_, w(), w0_, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), PV3DBase< T, PVType, FrameType >::z(), and ecaldqm::zside().

Referenced by analyze().

                                                                     {
  unsigned int nEC(0);
  std::map<unsigned int, EcalHit> hitMap;
  double etot(0), etotG(0);
  for (auto hit : hits) {
    double edep = hit.energy();
    double time = hit.time();
    unsigned int id_ = hit.id();
    double edepEM = hit.energyEM();
    double edepHad = hit.energyHad();
    if (indx == 0) {
      if ((hit.depth() == 1) || (hit.depth() == 2))
        continue;
    }
    if (time <= tmax_) {
      auto it = hitMap.find(id_);
      if (it == hitMap.end()) {
        uint16_t dep = hit.depth();
        hitMap[id_] = EcalHit(dep, time, edep);
      } else {
        (it->second).energy += edep;
      }
      etotG += edep;
      ++nEC;
    }
    time_[indx]->Fill(time);
    edep_[indx]->Fill(edep);
    edepEM_[indx]->Fill(edepEM);
    edepHad_[indx]->Fill(edepHad);
    etot += edep;
  }

  double edepM(0);
  unsigned int idM(0);
  uint16_t depM(0);
  for (auto it : hitMap) {
    if (it.second.energy > edepM) {
      idM = it.first;
      edepM = it.second.energy;
      depM = it.second.id;
    }
  }

  bool select(true);
  if (selX_ >= 0) {
    if ((depM & 0X4) != 0)
      select = (selX_ > 0);
    else
      select = (selX_ == 0);
  }
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HitStudy") << "EcalSimHitStudy::analyzeHits: Index " << indx << " Emax " << edepM << " IDMax "
                               << std::hex << idM << ":" << depM << std::dec << " Select " << select << ":" << selX_
                               << " Hits " << hits.size() << ":" << nEC << ":" << hitMap.size() << " ETotal " << etot
                               << ":" << etotG;
#endif
  if (select) {
    etot_[indx]->Fill(etot);
    etotg_[indx]->Fill(etotG);
    hit_[indx]->Fill(double(nEC));
    for (auto it : hitMap) {
      timeAll_[indx]->Fill((it.second).time);
      edepAll_[indx]->Fill((it.second).energy);
      DetId id(it.first);
      if (indx == 0) {
        if (EBDetId(id).zside() >= 0)
          poszp_[indx]->Fill(EBDetId(id).ietaAbs(), EBDetId(id).iphi());
        else
          poszn_[indx]->Fill(EBDetId(id).ietaAbs(), EBDetId(id).iphi());
      } else {
        GlobalPoint gpos = geometry_->getGeometry(id)->getPosition();
        if (EEDetId(id).zside() >= 0)
          poszp_[indx]->Fill(gpos.x(), gpos.y());
        else
          poszn_[indx]->Fill(gpos.x(), gpos.y());
      }
    }

    math::XYZVector meanPosition(0.0, 0.0, 0.0);
    std::vector<math::XYZVector> position;
    std::vector<double> energy;
    double e9(0), e25(0);
    for (auto it : hitMap) {
      DetId id(it.first);
      int deta(99), dphi(99), dz(0);
      if (indx == 0) {
        deta = std::abs(EBDetId(id).ietaAbs() - EBDetId(idM).ietaAbs());
        dphi = std::abs(EBDetId(id).iphi() - EBDetId(idM).iphi());
        if (dphi > 180)
          dphi = std::abs(dphi - 360);
        dz = std::abs(EBDetId(id).zside() - EBDetId(idM).zside());
      } else {
        deta = std::abs(EEDetId(id).ix() - EEDetId(idM).ix());
        dphi = std::abs(EEDetId(id).iy() - EEDetId(idM).iy());
        dz = std::abs(EEDetId(id).zside() - EEDetId(idM).zside());
      }
      if (deta <= 1 && dphi <= 1 && dz < 1)
        e9 += (it.second).energy;
      if (deta <= 2 && dphi <= 2 && dz < 1) {
        e25 += (it.second).energy;
        GlobalPoint gpos = geometry_->getGeometry(id)->getPosition();
        math::XYZVector pos(gpos.x(), gpos.y(), gpos.z());
        meanPosition += (it.second).energy * pos;
        position.push_back(pos);
        energy.push_back((it.second).energy);
      }
    }
    double r1by9 = (e9 > 0) ? (edepM / e9) : -1;
    double r1by25 = (e25 > 0) ? (edepM / e25) : -1;
    double r9by25 = (e25 > 0) ? (e9 / e25) : -1;

    meanPosition /= e25;
    double denom(0), numEtaEta(0), numEtaPhi(0), numPhiPhi(0);
    for (unsigned int k = 0; k < position.size(); ++k) {
      double dEta = position[k].eta() - meanPosition.eta();
      double dPhi = position[k].phi() - meanPosition.phi();
      if (dPhi > +M_PI) {
        dPhi = 2 * M_PI - dPhi;
      }
      if (dPhi < -M_PI) {
        dPhi = 2 * M_PI + dPhi;
      }

      double w = std::max(0.0, (w0_ + std::log(energy[k] / e25)));
      denom += w;
      numEtaEta += std::abs(w * dEta * dEta);
      numEtaPhi += std::abs(w * dEta * dPhi);
      numPhiPhi += std::abs(w * dPhi * dPhi);
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("HitStudy") << "[" << k << "] dEta " << dEta << " dPhi " << dPhi << " Wt " << energy[k] / e25
                                   << ":" << std::log(energy[k] / e25) << ":" << w;
#endif
    }
    double sEtaEta = (denom > 0) ? (numEtaEta / denom) : -1.0;
    double sEtaPhi = (denom > 0) ? (numEtaPhi / denom) : -1.0;
    double sPhiPhi = (denom > 0) ? (numPhiPhi / denom) : -1.0;

#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HitStudy") << "EcalSimHitStudy::Ratios " << r1by9 << " : " << r1by25 << " : " << r9by25
                                 << " Covariances " << sEtaEta << " : " << sEtaPhi << " : " << sPhiPhi;
#endif
    r1by9_[indx]->Fill(r1by9);
    r1by25_[indx]->Fill(r1by25);
    r9by25_[indx]->Fill(r9by25);
    sEtaEta_[indx]->Fill(sEtaEta);
    sEtaPhi_[indx]->Fill(sEtaPhi);
    sPhiPhi_[indx]->Fill(sPhiPhi);
  }
}

beginJob()

void EcalSimHitStudy::beginJob ( void  )

overrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 114 of file EcalSimHitStudy.cc.

References edep_, edepAll_, edepEM_, edepHad_, eneInc_, etaInc_, etot_, etotg_, hit_, mps_fire::i, maxEnergy_, Skims_PA_cff::name, ndets_, phiInc_, poszn_, poszp_, ptInc_, r1by25_, r1by9_, r9by25_, sEtaEta_, sEtaPhi_, sPhiPhi_, AlCaHLTBitMon_QueryRunRegistry::string, compare::tfile, time_, timeAll_, runGCPTkAlMap::title, and L1TOccupancyClient_cfi::ymax.

                               {
  edm::Service<TFileService> tfile;
  if (!tfile.isAvailable())
    throw cms::Exception("BadConfig") << "TFileService unavailable: "
                                      << "please add it to config file";
  char name[20], title[120];
  sprintf(title, "Incident PT (GeV)");
  ptInc_ = tfile->make<TH1F>("PtInc", title, 1000, 0., maxEnergy_);
  ptInc_->GetXaxis()->SetTitle(title);
  ptInc_->GetYaxis()->SetTitle("Events");
  sprintf(title, "Incident Energy (GeV)");
  eneInc_ = tfile->make<TH1F>("EneInc", title, 1000, 0., maxEnergy_);
  eneInc_->GetXaxis()->SetTitle(title);
  eneInc_->GetYaxis()->SetTitle("Events");
  sprintf(title, "Incident #eta");
  etaInc_ = tfile->make<TH1F>("EtaInc", title, 200, -5., 5.);
  etaInc_->GetXaxis()->SetTitle(title);
  etaInc_->GetYaxis()->SetTitle("Events");
  sprintf(title, "Incident #phi");
  phiInc_ = tfile->make<TH1F>("PhiInc", title, 200, -3.1415926, 3.1415926);
  phiInc_->GetXaxis()->SetTitle(title);
  phiInc_->GetYaxis()->SetTitle("Events");
  std::string dets[ndets_] = {"EB", "EE"};
  for (int i = 0; i < ndets_; i++) {
    sprintf(name, "Hit%d", i);
    sprintf(title, "Number of hits in %s", dets[i].c_str());
    hit_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 1000, 0., 20000.);
    hit_[i]->GetXaxis()->SetTitle(title);
    hit_[i]->GetYaxis()->SetTitle("Events");
    hit_[i]->Sumw2();
    sprintf(name, "Time%d", i);
    sprintf(title, "Time of the hit (ns) in %s", dets[i].c_str());
    time_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 1000, 0., 1000.);
    time_[i]->GetXaxis()->SetTitle(title);
    time_[i]->GetYaxis()->SetTitle("Hits");
    time_[i]->Sumw2();
    sprintf(name, "TimeAll%d", i);
    sprintf(title, "Hit time (ns) in %s (for first hit in the cell)", dets[i].c_str());
    timeAll_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 1000, 0., 1000.);
    timeAll_[i]->GetXaxis()->SetTitle(title);
    timeAll_[i]->GetYaxis()->SetTitle("Hits");
    timeAll_[i]->Sumw2();
    sprintf(name, "Edep%d", i);
    sprintf(title, "Energy deposit (GeV) in %s", dets[i].c_str());
    edep_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 5000, 0., maxEnergy_);
    edep_[i]->GetXaxis()->SetTitle(title);
    edep_[i]->GetYaxis()->SetTitle("Hits");
    edep_[i]->Sumw2();
    sprintf(name, "EdepAll%d", i);
    sprintf(title, "Total Energy deposit in the cell (GeV) in %s", dets[i].c_str());
    edepAll_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 5000, 0., maxEnergy_);
    edepAll_[i]->GetXaxis()->SetTitle(title);
    edepAll_[i]->GetYaxis()->SetTitle("Hits");
    edepAll_[i]->Sumw2();
    sprintf(name, "EdepEM%d", i);
    sprintf(title, "Energy deposit (GeV) by EM particles in %s", dets[i].c_str());
    edepEM_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 5000, 0., maxEnergy_);
    edepEM_[i]->GetXaxis()->SetTitle(title);
    edepEM_[i]->GetYaxis()->SetTitle("Hits");
    edepEM_[i]->Sumw2();
    sprintf(name, "EdepHad%d", i);
    sprintf(title, "Energy deposit (GeV) by hadrons in %s", dets[i].c_str());
    edepHad_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 5000, 0., maxEnergy_);
    edepHad_[i]->GetXaxis()->SetTitle(title);
    edepHad_[i]->GetYaxis()->SetTitle("Hits");
    edepHad_[i]->Sumw2();
    sprintf(name, "Etot%d", i);
    sprintf(title, "Total energy deposit (GeV) in %s", dets[i].c_str());
    etot_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 5000, 0., maxEnergy_);
    etot_[i]->GetXaxis()->SetTitle(title);
    etot_[i]->GetYaxis()->SetTitle("Events");
    etot_[i]->Sumw2();
    sprintf(name, "EtotG%d", i);
    sprintf(title, "Total energy deposit (GeV) in %s (t < 100 ns)", dets[i].c_str());
    etotg_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 5000, 0., maxEnergy_);
    etotg_[i]->GetXaxis()->SetTitle(title);
    etotg_[i]->GetYaxis()->SetTitle("Events");
    etotg_[i]->Sumw2();
    sprintf(name, "r1by9%d", i);
    sprintf(title, "E1/E9 in %s", dets[i].c_str());
    r1by9_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 100, 0.0, 1.0);
    r1by9_[i]->GetXaxis()->SetTitle(title);
    r1by9_[i]->GetYaxis()->SetTitle("Events");
    r1by9_[i]->Sumw2();
    sprintf(name, "r1by25%d", i);
    sprintf(title, "E1/E25 in %s", dets[i].c_str());
    r1by25_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 100, 0.0, 1.0);
    r1by25_[i]->GetXaxis()->SetTitle(title);
    r1by25_[i]->GetYaxis()->SetTitle("Events");
    r1by25_[i]->Sumw2();
    sprintf(name, "r9by25%d", i);
    sprintf(title, "E9/E25 in %s", dets[i].c_str());
    r9by25_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 100, 0.0, 1.0);
    r9by25_[i]->GetXaxis()->SetTitle(title);
    r9by25_[i]->GetYaxis()->SetTitle("Events");
    r9by25_[i]->Sumw2();
    double ymax = (i == 0) ? 0.0005 : 0.005;
    sprintf(name, "sEtaEta%d", i);
    sprintf(title, "Cov(#eta,#eta) in %s", dets[i].c_str());
    sEtaEta_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 1000, 0.0, ymax);
    sEtaEta_[i]->GetXaxis()->SetTitle(title);
    sEtaEta_[i]->GetYaxis()->SetTitle("Events");
    sEtaEta_[i]->Sumw2();
    sprintf(name, "sEtaPhi%d", i);
    sprintf(title, "Cov(#eta,#phi) in %s", dets[i].c_str());
    sEtaPhi_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 1000, 0.0, ymax);
    sEtaPhi_[i]->GetXaxis()->SetTitle(title);
    sEtaPhi_[i]->GetYaxis()->SetTitle("Events");
    sEtaPhi_[i]->Sumw2();
    ymax = (i == 0) ? 0.001 : 0.01;
    sprintf(name, "sPhiPhi%d", i);
    sprintf(title, "Cov(#phi,#phi) in %s", dets[i].c_str());
    sPhiPhi_[i] = tfile->make<TH1F>(name, dets[i].c_str(), 1000, 0.0, ymax);
    sPhiPhi_[i]->GetXaxis()->SetTitle(title);
    sPhiPhi_[i]->GetYaxis()->SetTitle("Events");
    sPhiPhi_[i]->Sumw2();
    if (i == 0) {
      sprintf(title, "%s+", dets[i].c_str());
      poszp_[i] = tfile->make<TH2F>("poszp0", title, 100, 0, 100, 360, 0, 360);
      poszp_[i]->GetXaxis()->SetTitle("i#eta");
      poszp_[i]->GetYaxis()->SetTitle("i#phi");
      sprintf(title, "%s-", dets[i].c_str());
      poszn_[i] = tfile->make<TH2F>("poszn0", title, 100, 0, 100, 360, 0, 360);
      poszn_[i]->GetXaxis()->SetTitle("i#eta");
      poszn_[i]->GetYaxis()->SetTitle("i#phi");
    } else {
      sprintf(title, "%s+", dets[i].c_str());
      poszp_[i] = tfile->make<TH2F>("poszp1", title, 100, -200, 200, 100, -200, 200);
      poszp_[i]->GetXaxis()->SetTitle("x (cm)");
      poszp_[i]->GetYaxis()->SetTitle("y (cm)");
      sprintf(title, "%s-", dets[i].c_str());
      poszn_[i] = tfile->make<TH2F>("poszn1", title, 100, -200, 200, 100, -200, 200);
      poszn_[i]->GetXaxis()->SetTitle("x (cm)");
      poszn_[i]->GetYaxis()->SetTitle("y (cm)");
    }
    poszp_[i]->GetYaxis()->SetTitleOffset(1.2);
    poszp_[i]->Sumw2();
    poszn_[i]->GetYaxis()->SetTitleOffset(1.2);
    poszn_[i]->Sumw2();
  }
}

beginRun()

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

inlineoverrideprotected

Definition at line 55 of file EcalSimHitStudy.cc.

{}

endRun()

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

inlineoverrideprotected

Definition at line 56 of file EcalSimHitStudy.cc.

{}

fillDescriptions()

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

static

Definition at line 102 of file EcalSimHitStudy.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                 {
  edm::ParameterSetDescription desc;
  std::vector<std::string> calonames;
  desc.addUntracked<std::string>("ModuleLabel", "g4SimHits");
  desc.addUntracked<std::vector<std::string>>("CaloCollection", calonames);
  desc.addUntracked<std::string>("SourceLabel", "VtxSmeared");
  desc.addUntracked<double>("MaxEnergy", 200.0);
  desc.addUntracked<double>("TimeCut", 100.0);
  desc.addUntracked<int>("SelectX", -1);
  descriptions.add("EcalSimHitStudy", desc);
}

Member Data Documentation

edep_

TH1F * EcalSimHitStudy::edep_[ndets_]

private

Definition at line 76 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

edepAll_

TH1F * EcalSimHitStudy::edepAll_[ndets_]

private

Definition at line 77 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

edepEM_

TH1F* EcalSimHitStudy::edepEM_[ndets_]

private

Definition at line 76 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

edepHad_

TH1F * EcalSimHitStudy::edepHad_[ndets_]

private

Definition at line 76 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

eneInc_

TH1F * EcalSimHitStudy::eneInc_

private

Definition at line 74 of file EcalSimHitStudy.cc.

Referenced by analyze(), and beginJob().

etaInc_

TH1F * EcalSimHitStudy::etaInc_

private

Definition at line 74 of file EcalSimHitStudy.cc.

Referenced by analyze(), and beginJob().

etot_

TH1F* EcalSimHitStudy::etot_[ndets_]

private

Definition at line 77 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

etotg_

TH1F * EcalSimHitStudy::etotg_[ndets_]

private

Definition at line 77 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

g4Label_

const std::string EcalSimHitStudy::g4Label_

private

Definition at line 67 of file EcalSimHitStudy.cc.

Referenced by EcalSimHitStudy().

geometry_

const CaloGeometry* EcalSimHitStudy::geometry_

private

Definition at line 73 of file EcalSimHitStudy.cc.

Referenced by analyze(), and analyzeHits().

hit_

TH1F* EcalSimHitStudy::hit_[ndets_]

private

Definition at line 75 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

hitLab_

const std::vector<std::string> EcalSimHitStudy::hitLab_

private

Definition at line 68 of file EcalSimHitStudy.cc.

Referenced by EcalSimHitStudy().

maxEnergy_

const double EcalSimHitStudy::maxEnergy_

private

Definition at line 69 of file EcalSimHitStudy.cc.

Referenced by beginJob().

ndets_

const int EcalSimHitStudy::ndets_ = 2

staticprivate

Definition at line 65 of file EcalSimHitStudy.cc.

Referenced by beginJob().

phiInc_

TH1F * EcalSimHitStudy::phiInc_

private

Definition at line 74 of file EcalSimHitStudy.cc.

Referenced by analyze(), and beginJob().

poszn_

TH2F * EcalSimHitStudy::poszn_[ndets_]

private

Definition at line 80 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

poszp_

TH2F* EcalSimHitStudy::poszp_[ndets_]

private

Definition at line 80 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

ptInc_

TH1F* EcalSimHitStudy::ptInc_

private

Definition at line 74 of file EcalSimHitStudy.cc.

Referenced by analyze(), and beginJob().

r1by25_

TH1F * EcalSimHitStudy::r1by25_[ndets_]

private

Definition at line 78 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

r1by9_

TH1F* EcalSimHitStudy::r1by9_[ndets_]

private

Definition at line 78 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

r9by25_

TH1F * EcalSimHitStudy::r9by25_[ndets_]

private

Definition at line 78 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

selX_

const int EcalSimHitStudy::selX_

private

Definition at line 70 of file EcalSimHitStudy.cc.

Referenced by analyzeHits().

sEtaEta_

TH1F* EcalSimHitStudy::sEtaEta_[ndets_]

private

Definition at line 79 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

sEtaPhi_

TH1F * EcalSimHitStudy::sEtaPhi_[ndets_]

private

Definition at line 79 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

sPhiPhi_

TH1F * EcalSimHitStudy::sPhiPhi_[ndets_]

private

Definition at line 79 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

time_

TH1F * EcalSimHitStudy::time_[ndets_]

private

Definition at line 75 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

timeAll_

TH1F * EcalSimHitStudy::timeAll_[ndets_]

private

Definition at line 75 of file EcalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

tmax_

const double EcalSimHitStudy::tmax_

private

Definition at line 69 of file EcalSimHitStudy.cc.

Referenced by analyzeHits().

tok_evt_

const edm::EDGetTokenT<edm::HepMCProduct> EcalSimHitStudy::tok_evt_

private

Definition at line 71 of file EcalSimHitStudy.cc.

Referenced by analyze().

tokGeom_

const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> EcalSimHitStudy::tokGeom_

private

Definition at line 66 of file EcalSimHitStudy.cc.

Referenced by analyze().

toks_calo_

const std::vector<edm::EDGetTokenT<edm::PCaloHitContainer> > EcalSimHitStudy::toks_calo_

private

Definition at line 72 of file EcalSimHitStudy.cc.

Referenced by analyze().

w0_

const double EcalSimHitStudy::w0_

private

Definition at line 69 of file EcalSimHitStudy.cc.

Referenced by analyzeHits().