#include <DQA_TOF.h>
Public Member Functions | |
| DQA_TOF (const std::string &name, ISvcLocator *pSvcLocator) | |
| DQA_TOF (const std::string &name, ISvcLocator *pSvcLocator) | |
| StatusCode | execute () |
| StatusCode | execute () |
| StatusCode | finalize () |
| StatusCode | finalize () |
| StatusCode | Fit_delT () |
| StatusCode | Fit_delT () |
| StatusCode | initialize () |
| StatusCode | initialize () |
| StatusCode | regHist () |
| StatusCode | regHist () |
Private Attributes | |
| TH2F * | B_path |
| TH2F * | B_path |
| TH1F * | Bt_delT |
| TH1F * | Bt_delT |
| TH2F * | cluster |
| TH2F * | cluster |
| TH2F * | counter |
| TH2F * | counter |
| TGraphErrors * | delT4 |
| TGraphErrors * | delT4 |
| TH2F * | delT_z1 |
| TH2F * | delT_z1 |
| TH2F * | delT_z2 |
| TH2F * | delT_z2 |
| TH2F * | delT_z3 |
| TH2F * | delT_z3 |
| TH2F * | delT_z4 |
| TH2F * | delT_z4 |
| TH2F * | E_delT |
| TH2F * | E_delT |
| TH2F * | E_delT_Q |
| TH2F * | E_delT_Q |
| TH2F * | E_path |
| TH2F * | E_path |
| TH2F * | EC_delT |
| TH2F * | EC_delT |
| TH1F * | Et_delT |
| TH1F * | Et_delT |
| TH1F * | h_Bzrhit |
| TH1F * | h_Bzrhit |
| TH1F * | h_Ezrhit |
| TH1F * | h_Ezrhit |
| TH1F * | h_path |
| TH1F * | h_path |
| TH1F * | h_ph |
| TH1F * | h_ph |
| ITHistSvc * | m_thsvc |
| ITHistSvc * | m_thsvc |
| TH1F * | S_BE |
| TH1F * | S_BE |
| TH1F * | S_cluster |
| TH1F * | S_cluster |
| TH1F * | S_counter |
| TH1F * | S_counter |
| TH1F * | S_delT_z1 |
| TH1F * | S_delT_z1 |
| TH1F * | S_delT_z2 |
| TH1F * | S_delT_z2 |
| TH1F * | S_delT_z3 |
| TH1F * | S_delT_z3 |
| TH1F * | S_delT_z4 |
| TH1F * | S_delT_z4 |
| TH1F * | S_E_delT |
| TH1F * | S_E_delT |
| TH1F * | S_E_delT_Q |
| TH1F * | S_E_delT_Q |
| TH1F * | S_EC_delT |
| TH1F * | S_EC_delT |
| TH1F * | S_W_delT |
| TH1F * | S_W_delT |
| TH1F * | S_W_delT_Q |
| TH1F * | S_W_delT_Q |
| TH1F * | SE_offset_Q |
| TH1F * | SE_offset_Q |
| TH1F * | SE_offset_z |
| TH1F * | SE_offset_z |
| TH1F * | SW_offset_Q |
| TH1F * | SW_offset_Q |
| TH1F * | SW_offset_z |
| TH1F * | SW_offset_z |
| TH2F * | W_delT |
| TH2F * | W_delT |
| TH2F * | W_delT_Q |
| TH2F * | W_delT_Q |
|
||||||||||||
|
00025 :
00026 Algorithm(name, pSvcLocator) {
00027
00028 //Declare the properties
00029 }
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
|
|
||||||||||||
|
|
|
|
|
|
|
00112 {
00113 MsgStream log(msgSvc(), name());
00114 log << MSG::INFO << "in execute()" << endreq;
00115
00116 SmartDataPtr<Event::EventHeader> eventHeader(eventSvc(),"/Event/EventHeader");
00117 int runNo=eventHeader->runNumber();
00118 int event=eventHeader->eventNumber();
00119 log<<MSG::DEBUG<<"run,evtnum="<<runNo<<","<<event<<endreq;
00120
00121 SmartDataPtr<DQAEvent::DQAEvent> dqaevt(eventSvc(), "/Event/DQATag");
00122 if ( !dqaevt ) {
00123 log << MSG::ERROR << "Error accessing DQAEvent" << endreq;
00124 return StatusCode::FAILURE;
00125 }
00126
00127 log << MSG::DEBUG << "event tag = " << dqaevt->EventTag() << endreq;
00128
00129 if ( !dqaevt->Bhabha() )return StatusCode::SUCCESS;
00130 SmartDataPtr<EvtRecEvent> evtRecEvent(eventSvc(), EventModel::EvtRec::EvtRecEvent);
00131 SmartDataPtr<EvtRecTrackCol> evtRecTrkCol(eventSvc(), EventModel::EvtRec::EvtRecTrackCol);
00132 for(int i = 0; i < evtRecEvent->totalCharged(); i++){
00133 EvtRecTrackIterator itTrk=evtRecTrkCol->begin() + i;
00134 if ( (*itTrk)->partId() != 1 ) continue; // only e+, e-
00135 if(!(*itTrk)->isTofTrackValid())continue;
00136 SmartRefVector<RecTofTrack> tofTrkCol = (*itTrk)->tofTrack();
00137 SmartRefVector<RecTofTrack>::iterator iter_tof = tofTrkCol.begin();
00138 for(;iter_tof != tofTrkCol.end(); iter_tof++ ) {
00139 TofHitStatus *status = new TofHitStatus;
00140 status->setStatus((*iter_tof)->status());
00141 int tofid=(*iter_tof)->tofID();
00142 double tof=(*iter_tof)->tof();
00143 double texpE=(*iter_tof)->texpElectron();
00144 double path=(*iter_tof)->path();
00145 double Q=(*iter_tof)->ph();
00146 double zrhit=(*iter_tof)->zrhit();
00147
00148 if(status->is_barrel()){
00149 h_Bzrhit->Fill(zrhit);
00150 h_path->Fill(path);
00151 //h_ph->Fill(Q);
00152 if(status->is_readout()){
00153 h_ph->Fill(Q);
00154 B_path->Fill(zrhit,path);
00155 if(status->is_east()){ //barrel readout east
00156 E_delT->Fill(tofid,tof-texpE);
00157 delT_z1->Fill(zrhit,tof-texpE);
00158 E_delT_Q->Fill(Q,tof-texpE);
00159 }
00160 else { //barrel readout west
00161 W_delT->Fill(tofid,tof-texpE);
00162 delT_z2->Fill(zrhit,tof-texpE);
00163 W_delT_Q->Fill(Q,tof-texpE);
00164 }
00165 } //end readout
00166 if(!status->is_readout()&&status->is_counter()){ //barrel counter
00167 counter->Fill(tofid,tof-texpE);
00168 delT_z3->Fill(zrhit,tof-texpE);
00169 }
00170 if(!status->is_counter()&&status->is_cluster()){ //barrel cluster
00171 cluster->Fill(tofid,tof-texpE);
00172 Bt_delT->Fill(tof-texpE);
00173 delT_z4->Fill(zrhit,tof-texpE);
00174 }
00175 }
00176 else{ //endcap
00177 E_path->Fill(zrhit,path);
00178 h_Ezrhit->Fill(zrhit);
00179 EC_delT->Fill(tofid,tof-texpE);
00180 Et_delT->Fill(tof-texpE);
00181 }
00182 } // loop 7 info fo tof
00183 } // loop every track
00184 return StatusCode::SUCCESS;
00185
00186 }
|
|
|
|
|
|
00188 {
00189 MsgStream log(msgSvc(), name());
00190 log << MSG::INFO << "in finalize()" << endmsg;
00191 //delT------->sigma
00192 DQA_TOF::Fit_delT();
00193
00194 //regHist
00195 DQA_TOF::regHist();
00196 return StatusCode::SUCCESS;
00197 }
|
|
|
|
|
|
00198 {
00199 TF1 *gauss1 = new TF1("gauss","gaus");
00200 for( int tofid=0; tofid<176; tofid++ ) {
00201 int binx = W_delT->GetXaxis()->FindBin( tofid );
00202 TH1D* hTmp = W_delT->ProjectionY("",binx,binx);
00203 hTmp->Fit(gauss1,"Q","",-1.5,1.5);
00204 S_W_delT->Fill(tofid+1,gauss1->GetParameter(2));
00205 S_W_delT->SetBinError(tofid+1,gauss1->GetParError(2));
00206 delete hTmp;
00207 }
00208 for( int tofid=0; tofid<176; tofid++ ) {
00209 int binx = E_delT->GetXaxis()->FindBin( tofid );
00210 TH1D* hTmp = E_delT->ProjectionY("",binx,binx);
00211 hTmp->Fit(gauss1,"Q","",-1.5,1.5);
00212 S_E_delT->Fill(tofid+1,gauss1->GetParameter(2));
00213 S_E_delT->SetBinError(tofid+1,gauss1->GetParError(2));
00214 delete hTmp;
00215 }
00216 for( int tofid=0; tofid<176; tofid++ ) {
00217 int binx = counter->GetXaxis()->FindBin( tofid );
00218 TH1D* hTmp = counter->ProjectionY("",binx,binx);
00219 hTmp->Fit(gauss1,"Q","",-1.5,1.5);
00220 S_counter->Fill(tofid+1,gauss1->GetParameter(2));
00221 S_counter->SetBinError(tofid+1,gauss1->GetParError(2));
00222 delete hTmp;
00223 }
00224 for( int tofid=0; tofid<88; tofid++ ) {
00225 int binx = cluster->GetXaxis()->FindBin( tofid );
00226 TH1D* hTmp = cluster->ProjectionY("",binx,binx);
00227 hTmp->Fit(gauss1,"Q","",-1.5,1.5);
00228 S_cluster->Fill(tofid+1,gauss1->GetParameter(2));
00229 S_cluster->SetBinError(tofid+1,gauss1->GetParError(2));
00230 delete hTmp;
00231 }
00232 for( int tofid=0; tofid<96; tofid++ ) {
00233 int binx = EC_delT->GetXaxis()->FindBin( tofid );
00234 TH1D* hTmp = EC_delT->ProjectionY("",binx,binx);
00235 hTmp->Fit(gauss1,"Q","",-1.5,1.5);
00236 S_EC_delT->Fill(tofid+1,gauss1->GetParameter(2));
00237 S_EC_delT->SetBinError(tofid+1,gauss1->GetParError(2));
00238 delete hTmp;
00239 }
00240
00241 Bt_delT->Fit(gauss1,"Q","",-1.5,1.5);
00242 S_BE->Fill(1,gauss1->GetParameter(2));
00243 S_BE->SetBinError(1,gauss1->GetParError(2));
00244 Et_delT->Fit(gauss1,"Q","",-1.5,1.5);
00245 S_BE->Fill(2,gauss1->GetParameter(2));
00246 S_BE->SetBinError(2,gauss1->GetParError(2));
00247
00248 // 4 delT : bin is base on event numbers;
00249 if(h_Bzrhit->GetEntries()!=0){
00250 int xfirst=h_Bzrhit->GetXaxis()->GetFirst();
00251 int xlast=h_Bzrhit->GetXaxis()->GetLast();
00252 int Nentries=h_Bzrhit->GetEntries();
00253 int ibin=Nentries/9;
00254 std::vector<int> Vbin;
00255 Vbin.push_back(xfirst);
00256 //divide delT_z4--> bins
00257 int num=1;
00258 int entries=0;
00259 for(int i=xfirst;i<xlast;i++){
00260 //for(int j=delT_z3->GetYaxis()->GetFirst();j<delT_z3->GetYaxis()->GetLast();j++){
00261 entries += h_Bzrhit->GetBinContent(i);
00262 //}
00263 if(entries>=num*ibin){
00264 std::cout<<"zhaohs is:"<<entries<<","<<i<<std::endl;
00265 num++;
00266 Vbin.push_back(i);
00267 }
00268 }
00269 Vbin.push_back(xlast);
00270 int Gr_index=0;
00271 for(int i=0;i<Vbin.size()-1;i++){
00272 //int nbin=240/6;
00273 //int binx1=delT_z1->GetXaxis()->FindBin(i*nbin);
00274 TH1D* htmp1=delT_z1->ProjectionY("",Vbin[i],Vbin[i+1]);
00275 htmp1->Fit(gauss1,"Q","",-1.5,1.5);
00276 int QQ=(Vbin[i]+Vbin[i+1])/2-120; //define bin: 0-240 ----z : -120 --- 120
00277 std::cout<<"S_delT_z:"<<QQ<<std::endl;
00278 S_delT_z1->Fill(QQ,gauss1->GetParameter(2));
00279 S_delT_z1->SetBinError(QQ,gauss1->GetParError(2));
00280 Gr_index++;
00281 delT4->SetPoint(Gr_index,(double)QQ,gauss1->GetParameter(2));
00282 delT4->SetPointError(Gr_index,0.0,gauss1->GetParError(2));
00283 SE_offset_z->Fill(QQ,gauss1->GetParameter(1));
00284 SE_offset_z->SetBinError(QQ,gauss1->GetParError(1));
00285
00286 //int binx2=delT_z2->GetXaxis()->FindBin(i*nbin);
00287 TH1D* htmp2=delT_z2->ProjectionY("",Vbin[i],Vbin[i+1]);
00288 htmp2->Fit(gauss1,"Q","",-1.5,1.5);
00289 S_delT_z2->Fill(QQ,gauss1->GetParameter(2));
00290 S_delT_z2->SetBinError(QQ,gauss1->GetParError(2));
00291 Gr_index++;
00292 delT4->SetPoint(Gr_index ,(double)QQ,gauss1->GetParameter(2));
00293 delT4->SetPointError(Gr_index,0.0,gauss1->GetParError(2));
00294 SW_offset_z->Fill(QQ,gauss1->GetParameter(1));
00295 SW_offset_z->SetBinError(QQ,gauss1->GetParError(1));
00296
00297 //int binx3=delT_z3->GetXaxis()->FindBin(i*nbin);
00298 TH1D* htmp3=delT_z3->ProjectionY("",Vbin[i],Vbin[i+1]);
00299 htmp3->Fit(gauss1,"Q","",-1.5,1.5);
00300 S_delT_z3->Fill(QQ,gauss1->GetParameter(2));
00301 S_delT_z3->SetBinError(QQ,gauss1->GetParError(2));
00302 Gr_index++;
00303 delT4->SetPoint(Gr_index ,(double)QQ,gauss1->GetParameter(2));
00304 delT4->SetPointError(Gr_index,0.0,gauss1->GetParError(2));
00305
00306 //int binx4=delT_z4->GetXaxis()->FindBin(i*nbin);
00307 TH1D* htmp4=delT_z4->ProjectionY("",Vbin[i],Vbin[i+1]);
00308 htmp4->Fit(gauss1,"Q","",-1.5,1.5);
00309 S_delT_z4->Fill(QQ,gauss1->GetParameter(2));
00310 S_delT_z4->SetBinError(QQ,gauss1->GetParError(2));
00311 Gr_index++;
00312 delT4->SetPoint(Gr_index ,(double)QQ,gauss1->GetParameter(2));
00313 delT4->SetPointError(Gr_index,0.0,gauss1->GetParError(2));
00314 //delete htmp1;delete htmp2; delete htmp3; delete htmp4;
00315 }
00316 }
00317 //divide Q base on entries;
00318 int Nbin=h_ph->GetEntries();
00319 if(Nbin!=0){
00320 int ibin=Nbin/8;//fill points
00321 std::cout<<"ibin:"<<ibin<<std::endl;
00322 int entries=0;
00323 int num=1;
00324 std::vector<int> Vbin;
00325 Vbin.push_back(h_ph->GetXaxis()->GetFirst());
00326 for(int xfirst=h_ph->GetXaxis()->GetFirst();xfirst<h_ph->GetXaxis()->GetLast();xfirst++){
00327 entries += h_ph->GetBinContent(xfirst);
00328 if(entries>=ibin*num){
00329 Vbin.push_back(xfirst);
00330 std::cout<<"xfirst:"<<xfirst<<std::endl;
00331 num++;
00332 }
00333
00334 }
00335 Vbin.push_back(h_ph->GetXaxis()->GetLast());
00336 for(int i=0;i<Vbin.size()-1;i++){
00337 //int nbin=9000/9;
00338 //int binx=W_delT_Q->GetXaxis()->FindBin(i*nbin);
00339 TH1D* htmp1=W_delT_Q->ProjectionY("",Vbin[i],Vbin[i+1]);
00340 std::cout<<"Vin is:"<<Vbin[i+1]<<std::endl;
00341 int nbin=(Vbin[i]+Vbin[i+1])*10/2;
00342 std::cout<<"nbin:"<<nbin<<std::endl;
00343
00344 htmp1->Fit(gauss1,"Q","",-1.5,1.5);
00345 S_W_delT_Q->Fill(nbin,gauss1->GetParameter(2));
00346 S_W_delT_Q->SetBinError(nbin,gauss1->GetParError(2));
00347 SW_offset_Q->Fill(nbin,gauss1->GetParameter(1));
00348 SW_offset_Q->SetBinError(nbin,gauss1->GetParError(1));
00349
00350 //int biny=E_delT_Q->GetXaxis()->FindBin(i*nbin);
00351 TH1D* htmp2=E_delT_Q->ProjectionY("",Vbin[i],Vbin[i+1]);
00352 htmp2->Fit(gauss1,"Q","",-1.5,1.5);
00353 S_E_delT_Q->Fill(nbin,gauss1->GetParameter(2));
00354 S_E_delT_Q->SetBinError(nbin,gauss1->GetParError(2));
00355 SE_offset_Q->Fill(nbin,gauss1->GetParameter(1));
00356 SE_offset_Q->SetBinError(nbin,gauss1->GetParError(1));
00357
00358 }
00359 }
00360 return StatusCode::SUCCESS;
00361 }
|
|
|
|
|
|
00031 {
00032 MsgStream log(msgSvc(), name());
00033
00034 log << MSG::INFO << "in initialize()" << endmsg;
00035 StatusCode status;
00036
00037 StatusCode sc=service("THistSvc", m_thsvc);
00038 if(sc.isFailure()) {
00039 log << MSG::ERROR << "Couldn't get THistSvc" << endreq;
00040 return StatusCode::FAILURE;
00041 }
00042
00043 h_path=new TH1F("h_path","barrel ",200,0,200);
00044 h_Bzrhit=new TH1F("h_Bzrhit","barrel z hitmap",240,-120,120);
00045 h_Ezrhit=new TH1F("h_Ezrhit","endcap z hitmap",50,40,90);
00046 h_ph=new TH1F("h_ph","barrel Q",900,0,9000);
00047
00048 W_delT=new TH2F("W_delT","barrel west PMT delT",176,0,176,300,-1.5,1.5);
00049 E_delT=new TH2F("E_delT","barrel east PMT delT",176,0,176,300,-1.5,1.5);
00050 counter=new TH2F("counter","barrel counter delT",176,0,176,300,-1.5,1.5);
00051 cluster=new TH2F("cluster","barrel cluster delT",88,0,88,300,-1.5,1.5);
00052 EC_delT=new TH2F("EC_delT","endcap delT",96,0,96,300,-1.5,1.5);
00053 Bt_delT=new TH1F("Bt_delT","barrel delT",300,-1.5,1.5);
00054 Et_delT=new TH1F("Et_delT","endcap delT",300,-1.5,1.5);
00055
00056 B_path=new TH2F("B_path","barrel flight distance vs z",240,-120,120,200,0.0,200.0);
00057 E_path=new TH2F("E_path","endcap path distance vs z",50,40,90,200,0.0,200.0);
00058
00059 delT_z1=new TH2F("delT_z1","barrel east delT vs Z",240,-120,120,300,-1.5,1.5);
00060 delT_z2=new TH2F("delT_z2","barrel west delT vs Z",240,-120,120,300,-1.5,1.5);
00061 delT_z3=new TH2F("delT_z3","barrel counter delT vs Z",240,-120,120,300,-1.5,1.5);
00062 delT_z4=new TH2F("delT_z4","barrel cluster delT vs Z",240,-120,120,300,-1.5,1.5);
00063
00064 W_delT_Q=new TH2F("W_delT_Q","west barrel delT vs Q",900,0,9000,300,-1.5,1.5);
00065 E_delT_Q=new TH2F("E_delT_Q","east barrel delT vs Q",900,0,9000,300,-1.5,1.5);
00066
00067 //EFA-----sigma
00068 S_W_delT=new TH1F("S_W_delT","barrel west PMT sigma",176,0.5,176.5);
00069 S_W_delT->GetYaxis()->SetRangeUser(0,0.5);
00070 S_E_delT=new TH1F("S_E_delT","barrel east PMT sigma",176,0.5,176.5);
00071 S_E_delT->GetYaxis()->SetRangeUser(0,0.5);
00072 S_counter=new TH1F("S_counter","barrel counter sigma",176,0.5,176.5);
00073 S_counter->GetYaxis()->SetRangeUser(0,0.5);
00074 S_cluster=new TH1F("S_cluster","barrel cluster sigma",88,0.5,88.5);
00075 S_cluster->GetYaxis()->SetRangeUser(0,0.5);
00076 S_EC_delT=new TH1F("S_EC_delT","endcap sigma",96,0.5,96.5);
00077 S_EC_delT->GetYaxis()->SetRangeUser(0,0.7);
00078 S_BE=new TH1F("S_BE","barrel and endcap sigma",2,0.5,2.5);
00079 S_delT_z1=new TH1F("S_delT_z1","barrel east sigma vs z",240,-120,120);
00080 S_delT_z1->GetYaxis()->SetRangeUser(0,0.5);
00081 S_delT_z2=new TH1F("S_delT_z2","barrel west sigma vs z",240,-120,120);
00082 S_delT_z2->GetYaxis()->SetRangeUser(0,0.5);
00083 S_delT_z3=new TH1F("S_delT_z3","barrel counter sigma vs z",240,-120,120);
00084 S_delT_z3->GetYaxis()->SetRangeUser(0,0.5);
00085 S_delT_z4=new TH1F("S_delT_z4","barrel cluster sigma vs z",240,-120,120);
00086 S_delT_z4->GetYaxis()->SetRangeUser(0,0.5);
00087
00088 SE_offset_z=new TH1F("SE_offset_z","east offset vs z",240,-120,120);
00089 SE_offset_z->GetYaxis()->SetRangeUser(-0.15,0.15);
00090 SW_offset_z=new TH1F("SW_offset_z","west offset vs z",240,-120,120);
00091 SW_offset_z->GetYaxis()->SetRangeUser(-0.15,0.15);
00092
00093 S_W_delT_Q=new TH1F("S_W_delT_Q","wset sigma vs Q",900,0,9000);
00094 S_W_delT_Q->GetYaxis()->SetRangeUser(0.,0.5);
00095 S_E_delT_Q=new TH1F("S_E_delT_Q","east sigma vs Q",900,0,9000);
00096 S_E_delT_Q->GetYaxis()->SetRangeUser(0.,0.5);
00097 SW_offset_Q=new TH1F("SW_offset_Q","wset offset vs Q",900,0,9000);
00098 SW_offset_Q->GetYaxis()->SetRangeUser(-0.15,0.15);
00099 SE_offset_Q=new TH1F("SE_offset_Q","east offset vs Q",900,0,9000);
00100 SE_offset_Q->GetYaxis()->SetRangeUser(-0.15,0.15);
00101
00102 delT4=new TGraphErrors();
00103 delT4->SetTitle("sigma");
00104 delT4->SetName("sigma4");
00105 delT4->SetMarkerColor(4);
00106 delT4->SetMarkerStyle(21);
00107
00108 log << MSG::INFO << "DQA_TOF successfully return from initialize()" <<endmsg;
00109 return StatusCode::SUCCESS;
00110 }
|
|
|
|
|
|
00363 {
00364 m_thsvc->regHist("/DQAHist/TOF/W_delT", W_delT);
00365 m_thsvc->regHist("/DQAHist/TOF/E_delT", E_delT);
00366 m_thsvc->regHist("/DQAHist/TOF/counter", counter);
00367 m_thsvc->regHist("/DQAHist/TOF/cluster", cluster);
00368 m_thsvc->regHist("/DQAHist/TOF/EC_delT", EC_delT);
00369 m_thsvc->regHist("/DQAHist/TOF/Bt_delT", Bt_delT);
00370 m_thsvc->regHist("/DQAHist/TOF/Et_delT", Et_delT);
00371
00372 m_thsvc->regHist("/DQAHist/TOF/B_path", B_path);
00373 m_thsvc->regHist("/DQAHist/TOF/E_path", E_path);
00374
00375 m_thsvc->regHist("/DQAHist/TOF/delT_z1",delT_z1);
00376 m_thsvc->regHist("/DQAHist/TOF/delT_z2",delT_z2);
00377 m_thsvc->regHist("/DQAHist/TOF/delT_z3",delT_z3);
00378 m_thsvc->regHist("/DQAHist/TOF/delT_z4",delT_z4);
00379
00380 m_thsvc->regHist("/DQAHist/TOF/W_delT_Q",W_delT_Q);
00381 m_thsvc->regHist("/DQAHist/TOF/E_delT_Q",E_delT_Q);
00382
00383 m_thsvc->regHist("/DQAHist/TOF/S_W_delT",S_W_delT);
00384 m_thsvc->regHist("/DQAHist/TOF/S_E_delT",S_E_delT);
00385 m_thsvc->regHist("/DQAHist/TOF/S_counter",S_counter);
00386 m_thsvc->regHist("/DQAHist/TOF/S_cluster",S_cluster);
00387 m_thsvc->regHist("/DQAHist/TOF/S_EC_delT",S_EC_delT);
00388 m_thsvc->regHist("/DQAHist/TOF/S_BE",S_BE);
00389 m_thsvc->regHist("/DQAHist/TOF/SE_offset_z",SE_offset_z);
00390 m_thsvc->regHist("/DQAHist/TOF/SW_offset_z",SW_offset_z);
00391 m_thsvc->regHist("/DQAHist/TOF/S_delT_z1",S_delT_z1);
00392 m_thsvc->regHist("/DQAHist/TOF/S_delT_z2",S_delT_z2);
00393 m_thsvc->regHist("/DQAHist/TOF/S_delT_z3",S_delT_z3);
00394 m_thsvc->regHist("/DQAHist/TOF/S_delT_z4",S_delT_z4);
00395 m_thsvc->regHist("/DQAHist/TOF/S_W_delT_Q",S_W_delT_Q);
00396 m_thsvc->regHist("/DQAHist/TOF/S_E_delT_Q",S_E_delT_Q);
00397 m_thsvc->regHist("/DQAHist/TOF/SW_offset_Q",SW_offset_Q);
00398 m_thsvc->regHist("/DQAHist/TOF/SE_offset_Q",SE_offset_Q);
00399
00400 m_thsvc->regHist("/DQAHist/TOF/h_Bzrhit",h_Bzrhit);
00401 m_thsvc->regHist("/DQAHist/TOF/h_Ezrhit",h_Ezrhit);
00402 m_thsvc->regHist("/DQAHist/TOF/h_path",h_path);
00403 m_thsvc->regHist("/DQAHist/TOF/h_ph",h_ph);
00404 m_thsvc->regGraph("/DQAHist/TOF/delT4",delT4);
00405 return StatusCode::SUCCESS;
00406
00407 }
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1.3.9.1