Inheritance diagram for HGCalTriggerNtupleHGCTriggerCells:

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

	HGCalTriggerNtupleHGCTriggerCells (const edm::ParameterSet &conf)

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

	~HGCalTriggerNtupleHGCTriggerCells () override

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

virtual	~HGCalTriggerNtupleBase ()

Private Member Functions
double	calibrate (double, int, int)

void	clear () final

void	simhits (const edm::Event &e, std::unordered_map< uint32_t, double > &simhits_ee, std::unordered_map< uint32_t, double > &simhits_fh, std::unordered_map< uint32_t, double > &simhits_bh)

Private Attributes
std::vector< double >	fcPerMip_

bool	fill_simenergy_

bool	filter_cells_in_multiclusters_

edm::ESHandle< HGCalTriggerGeometryBase >	geometry_

double	keV2fC_

std::vector< double >	layerWeights_

edm::EDGetToken	multiclusters_token_

edm::EDGetToken	simhits_bh_token_

edm::EDGetToken	simhits_ee_token_

edm::EDGetToken	simhits_fh_token_

std::vector< int >	tc_cell_

std::vector< uint32_t >	tc_cluster_id_

std::vector< uint32_t >	tc_data_

std::vector< float >	tc_energy_

std::vector< float >	tc_eta_

std::vector< uint32_t >	tc_id_

std::vector< int >	tc_layer_

std::vector< float >	tc_mipPt_

std::vector< uint32_t >	tc_multicluster_id_

std::vector< float >	tc_multicluster_pt_

int	tc_n_

std::vector< float >	tc_phi_

std::vector< float >	tc_pt_

std::vector< int >	tc_side_

std::vector< float >	tc_simenergy_

std::vector< int >	tc_subdet_

std::vector< int >	tc_wafer_

std::vector< int >	tc_wafertype_

std::vector< float >	tc_x_

std::vector< float >	tc_y_

std::vector< float >	tc_z_

std::vector< double >	thicknessCorrections_

edm::EDGetToken	trigger_cells_token_

HGCalTriggerTools	triggerTools_

Additional Inherited Members

Detailed Description

Definition at line 16 of file HGCalTriggerNtupleHGCTriggerCells.cc.

Constructor & Destructor Documentation

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

Definition at line 73 of file HGCalTriggerNtupleHGCTriggerCells.cc.

References fcPerMip_, fill_simenergy_, filter_cells_in_multiclusters_, edm::ParameterSet::getParameter(), initialize(), keV2fC_, layerWeights_, and thicknessCorrections_.

                                                               :HGCalTriggerNtupleBase(conf),
   fill_simenergy_(conf.getParameter<bool>("FillSimEnergy")),
   filter_cells_in_multiclusters_(conf.getParameter<bool>("FilterCellsInMulticlusters"))
 {
   fill_simenergy_ = conf.getParameter<bool>("FillSimEnergy");
   filter_cells_in_multiclusters_ = conf.getParameter<bool>("FilterCellsInMulticlusters");
   keV2fC_ = conf.getParameter<double>("keV2fC");
   fcPerMip_ = conf.getParameter<std::vector<double>>("fcPerMip");
   layerWeights_ = conf.getParameter<std::vector<double>>("layerWeights");
   thicknessCorrections_ = conf.getParameter<std::vector<double>>("thicknessCorrections");
 }

HGCalTriggerNtupleHGCTriggerCells::~HGCalTriggerNtupleHGCTriggerCells ( )

inlineoverride

Definition at line 21 of file HGCalTriggerNtupleHGCTriggerCells.cc.

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

21 {};

Member Function Documentation

double HGCalTriggerNtupleHGCTriggerCells::calibrate	(	double	energy,
		int	thickness,
		int	layer
	)

private

Definition at line 254 of file HGCalTriggerNtupleHGCTriggerCells.cc.

References hgcalLayerClusters_cfi::fcPerMip, fcPerMip_, keV2fC_, layerWeights_, simhits(), hgcalLayerClusters_cfi::thicknessCorrection, and thicknessCorrections_.

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

 {
   double fcPerMip = fcPerMip_[thickness];
   double thicknessCorrection = thicknessCorrections_[thickness];
   double layerWeight = layerWeights_[layer];
   double TeV2GeV = 1.e3;
   return energy*keV2fC_/fcPerMip*layerWeight*TeV2GeV/thicknessCorrection;
 }

void HGCalTriggerNtupleHGCTriggerCells::clear ( void )

finalprivatevirtual

Implements HGCalTriggerNtupleBase.

Definition at line 319 of file HGCalTriggerNtupleHGCTriggerCells.cc.

References tc_cell_, tc_cluster_id_, tc_data_, tc_energy_, tc_eta_, tc_id_, tc_layer_, tc_mipPt_, tc_multicluster_id_, tc_multicluster_pt_, tc_n_, tc_phi_, tc_pt_, tc_side_, tc_simenergy_, tc_subdet_, tc_wafer_, tc_wafertype_, tc_x_, tc_y_, and tc_z_.

Referenced by fill(), simhits(), and ~HGCalTriggerNtupleHGCTriggerCells().

 {
   tc_n_ = 0;
   tc_id_.clear();
   tc_subdet_.clear();
   tc_side_.clear();
   tc_layer_.clear();
   tc_wafer_.clear();
   tc_wafertype_.clear();
   tc_cell_.clear();
   tc_data_.clear();
   tc_mipPt_.clear();
   tc_pt_.clear();
   tc_energy_.clear();
   tc_simenergy_.clear();
   tc_eta_.clear();
   tc_phi_.clear();
   tc_x_.clear();
   tc_y_.clear();
   tc_z_.clear();
   tc_cluster_id_.clear();
   tc_multicluster_id_.clear();
   tc_multicluster_pt_.clear();
 }

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

finalvirtual

Implements HGCalTriggerNtupleBase.

Definition at line 125 of file HGCalTriggerNtupleHGCTriggerCells.cc.

References BXVector< T >::begin(), calibrate(), clear(), HGCalTopology::dddConstants(), HGCalTriggerGeometryBase::eeTopology(), BXVector< T >::end(), HGCalTriggerTools::eventSetup(), HGCalTriggerGeometryBase::fhTopology(), fill_simenergy_, filter_cells_in_multiclusters_, geometry_, edm::EventSetup::get(), edm::Event::getByToken(), HGCalTriggerGeometryBase::getCellsFromTriggerCell(), HGCEE, HGCHEB, HGCHEF, hcalTTPDigis_cfi::id, HGCalDetId::layer(), HGCalTriggerTools::layerWithOffset(), multiclusters_token_, simhits(), tc_cell_, tc_cluster_id_, tc_data_, tc_energy_, tc_eta_, tc_id_, tc_layer_, tc_mipPt_, tc_multicluster_id_, tc_multicluster_pt_, tc_n_, tc_phi_, tc_pt_, tc_side_, tc_simenergy_, tc_subdet_, tc_wafer_, tc_wafertype_, tc_x_, tc_y_, tc_z_, trigger_cells_token_, triggerTools_, HGCalDetId::wafer(), HGCalDDDConstants::waferTypeL(), and HcalDigiParam_cfi::zside.

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

 {
 
   // retrieve trigger cells
   edm::Handle<l1t::HGCalTriggerCellBxCollection> trigger_cells_h;
   e.getByToken(trigger_cells_token_, trigger_cells_h);
   const l1t::HGCalTriggerCellBxCollection& trigger_cells = *trigger_cells_h;
 
   // retrieve clusters
   edm::Handle<l1t::HGCalMulticlusterBxCollection> multiclusters_h;
   e.getByToken(multiclusters_token_, multiclusters_h);
   const l1t::HGCalMulticlusterBxCollection& multiclusters = *multiclusters_h;
 
   // retrieve geometry
   es.get<CaloGeometryRecord>().get(geometry_);
 
   // sim hit association
   std::unordered_map<uint32_t, double> simhits_ee;
   std::unordered_map<uint32_t, double> simhits_fh;  
   std::unordered_map<uint32_t, double> simhits_bh;  
   if(fill_simenergy_) simhits(e, simhits_ee, simhits_fh, simhits_bh);
 
   // Associate cells to clusters
   std::unordered_map<uint32_t, uint32_t> cell2cluster;
   std::unordered_map<uint32_t, l1t::HGCalMulticlusterBxCollection::const_iterator> cell2multicluster;
   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())
     {
       // loop on TC inside 2D clusters
       for(const auto& tc_ptr : cl_ptr.second->constituents())
       {
         cell2cluster.emplace(tc_ptr.second->detId(), cl_ptr.second->detId());
         cell2multicluster.emplace(tc_ptr.second->detId(), mcl_itr);
       }
     }
   }
 
   triggerTools_.eventSetup(es);
 
   clear();
   for(auto tc_itr=trigger_cells.begin(0); tc_itr!=trigger_cells.end(0); tc_itr++)
   {
     if(tc_itr->hwPt()>0)
     {
       auto cl_itr = cell2cluster.find(tc_itr->detId());
       auto mcl_itr = cell2multicluster.find(tc_itr->detId());
       uint32_t cl_id = (cl_itr!=cell2cluster.end() ? cl_itr->second : 0);
       uint32_t mcl_id = (mcl_itr!=cell2multicluster.end() ? mcl_itr->second->detId() : 0);
       float mcl_pt = (mcl_itr!=cell2multicluster.end() ? mcl_itr->second->pt() : 0.);
       // Filter cells not included in a multicluster, if requested
       if(filter_cells_in_multiclusters_ && mcl_id==0) continue;
       tc_n_++;
       // hardware data
       HGCalDetId id(tc_itr->detId());
       tc_id_.emplace_back(tc_itr->detId());
       tc_subdet_.emplace_back(id.subdetId());
       tc_side_.emplace_back(id.zside());
       tc_layer_.emplace_back(triggerTools_.layerWithOffset(id));
       tc_wafer_.emplace_back(id.wafer());
       tc_wafertype_.emplace_back(id.waferType());
       tc_cell_.emplace_back(id.cell());
       tc_data_.emplace_back(tc_itr->hwPt());
       tc_mipPt_.emplace_back(tc_itr->mipPt());
       // physical values 
       tc_pt_.emplace_back(tc_itr->pt());
       tc_energy_.emplace_back(tc_itr->energy());
       tc_eta_.emplace_back(tc_itr->eta());
       tc_phi_.emplace_back(tc_itr->phi());
       tc_x_.emplace_back(tc_itr->position().x());
       tc_y_.emplace_back(tc_itr->position().y());
       tc_z_.emplace_back(tc_itr->position().z());
       // Links between TC and clusters
       tc_cluster_id_.emplace_back(cl_id);
       tc_multicluster_id_.emplace_back(mcl_id);
       tc_multicluster_pt_.emplace_back(mcl_pt);
 
       if(fill_simenergy_)
       {
         double energy = 0;
         int subdet = id.subdetId();
         // search for simhit for all the cells inside the trigger cell
         for(uint32_t c_id : geometry_->getCellsFromTriggerCell(id))
         {
           switch(subdet)
           {
             case ForwardSubdetector::HGCEE:
               {
                 auto itr = simhits_ee.find(c_id);
                 if(itr!=simhits_ee.end())
                 {
                   HGCalDetId detid(c_id);
                   int thickness = geometry_->eeTopology().dddConstants().waferTypeL(detid.wafer())-1;
                   int layer = detid.layer();
                   energy += calibrate(itr->second, thickness, layer);
                 }
                 break;
               }
             case ForwardSubdetector::HGCHEF:
               {
                 auto itr = simhits_fh.find(c_id);
                 if(itr!=simhits_fh.end())
                 {
                   HGCalDetId detid(c_id);
                   int thickness = geometry_->fhTopology().dddConstants().waferTypeL(detid.wafer())-1;
                   int layer = detid.layer();
                   energy += calibrate(itr->second, thickness, layer);
                 }
                 break;
               }
             case ForwardSubdetector::HGCHEB:
               {
                 auto itr = simhits_bh.find(c_id);
                 if(itr!=simhits_bh.end()) energy += itr->second;
                 break;
               }
             default:
               break;
           }
         }
         tc_simenergy_.emplace_back(energy); 
       }
     }
   }
 }

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

finalvirtual

Implements HGCalTriggerNtupleBase.

Definition at line 87 of file HGCalTriggerNtupleHGCTriggerCells.cc.

References fill(), fill_simenergy_, edm::ParameterSet::getParameter(), multiclusters_token_, simhits_bh_token_, simhits_ee_token_, simhits_fh_token_, tc_cell_, tc_cluster_id_, tc_data_, tc_energy_, tc_eta_, tc_id_, tc_layer_, tc_mipPt_, tc_multicluster_id_, tc_multicluster_pt_, tc_n_, tc_phi_, tc_pt_, tc_side_, tc_simenergy_, tc_subdet_, tc_wafer_, tc_wafertype_, tc_x_, tc_y_, tc_z_, and trigger_cells_token_.

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

 {
   trigger_cells_token_ = collector.consumes<l1t::HGCalTriggerCellBxCollection>(conf.getParameter<edm::InputTag>("TriggerCells"));
   multiclusters_token_ = collector.consumes<l1t::HGCalMulticlusterBxCollection>(conf.getParameter<edm::InputTag>("Multiclusters"));
 
   if (fill_simenergy_) 
   {
     simhits_ee_token_ = collector.consumes<edm::PCaloHitContainer>(conf.getParameter<edm::InputTag>("eeSimHits"));
     simhits_fh_token_ = collector.consumes<edm::PCaloHitContainer>(conf.getParameter<edm::InputTag>("fhSimHits"));
     simhits_bh_token_ = collector.consumes<edm::PCaloHitContainer>(conf.getParameter<edm::InputTag>("bhSimHits"));
   }
 
   tree.Branch("tc_n", &tc_n_, "tc_n/I");
   tree.Branch("tc_id", &tc_id_);
   tree.Branch("tc_subdet", &tc_subdet_);
   tree.Branch("tc_zside", &tc_side_);
   tree.Branch("tc_layer", &tc_layer_);
   tree.Branch("tc_wafer", &tc_wafer_);
   tree.Branch("tc_wafertype", &tc_wafertype_);
   tree.Branch("tc_cell", &tc_cell_);
   tree.Branch("tc_data", &tc_data_);
   tree.Branch("tc_pt", &tc_pt_);
   tree.Branch("tc_mipPt", &tc_mipPt_);
   tree.Branch("tc_energy", &tc_energy_);
   if(fill_simenergy_) tree.Branch("tc_simenergy", &tc_simenergy_);
   tree.Branch("tc_eta", &tc_eta_);
   tree.Branch("tc_phi", &tc_phi_);
   tree.Branch("tc_x", &tc_x_);
   tree.Branch("tc_y", &tc_y_);
   tree.Branch("tc_z", &tc_z_);
   tree.Branch("tc_cluster_id", &tc_cluster_id_);
   tree.Branch("tc_multicluster_id", &tc_multicluster_id_);
   tree.Branch("tc_multicluster_pt", &tc_multicluster_pt_);
 
 }

void HGCalTriggerNtupleHGCTriggerCells::simhits	(	const edm::Event &	e,
		std::unordered_map< uint32_t, double > &	simhits_ee,
		std::unordered_map< uint32_t, double > &	simhits_fh,
		std::unordered_map< uint32_t, double > &	simhits_bh
	)

private

Definition at line 265 of file HGCalTriggerNtupleHGCTriggerCells.cc.

References HGCalTriggerGeometryBase::bhTopology(), clear(), HGCalTopology::dddConstants(), HcalTopology::dddConstants(), HGCalTopology::detectorType(), HGCalTriggerGeometryBase::eeTopology(), HGCalTriggerGeometryBase::fhTopology(), geometry_, edm::Event::getByToken(), HcalEndcap, HcalHitRelabeller::relabel(), simhits_bh_token_, simhits_ee_token_, simhits_fh_token_, HGCalDDDConstants::simToReco(), and HGCalTestNumbering::unpackHexagonIndex().

Referenced by calibrate(), fill(), and ~HGCalTriggerNtupleHGCTriggerCells().

 {
   edm::Handle<edm::PCaloHitContainer> ee_simhits_h;
   e.getByToken(simhits_ee_token_,ee_simhits_h);
   const edm::PCaloHitContainer& ee_simhits = *ee_simhits_h;
   edm::Handle<edm::PCaloHitContainer> fh_simhits_h;
   e.getByToken(simhits_fh_token_,fh_simhits_h);
   const edm::PCaloHitContainer& fh_simhits = *fh_simhits_h;
   edm::Handle<edm::PCaloHitContainer> bh_simhits_h;
   e.getByToken(simhits_bh_token_,bh_simhits_h);
   const edm::PCaloHitContainer& bh_simhits = *bh_simhits_h;
 
   //EE
   int layer=0,cell=0, sec=0, subsec=0, zp=0,subdet=0;
   ForwardSubdetector mysubdet;
   for( const auto& simhit : ee_simhits )
   { 
     HGCalTestNumbering::unpackHexagonIndex(simhit.id(), subdet, zp, layer, sec, subsec, cell); 
     mysubdet = (ForwardSubdetector)subdet;
     std::pair<int,int> recoLayerCell = geometry_->eeTopology().dddConstants().simToReco(cell,layer,sec,geometry_->eeTopology().detectorType());
     cell  = recoLayerCell.first;
     layer = recoLayerCell.second;
     if (layer<0 || cell<0) continue;
     auto itr_insert = simhits_ee.emplace(HGCalDetId(mysubdet,zp,layer,subsec,sec,cell), 0.);
     itr_insert.first->second += simhit.energy();
   }
 
   //  FH
   layer=0; cell=0; sec=0; subsec=0; zp=0; subdet=0;
   for( const auto& simhit : fh_simhits ) 
   { 
     HGCalTestNumbering::unpackHexagonIndex(simhit.id(), subdet, zp, layer, sec, subsec, cell); 
     mysubdet = (ForwardSubdetector)(subdet);
     std::pair<int,int> recoLayerCell = geometry_->fhTopology().dddConstants().simToReco(cell,layer,sec,geometry_->fhTopology().detectorType());
     cell  = recoLayerCell.first;
     layer = recoLayerCell.second;
     if (layer<0 || cell<0) continue;
     auto itr_insert = simhits_fh.emplace(HGCalDetId(mysubdet,zp,layer,subsec,sec,cell), 0.);
     itr_insert.first->second += simhit.energy();
   }      
 
   //  BH
   for( const auto& simhit : bh_simhits ) 
   { 
     HcalDetId id = HcalHitRelabeller::relabel(simhit.id(), geometry_->bhTopology().dddConstants());
     if (id.subdetId()!=HcalEndcap) continue;
     auto itr_insert = simhits_bh.emplace(id, 0.);
     itr_insert.first->second += simhit.energy();
   }      
 }