Inheritance diagram for TreeWriterForEcalCorrection:

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

	~TreeWriterForEcalCorrection ()

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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 &)

Private Member Functions
virtual 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 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 24 of file TreeWriterForEcalCorrection.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 40 of file TreeWriterForEcalCorrection.cc.

References file, TFileService::make(), tree, 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 ( )

inline

Definition at line 27 of file TreeWriterForEcalCorrection.cc.

27 {};

Member Function Documentation

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

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 49 of file TreeWriterForEcalCorrection.cc.

References reco::LeafCandidate::energy(), reco::LeafCandidate::eta(), Exception, relval_steps::gen(), edm::Event::getByLabel(), tree, tree_e, tree_eta, and tree_response.

 {
   // 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( "famosSimHits", "EcalHitsEB", SimHitsEB );
   if( isFastSim ) {
     iEvent.getByLabel( "famosSimHits", "EcalHitsEE", SimHitsEE );
     iEvent.getByLabel( "famosSimHits", "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 107 of file TreeWriterForEcalCorrection.cc.

References edm::ConfigurationDescriptions::add().

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