Inheritance diagram for HGCalTriggerNtupleHGCClusters:

Public Member Functions
void	fill (const edm::Event &e, const edm::EventSetup &es) final

	HGCalTriggerNtupleHGCClusters (const edm::ParameterSet &conf)

void	initialize (TTree &, const edm::ParameterSet &, edm::ConsumesCollector &&) final

	~HGCalTriggerNtupleHGCClusters () override

Public Member Functions inherited from HGCalTriggerNtupleBase
	HGCalTriggerNtupleBase (const edm::ParameterSet &conf)

virtual	~HGCalTriggerNtupleBase ()

Private Member Functions
void	clear () final

Private Attributes
std::vector< std::vector< uint32_t > >	cl_cells_id_

std::vector< int >	cl_cells_n_

std::vector< float >	cl_energy_

std::vector< float >	cl_eta_

std::vector< uint32_t >	cl_id_

std::vector< int >	cl_layer_

std::vector< float >	cl_mipPt_

std::vector< uint32_t >	cl_multicluster_id_

std::vector< float >	cl_multicluster_pt_

int	cl_n_

std::vector< float >	cl_phi_

std::vector< float >	cl_pt_

std::vector< int >	cl_subdet_

edm::EDGetToken	clusters_token_

bool	filter_clusters_in_multiclusters_

edm::EDGetToken	multiclusters_token_

HGCalTriggerTools	triggerTools_

Additional Inherited Members

Detailed Description

Definition at line 13 of file HGCalTriggerNtupleHGCClusters.cc.

Constructor & Destructor Documentation

HGCalTriggerNtupleHGCClusters::HGCalTriggerNtupleHGCClusters ( const edm::ParameterSet & conf )

Definition at line 54 of file HGCalTriggerNtupleHGCClusters.cc.

References initialize().

                                                           :HGCalTriggerNtupleBase(conf),
   filter_clusters_in_multiclusters_(conf.getParameter<bool>("FilterClustersInMulticlusters"))
 {
 }

HGCalTriggerNtupleHGCClusters::~HGCalTriggerNtupleHGCClusters ( )

inlineoverride

Definition at line 18 of file HGCalTriggerNtupleHGCClusters.cc.

References clear(), MillePedeFileConverter_cfg::e, fill(), and initialize().

18 {};

Member Function Documentation

void HGCalTriggerNtupleHGCClusters::clear ( void )

finalprivatevirtual

Implements HGCalTriggerNtupleBase.

Definition at line 145 of file HGCalTriggerNtupleHGCClusters.cc.

References cl_cells_id_, cl_cells_n_, cl_energy_, cl_eta_, cl_id_, cl_layer_, cl_mipPt_, cl_multicluster_id_, cl_multicluster_pt_, cl_n_, cl_phi_, cl_pt_, and cl_subdet_.

Referenced by fill(), and ~HGCalTriggerNtupleHGCClusters().

 {
   cl_n_ = 0;
   cl_id_.clear();
   cl_mipPt_.clear();
   cl_pt_.clear();
   cl_energy_.clear();
   cl_eta_.clear();
   cl_phi_.clear();
   cl_layer_.clear();
   cl_subdet_.clear();
   cl_cells_n_.clear();
   cl_cells_id_.clear();
   cl_multicluster_id_.clear();
   cl_multicluster_pt_.clear();
 }

void HGCalTriggerNtupleHGCClusters::fill	(	const edm::Event &	e,
		const edm::EventSetup &	es
	)

finalvirtual

Implements HGCalTriggerNtupleBase.

Definition at line 83 of file HGCalTriggerNtupleHGCClusters.cc.

References BXVector< T >::begin(), cl_cells_id_, cl_cells_n_, cl_energy_, cl_eta_, cl_id_, cl_layer_, cl_mipPt_, cl_multicluster_id_, cl_multicluster_pt_, cl_n_, cl_phi_, cl_pt_, cl_subdet_, clear(), fastPrimaryVertexProducer_cfi::clusters, clusters_token_, BXVector< T >::end(), HGCalTriggerTools::eventSetup(), filter_clusters_in_multiclusters_, geometry, edm::EventSetup::get(), edm::Event::getByToken(), HGCalTriggerTools::layerWithOffset(), multiclusters_token_, create_public_lumi_plots::transform, and triggerTools_.

Referenced by initialize(), and ~HGCalTriggerNtupleHGCClusters().

 {
 
   // retrieve clusters
   edm::Handle<l1t::HGCalClusterBxCollection> clusters_h;
   e.getByToken(clusters_token_, clusters_h);
   const l1t::HGCalClusterBxCollection& clusters = *clusters_h;
   edm::Handle<l1t::HGCalMulticlusterBxCollection> multiclusters_h;
   e.getByToken(multiclusters_token_, multiclusters_h);
   const l1t::HGCalMulticlusterBxCollection& multiclusters = *multiclusters_h;
 
   // retrieve geometry
   edm::ESHandle<HGCalTriggerGeometryBase> geometry;
   es.get<CaloGeometryRecord>().get(geometry);
 
   triggerTools_.eventSetup(es);
 
   // Associate cells to clusters
   std::unordered_map<uint32_t, l1t::HGCalMulticlusterBxCollection::const_iterator> cluster2multicluster;
   for(auto mcl_itr=multiclusters.begin(0); mcl_itr!=multiclusters.end(0); mcl_itr++)
   {
     // loop on 2D clusters inside 3D clusters
     for(const auto& cl_ptr : mcl_itr->constituents())
     {
       cluster2multicluster.emplace(cl_ptr.second->detId(), mcl_itr);
     }
   }
 
 
 
   clear();
   for(auto cl_itr=clusters.begin(0); cl_itr!=clusters.end(0); cl_itr++)
   {
     auto mcl_itr = cluster2multicluster.find(cl_itr->detId());
     uint32_t mcl_id = (mcl_itr!=cluster2multicluster.end() ? mcl_itr->second->detId() : 0);
     float mcl_pt = (mcl_itr!=cluster2multicluster.end() ? mcl_itr->second->pt() : 0.);
     if(filter_clusters_in_multiclusters_ && mcl_id==0) continue;
     cl_n_++;
     cl_mipPt_.emplace_back(cl_itr->mipPt());
     // physical values 
     cl_pt_.emplace_back(cl_itr->pt());
     cl_energy_.emplace_back(cl_itr->energy());
     cl_eta_.emplace_back(cl_itr->eta());
     cl_phi_.emplace_back(cl_itr->phi());
 
     cl_id_.emplace_back(cl_itr->detId());
     cl_layer_.emplace_back(triggerTools_.layerWithOffset(cl_itr->detId()));
     cl_subdet_.emplace_back(cl_itr->subdetId());
     cl_cells_n_.emplace_back(cl_itr->constituents().size());
     // Retrieve indices of trigger cells inside cluster
     cl_cells_id_.emplace_back(cl_itr->constituents().size());
     std::transform(cl_itr->constituents_begin(), cl_itr->constituents_end(),
         cl_cells_id_.back().begin(), [](const std::pair<uint32_t,edm::Ptr<l1t::HGCalTriggerCell>>& id_tc){return id_tc.second->detId();}
         );
     cl_multicluster_id_.emplace_back(mcl_id);
     cl_multicluster_pt_.emplace_back(mcl_pt);
   }
 }

void HGCalTriggerNtupleHGCClusters::initialize	(	TTree &	tree,
		const edm::ParameterSet &	conf,
		edm::ConsumesCollector &&	collector
	)

finalvirtual

Implements HGCalTriggerNtupleBase.

Definition at line 61 of file HGCalTriggerNtupleHGCClusters.cc.

References cl_cells_id_, cl_cells_n_, cl_energy_, cl_eta_, cl_id_, cl_layer_, cl_mipPt_, cl_multicluster_id_, cl_multicluster_pt_, cl_n_, cl_phi_, cl_pt_, cl_subdet_, clusters_token_, fill(), edm::ParameterSet::getParameter(), and multiclusters_token_.

Referenced by HGCalTriggerNtupleHGCClusters(), and ~HGCalTriggerNtupleHGCClusters().

 {
   clusters_token_ = collector.consumes<l1t::HGCalClusterBxCollection>(conf.getParameter<edm::InputTag>("Clusters"));
   multiclusters_token_ = collector.consumes<l1t::HGCalMulticlusterBxCollection>(conf.getParameter<edm::InputTag>("Multiclusters"));
 
   tree.Branch("cl_n", &cl_n_, "cl_n/I");
   tree.Branch("cl_id", &cl_id_);
   tree.Branch("cl_mipPt", &cl_mipPt_);
   tree.Branch("cl_pt", &cl_pt_);
   tree.Branch("cl_energy", &cl_energy_);
   tree.Branch("cl_eta", &cl_eta_);
   tree.Branch("cl_phi", &cl_phi_);
   tree.Branch("cl_layer", &cl_layer_);
   tree.Branch("cl_subdet", &cl_subdet_);
   tree.Branch("cl_cells_n", &cl_cells_n_);
   tree.Branch("cl_cells_id", &cl_cells_id_);
   tree.Branch("cl_multicluster_id", &cl_multicluster_id_);
   tree.Branch("cl_multicluster_pt", &cl_multicluster_pt_);
 }