ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitHzz4l/Util/src/ElectronEnergySmearingScaling.cc
(Generate patch)

Comparing UserCode/MitHzz4l/Util/src/ElectronEnergySmearingScaling.cc (file contents):
Revision 1.2 by anlevin, Wed Oct 17 01:29:56 2012 UTC vs.
Revision 1.6 by anlevin, Tue Nov 27 14:11:03 2012 UTC

# Line 98 | Line 98 | std::pair<float,float> apply_claudes_ele
98          if (run>=178424 && run<=180252) corr = 0.0047;
99        }
100      }
101 <    if (dataset_=="2012Jul13ReReco") {                    
101 >    else if (dataset_=="2012Jul13ReReco") {                    
102        // values from https://twiki.cern.ch/twiki/bin/view/CMS/ECALELF  
103        if (ele->IsEB() && fabs(ele->SCluster()->Eta())<1 && ele->SCluster()->R9()<0.94) {
104          if (run>=190645 && run<=190781) scale = 1.0020;
# Line 241 | Line 241 | std::pair<float,float> apply_claudes_ele
241      else
242        assert(0);
243    }
244 –  if(isMC_){
245 –    if (dataset_=="Summer12_DR53X_HCP2012") {
246 –      if (ele->IsEB() && fabs(ele->SCluster()->Eta())<1 && ele->SCluster()->R9()<0.94) dsigMC = 0.0103;
247 –      if (ele->IsEB() && fabs(ele->SCluster()->Eta())<1 && ele->SCluster()->R9()>=0.94) dsigMC = 0.0090;
248 –      if (ele->IsEB() && fabs(ele->SCluster()->Eta())>=1 && ele->SCluster()->R9()<0.94) dsigMC = 0.0190;
249 –      if (ele->IsEB() && fabs(ele->SCluster()->Eta())>=1 && ele->SCluster()->R9()>=0.94) dsigMC = 0.0156;
250 –      if (ele->IsEE() && fabs(ele->SCluster()->Eta())<2 && ele->SCluster()->R9()<0.94) dsigMC = 0.0269;
251 –      if (ele->IsEE() && fabs(ele->SCluster()->Eta())<2 && ele->SCluster()->R9()>=0.94) dsigMC = 0.0287;
252 –      if (ele->IsEE() && fabs(ele->SCluster()->Eta())>=2 && ele->SCluster()->R9()<0.94) dsigMC = 0.0364;
253 –      if (ele->IsEE() && fabs(ele->SCluster()->Eta())>=2 && ele->SCluster()->R9()>=0.94) dsigMC = 0.0321;  
254 –    }
255 –    else if (dataset_=="Fall11") { // values from https://hypernews.cern.ch/HyperNews/CMS/get/higgs2g/634.html, consistant with Jan16ReReco corrections
256 –      if (ele->IsEB() && fabs(ele->SCluster()->Eta())<1 && ele->SCluster()->R9()<0.94) dsigMC = 0.0096;
257 –      if (ele->IsEB() && fabs(ele->SCluster()->Eta())<1 && ele->SCluster()->R9()>=0.94) dsigMC = 0.0074;
258 –      if (ele->IsEB() && fabs(ele->SCluster()->Eta())>=1 && ele->SCluster()->R9()<0.94) dsigMC = 0.0196;
259 –      if (ele->IsEB() && fabs(ele->SCluster()->Eta())>=1 && ele->SCluster()->R9()>=0.94) dsigMC = 0.0141;
260 –      if (ele->IsEE() && fabs(ele->SCluster()->Eta())<2 && ele->SCluster()->R9()<0.94) dsigMC = 0.0279;
261 –      if (ele->IsEE() && fabs(ele->SCluster()->Eta())<2 && ele->SCluster()->R9()>=0.94) dsigMC = 0.0268;
262 –      if (ele->IsEE() && fabs(ele->SCluster()->Eta())>=2 && ele->SCluster()->R9()<0.94) dsigMC = 0.0301;
263 –      if (ele->IsEE() && fabs(ele->SCluster()->Eta())>=2 && ele->SCluster()->R9()>=0.94) dsigMC = 0.0293;  
264 –    }
244  
245 <    else
246 <      assert(0);
245 >  if (dataset_=="Summer12_DR53X_HCP2012" || dataset_ == "2012Jul13ReReco") {
246 >    if (ele->IsEB() && fabs(ele->SCluster()->Eta())<1 && ele->SCluster()->R9()<0.94) dsigMC = 0.0103;
247 >    if (ele->IsEB() && fabs(ele->SCluster()->Eta())<1 && ele->SCluster()->R9()>=0.94) dsigMC = 0.0090;
248 >    if (ele->IsEB() && fabs(ele->SCluster()->Eta())>=1 && ele->SCluster()->R9()<0.94) dsigMC = 0.0190;
249 >    if (ele->IsEB() && fabs(ele->SCluster()->Eta())>=1 && ele->SCluster()->R9()>=0.94) dsigMC = 0.0156;
250 >    if (ele->IsEE() && fabs(ele->SCluster()->Eta())<2 && ele->SCluster()->R9()<0.94) dsigMC = 0.0269;
251 >    if (ele->IsEE() && fabs(ele->SCluster()->Eta())<2 && ele->SCluster()->R9()>=0.94) dsigMC = 0.0287;
252 >    if (ele->IsEE() && fabs(ele->SCluster()->Eta())>=2 && ele->SCluster()->R9()<0.94) dsigMC = 0.0364;
253 >    if (ele->IsEE() && fabs(ele->SCluster()->Eta())>=2 && ele->SCluster()->R9()>=0.94) dsigMC = 0.0321;  
254    }
255 +  else if (dataset_=="Fall11" ||dataset_ == "Jan16ReReco") { // values from https://hypernews.cern.ch/HyperNews/CMS/get/higgs2g/634.html, consistant with Jan16ReReco corrections
256 +    if (ele->IsEB() && fabs(ele->SCluster()->Eta())<1 && ele->SCluster()->R9()<0.94) dsigMC = 0.0096;
257 +    if (ele->IsEB() && fabs(ele->SCluster()->Eta())<1 && ele->SCluster()->R9()>=0.94) dsigMC = 0.0074;
258 +    if (ele->IsEB() && fabs(ele->SCluster()->Eta())>=1 && ele->SCluster()->R9()<0.94) dsigMC = 0.0196;
259 +    if (ele->IsEB() && fabs(ele->SCluster()->Eta())>=1 && ele->SCluster()->R9()>=0.94) dsigMC = 0.0141;
260 +    if (ele->IsEE() && fabs(ele->SCluster()->Eta())<2 && ele->SCluster()->R9()<0.94) dsigMC = 0.0279;
261 +    if (ele->IsEE() && fabs(ele->SCluster()->Eta())<2 && ele->SCluster()->R9()>=0.94) dsigMC = 0.0268;
262 +    if (ele->IsEE() && fabs(ele->SCluster()->Eta())>=2 && ele->SCluster()->R9()<0.94) dsigMC = 0.0301;
263 +    if (ele->IsEE() && fabs(ele->SCluster()->Eta())>=2 && ele->SCluster()->R9()>=0.94) dsigMC = 0.0293;  
264 +  }
265 +  else
266 +    assert(!isMC_);
267  
268    float corrected_energy = -999;
269  
270    if (!isMC_){
271      if(dataset_ == "Jan16ReReco") corrected_energy = energy/(1+corr);
272 <    else if (!isMC_ && dataset_ == "2012Jul13ReReco") corrected_energy = energy*scale;
272 >    else if (dataset_ == "2012Jul13ReReco") corrected_energy = energy*scale;
273      else assert(0);
274    }
275    if (isMC_) {
276      TRandom1 * rand = new TRandom1();
277      rand->SetSeed(0);
278      corrMC = rand->Gaus(1,dsigMC);
279 <    corrMC = 1+dsigMC;
279 > //     corrMC = 1+dsigMC; //use this only for synchronization
280      if (debug_) std::cout << "[ElectronEnergyCalibrator] unsmeared energy " << energy << std::endl;
281      corrected_energy = energy*corrMC;  
282      if (debug_) std::cout << "[ElectronEnergyCalibrator] smeared energy " << corrected_energy << std::endl;

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines