CMS 3D CMS Logo
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
 
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 Member Functions

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

Private Attributes

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_
 
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_
 
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)
 
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 24 of file EcalEBTrigPrimAnalyzer.h.

Constructor & Destructor Documentation

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

Definition at line 35 of file EcalEBTrigPrimAnalyzer.cc.

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

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 
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  tokenEBdigi_=consumes<EBDigiCollection>(iConfig.getParameter<edm::InputTag>("barrelEcalDigis"));
66  nEvents_=0;
67 
68  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);
69  hAllTPperEvt_ = new TH1F("AllTPperEvt","TP per Event; N_{TP}; ", 100, 0., 20000.);
70  hTPperEvt_ = new TH1F("TPperEvt","N_{TP} per Event; N_{TP}; ", 100, 0., 500.);
71  hTP_iphiVsieta_= new TH2F("TP_iphiVsieta","TP i#phi vs i#eta ; i#eta(tp); i#phi(tp)",10,70,80,10,340,350);
72  hTP_iphiVsieta_fullrange_= new TH2F("TP_iphiVsieta_fullrange","TP i#phi vs i#eta ; i#eta(tp); i#phi(tp)",200,-100,100,350,0,350);
73 
74  if (recHits_) {
75  hTPvsTow_ieta_= new TH2F("TP_vs_Tow_ieta","TP vs Tow ieta ; i#eta(tow); i#eta(tp)",200,-100,100,200,-100,100);
76 
77  hTPvsRechit_= new TH2F("TP_vs_RecHit","TP vs rechit Et;E_{T}(rh) (GeV);E_{T}(tp) (GeV)",100,0,50,100,0,50);
78  hDeltaEt_ = new TH1F("DeltaEt","[Et(rh)-Et(TP)]/Et(rh); [E_{T}(rh)-E_{T}(tp)]/E_{T}(rh); Counts",200,-1,1);
79  hTPoverRechit_= new TH1F("TP_over_RecHit","Et(TP/rechit); E_{T}(tp)/E_{T}(rh); Counts",200,0,2);
80  hRechitEt_= new TH1F("RecHitEt","E_{T};E_{T}(rh) (GeV);Counts",100,0,50);
81  hTPEt_= new TH1F("TPEt","E_{T}{tp);E_{T}(rh) (GeV);Count",100,0,50);
82  hRatioEt_ = new TH1F("RatioTPoverRH","Et",100,0,50);
83  hAllRechitEt_= new TH1F("AllRecHit","Et",100,0,50);
84 
85  hRH_iphiVsieta_= new TH2F("RH_iphiVsieta","RH i#phi vs i#eta ; i#eta(rh); i#phi(rh)",10,70,80,10,340,350);
86  hRH_iphiVsieta_fullrange_= new TH2F("RH_iphiVsieta_fullrange","RH i#phi vs i#eta ; i#eta(rh); i#phi(rh)",200,-100,100,350,0,350);
87 
88 
89  }
90 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
mps_fire.i
i
Definition: mps_fire.py:341
EcalEBTrigPrimAnalyzer::rechits_labelEB_
edm::EDGetTokenT< EcalRecHitCollection > rechits_labelEB_
Definition: EcalEBTrigPrimAnalyzer.h:68
EcalEBTrigPrimAnalyzer::tokenEBdigi_
edm::EDGetTokenT< EBDigiCollection > tokenEBdigi_
Definition: EcalEBTrigPrimAnalyzer.h:69
overlapproblemtsosanalyzer_cfi.title
title
Definition: overlapproblemtsosanalyzer_cfi.py:7
EcalEBTrigPrimAnalyzer::ecal_parts_
std::vector< std::string > ecal_parts_
Definition: EcalEBTrigPrimAnalyzer.h:38
EcalEBTrigPrimAnalyzer::primToken_
edm::EDGetTokenT< EcalEBTrigPrimDigiCollection > primToken_
Definition: EcalEBTrigPrimAnalyzer.h:67
EcalEBTrigPrimAnalyzer::~EcalEBTrigPrimAnalyzer ( )
override

Definition at line 93 of file EcalEBTrigPrimAnalyzer.cc.

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

Member Function Documentation

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

Definition at line 114 of file EcalEBTrigPrimAnalyzer.cc.

References gather_cfg::cout, ztail::d, EcalEBTriggerPrimitiveDigi::encodedEt(), HCALHighEnergyHPDFilter_cfi::energy, edm::EventSetup::get(), edm::Event::getByToken(), CaloSubdetectorGeometry::getGeometry(), mps_fire::i, EcalEBTriggerPrimitiveDigi::id(), edm::EventBase::id(), EBDetId::ieta(), init, EBDetId::iphi(), dqmiolumiharvest::j, edm::Handle< T >::product(), edm::ESHandle< T >::product(), EcalEBTriggerPrimitiveDigi::sample(), EcalTPGScale::setEventSetup(), funct::sin(), EcalEBTriggerPrimitiveDigi::size(), edm::SortedCollection< T, SORT >::size(), and theta().

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

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 269 of file EcalEBTrigPrimAnalyzer.cc.

References mps_fire::i.

269  {
270  for (unsigned int i=0;i<ecal_parts_.size();++i) {
271  ecal_et_[i]->Write();
272  ecal_tt_[i]->Write();
273  ecal_fgvb_[i]->Write();
274  }
275 
276 
277  hAllTPperEvt_->Write();
278  hTPperEvt_->Write();
279 
280  if (recHits_) {
281  hTPvsTow_ieta_->Write();
282  hTPvsTow_eta_->Write();
283  hTPvsRechit_->Write();
284  hTPoverRechit_->Write();
285  hAllRechitEt_->Write();
286  hRechitEt_->Write();
287  hDeltaEt_ ->Write();
288  hTPEt_->Write();
289  hRatioEt_->Write();
290  hTP_iphiVsieta_->Write();
291  hRH_iphiVsieta_->Write();
292  hTP_iphiVsieta_fullrange_->Write();
293  hRH_iphiVsieta_fullrange_->Write();
294 
295  }
296 }
mps_fire.i
i
Definition: mps_fire.py:341
EcalEBTrigPrimAnalyzer::ecal_parts_
std::vector< std::string > ecal_parts_
Definition: EcalEBTrigPrimAnalyzer.h:38
int EcalEBTrigPrimAnalyzer::getIndex ( const EBDigiCollection *  ,
EcalTrigTowerDetId &  id 
)
inlineprivate

Definition at line 74 of file EcalEBTrigPrimAnalyzer.h.

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

Definition at line 104 of file EcalEBTrigPrimAnalyzer.cc.

References edm::EventSetup::get().

104  {
105  iSetup.get<IdealGeometryRecord>().get(eTTmap_);
106 }
edm::EventSetup::get
T get() const
Definition: EventSetup.h:73
EcalEBTrigPrimAnalyzer::eTTmap_
edm::ESHandle< EcalTrigTowerConstituentsMap > eTTmap_
Definition: EcalEBTrigPrimAnalyzer.h:72

Member Data Documentation

bool EcalEBTrigPrimAnalyzer::debug_
private

Definition at line 71 of file EcalEBTrigPrimAnalyzer.h.

TH1I* EcalEBTrigPrimAnalyzer::ecal_et_[2]
private

Definition at line 39 of file EcalEBTrigPrimAnalyzer.h.

TH1I* EcalEBTrigPrimAnalyzer::ecal_fgvb_[2]
private

Definition at line 41 of file EcalEBTrigPrimAnalyzer.h.

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

Definition at line 38 of file EcalEBTrigPrimAnalyzer.h.

TH1I* EcalEBTrigPrimAnalyzer::ecal_tt_[2]
private

Definition at line 40 of file EcalEBTrigPrimAnalyzer.h.

float EcalEBTrigPrimAnalyzer::eRec_
private

Definition at line 64 of file EcalEBTrigPrimAnalyzer.h.

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

Definition at line 72 of file EcalEBTrigPrimAnalyzer.h.

int EcalEBTrigPrimAnalyzer::fg_
private

Definition at line 62 of file EcalEBTrigPrimAnalyzer.h.

TH1F* EcalEBTrigPrimAnalyzer::hAllRechitEt_
private

Definition at line 49 of file EcalEBTrigPrimAnalyzer.h.

TH1F* EcalEBTrigPrimAnalyzer::hAllTPperEvt_
private

Definition at line 44 of file EcalEBTrigPrimAnalyzer.h.

TH1F* EcalEBTrigPrimAnalyzer::hDeltaEt_
private

Definition at line 48 of file EcalEBTrigPrimAnalyzer.h.

TH1I * EcalEBTrigPrimAnalyzer::histBar
private

Definition at line 42 of file EcalEBTrigPrimAnalyzer.h.

TH1I* EcalEBTrigPrimAnalyzer::histEndc
private

Definition at line 42 of file EcalEBTrigPrimAnalyzer.h.

TFile* EcalEBTrigPrimAnalyzer::histfile_
private

Definition at line 43 of file EcalEBTrigPrimAnalyzer.h.

TH1F* EcalEBTrigPrimAnalyzer::hRatioEt_
private

Definition at line 52 of file EcalEBTrigPrimAnalyzer.h.

TH1F* EcalEBTrigPrimAnalyzer::hRechitEt_
private

Definition at line 50 of file EcalEBTrigPrimAnalyzer.h.

TH2F* EcalEBTrigPrimAnalyzer::hRH_iphiVsieta_
private

Definition at line 56 of file EcalEBTrigPrimAnalyzer.h.

TH2F* EcalEBTrigPrimAnalyzer::hRH_iphiVsieta_fullrange_
private

Definition at line 58 of file EcalEBTrigPrimAnalyzer.h.

TH2F* EcalEBTrigPrimAnalyzer::hTP_iphiVsieta_
private

Definition at line 55 of file EcalEBTrigPrimAnalyzer.h.

TH2F* EcalEBTrigPrimAnalyzer::hTP_iphiVsieta_fullrange_
private

Definition at line 57 of file EcalEBTrigPrimAnalyzer.h.

TH1F* EcalEBTrigPrimAnalyzer::hTPEt_
private

Definition at line 51 of file EcalEBTrigPrimAnalyzer.h.

TH1F* EcalEBTrigPrimAnalyzer::hTPoverRechit_
private

Definition at line 47 of file EcalEBTrigPrimAnalyzer.h.

TH1F* EcalEBTrigPrimAnalyzer::hTPperEvt_
private

Definition at line 45 of file EcalEBTrigPrimAnalyzer.h.

TH2F* EcalEBTrigPrimAnalyzer::hTPvsRechit_
private

Definition at line 46 of file EcalEBTrigPrimAnalyzer.h.

TH2F* EcalEBTrigPrimAnalyzer::hTPvsTow_eta_
private

Definition at line 54 of file EcalEBTrigPrimAnalyzer.h.

TH2F* EcalEBTrigPrimAnalyzer::hTPvsTow_ieta_
private

Definition at line 53 of file EcalEBTrigPrimAnalyzer.h.

int EcalEBTrigPrimAnalyzer::nEvents_
private

Definition at line 34 of file EcalEBTrigPrimAnalyzer.h.

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

Definition at line 67 of file EcalEBTrigPrimAnalyzer.h.

bool EcalEBTrigPrimAnalyzer::recHits_
private

Definition at line 70 of file EcalEBTrigPrimAnalyzer.h.

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

Definition at line 68 of file EcalEBTrigPrimAnalyzer.h.

int EcalEBTrigPrimAnalyzer::rhIeta_
private

Definition at line 63 of file EcalEBTrigPrimAnalyzer.h.

int EcalEBTrigPrimAnalyzer::rhIphi_
private

Definition at line 63 of file EcalEBTrigPrimAnalyzer.h.

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

Definition at line 69 of file EcalEBTrigPrimAnalyzer.h.

int EcalEBTrigPrimAnalyzer::tpgADC_
private

Definition at line 62 of file EcalEBTrigPrimAnalyzer.h.

float EcalEBTrigPrimAnalyzer::tpgGeV_
private

Definition at line 64 of file EcalEBTrigPrimAnalyzer.h.

int EcalEBTrigPrimAnalyzer::tpIeta_
private

Definition at line 62 of file EcalEBTrigPrimAnalyzer.h.

int EcalEBTrigPrimAnalyzer::tpIphi_
private

Definition at line 62 of file EcalEBTrigPrimAnalyzer.h.

TTree* EcalEBTrigPrimAnalyzer::tree_
private

Definition at line 60 of file EcalEBTrigPrimAnalyzer.h.

int EcalEBTrigPrimAnalyzer::ttf_
private

Definition at line 62 of file EcalEBTrigPrimAnalyzer.h.

Generated for CMSSW Reference Manual by doxygen 1.8.11