CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
List of all members | Public Member Functions | Private Member Functions | Private Attributes
EcalEBTrigPrimAnalyzer Class Reference

#include <EcalEBTrigPrimAnalyzer.h>

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

Public Member Functions

void analyze (const edm::Event &, const edm::EventSetup &) override
 
 EcalEBTrigPrimAnalyzer (const edm::ParameterSet &)
 
void endJob () override
 
 ~EcalEBTrigPrimAnalyzer () override
 
- Public Member Functions inherited from edm::one::EDAnalyzer<>
 EDAnalyzer ()=default
 
 EDAnalyzer (const EDAnalyzer &)=delete
 
SerialTaskQueueglobalLuminosityBlocksQueue () final
 
SerialTaskQueueglobalRunsQueue () final
 
const EDAnalyzeroperator= (const EDAnalyzer &)=delete
 
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< 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
 
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::array< std::vector< ModuleDescription const * > *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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
 
void selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Private Member Functions

int getIndex (const EBDigiCollection *, EcalTrigTowerDetId &id)
 
void init (const edm::EventSetup &)
 

Private Attributes

edm::ESGetToken
< CaloSubdetectorGeometry,
EcalBarrelGeometryRecord
barrelGeomToken_
 
bool debug_
 
TH1I * ecal_et_ [2]
 
TH1I * ecal_fgvb_ [2]
 
std::vector< std::string > ecal_parts_
 
TH1I * ecal_tt_ [2]
 
float eRec_
 
edm::ESHandle
< EcalTrigTowerConstituentsMap
eTTmap_
 
edm::ESGetToken
< EcalTrigTowerConstituentsMap,
IdealGeometryRecord
eTTmapToken_
 
int fg_
 
TH1F * hAllRechitEt_
 
TH1F * hAllTPperEvt_
 
TH1F * hDeltaEt_
 
TH1I * histBar
 
TH1I * histEndc
 
TFile * histfile_
 
TH1F * hRatioEt_
 
TH1F * hRechitEt_
 
TH2F * hRH_iphiVsieta_
 
TH2F * hRH_iphiVsieta_fullrange_
 
TH2F * hTP_iphiVsieta_
 
TH2F * hTP_iphiVsieta_fullrange_
 
TH1F * hTPEt_
 
TH1F * hTPoverRechit_
 
TH1F * hTPperEvt_
 
TH2F * hTPvsRechit_
 
TH2F * hTPvsTow_eta_
 
TH2F * hTPvsTow_ieta_
 
int nEvents_
 
edm::EDGetTokenT
< EcalEBTrigPrimDigiCollection
primToken_
 
bool recHits_
 
edm::EDGetTokenT
< EcalRecHitCollection
rechits_labelEB_
 
int rhIeta_
 
int rhIphi_
 
edm::EDGetTokenT
< EBDigiCollection
tokenEBdigi_
 
EcalTPGScale::Tokens tokens_
 
int tpgADC_
 
float tpgGeV_
 
int tpIeta_
 
int tpIphi_
 
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)
 
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B > consumes (edm::InputTag tag) noexcept
 
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<Transition Tr = Transition::Event>
constexpr auto esConsumes () noexcept
 
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...
 
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)
 
void resetItemsToGetFrom (BranchType iType)
 

Detailed Description

Definition at line 29 of file EcalEBTrigPrimAnalyzer.h.

Constructor & Destructor Documentation

EcalEBTrigPrimAnalyzer::EcalEBTrigPrimAnalyzer ( const edm::ParameterSet iConfig)
explicit

Definition at line 34 of file EcalEBTrigPrimAnalyzer.cc.

References barrelGeomToken_, debug_, ecal_et_, ecal_fgvb_, ecal_parts_, ecal_tt_, eRec_, eTTmapToken_, fg_, edm::ParameterSet::getParameter(), hAllRechitEt_, hAllTPperEvt_, hDeltaEt_, histfile_, hRatioEt_, hRechitEt_, hRH_iphiVsieta_, hRH_iphiVsieta_fullrange_, hTP_iphiVsieta_, hTP_iphiVsieta_fullrange_, hTPEt_, hTPoverRechit_, hTPperEvt_, hTPvsRechit_, hTPvsTow_eta_, hTPvsTow_ieta_, mps_fire::i, nEvents_, primToken_, recHits_, rechits_labelEB_, rhIeta_, rhIphi_, runGCPTkAlMap::title, tokenEBdigi_, tpgADC_, tpgGeV_, tpIeta_, tpIphi_, tree_, and ttf_.

36 
37 {
38  ecal_parts_.push_back("Barrel");
39 
40  histfile_ = new TFile("histos.root", "RECREATE");
41  tree_ = new TTree("TPGtree", "TPGtree");
42  tree_->Branch("tpIphi", &tpIphi_, "tpIphi/I");
43  tree_->Branch("tpIeta", &tpIeta_, "tpIeta/I");
44  tree_->Branch("rhIphi", &rhIphi_, "rhIphi/I");
45  tree_->Branch("rhIeta", &rhIeta_, "rhIeta/I");
46  tree_->Branch("eRec", &eRec_, "eRec/F");
47  tree_->Branch("tpgADC", &tpgADC_, "tpgADC/I");
48  tree_->Branch("tpgGeV", &tpgGeV_, "tpgGeV/F");
49  tree_->Branch("ttf", &ttf_, "ttf/I");
50  tree_->Branch("fg", &fg_, "fg/I");
51  for (unsigned int i = 0; i < ecal_parts_.size(); ++i) {
52  char title[30];
53  sprintf(title, "%s_Et", ecal_parts_[i].c_str());
54  ecal_et_[i] = new TH1I(title, "Et", 255, 0, 255);
55  sprintf(title, "%s_ttf", ecal_parts_[i].c_str());
56  ecal_tt_[i] = new TH1I(title, "TTF", 10, 0, 10);
57  sprintf(title, "%s_fgvb", ecal_parts_[i].c_str());
58  ecal_fgvb_[i] = new TH1I(title, "FGVB", 10, 0, 10);
59  }
60  eTTmapToken_ = esConsumes<EcalTrigTowerConstituentsMap, IdealGeometryRecord>();
61  recHits_ = iConfig.getParameter<bool>("AnalyzeRecHits");
62  debug_ = iConfig.getParameter<bool>("Debug");
63  rechits_labelEB_ = consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("inputRecHitsEB"));
64  primToken_ = consumes<EcalEBTrigPrimDigiCollection>(iConfig.getParameter<edm::InputTag>("inputTP"));
65  barrelGeomToken_ = esConsumes<CaloSubdetectorGeometry, EcalBarrelGeometryRecord>(edm::ESInputTag("", "EcalBarrel"));
66  tokenEBdigi_ = consumes<EBDigiCollection>(iConfig.getParameter<edm::InputTag>("barrelEcalDigis"));
67  nEvents_ = 0;
68 
69  hTPvsTow_eta_ = new TH2F("TP_vs_Tow_eta", "TP vs Tow eta ; #eta(tow); #eta(tp)", 50, -2.5, 2.5, 50, -2.5, 2.5);
70  hAllTPperEvt_ = new TH1F("AllTPperEvt", "TP per Event; N_{TP}; ", 100, 0., 20000.);
71  hTPperEvt_ = new TH1F("TPperEvt", "N_{TP} per Event; N_{TP}; ", 100, 0., 500.);
72  hTP_iphiVsieta_ = new TH2F("TP_iphiVsieta", "TP i#phi vs i#eta ; i#eta(tp); i#phi(tp)", 10, 70, 80, 10, 340, 350);
74  new TH2F("TP_iphiVsieta_fullrange", "TP i#phi vs i#eta ; i#eta(tp); i#phi(tp)", 200, -100, 100, 350, 0, 350);
75 
76  if (recHits_) {
78  new TH2F("TP_vs_Tow_ieta", "TP vs Tow ieta ; i#eta(tow); i#eta(tp)", 200, -100, 100, 200, -100, 100);
79 
80  hTPvsRechit_ = new TH2F("TP_vs_RecHit", "TP vs rechit Et;E_{T}(rh) (GeV);E_{T}(tp) (GeV)", 100, 0, 50, 100, 0, 50);
81  hDeltaEt_ = new TH1F("DeltaEt", "[Et(rh)-Et(TP)]/Et(rh); [E_{T}(rh)-E_{T}(tp)]/E_{T}(rh); Counts", 200, -1, 1);
82  hTPoverRechit_ = new TH1F("TP_over_RecHit", "Et(TP/rechit); E_{T}(tp)/E_{T}(rh); Counts", 200, 0, 2);
83  hRechitEt_ = new TH1F("RecHitEt", "E_{T};E_{T}(rh) (GeV);Counts", 100, 0, 50);
84  hTPEt_ = new TH1F("TPEt", "E_{T}{tp);E_{T}(rh) (GeV);Count", 100, 0, 50);
85  hRatioEt_ = new TH1F("RatioTPoverRH", "Et", 100, 0, 50);
86  hAllRechitEt_ = new TH1F("AllRecHit", "Et", 100, 0, 50);
87 
88  hRH_iphiVsieta_ = new TH2F("RH_iphiVsieta", "RH i#phi vs i#eta ; i#eta(rh); i#phi(rh)", 10, 70, 80, 10, 340, 350);
90  new TH2F("RH_iphiVsieta_fullrange", "RH i#phi vs i#eta ; i#eta(rh); i#phi(rh)", 200, -100, 100, 350, 0, 350);
91  }
92 }
edm::EDGetTokenT< EcalRecHitCollection > rechits_labelEB_
edm::ESGetToken< CaloSubdetectorGeometry, EcalBarrelGeometryRecord > barrelGeomToken_
EcalTPGScale::Tokens tokens_
ConsumesCollector consumesCollector()
Use a ConsumesCollector to gather consumes information from helper functions.
edm::EDGetTokenT< EBDigiCollection > tokenEBdigi_
edm::ESGetToken< EcalTrigTowerConstituentsMap, IdealGeometryRecord > eTTmapToken_
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
std::vector< std::string > ecal_parts_
edm::EDGetTokenT< EcalEBTrigPrimDigiCollection > primToken_
EcalEBTrigPrimAnalyzer::~EcalEBTrigPrimAnalyzer ( )
override

Definition at line 94 of file EcalEBTrigPrimAnalyzer.cc.

References histfile_.

94  {
95  // do anything here that needs to be done at desctruction time
96  // (e.g. close files, deallocate resources etc.)
97 
98  histfile_->Write();
99  histfile_->Close();
100 }

Member Function Documentation

void EcalEBTrigPrimAnalyzer::analyze ( const edm::Event iEvent,
const edm::EventSetup iSetup 
)
overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 109 of file EcalEBTrigPrimAnalyzer.cc.

References barrelGeomToken_, gather_cfg::cout, ztail::d, debug_, ecal_et_, EcalEBTriggerPrimitiveDigi::encodedEt(), relval_parameters_module::energy, eRec_, edm::Event::getByToken(), edm::EventSetup::getHandle(), hAllRechitEt_, hAllTPperEvt_, hDeltaEt_, hRatioEt_, hRechitEt_, hRH_iphiVsieta_, hRH_iphiVsieta_fullrange_, hTP_iphiVsieta_, hTP_iphiVsieta_fullrange_, hTPEt_, hTPoverRechit_, hTPperEvt_, hTPvsRechit_, mps_fire::i, EcalEBTriggerPrimitiveDigi::id(), edm::EventBase::id(), EBDetId::ieta(), init(), EBDetId::iphi(), dqmiolumiharvest::j, nEvents_, primToken_, recHits_, rechits_labelEB_, rhIeta_, rhIphi_, EcalEBTriggerPrimitiveDigi::sample(), funct::sin(), EcalEBTriggerPrimitiveDigi::size(), theta(), tokens_, cmsswSequenceInfo::tp, tpgADC_, tpgGeV_, tpIeta_, tpIphi_, and tree_.

109  {
110  using namespace edm;
111  using namespace std;
112  nEvents_++;
113 
114  if (nEvents_ == 1)
115  this->init(iSetup);
116 
117  // Get input
119  iEvent.getByToken(primToken_, tp);
120  //
121  /*
122  edm::Handle<EBDigiCollection> barrelDigiHandle;
123  const EBDigiCollection *ebdigi=NULL;
124  iEvent.getByToken(tokenEBdigi_,barrelDigiHandle);
125  ebdigi=barrelDigiHandle.product();
126  */
127 
128  for (unsigned int i = 0; i < tp.product()->size(); i++) {
129  EcalEBTriggerPrimitiveDigi d = (*(tp.product()))[i];
130  int subdet = 0;
131  if (subdet == 0) {
132  ecal_et_[subdet]->Fill(d.encodedEt());
133  }
134  }
135 
136  // if (!recHits_) return;
137 
138  edm::Handle<EcalRecHitCollection> rechit_EB_col;
139  if (recHits_) {
140  // get the RecHits
141  iEvent.getByToken(rechits_labelEB_, rechit_EB_col);
142  }
143 
144  edm::ESHandle<CaloSubdetectorGeometry> theBarrelGeometry_handle = iSetup.getHandle(barrelGeomToken_);
145  const CaloSubdetectorGeometry* theBarrelGeometry = theBarrelGeometry_handle.product();
146 
147  if (debug_) {
148  std::cout << " TP analyzer =================> Treating event " << iEvent.id() << " Number of TPs "
149  << tp.product()->size() << std::endl;
150  if (recHits_)
151  std::cout << " Number of EB rechits " << rechit_EB_col.product()->size() << std::endl;
152  }
153  hAllTPperEvt_->Fill(float(tp.product()->size()));
154 
155  //if ( iEvent.id().event() != 648) return;
156 
157  //EcalEBTPGScale ecalScale ;
158  EcalTPGScale ecalScale(tokens_, iSetup);
159 
160  // for(unsigned int iDigi = 0; iDigi < ebdigi->size() ; ++iDigi) {
161  // EBDataFrame myFrame((*ebdigi)[iDigi]);
162  // const EBDetId & myId = myFrame.id();
163 
164  int nTP = 0;
165  for (unsigned int i = 0; i < tp.product()->size(); i++) {
166  EcalEBTriggerPrimitiveDigi d = (*(tp.product()))[i];
167  const EBDetId TPid = d.id();
168  // if ( myId != TPid ) continue;
169 
170  /*
171  int index=getIndex(ebdigi,coarser);
172  std::cout << " Same xTal " << myId << " " << TPid << " coarser " << coarser << " index " << index << std::endl;
173  double Et = ecalScale.getTPGInGeV(d.encodedEt(), coarser) ;
174  */
175  //this works if the energy is compressed into 8 bits float Et=d.compressedEt()/2.; // 2ADC counts/GeV
176  float Et = d.encodedEt() / 8.; // 8 ADCcounts/GeV
177  if (Et <= 5)
178  continue;
179  //if ( Et<= 0 ) continue;
180  nTP++;
181 
182  std::cout << " TP digi size " << d.size() << std::endl;
183  for (int iBx = 0; iBx < d.size(); iBx++) {
184  std::cout << " TP samples " << d.sample(iBx) << std::endl;
185  }
186 
187  // EcalTrigTowerDetId coarser=(*eTTmap_).towerOf(myId);
188  // does not work float etaTow = theBarrelGeometry->getGeometry(coarser)->getPosition().theta();
189  // float etaTP = theBarrelGeometry->getGeometry(TPid)->getPosition().eta();
190  // does not work hTPvsTow_eta_->Fill ( etaTow, etaTP );
191  // hTPvsTow_ieta_->Fill ( coarser.ieta(), TPid.ieta() );
192 
193  tpIphi_ = TPid.iphi();
194  tpIeta_ = TPid.ieta();
195  tpgADC_ = d.encodedEt();
196  tpgGeV_ = Et;
197 
198  hTP_iphiVsieta_->Fill(TPid.ieta(), TPid.iphi(), Et);
199  hTP_iphiVsieta_fullrange_->Fill(TPid.ieta(), TPid.iphi(), Et);
200 
201  if (recHits_) {
202  for (unsigned int j = 0; j < rechit_EB_col.product()->size(); j++) {
203  const EBDetId& myid1 = (*rechit_EB_col.product())[j].id();
204  float theta = theBarrelGeometry->getGeometry(myid1)->getPosition().theta();
205  float rhEt = ((*rechit_EB_col.product())[j].energy()) * sin(theta);
206  if (myid1 == TPid) {
207  if (debug_)
208  std::cout << " Analyzer same cristal " << myid1 << " " << TPid << std::endl;
209  // if ( rhEt < 1.5 && Et > 10 ) {
210  // std::cout << " TP analyzer =================> Treating event "<<iEvent.id()<< ", Number of EB rechits "<< rechit_EB_col.product()->size() << " Number of TPs " << tp.product()->size() << std::endl;
211  //std::cout << " TP compressed et " << d.encodedEt() << " Et in GeV " << Et << " RH Et " << rhEt << " Et/rhEt " << Et/rhEt << std::endl;
212  //}
213 
214  //std::cout << " TP out " << d << std::endl;
215 
216  // for (int isam=0;isam< d.size();++isam) {
217  // std::cout << " d[isam].raw() " << d[isam].raw() << std::endl;
218  //}
219 
220  rhIphi_ = myid1.iphi();
221  rhIeta_ = myid1.ieta();
222  hRH_iphiVsieta_->Fill(myid1.ieta(), myid1.iphi(), rhEt);
223  hRH_iphiVsieta_fullrange_->Fill(myid1.ieta(), myid1.iphi(), rhEt);
224 
225  hTPvsRechit_->Fill(rhEt, Et);
226  hTPoverRechit_->Fill(Et / rhEt);
227  hDeltaEt_->Fill((rhEt - Et) / rhEt);
228  if (debug_)
229  std::cout << " TP compressed et " << d.encodedEt() << " Et in GeV " << Et << " RH Et " << rhEt
230  << " Et/rhEt " << Et / rhEt << std::endl;
231  hRechitEt_->Fill(rhEt);
232  hTPEt_->Fill(Et);
233  if (rhEt < 1000000)
234  eRec_ = rhEt;
235  tree_->Fill();
236  }
237 
238  } // end loop of recHits
239  } // if recHits
240 
241  } // end loop over TP collection
242 
243  // } // end loop over digi collection
244 
245  hTPperEvt_->Fill(float(nTP));
246 
247  if (recHits_) {
248  hRatioEt_->Divide(hTPEt_, hRechitEt_);
249  for (unsigned int j = 0; j < rechit_EB_col.product()->size(); j++) {
250  const EBDetId& myid1 = (*rechit_EB_col.product())[j].id();
251  float theta = theBarrelGeometry->getGeometry(myid1)->getPosition().theta();
252  float rhEt = ((*rechit_EB_col.product())[j].energy()) * sin(theta);
253  if (rhEt > 0)
254  hAllRechitEt_->Fill(rhEt);
255  }
256  }
257 }
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:539
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
Geom::Theta< T > theta() const
edm::EDGetTokenT< EcalRecHitCollection > rechits_labelEB_
edm::ESGetToken< CaloSubdetectorGeometry, EcalBarrelGeometryRecord > barrelGeomToken_
int iphi() const
get the crystal iphi
Definition: EBDetId.h:51
tuple d
Definition: ztail.py:151
EcalTPGScale::Tokens tokens_
const EcalEBTriggerPrimitiveSample & sample(int i) const
int encodedEt() const
get the 10 bits Et of interesting sample
int ieta() const
get the crystal ieta
Definition: EBDetId.h:49
edm::EventID id() const
Definition: EventBase.h:59
tuple cout
Definition: gather_cfg.py:144
edm::EDGetTokenT< EcalEBTrigPrimDigiCollection > primToken_
ESHandle< T > getHandle(const ESGetToken< T, R > &iToken) const
Definition: EventSetup.h:157
void init(const edm::EventSetup &)
void EcalEBTrigPrimAnalyzer::endJob ( void  )
overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 259 of file EcalEBTrigPrimAnalyzer.cc.

References ecal_et_, ecal_fgvb_, ecal_parts_, ecal_tt_, hAllRechitEt_, hAllTPperEvt_, hDeltaEt_, hRatioEt_, hRechitEt_, hRH_iphiVsieta_, hRH_iphiVsieta_fullrange_, hTP_iphiVsieta_, hTP_iphiVsieta_fullrange_, hTPEt_, hTPoverRechit_, hTPperEvt_, hTPvsRechit_, hTPvsTow_eta_, hTPvsTow_ieta_, mps_fire::i, and recHits_.

259  {
260  for (unsigned int i = 0; i < ecal_parts_.size(); ++i) {
261  ecal_et_[i]->Write();
262  ecal_tt_[i]->Write();
263  ecal_fgvb_[i]->Write();
264  }
265 
266  hAllTPperEvt_->Write();
267  hTPperEvt_->Write();
268 
269  if (recHits_) {
270  hTPvsTow_ieta_->Write();
271  hTPvsTow_eta_->Write();
272  hTPvsRechit_->Write();
273  hTPoverRechit_->Write();
274  hAllRechitEt_->Write();
275  hRechitEt_->Write();
276  hDeltaEt_->Write();
277  hTPEt_->Write();
278  hRatioEt_->Write();
279  hTP_iphiVsieta_->Write();
280  hRH_iphiVsieta_->Write();
281  hTP_iphiVsieta_fullrange_->Write();
282  hRH_iphiVsieta_fullrange_->Write();
283  }
284 }
std::vector< std::string > ecal_parts_
int EcalEBTrigPrimAnalyzer::getIndex ( const EBDigiCollection ,
EcalTrigTowerDetId id 
)
inlineprivate

Definition at line 81 of file EcalEBTrigPrimAnalyzer.h.

81 { return id.hashedIndex(); }
void EcalEBTrigPrimAnalyzer::init ( const edm::EventSetup iSetup)
private

Definition at line 102 of file EcalEBTrigPrimAnalyzer.cc.

References eTTmap_, eTTmapToken_, and edm::EventSetup::getHandle().

Referenced by analyze().

102 { eTTmap_ = iSetup.getHandle(eTTmapToken_); }
edm::ESGetToken< EcalTrigTowerConstituentsMap, IdealGeometryRecord > eTTmapToken_
edm::ESHandle< EcalTrigTowerConstituentsMap > eTTmap_
ESHandle< T > getHandle(const ESGetToken< T, R > &iToken) const
Definition: EventSetup.h:157

Member Data Documentation

edm::ESGetToken<CaloSubdetectorGeometry, EcalBarrelGeometryRecord> EcalEBTrigPrimAnalyzer::barrelGeomToken_
private

Definition at line 78 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalEBTrigPrimAnalyzer().

bool EcalEBTrigPrimAnalyzer::debug_
private

Definition at line 75 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalEBTrigPrimAnalyzer().

TH1I* EcalEBTrigPrimAnalyzer::ecal_et_[2]
private

Definition at line 43 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), EcalEBTrigPrimAnalyzer(), and endJob().

TH1I* EcalEBTrigPrimAnalyzer::ecal_fgvb_[2]
private

Definition at line 45 of file EcalEBTrigPrimAnalyzer.h.

Referenced by EcalEBTrigPrimAnalyzer(), and endJob().

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

Definition at line 42 of file EcalEBTrigPrimAnalyzer.h.

Referenced by EcalEBTrigPrimAnalyzer(), and endJob().

TH1I* EcalEBTrigPrimAnalyzer::ecal_tt_[2]
private

Definition at line 44 of file EcalEBTrigPrimAnalyzer.h.

Referenced by EcalEBTrigPrimAnalyzer(), and endJob().

float EcalEBTrigPrimAnalyzer::eRec_
private

Definition at line 68 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalEBTrigPrimAnalyzer().

edm::ESHandle<EcalTrigTowerConstituentsMap> EcalEBTrigPrimAnalyzer::eTTmap_
private

Definition at line 76 of file EcalEBTrigPrimAnalyzer.h.

Referenced by init().

edm::ESGetToken<EcalTrigTowerConstituentsMap, IdealGeometryRecord> EcalEBTrigPrimAnalyzer::eTTmapToken_
private

Definition at line 77 of file EcalEBTrigPrimAnalyzer.h.

Referenced by EcalEBTrigPrimAnalyzer(), and init().

int EcalEBTrigPrimAnalyzer::fg_
private

Definition at line 66 of file EcalEBTrigPrimAnalyzer.h.

Referenced by EcalEBTrigPrimAnalyzer().

TH1F* EcalEBTrigPrimAnalyzer::hAllRechitEt_
private

Definition at line 53 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), EcalEBTrigPrimAnalyzer(), and endJob().

TH1F* EcalEBTrigPrimAnalyzer::hAllTPperEvt_
private

Definition at line 48 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), EcalEBTrigPrimAnalyzer(), and endJob().

TH1F* EcalEBTrigPrimAnalyzer::hDeltaEt_
private

Definition at line 52 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), EcalEBTrigPrimAnalyzer(), and endJob().

TH1I * EcalEBTrigPrimAnalyzer::histBar
private

Definition at line 46 of file EcalEBTrigPrimAnalyzer.h.

TH1I* EcalEBTrigPrimAnalyzer::histEndc
private

Definition at line 46 of file EcalEBTrigPrimAnalyzer.h.

TFile* EcalEBTrigPrimAnalyzer::histfile_
private

Definition at line 47 of file EcalEBTrigPrimAnalyzer.h.

Referenced by EcalEBTrigPrimAnalyzer(), and ~EcalEBTrigPrimAnalyzer().

TH1F* EcalEBTrigPrimAnalyzer::hRatioEt_
private

Definition at line 56 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), EcalEBTrigPrimAnalyzer(), and endJob().

TH1F* EcalEBTrigPrimAnalyzer::hRechitEt_
private

Definition at line 54 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), EcalEBTrigPrimAnalyzer(), and endJob().

TH2F* EcalEBTrigPrimAnalyzer::hRH_iphiVsieta_
private

Definition at line 60 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), EcalEBTrigPrimAnalyzer(), and endJob().

TH2F* EcalEBTrigPrimAnalyzer::hRH_iphiVsieta_fullrange_
private

Definition at line 62 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), EcalEBTrigPrimAnalyzer(), and endJob().

TH2F* EcalEBTrigPrimAnalyzer::hTP_iphiVsieta_
private

Definition at line 59 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), EcalEBTrigPrimAnalyzer(), and endJob().

TH2F* EcalEBTrigPrimAnalyzer::hTP_iphiVsieta_fullrange_
private

Definition at line 61 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), EcalEBTrigPrimAnalyzer(), and endJob().

TH1F* EcalEBTrigPrimAnalyzer::hTPEt_
private

Definition at line 55 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), EcalEBTrigPrimAnalyzer(), and endJob().

TH1F* EcalEBTrigPrimAnalyzer::hTPoverRechit_
private

Definition at line 51 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), EcalEBTrigPrimAnalyzer(), and endJob().

TH1F* EcalEBTrigPrimAnalyzer::hTPperEvt_
private

Definition at line 49 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), EcalEBTrigPrimAnalyzer(), and endJob().

TH2F* EcalEBTrigPrimAnalyzer::hTPvsRechit_
private

Definition at line 50 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), EcalEBTrigPrimAnalyzer(), and endJob().

TH2F* EcalEBTrigPrimAnalyzer::hTPvsTow_eta_
private

Definition at line 58 of file EcalEBTrigPrimAnalyzer.h.

Referenced by EcalEBTrigPrimAnalyzer(), and endJob().

TH2F* EcalEBTrigPrimAnalyzer::hTPvsTow_ieta_
private

Definition at line 57 of file EcalEBTrigPrimAnalyzer.h.

Referenced by EcalEBTrigPrimAnalyzer(), and endJob().

int EcalEBTrigPrimAnalyzer::nEvents_
private

Definition at line 38 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalEBTrigPrimAnalyzer().

edm::EDGetTokenT<EcalEBTrigPrimDigiCollection> EcalEBTrigPrimAnalyzer::primToken_
private

Definition at line 71 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalEBTrigPrimAnalyzer().

bool EcalEBTrigPrimAnalyzer::recHits_
private

Definition at line 74 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), EcalEBTrigPrimAnalyzer(), and endJob().

edm::EDGetTokenT<EcalRecHitCollection> EcalEBTrigPrimAnalyzer::rechits_labelEB_
private

Definition at line 72 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalEBTrigPrimAnalyzer().

int EcalEBTrigPrimAnalyzer::rhIeta_
private

Definition at line 67 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalEBTrigPrimAnalyzer().

int EcalEBTrigPrimAnalyzer::rhIphi_
private

Definition at line 67 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalEBTrigPrimAnalyzer().

edm::EDGetTokenT<EBDigiCollection> EcalEBTrigPrimAnalyzer::tokenEBdigi_
private

Definition at line 73 of file EcalEBTrigPrimAnalyzer.h.

Referenced by EcalEBTrigPrimAnalyzer().

EcalTPGScale::Tokens EcalEBTrigPrimAnalyzer::tokens_
private

Definition at line 79 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze().

int EcalEBTrigPrimAnalyzer::tpgADC_
private

Definition at line 66 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalEBTrigPrimAnalyzer().

float EcalEBTrigPrimAnalyzer::tpgGeV_
private

Definition at line 68 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalEBTrigPrimAnalyzer().

int EcalEBTrigPrimAnalyzer::tpIeta_
private

Definition at line 66 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalEBTrigPrimAnalyzer().

int EcalEBTrigPrimAnalyzer::tpIphi_
private

Definition at line 66 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalEBTrigPrimAnalyzer().

TTree* EcalEBTrigPrimAnalyzer::tree_
private

Definition at line 64 of file EcalEBTrigPrimAnalyzer.h.

Referenced by analyze(), and EcalEBTrigPrimAnalyzer().

int EcalEBTrigPrimAnalyzer::ttf_
private

Definition at line 66 of file EcalEBTrigPrimAnalyzer.h.

Referenced by EcalEBTrigPrimAnalyzer().