Inheritance diagram for HGCalTBAnalyzer:

Public Member Functions
	HGCalTBAnalyzer (edm::ParameterSet const &)

	~HGCalTBAnalyzer ()

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::SharedResources >
	EDAnalyzer ()=default

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	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 &	itemsToGetFromEvent () 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

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
virtual void	analyze (edm::Event const &, edm::EventSetup const &) override

template<class T1 >
void	analyzeDigi (int type, const T1 &detId, uint16_t adc)

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

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

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

virtual void	beginJob () override

virtual void	beginRun (edm::Run const &, edm::EventSetup const &) override

virtual void	endRun (edm::Run const &, edm::EventSetup const &) override

Private Attributes
std::string	detectorEE_

std::string	detectorHE_

bool	doDigis_

bool	doRecHits_

bool	doSimHits_

bool	doTree_

bool	doTreeCell_

edm::Service< TFileService >	fs_

TH1D *	hBeam_

TH1D *	hDigiADC_ [2]

TH1D *	hDigiLng_ [2]

TH2D *	hDigiOcc_ [2]

const HGCalDDDConstants *	hgcons_ [2]

const HGCalGeometry *	hgeom_ [2]

TH1D *	hRecHitE_ [2]

TProfile2D *	hRecHitLat_ [2]

TProfile *	hRecHitLng1_ [2]

TProfile *	hRecHitLng_ [2]

TH2D *	hRecHitOcc_ [2]

TH1D *	hSimHitE_ [2]

TH1D *	hSimHitEn_ [2]

TProfile2D *	hSimHitLat_ [2]

std::vector< TH1D * >	hSimHitLayEn1E_

std::vector< TH1D * >	hSimHitLayEn1H_

std::vector< TH1D * >	hSimHitLayEn2E_

std::vector< TH1D * >	hSimHitLayEn2H_

TProfile *	hSimHitLng1_ [2]

TProfile *	hSimHitLng2_ [2]

TProfile *	hSimHitLng_ [2]

TH1D *	hSimHitT_ [2]

bool	ifEE_

bool	ifHE_

double	pBeam

int	sampleIndex_

std::vector< float >	simHitCellEnE

std::vector< float >	simHitCellEnH

std::vector< uint32_t >	simHitCellIdE

std::vector< uint32_t >	simHitCellIdH

std::vector< float >	simHitLayEn1E

std::vector< float >	simHitLayEn1H

std::vector< float >	simHitLayEn2E

std::vector< float >	simHitLayEn2H

edm::EDGetToken	tok_digiEE_

edm::EDGetToken	tok_digiHE_

edm::EDGetTokenT < edm::HepMCProduct >	tok_hepMC_

edm::EDGetToken	tok_hitrEE_

edm::EDGetToken	tok_hitrHE_

edm::EDGetTokenT < edm::PCaloHitContainer >	tok_hitsEE_

edm::EDGetTokenT < edm::PCaloHitContainer >	tok_hitsHE_

edm::EDGetTokenT < edm::SimTrackContainer >	tok_simTk_

edm::EDGetTokenT < edm::SimVertexContainer >	tok_simVtx_

TTree *	tree_

double	xBeam

double	yBeam

double	zBeam

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 49 of file HGCalTBAnalyzer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 100 of file HGCalTBAnalyzer.cc.

References gather_cfg::cout, detectorEE_, detectorHE_, doDigis_, doRecHits_, doSimHits_, doTree_, doTreeCell_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), ifEE_, ifHE_, HLT_25ns10e33_v2_cff::InputTag, sampleIndex_, AlCaHLTBitMon_QueryRunRegistry::string, tok_digiEE_, tok_digiHE_, tok_hepMC_, tok_hitrEE_, tok_hitrHE_, tok_hitsEE_, tok_hitsHE_, tok_simTk_, tok_simVtx_, and usesResource().

                                                                {
 
   usesResource("TFileService");
 
   //now do whatever initialization is needed
   detectorEE_  = iConfig.getParameter<std::string>("DetectorEE");
   detectorHE_  = iConfig.getParameter<std::string>("DetectorHE");
   ifEE_        = iConfig.getParameter<bool>("UseEE");
   ifHE_        = iConfig.getParameter<bool>("UseHE");
   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);
 #ifdef DebugLog
   std::cout << "HGCalTBAnalyzer:: SimHits = " << doSimHits_ << " Digis = "
         << doDigis_ << ":" << sampleIndex_ << " RecHits = " << doRecHits_
         << " useDets " << ifEE_ << ":" << ifHE_ << std::endl;
 #endif
 
   edm::InputTag tmp0 = iConfig.getParameter<edm::InputTag>("GeneratorSrc");
   tok_hepMC_   = consumes<edm::HepMCProduct>(tmp0);
 #ifdef DebugLog
   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<HGCEEDigiCollection>(tmp2);
   edm::InputTag tmp3 = iConfig.getParameter<edm::InputTag>("RecHitSrcEE");
   tok_hitrEE_  = consumes<HGCRecHitCollection>(tmp3);
 #ifdef DebugLog
   if (ifEE_) {
     std::cout << "HGCalTBAnalyzer:: Detector " << detectorEE_ << " with tags "
           << tmp1 << ", " << tmp2 << ", " << tmp3 << std::endl;
   }
 #endif
   tmp1         = iConfig.getParameter<std::string>("CaloHitSrcHE");
   tok_hitsHE_  = consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits",tmp1));
   tmp2         = iConfig.getParameter<edm::InputTag>("DigiSrcHE");
   tok_digiHE_  = consumes<HGCHEDigiCollection>(tmp2);
   tmp3         = iConfig.getParameter<edm::InputTag>("RecHitSrcHE");
   tok_hitrHE_  = consumes<HGCRecHitCollection>(tmp3);
 #ifdef DebugLog
   if (ifHE_) {
     std::cout << "HGCalTBAnalyzer:: Detector " << detectorHE_ << " with tags "
           << tmp1 << ", " << tmp2 << ", " << tmp3 << std::endl;
   }
 #endif
 }

HGCalTBAnalyzer::~HGCalTBAnalyzer ( )

Definition at line 154 of file HGCalTBAnalyzer.cc.

154 {}

Member Function Documentation

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 313 of file HGCalTBAnalyzer.cc.

References ecalMGPA::adc(), analyzeDigi(), analyzeRecHits(), analyzeSimHits(), analyzeSimTracks(), edm::HandleBase::clear(), gather_cfg::cout, HGCSample::data(), detectorEE_, detectorHE_, doDigis_, doRecHits_, doSimHits_, doTree_, edm::Event::getByToken(), hBeam_, hgcons_, ifEE_, ifHE_, edm::HandleBase::isValid(), relval_2017::k, HGCalDDDConstants::layers(), AlCaHLTBitMon_ParallelJobs::p, sampleIndex_, simHitCellEnE, simHitCellEnH, simHitCellIdE, simHitCellIdH, simHitLayEn1E, simHitLayEn1H, simHitLayEn2E, simHitLayEn2H, tok_digiEE_, tok_digiHE_, tok_hepMC_, tok_hitrEE_, tok_hitrHE_, tok_hitsEE_, tok_hitsHE_, tok_simTk_, tok_simVtx_, and tree_.

                                                {
 
   //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 DebugLog
       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);
 
     simHitLayEn1E.clear(); simHitLayEn2E.clear();
     simHitLayEn1H.clear(); simHitLayEn2H.clear();
     simHitCellIdE.clear(); simHitCellEnE.clear();
     simHitCellIdH.clear(); simHitCellEnH.clear();
     edm::Handle<edm::PCaloHitContainer> theCaloHitContainers;
     std::vector<PCaloHit>               caloHits;
     if (ifEE_) {
       for (unsigned int k=0; k<hgcons_[0]->layers(false); ++k) 
     simHitLayEn1E.push_back(0);
       for (unsigned int k=0; k<hgcons_[0]->layers(true);  ++k) 
     simHitLayEn2E.push_back(0);
       iEvent.getByToken(tok_hitsEE_, theCaloHitContainers);
       if (theCaloHitContainers.isValid()) {
 #ifdef DebugLog
     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);
       } else {
 #ifdef DebugLog
     std::cout << "PCaloHitContainer does not exist for " << detectorEE_ 
           << " !!!" << std::endl;
 #endif
       }
     }
     if (ifHE_) {
       for (unsigned int k=0; k<hgcons_[1]->layers(false); ++k) 
     simHitLayEn1H.push_back(0);
       for (unsigned int k=0; k<hgcons_[1]->layers(true);  ++k) 
     simHitLayEn2H.push_back(0);
       iEvent.getByToken(tok_hitsHE_, theCaloHitContainers);
       if (theCaloHitContainers.isValid()) {
 #ifdef DebugLog
     std::cout << "PcalohitContainer for " << detectorHE_ << " has "
           << theCaloHitContainers->size() << " hits" << std::endl;
 #endif
     caloHits.clear();
     caloHits.insert(caloHits.end(), theCaloHitContainers->begin(), 
             theCaloHitContainers->end());
     analyzeSimHits(1, caloHits);
       } else {
 #ifdef DebugLog
     std::cout << "PCaloHitContainer does not exist for " << detectorHE_ 
           << " !!!" << std::endl;
 #endif
       }
     }
     if (doTree_) tree_->Fill();
   }
 
   //Now the Digis
   if (doDigis_) {
     if (ifEE_) {
       edm::Handle<HGCEEDigiCollection> theDigiContainers;
       iEvent.getByToken(tok_digiEE_, theDigiContainers);
       if (theDigiContainers.isValid()) {
 #ifdef DebugLog
     std::cout << "HGCDigiCintainer for " << detectorEE_ << " with " 
           << theDigiContainers->size() << " element(s)" << std::endl;
 #endif
     for (HGCEEDigiCollection::const_iterator it =theDigiContainers->begin();
          it !=theDigiContainers->end(); ++it) {
       HGCEEDetId detId     = (it->id());
       HGCSample  hgcSample = it->sample(sampleIndex_);
       uint16_t   adc       = hgcSample.data();
       analyzeDigi(0, detId, adc);
     }
       }
     }
     if (ifHE_) {
       edm::Handle<HGCHEDigiCollection> theDigiContainers;
       iEvent.getByToken(tok_digiHE_, theDigiContainers);
       if (theDigiContainers.isValid()) {
 #ifdef DebugLog
     std::cout << "HGCDigiContainer for " << detectorHE_ << " with " 
           << theDigiContainers->size() << " element(s)" << std::endl;
 #endif
     for (HGCHEDigiCollection::const_iterator it =theDigiContainers->begin();
          it !=theDigiContainers->end(); ++it) {
       HGCHEDetId detId     = (it->id());
       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 DebugLog
     std::cout << "HGCRecHitCollection for " << detectorEE_ << " has "
           << theCaloHitContainers->size() << " hits" << std::endl;
 #endif
     analyzeRecHits(0, theCaloHitContainers);
       } else {
 #ifdef DebugLog
     std::cout << "HGCRecHitCollection does not exist for " << detectorEE_ 
           << " !!!" << std::endl;
 #endif
       }
     }
     if (ifHE_) {
       iEvent.getByToken(tok_hitrHE_, theCaloHitContainers);
       if (theCaloHitContainers.isValid()) {
 #ifdef DebugLog
     std::cout << "HGCRecHitCollection for " << detectorHE_ << " has "
           << theCaloHitContainers->size() << " hits" << std::endl;
 #endif
     analyzeRecHits(1, theCaloHitContainers);
       } else {
 #ifdef DebugLog
     std::cout << "HGCRecHitCollection does not exist for " << detectorHE_ 
           << " !!!" << std::endl;
 #endif
       }
     }
   }
 }

template<class T1 >

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

private

Definition at line 634 of file HGCalTBAnalyzer.cc.

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

Referenced by analyze().

                                                                           {
 
   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::analyzeRecHits	(	int	type,
		edm::Handle< HGCRecHitCollection > &	hits
	)

private

Definition at line 643 of file HGCalTBAnalyzer.cc.

References HGCalDetId::cell(), gather_cfg::cout, relval_parameters_module::energy, 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().

                                                             {
  
   std::map<int,double>                   map_hitLayer;
   std::map<int,std::pair<DetId,double> > map_hitCell;
   for (HGCRecHitCollection::const_iterator it = hits->begin(); 
        it != hits->end(); ++it) {
     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 DebugLog
     std::cout << "RecHit: " << layer  << " " << global.x() << " " << global.y()
           << " " << global.z() << " " << energy << std::endl;
 #endif
   }
 
   for (std::map<int,double>::iterator itr = map_hitLayer.begin(); 
        itr != map_hitLayer.end(); ++itr)   {
     int    layer      = itr->first;
     double energy     = itr->second;
     double zp         = hgcons_[type]->waferZ(layer,true);
 #ifdef DebugLog
     std::cout << "SimHit:Layer " << layer << " " << zp << " " << energy 
           << std::endl;
 #endif
     hRecHitLng_[type]->Fill(zp,energy);
     hRecHitLng1_[type]->Fill(layer,energy);
   }
 
   for (std::map<int,std::pair<DetId,double> >::iterator itr = map_hitCell.begin(); 
        itr != map_hitCell.end(); ++itr) {
     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
	)

private

Definition at line 470 of file HGCalTBAnalyzer.cc.

References gather_cfg::cout, HLT_25ns10e33_v2_cff::depth, relval_parameters_module::energy, plotBeamSpotDB::first, hgcons_, hSimHitE_, hSimHitEn_, hSimHitLat_, hSimHitLayEn1E_, hSimHitLayEn1H_, hSimHitLayEn2E_, hSimHitLayEn2H_, hSimHitLng1_, hSimHitLng2_, hSimHitLng_, hSimHitT_, i, HGCalDDDConstants::locateCell(), HGCalTestNumbering::packHexagonIndex(), edm::second(), simHitCellEnE, simHitCellEnH, simHitCellIdE, simHitCellIdH, simHitLayEn1E, simHitLayEn1H, simHitLayEn2E, simHitLayEn2H, HGCalDDDConstants::simToReco(), HGCalTestNumbering::unpackHexagonIndex(), HGCalDDDConstants::waferZ(), create_public_lumi_plots::xy, and ecaldqm::zside().

Referenced by analyze().

                                                                          {
 
   std::map<uint32_t,double>                 map_hits, map_hitn;
   std::map<int,double>                      map_hitLayer, map_hitDepth;
   std::map<int,std::pair<uint32_t,double> > 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, cell;
     HGCalTestNumbering::unpackHexagonIndex(id, subdet, zside, layer, sector,
                        subsector, cell);
     int    depth       = hgcons_[type]->simToReco(cell,layer,sector,true).second;
 #ifdef DebugLog
     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) {
       map_hitLayer[layer] += energy;
     } else {
       map_hitLayer[layer]  = energy;
     }
     if (depth >= 0) {
       if (map_hitCell.count(cell) != 0) {
     double ee         = energy + map_hitCell[cell].second;
     map_hitCell[cell] = std::pair<uint32_t,double>(id,ee);
       } else {
     map_hitCell[cell] = 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  = 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 (std::map<uint32_t,double>::iterator itr = map_hits.begin() ; 
        itr != map_hits.end(); ++itr)   {
     hSimHitE_[type]->Fill(itr->second);
   }
 
   for (std::map<int,double>::iterator itr = map_hitLayer.begin(); 
        itr != map_hitLayer.end(); ++itr)   {
     int    layer      = itr->first - 1;
     double energy     = itr->second;
     double zp         = hgcons_[type]->waferZ(layer+1,false);
 #ifdef DebugLog
     std::cout << "SimHit:Layer " << layer+1 << " " << zp << " " << energy 
           << std::endl;
 #endif
     hSimHitLng_[type]->Fill(zp,energy);
     hSimHitLng2_[type]->Fill(layer+1,energy);
     if (type == 0) {
       if (layer < (int)(hSimHitLayEn1E_.size())) {
     simHitLayEn1E[layer] = energy;
     hSimHitLayEn1E_[layer]->Fill(energy);
       }
     } else {
       if (layer < (int)(hSimHitLayEn1H_.size())) {
     simHitLayEn1H[layer] = energy;
     hSimHitLayEn1H_[layer]->Fill(energy);
       }
     }
   }
   for (std::map<int,double>::iterator itr = map_hitDepth.begin(); 
        itr != map_hitDepth.end(); ++itr)   {
     int    layer      = itr->first - 1;
     double energy     = itr->second;
 #ifdef DebugLog
     std::cout << "SimHit:Layer " << layer+1 << " " << energy << std::endl;
 #endif
     hSimHitLng1_[type]->Fill(layer+1,energy);
     if (type == 0) {
       if (layer < (int)(hSimHitLayEn2E_.size())) {
     simHitLayEn2E[layer] = energy;
     hSimHitLayEn2E_[layer]->Fill(energy);
       }
     } else {
       if (layer < (int)(hSimHitLayEn2H_.size())) {
     simHitLayEn2H[layer] = energy;
     hSimHitLayEn2H_[layer]->Fill(energy);
       }
     }
   }
 
   for (std::map<int,std::pair<uint32_t,double> >::iterator itr = map_hitCell.begin(); 
        itr != map_hitCell.end(); ++itr) {
     uint32_t id       = ((itr->second).first);
     double   energy   = ((itr->second).second);
     int      subdet, zside, layer, sector, subsector, cell;
     HGCalTestNumbering::unpackHexagonIndex(id, subdet, zside, layer, sector,
                        subsector, cell);
     std::pair<float,float> xy = hgcons_[type]->locateCell(cell,layer,sector,false);
     double zp         = hgcons_[type]->waferZ(layer,false);
     double xx         = (zp < 0) ? -xy.first : xy.first;
     hSimHitLat_[type]->Fill(xx,xy.second,energy);
   }
 
   for (std::map<uint32_t,double>::iterator itr = map_hitn.begin(); 
        itr != map_hitn.end(); ++itr) {
     uint32_t id     = itr->first;
     double   energy = itr->second;
     if (type == 0) {
       simHitCellIdE.push_back(id); simHitCellEnE.push_back(energy);
     } else {
       simHitCellIdH.push_back(id); simHitCellEnH.push_back(energy);
     }
   }
 }

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

private

Definition at line 600 of file HGCalTBAnalyzer.cc.

References gather_cfg::cout, pBeam, xBeam, yBeam, and zBeam.

Referenced by analyze().

                                                                        {
 
   xBeam = yBeam = zBeam = pBeam = -1000000;
   int vertIndex(-1);
   for (edm::SimTrackContainer::const_iterator simTrkItr = SimTk->begin(); 
        simTrkItr!= SimTk->end(); simTrkItr++) {
 #ifdef DebugLog
     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 DebugLog
     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 164 of file HGCalTBAnalyzer.cc.

References detectorEE_, detectorHE_, doDigis_, doRecHits_, doSimHits_, doTree_, doTreeCell_, fs_, hBeam_, hDigiADC_, hDigiLng_, hDigiOcc_, hRecHitE_, hRecHitLat_, hRecHitLng1_, hRecHitLng_, hRecHitOcc_, hSimHitE_, hSimHitEn_, hSimHitLat_, hSimHitLng1_, hSimHitLng2_, hSimHitLng_, hSimHitT_, i, ifEE_, ifHE_, TFileService::make(), mergeVDriftHistosByStation::name, pBeam, simHitCellEnE, simHitCellEnH, simHitCellIdE, simHitCellIdH, simHitLayEn1E, simHitLayEn1H, simHitLayEn2E, simHitLayEn2H, AlCaHLTBitMon_QueryRunRegistry::string, SiPixelPhase1TrackClusters_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<2; ++i) {
     bool book = (i == 0) ? ifEE_ : ifHE_;
     std::string det = (i == 0) ? detectorEE_ : detectorHE_;
 
     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 (doSimHits_ && doTree_) {
     tree_ = fs_->make<TTree>("HGCTB","SimHitEnergy");
     tree_->Branch("simHitLayEn1E", &simHitLayEn1E);
     tree_->Branch("simHitLayEn2E", &simHitLayEn2E);
     tree_->Branch("simHitLayEn1H", &simHitLayEn1H);
     tree_->Branch("simHitLayEn2H", &simHitLayEn2H);
     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("simHitCellIdE", &simHitCellIdE);
       tree_->Branch("simHitCellEnE", &simHitCellEnE);
       tree_->Branch("simHitCellIdH", &simHitCellIdH);
       tree_->Branch("simHitCellEnH", &simHitCellEnH);
     }
   }
 }

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

overrideprivatevirtual

Definition at line 245 of file HGCalTBAnalyzer.cc.

References gather_cfg::cout, detectorEE_, detectorHE_, doDigis_, doRecHits_, fs_, relativeConstraints::geom, edm::EventSetup::get(), hgcons_, hgeom_, hSimHitLayEn1E_, hSimHitLayEn1H_, hSimHitLayEn2E_, ifEE_, ifHE_, cmsLHEtoEOSManager::l, HGCalDDDConstants::layers(), TFileService::make(), mergeVDriftHistosByStation::name, edm::ESHandle< class >::product(), and SiPixelPhase1TrackClusters_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] = 0;
     }
     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());
       hSimHitLayEn1E_.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());
     hSimHitLayEn2E_.push_back(fs_->make<TH1D>(name,title,100000,0.,0.2));
       }
     }
 #ifdef DebugLog
     std::cout << "HGCalTBAnalyzer::" << detectorEE_ << " defined with "
           << hgcons_[0]->layers(false) << " layers" << std::endl;
 #endif
   } else {
     hgcons_[0] = 0;
     hgeom_[0]  = 0;
   }
 
   if (ifHE_) {
     edm::ESHandle<HGCalDDDConstants>  pHGDC;
     iSetup.get<IdealGeometryRecord>().get(detectorHE_, pHGDC);
     hgcons_[1] = &(*pHGDC);
     if (doDigis_ || doRecHits_) {
       edm::ESHandle<HGCalGeometry> geom;
       iSetup.get<IdealGeometryRecord>().get(detectorHE_, geom);
       hgeom_[1] = geom.product();
     } else {
       hgeom_[1] = 0;
     }
     for (unsigned int l=0; l<hgcons_[1]->layers(false); ++l) {
       sprintf (name, "SimHitEnA%d%s", l, detectorHE_.c_str());
       sprintf (title,"Sim Hit Energy in layer %d for %s",l+1,
            detectorHE_.c_str());
       hSimHitLayEn1H_.push_back(fs_->make<TH1D>(name,title,100000,0.,0.2));
       if (l%3 == 0) {
     sprintf (name, "SimHitEnB%d%s", (l/3+1), detectorHE_.c_str());
     sprintf (title,"Sim Hit Energy in layer %d for %s",(l/3+1),
          detectorHE_.c_str());
     hSimHitLayEn2E_.push_back(fs_->make<TH1D>(name,title,100000,0.,0.2));
       }
     }
 #ifdef DebugLog
     std::cout << "HGCalTBAnalyzer::" << detectorHE_ << " defined with "
           << hgcons_[1]->layers(false) << " layers" << std::endl;
 #endif
   } else {
     hgcons_[1] = 0;
     hgeom_[1]  = 0;
   }
 }

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

inlineoverrideprivatevirtual

Definition at line 60 of file HGCalTBAnalyzer.cc.

60 {}

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

static

Definition at line 156 of file HGCalTBAnalyzer.cc.

References edm::ConfigurationDescriptions::addDefault(), and edm::ParameterSetDescription::setUnknown().

                                                                                  {
   //The following says we do not know what parameters are allowed so do no validation
   // Please change this to state exactly what you do use, even if it is no parameters
   edm::ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }

Member Data Documentation

std::string HGCalTBAnalyzer::detectorEE_

private

Definition at line 74 of file HGCalTBAnalyzer.cc.

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

std::string HGCalTBAnalyzer::detectorHE_

private

Definition at line 74 of file HGCalTBAnalyzer.cc.

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

bool HGCalTBAnalyzer::doDigis_

private

Definition at line 73 of file HGCalTBAnalyzer.cc.

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

bool HGCalTBAnalyzer::doRecHits_

private

Definition at line 73 of file HGCalTBAnalyzer.cc.

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

bool HGCalTBAnalyzer::doSimHits_

private

Definition at line 73 of file HGCalTBAnalyzer.cc.

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

bool HGCalTBAnalyzer::doTree_

private

Definition at line 72 of file HGCalTBAnalyzer.cc.

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

bool HGCalTBAnalyzer::doTreeCell_

private

Definition at line 72 of file HGCalTBAnalyzer.cc.

Referenced by beginJob(), and HGCalTBAnalyzer().

edm::Service<TFileService> HGCalTBAnalyzer::fs_

private

Definition at line 68 of file HGCalTBAnalyzer.cc.

Referenced by beginJob(), and beginRun().

TH1D * HGCalTBAnalyzer::hBeam_

private

Definition at line 85 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and beginJob().

TH1D* HGCalTBAnalyzer::hDigiADC_[2]

private

Definition at line 84 of file HGCalTBAnalyzer.cc.

Referenced by analyzeDigi(), and beginJob().

TH1D * HGCalTBAnalyzer::hDigiLng_[2]

private

Definition at line 84 of file HGCalTBAnalyzer.cc.

Referenced by analyzeDigi(), and beginJob().

TH2D* HGCalTBAnalyzer::hDigiOcc_[2]

private

Definition at line 86 of file HGCalTBAnalyzer.cc.

Referenced by analyzeDigi(), and beginJob().

const HGCalDDDConstants* HGCalTBAnalyzer::hgcons_[2]

private

Definition at line 69 of file HGCalTBAnalyzer.cc.

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

const HGCalGeometry* HGCalTBAnalyzer::hgeom_[2]

private

Definition at line 70 of file HGCalTBAnalyzer.cc.

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

TH1D* HGCalTBAnalyzer::hRecHitE_[2]

private

Definition at line 85 of file HGCalTBAnalyzer.cc.

Referenced by analyzeRecHits(), and beginJob().

TProfile2D * HGCalTBAnalyzer::hRecHitLat_[2]

private

Definition at line 90 of file HGCalTBAnalyzer.cc.

Referenced by analyzeRecHits(), and beginJob().

TProfile * HGCalTBAnalyzer::hRecHitLng1_[2]

private

Definition at line 89 of file HGCalTBAnalyzer.cc.

Referenced by analyzeRecHits(), and beginJob().

TProfile* HGCalTBAnalyzer::hRecHitLng_[2]

private

Definition at line 89 of file HGCalTBAnalyzer.cc.

Referenced by analyzeRecHits(), and beginJob().

TH2D * HGCalTBAnalyzer::hRecHitOcc_[2]

private

Definition at line 86 of file HGCalTBAnalyzer.cc.

Referenced by analyzeRecHits(), and beginJob().

TH1D* HGCalTBAnalyzer::hSimHitE_[2]

private

Definition at line 83 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

TH1D * HGCalTBAnalyzer::hSimHitEn_[2]

private

Definition at line 85 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

TProfile2D* HGCalTBAnalyzer::hSimHitLat_[2]

private

Definition at line 90 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

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

private

Definition at line 91 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

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

private

Definition at line 91 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

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

private

Definition at line 92 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

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

private

Definition at line 92 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits().

TProfile * HGCalTBAnalyzer::hSimHitLng1_[2]

private

Definition at line 87 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

TProfile* HGCalTBAnalyzer::hSimHitLng2_[2]

private

Definition at line 88 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

TProfile* HGCalTBAnalyzer::hSimHitLng_[2]

private

Definition at line 87 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

TH1D * HGCalTBAnalyzer::hSimHitT_[2]

private

Definition at line 83 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

bool HGCalTBAnalyzer::ifEE_

private

Definition at line 71 of file HGCalTBAnalyzer.cc.

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

bool HGCalTBAnalyzer::ifHE_

private

Definition at line 71 of file HGCalTBAnalyzer.cc.

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

double HGCalTBAnalyzer::pBeam

private

Definition at line 97 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

int HGCalTBAnalyzer::sampleIndex_

private

Definition at line 75 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

std::vector<float> HGCalTBAnalyzer::simHitCellEnE

private

Definition at line 96 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::simHitCellEnH

private

Definition at line 96 of file HGCalTBAnalyzer.cc.

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

std::vector<uint32_t> HGCalTBAnalyzer::simHitCellIdE

private

Definition at line 95 of file HGCalTBAnalyzer.cc.

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

std::vector<uint32_t> HGCalTBAnalyzer::simHitCellIdH

private

Definition at line 95 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::simHitLayEn1E

private

Definition at line 93 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::simHitLayEn1H

private

Definition at line 94 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::simHitLayEn2E

private

Definition at line 93 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::simHitLayEn2H

private

Definition at line 94 of file HGCalTBAnalyzer.cc.

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

edm::EDGetToken HGCalTBAnalyzer::tok_digiEE_

private

Definition at line 79 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

edm::EDGetToken HGCalTBAnalyzer::tok_digiHE_

private

Definition at line 79 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

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

private

Definition at line 81 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

edm::EDGetToken HGCalTBAnalyzer::tok_hitrEE_

private

Definition at line 80 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

edm::EDGetToken HGCalTBAnalyzer::tok_hitrHE_

private

Definition at line 80 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

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

private

Definition at line 76 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

edm::EDGetTokenT<edm::PCaloHitContainer> HGCalTBAnalyzer::tok_hitsHE_

private

Definition at line 76 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

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

private

Definition at line 77 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

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

private

Definition at line 78 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

TTree* HGCalTBAnalyzer::tree_

private

Definition at line 82 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HGCalTBAnalyzer::xBeam

private

Definition at line 97 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

double HGCalTBAnalyzer::yBeam

private

Definition at line 97 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

double HGCalTBAnalyzer::zBeam

private

Definition at line 97 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation