HGCalTBAnalyzer Class Reference

Inheritance diagram for HGCalTBAnalyzer:

Public Member Functions
	HGCalTBAnalyzer (edm::ParameterSet const &)
	~HGCalTBAnalyzer () override
Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::SharedResources >
	EDAnalyzer ()=default
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () 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
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::vector< ModuleDescription const * > &modules, 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
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
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)

Private Member Functions
void	analyze (edm::Event const &, edm::EventSetup const &) override
template<class T1 >
void	analyzeDigi (int type, const T1 &detId, uint16_t adc)
void	analyzePassiveHits (edm::Handle< edm::PassiveHitContainer > const &hgcPh, int subdet)
void	analyzeRecHits (int type, edm::Handle< HGCRecHitCollection > &hits)
void	analyzeSimHits (int type, std::vector< PCaloHit > &hits, double zFront)
void	analyzeSimTracks (edm::Handle< edm::SimTrackContainer > const &SimTk, edm::Handle< edm::SimVertexContainer > const &SimVtx)
void	beginJob () override
void	beginRun (edm::Run const &, edm::EventSetup const &) override
void	endRun (edm::Run const &, edm::EventSetup const &) override

Private Attributes
std::string	detectorBeam_
std::string	detectorBH_
std::string	detectorEE_
std::string	detectorFH_
bool	doDigis_
bool	doPassive_
bool	doRecHits_
bool	doSimHits_
bool	doTree_
bool	doTreeCell_
edm::Service< TFileService >	fs_
TH1D *	hBeam_
TH1D *	hDigiADC_ [3]
TH1D *	hDigiLng_ [2]
TH2D *	hDigiOcc_ [3]
const HGCalDDDConstants *	hgcons_ [2]
std::vector< float >	hgcPassiveBHEnergy_
std::vector< int >	hgcPassiveBHID_
std::vector< std::string >	hgcPassiveBHName_
std::vector< float >	hgcPassiveCMSEEnergy_
std::vector< int >	hgcPassiveCMSEID_
std::vector< std::string >	hgcPassiveCMSEName_
std::vector< float >	hgcPassiveEEEnergy_
std::vector< int >	hgcPassiveEEID_
std::vector< std::string >	hgcPassiveEEName_
std::vector< float >	hgcPassiveFHEnergy_
std::vector< int >	hgcPassiveFHID_
std::vector< std::string >	hgcPassiveFHName_
const HGCalGeometry *	hgeom_ [2]
TH1D *	hRecHitE_ [3]
TProfile2D *	hRecHitLat_ [3]
TProfile *	hRecHitLng1_ [3]
TProfile *	hRecHitLng_ [3]
TH2D *	hRecHitOcc_ [3]
TH1D *	hSimHitE_ [4]
TH1D *	hSimHitEn_ [4]
TProfile2D *	hSimHitLat_ [3]
std::vector< TH1D * >	hSimHitLayEn1BH_
std::vector< TH1D * >	hSimHitLayEn1EE_
std::vector< TH1D * >	hSimHitLayEn1FH_
std::vector< TH1D * >	hSimHitLayEn2BH_
std::vector< TH1D * >	hSimHitLayEn2EE_
std::vector< TH1D * >	hSimHitLayEn2FH_
std::vector< TH1D * >	hSimHitLayEnBeam_
TProfile *	hSimHitLng1_ [3]
TProfile *	hSimHitLng2_ [3]
TProfile *	hSimHitLng_ [3]
TH1D *	hSimHitT_ [4]
std::vector< int >	idBeams_
bool	ifBeam_
bool	ifBH_
bool	ifEE_
bool	ifFH_
double	pBeam_
int	sampleIndex_
std::vector< float >	simHitCellEnBeam_
std::vector< float >	simHitCellEnBH_
std::vector< float >	simHitCellEnEE_
std::vector< float >	simHitCellEnFH_
std::vector< uint32_t >	simHitCellIdBeam_
std::vector< uint32_t >	simHitCellIdBH_
std::vector< uint32_t >	simHitCellIdEE_
std::vector< uint32_t >	simHitCellIdFH_
std::vector< float >	simHitLayEn1BH_
std::vector< float >	simHitLayEn1EE_
std::vector< float >	simHitLayEn1FH_
std::vector< float >	simHitLayEn2BH_
std::vector< float >	simHitLayEn2EE_
std::vector< float >	simHitLayEn2FH_
std::vector< float >	simHitLayEnBeam_
edm::EDGetToken	tok_digiBH_
edm::EDGetToken	tok_digiEE_
edm::EDGetToken	tok_digiFH_
edm::EDGetTokenT< edm::HepMCProduct >	tok_hepMC_
edm::EDGetTokenT< edm::PassiveHitContainer >	tok_hgcPHBH_
edm::EDGetTokenT< edm::PassiveHitContainer >	tok_hgcPHCMSE_
edm::EDGetTokenT< edm::PassiveHitContainer >	tok_hgcPHEE_
edm::EDGetTokenT< edm::PassiveHitContainer >	tok_hgcPHFH_
edm::EDGetToken	tok_hitrBH_
edm::EDGetToken	tok_hitrEE_
edm::EDGetToken	tok_hitrFH_
edm::EDGetTokenT< edm::PCaloHitContainer >	tok_hitsBeam_
edm::EDGetTokenT< edm::PCaloHitContainer >	tok_hitsBH_
edm::EDGetTokenT< edm::PCaloHitContainer >	tok_hitsEE_
edm::EDGetTokenT< edm::PCaloHitContainer >	tok_hitsFH_
edm::EDGetTokenT< edm::SimTrackContainer >	tok_simTk_
edm::EDGetTokenT< edm::SimVertexContainer >	tok_simVtx_
TTree *	tree_
double	xBeam_
double	yBeam_
double	zBeam_
double	zFrontBH_
double	zFrontEE_
double	zFrontFH_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
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 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)

Detailed Description

Definition at line 53 of file HGCalTBAnalyzer.cc.

Constructor & Destructor Documentation

HGCalTBAnalyzer::HGCalTBAnalyzer ( edm::ParameterSet const &  )

explicit

Passive hits

Definition at line 125 of file HGCalTBAnalyzer.cc.

References gather_cfg::cout, detectorBeam_, detectorBH_, detectorEE_, detectorFH_, doDigis_, doPassive_, doRecHits_, doSimHits_, doTree_, doTreeCell_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), triggerObjects_cff::id, idBeams_, ifBeam_, ifBH_, ifEE_, ifFH_, sampleIndex_, AlCaHLTBitMon_QueryRunRegistry::string, tmp, tok_digiBH_, tok_digiEE_, tok_digiFH_, tok_hepMC_, tok_hgcPHBH_, tok_hgcPHCMSE_, tok_hgcPHEE_, tok_hgcPHFH_, tok_hitrBH_, tok_hitrEE_, tok_hitrFH_, tok_hitsBeam_, tok_hitsBH_, tok_hitsEE_, tok_hitsFH_, tok_simTk_, tok_simVtx_, zFrontBH_, zFrontEE_, and zFrontFH_.

                                                               {

  usesResource("TFileService");

  //now do whatever initialization is needed
  detectorEE_  = iConfig.getParameter<std::string>("DetectorEE");
  detectorFH_  = iConfig.getParameter<std::string>("DetectorFH");
  detectorBH_  = iConfig.getParameter<std::string>("DetectorBH");
  detectorBeam_= iConfig.getParameter<std::string>("DetectorBeam");
  ifEE_        = iConfig.getParameter<bool>("UseEE");
  ifFH_        = iConfig.getParameter<bool>("UseFH");
  ifBH_        = iConfig.getParameter<bool>("UseBH");
  ifBeam_      = iConfig.getParameter<bool>("UseBeam");
  zFrontEE_    = iConfig.getParameter<double>("ZFrontEE");
  zFrontFH_    = iConfig.getParameter<double>("ZFrontFH");
  zFrontBH_    = iConfig.getParameter<double>("ZFrontBH");
  idBeams_     = iConfig.getParameter<std::vector<int>>("IDBeams");  
  doSimHits_   = iConfig.getParameter<bool>("DoSimHits");
  doDigis_     = iConfig.getParameter<bool>("DoDigis");
  sampleIndex_ = iConfig.getParameter<int>("SampleIndex");
  doRecHits_   = iConfig.getParameter<bool>("DoRecHits");
  doTree_      = iConfig.getUntrackedParameter<bool>("DoTree",false);
  doTreeCell_  = iConfig.getUntrackedParameter<bool>("DoTreeCell",false);
  doPassive_   = iConfig.getUntrackedParameter<bool>("DoPassive",false);

#ifdef EDM_ML_DEBUG
  std::cout << "HGCalTBAnalyzer:: SimHits = " << doSimHits_ << " Digis = "
        << doDigis_ << ":" << sampleIndex_ << " RecHits = " << doRecHits_
        << " useDets " << ifEE_ << ":" << ifFH_ << ":" << ifBH_ << ":"
        << ifBeam_ << " zFront " << zFrontEE_ << ":" << zFrontFH_ << ":" 
        << zFrontBH_ << " IdBeam " << idBeams_.size() << ":";
  for (auto id : idBeams_) std::cout << " " << id;
  std::cout << std::endl;
#endif
  if (idBeams_.empty()) idBeams_.push_back(1001);

  edm::InputTag tmp0 = iConfig.getParameter<edm::InputTag>("GeneratorSrc");
  tok_hepMC_   = consumes<edm::HepMCProduct>(tmp0);

#ifdef EDM_ML_DEBUG
  std::cout << "HGCalTBAnalyzer:: GeneratorSource = " << tmp0 << std::endl;
#endif
  std::string   tmp1 = iConfig.getParameter<std::string>("CaloHitSrcEE");
  tok_hitsEE_  = consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits",tmp1));
  tok_simTk_   = consumes<edm::SimTrackContainer>(edm::InputTag("g4SimHits"));
  tok_simVtx_  = consumes<edm::SimVertexContainer>(edm::InputTag("g4SimHits"));
  edm::InputTag tmp2 = iConfig.getParameter<edm::InputTag>("DigiSrcEE");
  tok_digiEE_  = consumes<HGCalDigiCollection>(tmp2);
  edm::InputTag tmp3 = iConfig.getParameter<edm::InputTag>("RecHitSrcEE");
  tok_hitrEE_  = consumes<HGCRecHitCollection>(tmp3);
#ifdef EDM_ML_DEBUG
  if (ifEE_) {
    std::cout << "HGCalTBAnalyzer:: Detector " << detectorEE_ << " with tags "
          << tmp1 << ", " << tmp2 << ", " << tmp3 << std::endl;
  }
#endif
  tmp1         = iConfig.getParameter<std::string>("CaloHitSrcFH");
  tok_hitsFH_  = consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits",tmp1));
  tmp2         = iConfig.getParameter<edm::InputTag>("DigiSrcFH");
  tok_digiFH_  = consumes<HGCalDigiCollection>(tmp2);
  tmp3         = iConfig.getParameter<edm::InputTag>("RecHitSrcFH");
  tok_hitrFH_  = consumes<HGCRecHitCollection>(tmp3);
#ifdef EDM_ML_DEBUG
  if (ifFH_) {
    std::cout << "HGCalTBAnalyzer:: Detector " << detectorFH_ << " with tags "
          << tmp1 << ", " << tmp2 << ", " << tmp3 << std::endl;
  }
#endif
  tmp1         = iConfig.getParameter<std::string>("CaloHitSrcBH");
  tok_hitsBH_  = consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits",tmp1));
  tmp2         = iConfig.getParameter<edm::InputTag>("DigiSrcBH");
  tok_digiBH_  = consumes<HGCalDigiCollection>(tmp2);
  tmp3         = iConfig.getParameter<edm::InputTag>("RecHitSrcBH");
  tok_hitrBH_  = consumes<HGCRecHitCollection>(tmp3);

  edm::InputTag tmp = iConfig.getParameter<edm::InputTag>("HGCPassiveEE");
  tok_hgcPHEE_   = consumes<edm::PassiveHitContainer>(tmp);

  tmp = iConfig.getParameter<edm::InputTag>("HGCPassiveFH");
  tok_hgcPHFH_   = consumes<edm::PassiveHitContainer>(tmp);

  tmp = iConfig.getParameter<edm::InputTag>("HGCPassiveBH");
  tok_hgcPHBH_   = consumes<edm::PassiveHitContainer>(tmp);

  tmp = iConfig.getParameter<edm::InputTag>("HGCPassiveCMSE");
  tok_hgcPHCMSE_   = consumes<edm::PassiveHitContainer>(tmp);

#ifdef EDM_ML_DEBUG
  if (ifBH_) {
    std::cout << "HGCalTBAnalyzer:: Detector " << detectorBH_ << " with tags "
          << tmp1 << ", " << tmp2 << ", " << tmp3 << std::endl;
  }
#endif
  tmp1         = iConfig.getParameter<std::string>("CaloHitSrcBeam");
  tok_hitsBeam_= consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits",tmp1));
#ifdef EDM_ML_DEBUG
  if (ifBeam_) {
    std::cout << "HGCalTBAnalyzer:: Detector " << detectorBeam_ 
          << " with tags " << tmp1 << std::endl;
  }
#endif
}

HGCalTBAnalyzer::~HGCalTBAnalyzer ( )

override

Definition at line 229 of file HGCalTBAnalyzer.cc.

{}

Member Function Documentation

void HGCalTBAnalyzer::analyze ( edm::Event const &  , edm::EventSetup const &    )

overrideprivatevirtual

EE

FH

BH

CMSE

Implements edm::one::EDAnalyzerBase.

Definition at line 484 of file HGCalTBAnalyzer.cc.

References ecalMGPA::adc(), analyzeDigi(), analyzePassiveHits(), analyzeRecHits(), analyzeSimHits(), analyzeSimTracks(), edm::HandleBase::clear(), gather_cfg::cout, HGCSample::data(), detectorBeam_, detectorBH_, detectorEE_, detectorFH_, doDigis_, doPassive_, doRecHits_, doSimHits_, doTree_, edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), hBeam_, hgcons_, hgcPassiveBHEnergy_, hgcPassiveBHID_, hgcPassiveBHName_, hgcPassiveCMSEEnergy_, hgcPassiveCMSEID_, hgcPassiveCMSEName_, hgcPassiveEEEnergy_, hgcPassiveEEID_, hgcPassiveEEName_, hgcPassiveFHEnergy_, hgcPassiveFHID_, hgcPassiveFHName_, idBeams_, ifBeam_, ifBH_, ifEE_, ifFH_, edm::HandleBase::isValid(), gen::k, HGCalDDDConstants::layers(), AHCalDetId::MaxDepth, AlCaHLTBitMon_ParallelJobs::p, sampleIndex_, simHitCellEnBeam_, simHitCellEnBH_, simHitCellEnEE_, simHitCellEnFH_, simHitCellIdBeam_, simHitCellIdBH_, simHitCellIdEE_, simHitCellIdFH_, simHitLayEn1BH_, simHitLayEn1EE_, simHitLayEn1FH_, simHitLayEn2BH_, simHitLayEn2EE_, simHitLayEn2FH_, simHitLayEnBeam_, edm::SortedCollection< T, SORT >::size(), tok_digiEE_, tok_digiFH_, tok_hepMC_, tok_hgcPHBH_, tok_hgcPHCMSE_, tok_hgcPHEE_, tok_hgcPHFH_, tok_hitrEE_, tok_hitrFH_, tok_hitsBeam_, tok_hitsBH_, tok_hitsEE_, tok_hitsFH_, tok_simTk_, tok_simVtx_, tree_, zFrontBH_, zFrontEE_, and zFrontFH_.

Referenced by endRun().

                                               {

  //Generator input
  edm::Handle<edm::HepMCProduct> evtMC;
  iEvent.getByToken(tok_hepMC_,evtMC);
  if (!evtMC.isValid()) {
    edm::LogWarning("HGCal") << "no HepMCProduct found";
  } else { 
    const HepMC::GenEvent * myGenEvent = evtMC->GetEvent();
    unsigned int k(0);
    for (HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin();
         p != myGenEvent->particles_end(); ++p, ++k) {
      if (k == 0) hBeam_->Fill((*p)->momentum().rho());
#ifdef EDM_ML_DEBUG
      std::cout << "Particle[" << k << "] with p " << (*p)->momentum().rho() 
        << " theta " << (*p)->momentum().theta() << " phi "
        << (*p)->momentum().phi() << std::endl;
#endif
    }
  }

  //Now the Simhits
  if (doSimHits_) {
    edm::Handle<edm::SimTrackContainer>  SimTk;
    iEvent.getByToken(tok_simTk_, SimTk);
    edm::Handle<edm::SimVertexContainer> SimVtx;
    iEvent.getByToken(tok_simVtx_, SimVtx);
    analyzeSimTracks(SimTk, SimVtx);
  
    simHitLayEn1EE_.clear(); simHitLayEn2EE_.clear();
    simHitLayEn1FH_.clear(); simHitLayEn2FH_.clear();
    simHitLayEn1BH_.clear(); simHitLayEn2BH_.clear();
    simHitLayEnBeam_.clear();
    simHitCellIdEE_.clear(); simHitCellEnEE_.clear();
    simHitCellIdFH_.clear(); simHitCellEnFH_.clear();
    simHitCellIdBH_.clear(); simHitCellEnBH_.clear();
    simHitCellIdBeam_.clear(); simHitCellEnBeam_.clear();
    edm::Handle<edm::PCaloHitContainer> theCaloHitContainers;
    std::vector<PCaloHit>               caloHits;
    if (ifEE_) {
      simHitLayEn1EE_ = std::vector<float>(hgcons_[0]->layers(false),0);
      simHitLayEn2EE_ = std::vector<float>(hgcons_[0]->layers(true),0);
      iEvent.getByToken(tok_hitsEE_, theCaloHitContainers);
      if (theCaloHitContainers.isValid()) {
#ifdef EDM_ML_DEBUG
    std::cout << "PcalohitContainer for " << detectorEE_ << " has "
          << theCaloHitContainers->size() << " hits" << std::endl;
#endif
    caloHits.clear();
    caloHits.insert(caloHits.end(), theCaloHitContainers->begin(), 
            theCaloHitContainers->end());
    analyzeSimHits(0, caloHits, zFrontEE_);
      } else {
#ifdef EDM_ML_DEBUG
    std::cout << "PCaloHitContainer does not exist for " << detectorEE_ 
          << " !!!" << std::endl;
#endif
      }
    }
    if (ifFH_) {
      simHitLayEn1FH_ = std::vector<float>(hgcons_[1]->layers(false),0);
      simHitLayEn2FH_ = std::vector<float>(hgcons_[1]->layers(true),0);
      iEvent.getByToken(tok_hitsFH_, theCaloHitContainers);
      if (theCaloHitContainers.isValid()) {
#ifdef EDM_ML_DEBUG
    std::cout << "PcalohitContainer for " << detectorFH_ << " has "
          << theCaloHitContainers->size() << " hits" << std::endl;
#endif
    caloHits.clear();
    caloHits.insert(caloHits.end(), theCaloHitContainers->begin(), 
            theCaloHitContainers->end());
    analyzeSimHits(1, caloHits, zFrontFH_);
      } else {
#ifdef EDM_ML_DEBUG
    std::cout << "PCaloHitContainer does not exist for " << detectorFH_ 
          << " !!!" << std::endl;
#endif
      }
    }
    if (ifBH_) {
      simHitLayEn1BH_ = std::vector<float>(AHCalDetId::MaxDepth,0);
      simHitLayEn2BH_ = std::vector<float>(AHCalDetId::MaxDepth,0);
      iEvent.getByToken(tok_hitsBH_, theCaloHitContainers);
      if (theCaloHitContainers.isValid()) {
#ifdef EDM_ML_DEBUG
    std::cout << "PcalohitContainer for " << detectorBH_ << " has "
          << theCaloHitContainers->size() << " hits" << std::endl;
#endif
    caloHits.clear();
    caloHits.insert(caloHits.end(), theCaloHitContainers->begin(), 
            theCaloHitContainers->end());
    analyzeSimHits(2, caloHits, zFrontBH_);
      } else {
#ifdef EDM_ML_DEBUG
    std::cout << "PCaloHitContainer does not exist for " << detectorBH_ 
          << " !!!" << std::endl;
#endif
      }
    }
    if (ifBeam_) {
      simHitLayEnBeam_ = std::vector<float>(idBeams_.size(),0);
      iEvent.getByToken(tok_hitsBeam_, theCaloHitContainers);
      if (theCaloHitContainers.isValid()) {
#ifdef EDM_ML_DEBUG
    std::cout << "PcalohitContainer for " << detectorBeam_ << " has "
          << theCaloHitContainers->size() << " hits" << std::endl;
#endif
    caloHits.clear();
    caloHits.insert(caloHits.end(), theCaloHitContainers->begin(), 
            theCaloHitContainers->end());
    analyzeSimHits(3, caloHits, 0.0);
      } else {
#ifdef EDM_ML_DEBUG
    std::cout << "PCaloHitContainer does not exist for " << detectorBeam_ 
          << " !!!" << std::endl;
#endif
      }
    }
    //if (doTree_) tree_->Fill();
  }//if (doSimHits_)

  if (doPassive_) {
    hgcPassiveEEEnergy_.clear();     hgcPassiveFHEnergy_.clear();
    hgcPassiveBHEnergy_.clear();     hgcPassiveCMSEEnergy_.clear();
    hgcPassiveEEName_.clear();       hgcPassiveFHName_.clear();
    hgcPassiveBHName_.clear();       hgcPassiveCMSEName_.clear();
    hgcPassiveEEID_.clear();         hgcPassiveFHID_.clear();
    hgcPassiveBHID_.clear();         hgcPassiveCMSEID_.clear();
    
    edm::Handle<edm::PassiveHitContainer>  hgcPHEE;
    iEvent.getByToken(tok_hgcPHEE_,hgcPHEE);
    analyzePassiveHits(hgcPHEE, 1);

    edm::Handle<edm::PassiveHitContainer>  hgcPHFH;
    iEvent.getByToken(tok_hgcPHFH_,hgcPHFH);
    analyzePassiveHits(hgcPHFH, 2);

    edm::Handle<edm::PassiveHitContainer>  hgcPHBH;
    iEvent.getByToken(tok_hgcPHBH_,hgcPHBH);
    analyzePassiveHits(hgcPHBH, 3);

    edm::Handle<edm::PassiveHitContainer>  hgcPHCMSE;
    iEvent.getByToken(tok_hgcPHCMSE_,hgcPHCMSE);
    analyzePassiveHits(hgcPHCMSE, 4);
  }

  if ((doSimHits_ || doPassive_) && (doTree_)) tree_->Fill();

  //Now the Digis
  if (doDigis_) {
    if (ifEE_) {
      edm::Handle<HGCalDigiCollection> theDigiContainers;
      iEvent.getByToken(tok_digiEE_, theDigiContainers);
      if (theDigiContainers.isValid()) {
#ifdef EDM_ML_DEBUG
    std::cout << "HGCDigiCintainer for " << detectorEE_ << " with " 
          << theDigiContainers->size() << " element(s)" << std::endl;
#endif
    for (auto it : *theDigiContainers) {
      HGCalDetId detId     = (it.id());
      const HGCSample&  hgcSample = it.sample(sampleIndex_);
      uint16_t   adc       = hgcSample.data();
      analyzeDigi(0, detId, adc);
    }
      }
    }
    if (ifFH_) {
      edm::Handle<HGCalDigiCollection> theDigiContainers;
      iEvent.getByToken(tok_digiFH_, theDigiContainers);
      if (theDigiContainers.isValid()) {
#ifdef EDM_ML_DEBUG
    std::cout << "HGCDigiContainer for " << detectorFH_ << " with " 
          << theDigiContainers->size() << " element(s)" << std::endl;
#endif
    for (auto it : *theDigiContainers) {
      HGCalDetId detId     = (it.id());
      const HGCSample&  hgcSample = it.sample(sampleIndex_);
      uint16_t   adc       = hgcSample.data();
      analyzeDigi(1, detId, adc);
    }
      }
    }
  }

  //The Rechits
  if (doRecHits_) {
    edm::Handle<HGCRecHitCollection> theCaloHitContainers;
    if (ifEE_) {
      iEvent.getByToken(tok_hitrEE_, theCaloHitContainers);
      if (theCaloHitContainers.isValid()) {
#ifdef EDM_ML_DEBUG
    std::cout << "HGCRecHitCollection for " << detectorEE_ << " has "
          << theCaloHitContainers->size() << " hits" << std::endl;
#endif
    analyzeRecHits(0, theCaloHitContainers);
      } else {
#ifdef EDM_ML_DEBUG
    std::cout << "HGCRecHitCollection does not exist for " << detectorEE_ 
          << " !!!" << std::endl;
#endif
      }
    }
    if (ifFH_) {
      iEvent.getByToken(tok_hitrFH_, theCaloHitContainers);
      if (theCaloHitContainers.isValid()) {
#ifdef EDM_ML_DEBUG
    std::cout << "HGCRecHitCollection for " << detectorFH_ << " has "
          << theCaloHitContainers->size() << " hits" << std::endl;
#endif
    analyzeRecHits(1, theCaloHitContainers);
      } else {
#ifdef EDM_ML_DEBUG
    std::cout << "HGCRecHitCollection does not exist for " << detectorFH_ 
          << " !!!" << std::endl;
#endif
      }//else
    }//if (ifFH_)
  }//if (doRecHits_)

}//void HGCalTBAnalyzer::analyze

template<class T1 >

void HGCalTBAnalyzer::analyzeDigi ( int  type, const T1 &  detId, uint16_t  adc  )

private

Definition at line 923 of file HGCalTBAnalyzer.cc.

References HGCalGeometry::getPosition(), hDigiADC_, hDigiLng_, hDigiOcc_, hgeom_, and globals_cff::id1.

Referenced by analyze(), and endRun().

                                                                          {

  DetId id1 = DetId(detId.rawId());
  GlobalPoint global = hgeom_[type]->getPosition(id1);
  hDigiOcc_[type]->Fill(global.x(),global.y());
  hDigiLng_[type]->Fill(global.z());
  hDigiADC_[type]->Fill(adc);
}

void HGCalTBAnalyzer::analyzePassiveHits ( edm::Handle< edm::PassiveHitContainer > const &  hgcPh, int  subdet  )

private

Definition at line 982 of file HGCalTBAnalyzer.cc.

References gather_cfg::cout, DEFINE_FWK_MODULE, hgcPassiveBHEnergy_, hgcPassiveBHID_, hgcPassiveBHName_, hgcPassiveCMSEEnergy_, hgcPassiveCMSEID_, hgcPassiveCMSEName_, hgcPassiveEEEnergy_, hgcPassiveEEID_, hgcPassiveEEName_, hgcPassiveFHEnergy_, hgcPassiveFHID_, hgcPassiveFHName_, dataset::name, AlCaHLTBitMon_QueryRunRegistry::string, ntuplemaker::time, and findQualityFiles::v.

Referenced by analyze(), and endRun().

                                                                                                    {
  
  for (auto v : *hgcPH) {
    double       energy = v.energy();
    std::string  name   = v.vname();
    unsigned int id     = v.id();
#ifdef EDM_ML_DEBUG
    double       time   = v.time();
    std::cout << "HGCalTBAnalyzer::analyzePassiveHits:Energy:Time:Name:Id : "
          << energy << ":" << time << ":" << name << ":" << id 
          << std::endl;
#endif    

    if (subdet==1) {
      hgcPassiveEEEnergy_.push_back(energy);
      hgcPassiveEEName_.push_back(name);
      hgcPassiveEEID_.push_back(id);
    } else if (subdet==2) {
      hgcPassiveFHEnergy_.push_back(energy);
      hgcPassiveFHName_.push_back(name);
      hgcPassiveFHID_.push_back(id);
    } else if (subdet==3) {
      hgcPassiveBHEnergy_.push_back(energy);
      hgcPassiveBHName_.push_back(name);
      hgcPassiveBHID_.push_back(id);
    } else if (subdet==4) {
      hgcPassiveCMSEEnergy_.push_back(energy);
      hgcPassiveCMSEName_.push_back(name);
      hgcPassiveCMSEID_.push_back(id);
    }
  }
}

void HGCalTBAnalyzer::analyzeRecHits ( int  type, edm::Handle< HGCRecHitCollection > &  hits  )

private

Definition at line 932 of file HGCalTBAnalyzer.cc.

References HGCalDetId::cell(), gather_cfg::cout, plotBeamSpotDB::first, HGCalGeometry::getPosition(), hgcons_, hgeom_, hRecHitE_, hRecHitLat_, hRecHitLng1_, hRecHitLng_, hRecHitOcc_, HGCalDetId::layer(), edm::second(), HGCalDDDConstants::waferZ(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by analyze(), and endRun().

                                                            {
 
  std::map<int,double>                   map_hitLayer;
  std::map<int,std::pair<DetId,double> > map_hitCell;
  for (auto it : *hits) {
    DetId       detId  = it.id();
    GlobalPoint global = hgeom_[type]->getPosition(detId);
    double      energy = it.energy();
    int         layer  = HGCalDetId(detId).layer();
    int         cell   = HGCalDetId(detId).cell();
    hRecHitOcc_[type]->Fill(global.x(),global.y(),energy);
    hRecHitE_[type]->Fill(energy);
    if (map_hitLayer.count(layer) != 0) {
      map_hitLayer[layer] += energy;
    } else {
      map_hitLayer[layer]  = energy;
    }
    if (map_hitCell.count(cell) != 0) {
      double ee         = energy + map_hitCell[cell].second;
      map_hitCell[cell] = std::pair<uint32_t,double>(detId,ee);
    } else {
      map_hitCell[cell] = std::pair<uint32_t,double>(detId,energy);
    }
#ifdef EDM_ML_DEBUG
    std::cout << "RecHit: " << layer  << " " << global.x() << " " << global.y()
          << " " << global.z() << " " << energy << std::endl;
#endif
  }

  for (auto itr : map_hitLayer) {
    int    layer      = itr.first;
    double energy     = itr.second;
    double zp         = hgcons_[type]->waferZ(layer,true);
#ifdef EDM_ML_DEBUG
    std::cout << "SimHit:Layer " << layer << " " << zp << " " << energy 
          << std::endl;
#endif
    hRecHitLng_[type]->Fill(zp,energy);
    hRecHitLng1_[type]->Fill(layer,energy);
  }

  for (auto itr : map_hitCell) {
    DetId       detId  = ((itr.second).first);
    double      energy = ((itr.second).second);
    GlobalPoint global = hgeom_[type]->getPosition(detId);
    hRecHitLat_[type]->Fill(global.x(),global.y(),energy);
  }
}

void HGCalTBAnalyzer::analyzeSimHits ( int  type, std::vector< PCaloHit > &  hits, double  zFront  )

private

Definition at line 711 of file HGCalTBAnalyzer.cc.

References gather_cfg::cout, egammaForCoreTracking_cff::depth, AHCalDetId::depth(), plotBeamSpotDB::first, AHCalDetId::getXY(), AHCalDetId::getZ(), HcalOther, hgcons_, hSimHitE_, hSimHitEn_, hSimHitLat_, hSimHitLayEn1BH_, hSimHitLayEn1EE_, hSimHitLayEn1FH_, hSimHitLayEn2BH_, hSimHitLayEn2EE_, hSimHitLayEn2FH_, hSimHitLayEnBeam_, hSimHitLng1_, hSimHitLng2_, hSimHitLng_, hSimHitT_, mps_fire::i, AHCalDetId::icol(), AHCalDetId::icolAbs(), triggerObjects_cff::id, idBeams_, training_settings::idx, AHCalDetId::irow(), AHCalDetId::irowAbs(), gen::k, HGCalDDDConstants::locateCell(), HGCalTestNumbering::packHexagonIndex(), edm::second(), simHitCellEnBeam_, simHitCellEnBH_, simHitCellEnEE_, simHitCellEnFH_, simHitCellIdBeam_, simHitCellIdBH_, simHitCellIdEE_, simHitCellIdFH_, simHitLayEn1BH_, simHitLayEn1EE_, simHitLayEn1FH_, simHitLayEn2BH_, simHitLayEn2EE_, simHitLayEn2FH_, simHitLayEnBeam_, HGCalDDDConstants::simToReco(), HcalDetId::subdet(), ntuplemaker::time, HGCalTestNumbering::unpackHexagonIndex(), HcalTestBeamNumbering::unpackIndex(), HGCalDDDConstants::waferZ(), geometryCSVtoXML::xx, geometryCSVtoXML::xy, AHCalDetId::zside(), and ecaldqm::zside().

Referenced by analyze(), and endRun().

                                     {

  std::map<uint32_t,double>                 map_hits, map_hitn;
  std::map<int,double>                      map_hitDepth;
  std::map<int,std::pair<uint32_t,double> > map_hitLayer, map_hitCell;
  double                                    entot(0);
  for (unsigned int i=0; i<hits.size(); i++) {
    double energy      = hits[i].energy();
    double time        = hits[i].time();
    uint32_t id        = hits[i].id();
    entot             += energy;
    int      subdet, zside, layer, sector, subsector(0), cell, depth(0), idx(0);
    if (type == 2) {
      subdet           = HcalDetId(id).subdet();
      if (subdet != HcalOther) continue;
      AHCalDetId hid(id);
      layer = depth    = hid.depth();
      zside            = hid.zside();
      sector           = hid.irow();
      cell             = hid.icol();
      idx              = ((hid.irowAbs()*100) + (hid.icolAbs()));
    } else if (type == 3) {
      HcalTestBeamNumbering::unpackIndex(id, subdet, layer, sector, cell);
      depth = layer; zside = 1;
      idx              = subdet*1000 + layer;
      layer            = idx;
    } else {
      HGCalTestNumbering::unpackHexagonIndex(id, subdet, zside, layer, sector,
                         subsector, cell);
      depth           = hgcons_[type]->simToReco(cell,layer,sector,true).second;
      idx             = sector*1000+cell;
    }
#ifdef EDM_ML_DEBUG
    std::cout << "SimHit:Hit[" << i << "] Id " << subdet << ":" << zside << ":"
          << layer << ":" << sector << ":" << subsector << ":" << cell 
          << ":" << depth << " Energy " << energy << " Time " << time
          << std::endl;
#endif
    if (map_hits.count(id) != 0) {
      map_hits[id] += energy;
    } else {
      map_hits[id]  = energy;
    }
    if (map_hitLayer.count(layer) != 0) {
      double ee           = energy + map_hitLayer[layer].second;
      map_hitLayer[layer] = std::pair<uint32_t,double>(id,ee);
    } else {
      map_hitLayer[layer] = std::pair<uint32_t,double>(id,energy);
    }
    if (depth >= 0) {
      if (map_hitCell.count(idx) != 0) {
    double ee        = energy + map_hitCell[idx].second;
    map_hitCell[idx] = std::pair<uint32_t,double>(id,ee);
      } else {
    map_hitCell[idx] = std::pair<uint32_t,double>(id,energy);
      }
      if (map_hitDepth.count(depth) != 0) {
    map_hitDepth[depth] += energy;
      } else {
    map_hitDepth[depth]  = energy;
      }
      uint32_t idn  = (type >= 2) ? id : 
    HGCalTestNumbering::packHexagonIndex(subdet, zside, depth, sector, 
                         subsector, cell);
      if (map_hitn.count(idn) != 0) {
    map_hitn[idn] += energy;
      } else {
    map_hitn[idn]  = energy;
      }
    }
    hSimHitT_[type]->Fill(time,energy);
  }

  hSimHitEn_[type]->Fill(entot);
  for (auto itr : map_hits)  {
    hSimHitE_[type]->Fill(itr.second);
  }

  for (auto itr : map_hitLayer) {
    int    layer      = itr.first - 1;
    double energy     = (itr.second).second;
    double zp(0);
    if (type < 2)       zp = hgcons_[type]->waferZ(layer+1,false);
    else if (type == 2) zp = AHCalDetId((itr.second).first).getZ();
#ifdef EDM_ML_DEBUG
    std::cout << "SimHit:Layer " << layer+1 << " Z " << zp << ":" << zp-zFront
          << " E " << energy << std::endl;
#endif
    if (type < 3) {
      hSimHitLng_[type]->Fill(zp-zFront,energy);
      hSimHitLng2_[type]->Fill(layer+1,energy);
    }
    if (type == 0) {
      if (layer < (int)(hSimHitLayEn1EE_.size())) {
    simHitLayEn1EE_[layer] = energy;
    hSimHitLayEn1EE_[layer]->Fill(energy);
      }
    } else if (type == 1) {
      if (layer < (int)(hSimHitLayEn1FH_.size())) {
    simHitLayEn1FH_[layer] = energy;
    hSimHitLayEn1FH_[layer]->Fill(energy);
      }
    } else if (type == 2) {
      if (layer < (int)(hSimHitLayEn1BH_.size())) {
    simHitLayEn1BH_[layer] = energy;
    hSimHitLayEn1BH_[layer]->Fill(energy);
      }
    } else {
      for (unsigned int k=0; k<idBeams_.size(); ++k) {
    if (layer+1 == idBeams_[k]) {
      simHitLayEnBeam_[k] = energy;
      hSimHitLayEnBeam_[k]->Fill(energy);
      break;
    }
      }
    }
  }
  for (auto itr : map_hitDepth) {
    int    layer      = itr.first - 1;
    double energy     = itr.second;
#ifdef EDM_ML_DEBUG
    std::cout << "SimHit:Layer " << layer+1 << " " << energy << std::endl;
#endif
    hSimHitLng1_[type]->Fill(layer+1,energy);
    if (type == 0) {
      if (layer < (int)(hSimHitLayEn2EE_.size())) {
    simHitLayEn2EE_[layer] = energy;
    hSimHitLayEn2EE_[layer]->Fill(energy);
      }
    } else if (type == 1) {
      if (layer < (int)(hSimHitLayEn2FH_.size())) {
    simHitLayEn2FH_[layer] = energy;
    hSimHitLayEn2FH_[layer]->Fill(energy);
      }
    } else if (type == 2) {
      if (layer < (int)(hSimHitLayEn2BH_.size())) {
    simHitLayEn2BH_[layer] = energy;
    hSimHitLayEn2BH_[layer]->Fill(energy);
      }
    }
  }

  if (type < 3) {
    for (auto itr : map_hitCell) {
      uint32_t id       = ((itr.second).first);
      double   energy   = ((itr.second).second);
      std::pair<float,float> xy(0,0);
      double xx(0);
      if (type == 2) {
    xy = AHCalDetId(id).getXY();
    xx = xy.first;
      } else {
    int      subdet, zside, layer, sector, subsector, cell;
    HGCalTestNumbering::unpackHexagonIndex(id, subdet, zside, layer, sector,
                           subsector, cell);
    xy        = hgcons_[type]->locateCell(cell,layer,sector,false);
    double zp = hgcons_[type]->waferZ(layer,false);
    xx        = (zp < 0) ? -xy.first : xy.first;
      }
      hSimHitLat_[type]->Fill(xx,xy.second,energy);
    }
  }

  for (auto itr : map_hitn) {
    uint32_t id     = itr.first;
    double   energy = itr.second;
    if (type == 0) {
      simHitCellIdEE_.push_back(id); simHitCellEnEE_.push_back(energy);
    } else if (type == 1) {
      simHitCellIdFH_.push_back(id); simHitCellEnFH_.push_back(energy);
    } else if (type == 2) {
      simHitCellIdBH_.push_back(id); simHitCellEnBH_.push_back(energy);
    } else if (type == 3) {
      simHitCellIdBeam_.push_back(id); simHitCellEnBeam_.push_back(energy);
    }
  }
}

void HGCalTBAnalyzer::analyzeSimTracks ( edm::Handle< edm::SimTrackContainer > const &  SimTk, edm::Handle< edm::SimVertexContainer > const &  SimVtx  )

private

Definition at line 890 of file HGCalTBAnalyzer.cc.

References gather_cfg::cout, pBeam_, xBeam_, yBeam_, and zBeam_.

Referenced by analyze(), and endRun().

                                                                       {

  xBeam_ = yBeam_ = zBeam_ = pBeam_ = -1000000;
  int vertIndex(-1);
  for (auto simTrkItr : *SimTk) {
#ifdef EDM_ML_DEBUG
    std::cout << "Track " << simTrkItr.trackId() << " Vertex "
          << simTrkItr.vertIndex() << " Type " << simTrkItr.type()
          << " Charge " << simTrkItr.charge() << " momentum "
          << simTrkItr.momentum() << " " << simTrkItr.momentum().P()
          << std::endl;
#endif
    if (vertIndex == -1) {
      vertIndex = simTrkItr.vertIndex();
      pBeam_    = simTrkItr.momentum().P();
    }
  }
  if (vertIndex != -1 && vertIndex < (int)SimVtx->size()) {
    edm::SimVertexContainer::const_iterator simVtxItr= SimVtx->begin();
    for (int iv=0; iv<vertIndex; iv++) simVtxItr++;
#ifdef EDM_ML_DEBUG
    std::cout << "Vertex " << vertIndex << " position "
          << simVtxItr->position() << std::endl;
#endif
    xBeam_ = simVtxItr->position().X();
    yBeam_ = simVtxItr->position().Y();
    zBeam_ = simVtxItr->position().Z();
  }

}

void HGCalTBAnalyzer::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 274 of file HGCalTBAnalyzer.cc.

References detectorBeam_, detectorBH_, detectorEE_, detectorFH_, doDigis_, doPassive_, doRecHits_, doSimHits_, doTree_, doTreeCell_, fs_, hBeam_, hDigiADC_, hDigiLng_, hDigiOcc_, hgcPassiveBHEnergy_, hgcPassiveBHID_, hgcPassiveBHName_, hgcPassiveCMSEEnergy_, hgcPassiveCMSEID_, hgcPassiveCMSEName_, hgcPassiveEEEnergy_, hgcPassiveEEID_, hgcPassiveEEName_, hgcPassiveFHEnergy_, hgcPassiveFHID_, hgcPassiveFHName_, hRecHitE_, hRecHitLat_, hRecHitLng1_, hRecHitLng_, hRecHitOcc_, hSimHitE_, hSimHitEn_, hSimHitLat_, hSimHitLng1_, hSimHitLng2_, hSimHitLng_, hSimHitT_, mps_fire::i, ifBeam_, ifBH_, ifEE_, ifFH_, TFileService::make(), dataset::name, pBeam_, simHitCellEnBeam_, simHitCellEnBH_, simHitCellEnEE_, simHitCellEnFH_, simHitCellIdBeam_, simHitCellIdBH_, simHitCellIdEE_, simHitCellIdFH_, simHitLayEn1BH_, simHitLayEn1EE_, simHitLayEn1FH_, simHitLayEn2BH_, simHitLayEn2EE_, simHitLayEn2FH_, AlCaHLTBitMon_QueryRunRegistry::string, fftjetcommon_cfi::title, tree_, xBeam_, yBeam_, and zBeam_.

                               {

  char name[40], title[100];
  hBeam_ = fs_->make<TH1D>("BeamP", "Beam Momentum", 1000, 0, 1000.0);
  for (int i=0; i<3; ++i) {
    bool book(ifEE_);
    std::string det(detectorEE_); 
    if (i == 1) {
      book = ifFH_;
      det  = detectorFH_;
    } else if (i == 2) {
      book = ifBH_;
      det  = detectorBH_;
    }

    if (doSimHits_ && book) {
      sprintf (name, "SimHitEn%s", det.c_str());
      sprintf (title,"Sim Hit Energy for %s", det.c_str());
      hSimHitE_[i] = fs_->make<TH1D>(name,title,100000,0.,0.2);
      sprintf (name, "SimHitEnX%s", det.c_str());
      sprintf (title,"Sim Hit Energy for %s", det.c_str());
      hSimHitEn_[i] = fs_->make<TH1D>(name,title,100000,0.,0.2);
      sprintf (name, "SimHitTm%s", det.c_str());
      sprintf (title,"Sim Hit Timing for %s", det.c_str());
      hSimHitT_[i] = fs_->make<TH1D>(name,title,5000,0.,500.0);
      sprintf (name, "SimHitLat%s", det.c_str());
      sprintf (title,"Lateral Shower profile (Sim Hit) for %s", det.c_str());
      hSimHitLat_[i] = fs_->make<TProfile2D>(name,title,100,-100.,100.,100,-100.,100.);
      sprintf (name, "SimHitLng%s", det.c_str());
      sprintf (title,"Longitudinal Shower profile (Sim Hit) for %s",det.c_str());
      hSimHitLng_[i] = fs_->make<TProfile>(name,title,50,0.,100.);
      sprintf (name, "SimHitLng1%s", det.c_str());
      sprintf (title,"Longitudinal Shower profile (Layer) for %s",det.c_str());
      hSimHitLng1_[i] = fs_->make<TProfile>(name,title,200,0.,100.);
      sprintf (name, "SimHitLng2%s", det.c_str());
      sprintf (title,"Longitudinal Shower profile (Layer) for %s",det.c_str());
      hSimHitLng2_[i] = fs_->make<TProfile>(name,title,200,0.,100.);
    }

    if (doDigis_ && book) {
      sprintf (name, "DigiADC%s", det.c_str());
      sprintf (title,"ADC at Digi level for %s", det.c_str());
      hDigiADC_[i] = fs_->make<TH1D>(name,title,100,0.,100.0);
      sprintf (name, "DigiOcc%s", det.c_str());
      sprintf (title,"Occupancy (Digi)for %s", det.c_str());
      hDigiOcc_[i] = fs_->make<TH2D>(name,title,100,-10.,10.,100,-10.,10.);
      sprintf (name, "DigiLng%s", det.c_str());
      sprintf (title,"Longitudinal Shower profile (Digi) for %s",det.c_str());
      hDigiLng_[i] = fs_->make<TH1D>(name,title,100,0.,10.);
    }

    if (doRecHits_ && book) {
      sprintf (name, "RecHitEn%s", det.c_str());
      sprintf (title,"Rec Hit Energy for %s", det.c_str());
      hRecHitE_[i] = fs_->make<TH1D>(name,title,1000,0.,10.0);
      sprintf (name, "RecHitOcc%s", det.c_str());
      sprintf (title,"Occupancy (Rec Hit)for %s", det.c_str());
      hRecHitOcc_[i] = fs_->make<TH2D>(name,title,100,-10.,10.,100,-10.,10.);
      sprintf (name, "RecHitLat%s", det.c_str());
      sprintf (title,"Lateral Shower profile (Rec Hit) for %s", det.c_str());
      hRecHitLat_[i] = fs_->make<TProfile2D>(name,title,100,-10.,10.,100,-10.,10.);
      sprintf (name, "RecHitLng%s", det.c_str());
      sprintf (title,"Longitudinal Shower profile (Rec Hit) for %s",det.c_str());
      hRecHitLng_[i] = fs_->make<TProfile>(name,title,100,0.,10.);
      sprintf (name, "RecHitLng1%s", det.c_str());
      sprintf (title,"Longitudinal Shower profile vs Layer for %s",det.c_str());
      hRecHitLng1_[i] = fs_->make<TProfile>(name,title,120,0.,60.);
    }
  }
  if (ifBeam_ && doSimHits_) {
    sprintf (name, "SimHitEn%s", detectorBeam_.c_str());
    sprintf (title,"Sim Hit Energy for %s", detectorBeam_.c_str());
    hSimHitE_[3] = fs_->make<TH1D>(name,title,100000,0.,0.2);
    sprintf (name, "SimHitEnX%s", detectorBeam_.c_str());
    sprintf (title,"Sim Hit Energy for %s", detectorBeam_.c_str());
    hSimHitEn_[3] = fs_->make<TH1D>(name,title,100000,0.,0.2);
    sprintf (name, "SimHitTm%s", detectorBeam_.c_str());
    sprintf (title,"Sim Hit Timing for %s", detectorBeam_.c_str());
    hSimHitT_[3] = fs_->make<TH1D>(name,title,5000,0.,500.0);
  }
  if (doSimHits_ && doTree_) {
    tree_ = fs_->make<TTree>("HGCTB","SimHitEnergy");
    tree_->Branch("simHitLayEn1EE", &simHitLayEn1EE_);
    tree_->Branch("simHitLayEn2EE", &simHitLayEn2EE_);
    tree_->Branch("simHitLayEn1FH", &simHitLayEn1FH_);
    tree_->Branch("simHitLayEn2FH", &simHitLayEn2FH_);
    tree_->Branch("simHitLayEn1BH", &simHitLayEn1BH_);
    tree_->Branch("simHitLayEn2BH", &simHitLayEn2BH_);
    tree_->Branch("xBeam",          &xBeam_,           "xBeam/D");
    tree_->Branch("yBeam",          &yBeam_,           "yBeam/D");
    tree_->Branch("zBeam",          &zBeam_,           "zBeam/D");
    tree_->Branch("pBeam",          &pBeam_,           "pBeam/D");
    if (doTreeCell_) {
      tree_->Branch("simHitCellIdEE", &simHitCellIdEE_);
      tree_->Branch("simHitCellEnEE", &simHitCellEnEE_);
      tree_->Branch("simHitCellIdFH", &simHitCellIdFH_);
      tree_->Branch("simHitCellEnFH", &simHitCellEnFH_);
      tree_->Branch("simHitCellIdBH", &simHitCellIdBH_);
      tree_->Branch("simHitCellEnBH", &simHitCellEnBH_);
      tree_->Branch("simHitCellIdBeam", &simHitCellIdBeam_);
      tree_->Branch("simHitCellEnBeam", &simHitCellEnBeam_);
    }
  }

  if (doPassive_ && doTree_) {
    tree_->Branch("hgcPassiveEEEnergy",   &hgcPassiveEEEnergy_);
    tree_->Branch("hgcPassiveEEName",     &hgcPassiveEEName_);
    tree_->Branch("hgcPassiveEEID",       &hgcPassiveEEID_);
    tree_->Branch("hgcPassiveFHEnergy",   &hgcPassiveFHEnergy_);
    tree_->Branch("hgcPassiveFHName",     &hgcPassiveFHName_);
    tree_->Branch("hgcPassiveFHID",       &hgcPassiveFHID_);
    tree_->Branch("hgcPassiveBHEnergy",   &hgcPassiveBHEnergy_);
    tree_->Branch("hgcPassiveBHName",     &hgcPassiveBHName_);
    tree_->Branch("hgcPassiveBHID",       &hgcPassiveBHID_);
    tree_->Branch("hgcPassiveCMSEEnergy", &hgcPassiveCMSEEnergy_);
    tree_->Branch("hgcPassiveCMSEName",   &hgcPassiveCMSEName_);
    tree_->Branch("hgcPassiveCMSEID",     &hgcPassiveCMSEID_);
  }
}

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

overrideprivate

Definition at line 394 of file HGCalTBAnalyzer.cc.

References gather_cfg::cout, detectorBeam_, detectorBH_, detectorEE_, detectorFH_, doDigis_, doRecHits_, fs_, relativeConstraints::geom, edm::EventSetup::get(), hgcons_, hgeom_, hSimHitLayEn1BH_, hSimHitLayEn1EE_, hSimHitLayEn1FH_, hSimHitLayEn2BH_, hSimHitLayEn2EE_, hSimHitLayEn2FH_, hSimHitLayEnBeam_, idBeams_, ifBeam_, ifBH_, ifEE_, ifFH_, checklumidiff::l, HGCalDDDConstants::layers(), TFileService::make(), AHCalDetId::MaxDepth, dataset::name, edm::ESHandle< T >::product(), and fftjetcommon_cfi::title.

                                                                         {

  char name[40], title[100];
  if (ifEE_) {
    edm::ESHandle<HGCalDDDConstants>  pHGDC;
    iSetup.get<IdealGeometryRecord>().get(detectorEE_, pHGDC);
    hgcons_[0] = &(*pHGDC);
    if (doDigis_ || doRecHits_) {
      edm::ESHandle<HGCalGeometry> geom;
      iSetup.get<IdealGeometryRecord>().get(detectorEE_, geom);
      hgeom_[0] = geom.product();
    } else {
      hgeom_[0] = nullptr;
    }
    for (unsigned int l=0; l<hgcons_[0]->layers(false); ++l) {
      sprintf (name, "SimHitEnA%d%s", l, detectorEE_.c_str());
      sprintf (title,"Sim Hit Energy in SIM layer %d for %s",l+1,
           detectorEE_.c_str());
      hSimHitLayEn1EE_.push_back(fs_->make<TH1D>(name,title,100000,0.,0.2));
      if (l%3 == 0) {
    sprintf (name, "SimHitEnB%d%s", (l/3+1), detectorEE_.c_str());
    sprintf (title,"Sim Hit Energy in layer %d for %s",(l/3+1),
         detectorEE_.c_str());
    hSimHitLayEn2EE_.push_back(fs_->make<TH1D>(name,title,100000,0.,0.2));
      }
    }
#ifdef EDM_ML_DEBUG
    std::cout << "HGCalTBAnalyzer::" << detectorEE_ << " defined with "
          << hgcons_[0]->layers(false) << " layers" << std::endl;
#endif
  } else {
    hgcons_[0] = nullptr;
    hgeom_[0]  = nullptr;
  }

  if (ifFH_) {
    edm::ESHandle<HGCalDDDConstants>  pHGDC;
    iSetup.get<IdealGeometryRecord>().get(detectorFH_, pHGDC);
    hgcons_[1] = &(*pHGDC);
    if (doDigis_ || doRecHits_) {
      edm::ESHandle<HGCalGeometry> geom;
      iSetup.get<IdealGeometryRecord>().get(detectorFH_, geom);
      hgeom_[1] = geom.product();
    } else {
      hgeom_[1] = nullptr;
    }
    for (unsigned int l=0; l<hgcons_[1]->layers(false); ++l) {
      sprintf (name, "SimHitEnA%d%s", l, detectorFH_.c_str());
      sprintf (title,"Sim Hit Energy in layer %d for %s",l+1,
           detectorFH_.c_str());
      hSimHitLayEn1FH_.push_back(fs_->make<TH1D>(name,title,100000,0.,0.2));
      if (l%3 == 0) {
    sprintf (name, "SimHitEnB%d%s", (l/3+1), detectorFH_.c_str());
    sprintf (title,"Sim Hit Energy in layer %d for %s",(l/3+1),
         detectorFH_.c_str());
    hSimHitLayEn2FH_.push_back(fs_->make<TH1D>(name,title,100000,0.,0.2));
      }
    }
#ifdef EDM_ML_DEBUG
    std::cout << "HGCalTBAnalyzer::" << detectorFH_ << " defined with "
          << hgcons_[1]->layers(false) << " layers" << std::endl;
#endif
  } else {
    hgcons_[1] = nullptr;
    hgeom_[1]  = nullptr;
  }

  if (ifBH_) {
    for (int l=0; l<AHCalDetId::MaxDepth; ++l) {
      sprintf (name, "SimHitEnA%d%s", l, detectorBH_.c_str());
      sprintf (title,"Sim Hit Energy in layer %d for %s",l+1,
           detectorBH_.c_str());
      hSimHitLayEn1BH_.push_back(fs_->make<TH1D>(name,title,100000,0.,0.2));
      sprintf (name, "SimHitEnB%d%s", l, detectorBH_.c_str());
      sprintf (title,"Sim Hit Energy in layer %d for %s",l+1,
           detectorBH_.c_str());
      hSimHitLayEn2BH_.push_back(fs_->make<TH1D>(name,title,100000,0.,0.2));
    }
  }

  if (ifBeam_) {
    for (unsigned int l=0; l<idBeams_.size(); ++l) {
      sprintf (name, "SimHitEna%d%s", l, detectorBeam_.c_str());
      sprintf (title, "Sim Hit Energy in type %d for %s",idBeams_[l],
           detectorBeam_.c_str());
      hSimHitLayEnBeam_.push_back(fs_->make<TH1D>(name,title,100000,0.,0.2));
    }
  }
}

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

inlineoverrideprivate

Definition at line 64 of file HGCalTBAnalyzer.cc.

References ecalMGPA::adc(), analyze(), analyzeDigi(), analyzePassiveHits(), analyzeRecHits(), analyzeSimHits(), analyzeSimTracks(), and hfClusterShapes_cfi::hits.

{}

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

static

Definition at line 231 of file HGCalTBAnalyzer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addUntracked(), edm::ParameterSetDescription::setUnknown(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                 {
  edm::ParameterSetDescription desc;
  desc.setUnknown();
  desc.add<std::string>("DetectorEE","HGCalEESensitive");
  desc.add<bool>("UseEE",true);
  desc.add<double>("ZFrontEE",0.0);
  desc.add<std::string>("CaloHitSrcEE","HGCHitsEE");
  desc.add<edm::InputTag>("DigiSrcEE",edm::InputTag("hgcalDigis","EE"));
  desc.add<edm::InputTag>("RecHitSrcEE",edm::InputTag("HGCalRecHit","HGCEERecHits"));
  desc.add<std::string>("DetectorFH","HGCalHESiliconSensitive");
  desc.add<bool>("UseFH",false);
  desc.add<double>("ZFrontFH",0.0);
  desc.add<std::string>("CaloHitSrcFH","HGCHitsHEfront");
  desc.add<edm::InputTag>("DigiSrcFH",edm::InputTag("hgcalDigis","HEfront"));
  desc.add<edm::InputTag>("RecHitSrcFH",edm::InputTag("HGCalRecHit","HGCHEFRecHits"));
  desc.add<std::string>("DetectorBH","AHCal");
  desc.add<bool>("UseBH",false);
  desc.add<double>("ZFrontBH",0.0);
  desc.add<std::string>("CaloHitSrcBH","HcalHits");
  desc.add<edm::InputTag>("DigiSrcBH",edm::InputTag("hgcalDigis","HEback"));
  desc.add<edm::InputTag>("RecHitSrcBH",edm::InputTag("HGCalRecHit","HGCHEBRecHits"));
  desc.add<std::string>("DetectorBeam","HcalTB06BeamDetector");
  desc.add<bool>("UseBeam",false);
  desc.add<std::string>("CaloHitSrcBeam","HcalTB06BeamHits");
  std::vector<int> ids = {1000,1001,1002,1003,1004,1005,1006,1007,1008,1011,1012,1013,1014,2001,2002,2003,2004,2005};
  desc.add<std::vector<int>>("IDBeams",ids);
  desc.add<edm::InputTag>("GeneratorSrc",edm::InputTag("generatorSmeared"));
  desc.add<edm::InputTag>("HGCPassiveEE",edm::InputTag("g4SimHits","HGCalEEPassiveHits"));
  desc.add<edm::InputTag>("HGCPassiveFH",edm::InputTag("g4SimHits","HGCalHEPassiveHits"));
  desc.add<edm::InputTag>("HGCPassiveBH",edm::InputTag("g4SimHits","HGCalAHPassiveHits"));
  desc.add<edm::InputTag>("HGCPassiveCMSE",edm::InputTag("g4SimHits","CMSEPassiveHits"));

  desc.add<bool>("DoSimHits",true);
  desc.add<bool>("DoDigis",true);
  desc.add<bool>("DoRecHits",true);
  desc.add<int>("SampleIndex",0);
  desc.addUntracked<bool>("DoTree",true);
  desc.addUntracked<bool>("DoTreeCell",true);
  desc.addUntracked<bool>("DoPassive",false);

  descriptions.add("HGCalTBAnalyzer",desc);
}

Member Data Documentation

std::string HGCalTBAnalyzer::detectorBeam_

private

Definition at line 82 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), beginJob(), beginRun(), and HGCalTBAnalyzer().

std::string HGCalTBAnalyzer::detectorBH_

private

Definition at line 82 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), beginJob(), beginRun(), and HGCalTBAnalyzer().

std::string HGCalTBAnalyzer::detectorEE_

private

Definition at line 81 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), beginJob(), beginRun(), and HGCalTBAnalyzer().

std::string HGCalTBAnalyzer::detectorFH_

private

Definition at line 81 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), beginJob(), beginRun(), and HGCalTBAnalyzer().

bool HGCalTBAnalyzer::doDigis_

private

Definition at line 79 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), beginJob(), beginRun(), and HGCalTBAnalyzer().

bool HGCalTBAnalyzer::doPassive_

private

Definition at line 80 of file HGCalTBAnalyzer.cc.

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

bool HGCalTBAnalyzer::doRecHits_

private

Definition at line 79 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), beginJob(), beginRun(), and HGCalTBAnalyzer().

bool HGCalTBAnalyzer::doSimHits_

private

Definition at line 79 of file HGCalTBAnalyzer.cc.

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

bool HGCalTBAnalyzer::doTree_

private

Definition at line 78 of file HGCalTBAnalyzer.cc.

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

bool HGCalTBAnalyzer::doTreeCell_

private

Definition at line 78 of file HGCalTBAnalyzer.cc.

Referenced by beginJob(), and HGCalTBAnalyzer().

edm::Service<TFileService> HGCalTBAnalyzer::fs_

private

Definition at line 74 of file HGCalTBAnalyzer.cc.

Referenced by beginJob(), and beginRun().

TH1D * HGCalTBAnalyzer::hBeam_

private

Definition at line 99 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and beginJob().

TH1D* HGCalTBAnalyzer::hDigiADC_[3]

private

Definition at line 98 of file HGCalTBAnalyzer.cc.

Referenced by analyzeDigi(), and beginJob().

TH1D * HGCalTBAnalyzer::hDigiLng_[2]

private

Definition at line 98 of file HGCalTBAnalyzer.cc.

Referenced by analyzeDigi(), and beginJob().

TH2D* HGCalTBAnalyzer::hDigiOcc_[3]

private

Definition at line 100 of file HGCalTBAnalyzer.cc.

Referenced by analyzeDigi(), and beginJob().

const HGCalDDDConstants* HGCalTBAnalyzer::hgcons_[2]

private

Definition at line 75 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), analyzeRecHits(), analyzeSimHits(), and beginRun().

std::vector<float> HGCalTBAnalyzer::hgcPassiveBHEnergy_

private

Definition at line 118 of file HGCalTBAnalyzer.cc.

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

std::vector<int> HGCalTBAnalyzer::hgcPassiveBHID_

private

Definition at line 120 of file HGCalTBAnalyzer.cc.

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

std::vector<std::string> HGCalTBAnalyzer::hgcPassiveBHName_

private

Definition at line 119 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::hgcPassiveCMSEEnergy_

private

Definition at line 118 of file HGCalTBAnalyzer.cc.

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

std::vector<int> HGCalTBAnalyzer::hgcPassiveCMSEID_

private

Definition at line 120 of file HGCalTBAnalyzer.cc.

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

std::vector<std::string> HGCalTBAnalyzer::hgcPassiveCMSEName_

private

Definition at line 119 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::hgcPassiveEEEnergy_

private

Definition at line 118 of file HGCalTBAnalyzer.cc.

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

std::vector<int> HGCalTBAnalyzer::hgcPassiveEEID_

private

Definition at line 120 of file HGCalTBAnalyzer.cc.

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

std::vector<std::string> HGCalTBAnalyzer::hgcPassiveEEName_

private

Definition at line 119 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::hgcPassiveFHEnergy_

private

Definition at line 118 of file HGCalTBAnalyzer.cc.

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

std::vector<int> HGCalTBAnalyzer::hgcPassiveFHID_

private

Definition at line 120 of file HGCalTBAnalyzer.cc.

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

std::vector<std::string> HGCalTBAnalyzer::hgcPassiveFHName_

private

Definition at line 119 of file HGCalTBAnalyzer.cc.

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

const HGCalGeometry* HGCalTBAnalyzer::hgeom_[2]

private

Definition at line 76 of file HGCalTBAnalyzer.cc.

Referenced by analyzeDigi(), analyzeRecHits(), and beginRun().

TH1D* HGCalTBAnalyzer::hRecHitE_[3]

private

Definition at line 99 of file HGCalTBAnalyzer.cc.

Referenced by analyzeRecHits(), and beginJob().

TProfile2D * HGCalTBAnalyzer::hRecHitLat_[3]

private

Definition at line 104 of file HGCalTBAnalyzer.cc.

Referenced by analyzeRecHits(), and beginJob().

TProfile * HGCalTBAnalyzer::hRecHitLng1_[3]

private

Definition at line 103 of file HGCalTBAnalyzer.cc.

Referenced by analyzeRecHits(), and beginJob().

TProfile* HGCalTBAnalyzer::hRecHitLng_[3]

private

Definition at line 103 of file HGCalTBAnalyzer.cc.

Referenced by analyzeRecHits(), and beginJob().

TH2D * HGCalTBAnalyzer::hRecHitOcc_[3]

private

Definition at line 100 of file HGCalTBAnalyzer.cc.

Referenced by analyzeRecHits(), and beginJob().

TH1D* HGCalTBAnalyzer::hSimHitE_[4]

private

Definition at line 97 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

TH1D * HGCalTBAnalyzer::hSimHitEn_[4]

private

Definition at line 99 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

TProfile2D* HGCalTBAnalyzer::hSimHitLat_[3]

private

Definition at line 104 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

std::vector<TH1D*> HGCalTBAnalyzer::hSimHitLayEn1BH_

private

Definition at line 107 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

std::vector<TH1D*> HGCalTBAnalyzer::hSimHitLayEn1EE_

private

Definition at line 105 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

std::vector<TH1D*> HGCalTBAnalyzer::hSimHitLayEn1FH_

private

Definition at line 106 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

std::vector<TH1D*> HGCalTBAnalyzer::hSimHitLayEn2BH_

private

Definition at line 107 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

std::vector<TH1D*> HGCalTBAnalyzer::hSimHitLayEn2EE_

private

Definition at line 105 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

std::vector<TH1D*> HGCalTBAnalyzer::hSimHitLayEn2FH_

private

Definition at line 106 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

std::vector<TH1D*> HGCalTBAnalyzer::hSimHitLayEnBeam_

private

Definition at line 108 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

TProfile * HGCalTBAnalyzer::hSimHitLng1_[3]

private

Definition at line 101 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

TProfile* HGCalTBAnalyzer::hSimHitLng2_[3]

private

Definition at line 102 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

TProfile* HGCalTBAnalyzer::hSimHitLng_[3]

private

Definition at line 101 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

TH1D * HGCalTBAnalyzer::hSimHitT_[4]

private

Definition at line 97 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

std::vector<int> HGCalTBAnalyzer::idBeams_

private

Definition at line 85 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), analyzeSimHits(), beginRun(), and HGCalTBAnalyzer().

bool HGCalTBAnalyzer::ifBeam_

private

Definition at line 77 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), beginJob(), beginRun(), and HGCalTBAnalyzer().

bool HGCalTBAnalyzer::ifBH_

private

Definition at line 77 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), beginJob(), beginRun(), and HGCalTBAnalyzer().

bool HGCalTBAnalyzer::ifEE_

private

Definition at line 77 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), beginJob(), beginRun(), and HGCalTBAnalyzer().

bool HGCalTBAnalyzer::ifFH_

private

Definition at line 77 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), beginJob(), beginRun(), and HGCalTBAnalyzer().

double HGCalTBAnalyzer::pBeam_

private

Definition at line 122 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

int HGCalTBAnalyzer::sampleIndex_

private

Definition at line 84 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

std::vector<float> HGCalTBAnalyzer::simHitCellEnBeam_

private

Definition at line 116 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::simHitCellEnBH_

private

Definition at line 116 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::simHitCellEnEE_

private

Definition at line 115 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::simHitCellEnFH_

private

Definition at line 115 of file HGCalTBAnalyzer.cc.

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

std::vector<uint32_t> HGCalTBAnalyzer::simHitCellIdBeam_

private

Definition at line 114 of file HGCalTBAnalyzer.cc.

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

std::vector<uint32_t> HGCalTBAnalyzer::simHitCellIdBH_

private

Definition at line 114 of file HGCalTBAnalyzer.cc.

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

std::vector<uint32_t> HGCalTBAnalyzer::simHitCellIdEE_

private

Definition at line 113 of file HGCalTBAnalyzer.cc.

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

std::vector<uint32_t> HGCalTBAnalyzer::simHitCellIdFH_

private

Definition at line 113 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::simHitLayEn1BH_

private

Definition at line 111 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::simHitLayEn1EE_

private

Definition at line 109 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::simHitLayEn1FH_

private

Definition at line 110 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::simHitLayEn2BH_

private

Definition at line 111 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::simHitLayEn2EE_

private

Definition at line 109 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::simHitLayEn2FH_

private

Definition at line 110 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::simHitLayEnBeam_

private

Definition at line 112 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and analyzeSimHits().

edm::EDGetToken HGCalTBAnalyzer::tok_digiBH_

private

Definition at line 90 of file HGCalTBAnalyzer.cc.

Referenced by HGCalTBAnalyzer().

edm::EDGetToken HGCalTBAnalyzer::tok_digiEE_

private

Definition at line 90 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

edm::EDGetToken HGCalTBAnalyzer::tok_digiFH_

private

Definition at line 90 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

edm::EDGetTokenT<edm::HepMCProduct> HGCalTBAnalyzer::tok_hepMC_

private

Definition at line 92 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

edm::EDGetTokenT<edm::PassiveHitContainer> HGCalTBAnalyzer::tok_hgcPHBH_

private

Definition at line 94 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

edm::EDGetTokenT<edm::PassiveHitContainer> HGCalTBAnalyzer::tok_hgcPHCMSE_

private

Definition at line 94 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

edm::EDGetTokenT<edm::PassiveHitContainer> HGCalTBAnalyzer::tok_hgcPHEE_

private

Definition at line 93 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

edm::EDGetTokenT<edm::PassiveHitContainer> HGCalTBAnalyzer::tok_hgcPHFH_

private

Definition at line 93 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

edm::EDGetToken HGCalTBAnalyzer::tok_hitrBH_

private

Definition at line 91 of file HGCalTBAnalyzer.cc.

Referenced by HGCalTBAnalyzer().

edm::EDGetToken HGCalTBAnalyzer::tok_hitrEE_

private

Definition at line 91 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

edm::EDGetToken HGCalTBAnalyzer::tok_hitrFH_

private

Definition at line 91 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

edm::EDGetTokenT<edm::PCaloHitContainer> HGCalTBAnalyzer::tok_hitsBeam_

private

Definition at line 87 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

edm::EDGetTokenT<edm::PCaloHitContainer> HGCalTBAnalyzer::tok_hitsBH_

private

Definition at line 87 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

edm::EDGetTokenT<edm::PCaloHitContainer> HGCalTBAnalyzer::tok_hitsEE_

private

Definition at line 86 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

edm::EDGetTokenT<edm::PCaloHitContainer> HGCalTBAnalyzer::tok_hitsFH_

private

Definition at line 86 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

edm::EDGetTokenT<edm::SimTrackContainer> HGCalTBAnalyzer::tok_simTk_

private

Definition at line 88 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

edm::EDGetTokenT<edm::SimVertexContainer> HGCalTBAnalyzer::tok_simVtx_

private

Definition at line 89 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

TTree* HGCalTBAnalyzer::tree_

private

Definition at line 96 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HGCalTBAnalyzer::xBeam_

private

Definition at line 122 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

double HGCalTBAnalyzer::yBeam_

private

Definition at line 122 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

double HGCalTBAnalyzer::zBeam_

private

Definition at line 122 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

double HGCalTBAnalyzer::zFrontBH_

private

Definition at line 83 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

double HGCalTBAnalyzer::zFrontEE_

private

Definition at line 83 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

double HGCalTBAnalyzer::zFrontFH_

private

Definition at line 83 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().