Go to the documentation of this file.00001 #include "DataFormats/HcalDetId/interface/HcalDetId.h"
00002 #include "CalibCalorimetry/CaloMiscalibTools/interface/HcalRecHitRecalib.h"
00003
00004 #include "DataFormats/Common/interface/Handle.h"
00005 #include "FWCore/ParameterSet/interface/FileInPath.h"
00006 #include "FWCore/Framework/interface/ESHandle.h"
00007
00008 #include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"
00009 #include "FWCore/MessageLogger/interface/MessageLogger.h"
00010 #include "CalibCalorimetry/CaloMiscalibTools/interface/MiscalibReaderFromXMLHcal.h"
00011
00012 HcalRecHitRecalib::HcalRecHitRecalib(const edm::ParameterSet& iConfig)
00013 {
00014
00015 hbheLabel_ = iConfig.getParameter<edm::InputTag>("hbheInput");
00016 hoLabel_ = iConfig.getParameter<edm::InputTag>("hoInput");
00017 hfLabel_ = iConfig.getParameter<edm::InputTag>("hfInput");
00018
00019
00020
00021
00022
00023
00024
00025
00026 RecalibHBHEHits_ = iConfig.getParameter< std::string > ("RecalibHBHEHitCollection");
00027 RecalibHFHits_ = iConfig.getParameter< std::string > ("RecalibHFHitCollection");
00028 RecalibHOHits_ = iConfig.getParameter< std::string > ("RecalibHOHitCollection");
00029
00030 refactor_ = iConfig.getUntrackedParameter<double> ("Refactor",(double)1);
00031 refactor_mean_ = iConfig.getUntrackedParameter<double> ("Refactor_mean",(double)1);
00032
00033
00034 produces< HBHERecHitCollection >(RecalibHBHEHits_);
00035 produces< HFRecHitCollection >(RecalibHFHits_);
00036 produces< HORecHitCollection >(RecalibHOHits_);
00037
00038
00039
00040 hcalfileinpath_=iConfig.getUntrackedParameter<std::string> ("fileNameHcal","");
00041 edm::FileInPath hcalfiletmp("CalibCalorimetry/CaloMiscalibTools/data/"+hcalfileinpath_);
00042
00043 hcalfile_=hcalfiletmp.fullPath();
00044 }
00045
00046
00047 HcalRecHitRecalib::~HcalRecHitRecalib()
00048 {
00049
00050
00051 }
00052
00053 void
00054 HcalRecHitRecalib::beginRun(edm::Run&, edm::EventSetup const& iSetup)
00055 {
00056 edm::ESHandle<HcalTopology> topology;
00057 iSetup.get<IdealGeometryRecord>().get( topology );
00058
00059 mapHcal_.prefillMap(*topology);
00060
00061 MiscalibReaderFromXMLHcal hcalreader_(mapHcal_);
00062 if(!hcalfile_.empty()) hcalreader_.parseXMLMiscalibFile(hcalfile_);
00063 mapHcal_.print();
00064 }
00065
00066
00067 void
00068 HcalRecHitRecalib::produce(edm::Event& iEvent, const edm::EventSetup& iSetup)
00069 {
00070 using namespace edm;
00071 using namespace std;
00072
00073 Handle<HBHERecHitCollection> HBHERecHitsHandle;
00074 Handle<HFRecHitCollection> HFRecHitsHandle;
00075 Handle<HORecHitCollection> HORecHitsHandle;
00076
00077 const HBHERecHitCollection* HBHERecHits = 0;
00078 const HFRecHitCollection* HFRecHits = 0;
00079 const HORecHitCollection* HORecHits = 0;
00080
00081 iEvent.getByLabel(hbheLabel_,HBHERecHitsHandle);
00082 if (!HBHERecHitsHandle.isValid()) {
00083 LogDebug("") << "HcalREcHitRecalib: Error! can't get product!" << std::endl;
00084 } else {
00085 HBHERecHits = HBHERecHitsHandle.product();
00086 }
00087
00088 iEvent.getByLabel(hoLabel_,HORecHitsHandle);
00089 if (!HORecHitsHandle.isValid()) {
00090 LogDebug("") << "HcalREcHitRecalib: Error! can't get product!" << std::endl;
00091 } else {
00092 HORecHits = HORecHitsHandle.product();
00093 }
00094
00095 iEvent.getByLabel(hfLabel_,HFRecHitsHandle);
00096 if (!HFRecHitsHandle.isValid()) {
00097 LogDebug("") << "HcalREcHitRecalib: Error! can't get product!" << std::endl;
00098 } else {
00099 HFRecHits = HFRecHitsHandle.product();
00100 }
00101
00102
00103
00104
00105
00106
00107
00108
00109
00110
00111
00112
00113
00114 std::auto_ptr< HBHERecHitCollection > RecalibHBHERecHitCollection( new HBHERecHitCollection );
00115 std::auto_ptr< HFRecHitCollection > RecalibHFRecHitCollection( new HFRecHitCollection );
00116 std::auto_ptr< HORecHitCollection > RecalibHORecHitCollection( new HORecHitCollection );
00117
00118
00119
00120
00121
00122
00123 if(HBHERecHits)
00124 {
00125
00126
00127 HBHERecHitCollection::const_iterator itHBHE;
00128 for (itHBHE=HBHERecHits->begin(); itHBHE!=HBHERecHits->end(); itHBHE++) {
00129
00130
00131
00132
00133
00134
00135
00136
00137
00138
00139
00140
00141
00142
00143
00144 float icalconst=(mapHcal_.get().find(itHBHE->id().rawId()))->second;
00145
00146
00147 icalconst=refactor_mean_+(icalconst-refactor_mean_)*refactor_;
00148 HBHERecHit aHit(itHBHE->id(),itHBHE->energy()*icalconst,itHBHE->time());
00149
00150 RecalibHBHERecHitCollection->push_back( aHit);
00151 }
00152 }
00153
00154 if(HFRecHits)
00155 {
00156
00157
00158 HFRecHitCollection::const_iterator itHF;
00159 for (itHF=HFRecHits->begin(); itHF!=HFRecHits->end(); itHF++) {
00160
00161
00162
00163
00164
00165
00166
00167
00168
00169
00170
00171
00172
00173
00174
00175
00176
00177 float icalconst=(mapHcal_.get().find(itHF->id().rawId()))->second;
00178 icalconst=refactor_mean_+(icalconst-refactor_mean_)*refactor_;
00179 HFRecHit aHit(itHF->id(),itHF->energy()*icalconst,itHF->time());
00180
00181 RecalibHFRecHitCollection->push_back( aHit);
00182 }
00183 }
00184
00185 if(HORecHits)
00186 {
00187
00188
00189 HORecHitCollection::const_iterator itHO;
00190 for (itHO=HORecHits->begin(); itHO!=HORecHits->end(); itHO++) {
00191
00192
00193
00194
00195
00196
00197
00198
00199
00200
00201
00202
00203
00204
00205
00206
00207
00208 float icalconst=(mapHcal_.get().find(itHO->id().rawId()))->second;
00209 icalconst=refactor_mean_+(icalconst-refactor_mean_)*refactor_;
00210 HORecHit aHit(itHO->id(),itHO->energy()*icalconst,itHO->time());
00211
00212 RecalibHORecHitCollection->push_back( aHit);
00213 }
00214 }
00215
00216
00217
00218 iEvent.put( RecalibHBHERecHitCollection, RecalibHBHEHits_);
00219 iEvent.put( RecalibHFRecHitCollection, RecalibHFHits_);
00220 iEvent.put( RecalibHORecHitCollection, RecalibHOHits_);
00221 }
00222