XtalDedxAnalysis Class Reference

Inheritance diagram for XtalDedxAnalysis:

Public Member Functions
	XtalDedxAnalysis (const edm::ParameterSet &)
	~XtalDedxAnalysis () override
Public Member Functions inherited from edm::one::EDAnalyzer< 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 &&)=default
	EDConsumerBase (EDConsumerBase const &)=delete
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 &	operator= (EDConsumerBase &&)=default
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
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 (const edm::Event &, const edm::EventSetup &) override
void	analyzeHits (std::vector< PCaloHit > &, edm::Handle< edm::SimTrackContainer > &, edm::Handle< edm::SimVertexContainer > &)
void	beginJob () override
void	endJob () override
Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Private Attributes
edm::InputTag	caloHitSource_
TH1F *	meE1T0_ [4]
TH1F *	meE1T1_ [4]
TH1F *	meE9T0_ [4]
TH1F *	meE9T1_ [4]
TH1F *	meNHit_ [4]
TH1I *	mType_
std::string	simTkLabel_
edm::EDGetTokenT< edm::PCaloHitContainer >	tok_calo_
edm::EDGetTokenT< edm::SimTrackContainer >	tok_tk_
edm::EDGetTokenT< edm::SimVertexContainer >	tok_vtx_

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 36 of file XtalDedxAnalysis.cc.

Constructor & Destructor Documentation

XtalDedxAnalysis()

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

explicit

Definition at line 64 of file XtalDedxAnalysis.cc.

                                                            {
  usesResource(TFileService::kSharedResource);
  caloHitSource_ = ps.getParameter<edm::InputTag>("caloHitSource");
  simTkLabel_ = ps.getParameter<std::string>("moduleLabelTk");
  double energyMax = ps.getParameter<double>("energyMax");
  edm::LogVerbatim("CherenkovAnalysis") << "XtalDedxAnalysis::Source " << caloHitSource_ << " Track Label "
                                        << simTkLabel_ << " Energy Max " << energyMax;
  // register for data access
  tok_calo_ = consumes<edm::PCaloHitContainer>(caloHitSource_);
  tok_tk_ = consumes<edm::SimTrackContainer>(edm::InputTag(simTkLabel_));
  tok_vtx_ = consumes<edm::SimVertexContainer>(edm::InputTag(simTkLabel_));
 
  // Book histograms
  edm::Service<TFileService> tfile;
 
  if (!tfile.isAvailable())
    throw cms::Exception("BadConfig") << "TFileService unavailable: "
                                      << "please add it to config file";
  // Histograms for Hits
  std::string types[4] = {"total", "by dE/dx", "by delta-ray", "by bremms"};
  char name[20], title[80];
  for (int i = 0; i < 4; i++) {
    sprintf(name, "Hits%d", i);
    sprintf(title, "Number of hits (%s)", types[i].c_str());
    meNHit_[i] = tfile->make<TH1F>(name, title, 5000, 0., 5000.);
    meNHit_[i]->GetXaxis()->SetTitle(title);
    meNHit_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "E1T0%d", i);
    sprintf(title, "E1 (Loss %s) in GeV", types[i].c_str());
    meE1T0_[i] = tfile->make<TH1F>(name, title, 5000, 0, energyMax);
    meE1T0_[i]->GetXaxis()->SetTitle(title);
    meE1T0_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "E9T0%d", i);
    sprintf(title, "E9 (Loss %s) in GeV", types[i].c_str());
    meE9T0_[i] = tfile->make<TH1F>(name, title, 5000, 0, energyMax);
    meE9T0_[i]->GetXaxis()->SetTitle(title);
    meE9T0_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "E1T1%d", i);
    sprintf(title, "E1 (Loss %s with t < 400 ns) in GeV", types[i].c_str());
    meE1T1_[i] = tfile->make<TH1F>(name, title, 5000, 0, energyMax);
    meE1T1_[i]->GetXaxis()->SetTitle(title);
    meE1T1_[i]->GetYaxis()->SetTitle("Events");
    sprintf(name, "E9T1%d", i);
    sprintf(title, "E9 (Loss %s with t < 400 ns) in GeV", types[i].c_str());
    meE9T1_[i] = tfile->make<TH1F>(name, title, 5000, 0, energyMax);
    meE9T1_[i]->GetXaxis()->SetTitle(title);
    meE9T1_[i]->GetYaxis()->SetTitle("Events");
  }
  sprintf(name, "PDGType");
  sprintf(title, "PDG ID of first level secondary");
  mType_ = tfile->make<TH1I>(name, title, 5000, -2500, 2500);
  mType_->GetXaxis()->SetTitle(title);
  mType_->GetYaxis()->SetTitle("Tracks");
}

References caloHitSource_, edm::ParameterSet::getParameter(), mps_fire::i, TFileService::kSharedResource, meE1T0_, meE1T1_, meE9T0_, meE9T1_, meNHit_, mType_, Skims_PA_cff::name, simTkLabel_, AlCaHLTBitMon_QueryRunRegistry::string, compare::tfile, runGCPTkAlMap::title, tok_calo_, tok_tk_, and tok_vtx_.

~XtalDedxAnalysis()

XtalDedxAnalysis::~XtalDedxAnalysis ( )

inlineoverride

Definition at line 39 of file XtalDedxAnalysis.cc.

{}

Member Function Documentation

analyze()

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

overrideprotectedvirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 127 of file XtalDedxAnalysis.cc.

                                                                       {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("CherenkovAnalysis") << "XtalDedxAnalysis::Run = " << e.id().run() << " Event = " << e.id().event();
#endif
  std::vector<PCaloHit> caloHits;
  edm::Handle<edm::PCaloHitContainer> pCaloHits;
  e.getByToken(tok_calo_, pCaloHits);
 
  std::vector<SimTrack> theSimTracks;
  edm::Handle<edm::SimTrackContainer> simTk;
  e.getByToken(tok_tk_, simTk);
 
  std::vector<SimVertex> theSimVertex;
  edm::Handle<edm::SimVertexContainer> simVtx;
  e.getByToken(tok_vtx_, simVtx);
 
  if (pCaloHits.isValid()) {
    caloHits.insert(caloHits.end(), pCaloHits->begin(), pCaloHits->end());
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("CherenkovAnalysis") << "XtalDedxAnalysis: Hit buffer " << caloHits.size();
#endif
    analyzeHits(caloHits, simTk, simVtx);
  }
}

References analyzeHits(), MillePedeFileConverter_cfg::e, edm::HandleBase::isValid(), tok_calo_, tok_tk_, and tok_vtx_.

analyzeHits()

void XtalDedxAnalysis::analyzeHits ( std::vector< PCaloHit > &  hits, edm::Handle< edm::SimTrackContainer > &  SimTk, edm::Handle< edm::SimVertexContainer > &  SimVtx  )

protected

Definition at line 152 of file XtalDedxAnalysis.cc.

                                                                             {
  edm::SimTrackContainer::const_iterator simTrkItr;
  int nHit = hits.size();
  double e10[4], e90[4], e11[4], e91[4], hit[4];
  for (int i = 0; i < 4; i++)
    e10[i] = e90[i] = e11[i] = e91[i] = hit[i] = 0;
  for (int i = 0; i < nHit; i++) {
    double energy = hits[i].energy();
    double time = hits[i].time();
    unsigned int id_ = hits[i].id();
    int trackID = hits[i].geantTrackId();
    int type = 1;
    for (simTrkItr = SimTk->begin(); simTrkItr != SimTk->end(); simTrkItr++) {
      if (trackID == (int)(simTrkItr->trackId())) {
        int thePID = simTrkItr->type();
        if (thePID == 11)
          type = 2;
        else if (thePID != -13 && thePID != 13)
          type = 3;
        break;
      }
    }
    hit[0]++;
    hit[type]++;
    e90[0] += energy;
    e90[type] += energy;
    if (time < 400) {
      e91[0] += energy;
      e91[type] += energy;
    }
    if (id_ == 22) {
      e10[0] += energy;
      e10[type] += energy;
      if (time < 400) {
        e11[0] += energy;
        e11[type] += energy;
      }
    }
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("CherenkovAnalysis") << "XtalDedxAnalysis:Hit[" << i << "] ID " << id_ << " E " << energy
                                          << " time " << time << " track " << trackID << " type " << type;
#endif
  }
  for (int i = 0; i < 4; i++) {
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("CherenkovAnalysis") << "XtalDedxAnalysis:Type(" << i << ") Hit " << hit[i] << " E10 " << e10[i]
                                          << " E11 " << e11[i] << " E90 " << e90[i] << " E91 " << e91[i];
#endif
    meNHit_[i]->Fill(hit[i]);
    meE1T0_[i]->Fill(e10[i]);
    meE9T0_[i]->Fill(e90[i]);
    meE1T1_[i]->Fill(e11[i]);
    meE9T1_[i]->Fill(e91[i]);
  }
 
  // Type of the secondary (coming directly from a generator level track)
  int nvtx = 0, ntrk = 0, k1 = 0;
  edm::SimVertexContainer::const_iterator simVtxItr;
  for (simTrkItr = SimTk->begin(); simTrkItr != SimTk->end(); simTrkItr++)
    ntrk++;
  for (simVtxItr = SimVtx->begin(); simVtxItr != SimVtx->end(); simVtxItr++)
    nvtx++;
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("CherenkovAnalysis") << "XtalDedxAnalysis: " << ntrk << " tracks and " << nvtx << " vertices";
#endif
  for (simTrkItr = SimTk->begin(); simTrkItr != SimTk->end(); simTrkItr++, ++k1) {
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("CherenkovAnalysis") << "Track " << k1 << " PDGId " << simTrkItr->type() << " Vertex ID "
                                          << simTrkItr->vertIndex() << " Generator " << simTrkItr->noGenpart();
#endif
    if (simTrkItr->noGenpart()) {              // This is a secondary
      int vertIndex = simTrkItr->vertIndex();  // Vertex index of origin
      if (vertIndex >= 0 && vertIndex < nvtx) {
        simVtxItr = SimVtx->begin();
        for (int iv = 0; iv < vertIndex; iv++)
          simVtxItr++;
        int parent = simVtxItr->parentIndex(), k2 = 0;
        for (edm::SimTrackContainer::const_iterator trkItr = SimTk->begin(); trkItr != SimTk->end(); trkItr++, ++k2) {
#ifdef EDM_ML_DEBUG
          edm::LogVerbatim("CherenkovAnalysis") << "XtalDedxAnalysis::Track " << k2 << " ID " << trkItr->trackId()
                                                << " (" << parent << ")  Generator " << trkItr->noGenpart();
#endif
          if ((int)trkItr->trackId() == parent) {  // Parent track
            if (!trkItr->noGenpart()) {            // Generator level
#ifdef EDM_ML_DEBUG
              edm::LogVerbatim("CherenkovAnalysis") << "XtalDedxAnalysis::Track found with ID " << simTrkItr->type();
#endif
              mType_->Fill(simTrkItr->type());
            }
            break;
          }
        }
      }
    }
  }
}

References HCALHighEnergyHPDFilter_cfi::energy, hfClusterShapes_cfi::hits, mps_fire::i, gpuVertexFinder::iv, meE1T0_, meE1T1_, meE9T0_, meE9T1_, meNHit_, mType_, class-composition::parent, and protons_cff::time.

Referenced by analyze().

beginJob()

void XtalDedxAnalysis::beginJob ( void  )

inlineoverrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 44 of file XtalDedxAnalysis.cc.

{}

endJob()

void XtalDedxAnalysis::endJob ( void  )

inlineoverrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 46 of file XtalDedxAnalysis.cc.

{}

fillDescriptions()

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

static

Definition at line 119 of file XtalDedxAnalysis.cc.

                                                                                  {
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>("caloHitSource", edm::InputTag("g4SimHits", "EcalHitsEB"));
  desc.add<std::string>("moduleLabelTk", "g4SimHits");
  desc.add<double>("energyMax", 2.0);
  descriptions.add("xtalDedxAnalysis", desc);
}

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

Member Data Documentation

caloHitSource_

edm::InputTag XtalDedxAnalysis::caloHitSource_

private

Definition at line 53 of file XtalDedxAnalysis.cc.

Referenced by XtalDedxAnalysis().

meE1T0_

TH1F * XtalDedxAnalysis::meE1T0_[4]

private

Definition at line 60 of file XtalDedxAnalysis.cc.

Referenced by analyzeHits(), and XtalDedxAnalysis().

meE1T1_

TH1F * XtalDedxAnalysis::meE1T1_[4]

private

Definition at line 60 of file XtalDedxAnalysis.cc.

Referenced by analyzeHits(), and XtalDedxAnalysis().

meE9T0_

TH1F * XtalDedxAnalysis::meE9T0_[4]

private

Definition at line 60 of file XtalDedxAnalysis.cc.

Referenced by analyzeHits(), and XtalDedxAnalysis().

meE9T1_

TH1F * XtalDedxAnalysis::meE9T1_[4]

private

Definition at line 60 of file XtalDedxAnalysis.cc.

Referenced by analyzeHits(), and XtalDedxAnalysis().

meNHit_

TH1F* XtalDedxAnalysis::meNHit_[4]

private

Definition at line 60 of file XtalDedxAnalysis.cc.

Referenced by analyzeHits(), and XtalDedxAnalysis().

mType_

TH1I* XtalDedxAnalysis::mType_

private

Definition at line 61 of file XtalDedxAnalysis.cc.

Referenced by analyzeHits(), and XtalDedxAnalysis().

simTkLabel_

std::string XtalDedxAnalysis::simTkLabel_

private

Definition at line 54 of file XtalDedxAnalysis.cc.

Referenced by XtalDedxAnalysis().

tok_calo_

edm::EDGetTokenT<edm::PCaloHitContainer> XtalDedxAnalysis::tok_calo_

private

Definition at line 56 of file XtalDedxAnalysis.cc.

Referenced by analyze(), and XtalDedxAnalysis().

tok_tk_

edm::EDGetTokenT<edm::SimTrackContainer> XtalDedxAnalysis::tok_tk_

private

Definition at line 57 of file XtalDedxAnalysis.cc.

Referenced by analyze(), and XtalDedxAnalysis().

tok_vtx_

edm::EDGetTokenT<edm::SimVertexContainer> XtalDedxAnalysis::tok_vtx_

private

Definition at line 58 of file XtalDedxAnalysis.cc.

Referenced by analyze(), and XtalDedxAnalysis().