Inheritance diagram for HGCalTriggerNtupleHGCDigis:

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

	HGCalTriggerNtupleHGCDigis (const edm::ParameterSet &conf)

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

	~HGCalTriggerNtupleHGCDigis ()

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

virtual	~HGCalTriggerNtupleBase ()

Private Member Functions
virtual void	clear () override final

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

Private Attributes
edm::EDGetToken	ee_token_

edm::EDGetToken	fh_token_

edm::ESHandle< HGCalGeometry >	geom_ee

edm::ESHandle< HGCalGeometry >	geom_fh

std::vector< int >	hgcdigi_cell_

std::vector< uint32_t >	hgcdigi_data_

std::vector< float >	hgcdigi_eta_

std::vector< int >	hgcdigi_id_

std::vector< int >	hgcdigi_isadc_

std::vector< int >	hgcdigi_layer_

int	hgcdigi_n_

std::vector< float >	hgcdigi_phi_

std::vector< int >	hgcdigi_side_

std::vector< float >	hgcdigi_simenergy_

std::vector< int >	hgcdigi_subdet_

std::vector< int >	hgcdigi_wafer_

std::vector< int >	hgcdigi_wafertype_

std::vector< float >	hgcdigi_z_

bool	is_Simhit_comp_

edm::EDGetToken	SimHits_inputee_

edm::EDGetToken	SimHits_inputfh_

edm::ESHandle< HGCalTopology >	topo_ee

edm::ESHandle< HGCalTopology >	topo_fh

Additional Inherited Members

Detailed Description

Definition at line 12 of file HGCalTriggerNtupleHGCDigis.cc.

Constructor & Destructor Documentation

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

Definition at line 55 of file HGCalTriggerNtupleHGCDigis.cc.

References edm::ParameterSet::getParameter(), initialize(), and is_Simhit_comp_.

                                                        :HGCalTriggerNtupleBase(conf)
 {
     is_Simhit_comp_ = conf.getParameter<bool>("isSimhitComp");
 
 }

HGCalTriggerNtupleHGCDigis::~HGCalTriggerNtupleHGCDigis ( )

inline

Definition at line 17 of file HGCalTriggerNtupleHGCDigis.cc.

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

17 {};

Member Function Documentation

void HGCalTriggerNtupleHGCDigis::clear ( void )

finaloverrideprivatevirtual

Implements HGCalTriggerNtupleBase.

Definition at line 227 of file HGCalTriggerNtupleHGCDigis.cc.

References hgcdigi_cell_, hgcdigi_data_, hgcdigi_eta_, hgcdigi_id_, hgcdigi_isadc_, hgcdigi_layer_, hgcdigi_n_, hgcdigi_phi_, hgcdigi_side_, hgcdigi_simenergy_, hgcdigi_subdet_, hgcdigi_wafer_, hgcdigi_wafertype_, hgcdigi_z_, and is_Simhit_comp_.

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

 {
     hgcdigi_n_ = 0;
     hgcdigi_id_.clear();
     hgcdigi_subdet_.clear();
     hgcdigi_side_.clear();
     hgcdigi_layer_.clear();
     hgcdigi_wafer_.clear();
     hgcdigi_wafertype_.clear();
     hgcdigi_cell_.clear();
     hgcdigi_eta_.clear();
     hgcdigi_phi_.clear();
     hgcdigi_z_.clear();
     hgcdigi_data_.clear();
     hgcdigi_isadc_.clear();
     if  (is_Simhit_comp_) hgcdigi_simenergy_.clear();
 }

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

finaloverridevirtual

Implements HGCalTriggerNtupleBase.

Definition at line 90 of file HGCalTriggerNtupleHGCDigis.cc.

References clear(), data, ee_token_, PV3DBase< T, PVType, FrameType >::eta(), fh_token_, geom_ee, geom_fh, edm::EventSetup::get(), edm::Event::getByToken(), HGCalGeometry::getPosition(), hgcdigi_cell_, hgcdigi_data_, hgcdigi_eta_, hgcdigi_id_, hgcdigi_isadc_, hgcdigi_layer_, hgcdigi_n_, hgcdigi_phi_, hgcdigi_side_, hgcdigi_simenergy_, hgcdigi_subdet_, hgcdigi_wafer_, hgcdigi_wafertype_, hgcdigi_z_, HGCEE, HGCHEF, hcalTTPDigis_cfi::id, is_Simhit_comp_, ALCARECOPromptCalibProdSiPixelAli0T_cff::mode, PV3DBase< T, PVType, FrameType >::phi(), simhits(), edm::SortedCollection< T, SORT >::size(), topo_ee, topo_fh, PV3DBase< T, PVType, FrameType >::z(), and HcalDigiParam_cfi::zside.

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

 {
  
     es.get<IdealGeometryRecord>().get("HGCalEESensitive", geom_ee);
     es.get<IdealGeometryRecord>().get("HGCalHESiliconSensitive", geom_fh);
     es.get<IdealGeometryRecord>().get("HGCalEESensitive",topo_ee);
     es.get<IdealGeometryRecord>().get("HGCalHESiliconSensitive",topo_fh);
     edm::Handle<HGCEEDigiCollection> ee_digis_h;
     e.getByToken(ee_token_, ee_digis_h);
     const HGCEEDigiCollection& ee_digis = *ee_digis_h;
     edm::Handle<HGCHEDigiCollection> fh_digis_h;
     e.getByToken(fh_token_, fh_digis_h);
     const HGCHEDigiCollection& fh_digis = *fh_digis_h;
 
     // sim hit association
     std::unordered_map<uint32_t, double> simhits_ee;
     std::unordered_map<uint32_t, double> simhits_fh;  
     if (is_Simhit_comp_) simhits(e, simhits_ee, simhits_fh);
     
     clear();
     hgcdigi_n_ = ee_digis.size() + fh_digis.size();
     hgcdigi_id_.reserve(hgcdigi_n_);
     hgcdigi_subdet_.reserve(hgcdigi_n_);
     hgcdigi_side_.reserve(hgcdigi_n_);
     hgcdigi_layer_.reserve(hgcdigi_n_);
     hgcdigi_wafer_.reserve(hgcdigi_n_);
     hgcdigi_wafertype_.reserve(hgcdigi_n_);
     hgcdigi_cell_.reserve(hgcdigi_n_);
     hgcdigi_eta_.reserve(hgcdigi_n_);
     hgcdigi_phi_.reserve(hgcdigi_n_);
     hgcdigi_z_.reserve(hgcdigi_n_);
     hgcdigi_data_.reserve(hgcdigi_n_);
     hgcdigi_isadc_.reserve(hgcdigi_n_);
     if (is_Simhit_comp_) hgcdigi_simenergy_.reserve(hgcdigi_n_);
     
     for(const auto& digi : ee_digis)
       {
         const HGCalDetId id(digi.id());
         hgcdigi_id_.emplace_back(id.rawId());
         hgcdigi_subdet_.emplace_back(ForwardSubdetector::HGCEE);
         hgcdigi_side_.emplace_back(id.zside());
         hgcdigi_layer_.emplace_back(id.layer());
         hgcdigi_wafer_.emplace_back(id.wafer());
         hgcdigi_wafertype_.emplace_back(id.waferType());
         hgcdigi_cell_.emplace_back(id.cell());
         GlobalPoint cellpos = geom_ee->getPosition(id.rawId());
         hgcdigi_eta_.emplace_back(cellpos.eta());
         hgcdigi_phi_.emplace_back(cellpos.phi());
         hgcdigi_z_.emplace_back(cellpos.z());
         hgcdigi_data_.emplace_back(digi[2].data()); 
         int is_adc=0;
         if (!(digi[2].mode())) is_adc =1;
         hgcdigi_isadc_.emplace_back(is_adc);
         if (is_Simhit_comp_) {
           double hit_energy=0;
           auto itr = simhits_ee.find(id);
           if(itr!=simhits_ee.end())hit_energy = itr->second;
           hgcdigi_simenergy_.emplace_back(hit_energy); 
         }
       }
     
     for(const auto& digi : fh_digis)
       {
         const HGCalDetId id(digi.id());
         hgcdigi_id_.emplace_back(id.rawId());
         hgcdigi_subdet_.emplace_back(ForwardSubdetector::HGCHEF);
         hgcdigi_side_.emplace_back(id.zside());
         hgcdigi_layer_.emplace_back(id.layer());
         hgcdigi_wafer_.emplace_back(id.wafer());
         hgcdigi_wafertype_.emplace_back(id.waferType());
         hgcdigi_cell_.emplace_back(id.cell());
         GlobalPoint cellpos = geom_fh->getPosition(id.rawId());
         hgcdigi_eta_.emplace_back(cellpos.eta());
         hgcdigi_phi_.emplace_back(cellpos.phi());
         hgcdigi_z_.emplace_back(cellpos.z());
         hgcdigi_data_.emplace_back(digi[2].data()); 
         int is_adc=0;
         if (!(digi[2].mode())) is_adc =1;
         hgcdigi_isadc_.emplace_back(is_adc);
         if (is_Simhit_comp_) {
           double hit_energy=0;
           auto itr = simhits_fh.find(id);
           if(itr!=simhits_fh.end())hit_energy = itr->second;
           hgcdigi_simenergy_.emplace_back(hit_energy); 
         }
       }
 }

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

finaloverridevirtual

Implements HGCalTriggerNtupleBase.

Definition at line 63 of file HGCalTriggerNtupleHGCDigis.cc.

References ee_token_, fh_token_, fill(), edm::ParameterSet::getParameter(), hgcdigi_cell_, hgcdigi_data_, hgcdigi_eta_, hgcdigi_id_, hgcdigi_isadc_, hgcdigi_layer_, hgcdigi_n_, hgcdigi_phi_, hgcdigi_side_, hgcdigi_simenergy_, hgcdigi_subdet_, hgcdigi_wafer_, hgcdigi_wafertype_, hgcdigi_z_, is_Simhit_comp_, SimHits_inputee_, and SimHits_inputfh_.

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

 {
 
     ee_token_ = collector.consumes<HGCEEDigiCollection>(conf.getParameter<edm::InputTag>("HGCDigisEE")); 
     fh_token_ = collector.consumes<HGCHEDigiCollection>(conf.getParameter<edm::InputTag>("HGCDigisFH"));
     if (is_Simhit_comp_) {
       SimHits_inputee_ = collector.consumes<edm::PCaloHitContainer>(conf.getParameter<edm::InputTag>("eeSimHits"));
       SimHits_inputfh_ = collector.consumes<edm::PCaloHitContainer>(conf.getParameter<edm::InputTag>("fhSimHits"));
     }
     tree.Branch("hgcdigi_n", &hgcdigi_n_, "hgcdigi_n/I");
     tree.Branch("hgcdigi_id", &hgcdigi_id_);
     tree.Branch("hgcdigi_subdet", &hgcdigi_subdet_);
     tree.Branch("hgcdigi_zside", &hgcdigi_side_);
     tree.Branch("hgcdigi_layer", &hgcdigi_layer_);
     tree.Branch("hgcdigi_wafer", &hgcdigi_wafer_);
     tree.Branch("hgcdigi_wafertype", &hgcdigi_wafertype_);
     tree.Branch("hgcdigi_cell", &hgcdigi_cell_);    
     tree.Branch("hgcdigi_eta", &hgcdigi_eta_);    
     tree.Branch("hgcdigi_phi", &hgcdigi_phi_);    
     tree.Branch("hgcdigi_z", &hgcdigi_z_);    
     tree.Branch("hgcdigi_data", &hgcdigi_data_);
     tree.Branch("hgcdigi_isadc", &hgcdigi_isadc_);
     if (is_Simhit_comp_) tree.Branch("hgcdigi_simenergy", &hgcdigi_simenergy_);
 }

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

private

Definition at line 180 of file HGCalTriggerNtupleHGCDigis.cc.

References clear(), HGCalTopology::dddConstants(), HGCalTopology::detectorType(), edm::Event::getByToken(), SimHits_inputee_, SimHits_inputfh_, HGCalDDDConstants::simToReco(), topo_ee, topo_fh, and HGCalTestNumbering::unpackHexagonIndex().

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

 {
 
       edm::Handle<edm::PCaloHitContainer> ee_simhits_h;
       e.getByToken(SimHits_inputee_,ee_simhits_h);
       const edm::PCaloHitContainer& ee_simhits = *ee_simhits_h;
       edm::Handle<edm::PCaloHitContainer> fh_simhits_h;
       e.getByToken(SimHits_inputfh_,fh_simhits_h);
       const edm::PCaloHitContainer& fh_simhits = *fh_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 = topo_ee->dddConstants().simToReco(cell,layer,sec,topo_ee->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 = topo_fh->dddConstants().simToReco(cell,layer,sec,topo_fh->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();
       }      
 }