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

Comparing UserCode/yangyong/src/RecoAnalyzer.cc (file contents):
Revision 1.1 by yangyong, Mon Feb 14 12:29:43 2011 UTC vs.
Revision 1.2 by yangyong, Mon Feb 28 16:50:13 2011 UTC

# Line 193 | Line 193 | RecoAnalyzer::RecoAnalyzer(const edm::Pa
193      
194      selePtGammaPi0EndCap_region3_ = iConfig.getParameter<double> ("selePtGammaPi0EndCap_region3");  
195      selePtPi0EndCap_region3_ = iConfig.getParameter<double> ("selePtPi0EndCap_region3");
196 <        
196 >    selePtPi0MaxEndCap_region3_ = iConfig.getParameter<double> ("selePtPi0MaxEndCap_region3");
197      
198      seleS4S9GammaEndCap_ = iConfig.getParameter<double> ("seleS4S9GammaEndCap");  
199      seleMinvMaxPi0EndCap_ = iConfig.getParameter<double> ("seleMinvMaxPi0EndCap");  
# Line 251 | Line 251 | RecoAnalyzer::RecoAnalyzer(const edm::Pa
251      
252      selePtGammaEtaEndCap_region3_ = iConfig.getParameter<double> ("selePtGammaEtaEndCap_region3");  
253      selePtEtaEndCap_region3_ = iConfig.getParameter<double> ("selePtEtaEndCap_region3");
254 <        
254 >    selePtEtaMaxEndCap_region3_ = iConfig.getParameter<double> ("selePtEtaMaxEndCap_region3");
255      
256 +
257      seleS4S9GammaEtaEndCap_ = iConfig.getParameter<double> ("seleS4S9GammaEtaEndCap");  
258      seleS9S25GammaEtaEndCap_ = iConfig.getParameter<double> ("seleS9S25GammaEtaEndCap");  
259      
# Line 880 | Line 881 | RecoAnalyzer::analyze(const edm::Event&
881      
882      }catch(std::exception& ex ){
883        nErrorPrinted++;
884 <      if(nErrorPrinted< maxErrorToPrint) cout<<"encapRecHitsEta NA.."<<endl;
884 >      if(nErrorPrinted< maxErrorToPrint) cout<<"encapRecHitsPi0 NA.."<<endl;
885      }
886    }
887  
# Line 2000 | Line 2001 | void RecoAnalyzer::doSelectionAndFillTre
2001    double ptpairCut_endcap_region1 = 0 ;
2002    double ptpairCut_endcap_region2 = 0 ;
2003    double ptpairCut_endcap_region3 = 0 ;
2004 +  double ptpairMaxCut_endcap_region3 = 0 ;
2005 +
2006  
2007    double eta_endcap_region1 = 0 ;
2008    double eta_endcap_region2 = 0 ;
# Line 2016 | Line 2019 | void RecoAnalyzer::doSelectionAndFillTre
2019        ptpairCut_endcap_region1 = selePtPi0EndCap_region1_;
2020        ptpairCut_endcap_region2 = selePtPi0EndCap_region2_;
2021        ptpairCut_endcap_region3 = selePtPi0EndCap_region3_;
2022 <            
2022 >      ptpairMaxCut_endcap_region3 = selePtPi0MaxEndCap_region3_;
2023 >
2024 >
2025        eta_endcap_region1 = region1_Pi0EndCap_;
2026        eta_endcap_region2 = region2_Pi0EndCap_;
2027        
# Line 2029 | Line 2034 | void RecoAnalyzer::doSelectionAndFillTre
2034        ptpairCut_endcap_region1 = selePtEtaEndCap_region1_;
2035        ptpairCut_endcap_region2 = selePtEtaEndCap_region2_;
2036        ptpairCut_endcap_region3 = selePtEtaEndCap_region3_;
2037 <      
2037 >      ptpairMaxCut_endcap_region3 = selePtEtaMaxEndCap_region3_;
2038 >
2039        eta_endcap_region1 = region1_EtaEndCap_;
2040        eta_endcap_region2 = region2_EtaEndCap_;
2041        
# Line 2166 | Line 2172 | void RecoAnalyzer::doSelectionAndFillTre
2172            if(ptmin < ptminCut_endcap_region2 || ptpair < ptpairCut_endcap_region2) continue;
2173          }else{
2174            if(ptmin < ptminCut_endcap_region3 || ptpair < ptpairCut_endcap_region3) continue;
2175 +          if( ptpair > ptpairMaxCut_endcap_region3 ) continue;
2176          }
2177        }else{
2178          
# Line 2173 | Line 2180 | void RecoAnalyzer::doSelectionAndFillTre
2180          
2181        }
2182        
2183 <
2183 >      
2184        mpair = sqrt ( (eClus[i] + eClus[j])*(eClus[i] + eClus[j]) - (p0x+p1x)*(p0x+p1x) - (p0y+p1y)*(p0y+p1y) - (p0z+p1z)*(p0z+p1z) );  
2185  
2186        

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines