Inheritance diagram for TreeWriterForEcalCorrection:

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

	~TreeWriterForEcalCorrection () override

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

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () 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 const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (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 const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

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)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

Private Attributes
edm::Service< TFileService >	file

TTree *	tree

float	tree_e

float	tree_eta

float	tree_response

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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 23 of file TreeWriterForEcalCorrection.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 38 of file TreeWriterForEcalCorrection.cc.

References file, TFileService::make(), tree_e, tree_eta, and tree_response.

                                                                                        {
   tree = file->make<TTree>("responseTree", "same info as 3dhisto");
   tree->Branch("e", &tree_e, "e/F");
   tree->Branch("eta", &tree_eta, "eta/F");
   tree->Branch("r", &tree_response, "r/F");
 }

TreeWriterForEcalCorrection::~TreeWriterForEcalCorrection ( )

inlineoverride

Definition at line 26 of file TreeWriterForEcalCorrection.cc.

References analyze(), and fillDescriptions().

26 {};

Member Function Documentation

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

overrideprivate

Definition at line 45 of file TreeWriterForEcalCorrection.cc.

References reco::LeafCandidate::energy(), reco::LeafCandidate::eta(), Exception, heavyFlavorValidation_cfi::GenParticles, edm::Event::getByLabel(), tree_e, tree_eta, and tree_response.

Referenced by ~TreeWriterForEcalCorrection().

                                                                                              {
   // get generated particles
   edm::Handle<reco::GenParticleCollection> GenParticles;
   iEvent.getByLabel("genParticles", "", GenParticles);
 
   // As this module is intended for single particle guns, there should be
   // exactly one generated particle.
   if (GenParticles->size() != 1) {
     throw cms::Exception("MismatchedInputFiles") << "Intended for particle guns only\n";
   }
 
   // I assume here that the tracker simulation is disabled and no vertex
   // smearing is done, so that the generated particles is exactly the same
   // particle as the particle hitting the entrace of the ECAL.
   reco::GenParticle gen = GenParticles->at(0);
   float genE = gen.energy();
   float genEta = gen.eta();
   // genEta is positive for my photongun sample per definition.
   // If not, you could take the absolute value here.
 
   // get sim hits
   edm::Handle<edm::PCaloHitContainer> SimHitsEB;
   edm::Handle<edm::PCaloHitContainer> SimHitsEE;
   edm::Handle<edm::PCaloHitContainer> SimHitsES;
 
   // Finds out automatically, if this is fullsim or fastsim
   bool isFastSim = iEvent.getByLabel("fastSimProducer", "EcalHitsEB", SimHitsEB);
   if (isFastSim) {
     iEvent.getByLabel("fastSimProducer", "EcalHitsEE", SimHitsEE);
     iEvent.getByLabel("fastSimProducer", "EcalHitsES", SimHitsES);
   } else {
     iEvent.getByLabel("g4SimHits", "EcalHitsEB", SimHitsEB);
     iEvent.getByLabel("g4SimHits", "EcalHitsEE", SimHitsEE);
     iEvent.getByLabel("g4SimHits", "EcalHitsES", SimHitsES);
   }
 
   // merge them into one single vector
   auto SimHits = *SimHitsEB;
   SimHits.insert(SimHits.end(), SimHitsEE->begin(), SimHitsEE->end());
   SimHits.insert(SimHits.end(), SimHitsES->begin(), SimHitsES->end());
 
   // As we only had one generated particle (and hopefully no pileup),
   // the total energy is due to the generated particle only
   float energyTotal = 0;
   for (auto const& Hit : SimHits) {
     energyTotal += Hit.energy();
   }
 
   tree_e = genE;
   tree_eta = genEta;
   tree_response = energyTotal / genE;
   tree->Fill();
 }

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

static

Definition at line 99 of file TreeWriterForEcalCorrection.cc.

References edm::ConfigurationDescriptions::add(), and DEFINE_FWK_MODULE.

Referenced by ~TreeWriterForEcalCorrection().

                                                                                              {
   edm::ParameterSetDescription desc;
   descriptions.add("ecalScaleFactorCalculator", desc);
 }