SimplePi0DiscAnalyzer Class Reference

#include <SimplePi0DiscAnalyzer.h>

Inheritance diagram for SimplePi0DiscAnalyzer:

Public Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override
void	beginJob () override
void	endJob () override
	SimplePi0DiscAnalyzer (const edm::ParameterSet &conf)
	~SimplePi0DiscAnalyzer () override
Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
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 &&)=default
	EDConsumerBase (EDConsumerBase const &)=delete
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > 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::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase &	operator= (EDConsumerBase &&)=default
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Attributes
TH1F *	hAll_nnout_Assoc_
TH1F *	hAll_nnout_NoConv_Assoc_
TH1F *	hAll_nnout_NoConv_Assoc_R9_
TH1F *	hBarrel_nnout_Assoc_
TH1F *	hBarrel_nnout_NoConv_Assoc_
TH1F *	hBarrel_nnout_NoConv_Assoc_R9_
TH1F *	hConv_ntracks_
TH1F *	hEndcNoPresh_nnout_Assoc_
TH1F *	hEndcNoPresh_nnout_NoConv_Assoc_
TH1F *	hEndcNoPresh_nnout_NoConv_Assoc_R9_
TH1F *	hEndcWithPresh_nnout_Assoc_
TH1F *	hEndcWithPresh_nnout_NoConv_Assoc_
TH1F *	hEndcWithPresh_nnout_NoConv_Assoc_R9_
std::string	outputFile_
std::string	photonCollection_
std::string	photonCollectionProducer_
TFile *	rootFile_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Description: Pi0Disc analyzer using reco data

Implementation: <Notes on="" implementation>="">

Definition at line 40 of file SimplePi0DiscAnalyzer.h.

Constructor & Destructor Documentation

SimplePi0DiscAnalyzer()

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

explicit

Definition at line 49 of file SimplePi0DiscAnalyzer.cc.

                                                                        {
  outputFile_ = conf.getParameter<std::string>("outputFile");
  rootFile_ = new TFile(outputFile_.c_str(), "RECREATE");
 
  photonCollectionProducer_ = conf.getParameter<std::string>("phoProducer");
  photonCollection_ = conf.getParameter<std::string>("photonCollection");
}

References edm::ParameterSet::getParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

~SimplePi0DiscAnalyzer()

SimplePi0DiscAnalyzer::~SimplePi0DiscAnalyzer ( )

override

Definition at line 57 of file SimplePi0DiscAnalyzer.cc.

                                              {
  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)
  rootFile_->Write();
  rootFile_->Close();
}

Member Function Documentation

analyze()

void SimplePi0DiscAnalyzer::analyze ( const edm::Event &  e, const edm::EventSetup &  c  )

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 120 of file SimplePi0DiscAnalyzer.cc.

                                                                                       {
  std::cout << std::endl;
  std::cout << " -------------- NEW EVENT : Run, Event =  " << iEvent.id() << std::endl;
 
  edm::Handle<reco::PhotonCollection> PhotonHandle;
  iEvent.getByLabel(photonCollectionProducer_, photonCollection_, PhotonHandle);
  const reco::PhotonCollection photons = *(PhotonHandle.product());
 
  std::cout << "----> Photons size: " << photons.size() << std::endl;
 
  edm::Handle<reco::PhotonPi0DiscriminatorAssociationMap> map;
  iEvent.getByLabel("piZeroDiscriminators", "PhotonPi0DiscriminatorAssociationMap", map);
  reco::PhotonPi0DiscriminatorAssociationMap::const_iterator mapIter;
 
  //  int PhoInd = 0;
 
  for (reco::PhotonCollection::const_iterator iPho = photons.begin(); iPho != photons.end();
       iPho++) {  // Loop over Photons
 
    reco::Photon localPho(*iPho);
 
    float Photon_et = localPho.et();
    float Photon_eta = localPho.eta();
    float Photon_phi = localPho.phi();
    float Photon_r9 = localPho.r9();
    bool isPhotConv = localPho.hasConversionTracks();
    //    std::cout << "Photon Id = " << PhoInd
    std::cout << "Photon Id = " << iPho - photons.begin() << " with Et = " << Photon_et << " Eta = " << Photon_eta
              << " Phi = " << Photon_phi << " R9 = " << Photon_r9 << " and conv_id = " << isPhotConv << std::endl;
 
    SuperClusterRef it_super = localPho.superCluster();  // get the SC related to the Photon candidate
 
    //    hConv_ntracks_->Fill(Ntrk_conv);
 
    float nn = -10;
    //    mapIter = map->find(edm::Ref<reco::PhotonCollection>(PhotonHandle,PhoInd));
    mapIter = map->find(edm::Ref<reco::PhotonCollection>(PhotonHandle, iPho - photons.begin()));
    if (mapIter != map->end()) {
      nn = mapIter->val;
    }
    if (fabs(it_super->eta()) <= 1.442) {
      hBarrel_nnout_Assoc_->Fill(nn);
      hAll_nnout_Assoc_->Fill(nn);
      std::cout << "AssociationMap Barrel NN = " << nn << std::endl;
      if (!isPhotConv) {
        hBarrel_nnout_NoConv_Assoc_->Fill(nn);
        hAll_nnout_NoConv_Assoc_->Fill(nn);
      }
      if (Photon_r9 > 0.93) {
        hBarrel_nnout_NoConv_Assoc_R9_->Fill(nn);
        hAll_nnout_NoConv_Assoc_R9_->Fill(nn);
      }
    } else if ((fabs(it_super->eta()) >= 1.556 && fabs(it_super->eta()) < 1.65) || fabs(it_super->eta()) > 2.5) {
      hEndcNoPresh_nnout_Assoc_->Fill(nn);
      hAll_nnout_Assoc_->Fill(nn);
      std::cout << "AssociationMap EndcNoPresh NN = " << nn << std::endl;
      if (!isPhotConv) {
        hEndcNoPresh_nnout_NoConv_Assoc_->Fill(nn);
        hAll_nnout_NoConv_Assoc_->Fill(nn);
      }
      if (Photon_r9 > 0.93) {
        hEndcNoPresh_nnout_NoConv_Assoc_R9_->Fill(nn);
        hAll_nnout_NoConv_Assoc_R9_->Fill(nn);
      }
    } else if (fabs(it_super->eta()) >= 1.65 && fabs(it_super->eta()) <= 2.5) {
      hEndcWithPresh_nnout_Assoc_->Fill(nn);
      hAll_nnout_Assoc_->Fill(nn);
      std::cout << "AssociationMap EndcWithPresh NN = " << nn << std::endl;
      if (!isPhotConv) {
        hEndcWithPresh_nnout_NoConv_Assoc_->Fill(nn);
        hAll_nnout_NoConv_Assoc_->Fill(nn);
      }
      if (Photon_r9 > 0.93) {
        hEndcWithPresh_nnout_NoConv_Assoc_R9_->Fill(nn);
        hAll_nnout_NoConv_Assoc_R9_->Fill(nn);
      }
    }
 
    //    PhoInd++;
  }  // End Loop over Photons
}

References gather_cfg::cout, reco::LeafCandidate::et(), reco::LeafCandidate::eta(), reco::Photon::hasConversionTracks(), iEvent, genParticles_cff::map, groupFilesInBlocks::nn, reco::LeafCandidate::phi(), BPHMonitor_cfi::photons, edm::Handle< T >::product(), reco::Photon::r9(), reco::Photon::superCluster(), and edm::helpers::KeyVal< K, V >::val.

beginJob()

void SimplePi0DiscAnalyzer::beginJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 64 of file SimplePi0DiscAnalyzer.cc.

                                     {
  rootFile_->cd();
  std::cout << "beginJob() ->  Book the Histograms" << std::endl;
 
  hConv_ntracks_ = new TH1F("nConvTracks", "Number of tracks of converted Photons ", 10, 0., 10);
  hAll_nnout_Assoc_ = new TH1F("All_nnout_Assoc", "NNout for All Photons(AssociationMap)", 100, 0., 1.);
  hAll_nnout_NoConv_Assoc_ =
      new TH1F("All_nnout_NoConv_Assoc", "NNout for Unconverted Photons(AssociationMap)", 100, 0., 1.);
  hAll_nnout_NoConv_Assoc_R9_ =
      new TH1F("All_nnout_NoConv_Assoc_R9", "NNout for Unconverted Photons with R9>0.93 (AssociationMap)", 100, 0., 1.);
  hBarrel_nnout_Assoc_ = new TH1F("barrel_nnout_Assoc", "NNout for Barrel Photons(AssociationMap)", 100, 0., 1.);
  hBarrel_nnout_NoConv_Assoc_ =
      new TH1F("barrel_nnout_NoConv_Assoc", "NNout for Barrel Unconverted Photons(AssociationMap)", 100, 0., 1.);
  hBarrel_nnout_NoConv_Assoc_R9_ = new TH1F(
      "barrel_nnout_NoConv_Assoc_R9", "NNout for Barrel Unconverted Photons with R9>0.93 (AssociationMap)", 100, 0., 1.);
  hEndcNoPresh_nnout_Assoc_ =
      new TH1F("endcNoPresh_nnout_Assoc", "NNout for Endcap NoPresh Photons(AssociationMap)", 100, 0., 1.);
  hEndcNoPresh_nnout_NoConv_Assoc_ = new TH1F(
      "endcNoPresh_nnout_NoConv_Assoc", "NNout for Endcap Unconverted NoPresh Photons(AssociationMap)", 100, 0., 1.);
  hEndcNoPresh_nnout_NoConv_Assoc_R9_ =
      new TH1F("endcNoPresh_nnout_NoConv_Assoc_R9",
               "NNout for Endcap Unconverted NoPresh Photons with R9>0.93 (AssociationMap)",
               100,
               0.,
               1.);
  hEndcWithPresh_nnout_Assoc_ =
      new TH1F("endcWithPresh_nnout_Assoc", "NNout for Endcap WithPresh Photons(AssociationMap)", 100, 0., 1.);
  hEndcWithPresh_nnout_NoConv_Assoc_ = new TH1F(
      "endcWithPresh_nnout_NoConv_Assoc", "NNout for Endcap Unconverted WithPresh Photons(AssociationMap)", 100, 0., 1.);
  hEndcWithPresh_nnout_NoConv_Assoc_R9_ =
      new TH1F("endcWithPresh_nnout_NoConv_Assoc_R9",
               "NNout for Endcap Unconverted WithPresh Photons with R9>0.93 (AssociationMap)",
               100,
               0.,
               1.);
}

References gather_cfg::cout.

endJob()

void SimplePi0DiscAnalyzer::endJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 101 of file SimplePi0DiscAnalyzer.cc.

                                   {
  rootFile_->cd();
  std::cout << "endJob() ->  Write the Histograms" << std::endl;
  hConv_ntracks_->Write();
 
  hAll_nnout_Assoc_->Write();
  hAll_nnout_NoConv_Assoc_->Write();
  hAll_nnout_NoConv_Assoc_R9_->Write();
  hBarrel_nnout_Assoc_->Write();
  hBarrel_nnout_NoConv_Assoc_->Write();
  hBarrel_nnout_NoConv_Assoc_R9_->Write();
  hEndcNoPresh_nnout_Assoc_->Write();
  hEndcNoPresh_nnout_NoConv_Assoc_->Write();
  hEndcNoPresh_nnout_NoConv_Assoc_R9_->Write();
  hEndcWithPresh_nnout_Assoc_->Write();
  hEndcWithPresh_nnout_NoConv_Assoc_->Write();
  hEndcWithPresh_nnout_NoConv_Assoc_R9_->Write();
}

References gather_cfg::cout.

Member Data Documentation

hAll_nnout_Assoc_

TH1F* SimplePi0DiscAnalyzer::hAll_nnout_Assoc_

private

Definition at line 61 of file SimplePi0DiscAnalyzer.h.

hAll_nnout_NoConv_Assoc_

TH1F* SimplePi0DiscAnalyzer::hAll_nnout_NoConv_Assoc_

private

Definition at line 62 of file SimplePi0DiscAnalyzer.h.

hAll_nnout_NoConv_Assoc_R9_

TH1F* SimplePi0DiscAnalyzer::hAll_nnout_NoConv_Assoc_R9_

private

Definition at line 69 of file SimplePi0DiscAnalyzer.h.

hBarrel_nnout_Assoc_

TH1F* SimplePi0DiscAnalyzer::hBarrel_nnout_Assoc_

private

Definition at line 63 of file SimplePi0DiscAnalyzer.h.

hBarrel_nnout_NoConv_Assoc_

TH1F* SimplePi0DiscAnalyzer::hBarrel_nnout_NoConv_Assoc_

private

Definition at line 64 of file SimplePi0DiscAnalyzer.h.

hBarrel_nnout_NoConv_Assoc_R9_

TH1F* SimplePi0DiscAnalyzer::hBarrel_nnout_NoConv_Assoc_R9_

private

Definition at line 70 of file SimplePi0DiscAnalyzer.h.

hConv_ntracks_

TH1F* SimplePi0DiscAnalyzer::hConv_ntracks_

private

Definition at line 59 of file SimplePi0DiscAnalyzer.h.

hEndcNoPresh_nnout_Assoc_

TH1F* SimplePi0DiscAnalyzer::hEndcNoPresh_nnout_Assoc_

private

Definition at line 65 of file SimplePi0DiscAnalyzer.h.

hEndcNoPresh_nnout_NoConv_Assoc_

TH1F* SimplePi0DiscAnalyzer::hEndcNoPresh_nnout_NoConv_Assoc_

private

Definition at line 66 of file SimplePi0DiscAnalyzer.h.

hEndcNoPresh_nnout_NoConv_Assoc_R9_

TH1F* SimplePi0DiscAnalyzer::hEndcNoPresh_nnout_NoConv_Assoc_R9_

private

Definition at line 71 of file SimplePi0DiscAnalyzer.h.

hEndcWithPresh_nnout_Assoc_

TH1F* SimplePi0DiscAnalyzer::hEndcWithPresh_nnout_Assoc_

private

Definition at line 67 of file SimplePi0DiscAnalyzer.h.

hEndcWithPresh_nnout_NoConv_Assoc_

TH1F* SimplePi0DiscAnalyzer::hEndcWithPresh_nnout_NoConv_Assoc_

private

Definition at line 68 of file SimplePi0DiscAnalyzer.h.

hEndcWithPresh_nnout_NoConv_Assoc_R9_

TH1F* SimplePi0DiscAnalyzer::hEndcWithPresh_nnout_NoConv_Assoc_R9_

private

Definition at line 72 of file SimplePi0DiscAnalyzer.h.

outputFile_

std::string SimplePi0DiscAnalyzer::outputFile_

private

Definition at line 56 of file SimplePi0DiscAnalyzer.h.

photonCollection_

std::string SimplePi0DiscAnalyzer::photonCollection_

private

Definition at line 54 of file SimplePi0DiscAnalyzer.h.

photonCollectionProducer_

std::string SimplePi0DiscAnalyzer::photonCollectionProducer_

private

Definition at line 53 of file SimplePi0DiscAnalyzer.h.

rootFile_

TFile* SimplePi0DiscAnalyzer::rootFile_

private

Definition at line 57 of file SimplePi0DiscAnalyzer.h.