Inheritance diagram for HGCalTBAnalyzer:

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

	~HGCalTBAnalyzer () override

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

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () override

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

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
void	analyze (edm::Event const &, edm::EventSetup const &) override

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

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

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

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

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

void	beginJob () override

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

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

Static Private Member Functions
static bool	sortTime (const std::pair< double, double > &i, const std::pair< double, double > &j)

Private Attributes
const bool	addP_

std::unique_ptr< AHCalGeometry >	ahcalGeom_

const std::string	detectorBeam_

const std::string	detectorBH_

const std::string	detectorEE_

const std::string	detectorFH_

const bool	doBeam_

const bool	doDigis_

const bool	doPassive_

const bool	doPassiveBH_

const bool	doPassiveEE_

const bool	doPassiveHE_

const bool	doRecHits_

const bool	doSimHits_

const bool	doTree_

const bool	doTreeCell_

edm::Service< TFileService >	fs_

const double	gev2mip200_

const double	gev2mip300_

TH1D *	hBeam_

TH1D *	hDigiADC_ [3]

TH1D *	hDigiLng_ [2]

TH2D *	hDigiOcc_ [3]

const HGCalDDDConstants *	hgcons_ [2]

std::vector< float >	hgcPassiveBeamEnergy_

std::vector< int >	hgcPassiveBeamID_

std::vector< std::string >	hgcPassiveBeamName_

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_

const bool	ifBeam_

const bool	ifBH_

const bool	ifEE_

const bool	ifFH_

const edm::InputTag	labelDigiBH_

const edm::InputTag	labelDigiEE_

const edm::InputTag	labelDigiFH_

const edm::InputTag	labelGen_

const std::string	labelHitBeam_

const std::string	labelHitBH_

const std::string	labelHitEE_

const std::string	labelHitFH_

const edm::InputTag	labelPassiveBeam_

const edm::InputTag	labelPassiveBH_

const edm::InputTag	labelPassiveCMSE_

const edm::InputTag	labelPassiveEE_

const edm::InputTag	labelPassiveFH_

const edm::InputTag	labelRHitBH_

const edm::InputTag	labelRHitEE_

const edm::InputTag	labelRHitFH_

int	nBeamMC_

double	pBeam_

std::vector< float >	pBeamMC_

std::vector< int >	pdgIdBeamMC_

double	phiBeam_

std::vector< float >	pxBeamMC_

std::vector< float >	pyBeamMC_

std::vector< float >	pzBeamMC_

const int	sampleIndex_

std::vector< int >	simHitCellColBH_

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< int >	simHitCellLayerBH_

std::vector< int >	simHitCellRowBH_

std::vector< float >	simHitCellTime15MipBH_

std::vector< float >	simHitCellTime15MipEE_

std::vector< float >	simHitCellTime15MipFH_

std::vector< float >	simHitCellTimeFirstHitBH_

std::vector< float >	simHitCellTimeFirstHitEE_

std::vector< float >	simHitCellTimeFirstHitFH_

std::vector< float >	simHitCellTimeLastHitBH_

std::vector< float >	simHitCellTimeLastHitEE_

std::vector< float >	simHitCellTimeLastHitFH_

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_

const double	stoc_smear_time_200_

const double	stoc_smear_time_300_

double	thetaBeam_

const edm::EDGetTokenT< HGCalDigiCollection >	tok_digiBH_

const edm::EDGetTokenT< HGCalDigiCollection >	tok_digiEE_

const edm::EDGetTokenT< HGCalDigiCollection >	tok_digiFH_

const edm::EDGetTokenT< edm::HepMCProduct >	tok_hepMC_

const edm::EDGetTokenT< edm::PassiveHitContainer >	tok_hgcPHBeam_

const edm::EDGetTokenT< edm::PassiveHitContainer >	tok_hgcPHBH_

const edm::EDGetTokenT< edm::PassiveHitContainer >	tok_hgcPHCMSE_

const edm::EDGetTokenT< edm::PassiveHitContainer >	tok_hgcPHEE_

const edm::EDGetTokenT< edm::PassiveHitContainer >	tok_hgcPHFH_

const edm::EDGetTokenT< HGCRecHitCollection >	tok_hitrBH_

const edm::EDGetTokenT< HGCRecHitCollection >	tok_hitrEE_

const edm::EDGetTokenT< HGCRecHitCollection >	tok_hitrFH_

const edm::EDGetTokenT< edm::PCaloHitContainer >	tok_hitsBeam_

const edm::EDGetTokenT< edm::PCaloHitContainer >	tok_hitsBH_

const edm::EDGetTokenT< edm::PCaloHitContainer >	tok_hitsEE_

const edm::EDGetTokenT< edm::PCaloHitContainer >	tok_hitsFH_

const edm::EDGetTokenT< edm::SimTrackContainer >	tok_simTk_

const edm::EDGetTokenT< edm::SimVertexContainer >	tok_simVtx_

const edm::ESGetToken< HGCalDDDConstants, IdealGeometryRecord >	tokDDDEE_

const edm::ESGetToken< HGCalDDDConstants, IdealGeometryRecord >	tokDDDFH_

const edm::ESGetToken< HGCalGeometry, IdealGeometryRecord >	tokGeomEE_

const edm::ESGetToken< HGCalGeometry, IdealGeometryRecord >	tokGeomFH_

TTree *	tree_

double	xBeam_

std::vector< float >	xBeamMC_

double	yBeam_

std::vector< float >	yBeamMC_

double	zBeam_

std::vector< float >	zBeamMC_

const double	zFrontBH_

const double	zFrontEE_

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

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 54 of file HGCalTBAnalyzer.cc.

Constructor & Destructor Documentation

◆ HGCalTBAnalyzer()

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

explicit

TIME SMEARING Const = 50ps = 0.05ns Stochastic is 4ns/Q(fC)

Definition at line 150 of file HGCalTBAnalyzer.cc.

References addP_, ahcalGeom_, detectorBeam_, detectorBH_, detectorEE_, detectorFH_, doDigis_, doPassive_, doPassiveBH_, doPassiveEE_, doPassiveHE_, doRecHits_, doSimHits_, edm::ParameterSet::getParameter(), gev2mip200_, gev2mip300_, idBeams_, ifBeam_, ifBH_, ifEE_, ifFH_, dqmdumpme::k, labelDigiEE_, labelDigiFH_, labelGen_, labelHitBeam_, labelHitBH_, labelHitEE_, labelHitFH_, labelRHitBH_, labelRHitEE_, labelRHitFH_, sampleIndex_, stoc_smear_time_200_, stoc_smear_time_300_, zFrontBH_, zFrontEE_, and zFrontFH_.

     : ifEE_(iConfig.getParameter<bool>("useEE")),
       ifFH_(iConfig.getParameter<bool>("useFH")),
       ifBH_(iConfig.getParameter<bool>("useBH")),
       ifBeam_(iConfig.getParameter<bool>("useBeam")),
       doSimHits_(iConfig.getParameter<bool>("doSimHits")),
       doDigis_(iConfig.getParameter<bool>("doDigis")),
       doRecHits_(iConfig.getParameter<bool>("doRecHits")),
       doTree_(iConfig.getParameter<bool>("doTree")),
       doTreeCell_(iConfig.getParameter<bool>("doTreeCell")),
       doPassive_(iConfig.getParameter<bool>("doPassive")),
       doPassiveEE_(iConfig.getParameter<bool>("doPassiveEE")),
       doPassiveHE_(iConfig.getParameter<bool>("doPassiveHE")),
       doPassiveBH_(iConfig.getParameter<bool>("doPassiveBH")),
       addP_(iConfig.getParameter<bool>("addP")),
       doBeam_(iConfig.getParameter<bool>("doBeam")),
       detectorEE_(iConfig.getParameter<std::string>("detectorEE")),
       detectorFH_(iConfig.getParameter<std::string>("detectorFH")),
       detectorBH_(iConfig.getParameter<std::string>("detectorBH")),
       detectorBeam_(iConfig.getParameter<std::string>("detectorBeam")),
       zFrontEE_(iConfig.getParameter<double>("zFrontEE")),
       zFrontFH_(iConfig.getParameter<double>("zFrontFH")),
       zFrontBH_(iConfig.getParameter<double>("zFrontBH")),
       sampleIndex_(iConfig.getParameter<int>("sampleIndex")),
       gev2mip200_(iConfig.getUntrackedParameter<double>("gev2mip200", 57.0e-6)),
       gev2mip300_(iConfig.getUntrackedParameter<double>("gev2mip300", 85.5e-6)),
       stoc_smear_time_200_(iConfig.getUntrackedParameter<double>("stoc_smear_time_200", 10.24)),
       stoc_smear_time_300_(iConfig.getUntrackedParameter<double>("stoc_smear_time_300", 15.5)),
       labelGen_(iConfig.getParameter<edm::InputTag>("generatorSrc")),
       labelHitEE_(iConfig.getParameter<std::string>("caloHitSrcEE")),
       labelHitFH_(iConfig.getParameter<std::string>("caloHitSrcFH")),
       labelHitBH_(iConfig.getParameter<std::string>("caloHitSrcBH")),
       labelHitBeam_(iConfig.getParameter<std::string>("caloHitSrcBeam")),
       labelDigiEE_(iConfig.getParameter<edm::InputTag>("digiSrcEE")),
       labelDigiFH_(iConfig.getParameter<edm::InputTag>("digiSrcFH")),
       labelDigiBH_(iConfig.getParameter<edm::InputTag>("digiSrcBH")),
       labelRHitEE_(iConfig.getParameter<edm::InputTag>("recHitSrcEE")),
       labelRHitFH_(iConfig.getParameter<edm::InputTag>("recHitSrcFH")),
       labelRHitBH_(iConfig.getParameter<edm::InputTag>("recHitSrcBH")),
       labelPassiveEE_(iConfig.getParameter<edm::InputTag>("passiveEE")),
       labelPassiveFH_(iConfig.getParameter<edm::InputTag>("passiveFH")),
       labelPassiveBH_(iConfig.getParameter<edm::InputTag>("passiveBH")),
       labelPassiveCMSE_(iConfig.getParameter<edm::InputTag>("passiveCMSE")),
       labelPassiveBeam_(iConfig.getParameter<edm::InputTag>("passiveBeam")),
       tok_hepMC_(consumes<edm::HepMCProduct>(labelGen_)),
       tok_simTk_(consumes<edm::SimTrackContainer>(edm::InputTag("g4SimHits"))),
       tok_simVtx_(consumes<edm::SimVertexContainer>(edm::InputTag("g4SimHits"))),
       tok_hitsEE_(consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits", labelHitEE_))),
       tok_hitsFH_(consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits", labelHitFH_))),
       tok_hitsBH_(consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits", labelHitBH_))),
       tok_hitsBeam_(consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits", labelHitBeam_))),
       tok_digiEE_(consumes<HGCalDigiCollection>(labelDigiEE_)),
       tok_digiFH_(consumes<HGCalDigiCollection>(labelDigiFH_)),
       tok_digiBH_(consumes<HGCalDigiCollection>(labelDigiFH_)),
       tok_hitrEE_(consumes<HGCRecHitCollection>(labelRHitEE_)),
       tok_hitrFH_(consumes<HGCRecHitCollection>(labelRHitFH_)),
       tok_hitrBH_(consumes<HGCRecHitCollection>(labelRHitBH_)),
       tok_hgcPHEE_(consumes<edm::PassiveHitContainer>(labelPassiveEE_)),
       tok_hgcPHFH_(consumes<edm::PassiveHitContainer>(labelPassiveFH_)),
       tok_hgcPHBH_(consumes<edm::PassiveHitContainer>(labelPassiveBH_)),
       tok_hgcPHCMSE_(consumes<edm::PassiveHitContainer>(labelPassiveCMSE_)),
       tok_hgcPHBeam_(consumes<edm::PassiveHitContainer>(labelPassiveBeam_)),
       tokDDDEE_(esConsumes<HGCalDDDConstants, IdealGeometryRecord, edm::Transition::BeginRun>(
           edm::ESInputTag("", detectorEE_))),
       tokGeomEE_(
           esConsumes<HGCalGeometry, IdealGeometryRecord, edm::Transition::BeginRun>(edm::ESInputTag("", detectorEE_))),
       tokDDDFH_(esConsumes<HGCalDDDConstants, IdealGeometryRecord, edm::Transition::BeginRun>(
           edm::ESInputTag("", detectorFH_))),
       tokGeomFH_(
           esConsumes<HGCalGeometry, IdealGeometryRecord, edm::Transition::BeginRun>(edm::ESInputTag("", detectorFH_))) {
   usesResource("TFileService");
   ahcalGeom_ = std::make_unique<AHCalGeometry>(iConfig);
 
   // now do whatever initialization is needed
   idBeams_ = (iConfig.getParameter<std::vector<int>>("idBeams"));
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCSim") << "HGCalTBAnalyzer:: SimHits = " << doSimHits_ << " Digis = " << doDigis_ << ":"
                              << sampleIndex_ << " RecHits = " << doRecHits_ << " useDets " << ifEE_ << ":" << ifFH_
                              << ":" << ifBH_ << ":" << ifBeam_ << " zFront " << zFrontEE_ << ":" << zFrontFH_ << ":"
                              << zFrontBH_ << " IdBeam " << idBeams_.size() << ":";
   for (unsigned int k = 0; k < idBeams_.size(); ++k)
     edm::LogVerbatim("HGCSim") << " [" << k << "] " << idBeams_[k];
   edm::LogVerbatim("HGCSim") << "HGCalTBAnalyzer:: DoPassive " << doPassive_ << ":" << doPassiveEE_ << ":"
                              << doPassiveHE_ << ":" << doPassiveBH_;
   edm::LogVerbatim("HGCSim") << "HGCalTBAnalyzer:: MIP conversion factors " << gev2mip200_ << ":" << gev2mip300_
                              << " Time smearing " << stoc_smear_time_200_ << ":" << stoc_smear_time_300_ << " AddP "
                              << addP_;
 #endif
   if (idBeams_.empty())
     idBeams_.push_back(1001);
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCSim") << "HGCalTBAnalyzer:: GeneratorSource = " << labelGen_;
   if (ifEE_)
     edm::LogVerbatim("HGCSim") << "HGCalTBAnalyzer:: Detector " << detectorEE_ << " with tags " << labelHitEE_ << ", "
                                << labelDigiEE_ << ", " << labelRHitEE_;
   if (ifFH_)
     edm::LogVerbatim("HGCSim") << "HGCalTBAnalyzer:: Detector " << detectorFH_ << " with tags " << labelHitFH_ << ", "
                                << labelDigiFH_ << ", " << labelRHitFH_;
   if (ifBH_)
     edm::LogVerbatim("HGCSim") << "HGCalTBAnalyzer:: Detector " << detectorBH_ << " with tags " << labelHitBH_ << ", "
                                << labelDigiFH_ << ", " << labelRHitBH_;
   if (ifBeam_)
     edm::LogVerbatim("HGCSim") << "HGCalTBAnalyzer:: Detector " << detectorBeam_ << " with tags " << labelHitBeam_;
 #endif
 }

◆ ~HGCalTBAnalyzer()

HGCalTBAnalyzer::~HGCalTBAnalyzer ( )

override

Definition at line 262 of file HGCalTBAnalyzer.cc.

262 {}

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

EE

FH

BH

CMSE

Beam

Implements edm::one::EDAnalyzerBase.

Definition at line 549 of file HGCalTBAnalyzer.cc.

References gpuClustering::adc, addP_, ahcalGeom_, analyzeDigi(), analyzePassiveHits(), analyzeRecHits(), analyzeSimHits(), analyzeSimTracks(), HGCSample::data(), detectorBeam_, detectorBH_, detectorEE_, detectorFH_, doDigis_, doPassive_, doPassiveBH_, doPassiveEE_, doPassiveHE_, doRecHits_, doSimHits_, doTree_, edm::HepMCProduct::GetEvent(), hBeam_, hgcons_, hgcPassiveBeamEnergy_, hgcPassiveBeamID_, hgcPassiveBeamName_, hgcPassiveBHEnergy_, hgcPassiveBHID_, hgcPassiveBHName_, hgcPassiveCMSEEnergy_, hgcPassiveCMSEID_, hgcPassiveCMSEName_, hgcPassiveEEEnergy_, hgcPassiveEEID_, hgcPassiveEEName_, hgcPassiveFHEnergy_, hgcPassiveFHID_, hgcPassiveFHName_, idBeams_, iEvent, ifBeam_, ifBH_, ifEE_, ifFH_, edm::HandleBase::isValid(), dqmdumpme::k, HGCalDDDConstants::layers(), AlCaHLTBitMon_ParallelJobs::p, sampleIndex_, simHitCellColBH_, simHitCellEnBeam_, simHitCellEnBH_, simHitCellEnEE_, simHitCellEnFH_, simHitCellIdBeam_, simHitCellIdBH_, simHitCellIdEE_, simHitCellIdFH_, simHitCellLayerBH_, simHitCellRowBH_, simHitCellTime15MipEE_, simHitCellTime15MipFH_, simHitCellTimeFirstHitEE_, simHitCellTimeFirstHitFH_, simHitCellTimeLastHitEE_, simHitCellTimeLastHitFH_, simHitLayEn1BH_, simHitLayEn1EE_, simHitLayEn1FH_, simHitLayEn2BH_, simHitLayEn2EE_, simHitLayEn2FH_, simHitLayEnBeam_, edm::SortedCollection< T, SORT >::size(), tok_digiEE_, tok_digiFH_, tok_hepMC_, tok_hgcPHBeam_, 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_.

                                                                                  {
   // Generator input
   const edm::Handle<edm::HepMCProduct>& evtMC = iEvent.getHandle(tok_hepMC_);
   if (!evtMC.isValid()) {
     edm::LogWarning("HGCal") << "no HepMCProduct found";
   } else {
     const HepMC::GenEvent* myGenEvent = evtMC->GetEvent();
     unsigned int k(0);
     HepMC::FourVector pxyz(0, 0, 0, 0);
     for (HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin(); p != myGenEvent->particles_end();
          ++p, ++k) {
       edm::LogVerbatim("HGCSim") << "Particle [" << k << "] with p " << (*p)->momentum().rho() << " theta "
                                  << (*p)->momentum().theta() << " phi " << (*p)->momentum().phi() << " pxyz ("
                                  << (*p)->momentum().px() << ", " << (*p)->momentum().py() << ", "
                                  << (*p)->momentum().pz() << ")";
       if (addP_) {
         pxyz.setPx(pxyz.px() + (*p)->momentum().px());
         pxyz.setPy(pxyz.py() + (*p)->momentum().py());
         pxyz.setPz(pxyz.pz() + (*p)->momentum().pz());
         pxyz.setE(pxyz.e() + (*p)->momentum().e());
       } else if (!addP_ && (k == 0)) {
         pxyz = (*p)->momentum();
       }
     }
     hBeam_->Fill(pxyz.rho());
     edm::LogVerbatim("HGCSim") << "Particle with p " << pxyz.rho() << " theta " << pxyz.theta() << " phi "
                                << pxyz.phi();
   }
 
   // Now the Simhits
   if (doSimHits_) {
     const edm::Handle<edm::SimTrackContainer>& SimTk = iEvent.getHandle(tok_simTk_);
     const edm::Handle<edm::SimVertexContainer>& SimVtx = iEvent.getHandle(tok_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();
     simHitCellColBH_.clear();
     simHitCellRowBH_.clear();
     simHitCellLayerBH_.clear();
     simHitCellTimeFirstHitEE_.clear();
     simHitCellTime15MipEE_.clear();
     simHitCellTimeLastHitEE_.clear();
     simHitCellTimeFirstHitFH_.clear();
     simHitCellTime15MipFH_.clear();
     simHitCellTimeLastHitFH_.clear();
     std::vector<PCaloHit> caloHits;
     if (ifEE_) {
       simHitLayEn1EE_ = std::vector<float>(hgcons_[0]->layers(false), 0);
       simHitLayEn2EE_ = std::vector<float>(hgcons_[0]->layers(true), 0);
       const edm::Handle<edm::PCaloHitContainer>& theCaloHitContainers = iEvent.getHandle(tok_hitsEE_);
       if (theCaloHitContainers.isValid()) {
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HGCSim") << "PcalohitContainer for " << detectorEE_ << " has " << theCaloHitContainers->size()
                                    << " hits";
 #endif
         caloHits.clear();
         caloHits.insert(caloHits.end(), theCaloHitContainers->begin(), theCaloHitContainers->end());
         analyzeSimHits(0, caloHits, zFrontEE_);
       } else {
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HGCSim") << "PCaloHitContainer does not exist for " << detectorEE_ << " !!!";
 #endif
       }
     }
     if (ifFH_) {
       simHitLayEn1FH_ = std::vector<float>(hgcons_[1]->layers(false), 0);
       simHitLayEn2FH_ = std::vector<float>(hgcons_[1]->layers(true), 0);
       const edm::Handle<edm::PCaloHitContainer>& theCaloHitContainers = iEvent.getHandle(tok_hitsFH_);
       if (theCaloHitContainers.isValid()) {
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HGCSim") << "PcalohitContainer for " << detectorFH_ << " has " << theCaloHitContainers->size()
                                    << " hits";
 #endif
         caloHits.clear();
         caloHits.insert(caloHits.end(), theCaloHitContainers->begin(), theCaloHitContainers->end());
         analyzeSimHits(1, caloHits, zFrontFH_);
       } else {
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HGCSim") << "PCaloHitContainer does not exist for " << detectorFH_ << " !!!";
 #endif
       }
     }
     if (ifBH_) {
       simHitLayEn1BH_ = std::vector<float>(ahcalGeom_->maxDepth(), 0);
       simHitLayEn2BH_ = std::vector<float>(ahcalGeom_->maxDepth(), 0);
       const edm::Handle<edm::PCaloHitContainer>& theCaloHitContainers = iEvent.getHandle(tok_hitsBH_);
       if (theCaloHitContainers.isValid()) {
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HGCSim") << "PcalohitContainer for " << detectorBH_ << " has " << theCaloHitContainers->size()
                                    << " hits";
 #endif
         caloHits.clear();
         caloHits.insert(caloHits.end(), theCaloHitContainers->begin(), theCaloHitContainers->end());
         analyzeSimHits(2, caloHits, zFrontBH_);
       } else {
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HGCSim") << "PCaloHitContainer does not exist for " << detectorBH_ << " !!!";
 #endif
       }
     }
     if (ifBeam_) {
       simHitLayEnBeam_ = std::vector<float>(idBeams_.size(), 0);
       const edm::Handle<edm::PCaloHitContainer>& theCaloHitContainers = iEvent.getHandle(tok_hitsBeam_);
       if (theCaloHitContainers.isValid()) {
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HGCSim") << "PcalohitContainer for " << detectorBeam_ << " has "
                                    << theCaloHitContainers->size() << " hits";
 #endif
         caloHits.clear();
         caloHits.insert(caloHits.end(), theCaloHitContainers->begin(), theCaloHitContainers->end());
         analyzeSimHits(3, caloHits, 0.0);
       } else {
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HGCSim") << "PCaloHitContainer does not exist for " << detectorBeam_ << " !!!";
 #endif
       }
     }
   }  // if (doSimHits_)
 
   if (doPassive_) {
     hgcPassiveEEEnergy_.clear();
     hgcPassiveEEName_.clear();
     hgcPassiveEEID_.clear();
     if (doPassiveEE_) {
       const edm::Handle<edm::PassiveHitContainer>& hgcPHEE = iEvent.getHandle(tok_hgcPHEE_);
       analyzePassiveHits(hgcPHEE, 1);
     }
     hgcPassiveFHEnergy_.clear();
     hgcPassiveFHName_.clear();
     hgcPassiveFHID_.clear();
     if (doPassiveHE_) {
       const edm::Handle<edm::PassiveHitContainer>& hgcPHFH = iEvent.getHandle(tok_hgcPHFH_);
       analyzePassiveHits(hgcPHFH, 2);
     }
     hgcPassiveBHEnergy_.clear();
     hgcPassiveBHName_.clear();
     hgcPassiveBHID_.clear();
     if (doPassiveBH_) {
       const edm::Handle<edm::PassiveHitContainer>& hgcPHBH = iEvent.getHandle(tok_hgcPHBH_);
       analyzePassiveHits(hgcPHBH, 3);
     }
     hgcPassiveCMSEEnergy_.clear();
     hgcPassiveCMSEName_.clear();
     hgcPassiveCMSEID_.clear();
     const edm::Handle<edm::PassiveHitContainer>& hgcPHCMSE = iEvent.getHandle(tok_hgcPHCMSE_);
     analyzePassiveHits(hgcPHCMSE, 4);
     hgcPassiveBeamName_.clear();
     hgcPassiveBeamEnergy_.clear();
     hgcPassiveBeamID_.clear();
     const edm::Handle<edm::PassiveHitContainer>& hgcPHBeam = iEvent.getHandle(tok_hgcPHBeam_);
     analyzePassiveHits(hgcPHBeam, 5);
   }
 
   if ((doSimHits_ || doPassive_) && (doTree_))
     tree_->Fill();
 
   // Now the Digis
   if (doDigis_) {
     if (ifEE_) {
       const edm::Handle<HGCalDigiCollection>& theDigiContainers = iEvent.getHandle(tok_digiEE_);
       if (theDigiContainers.isValid()) {
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HGCSim") << "HGCDigiCintainer for " << detectorEE_ << " with " << theDigiContainers->size()
                                    << " element(s)";
 #endif
         for (const auto& it : *theDigiContainers) {
           HGCalDetId detId = (it.id());
           const HGCSample& hgcSample = it.sample(sampleIndex_);
           uint16_t adc = hgcSample.data();
           analyzeDigi(0, detId, adc);
         }
       }
     }
     if (ifFH_) {
       const edm::Handle<HGCalDigiCollection>& theDigiContainers = iEvent.getHandle(tok_digiFH_);
       if (theDigiContainers.isValid()) {
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HGCSim") << "HGCDigiContainer for " << detectorFH_ << " with " << theDigiContainers->size()
                                    << " element(s)";
 #endif
         for (const auto& it : *theDigiContainers) {
           HGCalDetId detId = (it.id());
           const HGCSample& hgcSample = it.sample(sampleIndex_);
           uint16_t adc = hgcSample.data();
           analyzeDigi(1, detId, adc);
         }
       }
     }
   }
 
   // The Rechits
   if (doRecHits_) {
     if (ifEE_) {
       const edm::Handle<HGCRecHitCollection>& theCaloHitContainers = iEvent.getHandle(tok_hitrEE_);
       if (theCaloHitContainers.isValid()) {
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HGCSim") << "HGCRecHitCollection for " << detectorEE_ << " has "
                                    << theCaloHitContainers->size() << " hits";
 #endif
         analyzeRecHits(0, theCaloHitContainers);
       } else {
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HGCSim") << "HGCRecHitCollection does not exist for " << detectorEE_ << " !!!";
 #endif
       }
     }
     if (ifFH_) {
       const edm::Handle<HGCRecHitCollection>& theCaloHitContainers = iEvent.getHandle(tok_hitrFH_);
       if (theCaloHitContainers.isValid()) {
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HGCSim") << "HGCRecHitCollection for " << detectorFH_ << " has "
                                    << theCaloHitContainers->size() << " hits";
 #endif
         analyzeRecHits(1, theCaloHitContainers);
       } else {
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HGCSim") << "HGCRecHitCollection does not exist for " << detectorFH_ << " !!!";
 #endif
       }  // else
     }    // if (ifFH_)
   }      // if (doRecHits_)
 
 }  // void HGCalTBAnalyzer::analyze

◆ analyzeDigi()

template<class T1 >

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

private

Definition at line 1197 of file HGCalTBAnalyzer.cc.

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

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

◆ analyzePassiveHits()

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

private

Definition at line 1255 of file HGCalTBAnalyzer.cc.

References HCALHighEnergyHPDFilter_cfi::energy, hgcPassiveBeamEnergy_, hgcPassiveBeamID_, hgcPassiveBeamName_, hgcPassiveBHEnergy_, hgcPassiveBHID_, hgcPassiveBHName_, hgcPassiveCMSEEnergy_, hgcPassiveCMSEID_, hgcPassiveCMSEName_, hgcPassiveEEEnergy_, hgcPassiveEEID_, hgcPassiveEEName_, hgcPassiveFHEnergy_, hgcPassiveFHID_, hgcPassiveFHName_, triggerObjects_cff::id, Skims_PA_cff::name, AlCaHLTBitMon_QueryRunRegistry::string, protons_cff::time, and findQualityFiles::v.

Referenced by analyze().

                                                                                                    {
   for (const 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();
     edm::LogVerbatim("HGCSim") << "HGCalTBAnalyzer::analyzePassiveHits:Energy:"
                                << "Time:Name:Id : " << energy << ":" << time << ":" << name << ":" << id;
 #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);
     } else if (subdet == 5) {
       hgcPassiveBeamEnergy_.push_back(energy);
       hgcPassiveBeamName_.push_back(name);
       hgcPassiveBeamID_.push_back(id);
     }
   }
 }

◆ analyzeRecHits()

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

private

Definition at line 1205 of file HGCalTBAnalyzer.cc.

References HGCalDetId::cell(), HCALHighEnergyHPDFilter_cfi::energy, first, HGCalGeometry::getPosition(), hgcons_, hgeom_, hfClusterShapes_cfi::hits, hRecHitE_, hRecHitLat_, hRecHitLng1_, hRecHitLng_, hRecHitOcc_, HGCalDetId::layer(), phase1PixelTopology::layer, edm::second(), HGCalDDDConstants::waferTypeL(), 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 (const 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();
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCSim") << "Layer thickness " << hgcons_[type]->waferTypeL(HGCalDetId(detId).wafer());
 #endif
     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
     edm::LogVerbatim("HGCSim") << "RecHit: " << layer << " " << global.x() << " " << global.y() << " " << global.z()
                                << " " << energy;
 #endif
   }
 
   for (const auto& itr : map_hitLayer) {
     int layer = itr.first;
     double energy = itr.second;
     double zp = hgcons_[type]->waferZ(layer, true);
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCSim") << "SimHit:Layer " << layer << " " << zp << " " << energy;
 #endif
     hRecHitLng_[type]->Fill(zp, energy);
     hRecHitLng1_[type]->Fill(layer, energy);
   }
 
   for (const 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);
   }
 }

◆ analyzeSimHits()

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

private

storing time of each pcalohit in each cell ///SJ

now sort the vector of each cell hits

id: reco and ID[id]: sim ID

once it is sorted now start adding the energy

Definition at line 793 of file HGCalTBAnalyzer.cc.

References ahcalGeom_, HGCalDDDConstants::cellThickness(), cuy::col, debug, TauDecayModes::dec, electronHcalTowerIsolationLcone_cfi::Depth, AHCalDetId::depth(), LEDCalibrationChannels::depth, relativeConstraints::empty, mps_fire::end, HCALHighEnergyHPDFilter_cfi::energy, first, gev2mip200_, gev2mip300_, HcalOther, hgcons_, hfClusterShapes_cfi::hits, hSimHitE_, hSimHitEn_, hSimHitLat_, hSimHitLayEn1BH_, hSimHitLayEn1EE_, hSimHitLayEn1FH_, hSimHitLayEn2BH_, hSimHitLayEn2EE_, hSimHitLayEn2FH_, hSimHitLayEnBeam_, hSimHitLng1_, hSimHitLng2_, hSimHitLng_, hSimHitT_, mps_fire::i, AHCalDetId::icol(), AHCalDetId::icolAbs(), triggerObjects_cff::id, idBeams_, heavyIonCSV_trainingSettings::idx, AHCalDetId::irow(), AHCalDetId::irowAbs(), dqmdumpme::k, HGCalDetId::layer(), phase1PixelTopology::layer, HGCalDDDConstants::locateCell(), nhits, HGCalTestNumbering::packHexagonIndex(), edm::second(), simHitCellColBH_, simHitCellEnBeam_, simHitCellEnBH_, simHitCellEnEE_, simHitCellEnFH_, simHitCellIdBeam_, simHitCellIdBH_, simHitCellIdEE_, simHitCellIdFH_, simHitCellLayerBH_, simHitCellRowBH_, simHitCellTime15MipEE_, simHitCellTime15MipFH_, simHitCellTimeFirstHitEE_, simHitCellTimeFirstHitFH_, simHitCellTimeLastHitEE_, simHitCellTimeLastHitFH_, simHitLayEn1BH_, simHitLayEn1EE_, simHitLayEn1FH_, simHitLayEn2BH_, simHitLayEn2EE_, simHitLayEn2FH_, simHitLayEnBeam_, HGCalDDDConstants::simToReco(), findQualityFiles::size, jetUpdater_cfi::sort, sortTime(), stoc_smear_time_200_, stoc_smear_time_300_, HcalDetId::subdet(), Calorimetry_cff::thickness, remoteMonitoring_LASER_era2018_cfg::threshold, protons_cff::time, HGCalTestNumbering::unpackHexagonIndex(), HcalTestBeamNumbering::unpackIndex(), HGCalDetId::wafer(), HGCalDDDConstants::waferZ(), geometryCSVtoXML::xx, geometryCSVtoXML::xy, AHCalDetId::zside(), and ecaldqm::zside().

Referenced by analyze().

                                                                                        {
   std::map<uint32_t, double> map_hits, map_hitn;
   std::map<uint32_t, double> map_hittime_firsthit, map_hittime_lasthit, map_hittime_15Mip;
   std::map<int, double> map_hitDepth, map_hitWafer;
   std::map<int, std::pair<uint32_t, double>> map_hitLayer, map_hitCell;
   double entot(0);
   std::map<uint32_t, double> nhits;
   std::map<uint32_t, int> ID, Depth;
   std::map<uint32_t, double> GeV2Mip;
   std::map<uint32_t, double> StochTermTime;
   std::vector<int> nSimLayers;
   std::map<uint32_t, std::vector<std::pair<double, double>>> map_hitTimeEn;
   //bool debug = true;
   bool debug = false;
   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()));
       if (debug)
         edm::LogVerbatim("HGCSim") << "depth, sector, cell " << depth << ":" << sector << ":" << cell;
     } 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::pair<float, float> xy = hgcons_[type]->locateCell(cell, layer, sector, false);
       edm::LogVerbatim("HGCSim") << "HGCalTBAnalyzer::detId " << std::hex << id << std::dec << " Layer:Wafer:Cell "
                                  << layer << ":" << sector << ":" << cell << " Position " << xy.first << ":"
                                  << xy.second << ":" << hgcons_[type]->waferZ(layer, false);
 #endif
     }
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCSim") << "SimHit:Hit[" << i << "] Id " << subdet << ":" << zside << ":" << layer << ":"
                                << sector << ":" << subsector << ":" << cell << ":" << depth << " Energy " << energy
                                << " Time " << time;
 #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::make_pair(id, ee);
     } else {
       map_hitLayer[layer] = std::make_pair(id, energy);
     }
     if (map_hitWafer.count(sector) != 0)
       map_hitWafer[sector] += energy;
     else
       map_hitWafer[sector] = energy;
     if (depth >= 0) {
       if (map_hitCell.count(idx) != 0) {
         double ee = energy + map_hitCell[idx].second;
         map_hitCell[idx] = std::make_pair(id, ee);
       } else {
         map_hitCell[idx] = std::make_pair(id, energy);
       }
       if (debug) {
         if (type == 1)
           edm::LogVerbatim("HGCSim") << "EE, depth is and map_hitDepth[depth] " << depth << " " << map_hitDepth[depth];
         if (type == 2)
           edm::LogVerbatim("HGCSim") << "BH, depth is and map_hitDepth[depth] " << depth << " " << map_hitDepth[depth];
       }
       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);
       map_hitTimeEn[idn].push_back(std::make_pair(time, energy));
       GeV2Mip[idn] = gev2mip300_;
       StochTermTime[idn] = stoc_smear_time_300_;
       ID[idn] = id;
       Depth[idn] = depth;
       if (map_hitn.count(idn) != 0) {
         map_hitn[idn] += energy;
         ++nhits[idn];
       } else {
         map_hitn[idn] = energy;
         nhits[idn] = 1;
       }
     }
     hSimHitT_[type]->Fill(time, energy);
   }
 
   if (type < 2) {  //store only for EE and FH
     edm::LogVerbatim("HGCSim") << "HGCalTAnalyzer:: " << map_hitWafer.size() << " wafers are hit in type " << type;
     for (auto itr = map_hitWafer.begin(); itr != map_hitWafer.end(); ++itr)
       edm::LogVerbatim("HGCSim") << "Wafer: " << itr->first << " Deposited Energy " << itr->second;
     for (const auto& itr : map_hitTimeEn) {
       uint32_t id = itr.first;
       int wafer = -99;
       wafer = HGCalDetId(ID[id]).wafer();
       double layer = HGCalDetId(id).layer();
       double thickness = hgcons_[type]->cellThickness(layer, wafer, 0);
       if (debug)
         edm::LogVerbatim("HGCSim") << "wafer is : depth (reco) " << wafer << " " << Depth[id]
                                    << "\ntype : layer : wafer thickness " << type << " " << layer << " " << thickness
                                    << "\nID(sim) and id(reco) " << std::hex << ID[id] << " " << id << std::dec;
       if (thickness == 300) {
         GeV2Mip[id] = gev2mip300_;
         StochTermTime[id] = stoc_smear_time_300_;
       } else if (thickness == 200) {
         GeV2Mip[id] = gev2mip200_;
         StochTermTime[id] = stoc_smear_time_200_;
       }
       //first sort
       std::sort(map_hitTimeEn[id].begin(), map_hitTimeEn[id].end(), sortTime);
       double threshold = 15.;
       double totE = 0.;
       double totEbeforeThreshold = 0.;
       double timebeforeThreshold = 0.;
       double timeAtThresohld = 0.;
       for (unsigned int ihit = 0; ihit < map_hitTimeEn[id].size(); ihit++) {
         double energy = (map_hitTimeEn[id].at(ihit)).second / GeV2Mip[id];
         totE += energy;
         double time = (map_hitTimeEn[id].at(ihit)).first;
         if (debug)
           edm::LogVerbatim("HGCSim") << "Tot E till now : time of that E : GeV2Mip[id] is " << totE << " " << time
                                      << " " << GeV2Mip[id];
         if (totE < threshold) {
           totEbeforeThreshold = totE;
           timebeforeThreshold = time;
         } else {
           timeAtThresohld =
               (threshold - totEbeforeThreshold) * (time - timebeforeThreshold) / (totE - totEbeforeThreshold) +
               timebeforeThreshold;
           map_hittime_15Mip[id] = timeAtThresohld;
           if (debug)
             edm::LogVerbatim("HGCSim") << "ihit : energyBefore : timeBefore : energyTot : timeTot : timeAt15MIP  "
                                        << ihit << " " << totEbeforeThreshold << " " << timebeforeThreshold << " "
                                        << totE << " " << time << " " << map_hittime_15Mip[id];
           break;
         }
       }
       if (!map_hitTimeEn[id].empty()) {
         map_hittime_firsthit[id] = (map_hitTimeEn[id].at(0)).first;
         map_hittime_lasthit[id] = (map_hitTimeEn[id].at(map_hitTimeEn[id].size() - 1)).first;
         if (map_hittime_15Mip[id] < map_hittime_firsthit[id])
           map_hittime_15Mip[id] = map_hittime_firsthit[id];
         /*
         double firsthit_sm = ran3.Gaus(map_hittime_firsthit_vtxCorr[id], StochTermTime[id]);
         double lasthit_sm = ran3.Gaus(map_hittime_lasthit_vtxCorr[id], StochTermTime[id]);
         double threshmiphit_sm = ran3.Gaus(map_hittime_15Mip_vtxCorr[id], StochTermTime[id]);
         */
       }
 
       if (totE < threshold)
         map_hittime_15Mip[id] = -99;
       if (debug)
         edm::LogVerbatim("HGCSim") << "id : first hit time :  last hit time " << id << " " << map_hittime_firsthit[id]
                                    << " " << map_hittime_lasthit[id] << "\nFinally for this cell, time is "
                                    << map_hittime_15Mip[id];
     }
   }
 
   hSimHitEn_[type]->Fill(entot);
   for (const auto& itr : map_hits) {
     hSimHitE_[type]->Fill(itr.second);
   }
 
   if (debug)
     edm::LogVerbatim("HGCSim") << "Now looping over map_hitLayer";
   for (const auto& itr : map_hitLayer) {
     int layer = (type == 2) ? itr.first : (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 = ahcalGeom_->getZ(AHCalDetId((itr.second).first));
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCSim") << "SimHit:Layer " << layer + 1 << " Z " << zp << ":" << zp - zFront << " E " << energy;
 #endif
     if (type < 3) {
       hSimHitLng_[type]->Fill(zp - zFront, energy);
       hSimHitLng2_[type]->Fill(layer + 1, energy);
     }
     if (type == 0) {
       if (layer < static_cast<int>(hSimHitLayEn1EE_.size())) {
         simHitLayEn1EE_[layer] = energy;
         hSimHitLayEn1EE_[layer]->Fill(energy);
       }
     } else if (type == 1) {
       if (layer < static_cast<int>(hSimHitLayEn1FH_.size())) {
         simHitLayEn1FH_[layer] = energy;
         hSimHitLayEn1FH_[layer]->Fill(energy);
       }
     } else if (type == 2) {
       if (debug)
         edm::LogVerbatim("HGCSim") << "layer < hSimHitLayEn1BH_.size()";
       if (layer < static_cast<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 (const auto& itr : map_hitDepth) {
     int layer = (type == 2) ? itr.first : (itr.first - 1);
     double energy = itr.second;
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCSim") << "SimHit:Layer " << layer + 1 << " " << energy;
 #endif
     hSimHitLng1_[type]->Fill(layer + 1, energy);
     if (type == 0) {
       if (layer < static_cast<int>(hSimHitLayEn2EE_.size())) {
         simHitLayEn2EE_[layer] = energy;
         hSimHitLayEn2EE_[layer]->Fill(energy);
       }
     } else if (type == 1) {
       if (layer < static_cast<int>(hSimHitLayEn2FH_.size())) {
         simHitLayEn2FH_[layer] = energy;
         hSimHitLayEn2FH_[layer]->Fill(energy);
       }
     } else if (type == 2) {
       if (debug)
         edm::LogVerbatim("HGCSim") << "Inside map_hitDepth, layer no. " << layer;
       if (layer < static_cast<int>(hSimHitLayEn2BH_.size())) {
         simHitLayEn2BH_[layer] = energy;
         hSimHitLayEn2BH_[layer]->Fill(energy);
       }
     }
   }
 
   if (debug)
     edm::LogVerbatim("HGCSim") << "Now looping over map_hitCell";
   if (type < 3) {
     for (const 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 = ahcalGeom_->getXY(AHCalDetId(id));
         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 (const auto& itr : map_hitn) {
     uint32_t id = itr.first;
     double energy = itr.second;
     if (type == 0) {
       double time_firsthit = map_hittime_firsthit[id];
       double time15Mip = map_hittime_15Mip[id];
       double time_lasthit = map_hittime_lasthit[id];
       simHitCellIdEE_.push_back(id);
       simHitCellEnEE_.push_back(energy);
       simHitCellTimeFirstHitEE_.push_back(time_firsthit);
       simHitCellTime15MipEE_.push_back(time15Mip);
       simHitCellTimeLastHitEE_.push_back(time_lasthit);
       if (debug && (energy / GeV2Mip[id] < 15) && (map_hittime_15Mip[id] > 0))
         edm::LogVerbatim("HGCSim") << "FOUND!!!!rechit energy : Finally for this cell, time is " << energy / GeV2Mip[id]
                                    << " " << map_hittime_15Mip[id];
     } else if (type == 1) {
       double time_firsthit = map_hittime_firsthit[id];
       double time15Mip = map_hittime_15Mip[id];
       double time_lasthit = map_hittime_lasthit[id];
       simHitCellIdFH_.push_back(id);
       simHitCellEnFH_.push_back(energy);
       simHitCellTimeFirstHitFH_.push_back(time_firsthit);
       simHitCellTime15MipFH_.push_back(time15Mip);
       simHitCellTimeLastHitFH_.push_back(time_lasthit);
     } else if (type == 2) {
       simHitCellIdBH_.push_back(id);
       simHitCellEnBH_.push_back(energy);
       AHCalDetId hid(id);
       int row = hid.irow();
       int col = hid.icol();
       int layer = hid.depth();
       simHitCellColBH_.push_back(col);
       simHitCellRowBH_.push_back(row);
       simHitCellLayerBH_.push_back(layer);
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HGCSim") << "ID: " << std::hex << id << std::dec << " Layer: " << layer << " col: " << col
                                  << " row: " << row;
 #endif
     } else if (type == 3) {
       simHitCellIdBeam_.push_back(id);
       simHitCellEnBeam_.push_back(energy);
     }
   }
 }

◆ analyzeSimTracks()

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

private

Definition at line 1117 of file HGCalTBAnalyzer.cc.

References addP_, doBeam_, gpuVertexFinder::iv, M_PI, nBeamMC_, pBeam_, pBeamMC_, pdgIdBeamMC_, phiBeam_, pxBeamMC_, pyBeamMC_, pzBeamMC_, thetaBeam_, xBeam_, xBeamMC_, yBeam_, yBeamMC_, zBeam_, and zBeamMC_.

Referenced by analyze().

                                                                                        {
   xBeam_ = yBeam_ = zBeam_ = pBeam_ = -9999;
   nBeamMC_ = thetaBeam_ = phiBeam_ = -9999;
   int nParBeam = 0;
   int vertIndex(-1);
   if (doBeam_) {
     pdgIdBeamMC_.clear();
     xBeamMC_.clear();
     yBeamMC_.clear();
     zBeamMC_.clear();
     pxBeamMC_.clear();
     pyBeamMC_.clear();
     pzBeamMC_.clear();
     pBeamMC_.clear();
   }
   std::vector<float> verX, verY, verZ;
   verX.clear();
   verY.clear();
   verZ.clear();
   for (const auto& simVtxItr : *SimVtx) {
     verX.push_back(simVtxItr.position().X());
     verY.push_back(simVtxItr.position().Y());
     verZ.push_back(simVtxItr.position().Z());
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCSim") << "Size of track " << SimTk->size();
 #endif
   HepMC::FourVector pxyz(0, 0, 0, 0);
   for (const auto& simTrkItr : *SimTk) {
     if (addP_ && !(simTrkItr.noGenpart())) {
       pxyz.setPx(pxyz.px() + simTrkItr.momentum().px());
       pxyz.setPy(pxyz.py() + simTrkItr.momentum().py());
       pxyz.setPz(pxyz.pz() + simTrkItr.momentum().pz());
       pxyz.setE(pxyz.e() + simTrkItr.momentum().e());
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HGCSim") << "Track " << simTrkItr.trackId() << " Vertex " << simTrkItr.vertIndex() << " Type "
                                  << simTrkItr.type() << " Charge " << simTrkItr.charge() << " px "
                                  << simTrkItr.momentum().px() << " py " << simTrkItr.momentum().py() << " pz "
                                  << simTrkItr.momentum().pz() << " P " << simTrkItr.momentum().P() << " GenIndex "
                                  << simTrkItr.genpartIndex();
       edm::LogVerbatim("HGCSim") << "Vertex " << simTrkItr.vertIndex()
                                  << " position-> X: " << verX[simTrkItr.vertIndex()]
                                  << " Y: " << verY[simTrkItr.vertIndex()] << " Z: " << verZ[simTrkItr.vertIndex()];
 #endif
     }
     if (doBeam_ && !(simTrkItr.noGenpart())) {
       nParBeam++;
       pdgIdBeamMC_.push_back(simTrkItr.type());
       xBeamMC_.push_back(verX[simTrkItr.vertIndex()]);
       yBeamMC_.push_back(verY[simTrkItr.vertIndex()]);
       zBeamMC_.push_back(verZ[simTrkItr.vertIndex()]);
       pxBeamMC_.push_back(simTrkItr.momentum().px());
       pyBeamMC_.push_back(simTrkItr.momentum().py());
       pzBeamMC_.push_back(simTrkItr.momentum().pz());
       pBeamMC_.push_back(simTrkItr.momentum().P());
     } else if (!addP_ && (vertIndex == -1)) {
       pxyz = simTrkItr.momentum();
     }
     if (vertIndex == -1)
       vertIndex = simTrkItr.vertIndex();
   }
   nBeamMC_ = nParBeam;
   pBeam_ = pxyz.rho();
   thetaBeam_ = pxyz.theta();
   phiBeam_ = pxyz.phi();
   if (phiBeam_ < 0)
     phiBeam_ += (2 * M_PI);
   if (vertIndex != -1 && vertIndex < static_cast<int>(SimVtx->size())) {
     edm::SimVertexContainer::const_iterator simVtxItr = SimVtx->begin();
     for (int iv = 0; iv < vertIndex; iv++)
       simVtxItr++;
     edm::LogVerbatim("HGCSim") << "Vertex " << vertIndex << " position " << simVtxItr->position();
     xBeam_ = verX[0];
     yBeam_ = verY[0];
     zBeam_ = verZ[0];
   }
 }

◆ beginJob()

void HGCalTBAnalyzer::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 321 of file HGCalTBAnalyzer.cc.

References detectorBeam_, detectorBH_, detectorEE_, detectorFH_, doBeam_, doDigis_, doPassive_, doPassiveBH_, doPassiveEE_, doPassiveHE_, doRecHits_, doSimHits_, doTree_, doTreeCell_, fs_, hBeam_, hDigiADC_, hDigiLng_, hDigiOcc_, hgcPassiveBeamEnergy_, hgcPassiveBeamID_, hgcPassiveBeamName_, 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(), Skims_PA_cff::name, nBeamMC_, pBeam_, pBeamMC_, pdgIdBeamMC_, phiBeam_, pxBeamMC_, pyBeamMC_, pzBeamMC_, simHitCellColBH_, simHitCellEnBeam_, simHitCellEnBH_, simHitCellEnEE_, simHitCellEnFH_, simHitCellIdBeam_, simHitCellIdBH_, simHitCellIdEE_, simHitCellIdFH_, simHitCellLayerBH_, simHitCellRowBH_, simHitCellTime15MipEE_, simHitCellTime15MipFH_, simHitCellTimeFirstHitEE_, simHitCellTimeFirstHitFH_, simHitCellTimeLastHitEE_, simHitCellTimeLastHitFH_, simHitLayEn1BH_, simHitLayEn1EE_, simHitLayEn1FH_, simHitLayEn2BH_, simHitLayEn2EE_, simHitLayEn2FH_, AlCaHLTBitMon_QueryRunRegistry::string, thetaBeam_, runGCPTkAlMap::title, tree_, xBeam_, xBeamMC_, yBeam_, yBeamMC_, zBeam_, and zBeamMC_.

                                {
   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");
     tree_->Branch("thetaBeam", &thetaBeam_, "thetaBeam/D");
     tree_->Branch("phiBeam", &phiBeam_, "phiBeam/D");
     if (doBeam_) {
       tree_->Branch("nBeamMC", &nBeamMC_, "nBeamMC/I");
       tree_->Branch("pdgIdBeamMC", &pdgIdBeamMC_);
       tree_->Branch("xBeamMC", &xBeamMC_);
       tree_->Branch("yBeamMC", &yBeamMC_);
       tree_->Branch("zBeamMC", &zBeamMC_);
       tree_->Branch("pxBeamMC", &pxBeamMC_);
       tree_->Branch("pyBeamMC", &pyBeamMC_);
       tree_->Branch("pzBeamMC", &pzBeamMC_);
       tree_->Branch("pBeamMC", &pBeamMC_);
     }
     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_);
 
       tree_->Branch("simHitCellColBH", &simHitCellColBH_);
       tree_->Branch("simHitCellRowBH", &simHitCellRowBH_);
       tree_->Branch("simHitCellLayerBH", &simHitCellLayerBH_);
       tree_->Branch("simHitCellTimeFirstHitEE", &simHitCellTimeFirstHitEE_);
       tree_->Branch("simHitCellTimeFirstHitFH", &simHitCellTimeFirstHitFH_);
       //tree_->Branch("simHitCellTimeFirstHitBH", &simHitCellTimeFirstHitBH_);
       tree_->Branch("simHitCellTime15MipEE", &simHitCellTime15MipEE_);
       tree_->Branch("simHitCellTime15MipFH", &simHitCellTime15MipFH_);
       //tree_->Branch("simHitCellTime15MipBH", &simHitCellTime15MipBH_);
       tree_->Branch("simHitCellTimeLastHitEE", &simHitCellTimeLastHitEE_);
       tree_->Branch("simHitCellTimeLastHitFH", &simHitCellTimeLastHitFH_);
       //tree_->Branch("simHitCellTimeLastHitBH", &simHitCellTimeLastHitBH_);
     }
   }
 
   if (doPassive_ && doTree_) {
     if (doPassiveEE_) {
       tree_->Branch("hgcPassiveEEEnergy", &hgcPassiveEEEnergy_);
       tree_->Branch("hgcPassiveEEName", &hgcPassiveEEName_);
       tree_->Branch("hgcPassiveEEID", &hgcPassiveEEID_);
     }
     if (doPassiveHE_) {
       tree_->Branch("hgcPassiveFHEnergy", &hgcPassiveFHEnergy_);
       tree_->Branch("hgcPassiveFHName", &hgcPassiveFHName_);
       tree_->Branch("hgcPassiveFHID", &hgcPassiveFHID_);
     }
     if (doPassiveBH_) {
       tree_->Branch("hgcPassiveBHEnergy", &hgcPassiveBHEnergy_);
       tree_->Branch("hgcPassiveBHName", &hgcPassiveBHName_);
       tree_->Branch("hgcPassiveBHID", &hgcPassiveBHID_);
     }
     tree_->Branch("hgcPassiveCMSEEnergy", &hgcPassiveCMSEEnergy_);
     tree_->Branch("hgcPassiveCMSEName", &hgcPassiveCMSEName_);
     tree_->Branch("hgcPassiveCMSEID", &hgcPassiveCMSEID_);
     tree_->Branch("hgcPassiveBeamEnergy", &hgcPassiveBeamEnergy_);
     tree_->Branch("hgcPassiveBeamName", &hgcPassiveBeamName_);
     tree_->Branch("hgcPassiveBeamID", &hgcPassiveBeamID_);
   }
 }

◆ beginRun()

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

overrideprivate

Definition at line 475 of file HGCalTBAnalyzer.cc.

References ahcalGeom_, detectorBeam_, detectorBH_, detectorEE_, detectorFH_, doDigis_, doRecHits_, fs_, edm::EventSetup::getData(), hgcons_, hgeom_, hSimHitLayEn1BH_, hSimHitLayEn1EE_, hSimHitLayEn1FH_, hSimHitLayEn2BH_, hSimHitLayEn2EE_, hSimHitLayEn2FH_, hSimHitLayEnBeam_, idBeams_, ifBeam_, ifBH_, ifEE_, ifFH_, cmsLHEtoEOSManager::l, HGCalDDDConstants::layers(), TFileService::make(), Skims_PA_cff::name, runGCPTkAlMap::title, tokDDDEE_, tokDDDFH_, tokGeomEE_, and tokGeomFH_.

                                                                          {
   char name[40], title[100];
   if (ifEE_) {
     hgcons_[0] = &iSetup.getData(tokDDDEE_);
     if (doDigis_ || doRecHits_) {
       hgeom_[0] = &iSetup.getData(tokGeomEE_);
     } else {
       hgeom_[0] = nullptr;
     }
     for (unsigned int l = 0; l < hgcons_[0]->layers(false); ++l) {
       sprintf(name, "SimHitEnA%d%s", l, detectorEE_.c_str());
       sprintf(title, "Sim Hit Energy in SIM layer %d for %s", l + 1, detectorEE_.c_str());
       hSimHitLayEn1EE_.push_back(fs_->make<TH1D>(name, title, 100000, 0., 0.2));
       if (l % 3 == 0) {
         sprintf(name, "SimHitEnB%d%s", (l / 3 + 1), detectorEE_.c_str());
         sprintf(title, "Sim Hit Energy in layer %d for %s", (l / 3 + 1), detectorEE_.c_str());
         hSimHitLayEn2EE_.push_back(fs_->make<TH1D>(name, title, 100000, 0., 0.2));
       }
     }
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCSim") << "HGCalTBAnalyzer::" << detectorEE_ << " defined with " << hgcons_[0]->layers(false)
                                << " layers";
 #endif
   } else {
     hgcons_[0] = nullptr;
     hgeom_[0] = nullptr;
   }
 
   if (ifFH_) {
     hgcons_[1] = &iSetup.getData(tokDDDFH_);
     if (doDigis_ || doRecHits_) {
       hgeom_[1] = &iSetup.getData(tokGeomFH_);
     } else {
       hgeom_[1] = nullptr;
     }
     for (unsigned int l = 0; l < hgcons_[1]->layers(false); ++l) {
       sprintf(name, "SimHitEnA%d%s", l, detectorFH_.c_str());
       sprintf(title, "Sim Hit Energy in layer %d for %s", l + 1, detectorFH_.c_str());
       hSimHitLayEn1FH_.push_back(fs_->make<TH1D>(name, title, 100000, 0., 0.2));
       if (l % 3 == 0) {
         sprintf(name, "SimHitEnB%d%s", (l / 3 + 1), detectorFH_.c_str());
         sprintf(title, "Sim Hit Energy in layer %d for %s", (l / 3 + 1), detectorFH_.c_str());
         hSimHitLayEn2FH_.push_back(fs_->make<TH1D>(name, title, 100000, 0., 0.2));
       }
     }
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCSim") << "HGCalTBAnalyzer::" << detectorFH_ << " defined with " << hgcons_[1]->layers(false)
                                << " layers";
 #endif
   } else {
     hgcons_[1] = nullptr;
     hgeom_[1] = nullptr;
   }
 
   if (ifBH_) {
     for (int l = 0; l < ahcalGeom_->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));
     }
   }
 }

◆ endRun()

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

inlineoverrideprivate

Definition at line 64 of file HGCalTBAnalyzer.cc.

64 {}

◆ fillDescriptions()

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

static

Definition at line 264 of file HGCalTBAnalyzer.cc.

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

                                                                                  {
   edm::ParameterSetDescription desc;
   desc.add<std::string>("detectorEE", "HGCalEESensitive");
   desc.add<bool>("useEE", true);
   desc.add<double>("zFrontEE", 0.0);
   desc.add<std::string>("caloHitSrcEE", "HGCHitsEE");
   desc.add<edm::InputTag>("digiSrcEE", edm::InputTag("hgcalDigis", "EE"));
   desc.add<edm::InputTag>("recHitSrcEE", edm::InputTag("HGCalRecHit", "HGCEERecHits"));
   desc.add<std::string>("detectorFH", "HGCalHESiliconSensitive");
   desc.add<bool>("useFH", false);
   desc.add<double>("zFrontFH", 0.0);
   desc.add<std::string>("caloHitSrcFH", "HGCHitsHEfront");
   desc.add<edm::InputTag>("digiSrcFH", edm::InputTag("hgcalDigis", "HEfront"));
   desc.add<edm::InputTag>("recHitSrcFH", edm::InputTag("HGCalRecHit", "HGCHEFRecHits"));
   desc.add<std::string>("detectorBH", "AHCal");
   desc.add<bool>("useBH", false);
   desc.add<double>("zFrontBH", 0.0);
   desc.add<std::string>("caloHitSrcBH", "HcalHits");
   desc.add<edm::InputTag>("digiSrcBH", edm::InputTag("hgcalDigis", "HEback"));
   desc.add<edm::InputTag>("recHitSrcBH", edm::InputTag("HGCalRecHit", "HGCHEBRecHits"));
   desc.add<std::string>("detectorBeam", "HcalTB06BeamDetector");
   desc.add<bool>("useBeam", false);
   desc.add<std::string>("caloHitSrcBeam", "HcalTB06BeamHits");
   std::vector<int> ids = {
       1000, 1001, 1002, 1003, 1004, 1005, 1006, 1007, 1008, 1011, 1012, 1013, 1014, 2001, 2002, 2003, 2004, 2005};
   desc.add<std::vector<int>>("idBeams", ids);
   desc.add<edm::InputTag>("generatorSrc", edm::InputTag("generatorSmeared"));
   desc.add<edm::InputTag>("passiveEE", edm::InputTag("g4SimHits", "HGCalEEPassiveHits"));
   desc.add<edm::InputTag>("passiveFH", edm::InputTag("g4SimHits", "HGCalHEPassiveHits"));
   desc.add<edm::InputTag>("passiveBH", edm::InputTag("g4SimHits", "HGCalAHPassiveHits"));
   desc.add<edm::InputTag>("passiveCMSE", edm::InputTag("g4SimHits", "CMSEPassiveHits"));
   desc.add<edm::InputTag>("passiveBeam", edm::InputTag("g4SimHits", "HGCalBeamPassiveHits"));
 
   desc.add<bool>("doSimHits", true);
   desc.add<bool>("doDigis", true);
   desc.add<bool>("doRecHits", true);
   desc.add<int>("sampleIndex", 0);
   desc.add<bool>("doTree", true);
   desc.add<bool>("doTreeCell", true);
   desc.add<bool>("doPassive", false);
   desc.add<bool>("doPassiveEE", false);
   desc.add<bool>("doPassiveHE", false);
   desc.add<bool>("doPassiveBH", false);
   desc.add<bool>("addP", false);
   desc.add<bool>("doBeam", false);
   desc.addUntracked<double>("gev2mip200", 57.0e-6);
   desc.addUntracked<double>("gev2mip300", 85.5e-6);
   desc.addUntracked<double>("stoc_smear_time_200", 10.24);
   desc.addUntracked<double>("stoc_smear_time_300", 15.5);
   desc.addUntracked<int>("maxDepth", 12);
   desc.addUntracked<double>("deltaX", 30.0);  // Size of tile along X
   desc.addUntracked<double>("deltaY", 30.0);  // Size of tile along Y
   desc.addUntracked<double>("deltaZ", 81.0);  // Thickness of a single layer
   desc.addUntracked<double>("zFirst", 17.6);  // Position of the center
   descriptions.add("HGCalTBAnalyzer", desc);
 }

◆ sortTime()

bool HGCalTBAnalyzer::sortTime	(	const std::pair< double, double > &	i,
		const std::pair< double, double > &	j
	)

staticprivate

Definition at line 1290 of file HGCalTBAnalyzer.cc.

References mps_fire::i, and dqmiolumiharvest::j.

Referenced by analyzeSimHits().

                                                                                                  {
   return i.first < j.first;
 }

Member Data Documentation

◆ addP_

const bool HGCalTBAnalyzer::addP_

private

Definition at line 82 of file HGCalTBAnalyzer.cc.

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

◆ ahcalGeom_

std::unique_ptr<AHCalGeometry> HGCalTBAnalyzer::ahcalGeom_

private

Definition at line 76 of file HGCalTBAnalyzer.cc.

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

◆ detectorBeam_

const std::string HGCalTBAnalyzer::detectorBeam_

private

Definition at line 84 of file HGCalTBAnalyzer.cc.

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

◆ detectorBH_

const std::string HGCalTBAnalyzer::detectorBH_

private

Definition at line 84 of file HGCalTBAnalyzer.cc.

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

◆ detectorEE_

const std::string HGCalTBAnalyzer::detectorEE_

private

Definition at line 83 of file HGCalTBAnalyzer.cc.

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

◆ detectorFH_

const std::string HGCalTBAnalyzer::detectorFH_

private

Definition at line 83 of file HGCalTBAnalyzer.cc.

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

◆ doBeam_

const bool HGCalTBAnalyzer::doBeam_

private

Definition at line 82 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ doDigis_

const bool HGCalTBAnalyzer::doDigis_

private

Definition at line 80 of file HGCalTBAnalyzer.cc.

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

◆ doPassive_

const bool HGCalTBAnalyzer::doPassive_

private

Definition at line 82 of file HGCalTBAnalyzer.cc.

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

◆ doPassiveBH_

const bool HGCalTBAnalyzer::doPassiveBH_

private

Definition at line 82 of file HGCalTBAnalyzer.cc.

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

◆ doPassiveEE_

const bool HGCalTBAnalyzer::doPassiveEE_

private

Definition at line 82 of file HGCalTBAnalyzer.cc.

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

◆ doPassiveHE_

const bool HGCalTBAnalyzer::doPassiveHE_

private

Definition at line 82 of file HGCalTBAnalyzer.cc.

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

◆ doRecHits_

const bool HGCalTBAnalyzer::doRecHits_

private

Definition at line 80 of file HGCalTBAnalyzer.cc.

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

◆ doSimHits_

const bool HGCalTBAnalyzer::doSimHits_

private

Definition at line 80 of file HGCalTBAnalyzer.cc.

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

◆ doTree_

const bool HGCalTBAnalyzer::doTree_

private

Definition at line 81 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ doTreeCell_

const bool HGCalTBAnalyzer::doTreeCell_

private

Definition at line 81 of file HGCalTBAnalyzer.cc.

Referenced by beginJob().

◆ fs_

edm::Service<TFileService> HGCalTBAnalyzer::fs_

private

Definition at line 75 of file HGCalTBAnalyzer.cc.

Referenced by beginJob(), and beginRun().

◆ gev2mip200_

const double HGCalTBAnalyzer::gev2mip200_

private

Definition at line 87 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and HGCalTBAnalyzer().

◆ gev2mip300_

const double HGCalTBAnalyzer::gev2mip300_

private

Definition at line 87 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and HGCalTBAnalyzer().

◆ hBeam_

TH1D * HGCalTBAnalyzer::hBeam_

private

Definition at line 113 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ hDigiADC_

TH1D* HGCalTBAnalyzer::hDigiADC_[3]

private

Definition at line 112 of file HGCalTBAnalyzer.cc.

Referenced by analyzeDigi(), and beginJob().

◆ hDigiLng_

TH1D * HGCalTBAnalyzer::hDigiLng_[2]

private

Definition at line 112 of file HGCalTBAnalyzer.cc.

Referenced by analyzeDigi(), and beginJob().

◆ hDigiOcc_

TH2D* HGCalTBAnalyzer::hDigiOcc_[3]

private

Definition at line 114 of file HGCalTBAnalyzer.cc.

Referenced by analyzeDigi(), and beginJob().

◆ hgcons_

const HGCalDDDConstants* HGCalTBAnalyzer::hgcons_[2]

private

Definition at line 77 of file HGCalTBAnalyzer.cc.

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

◆ hgcPassiveBeamEnergy_

std::vector<float> HGCalTBAnalyzer::hgcPassiveBeamEnergy_

private

Definition at line 137 of file HGCalTBAnalyzer.cc.

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

◆ hgcPassiveBeamID_

std::vector<int> HGCalTBAnalyzer::hgcPassiveBeamID_

private

Definition at line 140 of file HGCalTBAnalyzer.cc.

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

◆ hgcPassiveBeamName_

std::vector<std::string> HGCalTBAnalyzer::hgcPassiveBeamName_

private

Definition at line 139 of file HGCalTBAnalyzer.cc.

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

◆ hgcPassiveBHEnergy_

std::vector<float> HGCalTBAnalyzer::hgcPassiveBHEnergy_

private

Definition at line 136 of file HGCalTBAnalyzer.cc.

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

◆ hgcPassiveBHID_

std::vector<int> HGCalTBAnalyzer::hgcPassiveBHID_

private

Definition at line 140 of file HGCalTBAnalyzer.cc.

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

◆ hgcPassiveBHName_

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

private

Definition at line 138 of file HGCalTBAnalyzer.cc.

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

◆ hgcPassiveCMSEEnergy_

std::vector<float> HGCalTBAnalyzer::hgcPassiveCMSEEnergy_

private

Definition at line 136 of file HGCalTBAnalyzer.cc.

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

◆ hgcPassiveCMSEID_

std::vector<int> HGCalTBAnalyzer::hgcPassiveCMSEID_

private

Definition at line 140 of file HGCalTBAnalyzer.cc.

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

◆ hgcPassiveCMSEName_

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

private

Definition at line 138 of file HGCalTBAnalyzer.cc.

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

◆ hgcPassiveEEEnergy_

std::vector<float> HGCalTBAnalyzer::hgcPassiveEEEnergy_

private

Definition at line 136 of file HGCalTBAnalyzer.cc.

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

◆ hgcPassiveEEID_

std::vector<int> HGCalTBAnalyzer::hgcPassiveEEID_

private

Definition at line 140 of file HGCalTBAnalyzer.cc.

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

◆ hgcPassiveEEName_

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

private

Definition at line 138 of file HGCalTBAnalyzer.cc.

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

◆ hgcPassiveFHEnergy_

std::vector<float> HGCalTBAnalyzer::hgcPassiveFHEnergy_

private

Definition at line 136 of file HGCalTBAnalyzer.cc.

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

◆ hgcPassiveFHID_

std::vector<int> HGCalTBAnalyzer::hgcPassiveFHID_

private

Definition at line 140 of file HGCalTBAnalyzer.cc.

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

◆ hgcPassiveFHName_

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

private

Definition at line 138 of file HGCalTBAnalyzer.cc.

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

◆ hgeom_

const HGCalGeometry* HGCalTBAnalyzer::hgeom_[2]

private

Definition at line 78 of file HGCalTBAnalyzer.cc.

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

◆ hRecHitE_

TH1D* HGCalTBAnalyzer::hRecHitE_[3]

private

Definition at line 113 of file HGCalTBAnalyzer.cc.

Referenced by analyzeRecHits(), and beginJob().

◆ hRecHitLat_

TProfile2D * HGCalTBAnalyzer::hRecHitLat_[3]

private

Definition at line 118 of file HGCalTBAnalyzer.cc.

Referenced by analyzeRecHits(), and beginJob().

◆ hRecHitLng1_

TProfile * HGCalTBAnalyzer::hRecHitLng1_[3]

private

Definition at line 117 of file HGCalTBAnalyzer.cc.

Referenced by analyzeRecHits(), and beginJob().

◆ hRecHitLng_

TProfile* HGCalTBAnalyzer::hRecHitLng_[3]

private

Definition at line 117 of file HGCalTBAnalyzer.cc.

Referenced by analyzeRecHits(), and beginJob().

◆ hRecHitOcc_

TH2D * HGCalTBAnalyzer::hRecHitOcc_[3]

private

Definition at line 114 of file HGCalTBAnalyzer.cc.

Referenced by analyzeRecHits(), and beginJob().

◆ hSimHitE_

TH1D* HGCalTBAnalyzer::hSimHitE_[4]

private

Definition at line 111 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

◆ hSimHitEn_

TH1D * HGCalTBAnalyzer::hSimHitEn_[4]

private

Definition at line 113 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

◆ hSimHitLat_

TProfile2D* HGCalTBAnalyzer::hSimHitLat_[3]

private

Definition at line 118 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

◆ hSimHitLayEn1BH_

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

private

Definition at line 121 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

◆ hSimHitLayEn1EE_

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

private

Definition at line 119 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

◆ hSimHitLayEn1FH_

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

private

Definition at line 120 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

◆ hSimHitLayEn2BH_

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

private

Definition at line 121 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

◆ hSimHitLayEn2EE_

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

private

Definition at line 119 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

◆ hSimHitLayEn2FH_

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

private

Definition at line 120 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

◆ hSimHitLayEnBeam_

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

private

Definition at line 122 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

◆ hSimHitLng1_

TProfile * HGCalTBAnalyzer::hSimHitLng1_[3]

private

Definition at line 115 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

◆ hSimHitLng2_

TProfile* HGCalTBAnalyzer::hSimHitLng2_[3]

private

Definition at line 116 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

◆ hSimHitLng_

TProfile* HGCalTBAnalyzer::hSimHitLng_[3]

private

Definition at line 115 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

◆ hSimHitT_

TH1D * HGCalTBAnalyzer::hSimHitT_[4]

private

Definition at line 111 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

◆ idBeams_

std::vector<int> HGCalTBAnalyzer::idBeams_

private

Definition at line 88 of file HGCalTBAnalyzer.cc.

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

◆ ifBeam_

const bool HGCalTBAnalyzer::ifBeam_

private

Definition at line 79 of file HGCalTBAnalyzer.cc.

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

◆ ifBH_

const bool HGCalTBAnalyzer::ifBH_

private

Definition at line 79 of file HGCalTBAnalyzer.cc.

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

◆ ifEE_

const bool HGCalTBAnalyzer::ifEE_

private

Definition at line 79 of file HGCalTBAnalyzer.cc.

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

◆ ifFH_

const bool HGCalTBAnalyzer::ifFH_

private

Definition at line 79 of file HGCalTBAnalyzer.cc.

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

◆ labelDigiBH_

const edm::InputTag HGCalTBAnalyzer::labelDigiBH_

private

Definition at line 91 of file HGCalTBAnalyzer.cc.

◆ labelDigiEE_

const edm::InputTag HGCalTBAnalyzer::labelDigiEE_

private

Definition at line 91 of file HGCalTBAnalyzer.cc.

Referenced by HGCalTBAnalyzer().

◆ labelDigiFH_

const edm::InputTag HGCalTBAnalyzer::labelDigiFH_

private

Definition at line 91 of file HGCalTBAnalyzer.cc.

Referenced by HGCalTBAnalyzer().

◆ labelGen_

const edm::InputTag HGCalTBAnalyzer::labelGen_

private

Definition at line 89 of file HGCalTBAnalyzer.cc.

Referenced by HGCalTBAnalyzer().

◆ labelHitBeam_

const std::string HGCalTBAnalyzer::labelHitBeam_

private

Definition at line 90 of file HGCalTBAnalyzer.cc.

Referenced by HGCalTBAnalyzer().

◆ labelHitBH_

const std::string HGCalTBAnalyzer::labelHitBH_

private

Definition at line 90 of file HGCalTBAnalyzer.cc.

Referenced by HGCalTBAnalyzer().

◆ labelHitEE_

const std::string HGCalTBAnalyzer::labelHitEE_

private

Definition at line 90 of file HGCalTBAnalyzer.cc.

Referenced by HGCalTBAnalyzer().

◆ labelHitFH_

const std::string HGCalTBAnalyzer::labelHitFH_

private

Definition at line 90 of file HGCalTBAnalyzer.cc.

Referenced by HGCalTBAnalyzer().

◆ labelPassiveBeam_

const edm::InputTag HGCalTBAnalyzer::labelPassiveBeam_

private

Definition at line 94 of file HGCalTBAnalyzer.cc.

◆ labelPassiveBH_

const edm::InputTag HGCalTBAnalyzer::labelPassiveBH_

private

Definition at line 93 of file HGCalTBAnalyzer.cc.

◆ labelPassiveCMSE_

const edm::InputTag HGCalTBAnalyzer::labelPassiveCMSE_

private

Definition at line 94 of file HGCalTBAnalyzer.cc.

◆ labelPassiveEE_

const edm::InputTag HGCalTBAnalyzer::labelPassiveEE_

private

Definition at line 93 of file HGCalTBAnalyzer.cc.

◆ labelPassiveFH_

const edm::InputTag HGCalTBAnalyzer::labelPassiveFH_

private

Definition at line 93 of file HGCalTBAnalyzer.cc.

◆ labelRHitBH_

const edm::InputTag HGCalTBAnalyzer::labelRHitBH_

private

Definition at line 92 of file HGCalTBAnalyzer.cc.

Referenced by HGCalTBAnalyzer().

◆ labelRHitEE_

const edm::InputTag HGCalTBAnalyzer::labelRHitEE_

private

Definition at line 92 of file HGCalTBAnalyzer.cc.

Referenced by HGCalTBAnalyzer().

◆ labelRHitFH_

const edm::InputTag HGCalTBAnalyzer::labelRHitFH_

private

Definition at line 92 of file HGCalTBAnalyzer.cc.

Referenced by HGCalTBAnalyzer().

◆ nBeamMC_

int HGCalTBAnalyzer::nBeamMC_

private

Definition at line 144 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ pBeam_

double HGCalTBAnalyzer::pBeam_

private

Definition at line 142 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ pBeamMC_

std::vector<float> HGCalTBAnalyzer::pBeamMC_

private

Definition at line 147 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ pdgIdBeamMC_

std::vector<int> HGCalTBAnalyzer::pdgIdBeamMC_

private

Definition at line 145 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ phiBeam_

double HGCalTBAnalyzer::phiBeam_

private

Definition at line 143 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ pxBeamMC_

std::vector<float> HGCalTBAnalyzer::pxBeamMC_

private

Definition at line 147 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ pyBeamMC_

std::vector<float> HGCalTBAnalyzer::pyBeamMC_

private

Definition at line 147 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ pzBeamMC_

std::vector<float> HGCalTBAnalyzer::pzBeamMC_

private

Definition at line 147 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ sampleIndex_

const int HGCalTBAnalyzer::sampleIndex_

private

Definition at line 86 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

◆ simHitCellColBH_

std::vector<int> HGCalTBAnalyzer::simHitCellColBH_

private

Definition at line 131 of file HGCalTBAnalyzer.cc.

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

◆ simHitCellEnBeam_

std::vector<float> HGCalTBAnalyzer::simHitCellEnBeam_

private

Definition at line 130 of file HGCalTBAnalyzer.cc.

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

◆ simHitCellEnBH_

std::vector<float> HGCalTBAnalyzer::simHitCellEnBH_

private

Definition at line 130 of file HGCalTBAnalyzer.cc.

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

◆ simHitCellEnEE_

std::vector<float> HGCalTBAnalyzer::simHitCellEnEE_

private

Definition at line 129 of file HGCalTBAnalyzer.cc.

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

◆ simHitCellEnFH_

std::vector<float> HGCalTBAnalyzer::simHitCellEnFH_

private

Definition at line 129 of file HGCalTBAnalyzer.cc.

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

◆ simHitCellIdBeam_

std::vector<uint32_t> HGCalTBAnalyzer::simHitCellIdBeam_

private

Definition at line 128 of file HGCalTBAnalyzer.cc.

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

◆ simHitCellIdBH_

std::vector<uint32_t> HGCalTBAnalyzer::simHitCellIdBH_

private

Definition at line 128 of file HGCalTBAnalyzer.cc.

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

◆ simHitCellIdEE_

std::vector<uint32_t> HGCalTBAnalyzer::simHitCellIdEE_

private

Definition at line 127 of file HGCalTBAnalyzer.cc.

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

◆ simHitCellIdFH_

std::vector<uint32_t> HGCalTBAnalyzer::simHitCellIdFH_

private

Definition at line 127 of file HGCalTBAnalyzer.cc.

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

◆ simHitCellLayerBH_

std::vector<int> HGCalTBAnalyzer::simHitCellLayerBH_

private

Definition at line 131 of file HGCalTBAnalyzer.cc.

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

◆ simHitCellRowBH_

std::vector<int> HGCalTBAnalyzer::simHitCellRowBH_

private

Definition at line 131 of file HGCalTBAnalyzer.cc.

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

◆ simHitCellTime15MipBH_

std::vector<float> HGCalTBAnalyzer::simHitCellTime15MipBH_

private

Definition at line 133 of file HGCalTBAnalyzer.cc.

◆ simHitCellTime15MipEE_

std::vector<float> HGCalTBAnalyzer::simHitCellTime15MipEE_

private

Definition at line 133 of file HGCalTBAnalyzer.cc.

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

◆ simHitCellTime15MipFH_

std::vector<float> HGCalTBAnalyzer::simHitCellTime15MipFH_

private

Definition at line 133 of file HGCalTBAnalyzer.cc.

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

◆ simHitCellTimeFirstHitBH_

std::vector<float> HGCalTBAnalyzer::simHitCellTimeFirstHitBH_

private

Definition at line 132 of file HGCalTBAnalyzer.cc.

◆ simHitCellTimeFirstHitEE_

std::vector<float> HGCalTBAnalyzer::simHitCellTimeFirstHitEE_

private

Definition at line 132 of file HGCalTBAnalyzer.cc.

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

◆ simHitCellTimeFirstHitFH_

std::vector<float> HGCalTBAnalyzer::simHitCellTimeFirstHitFH_

private

Definition at line 132 of file HGCalTBAnalyzer.cc.

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

◆ simHitCellTimeLastHitBH_

std::vector<float> HGCalTBAnalyzer::simHitCellTimeLastHitBH_

private

Definition at line 134 of file HGCalTBAnalyzer.cc.

◆ simHitCellTimeLastHitEE_

std::vector<float> HGCalTBAnalyzer::simHitCellTimeLastHitEE_

private

Definition at line 134 of file HGCalTBAnalyzer.cc.

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

◆ simHitCellTimeLastHitFH_

std::vector<float> HGCalTBAnalyzer::simHitCellTimeLastHitFH_

private

Definition at line 134 of file HGCalTBAnalyzer.cc.

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

◆ simHitLayEn1BH_

std::vector<float> HGCalTBAnalyzer::simHitLayEn1BH_

private

Definition at line 125 of file HGCalTBAnalyzer.cc.

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

◆ simHitLayEn1EE_

std::vector<float> HGCalTBAnalyzer::simHitLayEn1EE_

private

Definition at line 123 of file HGCalTBAnalyzer.cc.

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

◆ simHitLayEn1FH_

std::vector<float> HGCalTBAnalyzer::simHitLayEn1FH_

private

Definition at line 124 of file HGCalTBAnalyzer.cc.

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

◆ simHitLayEn2BH_

std::vector<float> HGCalTBAnalyzer::simHitLayEn2BH_

private

Definition at line 125 of file HGCalTBAnalyzer.cc.

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

◆ simHitLayEn2EE_

std::vector<float> HGCalTBAnalyzer::simHitLayEn2EE_

private

Definition at line 123 of file HGCalTBAnalyzer.cc.

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

◆ simHitLayEn2FH_

std::vector<float> HGCalTBAnalyzer::simHitLayEn2FH_

private

Definition at line 124 of file HGCalTBAnalyzer.cc.

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

◆ simHitLayEnBeam_

std::vector<float> HGCalTBAnalyzer::simHitLayEnBeam_

private

Definition at line 126 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and analyzeSimHits().

◆ stoc_smear_time_200_

const double HGCalTBAnalyzer::stoc_smear_time_200_

private

Definition at line 87 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and HGCalTBAnalyzer().

◆ stoc_smear_time_300_

const double HGCalTBAnalyzer::stoc_smear_time_300_

private

Definition at line 87 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimHits(), and HGCalTBAnalyzer().

◆ thetaBeam_

double HGCalTBAnalyzer::thetaBeam_

private

Definition at line 143 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ tok_digiBH_

const edm::EDGetTokenT<HGCalDigiCollection> HGCalTBAnalyzer::tok_digiBH_

private

Definition at line 101 of file HGCalTBAnalyzer.cc.

◆ tok_digiEE_

const edm::EDGetTokenT<HGCalDigiCollection> HGCalTBAnalyzer::tok_digiEE_

private

Definition at line 101 of file HGCalTBAnalyzer.cc.

Referenced by analyze().

◆ tok_digiFH_

const edm::EDGetTokenT<HGCalDigiCollection> HGCalTBAnalyzer::tok_digiFH_

private

Definition at line 101 of file HGCalTBAnalyzer.cc.

Referenced by analyze().

◆ tok_hepMC_

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

private

Definition at line 96 of file HGCalTBAnalyzer.cc.

Referenced by analyze().

◆ tok_hgcPHBeam_

const edm::EDGetTokenT<edm::PassiveHitContainer> HGCalTBAnalyzer::tok_hgcPHBeam_

private

Definition at line 104 of file HGCalTBAnalyzer.cc.

Referenced by analyze().

◆ tok_hgcPHBH_

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

private

Definition at line 104 of file HGCalTBAnalyzer.cc.

Referenced by analyze().

◆ tok_hgcPHCMSE_

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

private

Definition at line 104 of file HGCalTBAnalyzer.cc.

Referenced by analyze().

◆ tok_hgcPHEE_

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

private

Definition at line 103 of file HGCalTBAnalyzer.cc.

Referenced by analyze().

◆ tok_hgcPHFH_

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

private

Definition at line 103 of file HGCalTBAnalyzer.cc.

Referenced by analyze().

◆ tok_hitrBH_

const edm::EDGetTokenT<HGCRecHitCollection> HGCalTBAnalyzer::tok_hitrBH_

private

Definition at line 102 of file HGCalTBAnalyzer.cc.

◆ tok_hitrEE_

const edm::EDGetTokenT<HGCRecHitCollection> HGCalTBAnalyzer::tok_hitrEE_

private

Definition at line 102 of file HGCalTBAnalyzer.cc.

Referenced by analyze().

◆ tok_hitrFH_

const edm::EDGetTokenT<HGCRecHitCollection> HGCalTBAnalyzer::tok_hitrFH_

private

Definition at line 102 of file HGCalTBAnalyzer.cc.

Referenced by analyze().

◆ tok_hitsBeam_

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

private

Definition at line 100 of file HGCalTBAnalyzer.cc.

Referenced by analyze().

◆ tok_hitsBH_

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

private

Definition at line 100 of file HGCalTBAnalyzer.cc.

Referenced by analyze().

◆ tok_hitsEE_

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

private

Definition at line 99 of file HGCalTBAnalyzer.cc.

Referenced by analyze().

◆ tok_hitsFH_

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

private

Definition at line 99 of file HGCalTBAnalyzer.cc.

Referenced by analyze().

◆ tok_simTk_

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

private

Definition at line 97 of file HGCalTBAnalyzer.cc.

Referenced by analyze().

◆ tok_simVtx_

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

private

Definition at line 98 of file HGCalTBAnalyzer.cc.

Referenced by analyze().

◆ tokDDDEE_

const edm::ESGetToken<HGCalDDDConstants, IdealGeometryRecord> HGCalTBAnalyzer::tokDDDEE_

private

Definition at line 105 of file HGCalTBAnalyzer.cc.

Referenced by beginRun().

◆ tokDDDFH_

const edm::ESGetToken<HGCalDDDConstants, IdealGeometryRecord> HGCalTBAnalyzer::tokDDDFH_

private

Definition at line 107 of file HGCalTBAnalyzer.cc.

Referenced by beginRun().

◆ tokGeomEE_

const edm::ESGetToken<HGCalGeometry, IdealGeometryRecord> HGCalTBAnalyzer::tokGeomEE_

private

Definition at line 106 of file HGCalTBAnalyzer.cc.

Referenced by beginRun().

◆ tokGeomFH_

const edm::ESGetToken<HGCalGeometry, IdealGeometryRecord> HGCalTBAnalyzer::tokGeomFH_

private

Definition at line 108 of file HGCalTBAnalyzer.cc.

Referenced by beginRun().

◆ tree_

TTree* HGCalTBAnalyzer::tree_

private

Definition at line 110 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ xBeam_

double HGCalTBAnalyzer::xBeam_

private

Definition at line 142 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ xBeamMC_

std::vector<float> HGCalTBAnalyzer::xBeamMC_

private

Definition at line 146 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ yBeam_

double HGCalTBAnalyzer::yBeam_

private

Definition at line 142 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ yBeamMC_

std::vector<float> HGCalTBAnalyzer::yBeamMC_

private

Definition at line 146 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ zBeam_

double HGCalTBAnalyzer::zBeam_

private

Definition at line 142 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ zBeamMC_

std::vector<float> HGCalTBAnalyzer::zBeamMC_

private

Definition at line 146 of file HGCalTBAnalyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ zFrontBH_

const double HGCalTBAnalyzer::zFrontBH_

private

Definition at line 85 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

◆ zFrontEE_

const double HGCalTBAnalyzer::zFrontEE_

private

Definition at line 85 of file HGCalTBAnalyzer.cc.

Referenced by analyze(), and HGCalTBAnalyzer().

◆ zFrontFH_

const 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

Static Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HGCalTBAnalyzer()

◆ ~HGCalTBAnalyzer()

Member Function Documentation

◆ analyze()

◆ analyzeDigi()

◆ analyzePassiveHits()

◆ analyzeRecHits()

◆ analyzeSimHits()

◆ analyzeSimTracks()

◆ beginJob()

◆ beginRun()

◆ endRun()

◆ fillDescriptions()

◆ sortTime()

Member Data Documentation

◆ addP_

◆ ahcalGeom_

◆ detectorBeam_

◆ detectorBH_

◆ detectorEE_

◆ detectorFH_

◆ doBeam_

◆ doDigis_

◆ doPassive_

◆ doPassiveBH_

◆ doPassiveEE_

◆ doPassiveHE_

◆ doRecHits_

◆ doSimHits_

◆ doTree_

◆ doTreeCell_

◆ fs_

◆ gev2mip200_

◆ gev2mip300_

◆ hBeam_

◆ hDigiADC_

◆ hDigiLng_

◆ hDigiOcc_

◆ hgcons_

◆ hgcPassiveBeamEnergy_

◆ hgcPassiveBeamID_

◆ hgcPassiveBeamName_

◆ hgcPassiveBHEnergy_

◆ hgcPassiveBHID_

◆ hgcPassiveBHName_

◆ hgcPassiveCMSEEnergy_

◆ hgcPassiveCMSEID_

◆ hgcPassiveCMSEName_

◆ hgcPassiveEEEnergy_

◆ hgcPassiveEEID_

◆ hgcPassiveEEName_

◆ hgcPassiveFHEnergy_

◆ hgcPassiveFHID_

◆ hgcPassiveFHName_

◆ hgeom_

◆ hRecHitE_

◆ hRecHitLat_

◆ hRecHitLng1_

◆ hRecHitLng_

◆ hRecHitOcc_

◆ hSimHitE_

◆ hSimHitEn_

◆ hSimHitLat_

◆ hSimHitLayEn1BH_

◆ hSimHitLayEn1EE_

◆ hSimHitLayEn1FH_

◆ hSimHitLayEn2BH_

◆ hSimHitLayEn2EE_

◆ hSimHitLayEn2FH_

◆ hSimHitLayEnBeam_

◆ hSimHitLng1_

◆ hSimHitLng2_