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

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

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

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

Constructor & Destructor Documentation

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

explicit

Passive hits

Definition at line 127 of file HGCalTBAnalyzer.cc.

References gather_cfg::cout, detectorBeam_, detectorBH_, detectorEE_, detectorFH_, doDigis_, doPassive_, doRecHits_, doSimHits_, doTree_, doTreeCell_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hcalTTPDigis_cfi::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_.size() == 0) 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<HGCEEDigiCollection>(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<HGCHEDigiCollection>(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<HGCBHDigiCollection>(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 ( )

Definition at line 231 of file HGCalTBAnalyzer.cc.

231 {}

Member Function Documentation

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

overrideprivatevirtual

EE

FH

BH

CMSE

Implements edm::one::EDAnalyzerBase.

Definition at line 486 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<HGCEEDigiCollection> 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) {
       HGCEEDetId detId     = (it.id());
       HGCSample  hgcSample = it.sample(sampleIndex_);
       uint16_t   adc       = hgcSample.data();
       analyzeDigi(0, detId, adc);
     }
       }
     }
     if (ifFH_) {
       edm::Handle<HGCHEDigiCollection> 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) {
       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 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 925 of file HGCalTBAnalyzer.cc.

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

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 984 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 934 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 713 of file HGCalTBAnalyzer.cc.

References gather_cfg::cout, particleFlowClusterECALTimeSelected_cfi::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(), hcalTTPDigis_cfi::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, HcalDigiParam_cfi::zside, and AHCalDetId::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 892 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 276 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 &
	)

overrideprivatevirtual

Definition at line 396 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] = 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());
       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] = 0;
     hgeom_[0]  = 0;
   }
 
   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] = 0;
     }
     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] = 0;
     hgeom_[1]  = 0;
   }
 
   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));
     }
   }
 }

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

inlineoverrideprivatevirtual

Definition at line 66 of file HGCalTBAnalyzer.cc.

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

66 {}

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

static

Definition at line 233 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("mix","HGCDigisEE"));
   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("mix","HGCDigisHEfront"));
   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("mix","HGCDigisHEback"));
   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 84 of file HGCalTBAnalyzer.cc.

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

std::string HGCalTBAnalyzer::detectorBH_

private

Definition at line 84 of file HGCalTBAnalyzer.cc.

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

std::string HGCalTBAnalyzer::detectorEE_

private

Definition at line 83 of file HGCalTBAnalyzer.cc.

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

std::string HGCalTBAnalyzer::detectorFH_

private

Definition at line 83 of file HGCalTBAnalyzer.cc.

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

bool HGCalTBAnalyzer::doDigis_

private

Definition at line 81 of file HGCalTBAnalyzer.cc.

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

bool HGCalTBAnalyzer::doPassive_

private

Definition at line 82 of file HGCalTBAnalyzer.cc.

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

bool HGCalTBAnalyzer::doRecHits_

private

Definition at line 81 of file HGCalTBAnalyzer.cc.

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

bool HGCalTBAnalyzer::doSimHits_

private

Definition at line 81 of file HGCalTBAnalyzer.cc.

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

bool HGCalTBAnalyzer::doTree_

private

Definition at line 80 of file HGCalTBAnalyzer.cc.

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

bool HGCalTBAnalyzer::doTreeCell_

private

Definition at line 80 of file HGCalTBAnalyzer.cc.

Referenced by beginJob(), and HGCalTBAnalyzer().

edm::Service<TFileService> HGCalTBAnalyzer::fs_

private

Definition at line 76 of file HGCalTBAnalyzer.cc.

Referenced by beginJob(), and beginRun().

TH1D * HGCalTBAnalyzer::hBeam_

private

Definition at line 101 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and beginJob().

TH1D* HGCalTBAnalyzer::hDigiADC_[3]

private

Definition at line 100 of file HGCalTBAnalyzer.cc.

Referenced by analyzeDigi(), and beginJob().

TH1D * HGCalTBAnalyzer::hDigiLng_[2]

private

Definition at line 100 of file HGCalTBAnalyzer.cc.

Referenced by analyzeDigi(), and beginJob().

TH2D* HGCalTBAnalyzer::hDigiOcc_[3]

private

Definition at line 102 of file HGCalTBAnalyzer.cc.

Referenced by analyzeDigi(), and beginJob().

const HGCalDDDConstants* HGCalTBAnalyzer::hgcons_[2]

private

Definition at line 77 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::hgcPassiveBHEnergy_

private

Definition at line 120 of file HGCalTBAnalyzer.cc.

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

std::vector<int> HGCalTBAnalyzer::hgcPassiveBHID_

private

Definition at line 122 of file HGCalTBAnalyzer.cc.

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

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

private

Definition at line 121 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::hgcPassiveCMSEEnergy_

private

Definition at line 120 of file HGCalTBAnalyzer.cc.

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

std::vector<int> HGCalTBAnalyzer::hgcPassiveCMSEID_

private

Definition at line 122 of file HGCalTBAnalyzer.cc.

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

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

private

Definition at line 121 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::hgcPassiveEEEnergy_

private

Definition at line 120 of file HGCalTBAnalyzer.cc.

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

std::vector<int> HGCalTBAnalyzer::hgcPassiveEEID_

private

Definition at line 122 of file HGCalTBAnalyzer.cc.

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

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

private

Definition at line 121 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::hgcPassiveFHEnergy_

private

Definition at line 120 of file HGCalTBAnalyzer.cc.

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

std::vector<int> HGCalTBAnalyzer::hgcPassiveFHID_

private

Definition at line 122 of file HGCalTBAnalyzer.cc.

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

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

private

Definition at line 121 of file HGCalTBAnalyzer.cc.

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

const HGCalGeometry* HGCalTBAnalyzer::hgeom_[2]

private

Definition at line 78 of file HGCalTBAnalyzer.cc.

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

TH1D* HGCalTBAnalyzer::hRecHitE_[3]

private

Definition at line 101 of file HGCalTBAnalyzer.cc.

Referenced by analyzeRecHits(), and beginJob().

TProfile2D * HGCalTBAnalyzer::hRecHitLat_[3]

private

Definition at line 106 of file HGCalTBAnalyzer.cc.

Referenced by analyzeRecHits(), and beginJob().

TProfile * HGCalTBAnalyzer::hRecHitLng1_[3]

private

Definition at line 105 of file HGCalTBAnalyzer.cc.

Referenced by analyzeRecHits(), and beginJob().

TProfile* HGCalTBAnalyzer::hRecHitLng_[3]

private

Definition at line 105 of file HGCalTBAnalyzer.cc.

Referenced by analyzeRecHits(), and beginJob().

TH2D * HGCalTBAnalyzer::hRecHitOcc_[3]

private

Definition at line 102 of file HGCalTBAnalyzer.cc.

Referenced by analyzeRecHits(), and beginJob().

TH1D* HGCalTBAnalyzer::hSimHitE_[4]

private

Definition at line 99 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

TH1D * HGCalTBAnalyzer::hSimHitEn_[4]

private

Definition at line 101 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

TProfile2D* HGCalTBAnalyzer::hSimHitLat_[3]

private

Definition at line 106 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

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

private

Definition at line 109 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

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

private

Definition at line 107 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

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

private

Definition at line 108 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

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

private

Definition at line 109 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

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

private

Definition at line 107 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

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

private

Definition at line 108 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

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

private

Definition at line 110 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

TProfile * HGCalTBAnalyzer::hSimHitLng1_[3]

private

Definition at line 103 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

TProfile* HGCalTBAnalyzer::hSimHitLng2_[3]

private

Definition at line 104 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

TProfile* HGCalTBAnalyzer::hSimHitLng_[3]

private

Definition at line 103 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

TH1D * HGCalTBAnalyzer::hSimHitT_[4]

private

Definition at line 99 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

std::vector<int> HGCalTBAnalyzer::idBeams_

private

Definition at line 87 of file HGCalTBAnalyzer.cc.

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

bool HGCalTBAnalyzer::ifBeam_

private

Definition at line 79 of file HGCalTBAnalyzer.cc.

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

bool HGCalTBAnalyzer::ifBH_

private

Definition at line 79 of file HGCalTBAnalyzer.cc.

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

bool HGCalTBAnalyzer::ifEE_

private

Definition at line 79 of file HGCalTBAnalyzer.cc.

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

bool HGCalTBAnalyzer::ifFH_

private

Definition at line 79 of file HGCalTBAnalyzer.cc.

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

double HGCalTBAnalyzer::pBeam_

private

Definition at line 124 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

int HGCalTBAnalyzer::sampleIndex_

private

Definition at line 86 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

std::vector<float> HGCalTBAnalyzer::simHitCellEnBeam_

private

Definition at line 118 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::simHitCellEnBH_

private

Definition at line 118 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::simHitCellEnEE_

private

Definition at line 117 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::simHitCellEnFH_

private

Definition at line 117 of file HGCalTBAnalyzer.cc.

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

std::vector<uint32_t> HGCalTBAnalyzer::simHitCellIdBeam_

private

Definition at line 116 of file HGCalTBAnalyzer.cc.

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

std::vector<uint32_t> HGCalTBAnalyzer::simHitCellIdBH_

private

Definition at line 116 of file HGCalTBAnalyzer.cc.

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

std::vector<uint32_t> HGCalTBAnalyzer::simHitCellIdEE_

private

Definition at line 115 of file HGCalTBAnalyzer.cc.

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

std::vector<uint32_t> HGCalTBAnalyzer::simHitCellIdFH_

private

Definition at line 115 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::simHitLayEn1BH_

private

Definition at line 113 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::simHitLayEn1EE_

private

Definition at line 111 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::simHitLayEn1FH_

private

Definition at line 112 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::simHitLayEn2BH_

private

Definition at line 113 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::simHitLayEn2EE_

private

Definition at line 111 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::simHitLayEn2FH_

private

Definition at line 112 of file HGCalTBAnalyzer.cc.

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

std::vector<float> HGCalTBAnalyzer::simHitLayEnBeam_

private

Definition at line 114 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and analyzeSimHits().

edm::EDGetToken HGCalTBAnalyzer::tok_digiBH_

private

Definition at line 92 of file HGCalTBAnalyzer.cc.

Referenced by HGCalTBAnalyzer().

edm::EDGetToken HGCalTBAnalyzer::tok_digiEE_

private

Definition at line 92 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

edm::EDGetToken HGCalTBAnalyzer::tok_digiFH_

private

Definition at line 92 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

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

private

Definition at line 94 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

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

private

Definition at line 96 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

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

private

Definition at line 96 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

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

private

Definition at line 95 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

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

private

Definition at line 95 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

edm::EDGetToken HGCalTBAnalyzer::tok_hitrBH_

private

Definition at line 93 of file HGCalTBAnalyzer.cc.

Referenced by HGCalTBAnalyzer().

edm::EDGetToken HGCalTBAnalyzer::tok_hitrEE_

private

Definition at line 93 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

edm::EDGetToken HGCalTBAnalyzer::tok_hitrFH_

private

Definition at line 93 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

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

private

Definition at line 89 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

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

private

Definition at line 89 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

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

private

Definition at line 88 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

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

private

Definition at line 88 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

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

private

Definition at line 90 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

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

private

Definition at line 91 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

TTree* HGCalTBAnalyzer::tree_

private

Definition at line 98 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HGCalTBAnalyzer::xBeam_

private

Definition at line 124 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

double HGCalTBAnalyzer::yBeam_

private

Definition at line 124 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

double HGCalTBAnalyzer::zBeam_

private

Definition at line 124 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

double HGCalTBAnalyzer::zFrontBH_

private

Definition at line 85 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

double HGCalTBAnalyzer::zFrontEE_

private

Definition at line 85 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

double HGCalTBAnalyzer::zFrontFH_

private

Definition at line 85 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

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