Status report. Wataru OKAMURA 06/02/2012 ATLAS Osaka meeting

Transcription

1 Status report Wataru OKAMURA // ATLAS Osaka meeting

2 Noise occupancy test problem Boosted ZH to llbb analysis I m cleaning up my source codes and implementing any collections / smearing.

3 Noise occupancy test problem Threshold against channel X axis : channel, Y axis : Threshold [fc] 9..57_link..9 Triggers Channel _link Channel.

4 as a function of threshold [q] X axis : Threshold [fc], Y axis : Chip.... Chip Chip Chip Chip Chip Chip....7 Chip Chip Chip Chip.... Chip

5 ln() as a function of threshold The plots are fitted with a linear function Fit function : y = p + p * x ENC [e] is defined as 5*sqrt{/(*p)} Ln.e Threshold^ (fc^) Ln.e Threshold^ (fc^) Ln 9.e..... Threshold^ (fc^) Ln.e Threshold^ (fc^).7e7.e 7.e.e Noise occupancy test problem Threshold^ (fc^) Threshold^ (fc^) Threshold^ (fc^) Threshold^ (fc^) Ln Ln Ln Ln Ln.e Threshold^ (fc^) Ln.e..... Threshold^ (fc^) Ln.e..... Threshold^ (fc^) Ln.e..... Threshold^ (fc^) 5

6 ENC in Online DB vs ENC obtained by using the latest version from SVN for all modules h_estenc_bec h_estenc_bec Entries 7 9 Mean 9 ±..77 RMS Underflow 5 Overflow 79 : Online DB : SVN Integral.5e+.e+ These are not agreement.

7 auto_ptr<tf> f = NoiseTestResult::createFitFunction(); //Check for a dead chip (which hasnt been caused by vetoing due to other serious defect!) if (nbins == && defects.defectseverityencompassingelement(moduleelement::chip(ichip)) < SERIOUS ) { defects.adddefect(defect(standarddefects::dead, ModuleElement::Chip(iChip))); boost::shared_ptr<tgraphasymmerrors> g( new TGraphAsymmErrors()); std::cout<<"returning empty chip result for chip "<< ichip <<std::endl; total_dead_chips++; return ChipNOResult(g, boost::shared_ptr<tf>(f),.,.,.,.,.); } boost::shared_ptr<tgraphasymmerrors> g( new TGraphAsymmErrors(nBins, x, y, xlowererr, xuppererr, ylowererr, yuppererr)); // begin ingo float range_min =.; float range_max =.5; if(m_scanvariable==){ f>setrange(range_min,range_max); // for barrel if ( m_mur_type==sctconfiguration::endcap && m_endcap_type==sctconfiguration::endcap_outer) f>setrange(.,range_max); else if ( m_mur_type==sctconfiguration::endcap && m_endcap_type==sctconfiguration::endcap_short_middle) f>setrange(.5,range_max); else if ( m_mur_type==sctconfiguration::endcap && m_endcap_type==sctconfiguration::endcap_middle) f>setrange(.7,range_max); else if ( m_mur_type==sctconfiguration::endcap && m_endcap_type==sctconfiguration::endcap_inner) f>setrange(.7,range_max); } // end ingo These lines don t work with my script to reanalyze raw data? I make sure that in next pages. //We don't usually use FitStrategy because there is no asymm errors method. However since ROOT v./ seems to core dump with this method, we use it AJB AnalysisService::instance().getFitStrategy().fitTGraphAsymmErrors(*g, *f); double noise =.; if (f>getparameter() <.) noise = 5.*sqrt(/(*f>GetParameter())); double chisq =.; if (f>getndf() >.) { chisq = f>getchisquare() / (float)(f>getndf()); } 7

8 ENC which is obtained by the following different ways are compared. Case : Online DB Case : Fit range [. to.5] for Barrel modules Case : Fit range [. to.5] for EndCap outer modules Case : Fit range [.5 to.5] for EndCap shortmiddle modules Case : Fit range [.7 to.5] for EndCap middle and inner modules

9 ENC for Barrel modules h_estenc_b h_estenc_b Entries 5 55 Mean ± RMS Underflow Overflow 97 Integral.e+.7e+.9e+.5e+ : case : case : case : case : case Case and case are agreement. 9

10 ENC for EndCap outer modules h_estenc_ecouter h_estenc_ecouter : case Entries Mean ± RMS Underflow Overflow Integral 9 : case : case : case : case Case and case are agreement.

11 ENC for EndCap shortmiddle modules h_estenc_ecshortmiddle h_estenc_ecshortmiddle Entries 7 Mean ± RMS Underflow Overflow 5 Integral : case : case : case : case : case Case and case are agreement.

12 ENC for EndCap middle and inner modules h_estenc_ecinner h_estenc_ecinner Entries 55 Mean 5 ±.55 RMS 5. Underflow Overflow 57 Integral 57 h_estenc_ecmiddle : case : case : case : case : case h_estenc_ecmiddle Entries 57 Mean 5 5 ± RMS Underflow Overflow Integral Case and case are agreement.

13 Summary The lines which shown in page don t work with my script to reanalyze raw data. As possible causes... The reanalysis script which given from Per has some issue?? The latest version Software from SVN and one which is used now are different?? To do Compare mean occupancy and occupancy obtained by fitting.

14 backup

15 for all modules : Online DB : SVN h_meanocc_bec h_meanocc_bec Entries 7 9 Mean RMS.9e5.9e5 ±.95e7.e7.e5.e5 Underflow Overflow 7 Integral.9e+.e+!

16 for Barrel modules : Online DB : SVN h_meanocc_b h_meanocc_b Entries 5 55 Mean RMS.e5.7e5 ± 5.9e 5.9e 9.e 9.e Underflow Overflow Integral.e+.5e+!.....5

17 for EndCap modules h_meanocc_ecouter h_meanocc_ecouter Entries h_meanocc_ecmiddle : Online DB : SVN h_meanocc_ecmiddle Entries 57 Mean RMS.e5 ±.5e7 5.9e5 Mean RMS.e5 ±.e7.e5 Underflow Underflow Overflow Overflow Integral Integral 57!.....5! h_meanocc_ecshortmiddle h_meanocc_ecshortmiddle h_meanocc_ecinner h_meanocc_ecinner Entries 7 Entries 55 Mean RMS.777e.5e7 ±.e7.9e7.e5 5.7e Mean RMS.e.7e7 ± 5.e7.e7.9e5.7e Underflow Underflow Overflow Overflow 5 Integral 7 7 Integral 5 55!.....5!

18 h_etaphi_b h_etaphi_b h_etaphi_b h_etaphi_b

19 h_etaphi_eca h_etaphi_eca h_etaphi_eca h_etaphi_eca h_etaphi_eca h_etaphi_eca h_etaphi_eca h_etaphi_eca7 h_etaphi_eca

20 h_etaphi_ecc h_etaphi_ecc h_etaphi_ecc h_etaphi_ecc h_etaphi_ecc h_etaphi_ecc h_etaphi_ecc h_etaphi_ecc7 h_etaphi_ecc