BOSS_7.0.2: /home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/DQA/DQA_TOF/DQA_TOF-00-00-07/src/DQA

00001 #include "GaudiKernel/MsgStream.h"
00002 #include "GaudiKernel/SmartDataPtr.h"
00003 #include "GaudiKernel/PropertyMgr.h"
00004 #include "GaudiKernel/Bootstrap.h"
00005 
00006 #include "GaudiKernel/INTupleSvc.h"
00007 #include "GaudiKernel/NTuple.h"
00008 #include "GaudiKernel/ITHistSvc.h"
00009 
00010 #include "CLHEP/Vector/ThreeVector.h"
00011 #include "CLHEP/Vector/LorentzVector.h"
00012 
00013 #include "EventModel/EventModel.h"
00014 #include "EventModel/Event.h"
00015 
00016 #include "EvtRecEvent/EvtRecEvent.h"
00017 #include "EvtRecEvent/EvtRecTrack.h"
00018 #include "DstEvent/TofHitStatus.h"
00019 #include "EventModel/EventHeader.h"
00020 #include "ParticleID/ParticleID.h"
00021 #include "DQAEvent/DQAEvent.h"
00022 #include "DQA_TOF/DQA_TOF.h"
00023 #include <vector>
00025 DQA_TOF::DQA_TOF(const std::string& name, ISvcLocator* pSvcLocator) :
00026         Algorithm(name, pSvcLocator) {
00027 
00028                 //Declare the properties  
00029         }
00030 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 
00031 StatusCode DQA_TOF::initialize(){
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         if(m_thsvc->regHist("/DQAHist/TOF/h_path",h_path).isFailure()){
00045             log << MSG::ERROR << "Couldn't register h_path" << endreq;
00046         }
00047         h_Bzrhit=new TH1F("h_Bzrhit","barrel z hitmap",240,-120,120);
00048         m_thsvc->regHist("/DQAHist/TOF/h_Bzrhit", h_Bzrhit);
00049         h_Ezrhit=new TH1F("h_Ezrhit","endcap z hitmap",50,40,90);
00050         m_thsvc->regHist("/DQAHist/TOF/h_Ezrhit", h_Ezrhit);
00051         h_ph=new TH1F("h_ph","barrel Q",900,0,9000);
00052         m_thsvc->regHist("/DQAHist/TOF/h_ph", h_ph);
00053         
00054         W_delT=new TH2F("W_delT","barrel west PMT delT",176,0,176,300,-1.5,1.5);
00055         m_thsvc->regHist("/DQAHist/TOF/W_delT", W_delT);
00056         E_delT=new TH2F("E_delT","barrel east PMT delT",176,0,176,300,-1.5,1.5);
00057         m_thsvc->regHist("/DQAHist/TOF/E_delT", E_delT);
00058         counter=new TH2F("counter","barrel counter delT",176,0,176,300,-1.5,1.5);
00059         m_thsvc->regHist("/DQAHist/TOF/counter", counter);
00060         cluster=new TH2F("cluster","barrel cluster delT",88,0,88,300,-1.5,1.5);
00061         m_thsvc->regHist("/DQAHist/TOF/cluster", cluster);
00062         EC_delT=new TH2F("EC_delT","endcap delT",96,0,96,300,-1.5,1.5);
00063         m_thsvc->regHist("/DQAHist/TOF/EC_delT", EC_delT);
00064         Bt_delT=new TH1F("Bt_delT","barrel delT",300,-1.5,1.5);
00065         m_thsvc->regHist("/DQAHist/TOF/Bt_delT", Bt_delT);
00066         Et_delT=new TH1F("Et_delT","endcap delT",300,-1.5,1.5);
00067         m_thsvc->regHist("/DQAHist/TOF/Et_delT", Et_delT);
00068         
00069         B_path=new TH2F("B_path","barrel flight distance vs z",240,-120,120,200,0.0,200.0);
00070         m_thsvc->regHist("/DQAHist/TOF/B_path", B_path);
00071         E_path=new TH2F("E_path","endcap path distance vs z",50,40,90,200,0.0,200.0);
00072         m_thsvc->regHist("/DQAHist/TOF/E_path", E_path);
00073         
00074         delT_z1=new TH2F("delT_z1","barrel east delT vs Z",240,-120,120,300,-1.5,1.5);
00075         m_thsvc->regHist("/DQAHist/TOF/delT_z1", delT_z1);
00076         delT_z2=new TH2F("delT_z2","barrel west delT vs Z",240,-120,120,300,-1.5,1.5);
00077         m_thsvc->regHist("/DQAHist/TOF/delT_z2", delT_z2);
00078         delT_z3=new TH2F("delT_z3","barrel counter delT vs Z",240,-120,120,300,-1.5,1.5);
00079         m_thsvc->regHist("/DQAHist/TOF/delT_z3", delT_z3);
00080         delT_z4=new TH2F("delT_z4","barrel cluster delT vs Z",240,-120,120,300,-1.5,1.5);
00081         m_thsvc->regHist("/DQAHist/TOF/delT_z4", delT_z4);
00082         
00083         W_delT_Q=new TH2F("W_delT_Q","west barrel delT vs Q",900,0,9000,300,-1.5,1.5);
00084         m_thsvc->regHist("/DQAHist/TOF/W_delT_Q", W_delT_Q);
00085         E_delT_Q=new TH2F("E_delT_Q","east barrel delT vs Q",900,0,9000,300,-1.5,1.5);
00086         m_thsvc->regHist("/DQAHist/TOF/E_delT_Q", E_delT_Q);
00087         
00088         delT_pp=new TH1F("delT_pp","proton delT",300,-1.5,1.5);
00089         m_thsvc->regHist("/DQAHist/TOF/delT_pp", delT_pp);
00090         delT_pm=new TH1F("delT_Pm","anti-proton delT",300,-1.5,1.5);
00091         m_thsvc->regHist("/DQAHist/TOF/delT_pm", delT_pm);
00092         delT_pi=new TH1F("delT_pi","pi delT",300,-1.5,1.5);
00093         m_thsvc->regHist("/DQAHist/TOF/delT_pi", delT_pi);
00094         delT_k=new TH1F("delT_k","k delT",300,-1.5,1.5);
00095         m_thsvc->regHist("/DQAHist/TOF/delT_k", delT_k);
00096         log << MSG::INFO << "DQA_TOF successfully return from initialize()" <<endmsg;
00097         return StatusCode::SUCCESS;
00098 }
00099 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 
00100 StatusCode DQA_TOF::execute() {
00101         MsgStream log(msgSvc(), name());
00102         log << MSG::INFO << "in execute()" << endreq;
00103 
00104         SmartDataPtr<Event::EventHeader> eventHeader(eventSvc(),"/Event/EventHeader");
00105         int runNo=eventHeader->runNumber();
00106         int event=eventHeader->eventNumber();
00107         log<<MSG::DEBUG<<"run,evtnum="<<runNo<<","<<event<<endreq;
00108 
00109         SmartDataPtr<DQAEvent::DQAEvent> dqaevt(eventSvc(), "/Event/DQATag");
00110         if ( !dqaevt ) {
00111                 log << MSG::ERROR << "Error accessing DQAEvent" << endreq;
00112                 return StatusCode::FAILURE;
00113         }
00114 
00115         log << MSG::DEBUG << "event tag = " << dqaevt->EventTag() << endreq;
00116 
00117         //if ( !dqaevt->Bhabha() )return StatusCode::SUCCESS; 
00118         SmartDataPtr<EvtRecEvent> evtRecEvent(eventSvc(), EventModel::EvtRec::EvtRecEvent);     
00119         SmartDataPtr<EvtRecTrackCol> evtRecTrkCol(eventSvc(),  EventModel::EvtRec::EvtRecTrackCol);
00120         if(dqaevt->Bhabha()){
00121                 for(int i = 0; i < evtRecEvent->totalCharged(); i++){
00122                         EvtRecTrackIterator itTrk=evtRecTrkCol->begin() + i;
00123                         //if ( (*itTrk)->partId() != 1 ) continue; // only e+, e-
00124                         if(!(*itTrk)->isElectron())continue;
00125                         if(!(*itTrk)->isTofTrackValid())continue;
00126                         SmartRefVector<RecTofTrack> tofTrkCol = (*itTrk)->tofTrack();
00127                         SmartRefVector<RecTofTrack>::iterator iter_tof = tofTrkCol.begin();
00128                         for(;iter_tof != tofTrkCol.end(); iter_tof++ ) {
00129                                 TofHitStatus *status = new TofHitStatus;
00130                                 status->setStatus((*iter_tof)->status());
00131                                 int tofid=(*iter_tof)->tofID();
00132                                 double tof=(*iter_tof)->tof();
00133                                 double texpE=(*iter_tof)->texpElectron();
00134                                 double path=(*iter_tof)->path();
00135                                 double Q=(*iter_tof)->ph();
00136                                 double zrhit=(*iter_tof)->zrhit();
00137 
00138                                 if(status->is_barrel()){
00139                                         h_Bzrhit->Fill(zrhit);
00140                                         h_path->Fill(path);
00141                                         //h_ph->Fill(Q);
00142                                         if(status->is_readout()){
00143                                                 h_ph->Fill(Q);
00144                                                 B_path->Fill(zrhit,path);
00145                                                 if(status->is_east()){                    //barrel readout east
00146                                                         E_delT->Fill(tofid,tof-texpE);
00147                                                         delT_z1->Fill(zrhit,tof-texpE);
00148                                                         E_delT_Q->Fill(Q,tof-texpE);
00149                                                 }
00150                                                 else {                                    //barrel readout west
00151                                                         W_delT->Fill(tofid,tof-texpE);
00152                                                         delT_z2->Fill(zrhit,tof-texpE);
00153                                                         W_delT_Q->Fill(Q,tof-texpE);
00154                                                 }
00155                                         }                                                 //end readout
00156                                         if(!status->is_readout()&&status->is_counter()){  //barrel counter
00157                                                 counter->Fill(tofid,tof-texpE);
00158                                                 delT_z3->Fill(zrhit,tof-texpE);
00159                                         }
00160                                         if(!status->is_counter()&&status->is_cluster()){  //barrel cluster
00161                                                 cluster->Fill(tofid,tof-texpE);
00162                                                 Bt_delT->Fill(tof-texpE);
00163                                                 delT_z4->Fill(zrhit,tof-texpE);
00164                                         }
00165                                 }
00166                                 else{                                                     //endcap
00167                                         E_path->Fill(zrhit,path);
00168                                         h_Ezrhit->Fill(zrhit);
00169                                         EC_delT->Fill(tofid,tof-texpE);
00170                                         Et_delT->Fill(tof-texpE);
00171                                 }
00172                         }                                                                 // loop 7 info fo tof
00173                 }                                                                         // loop every track
00174         }
00175         else{
00176                 for(int i = 0; i < evtRecEvent->totalCharged(); i++){
00177                         EvtRecTrackIterator itTrk=evtRecTrkCol->begin() + i;
00178                         //if ( (*itTrk)->partId() != 1 ) continue; // only e+, e-
00179                         if(!(*itTrk)->isTofTrackValid())continue;
00180                         SmartRefVector<RecTofTrack> tofTrkCol = (*itTrk)->tofTrack();
00181                         SmartRefVector<RecTofTrack>::iterator iter_tof = tofTrkCol.begin();
00182                         for(;iter_tof != tofTrkCol.end(); iter_tof++ ) {
00183                                 TofHitStatus *status = new TofHitStatus;
00184                                 status->setStatus((*iter_tof)->status());
00185                                 if(!status->is_cluster())continue;
00186                                 double tof=(*iter_tof)->tof();
00187                                 double texpPi=(*iter_tof)->texpPion();
00188                                 double texpK=(*iter_tof)->texpKaon();
00189                                 double texpP=(*iter_tof)->texpProton();
00190                                 //if((*itTrk)->partId()==3){  //pi
00191                                 if((*itTrk)->isPion()){
00192                                         double texpPi=(*iter_tof)->texpPion();
00193                                         delT_pi->Fill(tof-texpPi);
00194                                 }
00195                                 //else if((*itTrk)->partId()==4){ //k
00196                                 else if((*itTrk)->isKaon()){
00197                                         double texpK=(*iter_tof)->texpKaon();
00198                                         delT_k->Fill(tof-texpK);
00199                                 }
00200                                 //else if((*itTrk)->partId()==5){ //p
00201                                 else if((*itTrk)->isProton()){
00202                                         double texpP=(*iter_tof)->texpProton();
00203                                         if(!(*itTrk)->isMdcTrackValid())continue;
00204                                         RecMdcTrack *mdcTrk = (*itTrk)->mdcTrack();
00205                                         if(mdcTrk->charge()>0)delT_pp->Fill(tof-texpP);
00206                                         else delT_pm->Fill(tof-texpP);
00207                                 }
00208                         }
00209                 }
00210         }
00211         return StatusCode::SUCCESS;
00212 
00213 }
00214 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
00215 StatusCode DQA_TOF::finalize() {
00216         MsgStream log(msgSvc(), name());
00217         log << MSG::INFO << "in finalize()" << endmsg;
00218 
00219         return StatusCode::SUCCESS;
00220 }
00221 
00222
/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/DQA/DQA_TOF/DQA_TOF-00-00-07/src/DQA_TOF.cxx
