CMS 3D CMS Logo
List of all members | Public Member Functions | Private Attributes
EcalTrigPrimAnalyzer Class Reference

#include <EcalTrigPrimAnalyzer.h>

Inheritance diagram for EcalTrigPrimAnalyzer:
edm::one::EDAnalyzer<> edm::one::EDAnalyzerBase edm::EDConsumerBase

Public Member Functions

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

Private Attributes

TH1I * ecal_et_ [2]
 
TH1I * ecal_fgvb_ [2]
 
std::vector< std::string > ecal_parts_
 
TH1I * ecal_tt_ [2]
 
float eRec_
 
int fg_
 
TH1I * histBar
 
TH1I * histEndc
 
TFile * histfile_
 
TH1F * hTPoverRechit_
 
TH2F * hTPvsRechit_
 
int ieta_
 
int iphi_
 
edm::InputTag label_
 
bool recHits_
 
edm::InputTag rechits_labelEB_
 
edm::InputTag rechits_labelEE_
 
int tpgADC_
 
float tpgGeV_
 
TTree * tree_
 
int ttf_
 

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

Description: test of the output of EcalTrigPrimProducer

Description: rereads the result of the EcalTrigPrimProducer

Definition at line 37 of file EcalTrigPrimAnalyzer.h.

Constructor & Destructor Documentation

EcalTrigPrimAnalyzer::EcalTrigPrimAnalyzer ( const edm::ParameterSet &  iConfig)
explicit

Definition at line 50 of file EcalTrigPrimAnalyzer.cc.

References edm::ParameterSet::getParameter(), mps_fire::i, and fftjetcommon_cfi::title.

52 {
53  ecal_parts_.push_back("Barrel");
54  ecal_parts_.push_back("Endcap");
55 
56  histfile_ = new TFile("histos.root", "RECREATE");
57  tree_ = new TTree("TPGtree", "TPGtree");
58  tree_->Branch("iphi", &iphi_, "iphi/I");
59  tree_->Branch("ieta", &ieta_, "ieta/I");
60  tree_->Branch("eRec", &eRec_, "eRec/F");
61  tree_->Branch("tpgADC", &tpgADC_, "tpgADC/I");
62  tree_->Branch("tpgGeV", &tpgGeV_, "tpgGeV/F");
63  tree_->Branch("ttf", &ttf_, "ttf/I");
64  tree_->Branch("fg", &fg_, "fg/I");
65  for (unsigned int i = 0; i < 2; ++i) {
66  ecal_et_[i] = new TH1I(ecal_parts_[i].c_str(), "Et", 255, 0, 255);
67  char title[30];
68  sprintf(title, "%s_ttf", ecal_parts_[i].c_str());
69  ecal_tt_[i] = new TH1I(title, "TTF", 10, 0, 10);
70  sprintf(title, "%s_fgvb", ecal_parts_[i].c_str());
71  ecal_fgvb_[i] = new TH1I(title, "FGVB", 10, 0, 10);
72  }
73 
74  recHits_ = iConfig.getParameter<bool>("AnalyzeRecHits");
75  label_ = iConfig.getParameter<edm::InputTag>("inputTP");
76  if (recHits_) {
77  hTPvsRechit_ =
78  new TH2F("TP_vs_RecHit", "TP vs rechit", 256, -1, 255, 255, 0, 255);
79  hTPoverRechit_ = new TH1F("TP_over_RecHit", "TP over rechit", 500, 0, 4);
80  rechits_labelEB_ = iConfig.getParameter<edm::InputTag>("inputRecHitsEB");
81  rechits_labelEE_ = iConfig.getParameter<edm::InputTag>("inputRecHitsEE");
82  }
83 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
mps_fire.i
i
Definition: mps_fire.py:338
EcalTrigPrimAnalyzer::ecal_parts_
std::vector< std::string > ecal_parts_
Definition: EcalTrigPrimAnalyzer.h:47
EcalTrigPrimAnalyzer::rechits_labelEE_
edm::InputTag rechits_labelEE_
Definition: EcalTrigPrimAnalyzer.h:63
EcalTrigPrimAnalyzer::rechits_labelEB_
edm::InputTag rechits_labelEB_
Definition: EcalTrigPrimAnalyzer.h:62
fftjetcommon_cfi.title
title
Definition: fftjetcommon_cfi.py:32
EcalTrigPrimAnalyzer::~EcalTrigPrimAnalyzer ( )
override

Definition at line 85 of file EcalTrigPrimAnalyzer.cc.

85  {
86 
87  // do anything here that needs to be done at desctruction time
88  // (e.g. close files, deallocate resources etc.)
89 
90  histfile_->Write();
91  histfile_->Close();
92 }

Member Function Documentation

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

Definition at line 99 of file EcalTrigPrimAnalyzer.cc.

References EcalTriggerPrimitiveDigi::compressedEt(), KineDebug3::count(), edmIntegrityCheck::d, EcalTriggerPrimitiveDigi::fineGrain(), edm::EventSetup::get(), edm::Event::getByLabel(), CaloSubdetectorGeometry::getGeometry(), EcalTPGScale::getTPGInGeV(), mps_fire::i, EcalTriggerPrimitiveDigi::id(), EcalTrigTowerDetId::ieta(), EcalTrigTowerDetId::ietaAbs(), EcalTrigTowerDetId::iphi(), edm::Handle< T >::product(), edm::ESHandle< T >::product(), EcalTPGScale::setEventSetup(), funct::sin(), edm::SortedCollection< T, SORT >::size(), EcalTrigTowerDetId::subDet(), theta(), EBDetId::tower(), and EcalTriggerPrimitiveDigi::ttFlag().

100  {
101  using namespace edm;
102  using namespace std;
103 
104  // Get input
105  edm::Handle<EcalTrigPrimDigiCollection> tp;
106  iEvent.getByLabel(label_, tp);
107  for (unsigned int i = 0; i < tp.product()->size(); i++) {
108  EcalTriggerPrimitiveDigi d = (*(tp.product()))[i];
109  int subdet = d.id().subDet() - 1;
110  if (subdet == 0) {
111  ecal_et_[subdet]->Fill(d.compressedEt());
112  } else {
113  if (d.id().ietaAbs() == 27 || d.id().ietaAbs() == 28) {
114  if (i % 2)
115  ecal_et_[subdet]->Fill(d.compressedEt() * 2.);
116  } else
117  ecal_et_[subdet]->Fill(d.compressedEt());
118  }
119  ecal_tt_[subdet]->Fill(d.ttFlag());
120  ecal_fgvb_[subdet]->Fill(d.fineGrain());
121  }
122  if (!recHits_)
123  return;
124 
125  // comparison with RecHits
126  edm::Handle<EcalRecHitCollection> rechit_EB_col;
127  iEvent.getByLabel(rechits_labelEB_, rechit_EB_col);
128 
129  edm::Handle<EcalRecHitCollection> rechit_EE_col;
130  iEvent.getByLabel(rechits_labelEE_, rechit_EE_col);
131 
132  edm::ESHandle<CaloGeometry> theGeometry;
133  edm::ESHandle<CaloSubdetectorGeometry> theBarrelGeometry_handle;
134  edm::ESHandle<CaloSubdetectorGeometry> theEndcapGeometry_handle;
135  iSetup.get<CaloGeometryRecord>().get(theGeometry);
136  iSetup.get<EcalEndcapGeometryRecord>().get("EcalEndcap",
137  theEndcapGeometry_handle);
138  iSetup.get<EcalBarrelGeometryRecord>().get("EcalBarrel",
139  theBarrelGeometry_handle);
140 
141  const CaloSubdetectorGeometry *theEndcapGeometry =
142  theEndcapGeometry_handle.product();
143  const CaloSubdetectorGeometry *theBarrelGeometry =
144  theBarrelGeometry_handle.product();
145  edm::ESHandle<EcalTrigTowerConstituentsMap> eTTmap_;
146  iSetup.get<IdealGeometryRecord>().get(eTTmap_);
147 
148  map<EcalTrigTowerDetId, float> mapTow_Et;
149 
150  for (unsigned int i = 0; i < rechit_EB_col.product()->size(); i++) {
151  const EBDetId &myid1 = (*rechit_EB_col.product())[i].id();
152  EcalTrigTowerDetId towid1 = myid1.tower();
153  float theta = theBarrelGeometry->getGeometry(myid1)->getPosition().theta();
154  float Etsum = ((*rechit_EB_col.product())[i].energy()) * sin(theta);
155  bool test_alreadyin = false;
156  map<EcalTrigTowerDetId, float>::iterator ittest = mapTow_Et.find(towid1);
157  if (ittest != mapTow_Et.end())
158  test_alreadyin = true;
159  if (test_alreadyin)
160  continue;
161  unsigned int j = i + 1;
162  bool loopend = false;
163  unsigned int count = 0;
164  while (j < rechit_EB_col.product()->size() && !loopend) {
165  count++;
166  const EBDetId &myid2 = (*rechit_EB_col.product())[j].id();
167  EcalTrigTowerDetId towid2 = myid2.tower();
168  if (towid1 == towid2) {
169  float theta =
170  theBarrelGeometry->getGeometry(myid2)->getPosition().theta();
171  Etsum += (*rechit_EB_col.product())[j].energy() * sin(theta);
172  }
173  j++;
174  if (count > 1800)
175  loopend = true;
176  }
177  mapTow_Et.insert(pair<EcalTrigTowerDetId, float>(towid1, Etsum));
178  }
179 
180  for (unsigned int i = 0; i < rechit_EE_col.product()->size(); i++) {
181  const EEDetId &myid1 = (*rechit_EE_col.product())[i].id();
182  EcalTrigTowerDetId towid1 = (*eTTmap_).towerOf(myid1);
183  float theta = theEndcapGeometry->getGeometry(myid1)->getPosition().theta();
184  float Etsum = (*rechit_EE_col.product())[i].energy() * sin(theta);
185  bool test_alreadyin = false;
186  map<EcalTrigTowerDetId, float>::iterator ittest = mapTow_Et.find(towid1);
187  if (ittest != mapTow_Et.end())
188  test_alreadyin = true;
189  if (test_alreadyin)
190  continue;
191  unsigned int j = i + 1;
192  bool loopend = false;
193  unsigned int count = 0;
194  while (j < rechit_EE_col.product()->size() && !loopend) {
195  const EEDetId &myid2 = (*rechit_EE_col.product())[j].id();
196  EcalTrigTowerDetId towid2 = (*eTTmap_).towerOf(myid2);
197  if (towid1 == towid2) {
198  float theta =
199  theEndcapGeometry->getGeometry(myid2)->getPosition().theta();
200  Etsum += (*rechit_EE_col.product())[j].energy() * sin(theta);
201  }
202  // else loopend=true;
203  j++;
204  if (count > 500)
205  loopend = true;
206  }
207  // alreadyin_EE.push_back(towid1);
208  mapTow_Et.insert(pair<EcalTrigTowerDetId, float>(towid1, Etsum));
209  }
210 
211  EcalTPGScale ecalScale;
212  ecalScale.setEventSetup(iSetup);
213  for (unsigned int i = 0; i < tp.product()->size(); i++) {
214  EcalTriggerPrimitiveDigi d = (*(tp.product()))[i];
215  const EcalTrigTowerDetId TPtowid = d.id();
216  map<EcalTrigTowerDetId, float>::iterator it = mapTow_Et.find(TPtowid);
217  float Et = ecalScale.getTPGInGeV(d.compressedEt(), TPtowid);
218  if (d.id().ietaAbs() == 27 || d.id().ietaAbs() == 28)
219  Et *= 2;
220  iphi_ = TPtowid.iphi();
221  ieta_ = TPtowid.ieta();
222  tpgADC_ = d.compressedEt();
223  tpgGeV_ = Et;
224  ttf_ = d.ttFlag();
225  fg_ = d.fineGrain();
226  if (it != mapTow_Et.end()) {
227  hTPvsRechit_->Fill(it->second, Et);
228  hTPoverRechit_->Fill(Et / it->second);
229  eRec_ = it->second;
230  }
231  tree_->Fill();
232  }
233 }
mps_fire.i
i
Definition: mps_fire.py:338
EcalTPGScale::setEventSetup
void setEventSetup(const edm::EventSetup &evtSetup)
Definition: EcalTPGScale.cc:19
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
theta
Geom::Theta< T > theta() const
Definition: Basic3DVectorLD.h:179
EcalTPGScale::getTPGInGeV
double getTPGInGeV(const EcalTriggerPrimitiveDigi &tpDigi)
Definition: EcalTPGScale.cc:24
EcalTrigTowerDetId::ieta
int ieta() const
get the tower ieta
Definition: EcalTrigTowerDetId.h:44
EcalTriggerPrimitiveDigi::compressedEt
int compressedEt() const
get the encoded/compressed Et of interesting sample
Definition: EcalTriggerPrimitiveDigi.cc:24
EcalTrigTowerDetId::ietaAbs
int ietaAbs() const
get the absolute value of the tower ieta
Definition: EcalTrigTowerDetId.h:36
EBDetId::tower
EcalTrigTowerDetId tower() const
get the HCAL/trigger iphi of this crystal
Definition: EBDetId.h:57
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:480
EcalTrigTowerDetId::iphi
int iphi() const
get the tower iphi
Definition: EcalTrigTowerDetId.h:52
EcalTriggerPrimitiveDigi::id
const EcalTrigTowerDetId & id() const
Definition: EcalTriggerPrimitiveDigi.h:35
edm::Handle::product
T const * product() const
Definition: Handle.h:74
CaloSubdetectorGeometry::getGeometry
virtual std::shared_ptr< const CaloCellGeometry > getGeometry(const DetId &id) const
Get the cell geometry of a given detector id. Should return false if not found.
Definition: CaloSubdetectorGeometry.cc:46
EcalTrigPrimAnalyzer::rechits_labelEE_
edm::InputTag rechits_labelEE_
Definition: EcalTrigPrimAnalyzer.h:63
EcalTrigPrimAnalyzer::rechits_labelEB_
edm::InputTag rechits_labelEB_
Definition: EcalTrigPrimAnalyzer.h:62
EcalTrigTowerDetId::subDet
EcalSubdetector subDet() const
get the subDetector associated to the Trigger Tower
Definition: EcalTrigTowerDetId.h:33
edm
HLT enums.
Definition: AlignableModifier.h:17
edm::SortedCollection::size
size_type size() const
Definition: SortedCollection.h:236
edm::EventSetup::get
T get() const
Definition: EventSetup.h:71
EcalTriggerPrimitiveDigi::fineGrain
bool fineGrain() const
get the fine-grain bit of interesting sample
Definition: EcalTriggerPrimitiveDigi.cc:34
EcalTriggerPrimitiveDigi::ttFlag
int ttFlag() const
get the Trigger tower Flag of interesting sample
Definition: EcalTriggerPrimitiveDigi.cc:43
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
void EcalTrigPrimAnalyzer::endJob ( void  )
overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 235 of file EcalTrigPrimAnalyzer.cc.

References mps_fire::i.

235  {
236  for (unsigned int i = 0; i < 2; ++i) {
237  ecal_et_[i]->Write();
238  ecal_tt_[i]->Write();
239  ecal_fgvb_[i]->Write();
240  }
241  if (recHits_) {
242  hTPvsRechit_->Write();
243  hTPoverRechit_->Write();
244  }
245 }
mps_fire.i
i
Definition: mps_fire.py:338

Member Data Documentation

TH1I* EcalTrigPrimAnalyzer::ecal_et_[2]
private

Definition at line 48 of file EcalTrigPrimAnalyzer.h.

TH1I* EcalTrigPrimAnalyzer::ecal_fgvb_[2]
private

Definition at line 50 of file EcalTrigPrimAnalyzer.h.

std::vector<std::string> EcalTrigPrimAnalyzer::ecal_parts_
private

Definition at line 47 of file EcalTrigPrimAnalyzer.h.

TH1I* EcalTrigPrimAnalyzer::ecal_tt_[2]
private

Definition at line 49 of file EcalTrigPrimAnalyzer.h.

float EcalTrigPrimAnalyzer::eRec_
private

Definition at line 58 of file EcalTrigPrimAnalyzer.h.

int EcalTrigPrimAnalyzer::fg_
private

Definition at line 57 of file EcalTrigPrimAnalyzer.h.

TH1I * EcalTrigPrimAnalyzer::histBar
private

Definition at line 51 of file EcalTrigPrimAnalyzer.h.

TH1I* EcalTrigPrimAnalyzer::histEndc
private

Definition at line 51 of file EcalTrigPrimAnalyzer.h.

TFile* EcalTrigPrimAnalyzer::histfile_
private

Definition at line 52 of file EcalTrigPrimAnalyzer.h.

TH1F* EcalTrigPrimAnalyzer::hTPoverRechit_
private

Definition at line 54 of file EcalTrigPrimAnalyzer.h.

TH2F* EcalTrigPrimAnalyzer::hTPvsRechit_
private

Definition at line 53 of file EcalTrigPrimAnalyzer.h.

int EcalTrigPrimAnalyzer::ieta_
private

Definition at line 57 of file EcalTrigPrimAnalyzer.h.

int EcalTrigPrimAnalyzer::iphi_
private

Definition at line 57 of file EcalTrigPrimAnalyzer.h.

edm::InputTag EcalTrigPrimAnalyzer::label_
private

Definition at line 60 of file EcalTrigPrimAnalyzer.h.

Referenced by SequenceTypes.DummyModule::__repr__(), Mixins._Labelable::_findDependencies(), and Mixins._Labelable::setLabel().

bool EcalTrigPrimAnalyzer::recHits_
private

Definition at line 65 of file EcalTrigPrimAnalyzer.h.

edm::InputTag EcalTrigPrimAnalyzer::rechits_labelEB_
private

Definition at line 62 of file EcalTrigPrimAnalyzer.h.

edm::InputTag EcalTrigPrimAnalyzer::rechits_labelEE_
private

Definition at line 63 of file EcalTrigPrimAnalyzer.h.

int EcalTrigPrimAnalyzer::tpgADC_
private

Definition at line 57 of file EcalTrigPrimAnalyzer.h.

float EcalTrigPrimAnalyzer::tpgGeV_
private

Definition at line 58 of file EcalTrigPrimAnalyzer.h.

TTree* EcalTrigPrimAnalyzer::tree_
private

Definition at line 55 of file EcalTrigPrimAnalyzer.h.

int EcalTrigPrimAnalyzer::ttf_
private

Definition at line 57 of file EcalTrigPrimAnalyzer.h.

Generated for CMSSW Reference Manual by doxygen 1.8.11