Inheritance diagram for HGCalTB23Analyzer:

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

	~HGCalTB23Analyzer () override=default

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

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

std::vector< ESResolverIndex > 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

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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices 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

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

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	doPassive_

const bool	doPassiveBH_

const bool	doPassiveEE_

const bool	doPassiveHE_

const bool	doSimHits_

const bool	doTree_

const bool	doTreeCell_

edm::Service< TFileService >	fs_

const double	gev2mip200_

const double	gev2mip300_

TH1D *	hBeam_

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_

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

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_

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

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 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 51 of file HGCalTB23Analyzer.cc.

Constructor & Destructor Documentation

◆ HGCalTB23Analyzer()

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

explicit

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

Definition at line 131 of file HGCalTB23Analyzer.cc.

References addP_, ahcalGeom_, detectorBeam_, detectorBH_, detectorEE_, detectorFH_, doPassive_, doPassiveBH_, doPassiveEE_, doPassiveHE_, doSimHits_, edm::ParameterSet::getParameter(), gev2mip200_, gev2mip300_, idBeams_, ifBeam_, ifBH_, ifEE_, ifFH_, dqmdumpme::k, labelGen_, labelHitBeam_, labelHitBH_, labelHitEE_, labelHitFH_, 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")),
       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")),
       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")),
       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_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_))),
       tokDDDFH_(esConsumes<HGCalDDDConstants, 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") << "HGCalTB23Analyzer:: SimHits = " << doSimHits_ << " 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") << "HGCalTB23Analyzer:: DoPassive " << doPassive_ << ":" << doPassiveEE_ << ":"
                              << doPassiveHE_ << ":" << doPassiveBH_;
   edm::LogVerbatim("HGCSim") << "HGCalTB23Analyzer:: 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") << "HGCalTB23Analyzer:: GeneratorSource = " << labelGen_;
   if (ifEE_)
     edm::LogVerbatim("HGCSim") << "HGCalTB23Analyzer:: Detector " << detectorEE_ << " with tags " << labelHitEE_;
   if (ifFH_)
     edm::LogVerbatim("HGCSim") << "HGCalTB23Analyzer:: Detector " << detectorFH_ << " with tags " << labelHitFH_;
   if (ifBH_)
     edm::LogVerbatim("HGCSim") << "HGCalTB23Analyzer:: Detector " << detectorBH_ << " with tags " << labelHitBH_;
   if (ifBeam_)
     edm::LogVerbatim("HGCSim") << "HGCalTB23Analyzer:: Detector " << detectorBeam_ << " with tags " << labelHitBeam_;
 #endif
 }

◆ ~HGCalTB23Analyzer()

HGCalTB23Analyzer::~HGCalTB23Analyzer ( )

overridedefault

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

EE

FH

BH

CMSE

Beam

Implements edm::one::EDAnalyzerBase.

Definition at line 454 of file HGCalTB23Analyzer.cc.

References addP_, ahcalGeom_, analyzePassiveHits(), analyzeSimHits(), analyzeSimTracks(), detectorBeam_, detectorBH_, detectorEE_, detectorFH_, doPassive_, doPassiveBH_, doPassiveEE_, doPassiveHE_, 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, simHitCellColBH_, simHitCellEnBeam_, simHitCellEnBH_, simHitCellEnEE_, simHitCellEnFH_, simHitCellIdBeam_, simHitCellIdBH_, simHitCellIdEE_, simHitCellIdFH_, simHitCellLayerBH_, simHitCellRowBH_, simHitCellTime15MipEE_, simHitCellTime15MipFH_, simHitCellTimeFirstHitEE_, simHitCellTimeFirstHitFH_, simHitCellTimeLastHitEE_, simHitCellTimeLastHitFH_, simHitLayEn1BH_, simHitLayEn1EE_, simHitLayEn1FH_, simHitLayEn2BH_, simHitLayEn2EE_, simHitLayEn2FH_, simHitLayEnBeam_, tok_hepMC_, tok_hgcPHBeam_, tok_hgcPHBH_, tok_hgcPHCMSE_, tok_hgcPHEE_, tok_hgcPHFH_, 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();
 }  // void HGCalTB23Analyzer::analyze

◆ analyzePassiveHits()

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

private

Definition at line 1035 of file HGCalTB23Analyzer.cc.

References hcalRecHitTable_cff::energy, hgcPassiveBeamEnergy_, hgcPassiveBeamID_, hgcPassiveBeamName_, hgcPassiveBHEnergy_, hgcPassiveBHID_, hgcPassiveBHName_, hgcPassiveCMSEEnergy_, hgcPassiveCMSEID_, hgcPassiveCMSEName_, hgcPassiveEEEnergy_, hgcPassiveEEID_, hgcPassiveEEName_, hgcPassiveFHEnergy_, hgcPassiveFHID_, hgcPassiveFHName_, l1ctLayer2EG_cff::id, Skims_PA_cff::name, AlCaHLTBitMon_QueryRunRegistry::string, hcalRecHitTable_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") << "HGCalTB23Analyzer::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);
     }
   }
 }

◆ analyzeSimHits()

void HGCalTB23Analyzer::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 631 of file HGCalTB23Analyzer.cc.

References ahcalGeom_, HGCalDDDConstants::cellThickness(), cuy::col, debug, TauDecayModes::dec, electronHcalTowerIsolationLcone_cfi::Depth, hcalRecHitTable_cff::depth, AHCalDetId::depth(), relativeConstraints::empty, mps_fire::end, hcalRecHitTable_cff::energy, dqmdumpme::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(), l1ctLayer2EG_cff::id, idBeams_, heavyIonCSV_trainingSettings::idx, AHCalDetId::irow(), AHCalDetId::irowAbs(), dqmdumpme::k, nano_mu_digi_cff::layer, HGCSiliconDetId::layer(), HGCalDDDConstants::locateCell(), nhits, HGCalTestNumbering::packHexagonIndex(), edm::second(), nano_mu_digi_cff::sector, 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(), nano_mu_digi_cff::subsector, Calorimetry_cff::thickness, remoteMonitoring_LASER_era2018_cfg::threshold, hcalRecHitTable_cff::time, HGCalTestNumbering::unpackHexagonIndex(), HcalTestBeamNumbering::unpackIndex(), HGCalWaferIndex::waferU(), HGCSiliconDetId::waferU(), HGCalWaferIndex::waferV(), HGCSiliconDetId::waferV(), 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") << "HGCalTB23Analyzer::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") << "HGCalTB23Analyzer:: " << 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 waferU(-99), waferV(-99);
       waferU = HGCSiliconDetId(ID[id]).waferU();
       waferV = HGCSiliconDetId(ID[id]).waferV();
       double layer = HGCSiliconDetId(id).layer();
       double thickness = hgcons_[type]->cellThickness(layer, waferU, waferV);
       if (debug)
         edm::LogVerbatim("HGCSim") << "wafer is : depth (reco) " << waferU << ":" << waferV << " " << 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 HGCalTB23Analyzer::analyzeSimTracks	(	edm::Handle< edm::SimTrackContainer > const &	SimTk,
		edm::Handle< edm::SimVertexContainer > const &	SimVtx
	)

private

Definition at line 956 of file HGCalTB23Analyzer.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 HGCalTB23Analyzer::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 268 of file HGCalTB23Analyzer.cc.

References detectorBeam_, detectorBH_, detectorEE_, detectorFH_, doBeam_, doPassive_, doPassiveBH_, doPassiveEE_, doPassiveHE_, doSimHits_, doTree_, doTreeCell_, fs_, hBeam_, hgcPassiveBeamEnergy_, hgcPassiveBeamID_, hgcPassiveBeamName_, hgcPassiveBHEnergy_, hgcPassiveBHID_, hgcPassiveBHName_, hgcPassiveCMSEEnergy_, hgcPassiveCMSEID_, hgcPassiveCMSEName_, hgcPassiveEEEnergy_, hgcPassiveEEID_, hgcPassiveEEName_, hgcPassiveFHEnergy_, hgcPassiveFHID_, hgcPassiveFHName_, 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 (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 HGCalTB23Analyzer::beginRun	(	edm::Run const &	,
		edm::EventSetup const &	iSetup
	)

overrideprivate

Definition at line 392 of file HGCalTB23Analyzer.cc.

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

                                                                            {
   char name[40], title[100];
   if (ifEE_) {
     hgcons_[0] = &iSetup.getData(tokDDDEE_);
     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") << "HGCalTB23Analyzer::" << detectorEE_ << " defined with " << hgcons_[0]->layers(false)
                                << " layers";
 #endif
   } else {
     hgcons_[0] = nullptr;
   }
 
   if (ifFH_) {
     hgcons_[1] = &iSetup.getData(tokDDDFH_);
     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") << "HGCalTB23Analyzer::" << detectorFH_ << " defined with " << hgcons_[1]->layers(false)
                                << " layers";
 #endif
   } else {
     hgcons_[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 HGCalTB23Analyzer::endRun	(	edm::Run const &	,
		edm::EventSetup const &
	)

inlineoverrideprivate

Definition at line 61 of file HGCalTB23Analyzer.cc.

61 {}

◆ fillDescriptions()

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

static

Definition at line 220 of file HGCalTB23Analyzer.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, ProducerED_cfi::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<std::string>("detectorFH", "HGCalHESiliconSensitive");
   desc.add<bool>("useFH", false);
   desc.add<double>("zFrontFH", 0.0);
   desc.add<std::string>("caloHitSrcFH", "HGCHitsHEfront");
   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<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>("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("HGCalTB23Analyzer", desc);
 }

◆ sortTime()

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

staticprivate

Definition at line 1070 of file HGCalTB23Analyzer.cc.

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

Referenced by analyzeSimHits().

                                                                                                    {
   return i.first < j.first;
 }

Member Data Documentation

◆ addP_

const bool HGCalTB23Analyzer::addP_

private

Definition at line 74 of file HGCalTB23Analyzer.cc.

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

◆ ahcalGeom_

std::unique_ptr<AHCalGeometry> HGCalTB23Analyzer::ahcalGeom_

private

Definition at line 70 of file HGCalTB23Analyzer.cc.

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

◆ detectorBeam_

const std::string HGCalTB23Analyzer::detectorBeam_

private

Definition at line 76 of file HGCalTB23Analyzer.cc.

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

◆ detectorBH_

const std::string HGCalTB23Analyzer::detectorBH_

private

Definition at line 76 of file HGCalTB23Analyzer.cc.

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

◆ detectorEE_

const std::string HGCalTB23Analyzer::detectorEE_

private

Definition at line 75 of file HGCalTB23Analyzer.cc.

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

◆ detectorFH_

const std::string HGCalTB23Analyzer::detectorFH_

private

Definition at line 75 of file HGCalTB23Analyzer.cc.

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

◆ doBeam_

const bool HGCalTB23Analyzer::doBeam_

private

Definition at line 74 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ doPassive_

const bool HGCalTB23Analyzer::doPassive_

private

Definition at line 74 of file HGCalTB23Analyzer.cc.

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

◆ doPassiveBH_

const bool HGCalTB23Analyzer::doPassiveBH_

private

Definition at line 74 of file HGCalTB23Analyzer.cc.

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

◆ doPassiveEE_

const bool HGCalTB23Analyzer::doPassiveEE_

private

Definition at line 74 of file HGCalTB23Analyzer.cc.

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

◆ doPassiveHE_

const bool HGCalTB23Analyzer::doPassiveHE_

private

Definition at line 74 of file HGCalTB23Analyzer.cc.

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

◆ doSimHits_

const bool HGCalTB23Analyzer::doSimHits_

private

Definition at line 73 of file HGCalTB23Analyzer.cc.

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

◆ doTree_

const bool HGCalTB23Analyzer::doTree_

private

Definition at line 73 of file HGCalTB23Analyzer.cc.

Referenced by analyze(), and beginJob().

◆ doTreeCell_

const bool HGCalTB23Analyzer::doTreeCell_

private

Definition at line 73 of file HGCalTB23Analyzer.cc.

Referenced by beginJob().

◆ fs_

edm::Service<TFileService> HGCalTB23Analyzer::fs_

private

Definition at line 69 of file HGCalTB23Analyzer.cc.

Referenced by beginJob(), and beginRun().

◆ gev2mip200_

const double HGCalTB23Analyzer::gev2mip200_

private

Definition at line 78 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimHits(), and HGCalTB23Analyzer().

◆ gev2mip300_

const double HGCalTB23Analyzer::gev2mip300_

private

Definition at line 78 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimHits(), and HGCalTB23Analyzer().

◆ hBeam_

TH1D * HGCalTB23Analyzer::hBeam_

private

Definition at line 96 of file HGCalTB23Analyzer.cc.

Referenced by analyze(), and beginJob().

◆ hgcons_

const HGCalDDDConstants* HGCalTB23Analyzer::hgcons_[2]

private

Definition at line 71 of file HGCalTB23Analyzer.cc.

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

◆ hgcPassiveBeamEnergy_

std::vector<float> HGCalTB23Analyzer::hgcPassiveBeamEnergy_

private

Definition at line 118 of file HGCalTB23Analyzer.cc.

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

◆ hgcPassiveBeamID_

std::vector<int> HGCalTB23Analyzer::hgcPassiveBeamID_

private

Definition at line 121 of file HGCalTB23Analyzer.cc.

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

◆ hgcPassiveBeamName_

std::vector<std::string> HGCalTB23Analyzer::hgcPassiveBeamName_

private

Definition at line 120 of file HGCalTB23Analyzer.cc.

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

◆ hgcPassiveBHEnergy_

std::vector<float> HGCalTB23Analyzer::hgcPassiveBHEnergy_

private

Definition at line 117 of file HGCalTB23Analyzer.cc.

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

◆ hgcPassiveBHID_

std::vector<int> HGCalTB23Analyzer::hgcPassiveBHID_

private

Definition at line 121 of file HGCalTB23Analyzer.cc.

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

◆ hgcPassiveBHName_

std::vector<std::string> HGCalTB23Analyzer::hgcPassiveBHName_

private

Definition at line 119 of file HGCalTB23Analyzer.cc.

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

◆ hgcPassiveCMSEEnergy_

std::vector<float> HGCalTB23Analyzer::hgcPassiveCMSEEnergy_

private

Definition at line 117 of file HGCalTB23Analyzer.cc.

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

◆ hgcPassiveCMSEID_

std::vector<int> HGCalTB23Analyzer::hgcPassiveCMSEID_

private

Definition at line 121 of file HGCalTB23Analyzer.cc.

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

◆ hgcPassiveCMSEName_

std::vector<std::string> HGCalTB23Analyzer::hgcPassiveCMSEName_

private

Definition at line 119 of file HGCalTB23Analyzer.cc.

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

◆ hgcPassiveEEEnergy_

std::vector<float> HGCalTB23Analyzer::hgcPassiveEEEnergy_

private

Definition at line 117 of file HGCalTB23Analyzer.cc.

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

◆ hgcPassiveEEID_

std::vector<int> HGCalTB23Analyzer::hgcPassiveEEID_

private

Definition at line 121 of file HGCalTB23Analyzer.cc.

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

◆ hgcPassiveEEName_

std::vector<std::string> HGCalTB23Analyzer::hgcPassiveEEName_

private

Definition at line 119 of file HGCalTB23Analyzer.cc.

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

◆ hgcPassiveFHEnergy_

std::vector<float> HGCalTB23Analyzer::hgcPassiveFHEnergy_

private

Definition at line 117 of file HGCalTB23Analyzer.cc.

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

◆ hgcPassiveFHID_

std::vector<int> HGCalTB23Analyzer::hgcPassiveFHID_

private

Definition at line 121 of file HGCalTB23Analyzer.cc.

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

◆ hgcPassiveFHName_

std::vector<std::string> HGCalTB23Analyzer::hgcPassiveFHName_

private

Definition at line 119 of file HGCalTB23Analyzer.cc.

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

◆ hSimHitE_

TH1D* HGCalTB23Analyzer::hSimHitE_[4]

private

Definition at line 96 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

◆ hSimHitEn_

TH1D * HGCalTB23Analyzer::hSimHitEn_[4]

private

Definition at line 96 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

◆ hSimHitLat_

TProfile2D* HGCalTB23Analyzer::hSimHitLat_[3]

private

Definition at line 99 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

◆ hSimHitLayEn1BH_

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

private

Definition at line 102 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

◆ hSimHitLayEn1EE_

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

private

Definition at line 100 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

◆ hSimHitLayEn1FH_

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

private

Definition at line 101 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

◆ hSimHitLayEn2BH_

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

private

Definition at line 102 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

◆ hSimHitLayEn2EE_

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

private

Definition at line 100 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

◆ hSimHitLayEn2FH_

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

private

Definition at line 101 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

◆ hSimHitLayEnBeam_

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

private

Definition at line 103 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimHits(), and beginRun().

◆ hSimHitLng1_

TProfile * HGCalTB23Analyzer::hSimHitLng1_[3]

private

Definition at line 97 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

◆ hSimHitLng2_

TProfile* HGCalTB23Analyzer::hSimHitLng2_[3]

private

Definition at line 98 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

◆ hSimHitLng_

TProfile* HGCalTB23Analyzer::hSimHitLng_[3]

private

Definition at line 97 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

◆ hSimHitT_

TH1D * HGCalTB23Analyzer::hSimHitT_[4]

private

Definition at line 96 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimHits(), and beginJob().

◆ idBeams_

std::vector<int> HGCalTB23Analyzer::idBeams_

private

Definition at line 79 of file HGCalTB23Analyzer.cc.

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

◆ ifBeam_

const bool HGCalTB23Analyzer::ifBeam_

private

Definition at line 72 of file HGCalTB23Analyzer.cc.

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

◆ ifBH_

const bool HGCalTB23Analyzer::ifBH_

private

Definition at line 72 of file HGCalTB23Analyzer.cc.

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

◆ ifEE_

const bool HGCalTB23Analyzer::ifEE_

private

Definition at line 72 of file HGCalTB23Analyzer.cc.

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

◆ ifFH_

const bool HGCalTB23Analyzer::ifFH_

private

Definition at line 72 of file HGCalTB23Analyzer.cc.

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

◆ labelGen_

const edm::InputTag HGCalTB23Analyzer::labelGen_

private

Definition at line 80 of file HGCalTB23Analyzer.cc.

Referenced by HGCalTB23Analyzer().

◆ labelHitBeam_

const std::string HGCalTB23Analyzer::labelHitBeam_

private

Definition at line 81 of file HGCalTB23Analyzer.cc.

Referenced by HGCalTB23Analyzer().

◆ labelHitBH_

const std::string HGCalTB23Analyzer::labelHitBH_

private

Definition at line 81 of file HGCalTB23Analyzer.cc.

Referenced by HGCalTB23Analyzer().

◆ labelHitEE_

const std::string HGCalTB23Analyzer::labelHitEE_

private

Definition at line 81 of file HGCalTB23Analyzer.cc.

Referenced by HGCalTB23Analyzer().

◆ labelHitFH_

const std::string HGCalTB23Analyzer::labelHitFH_

private

Definition at line 81 of file HGCalTB23Analyzer.cc.

Referenced by HGCalTB23Analyzer().

◆ labelPassiveBeam_

const edm::InputTag HGCalTB23Analyzer::labelPassiveBeam_

private

Definition at line 83 of file HGCalTB23Analyzer.cc.

◆ labelPassiveBH_

const edm::InputTag HGCalTB23Analyzer::labelPassiveBH_

private

Definition at line 82 of file HGCalTB23Analyzer.cc.

◆ labelPassiveCMSE_

const edm::InputTag HGCalTB23Analyzer::labelPassiveCMSE_

private

Definition at line 83 of file HGCalTB23Analyzer.cc.

◆ labelPassiveEE_

const edm::InputTag HGCalTB23Analyzer::labelPassiveEE_

private

Definition at line 82 of file HGCalTB23Analyzer.cc.

◆ labelPassiveFH_

const edm::InputTag HGCalTB23Analyzer::labelPassiveFH_

private

Definition at line 82 of file HGCalTB23Analyzer.cc.

◆ nBeamMC_

int HGCalTB23Analyzer::nBeamMC_

private

Definition at line 125 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ pBeam_

double HGCalTB23Analyzer::pBeam_

private

Definition at line 123 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ pBeamMC_

std::vector<float> HGCalTB23Analyzer::pBeamMC_

private

Definition at line 128 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ pdgIdBeamMC_

std::vector<int> HGCalTB23Analyzer::pdgIdBeamMC_

private

Definition at line 126 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ phiBeam_

double HGCalTB23Analyzer::phiBeam_

private

Definition at line 124 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ pxBeamMC_

std::vector<float> HGCalTB23Analyzer::pxBeamMC_

private

Definition at line 128 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ pyBeamMC_

std::vector<float> HGCalTB23Analyzer::pyBeamMC_

private

Definition at line 128 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ pzBeamMC_

std::vector<float> HGCalTB23Analyzer::pzBeamMC_

private

Definition at line 128 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ simHitCellColBH_

std::vector<int> HGCalTB23Analyzer::simHitCellColBH_

private

Definition at line 112 of file HGCalTB23Analyzer.cc.

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

◆ simHitCellEnBeam_

std::vector<float> HGCalTB23Analyzer::simHitCellEnBeam_

private

Definition at line 111 of file HGCalTB23Analyzer.cc.

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

◆ simHitCellEnBH_

std::vector<float> HGCalTB23Analyzer::simHitCellEnBH_

private

Definition at line 111 of file HGCalTB23Analyzer.cc.

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

◆ simHitCellEnEE_

std::vector<float> HGCalTB23Analyzer::simHitCellEnEE_

private

Definition at line 110 of file HGCalTB23Analyzer.cc.

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

◆ simHitCellEnFH_

std::vector<float> HGCalTB23Analyzer::simHitCellEnFH_

private

Definition at line 110 of file HGCalTB23Analyzer.cc.

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

◆ simHitCellIdBeam_

std::vector<uint32_t> HGCalTB23Analyzer::simHitCellIdBeam_

private

Definition at line 109 of file HGCalTB23Analyzer.cc.

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

◆ simHitCellIdBH_

std::vector<uint32_t> HGCalTB23Analyzer::simHitCellIdBH_

private

Definition at line 109 of file HGCalTB23Analyzer.cc.

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

◆ simHitCellIdEE_

std::vector<uint32_t> HGCalTB23Analyzer::simHitCellIdEE_

private

Definition at line 108 of file HGCalTB23Analyzer.cc.

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

◆ simHitCellIdFH_

std::vector<uint32_t> HGCalTB23Analyzer::simHitCellIdFH_

private

Definition at line 108 of file HGCalTB23Analyzer.cc.

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

◆ simHitCellLayerBH_

std::vector<int> HGCalTB23Analyzer::simHitCellLayerBH_

private

Definition at line 112 of file HGCalTB23Analyzer.cc.

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

◆ simHitCellRowBH_

std::vector<int> HGCalTB23Analyzer::simHitCellRowBH_

private

Definition at line 112 of file HGCalTB23Analyzer.cc.

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

◆ simHitCellTime15MipBH_

std::vector<float> HGCalTB23Analyzer::simHitCellTime15MipBH_

private

Definition at line 114 of file HGCalTB23Analyzer.cc.

◆ simHitCellTime15MipEE_

std::vector<float> HGCalTB23Analyzer::simHitCellTime15MipEE_

private

Definition at line 114 of file HGCalTB23Analyzer.cc.

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

◆ simHitCellTime15MipFH_

std::vector<float> HGCalTB23Analyzer::simHitCellTime15MipFH_

private

Definition at line 114 of file HGCalTB23Analyzer.cc.

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

◆ simHitCellTimeFirstHitBH_

std::vector<float> HGCalTB23Analyzer::simHitCellTimeFirstHitBH_

private

Definition at line 113 of file HGCalTB23Analyzer.cc.

◆ simHitCellTimeFirstHitEE_

std::vector<float> HGCalTB23Analyzer::simHitCellTimeFirstHitEE_

private

Definition at line 113 of file HGCalTB23Analyzer.cc.

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

◆ simHitCellTimeFirstHitFH_

std::vector<float> HGCalTB23Analyzer::simHitCellTimeFirstHitFH_

private

Definition at line 113 of file HGCalTB23Analyzer.cc.

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

◆ simHitCellTimeLastHitBH_

std::vector<float> HGCalTB23Analyzer::simHitCellTimeLastHitBH_

private

Definition at line 115 of file HGCalTB23Analyzer.cc.

◆ simHitCellTimeLastHitEE_

std::vector<float> HGCalTB23Analyzer::simHitCellTimeLastHitEE_

private

Definition at line 115 of file HGCalTB23Analyzer.cc.

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

◆ simHitCellTimeLastHitFH_

std::vector<float> HGCalTB23Analyzer::simHitCellTimeLastHitFH_

private

Definition at line 115 of file HGCalTB23Analyzer.cc.

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

◆ simHitLayEn1BH_

std::vector<float> HGCalTB23Analyzer::simHitLayEn1BH_

private

Definition at line 106 of file HGCalTB23Analyzer.cc.

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

◆ simHitLayEn1EE_

std::vector<float> HGCalTB23Analyzer::simHitLayEn1EE_

private

Definition at line 104 of file HGCalTB23Analyzer.cc.

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

◆ simHitLayEn1FH_

std::vector<float> HGCalTB23Analyzer::simHitLayEn1FH_

private

Definition at line 105 of file HGCalTB23Analyzer.cc.

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

◆ simHitLayEn2BH_

std::vector<float> HGCalTB23Analyzer::simHitLayEn2BH_

private

Definition at line 106 of file HGCalTB23Analyzer.cc.

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

◆ simHitLayEn2EE_

std::vector<float> HGCalTB23Analyzer::simHitLayEn2EE_

private

Definition at line 104 of file HGCalTB23Analyzer.cc.

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

◆ simHitLayEn2FH_

std::vector<float> HGCalTB23Analyzer::simHitLayEn2FH_

private

Definition at line 105 of file HGCalTB23Analyzer.cc.

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

◆ simHitLayEnBeam_

std::vector<float> HGCalTB23Analyzer::simHitLayEnBeam_

private

Definition at line 107 of file HGCalTB23Analyzer.cc.

Referenced by analyze(), and analyzeSimHits().

◆ stoc_smear_time_200_

const double HGCalTB23Analyzer::stoc_smear_time_200_

private

Definition at line 78 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimHits(), and HGCalTB23Analyzer().

◆ stoc_smear_time_300_

const double HGCalTB23Analyzer::stoc_smear_time_300_

private

Definition at line 78 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimHits(), and HGCalTB23Analyzer().

◆ thetaBeam_

double HGCalTB23Analyzer::thetaBeam_

private

Definition at line 124 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ tok_hepMC_

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

private

Definition at line 85 of file HGCalTB23Analyzer.cc.

Referenced by analyze().

◆ tok_hgcPHBeam_

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

private

Definition at line 91 of file HGCalTB23Analyzer.cc.

Referenced by analyze().

◆ tok_hgcPHBH_

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

private

Definition at line 91 of file HGCalTB23Analyzer.cc.

Referenced by analyze().

◆ tok_hgcPHCMSE_

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

private

Definition at line 91 of file HGCalTB23Analyzer.cc.

Referenced by analyze().

◆ tok_hgcPHEE_

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

private

Definition at line 90 of file HGCalTB23Analyzer.cc.

Referenced by analyze().

◆ tok_hgcPHFH_

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

private

Definition at line 90 of file HGCalTB23Analyzer.cc.

Referenced by analyze().

◆ tok_hitsBeam_

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

private

Definition at line 89 of file HGCalTB23Analyzer.cc.

Referenced by analyze().

◆ tok_hitsBH_

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

private

Definition at line 89 of file HGCalTB23Analyzer.cc.

Referenced by analyze().

◆ tok_hitsEE_

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

private

Definition at line 88 of file HGCalTB23Analyzer.cc.

Referenced by analyze().

◆ tok_hitsFH_

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

private

Definition at line 88 of file HGCalTB23Analyzer.cc.

Referenced by analyze().

◆ tok_simTk_

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

private

Definition at line 86 of file HGCalTB23Analyzer.cc.

Referenced by analyze().

◆ tok_simVtx_

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

private

Definition at line 87 of file HGCalTB23Analyzer.cc.

Referenced by analyze().

◆ tokDDDEE_

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

private

Definition at line 92 of file HGCalTB23Analyzer.cc.

Referenced by beginRun().

◆ tokDDDFH_

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

private

Definition at line 93 of file HGCalTB23Analyzer.cc.

Referenced by beginRun().

◆ tree_

TTree* HGCalTB23Analyzer::tree_

private

Definition at line 95 of file HGCalTB23Analyzer.cc.

Referenced by analyze(), and beginJob().

◆ xBeam_

double HGCalTB23Analyzer::xBeam_

private

Definition at line 123 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ xBeamMC_

std::vector<float> HGCalTB23Analyzer::xBeamMC_

private

Definition at line 127 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ yBeam_

double HGCalTB23Analyzer::yBeam_

private

Definition at line 123 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ yBeamMC_

std::vector<float> HGCalTB23Analyzer::yBeamMC_

private

Definition at line 127 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ zBeam_

double HGCalTB23Analyzer::zBeam_

private

Definition at line 123 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ zBeamMC_

std::vector<float> HGCalTB23Analyzer::zBeamMC_

private

Definition at line 127 of file HGCalTB23Analyzer.cc.

Referenced by analyzeSimTracks(), and beginJob().

◆ zFrontBH_

const double HGCalTB23Analyzer::zFrontBH_

private

Definition at line 77 of file HGCalTB23Analyzer.cc.

Referenced by analyze(), and HGCalTB23Analyzer().

◆ zFrontEE_

const double HGCalTB23Analyzer::zFrontEE_

private

Definition at line 77 of file HGCalTB23Analyzer.cc.

Referenced by analyze(), and HGCalTB23Analyzer().

◆ zFrontFH_

const double HGCalTB23Analyzer::zFrontFH_

private

Definition at line 77 of file HGCalTB23Analyzer.cc.

Referenced by analyze(), and HGCalTB23Analyzer().

Public Member Functions

Static Public Member Functions

Private Member Functions

Static Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HGCalTB23Analyzer()

◆ ~HGCalTB23Analyzer()

Member Function Documentation

◆ analyze()

◆ analyzePassiveHits()

◆ analyzeSimHits()

◆ analyzeSimTracks()

◆ beginJob()

◆ beginRun()

◆ endRun()

◆ fillDescriptions()

◆ sortTime()

Member Data Documentation

◆ addP_

◆ ahcalGeom_

◆ detectorBeam_

◆ detectorBH_

◆ detectorEE_

◆ detectorFH_

◆ doBeam_

◆ doPassive_

◆ doPassiveBH_

◆ doPassiveEE_

◆ doPassiveHE_

◆ doSimHits_

◆ doTree_

◆ doTreeCell_

◆ fs_

◆ gev2mip200_

◆ gev2mip300_

◆ hBeam_

◆ hgcons_

◆ hgcPassiveBeamEnergy_

◆ hgcPassiveBeamID_

◆ hgcPassiveBeamName_

◆ hgcPassiveBHEnergy_

◆ hgcPassiveBHID_

◆ hgcPassiveBHName_

◆ hgcPassiveCMSEEnergy_

◆ hgcPassiveCMSEID_

◆ hgcPassiveCMSEName_

◆ hgcPassiveEEEnergy_

◆ hgcPassiveEEID_

◆ hgcPassiveEEName_

◆ hgcPassiveFHEnergy_

◆ hgcPassiveFHID_

◆ hgcPassiveFHName_

◆ hSimHitE_

◆ hSimHitEn_

◆ hSimHitLat_

◆ hSimHitLayEn1BH_

◆ hSimHitLayEn1EE_

◆ hSimHitLayEn1FH_

◆ hSimHitLayEn2BH_

◆ hSimHitLayEn2EE_

◆ hSimHitLayEn2FH_

◆ hSimHitLayEnBeam_

◆ hSimHitLng1_

◆ hSimHitLng2_

◆ hSimHitLng_

◆ hSimHitT_

◆ idBeams_

◆ ifBeam_

◆ ifBH_

◆ ifEE_

◆ ifFH_

◆ labelGen_

◆ labelHitBeam_

◆ labelHitBH_

◆ labelHitEE_

◆ labelHitFH_

◆ labelPassiveBeam_