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
 
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)
 
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 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 38 of file EcalTrigPrimAnalyzer.h.

Constructor & Destructor Documentation

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

Definition at line 51 of file EcalTrigPrimAnalyzer.cc.

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

53 {
54  ecal_parts_.push_back("Barrel");
55  ecal_parts_.push_back("Endcap");
56 
57  histfile_=new TFile("histos.root","RECREATE");
58  tree_ = new TTree("TPGtree","TPGtree");
59  tree_->Branch("iphi",&iphi_,"iphi/I");
60  tree_->Branch("ieta",&ieta_,"ieta/I");
61  tree_->Branch("eRec",&eRec_,"eRec/F");
62  tree_->Branch("tpgADC",&tpgADC_,"tpgADC/I");
63  tree_->Branch("tpgGeV",&tpgGeV_,"tpgGeV/F");
64  tree_->Branch("ttf",&ttf_,"ttf/I");
65  tree_->Branch("fg",&fg_,"fg/I");
66  for (unsigned int i=0;i<2;++i) {
67  ecal_et_[i]=new TH1I(ecal_parts_[i].c_str(),"Et",255,0,255);
68  char title[30];
69  sprintf(title,"%s_ttf",ecal_parts_[i].c_str());
70  ecal_tt_[i]=new TH1I(title,"TTF",10,0,10);
71  sprintf(title,"%s_fgvb",ecal_parts_[i].c_str());
72  ecal_fgvb_[i]=new TH1I(title,"FGVB",10,0,10);
73  }
74 
75  recHits_= iConfig.getParameter<bool>("AnalyzeRecHits");
76  label_=iConfig.getParameter<edm::InputTag>("inputTP");
77  if (recHits_) {
78  hTPvsRechit_= 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:277
EcalTrigPrimAnalyzer::ecal_parts_
std::vector< std::string > ecal_parts_
Definition: EcalTrigPrimAnalyzer.h:49
EcalTrigPrimAnalyzer::rechits_labelEE_
edm::InputTag rechits_labelEE_
Definition: EcalTrigPrimAnalyzer.h:65
EcalTrigPrimAnalyzer::rechits_labelEB_
edm::InputTag rechits_labelEB_
Definition: EcalTrigPrimAnalyzer.h:64
fftjetcommon_cfi.title
title
Definition: fftjetcommon_cfi.py:32
EcalTrigPrimAnalyzer::~EcalTrigPrimAnalyzer ( )
override

Definition at line 86 of file EcalTrigPrimAnalyzer.cc.

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

Member Function Documentation

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

Definition at line 104 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().

105 {
106  using namespace edm;
107  using namespace std;
108 
109  // Get input
110  edm::Handle<EcalTrigPrimDigiCollection> tp;
111  iEvent.getByLabel(label_,tp);
112  for (unsigned int i=0;i<tp.product()->size();i++) {
113  EcalTriggerPrimitiveDigi d=(*(tp.product()))[i];
114  int subdet=d.id().subDet()-1;
115  if (subdet==0) {
116  ecal_et_[subdet]->Fill(d.compressedEt());
117  }
118  else {
119  if (d.id().ietaAbs()==27 || d.id().ietaAbs()==28) {
120  if (i%2) ecal_et_[subdet]->Fill(d.compressedEt()*2.);
121  }
122  else ecal_et_[subdet]->Fill(d.compressedEt());
123  }
124  ecal_tt_[subdet]->Fill(d.ttFlag());
125  ecal_fgvb_[subdet]->Fill(d.fineGrain());
126 
127  }
128  if (!recHits_) return;
129 
130  // comparison with RecHits
131  edm::Handle<EcalRecHitCollection> rechit_EB_col;
132  iEvent.getByLabel(rechits_labelEB_,rechit_EB_col);
133 
134  edm::Handle<EcalRecHitCollection> rechit_EE_col;
135  iEvent.getByLabel(rechits_labelEE_,rechit_EE_col);
136 
137 
138  edm::ESHandle<CaloGeometry> theGeometry;
139  edm::ESHandle<CaloSubdetectorGeometry> theBarrelGeometry_handle;
140  edm::ESHandle<CaloSubdetectorGeometry> theEndcapGeometry_handle;
141  iSetup.get<CaloGeometryRecord>().get( theGeometry );
142  iSetup.get<EcalEndcapGeometryRecord>().get("EcalEndcap",theEndcapGeometry_handle);
143  iSetup.get<EcalBarrelGeometryRecord>().get("EcalBarrel",theBarrelGeometry_handle);
144 
145  const CaloSubdetectorGeometry* theEndcapGeometry = theEndcapGeometry_handle.product();
146  const CaloSubdetectorGeometry* theBarrelGeometry = theBarrelGeometry_handle.product();
147  edm::ESHandle<EcalTrigTowerConstituentsMap> eTTmap_;
148  iSetup.get<IdealGeometryRecord>().get(eTTmap_);
149 
150  map<EcalTrigTowerDetId, float> mapTow_Et;
151 
152 
153  for (unsigned int i=0;i<rechit_EB_col.product()->size();i++) {
154  const EBDetId & myid1=(*rechit_EB_col.product())[i].id();
155  EcalTrigTowerDetId towid1= myid1.tower();
156  float theta = theBarrelGeometry->getGeometry(myid1)->getPosition().theta();
157  float Etsum=((*rechit_EB_col.product())[i].energy())*sin(theta);
158  bool test_alreadyin= false;
159  map<EcalTrigTowerDetId, float>::iterator ittest= mapTow_Et.find(towid1);
160  if (ittest!= mapTow_Et.end()) test_alreadyin=true;
161  if (test_alreadyin) continue;
162  unsigned int j=i+1;
163  bool loopend=false;
164  unsigned int count=0;
165  while( j<rechit_EB_col.product()->size() && !loopend){
166  count++;
167  const EBDetId & myid2=(*rechit_EB_col.product())[j].id();
168  EcalTrigTowerDetId towid2= myid2.tower();
169  if( towid1==towid2 ) {
170  float theta=theBarrelGeometry->getGeometry(myid2)->getPosition().theta();
171  Etsum += (*rechit_EB_col.product())[j].energy()*sin(theta);
172  }
173  j++;
174  if (count>1800) loopend=true;
175  }
176  mapTow_Et.insert(pair<EcalTrigTowerDetId,float>(towid1, Etsum));
177  }
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()) test_alreadyin=true;
188  if (test_alreadyin) continue;
189  unsigned int j=i+1;
190  bool loopend=false;
191  unsigned int count=0;
192  while( j<rechit_EE_col.product()->size() && !loopend){
193  const EEDetId & myid2=(*rechit_EE_col.product())[j].id();
194  EcalTrigTowerDetId towid2= (*eTTmap_).towerOf(myid2);
195  if( towid1==towid2 ) {
196  float theta=theEndcapGeometry->getGeometry(myid2)->getPosition().theta();
197  Etsum += (*rechit_EE_col.product())[j].energy()*sin(theta);
198  }
199  // else loopend=true;
200  j++;
201  if (count>500) loopend=true;
202  }
203  // alreadyin_EE.push_back(towid1);
204  mapTow_Et.insert(pair<EcalTrigTowerDetId,float>(towid1, Etsum));
205  }
206 
207 
208  EcalTPGScale ecalScale ;
209  ecalScale.setEventSetup(iSetup) ;
210  for (unsigned int i=0;i<tp.product()->size();i++) {
211  EcalTriggerPrimitiveDigi d=(*(tp.product()))[i];
212  const EcalTrigTowerDetId TPtowid= d.id();
213  map<EcalTrigTowerDetId, float>::iterator it= mapTow_Et.find(TPtowid);
214  float Et = ecalScale.getTPGInGeV(d.compressedEt(), TPtowid) ;
215  if (d.id().ietaAbs()==27 || d.id().ietaAbs()==28) Et*=2;
216  iphi_ = TPtowid.iphi() ;
217  ieta_ = TPtowid.ieta() ;
218  tpgADC_ = d.compressedEt() ;
219  tpgGeV_ = Et ;
220  ttf_ = d.ttFlag() ;
221  fg_ = d.fineGrain() ;
222  if (it!= mapTow_Et.end()) {
223  hTPvsRechit_->Fill(it->second,Et);
224  hTPoverRechit_->Fill(Et/it->second);
225  eRec_ = it->second ;
226  }
227  tree_->Fill() ;
228  }
229 
230 }
mps_fire.i
i
Definition: mps_fire.py:277
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:48
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:39
EBDetId::tower
EcalTrigTowerDetId tower() const
get the HCAL/trigger iphi of this crystal
Definition: EBDetId.h:59
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:535
EcalTrigTowerDetId::iphi
int iphi() const
get the tower iphi
Definition: EcalTrigTowerDetId.h:57
EcalTriggerPrimitiveDigi::id
const EcalTrigTowerDetId & id() const
Definition: EcalTriggerPrimitiveDigi.h:35
edm::Handle::product
T const * product() const
Definition: Handle.h:81
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:65
EcalTrigPrimAnalyzer::rechits_labelEB_
edm::InputTag rechits_labelEB_
Definition: EcalTrigPrimAnalyzer.h:64
EcalTrigTowerDetId::subDet
EcalSubdetector subDet() const
get the subDetector associated to the Trigger Tower
Definition: EcalTrigTowerDetId.h:36
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:62
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 233 of file EcalTrigPrimAnalyzer.cc.

References mps_fire::i.

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

Member Data Documentation

TH1I* EcalTrigPrimAnalyzer::ecal_et_[2]
private

Definition at line 50 of file EcalTrigPrimAnalyzer.h.

TH1I* EcalTrigPrimAnalyzer::ecal_fgvb_[2]
private

Definition at line 52 of file EcalTrigPrimAnalyzer.h.

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

Definition at line 49 of file EcalTrigPrimAnalyzer.h.

TH1I* EcalTrigPrimAnalyzer::ecal_tt_[2]
private

Definition at line 51 of file EcalTrigPrimAnalyzer.h.

float EcalTrigPrimAnalyzer::eRec_
private

Definition at line 60 of file EcalTrigPrimAnalyzer.h.

int EcalTrigPrimAnalyzer::fg_
private

Definition at line 59 of file EcalTrigPrimAnalyzer.h.

TH1I * EcalTrigPrimAnalyzer::histBar
private

Definition at line 53 of file EcalTrigPrimAnalyzer.h.

TH1I* EcalTrigPrimAnalyzer::histEndc
private

Definition at line 53 of file EcalTrigPrimAnalyzer.h.

TFile* EcalTrigPrimAnalyzer::histfile_
private

Definition at line 54 of file EcalTrigPrimAnalyzer.h.

TH1F* EcalTrigPrimAnalyzer::hTPoverRechit_
private

Definition at line 56 of file EcalTrigPrimAnalyzer.h.

TH2F* EcalTrigPrimAnalyzer::hTPvsRechit_
private

Definition at line 55 of file EcalTrigPrimAnalyzer.h.

int EcalTrigPrimAnalyzer::ieta_
private

Definition at line 59 of file EcalTrigPrimAnalyzer.h.

int EcalTrigPrimAnalyzer::iphi_
private

Definition at line 59 of file EcalTrigPrimAnalyzer.h.

edm::InputTag EcalTrigPrimAnalyzer::label_
private

Definition at line 62 of file EcalTrigPrimAnalyzer.h.

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

bool EcalTrigPrimAnalyzer::recHits_
private

Definition at line 67 of file EcalTrigPrimAnalyzer.h.

edm::InputTag EcalTrigPrimAnalyzer::rechits_labelEB_
private

Definition at line 64 of file EcalTrigPrimAnalyzer.h.

edm::InputTag EcalTrigPrimAnalyzer::rechits_labelEE_
private

Definition at line 65 of file EcalTrigPrimAnalyzer.h.

int EcalTrigPrimAnalyzer::tpgADC_
private

Definition at line 59 of file EcalTrigPrimAnalyzer.h.

float EcalTrigPrimAnalyzer::tpgGeV_
private

Definition at line 60 of file EcalTrigPrimAnalyzer.h.

TTree* EcalTrigPrimAnalyzer::tree_
private

Definition at line 57 of file EcalTrigPrimAnalyzer.h.

int EcalTrigPrimAnalyzer::ttf_
private

Definition at line 59 of file EcalTrigPrimAnalyzer.h.

Generated for CMSSW Reference Manual by doxygen 1.8.11