CMS 3D CMS Logo

List of all members | Classes | Public Member Functions | Private Member Functions | Private Attributes
EnergyScaleAnalyzer Class Reference

#include <Validation/EcalClusters/src/EnergyScaleAnalyzer.cc>

Inheritance diagram for EnergyScaleAnalyzer:
edm::EDAnalyzer edm::EDConsumerBase

Classes

struct  tree_structure_
 

Public Member Functions

void analyze (const edm::Event &, const edm::EventSetup &) override
 
void beginJob () override
 
void endJob () override
 
 EnergyScaleAnalyzer (const edm::ParameterSet &)
 
 ~EnergyScaleAnalyzer () override
 
- Public Member Functions inherited from edm::EDAnalyzer
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzer ()
 
SerialTaskQueueglobalLuminosityBlocksQueue ()
 
SerialTaskQueueglobalRunsQueue ()
 
ModuleDescription const & moduleDescription () const
 
std::string workerType () const
 
 ~EDAnalyzer () override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfoconsumesInfo () 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
 
EDConsumerBaseoperator= (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)
 
virtual ~EDConsumerBase () noexcept(false)
 

Private Member Functions

void fillTree (const reco::SuperClusterCollection *scColl, const reco::SuperClusterCollection *corrSCColl, HepMC::GenParticle *mc, tree_structure_ &tree_, float xV, float yV, float zV, int scType)
 

Private Attributes

edm::EDGetTokenT< reco::SuperClusterCollectioncorrectedDynamicHybridSuperClusters_token
 
edm::EDGetTokenT< reco::SuperClusterCollectioncorrectedFixedMatrixSuperClustersWithPreshower_token
 
edm::EDGetTokenT< reco::SuperClusterCollectioncorrectedHybridSuperClusters_token
 
edm::EDGetTokenT< reco::SuperClusterCollectiondynamicHybridSuperClusters_token
 
float energyMax_
 
float etaMax_
 
float etaMaxVtx_
 
float eTMax_
 
float eTMaxVtx_
 
int evtN
 
edm::EDGetTokenT< reco::SuperClusterCollectionfixedMatrixSuperClustersWithPreshower_token
 
edm::EDGetTokenT< edm::HepMCProducthepMCLabel_
 
edm::EDGetTokenT< reco::SuperClusterCollectionhybridSuperClusters_token
 
TTree * mytree_
 
std::string outputFile_
 
float phiMax_
 
float phiMaxVtx_
 
float rClust_vtx_
 
TFile * rootFile_
 
float thetaMax_
 
float thetaMaxVtx_
 
tree_structure_ tree_
 
float xClust_vtx_
 
float xClust_zero_
 
float xVtx_
 
float yClust_vtx_
 
float yClust_zero_
 
float yVtx_
 
float zClust_vtx_
 
float zClust_zero_
 
float zVtx_
 

Additional Inherited Members

- Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- 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 ()
 
- 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

Description: <one line="" class="" summary>="">

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

Definition at line 48 of file EnergyScaleAnalyzer.h.

Constructor & Destructor Documentation

EnergyScaleAnalyzer::EnergyScaleAnalyzer ( const edm::ParameterSet ps)
explicit

Definition at line 73 of file EnergyScaleAnalyzer.cc.

References correctedDynamicHybridSuperClusters_token, correctedFixedMatrixSuperClustersWithPreshower_token, correctedHybridSuperClusters_token, dynamicHybridSuperClusters_token, evtN, fixedMatrixSuperClustersWithPreshower_token, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hepMCLabel_, hybridSuperClusters_token, outputFile_, rootFile_, and AlCaHLTBitMon_QueryRunRegistry::string.

75 {
76 
77  hepMCLabel_ = consumes<edm::HepMCProduct>(ps.getParameter<std::string>("hepMCLabel"));
78  hybridSuperClusters_token = consumes<reco::SuperClusterCollection>(ps.getUntrackedParameter<edm::InputTag>("hybridSuperClusters",edm::InputTag("hybridSuperClusters")));
79  dynamicHybridSuperClusters_token = consumes<reco::SuperClusterCollection>(ps.getUntrackedParameter<edm::InputTag>("dynamicHybridSuperClusters",edm::InputTag("dynamicHybridSuperClusters")));
80  correctedHybridSuperClusters_token = consumes<reco::SuperClusterCollection>(ps.getUntrackedParameter<edm::InputTag>("correctedHybridSuperClusters",edm::InputTag("correctedHybridSuperClusters")));
81  correctedDynamicHybridSuperClusters_token = consumes<reco::SuperClusterCollection>(ps.getUntrackedParameter<edm::InputTag>("correctedDynamicHybridSuperClusters",edm::InputTag("correctedDynamicHybridSuperClusters")));
82  correctedFixedMatrixSuperClustersWithPreshower_token = consumes<reco::SuperClusterCollection>(ps.getUntrackedParameter<edm::InputTag>("correctedFixedMatrixSuperClustersWithPreshower",edm::InputTag("correctedFixedMatrixSuperClustersWithPreshower")));
83  fixedMatrixSuperClustersWithPreshower_token = consumes<reco::SuperClusterCollection>(ps.getUntrackedParameter<edm::InputTag>("fixedMatrixSuperClustersWithPreshower",edm::InputTag("fixedMatrixSuperClustersWithPreshower")));
84 
85  outputFile_ = ps.getParameter<std::string>("outputFile");
86  rootFile_ = TFile::Open(outputFile_.c_str(),"RECREATE"); // open output file to store histograms
87 
88  evtN = 0;
89 }
edm::EDGetTokenT< reco::SuperClusterCollection > dynamicHybridSuperClusters_token
T getParameter(std::string const &) const
T getUntrackedParameter(std::string const &, T const &) const
edm::EDGetTokenT< reco::SuperClusterCollection > fixedMatrixSuperClustersWithPreshower_token
edm::EDGetTokenT< edm::HepMCProduct > hepMCLabel_
edm::EDGetTokenT< reco::SuperClusterCollection > correctedFixedMatrixSuperClustersWithPreshower_token
edm::EDGetTokenT< reco::SuperClusterCollection > hybridSuperClusters_token
edm::EDGetTokenT< reco::SuperClusterCollection > correctedDynamicHybridSuperClusters_token
edm::EDGetTokenT< reco::SuperClusterCollection > correctedHybridSuperClusters_token
EnergyScaleAnalyzer::~EnergyScaleAnalyzer ( )
override

Definition at line 93 of file EnergyScaleAnalyzer.cc.

References rootFile_.

95 {
96  delete rootFile_;
97 }

Member Function Documentation

void EnergyScaleAnalyzer::analyze ( const edm::Event evt,
const edm::EventSetup es 
)
override

Definition at line 121 of file EnergyScaleAnalyzer.cc.

References funct::abs(), correctedDynamicHybridSuperClusters_token, correctedFixedMatrixSuperClustersWithPreshower_token, correctedHybridSuperClusters_token, kinem::delta_R(), dynamicHybridSuperClusters_cfi::dynamicHybridSuperClusters, dynamicHybridSuperClusters_token, fillTree(), fixedMatrixSuperClustersWithPreshower_token, nano_cff::genEvent, GenParticle::GenParticle, edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), edm::Event::getManyByType(), MuonClassifier_cff::hepMC, hepMCLabel_, unifiedSCCollection_cfi::hybridSuperClusters, hybridSuperClusters_token, mps_fire::i, edm::HandleBase::isValid(), checklumidiff::l, edm::EDConsumerBase::labelsForToken(), hiMultiTrackSelector_cfi::max_eta, CaloTowersParam_cfi::mc, EnergyScaleAnalyzer::tree_structure_::mc_e, EnergyScaleAnalyzer::tree_structure_::mc_et, EnergyScaleAnalyzer::tree_structure_::mc_eta, EnergyScaleAnalyzer::tree_structure_::mc_npar, EnergyScaleAnalyzer::tree_structure_::mc_phi, EnergyScaleAnalyzer::tree_structure_::mc_sep, EnergyScaleAnalyzer::tree_structure_::mc_theta, edm::ProductLabels::module, AlCaHLTBitMon_ParallelJobs::p, EnergyScaleAnalyzer::tree_structure_::parID, edm::Handle< T >::product(), funct::sin(), AlCaHLTBitMon_QueryRunRegistry::string, tree_, badGlobalMuonTaggersAOD_cff::vtx, xVtx_, yVtx_, and zVtx_.

121  {
122  using namespace edm; // needed for all fwk related classes
123  using namespace std;
124 
125  // std::cout << "Proceccing event # " << ++evtN << std::endl;
126 
127  //Get containers for MC truth, SC etc. ===================================================
128  // =======================================================================================
129  // =======================================================================================
131  evt.getByToken( hepMCLabel_, hepMC ) ;
132 
133  Labels l;
135 
136  const HepMC::GenEvent* genEvent = hepMC->GetEvent();
137  if ( !(hepMC.isValid())) {
138  LogInfo("EnergyScaleAnalyzer") << "Could not get MC Product!";
139  return;
140  }
141 
142 
143  //=======================For Vertex correction
144  std::vector< Handle< HepMCProduct > > evtHandles ;
145  evt.getManyByType( evtHandles ) ;
146 
147  for ( unsigned int i=0; i<evtHandles.size(); ++i) {
148  if ( evtHandles[i].isValid() ) {
149  const HepMC::GenEvent* evt = evtHandles[i]->GetEvent() ;
150 
151  // take only 1st vertex for now - it's been tested only of PGuns...
152  //
153  HepMC::GenEvent::vertex_const_iterator vtx = evt->vertices_begin() ;
154  if ( evtHandles[i].provenance()->moduleLabel() == std::string(l.module) ) {
155  //Corrdinates of Vertex w.r.o. the point (0,0,0)
156  xVtx_ = 0.1*(*vtx)->position().x();
157  yVtx_ = 0.1*(*vtx)->position().y();
158  zVtx_ = 0.1*(*vtx)->position().z();
159  }
160  }
161  }
162  //==============================================================================
163  //Get handle to SC collections
164 
166  try {
167  evt.getByToken(hybridSuperClusters_token,hybridSuperClusters);
168  }catch (cms::Exception& ex) {
169  edm::LogError("EnergyScaleAnalyzer") << "Can't get collection with producer hybridSuperClusters.";
170  }
171 
173  try {
174  evt.getByToken(dynamicHybridSuperClusters_token,dynamicHybridSuperClusters);
175  }catch (cms::Exception& ex) {
176  edm::LogError("EnergyScaleAnalyzer") << "Can't get collection with producer dynamicHybridSuperClusters.";
177  }
178 
179  Handle<reco::SuperClusterCollection> fixedMatrixSuperClustersWithPS;
180  try {
181  evt.getByToken(fixedMatrixSuperClustersWithPreshower_token,fixedMatrixSuperClustersWithPS);
182  }catch (cms::Exception& ex) {
183  edm::LogError("EnergyScaleAnalyzer") << "Can't get collection with producer fixedMatrixSuperClustersWithPreshower.";
184  }
185 
186  //Corrected collections
187  Handle<reco::SuperClusterCollection> correctedHybridSC;
188  try {
189  evt.getByToken(correctedHybridSuperClusters_token,correctedHybridSC);
190  }catch (cms::Exception& ex) {
191  edm::LogError("EnergyScaleAnalyzer") << "Can't get collection with producer correctedHybridSuperClusters.";
192  }
193 
194  Handle<reco::SuperClusterCollection> correctedDynamicHybridSC;
195  try{
196  evt.getByToken(correctedDynamicHybridSuperClusters_token,correctedDynamicHybridSC);
197  }catch (cms::Exception& ex) {
198  edm::LogError("EnergyScaleAnalyzer") << "Can't get collection with producer correctedDynamicHybridSuperClusters.";
199  }
200 
201  Handle<reco::SuperClusterCollection> correctedFixedMatrixSCWithPS;
202  try {
203  evt.getByToken(correctedFixedMatrixSuperClustersWithPreshower_token,correctedFixedMatrixSCWithPS);
204  }catch (cms::Exception& ex ) {
205  edm::LogError("EnergyScaleAnalyzer") << "Can't get collection with producer correctedFixedMatrixSuperClustersWithPreshower.";
206  }
207 
208  const reco::SuperClusterCollection* dSC = dynamicHybridSuperClusters.product();
209  const reco::SuperClusterCollection* hSC = hybridSuperClusters.product();
210  const reco::SuperClusterCollection* fmSC = fixedMatrixSuperClustersWithPS.product();
211  const reco::SuperClusterCollection* chSC = correctedHybridSC.product();
212  const reco::SuperClusterCollection* cdSC = correctedDynamicHybridSC.product();
213  const reco::SuperClusterCollection* cfmSC = correctedFixedMatrixSCWithPS.product();
214 
215  // ----------------------- Print outs for debugging
216  /*
217  std::cout << "MC truth" << std::endl;
218  int counterI = 0;
219  for( HepMC::GenEvent::particle_const_iterator p = genEvent->particles_begin();
220  p != genEvent->particles_end(); ++p ) {
221  if ( fabs((*p)->momentum().eta()) < 1.5 ) {
222  std::cout << ++counterI << " " << (*p)->momentum().e() << " "
223  << (*p)->momentum().phi() << " " << (*p)->momentum().eta() << std::endl;
224  }
225  }
226 
227  std::cout << "Standard clusters" << std::endl;
228  counterI = 0;
229  for(reco::SuperClusterCollection::const_iterator em = hSC->begin();
230  em != hSC->end(); ++em )
231  std::cout << ++counterI << " " << em->energy() << " " << em->position().phi() << " " << em->position().eta() << std::endl;
232 
233  std::cout << "Dynamic clusters" << std::endl;
234  counterI = 0;
235  for(reco::SuperClusterCollection::const_iterator em = dSC->begin();
236  em != dSC->end(); ++em )
237  std::cout << ++counterI << " " << em->energy() << " " << em->position().phi() << " " << em->position().eta() << std::endl;
238 
239  std::cout << "FixedMatrix clusters with PS" << std::endl;
240  counterI = 0;
241  for(reco::SuperClusterCollection::const_iterator em = fmSC->begin();
242  em != fmSC->end(); ++em )
243  std::cout << ++counterI << " " << em->energy() << " " << em->position().phi() << " " << em->position().eta() << std::endl;
244  */
245  // -----------------------------
246  //=====================================================================
247  // All containers are loaded, perform the analysis
248  //====================================================================
249 
250  // --------------------------- Store MC particles
251  HepMC::GenEvent::particle_const_iterator p = genEvent->particles_begin();
252 
253  // Search for MC electrons or photons that satisfy the criteria
254  float min_eT = 5;
255  float max_eta = 2.5;
256 
257  std::vector<HepMC::GenParticle* > mcParticles;
258  //int counter = 0;
259  for ( HepMC::GenEvent::particle_const_iterator p = genEvent->particles_begin();
260  p != genEvent->particles_end();
261  ++p ) {
262  //LogInfo("EnergyScaleAnalyzer") << "Particle " << ++counter
263  //<< " PDG ID = " << (*p)->pdg_id() << " pT = " << (*p)->momentum().perp();
264  // require photon or electron
265  if ( (*p)->pdg_id() != 22 && abs((*p)->pdg_id()) != 11 ) continue;
266 
267  // require selection criteria
268  bool satisfySelectionCriteria =
269  (*p)->momentum().perp() > min_eT &&
270  fabs((*p)->momentum().eta()) < max_eta;
271 
272  if ( ! satisfySelectionCriteria ) continue;
273 
274  // EM MC particle is found, save it in the vector
275  mcParticles.push_back(*p);
276  }
277  // separation in dR between 2 first MC particles
278  // should not be used for MC samples with > 2 em objects generated!
279  if ( mcParticles.size() == 2 ) {
280  HepMC::GenParticle* mc1 = mcParticles[0];
281  HepMC::GenParticle* mc2 = mcParticles[1];
282  tree_.mc_sep = kinem::delta_R(mc1->momentum().eta(), mc1->momentum().phi(),
283  mc2->momentum().eta(), mc2->momentum().phi());
284  } else
285  tree_.mc_sep = -100;
286 
287 
288  // now loop over MC particles, find the match with SC and do everything we need
289  // then save info in the tree for every MC particle
290  for(std::vector<HepMC::GenParticle* >::const_iterator p = mcParticles.begin();
291  p != mcParticles.end(); ++p) {
292  HepMC::GenParticle* mc = *p;
293 
294  // Fill MC information
295  tree_.mc_npar = mcParticles.size();
296  tree_.parID = mc->pdg_id();
297  tree_.mc_e = mc->momentum().e();
298  tree_.mc_et = mc->momentum().e()*sin(mc->momentum().theta());
299  tree_.mc_phi = mc->momentum().phi();
300  tree_.mc_eta = mc->momentum().eta();
301  tree_.mc_theta = mc->momentum().theta();
302 
303 
304  //Call function to fill tree
305  // scType coprreponds:
306  // HybridSuperCluster -- 1
307  // DynamicHybridSuperCluster -- 2
308  // FixedMatrixSuperClustersWithPreshower -- 3
309 
310  fillTree( hSC, chSC, mc, tree_, xVtx_, yVtx_, zVtx_, 1);
311  // std::cout << " TYPE " << 1 << " : " << tree_.em_e << " : " << tree_.em_phi << " : " << tree_.em_eta << std::endl;
312 
313  fillTree( dSC, cdSC, mc, tree_, xVtx_, yVtx_, zVtx_, 2);
314  // std::cout << " TYPE " << 2 << " : " << tree_.em_e << " : " << tree_.em_phi << " : " << tree_.em_eta << std::endl;
315 
316  fillTree( fmSC, cfmSC, mc, tree_, xVtx_, yVtx_, zVtx_, 3);
317  // std::cout << " TYPE " << 3 << " : " << tree_.em_e << " : " << tree_.em_phi << " : " << tree_.em_eta << std::endl;
318 
319  // mytree_->Fill();
320  } // loop over particles
321 }
void getManyByType(std::vector< Handle< PROD > > &results) const
Definition: Event.h:509
edm::EDGetTokenT< reco::SuperClusterCollection > dynamicHybridSuperClusters_token
edm::EDGetTokenT< reco::SuperClusterCollection > fixedMatrixSuperClustersWithPreshower_token
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:519
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
edm::EDGetTokenT< edm::HepMCProduct > hepMCLabel_
edm::EDGetTokenT< reco::SuperClusterCollection > correctedFixedMatrixSuperClustersWithPreshower_token
edm::EDGetTokenT< reco::SuperClusterCollection > hybridSuperClusters_token
void fillTree(const reco::SuperClusterCollection *scColl, const reco::SuperClusterCollection *corrSCColl, HepMC::GenParticle *mc, tree_structure_ &tree_, float xV, float yV, float zV, int scType)
std::vector< SuperCluster > SuperClusterCollection
collection of SuperCluser objectr
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
double delta_R(double eta1, double phi1, double eta2, double phi2)
Definition: AnglesUtil.h:116
bool isValid() const
Definition: HandleBase.h:74
edm::EDGetTokenT< reco::SuperClusterCollection > correctedDynamicHybridSuperClusters_token
edm::EDGetTokenT< reco::SuperClusterCollection > correctedHybridSuperClusters_token
ProductLabels Labels
const HepMC::GenEvent * GetEvent() const
Definition: HepMCProduct.h:38
T const * product() const
Definition: Handle.h:81
void labelsForToken(EDGetToken iToken, Labels &oLabels) const
HLT enums.
void EnergyScaleAnalyzer::beginJob ( void  )
overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 101 of file EnergyScaleAnalyzer.cc.

References EnergyScaleAnalyzer::tree_structure_::mc_npar, mytree_, and tree_.

101  {
102 //========================================================================
103 
104  mytree_ = new TTree("energyScale","");
105  TString treeVariables = "mc_npar/I:parID:mc_sep/F:mc_e:mc_et:mc_phi:mc_eta:mc_theta:"; // MC information
106  treeVariables += "em_dR/F:"; // MC <-> EM matching information
107  treeVariables += "em_isInCrack/I:em_scType:em_e/F:em_et:em_phi:em_eta:em_theta:em_nCell/I:em_nBC:"; // EM SC info
108  treeVariables += "em_pet/F:em_pe:em_peta:em_ptheta:" ; // EM SC physics (eta corrected information)
109 
110  treeVariables += "emCorr_e/F:emCorr_et:emCorr_eta:emCorr_phi:emCorr_theta:";// CMSSW standard corrections
111  treeVariables += "emCorr_pet/F:emCorr_peta:emCorr_ptheta:" ;// CMSSW standard physics
112 
113  treeVariables += "em_pw/F:em_ew:em_br" ; // EM widths pw -- phiWidth, ew -- etaWidth, ratios of pw/ew
114 
115  mytree_->Branch("energyScale",&(tree_.mc_npar),treeVariables);
116 
117 }
void EnergyScaleAnalyzer::endJob ( void  )
overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 450 of file EnergyScaleAnalyzer.cc.

References DEFINE_FWK_MODULE, and rootFile_.

450  {
451  //========================================================================
452  //Fill ROOT tree
453  rootFile_->Write();
454 }
void EnergyScaleAnalyzer::fillTree ( const reco::SuperClusterCollection scColl,
const reco::SuperClusterCollection corrSCColl,
HepMC::GenParticle *  mc,
tree_structure_ tree_,
float  xV,
float  yV,
float  zV,
int  scType 
)
private

fill tree with kinematic variables of corrected Super Cluster

Definition at line 324 of file EnergyScaleAnalyzer.cc.

References kinem::delta_R(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, EnergyScaleAnalyzer::tree_structure_::em_br, EnergyScaleAnalyzer::tree_structure_::em_dR, EnergyScaleAnalyzer::tree_structure_::em_e, EnergyScaleAnalyzer::tree_structure_::em_et, EnergyScaleAnalyzer::tree_structure_::em_eta, EnergyScaleAnalyzer::tree_structure_::em_ew, EnergyScaleAnalyzer::tree_structure_::em_isInCrack, EnergyScaleAnalyzer::tree_structure_::em_nBC, EnergyScaleAnalyzer::tree_structure_::em_nCell, EnergyScaleAnalyzer::tree_structure_::em_pe, EnergyScaleAnalyzer::tree_structure_::em_pet, EnergyScaleAnalyzer::tree_structure_::em_peta, EnergyScaleAnalyzer::tree_structure_::em_phi, EnergyScaleAnalyzer::tree_structure_::em_ptheta, EnergyScaleAnalyzer::tree_structure_::em_pw, EnergyScaleAnalyzer::tree_structure_::em_scType, EnergyScaleAnalyzer::tree_structure_::em_theta, EnergyScaleAnalyzer::tree_structure_::emCorr_e, EnergyScaleAnalyzer::tree_structure_::emCorr_et, EnergyScaleAnalyzer::tree_structure_::emCorr_eta, EnergyScaleAnalyzer::tree_structure_::emCorr_pet, EnergyScaleAnalyzer::tree_structure_::emCorr_peta, EnergyScaleAnalyzer::tree_structure_::emCorr_phi, EnergyScaleAnalyzer::tree_structure_::emCorr_ptheta, EnergyScaleAnalyzer::tree_structure_::emCorr_theta, energyMax_, PVValHelper::eta, etaMax_, etaMaxVtx_, eTMax_, eTMaxVtx_, cmsBatch::log, mytree_, phi, phiMax_, phiMaxVtx_, rClust_vtx_, funct::sin(), mathSSE::sqrt(), funct::tan(), thetaMax_, thetaMaxVtx_, xClust_vtx_, xClust_zero_, yClust_vtx_, yClust_zero_, zClust_vtx_, and zClust_zero_.

Referenced by analyze().

325  {
326 
327  // ----------------------------- SuperClusters before energy correction
328  reco::SuperClusterCollection::const_iterator em = scColl->end();
329  float energyMax = -100.0; // dummy energy of the matched SC
330  for(reco::SuperClusterCollection::const_iterator aClus = scColl->begin();
331  aClus != scColl->end(); ++aClus) {
332  // check the matching
333  float dR = kinem::delta_R(mc ->momentum().eta(), mc ->momentum().phi(),
334  aClus->position().eta(), aClus->position().phi());
335  if (dR < 0.4) { // a rather loose matching cut
336  float energy = aClus->energy();
337  if ( energy < energyMax ) continue;
338  energyMax = energy;
339  em = aClus;
340  }
341  }
342 
343  if (em == scColl->end() ) {
344  // std::cout << "No matching SC with type " << scType << " was found for MC particle! " << std::endl;
345  // std::cout << "Going to next type of SC. " << std::endl;
346  return;
347  }
348  // ------------
349 
350  tree_.em_scType = scType;
351 
352  tree_.em_isInCrack = 0;
353  double emAbsEta = fabs(em->position().eta());
354  // copied from RecoEgama/EgammaElectronAlgos/src/EgammaElectronClassification.cc
355  if ( emAbsEta < 0.018 ||
356  (emAbsEta > 0.423 && emAbsEta < 0.461) ||
357  (emAbsEta > 0.770 && emAbsEta < 0.806) ||
358  (emAbsEta > 1.127 && emAbsEta < 1.163) ||
359  (emAbsEta > 1.460 && emAbsEta < 1.558) )
360  tree_.em_isInCrack = 1;
361 
362  tree_.em_dR = kinem::delta_R(mc->momentum().eta(), mc->momentum().phi(),
363  em->position().eta(), em->position().phi());
364  tree_.em_e = em->energy();
365  tree_.em_et = em->energy() * sin(em->position().theta());
366  tree_.em_phi = em->position().phi();
367  tree_.em_eta = em->position().eta();
368  tree_.em_theta = em->position().theta();
369  tree_.em_nCell = em->size();
370  tree_.em_nBC = em->clustersSize();
371 
372  //Get physics e, et etc:
373  //Coordinates of EM object with respect of the point (0,0,0)
374  xClust_zero_ = em->position().x();
375  yClust_zero_ = em->position().y();
376  zClust_zero_ = em->position().z();
377  //Coordinates of EM object w.r.o. the Vertex position
378  xClust_vtx_ = xClust_zero_ - xV;
379  yClust_vtx_ = yClust_zero_ - yV;
380  zClust_vtx_ = zClust_zero_ - zV;
381 
382  energyMax_ = em->energy();
383  thetaMax_ = em->position().theta();
384  etaMax_ = em->position().eta();
385  phiMax_ = em->position().phi();
387  if ( phiMax_ < 0) phiMax_ += 2*3.14159;
388 
391  etaMaxVtx_ = - log(tan(thetaMaxVtx_/2));
394  if ( phiMaxVtx_ < 0) phiMaxVtx_ += 2*3.14159;
395  //=============================
396  //parametres of EM object after vertex correction
401 
402 
403  //------------------------------- Get SC after energy correction
404  em = corrSCColl->end();
405  energyMax = -100.0; // dummy energy of the matched SC
406  for(reco::SuperClusterCollection::const_iterator aClus = corrSCColl->begin();
407  aClus != corrSCColl->end(); ++aClus) {
408  // check the matching
409  float dR = kinem::delta_R(mc ->momentum().eta(), mc ->momentum().phi(),
410  aClus->position().eta(), aClus->position().phi());
411  if (dR < 0.4) {
412  float energy = aClus->energy();
413  if ( energy < energyMax ) continue;
414  energyMax = energy;
415  em = aClus;
416  }
417  }
418 
419  if (em == corrSCColl->end() ) {
420  // std::cout << "No matching corrected SC with type " << scType << " was found for MC particle! " << std::endl;
421  // std::cout << "Going to next type of SC. " << std::endl;
422  return;
423  }
424  // ------------
425 
427  tree_.emCorr_e = em->energy();
428  tree_.emCorr_et = em->energy() * sin(em->position().theta());
429  tree_.emCorr_phi = em->position().phi();
430  tree_.emCorr_eta = em->position().eta();
431  tree_.emCorr_theta = em->position().theta();
432 
433  // =========== Eta and Theta wrt Vertex does not change after energy corrections are applied
434  // =========== So, no need to calculate them again
435 
439 
440  tree_.em_pw = em->phiWidth();
441  tree_.em_ew = em->etaWidth();
443 
444  mytree_->Fill();
445 
446 }
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
T sqrt(T t)
Definition: SSEVec.h:18
Tan< T >::type tan(const T &t)
Definition: Tan.h:22
double delta_R(double eta1, double phi1, double eta2, double phi2)
Definition: AnglesUtil.h:116

Member Data Documentation

edm::EDGetTokenT<reco::SuperClusterCollection> EnergyScaleAnalyzer::correctedDynamicHybridSuperClusters_token
private

Definition at line 67 of file EnergyScaleAnalyzer.h.

Referenced by analyze(), and EnergyScaleAnalyzer().

edm::EDGetTokenT<reco::SuperClusterCollection> EnergyScaleAnalyzer::correctedFixedMatrixSuperClustersWithPreshower_token
private

Definition at line 68 of file EnergyScaleAnalyzer.h.

Referenced by analyze(), and EnergyScaleAnalyzer().

edm::EDGetTokenT<reco::SuperClusterCollection> EnergyScaleAnalyzer::correctedHybridSuperClusters_token
private

Definition at line 66 of file EnergyScaleAnalyzer.h.

Referenced by analyze(), and EnergyScaleAnalyzer().

edm::EDGetTokenT<reco::SuperClusterCollection> EnergyScaleAnalyzer::dynamicHybridSuperClusters_token
private

Definition at line 65 of file EnergyScaleAnalyzer.h.

Referenced by analyze(), and EnergyScaleAnalyzer().

float EnergyScaleAnalyzer::energyMax_
private

Definition at line 140 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

float EnergyScaleAnalyzer::etaMax_
private

Definition at line 143 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

float EnergyScaleAnalyzer::etaMaxVtx_
private

Definition at line 144 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

float EnergyScaleAnalyzer::eTMax_
private

Definition at line 141 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

float EnergyScaleAnalyzer::eTMaxVtx_
private

Definition at line 142 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

int EnergyScaleAnalyzer::evtN
private

Definition at line 151 of file EnergyScaleAnalyzer.h.

Referenced by EnergyScaleAnalyzer().

edm::EDGetTokenT<reco::SuperClusterCollection> EnergyScaleAnalyzer::fixedMatrixSuperClustersWithPreshower_token
private

Definition at line 69 of file EnergyScaleAnalyzer.h.

Referenced by analyze(), and EnergyScaleAnalyzer().

edm::EDGetTokenT<edm::HepMCProduct> EnergyScaleAnalyzer::hepMCLabel_
private

Definition at line 63 of file EnergyScaleAnalyzer.h.

Referenced by analyze(), and EnergyScaleAnalyzer().

edm::EDGetTokenT<reco::SuperClusterCollection> EnergyScaleAnalyzer::hybridSuperClusters_token
private

Definition at line 64 of file EnergyScaleAnalyzer.h.

Referenced by analyze(), and EnergyScaleAnalyzer().

TTree* EnergyScaleAnalyzer::mytree_
private

Definition at line 75 of file EnergyScaleAnalyzer.h.

Referenced by beginJob(), and fillTree().

std::string EnergyScaleAnalyzer::outputFile_
private

Definition at line 61 of file EnergyScaleAnalyzer.h.

Referenced by EnergyScaleAnalyzer().

float EnergyScaleAnalyzer::phiMax_
private

Definition at line 145 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

float EnergyScaleAnalyzer::phiMaxVtx_
private

Definition at line 146 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

float EnergyScaleAnalyzer::rClust_vtx_
private

Definition at line 138 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

TFile* EnergyScaleAnalyzer::rootFile_
private

Definition at line 72 of file EnergyScaleAnalyzer.h.

Referenced by endJob(), EnergyScaleAnalyzer(), and ~EnergyScaleAnalyzer().

float EnergyScaleAnalyzer::thetaMax_
private

Definition at line 147 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

float EnergyScaleAnalyzer::thetaMaxVtx_
private

Definition at line 148 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

tree_structure_ EnergyScaleAnalyzer::tree_
private

Definition at line 123 of file EnergyScaleAnalyzer.h.

Referenced by analyze(), and beginJob().

float EnergyScaleAnalyzer::xClust_vtx_
private

Definition at line 134 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

float EnergyScaleAnalyzer::xClust_zero_
private

Definition at line 130 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

float EnergyScaleAnalyzer::xVtx_
private

Definition at line 126 of file EnergyScaleAnalyzer.h.

Referenced by analyze().

float EnergyScaleAnalyzer::yClust_vtx_
private

Definition at line 135 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

float EnergyScaleAnalyzer::yClust_zero_
private

Definition at line 131 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

float EnergyScaleAnalyzer::yVtx_
private

Definition at line 127 of file EnergyScaleAnalyzer.h.

Referenced by analyze().

float EnergyScaleAnalyzer::zClust_vtx_
private

Definition at line 136 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

float EnergyScaleAnalyzer::zClust_zero_
private

Definition at line 132 of file EnergyScaleAnalyzer.h.

Referenced by fillTree().

float EnergyScaleAnalyzer::zVtx_
private

Definition at line 128 of file EnergyScaleAnalyzer.h.

Referenced by analyze().