Inheritance diagram for HGCalHitCalibration:

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

	~HGCalHitCalibration ()

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

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

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

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::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

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

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

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

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
virtual void	analyze (const edm::Event &, const edm::EventSetup &) override

virtual void	beginJob () override

virtual void	endJob () override

Private Attributes
edm::EDGetTokenT< std::vector< CaloParticle > >	_caloParticles

edm::EDGetTokenT< HGCRecHitCollection >	_recHitsBH

edm::EDGetTokenT< HGCRecHitCollection >	_recHitsEE

edm::EDGetTokenT< HGCRecHitCollection >	_recHitsFH

int	algo_

bool	cutOnPt_

double	cutValue_

std::vector< float >	Energy_layer_calib

std::vector< float >	Energy_layer_calib_fraction

TH1F *	h_EoP_CPene_100_calib_fraction_ [6]

TH1F *	h_EoP_CPene_200_calib_fraction_ [6]

TH1F *	h_EoP_CPene_300_calib_fraction_ [6]

TH1F *	h_LayerOccupancy_ [6]

TH1F *	h_MissingHit_ [3]

bool	rawRecHits_

hgcal::RecHitTools	recHitTools

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)

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

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

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

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

void	consumesMany (const TypeToGet &id)

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

template<typename 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)

Detailed Description

Definition at line 44 of file HGCalHitCalibration.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 76 of file HGCalHitCalibration.cc.

References _caloParticles, _recHitsBH, _recHitsEE, _recHitsFH, algo_, gamEcalExtractorBlocks_cff::detector, edm::ParameterSet::getParameter(), h_EoP_CPene_100_calib_fraction_, h_EoP_CPene_200_calib_fraction_, h_EoP_CPene_300_calib_fraction_, h_LayerOccupancy_, h_MissingHit_, gen::k, TFileService::kSharedResource, dataset::name, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                        :
   rawRecHits_(iConfig.getParameter<bool>("rawRecHits")),
   cutOnPt_(iConfig.getParameter<bool>("cutOnPt")),
   cutValue_(iConfig.getParameter<double>("cutValue")) {
 
   usesResource(TFileService::kSharedResource);
   std::string detector = iConfig.getParameter<std::string >("detector");
   _recHitsEE = consumes<HGCRecHitCollection>(edm::InputTag("HGCalRecHit","HGCEERecHits"));
   _recHitsFH = consumes<HGCRecHitCollection>(edm::InputTag("HGCalRecHit","HGCHEFRecHits"));
   _recHitsBH = consumes<HGCRecHitCollection>(edm::InputTag("HGCalRecHit","HGCHEBRecHits"));
   if (detector=="all") {
     algo_ = 1;
   } else if (detector=="EM") {
     algo_ = 2;
   } else if (detector=="HAD") {
     algo_ = 3;
   }
   _caloParticles = consumes<std::vector<CaloParticle> >(edm::InputTag("mix","MergedCaloTruth"));
 
   edm::Service<TFileService> fs;
   std::string particle[6] = {"All", "Electron", "Muon", "Photon", "ChgHad", "NeutHad"};
   char name[100];
   for (int k=0; k<6; ++k) {
     sprintf(name,"h_EoP_CPene_100_calib_fraction_%s",particle[k].c_str());
     h_EoP_CPene_100_calib_fraction_[k] = fs->make<TH1F>(name,"",1000,-0.5,2.5);
     sprintf(name,"h_EoP_CPene_200_calib_fraction_%s",particle[k].c_str());
     h_EoP_CPene_200_calib_fraction_[k] = fs->make<TH1F>(name,"",1000,-0.5,2.5);
     sprintf(name,"h_EoP_CPene_300_calib_fraction_%s",particle[k].c_str());
     h_EoP_CPene_300_calib_fraction_[k] = fs->make<TH1F>(name,"",1000,-0.5,2.5);
     sprintf(name,"h_LayerOccupancy_%s",particle[k].c_str());
     h_LayerOccupancy_[k] = fs->make<TH1F>(name, "", 60, 0., 60.);
   }
   std::string dets[3] = {"EE", "FH", "BH"};
   for (int k=0; k<3; ++k) {
     sprintf(name,"h_missHit_%s",dets[k].c_str());
     h_MissingHit_[k] = fs->make<TH1F>(name, "", 200, 0., 200.);
   }
 }

HGCalHitCalibration::~HGCalHitCalibration ( )

Definition at line 115 of file HGCalHitCalibration.cc.

115 { }

Member Function Documentation

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

overrideprivatevirtual

Definition at line 117 of file HGCalHitCalibration.cc.

References _caloParticles, _recHitsBH, _recHitsEE, _recHitsFH, funct::abs(), algo_, edm::SortedCollection< T, SORT >::begin(), caloTruthProducer_cfi::caloParticles, gather_cfg::cout, TkAlMuonSelectors_cfi::cut, cutOnPt_, cutValue_, DetId::det(), edm::SortedCollection< T, SORT >::end(), Energy_layer_calib, Energy_layer_calib_fraction, HcalObjRepresent::Fill(), DetId::Forward, runEdmFileComparison::found, edm::Event::getByToken(), hgcal::RecHitTools::getEventSetup(), hgcal::RecHitTools::getLayerWithOffset(), hgcal::RecHitTools::getSiThickness(), h_EoP_CPene_100_calib_fraction_, h_EoP_CPene_200_calib_fraction_, h_EoP_CPene_300_calib_fraction_, h_LayerOccupancy_, h_MissingHit_, DetId::Hcal, HcalEndcap, HGCEE, HGCHEF, gen::k, cosmictrackingParticleSelector_cfi::pdgId, rawRecHits_, recHitTools, and DetId::subdetId().

Referenced by beginJob().

                                                                                      {
 
   recHitTools.getEventSetup(iSetup);
 
   edm::Handle<HGCRecHitCollection> recHitHandleEE;
   edm::Handle<HGCRecHitCollection> recHitHandleFH;
   edm::Handle<HGCRecHitCollection> recHitHandleBH;
   iEvent.getByToken(_recHitsEE,recHitHandleEE);
   iEvent.getByToken(_recHitsFH,recHitHandleFH);
   iEvent.getByToken(_recHitsBH,recHitHandleBH);
   const HGCRecHitCollection& rechitsEE = *recHitHandleEE;
   const HGCRecHitCollection& rechitsFH = *recHitHandleFH;
   const HGCRecHitCollection& rechitsBH = *recHitHandleBH;
   
   edm::Handle<std::vector<CaloParticle> > caloParticleHandle;
   iEvent.getByToken(_caloParticles, caloParticleHandle);
   const std::vector<CaloParticle>& caloParticles = *caloParticleHandle;
   const int pdgId[3] = {11, 13, 22};
 
   // loop over caloParticles
   int mhit[3] = {0,0,0};
   for (auto it_caloPart = caloParticles.begin();
        it_caloPart != caloParticles.end(); ++it_caloPart) {
     double cut = (cutOnPt_) ? it_caloPart->pt() : it_caloPart->energy();
     if (cut > cutValue_) {
       int type(5);
       if      (std::abs(it_caloPart->pdgId()) == pdgId[0]) type = 1;
       else if (std::abs(it_caloPart->pdgId()) == pdgId[1]) type = 2;
       else if (it_caloPart->pdgId()           == pdgId[2]) type = 3;
       else if (it_caloPart->threeCharge()     != 0)        type = 4;
       const SimClusterRefVector simClusterRefVector = it_caloPart->simClusters();
 
       //size should be HGC layers 52 is enough
       Energy_layer_calib.assign(60,0.);
       Energy_layer_calib_fraction.assign(60,0.);
     
       int seedDet = 0;
       float seedEnergy = 0.;
       int simClusterCount = 0;
     
       for ( const auto & simCluster : simClusterRefVector) {
     ++simClusterCount;
 #ifdef EDM_ML_DEBUG
     std::cout << ">>> simCluster.energy() = " << simCluster->energy() << std::endl;
 #endif
     const std::vector<std::pair<uint32_t,float> > hits_and_fractions = simCluster->hits_and_fractions();
 
     //loop over hits      
     for (auto it_haf = hits_and_fractions.begin(); 
          it_haf != hits_and_fractions.end(); ++it_haf) {
       unsigned int hitlayer = recHitTools.getLayerWithOffset(it_haf->first);
       DetId hitid = (it_haf->first); 
 
       // dump raw RecHits and match
       if (rawRecHits_) {
         if (hitid.det() == DetId::Forward && hitid.subdetId() == HGCEE &&
         (algo_ == 1 || algo_ == 2)) {
           // loop over EE RecHits
           bool found(false);
           for (auto it_hit = rechitsEE.begin(); 
            it_hit < rechitsEE.end(); ++it_hit) {
         const HGCalDetId detid = it_hit->detid();
         unsigned int layer = recHitTools.getLayerWithOffset(detid);
         if (detid == hitid) {
           found = true;
           if (hitlayer != layer) {
 #ifdef EDM_ML_DEBUG
             std::cout << " recHit ID problem EE " << std::endl;
 #endif
             return;
           }
           Energy_layer_calib_fraction[layer] += it_hit->energy()*it_haf->second;
           h_LayerOccupancy_[0]->Fill(layer);
           h_LayerOccupancy_[type]->Fill(layer);
           if(seedEnergy < it_hit->energy()){
             seedEnergy = it_hit->energy();
             seedDet = recHitTools.getSiThickness(detid);
           }
           break;
         }
           }
           if (!found) ++mhit[0];
         }
         if (hitid.det() == DetId::Forward && hitid.subdetId() == HGCHEF &&
         (algo_ == 1 || algo_ == 3)) {
           // loop over HEF RecHits
           bool found(false);
           for (auto it_hit = rechitsFH.begin(); 
            it_hit < rechitsFH.end(); ++it_hit) {
         const HGCalDetId detid = it_hit->detid();
         unsigned int layer = recHitTools.getLayerWithOffset(detid);
         if (detid == hitid) {
           found = true;
           if (hitlayer != layer) {
 #ifdef EDM_ML_DEBUG
             std::cout << " recHit ID problem FH " << std::endl; 
 #endif
             return;
           }
           Energy_layer_calib_fraction[layer] += it_hit->energy()*it_haf->second;
           h_LayerOccupancy_[0]->Fill(layer);
           h_LayerOccupancy_[type]->Fill(layer);
           if(seedEnergy < it_hit->energy()){
             seedEnergy = it_hit->energy();
             seedDet = recHitTools.getSiThickness(detid);
           }
           break;
         }
           }
           if (!found) ++mhit[1];
         }
         if (hitid.det() == DetId::Hcal && hitid.subdetId() == HcalEndcap &&
         (algo_ == 1 || algo_ == 3)) {
           // loop over BH RecHits
           bool found(false);
           for (auto it_hit = rechitsBH.begin();
            it_hit < rechitsBH.end(); ++it_hit) {
         const HcalDetId detid = it_hit->detid();
         unsigned int layer = recHitTools.getLayerWithOffset(detid);
         if (detid == hitid) {
           found = true;
           if (hitlayer != layer) {
 #ifdef EDM_ML_DEBUG
             std::cout << " recHit ID problem BH " << std::endl; 
 #endif
             return;
           }
           Energy_layer_calib_fraction[layer] += it_hit->energy()*it_haf->second;
           h_LayerOccupancy_[0]->Fill(layer);
           h_LayerOccupancy_[type]->Fill(layer);
           if (seedEnergy < it_hit->energy()){
             seedEnergy = it_hit->energy();
             seedDet = 300; // recHitTools.getSiThickness(detid);
           }
           break;
         }
           }
           if (!found) ++mhit[2];
         }
       
       }//end recHits
     }// end simHit
       }//end simCluster
 
       float sumCalibRecHitCalib_fraction = 0;
       for(unsigned int iL=0; iL<Energy_layer_calib_fraction.size(); ++iL){
     sumCalibRecHitCalib_fraction += Energy_layer_calib_fraction[iL];
       }
     
       double ebyp = sumCalibRecHitCalib_fraction / it_caloPart->energy();
       if(seedDet == 100) {
     h_EoP_CPene_100_calib_fraction_[0]->Fill(ebyp);
     h_EoP_CPene_100_calib_fraction_[type]->Fill(ebyp);
       } else if (seedDet == 200) {
     h_EoP_CPene_200_calib_fraction_[0]->Fill(ebyp);
     h_EoP_CPene_200_calib_fraction_[type]->Fill(ebyp);
       } else if (seedDet == 300){
     h_EoP_CPene_300_calib_fraction_[0]->Fill(ebyp);
     h_EoP_CPene_300_calib_fraction_[type]->Fill(ebyp);
       }
     }
   }//end caloparticle
 
   for (int k=0; k<3; ++k) h_MissingHit_[k]->Fill(mhit[k]);
 }

virtual void HGCalHitCalibration::beginJob ( void )

inlineoverrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 52 of file HGCalHitCalibration.cc.

References analyze().

52 {}

virtual void HGCalHitCalibration::endJob ( void )

inlineoverrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 54 of file HGCalHitCalibration.cc.

54 {}

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

static

Definition at line 285 of file HGCalHitCalibration.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), DEFINE_FWK_MODULE, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                 {
   //The following says we do not know what parameters are allowed so do no validation
   // Please change this to state exactly what you do use, even if it is no parameters
   edm::ParameterSetDescription desc;
   desc.add<std::string>("detector","all");
   desc.add<bool>("rawRecHits",true);
   desc.add<bool>("cutOnPt",true);
   desc.add<double>("cutValue",10.0);
   descriptions.add("hgcalHitCalibration",desc);
 }