BOSS6.5.5: DQADtag Class Reference

00864                             {
00865 
00866   //  std::cout << "execute()" << std::endl;
00867 
00868   MsgStream log(msgSvc(), name());
00869   log << MSG::INFO << "in execute()" << endreq;
00870 
00871   setFilterPassed(false);
00872 
00873   SmartDataPtr<Event::EventHeader> eventHeader(eventSvc(),"/Event/EventHeader");
00874   int run=eventHeader->runNumber();
00875   int event=eventHeader->eventNumber();
00876   log << MSG::DEBUG <<"run, evtnum = "
00877         << run << " , "
00878         << event <<endreq;
00879 
00880   m_run = eventHeader->runNumber();  
00881   m_rec = eventHeader->eventNumber();
00882   m_evtrun = eventHeader->runNumber();  
00883   m_evtrec = eventHeader->eventNumber();
00884 
00886   // get beam energy
00888   if(m_ReadBeamEFromDB&&m_irun!=run){ 
00889         m_irun=run;
00890         if(m_readDb.isRunValid(m_irun))   
00891           m_beamE=m_readDb.getbeamE(m_irun);
00892   }
00893 
00894   double beamE=m_beamE;
00895   m_evtbeamE = m_beamE;
00896 
00897   if (m_debug) 
00898         cout<<"event "<<event<<endl;
00899   Ncut[0]++;
00900   if (Ncut[0]%10000 == 0)
00901   {
00902         cout<<"Ncut[0]: "<<Ncut[0]<<endl;
00903   }
00904   //if (Ncut[0]<120000)
00905   //return StatusCode::SUCCESS;
00906   SmartDataPtr<EvtRecEvent> evtRecEvent(eventSvc(), EventModel::EvtRec::EvtRecEvent);
00907   log << MSG::INFO << "get event tag OK" << endreq;
00908   log << MSG::DEBUG <<"ncharg, nneu, tottks = " 
00909         << evtRecEvent->totalCharged() << " , "
00910         << evtRecEvent->totalNeutral() << " , "
00911         << evtRecEvent->totalTracks() <<endreq;
00912 
00913   m_nch  = evtRecEvent->totalCharged();
00914   m_nneu = evtRecEvent->totalNeutral();
00915   m_evtnch  = evtRecEvent->totalCharged();
00916   m_evtnneu = evtRecEvent->totalNeutral();
00917 
00918   SmartDataPtr<EvtRecTrackCol> evtRecTrkCol(eventSvc(), EventModel::EvtRec::EvtRecTrackCol);
00919 
00920   //get primay vetex from db
00921   Hep3Vector xorigin(0,0,0);
00922   IVertexDbSvc*  vtxsvc;
00923   HepPoint3D vx(0., 0., 0.);
00924   HepSymMatrix Evx(3, 0);
00925   Gaudi::svcLocator()->service("VertexDbSvc", vtxsvc);
00926   if (vtxsvc->isVertexValid()) {
00927 
00928         double* vertex = vtxsvc->PrimaryVertex();     //vertex[0] = vx; vertex[1]= vy; vertex[2] = vz;
00929         double* vertexsigma = vtxsvc->SigmaPrimaryVertex();     
00930         xorigin.setX(vertex[0]);
00931         xorigin.setY(vertex[1]);
00932         xorigin.setZ(vertex[2]);
00933         vx[0] = vertex[0];
00934         vx[1] = vertex[1];
00935         vx[2] = vertex[2];
00936         Evx[0][0] = vertexsigma[0]*vertexsigma[0];
00937         Evx[1][1] = vertexsigma[1]*vertexsigma[1];
00938         Evx[2][2] = vertexsigma[2]*vertexsigma[2];
00939         m_privxpar.setVx(vx);
00940         m_privxpar.setEvx(Evx);
00941   }
00942 
00943   Vp4 pDtagCand_cms;
00944   pDtagCand_cms.clear();
00945   int numDtagCand = 0;
00946   vector<double> vmBC_D0, vmDCand_D0, vdeltE_D0, vpDCand_D0, vmpi0, vvmks, vmksp4,\
00947         vvlks, vmBC_DC, vmDCand_DC, vdeltE_DC, vpDCand_DC;
00948   vmBC_D0.clear();
00949   vmDCand_D0.clear();
00950   vdeltE_D0.clear();
00951   vpDCand_D0.clear();
00952   vmBC_DC.clear();
00953   vmDCand_DC.clear();
00954   vdeltE_DC.clear();
00955   vpDCand_DC.clear();
00956   vmpi0.clear();
00957   vvmks.clear();
00958   vmksp4.clear();
00959   vvlks.clear();
00960 
00961   SmartDataPtr<EvtRecDTagCol> evtRecDTagCol(eventSvc(), EventModel::EvtRec::EvtRecDTagCol);
00962   if ( ! evtRecDTagCol ) {
00963         cout << "Could not find EvtRecDTagCol" << endl;
00964         return StatusCode::SUCCESS;
00965   }
00966   //m_numAllMode = evtRecDTagCol->size();
00967 
00968   //use util tool to get dtaglist
00969   DTagTool dtagTool;
00970   if( dtagTool.isDTagListEmpty() ){
00971         if (m_debug) 
00972           cout<<"no D candidates found"<<endl;
00973         return StatusCode::SUCCESS;
00974   }
00975 
00976   DTagToolIterator iter_begin = dtagTool.modes_begin();
00977   DTagToolIterator iter_end =dtagTool.modes_end();
00978 
00979 
00980   int nCharge = 0;
00981 
00982   if (m_debug)
00983   {
00984         //test to get one mode only
00985         cout<<"test single mode search :"<< EvtRecDTag::kD0toKPi<<endl;  
00986 
00987         vector<int> mode = dtagTool.mode( EvtRecDTag::kD0toKPi );
00988         cout<<" there are "<< mode.size() <<" candidates for this  mode" <<endl;
00989 
00990         for( int i=0; i < mode.size(); i++){
00991 
00992           DTagToolIterator iter= dtagTool.modes_begin()+ mode[i];
00993           cout<<"No."<<i+1<<" candidate deltaE is : "<< (*iter)->deltaE()<<endl; 
00994 
00995         }
00996   }
00997 
00998   //loop the whole DTag lists, do all kinds of fit in side the loop
00999   const int numMode = 11;
01000   int indexD[numMode], indexTof[numMode], indexDedx[numMode],numDCand[numMode];
01001   for (int i=0;i!=numMode;i++) 
01002   {
01003         indexD[i] = 0;
01004         indexTof[i] = 0;
01005         indexDedx[i] = 0;
01006         numDCand[i] = 0;
01007   }
01008   int indexDAll = 0;
01009   int indexD0All = 0;
01010   int indexDCAll = 0;
01011 
01012   for (DTagToolIterator iter_dtag=iter_begin; iter_dtag != iter_end; iter_dtag++)
01013   {
01014         bool iswritten = 1;
01015         //decay mode name
01016         if (m_debug)
01017         {
01018           cout<<"***************"<<endl;
01019           cout<<"***************"<<endl;
01020           dtagTool<< iter_dtag;
01021           cout<<"***************"<<endl;
01022           cout<<"***************"<<endl;
01023         }
01024 
01025         // ***********************************************************************
01026         // m_checkDecayMode == 0 : allDecayMode         donot check pi0 Ks
01027         // m_checkDecayMode == 1 : singleDecayMode  D0 to KPi
01028         // m_checkDecayMode == 2 : fewDecayMode         10 decay modes  check pi0 Ks    
01029         // ***********************************************************************
01030         if (m_checkDecayMode < 3)               
01031         {
01032           //    general variables
01033           //HepLorentzVector p4=(*iter_dtag)->p4();
01034           //p4.boost(-0.011,0,0);
01035 
01036           //Hep3Vector p3=p4.v();
01037           //double deltaEDCand = (*iter_dtag)->deltaE();
01038           //double mBCDCand =(*iter_dtag)->mBC();  
01039           //if(fabs(deltaEDCand)>=0.1) continue;
01040           //if(mBCDCand<=1.8 || mBCDCand>=2.0) continue;
01041           m_dtagbeamE[indexDAll] = m_beamE;
01042           m_mode[indexDAll]=(*iter_dtag)->decayMode();
01043           m_charge[indexDAll]=(*iter_dtag)->charge();
01044           m_charm[indexDAll]=(*iter_dtag)->charm();
01045           m_numofchildren[indexDAll]=(*iter_dtag)->numOfChildren();
01046           m_mass[indexDAll]=(*iter_dtag)->mass();
01047           m_mBC[indexDAll]=(*iter_dtag)->mBC();
01048           m_e[indexDAll]=(*iter_dtag)->beamE();
01049           m_deltae[indexDAll]=(*iter_dtag)->deltaE();
01050           pDtagCand_cms.push_back((*iter_dtag)->p4());  
01051           SmartRefVector<EvtRecTrack> tracks=(*iter_dtag)->tracks();
01052           SmartRefVector<EvtRecTrack> othertracks=(*iter_dtag)->otherTracks();
01053           SmartRefVector<EvtRecTrack> othershowers=(*iter_dtag)->otherShowers();
01054           SmartRefVector<EvtRecTrack> showers=(*iter_dtag)->showers();
01055           m_nshower[indexDAll] = showers.size();
01056           m_ntrk[indexDAll]=tracks.size();
01057           m_nothertrk[indexDAll]=othertracks.size();
01058           m_nothershower[indexDAll] = othershowers.size();
01059 
01060           indexDAll++;
01061           if (m_checkDecayMode == 0 ) continue;
01062           //    end general variables
01063           if((*iter_dtag)->decayMode()==EvtRecDTag::kD0toKPi) {
01064                 const int MODE = 0;
01065                 fillDecayModeItem(iter_dtag,"D02KPi",indexD[MODE]);
01066                 fillTrackItem(tracks[0],"D02KPi_kaon",3,indexD[MODE],xorigin);
01067                 fillTrackItem(tracks[1],"D02KPi_pion",2,indexD[MODE],xorigin);
01068                 indexD[MODE]++;
01069                 if (m_checkDedx)
01070                 {
01071                   //if (!tracks[0]->isMdcDedxValid()) continue;
01072                   //if (!tracks[1]->isMdcDedxValid()) continue;
01073                   fillTrackDedxItem(tracks[0],"D02KPi_kaon",3,indexD[MODE]-1);
01074                   fillTrackDedxItem(tracks[1],"D02KPi_pion",2,indexD[MODE]-1);
01075                   //indexDedx[MODE]++;
01076                 }
01077                 if (m_checkTof)
01078                 {
01079                   //if(!tracks[0]->isTofTrackValid()) continue;
01080                   //if(!tracks[1]->isTofTrackValid()) continue;
01081                   fillTrackTofItem(tracks[0],"D02KPi_kaon",3,indexD[MODE]-1);
01082                   fillTrackTofItem(tracks[1],"D02KPi_pion",2,indexD[MODE]-1);
01083                   //indexTof[MODE]++;
01084                 }
01085           } //end of D0 to K Pi mode
01086           if (m_checkDecayMode == 1 ) continue;
01087           if((*iter_dtag)->decayMode()==EvtRecDTag::kD0toKPiPi0) {
01088                 //pi0 mass
01089                 const int MODE = 1;
01090                 DQAPi0Info pi0info;
01091                 pi0info.setchilds(showers[0],showers[1]);
01092                 if (!pi0info.calculate()) continue;
01093                 if (pi0info.m() >m_mpi0Max || pi0info.m() < m_mpi0Min) 
01094                   iswritten = 0;
01095                 vmpi0.push_back(pi0info.m());
01096                 fillPi0Item(showers,"D02KPiPi0",indexD[MODE],&pi0info );
01097                 fillDecayModeItem(iter_dtag,"D02KPiPi0",indexD[MODE]);
01098                 fillTrackItem(tracks[0],"D02KPiPi0_kaon",3,indexD[MODE],xorigin);
01099                 fillTrackItem(tracks[1],"D02KPiPi0_pion",2,indexD[MODE],xorigin);
01100                 fillPhotonItem(showers[0],"D02KPiPi0_photon1",indexD[MODE]);
01101                 fillPhotonItem(showers[1],"D02KPiPi0_photon2",indexD[MODE]);
01102                 indexD[MODE]++;
01103                 if (m_checkDedx)
01104                 {
01105                   //if (!tracks[0]->isMdcDedxValid()) continue;
01106                   //if (!tracks[1]->isMdcDedxValid()) continue;
01107                   fillTrackDedxItem(tracks[0],"D02KPiPi0_kaon",3,indexD[MODE]-1);
01108                   fillTrackDedxItem(tracks[1],"D02KPiPi0_pion",2,indexD[MODE]-1);
01109                   //indexDedx[MODE]++;
01110                 }
01111                 if (m_checkTof)
01112                 {
01113                   //if(!tracks[0]->isTofTrackValid()) continue;
01114                   //if(!tracks[1]->isTofTrackValid()) continue;
01115                   fillTrackTofItem(tracks[0],"D02KPiPi0_kaon",3,indexD[MODE]-1);
01116                   fillTrackTofItem(tracks[1],"D02KPiPi0_pion",2,indexD[MODE]-1);
01117                   //indexTof[MODE]++;
01118                 }
01119           }// end of D0 to KPiPi0
01120 
01121           //****************************************************
01122           //mode D0 to KPiPiPi only
01123           //****************************************************
01124           if((*iter_dtag)->decayMode()==EvtRecDTag::kD0toKPiPiPi) {
01125                 const int MODE = 2;
01126                 fillDecayModeItem(iter_dtag,"D02KPiPiPi",indexD[MODE]);
01127 
01128                 RecMdcKalTrack* mdcKalTrk1 = tracks[0]->mdcKalTrack();
01129                 RecMdcKalTrack* mdcKalTrk2 = tracks[1]->mdcKalTrack();
01130                 RecMdcKalTrack* mdcKalTrk3 = tracks[2]->mdcKalTrack();
01131                 RecMdcKalTrack* mdcKalTrk4 = tracks[3]->mdcKalTrack();
01132                 HepLorentzVector p4kaon = mdcKalTrk1->p4(xmass[3]);
01133                 HepLorentzVector p4pion1 = mdcKalTrk2->p4(xmass[2]);
01134                 HepLorentzVector p4pion2 = mdcKalTrk3->p4(xmass[2]);
01135                 HepLorentzVector p4pion3 = mdcKalTrk4->p4(xmass[2]);
01136                 m_D02KPiPiPi_kaon_pp[indexD[MODE]] = p4kaon.rho();
01137                 m_D02KPiPiPi_pion1_pp[indexD[MODE]] = p4pion1.rho();
01138                 m_D02KPiPiPi_pion2_pp[indexD[MODE]] = p4pion2.rho();
01139                 m_D02KPiPiPi_pion3_pp[indexD[MODE]] = p4pion3.rho();
01140                 m_D02KPiPiPi_kaon_Id[indexD[MODE]] = tracks[0]->trackId(); 
01141                 m_D02KPiPiPi_pion1_Id[indexD[MODE]] = tracks[1]->trackId();
01142                 m_D02KPiPiPi_pion2_Id[indexD[MODE]] = tracks[2]->trackId();
01143                 m_D02KPiPiPi_pion3_Id[indexD[MODE]] = tracks[3]->trackId();
01144                 indexD[MODE]++;
01145           }
01146           //****************************************************
01147           //mode D0 to KPiPiPiPi0 only
01148           //****************************************************
01149           if((*iter_dtag)->decayMode()==EvtRecDTag::kD0toKPiPiPiPi0) {
01150                 const int MODE = 3;
01151                 DQAPi0Info pi0info;
01152                 pi0info.setchilds(showers[0],showers[1]);
01153                 if (!pi0info.calculate()) continue;
01154                 if (pi0info.m() >m_mpi0Max || pi0info.m() < m_mpi0Min) 
01155                   iswritten = 0;
01156                 fillPi0Item(showers,"D02KPiPiPiPi0",indexD[MODE],&pi0info );
01157                 vmpi0.push_back(pi0info.m());
01158 
01159                 fillDecayModeItem(iter_dtag,"D02KPiPiPiPi0",indexD[MODE]);
01160 
01161                 RecMdcKalTrack* mdcKalTrk1 = tracks[0]->mdcKalTrack();
01162                 RecMdcKalTrack* mdcKalTrk2 = tracks[1]->mdcKalTrack();
01163                 RecMdcKalTrack* mdcKalTrk3 = tracks[2]->mdcKalTrack();
01164                 RecMdcKalTrack* mdcKalTrk4 = tracks[3]->mdcKalTrack();
01165                 HepLorentzVector p4kaon = mdcKalTrk1->p4(xmass[3]);
01166                 HepLorentzVector p4pion1 = mdcKalTrk2->p4(xmass[2]);
01167                 HepLorentzVector p4pion2 = mdcKalTrk3->p4(xmass[2]);
01168                 HepLorentzVector p4pion3 = mdcKalTrk4->p4(xmass[2]);
01169                 m_D02KPiPiPiPi0_kaon_pp[indexD[MODE]]  = p4kaon.rho();
01170                 m_D02KPiPiPiPi0_pion1_pp[indexD[MODE]] = p4pion1.rho();
01171                 m_D02KPiPiPiPi0_pion2_pp[indexD[MODE]] = p4pion2.rho();
01172                 m_D02KPiPiPiPi0_pion3_pp[indexD[MODE]] = p4pion3.rho();
01173                 m_D02KPiPiPiPi0_kaon_Id[indexD[MODE]] = tracks[0]->trackId(); 
01174                 m_D02KPiPiPiPi0_pion1_Id[indexD[MODE]] = tracks[1]->trackId();
01175                 m_D02KPiPiPiPi0_pion2_Id[indexD[MODE]] = tracks[2]->trackId();
01176                 m_D02KPiPiPiPi0_pion3_Id[indexD[MODE]] = tracks[3]->trackId();
01177                 indexD[MODE]++;
01178           }
01179           //****************************************************
01180           //mode D0 to KsPiPi only
01181           //****************************************************
01182           if((*iter_dtag)->decayMode()==EvtRecDTag::kD0toKsPiPi) {
01183                 const int MODE = 4;
01184                 DQAKsInfo ksinfo;
01185                 ksinfo.setchilds(tracks[0],tracks[1]);
01186                 ksinfo.setVtxPar(m_vxpar);
01187                 ksinfo.setPrimaryVtxPar(m_privxpar);
01188                 if (!ksinfo.calculate()) continue;
01189 
01190                 fillDecayModeItem(iter_dtag,"D02KsPiPi",indexD[MODE]);
01191                 fillTrackItem(tracks[0],"D02KsPiPi_pion1",2,indexD[MODE],xorigin);
01192                 fillTrackItem(tracks[1],"D02KsPiPi_pion2",2,indexD[MODE],xorigin);
01193                 fillTrackItem(tracks[2],"D02KsPiPi_pion3",2,indexD[MODE],xorigin);
01194                 fillTrackItem(tracks[3],"D02KsPiPi_pion4",2,indexD[MODE],xorigin);
01195                 fillKsItem(tracks,"D02KsPiPi",indexD[MODE],&ksinfo);
01196                 indexD[MODE]++;
01197 
01198                 if (iswritten) vvmks.push_back(ksinfo.m());
01199                 if (iswritten) vmksp4.push_back(ksinfo.mksp4());
01200                 if (iswritten) vvlks.push_back(ksinfo.decayLengthRatio());
01201 
01202           }// end of D0 to KsPiPi
01203           //****************************************************
01204           //mode D0 to KsPiPiPi0 only
01205           //****************************************************
01206           if((*iter_dtag)->decayMode()==EvtRecDTag::kD0toKsPiPiPi0) {
01207                 const int MODE = 5;
01208                 DQAKsInfo ksinfo;
01209                 ksinfo.setchilds(tracks[0],tracks[1]);
01210                 ksinfo.setVtxPar(m_vxpar);
01211                 ksinfo.setPrimaryVtxPar(m_privxpar);
01212                 if (!ksinfo.calculate()) continue;
01213                 DQAPi0Info pi0info;
01214                 pi0info.setchilds(showers[0],showers[1]);
01215                 if (!pi0info.calculate()) continue;
01216                 if (pi0info.m() >m_mpi0Max || pi0info.m() < m_mpi0Min) 
01217                   iswritten = 0;
01218                 fillPi0Item(showers,"D02KsPiPiPi0",indexD[MODE],&pi0info );
01219                 vmpi0.push_back(pi0info.m());
01220                 fillKsItem(tracks,"D02KsPiPiPi0",indexD[MODE],&ksinfo);
01221                 if (iswritten) vvmks.push_back(ksinfo.m());
01222                 if (iswritten) vmksp4.push_back(ksinfo.mksp4());
01223                 if (iswritten) vvlks.push_back(ksinfo.decayLengthRatio());
01224                 fillDecayModeItem(iter_dtag,"D02KsPiPiPi0",indexD[MODE]);
01225 
01226                 RecMdcKalTrack* mdcKalTrk1 = tracks[0]->mdcKalTrack();
01227                 RecMdcKalTrack* mdcKalTrk2 = tracks[1]->mdcKalTrack();
01228                 RecMdcKalTrack* mdcKalTrk3 = tracks[2]->mdcKalTrack();
01229                 RecMdcKalTrack* mdcKalTrk4 = tracks[3]->mdcKalTrack();
01230                 HepLorentzVector p4pion1 = mdcKalTrk1->p4(xmass[2]);
01231                 HepLorentzVector p4pion2 = mdcKalTrk2->p4(xmass[2]);
01232                 HepLorentzVector p4pion3 = mdcKalTrk3->p4(xmass[2]);
01233                 HepLorentzVector p4pion4 = mdcKalTrk4->p4(xmass[2]);
01234                 m_D02KsPiPiPi0_pion1_pp[indexD[MODE]] = p4pion1.rho();
01235                 m_D02KsPiPiPi0_pion2_pp[indexD[MODE]] = p4pion2.rho();
01236                 m_D02KsPiPiPi0_pion3_pp[indexD[MODE]] = p4pion3.rho();
01237                 m_D02KsPiPiPi0_pion4_pp[indexD[MODE]] = p4pion4.rho();
01238                 m_D02KsPiPiPi0_pion1_Id[indexD[MODE]] = tracks[0]->trackId();
01239                 m_D02KsPiPiPi0_pion2_Id[indexD[MODE]] = tracks[1]->trackId();
01240                 m_D02KsPiPiPi0_pion3_Id[indexD[MODE]] = tracks[2]->trackId();
01241                 m_D02KsPiPiPi0_pion4_Id[indexD[MODE]] = tracks[3]->trackId();
01242                 indexD[MODE]++;
01243           }// end of D0 to KsPiPiPi0
01244           //****************************************************
01245           //mode Dp to KPiPi only
01246           //****************************************************
01247           if((*iter_dtag)->decayMode()==EvtRecDTag::kDptoKPiPi) {
01248                 const int MODE = 6;
01249                 fillDecayModeItem(iter_dtag,"DC2KPiPi",indexD[MODE]);
01250                 fillTrackItem(tracks[0],"DC2KPiPi_kaon",3,indexD[MODE],xorigin);
01251                 fillTrackItem(tracks[1],"DC2KPiPi_pion1",2,indexD[MODE],xorigin);
01252                 fillTrackItem(tracks[2],"DC2KPiPi_pion2",2,indexD[MODE],xorigin);
01253                 indexD[MODE]++;
01254           }
01255 
01256           //****************************************************
01257           //mode Dp to KPiPiPi0 only
01258           //****************************************************
01259           if((*iter_dtag)->decayMode()==EvtRecDTag::kDptoKPiPiPi0) {
01260                 const int MODE = 7;
01261                 DQAPi0Info pi0info;
01262                 pi0info.setchilds(showers[0],showers[1]);
01263                 if (!pi0info.calculate()) continue;
01264                 if (pi0info.m() >m_mpi0Max || pi0info.m() < m_mpi0Min) 
01265                   iswritten = 0;
01266                 fillPi0Item(showers,"DC2KPiPiPi0",indexD[MODE],&pi0info );
01267                 vmpi0.push_back(pi0info.m());
01268                 fillDecayModeItem(iter_dtag,"DC2KPiPiPi0",indexD[MODE]);
01269 
01270                 fillTrackItem(tracks[0],"DC2KPiPiPi0_kaon",3,indexD[MODE],xorigin);
01271                 fillTrackItem(tracks[1],"DC2KPiPiPi0_pion1",2,indexD[MODE],xorigin);
01272                 fillTrackItem(tracks[2],"DC2KPiPiPi0_pion2",2,indexD[MODE],xorigin);
01273                 indexD[MODE]++;
01274           }
01275           //****************************************************
01276           //mode Dp to KsPiPi0 only
01277           //****************************************************
01278           if((*iter_dtag)->decayMode()==EvtRecDTag::kDptoKsPiPi0) {
01279                 const int MODE = 8;
01280                 DQAPi0Info pi0info;
01281                 pi0info.setchilds(showers[0],showers[1]);
01282                 if (!pi0info.calculate()) continue;
01283                 DQAKsInfo ksinfo;
01284                 ksinfo.setchilds(tracks[0],tracks[1]);
01285                 ksinfo.setVtxPar(m_vxpar);
01286                 ksinfo.setPrimaryVtxPar(m_privxpar);
01287                 if (!ksinfo.calculate()) continue;
01288                 if (pi0info.m() >m_mpi0Max || pi0info.m() < m_mpi0Min) 
01289                   iswritten = 0;
01290                 fillPi0Item(showers,"DC2KsPiPi0",indexD[MODE],&pi0info );
01291                 fillDecayModeItem(iter_dtag,"DC2KsPiPi0",indexD[MODE]);
01292                 fillTrackItem(tracks[0],"DC2KsPiPi0_pion1",2,indexD[MODE],xorigin);
01293                 fillTrackItem(tracks[1],"DC2KsPiPi0_pion2",2,indexD[MODE],xorigin);
01294                 fillTrackItem(tracks[2],"DC2KsPiPi0_pion3",2,indexD[MODE],xorigin);
01295                 fillKsItem(tracks,"DC2KsPiPi0",indexD[MODE],&ksinfo);
01296                 indexD[MODE]++;
01297                 if (iswritten) vvmks.push_back(ksinfo.m());
01298                 if (iswritten) vmksp4.push_back(ksinfo.mksp4());
01299                 if (iswritten) vvlks.push_back(ksinfo.decayLengthRatio());
01300           }// end of Dp to Kspipi0
01301           //****************************************************
01302           //mode Dp to KKPi only
01303           //****************************************************
01304           if((*iter_dtag)->decayMode()==EvtRecDTag::kDptoKsPiPiPi) {
01305                 const int MODE = 9;
01306                 DQAKsInfo ksinfo;
01307                 ksinfo.setchilds(tracks[0],tracks[1]);
01308                 ksinfo.setVtxPar(m_vxpar);
01309                 ksinfo.setPrimaryVtxPar(m_privxpar);
01310                 if (!ksinfo.calculate()) continue;
01311 
01312                 fillKsItem(tracks,"DC2KsPiPiPi",indexD[MODE],&ksinfo);
01313                 if (iswritten) vvmks.push_back(ksinfo.m());
01314                 if (iswritten) vmksp4.push_back(ksinfo.mksp4());
01315                 if (iswritten) vvlks.push_back(ksinfo.decayLengthRatio());
01316                 fillDecayModeItem(iter_dtag,"DC2KsPiPiPi",indexD[MODE]);
01317 
01318                 RecMdcKalTrack* mdcKalTrk1 = tracks[0]->mdcKalTrack();
01319                 RecMdcKalTrack* mdcKalTrk2 = tracks[1]->mdcKalTrack();
01320                 RecMdcKalTrack* mdcKalTrk3 = tracks[2]->mdcKalTrack();
01321                 RecMdcKalTrack* mdcKalTrk4 = tracks[3]->mdcKalTrack();
01322                 RecMdcKalTrack* mdcKalTrk5 = tracks[4]->mdcKalTrack();
01323                 HepLorentzVector p4pion1 = mdcKalTrk1->p4(xmass[2]);
01324                 HepLorentzVector p4pion2 = mdcKalTrk2->p4(xmass[2]);
01325                 HepLorentzVector p4pion3 = mdcKalTrk3->p4(xmass[2]);
01326                 HepLorentzVector p4pion4 = mdcKalTrk4->p4(xmass[2]);
01327                 HepLorentzVector p4pion5 = mdcKalTrk5->p4(xmass[2]);
01328                 m_DC2KsPiPiPi_pion1_pp[indexD[MODE]] = p4pion1.rho();
01329                 m_DC2KsPiPiPi_pion2_pp[indexD[MODE]] = p4pion2.rho();
01330                 m_DC2KsPiPiPi_pion3_pp[indexD[MODE]] = p4pion3.rho();
01331                 m_DC2KsPiPiPi_pion4_pp[indexD[MODE]] = p4pion4.rho();
01332                 m_DC2KsPiPiPi_pion5_pp[indexD[MODE]] = p4pion5.rho();
01333                 m_DC2KsPiPiPi_pion1_Id[indexD[MODE]] = tracks[0]->trackId();
01334                 m_DC2KsPiPiPi_pion2_Id[indexD[MODE]] = tracks[1]->trackId();
01335                 m_DC2KsPiPiPi_pion3_Id[indexD[MODE]] = tracks[2]->trackId();
01336                 m_DC2KsPiPiPi_pion4_Id[indexD[MODE]] = tracks[3]->trackId();
01337                 m_DC2KsPiPiPi_pion5_Id[indexD[MODE]] = tracks[4]->trackId();
01338                 indexD[MODE]++;
01339           }//   end of kDptoKsPiPiPi
01340           if((*iter_dtag)->decayMode()==EvtRecDTag::kDptoKKPi) {
01341                 const int MODE = 10;
01342                 fillDecayModeItem(iter_dtag,"DC2KKPi",indexD[MODE]);
01343 
01344                 RecMdcKalTrack* mdcKalTrk1 = tracks[0]->mdcKalTrack();
01345                 RecMdcKalTrack* mdcKalTrk2 = tracks[1]->mdcKalTrack();
01346                 RecMdcKalTrack* mdcKalTrk3 = tracks[2]->mdcKalTrack();
01347                 HepLorentzVector p4kaon1 = mdcKalTrk1->p4(xmass[3]);
01348                 HepLorentzVector p4kaon2 = mdcKalTrk2->p4(xmass[3]);
01349                 HepLorentzVector p4pion = mdcKalTrk3->p4(xmass[2]);
01350                 m_DC2KKPi_kaon1_pp[indexD[MODE]] = p4kaon1.rho();
01351                 m_DC2KKPi_kaon2_pp[indexD[MODE]] = p4kaon2.rho();
01352                 m_DC2KKPi_pion_pp[indexD[MODE]] = p4pion.rho();
01353                 m_DC2KKPi_kaon1_Id[indexD[MODE]] = tracks[0]->trackId();
01354                 m_DC2KKPi_kaon2_Id[indexD[MODE]] = tracks[1]->trackId();
01355                 m_DC2KKPi_pion_Id[indexD[MODE]] = tracks[2]->trackId();
01356                 indexD[MODE]++;
01357           }//   end of kDptoKKPi
01358           if (iswritten && (*iter_dtag)->charge() == 0)
01359           {
01360                 m_D0_mBC[indexD0All]=(*iter_dtag)->mBC();
01361                 m_D0_deltE[indexD0All]=(*iter_dtag)->deltaE();
01362                 m_D0_mass[indexD0All]=(*iter_dtag)->mass();
01363                 indexD0All++;
01364           }
01365           if (iswritten && (*iter_dtag)->charge() != 0)
01366           {
01367                 m_DC_mBC[indexDCAll]=(*iter_dtag)->mBC();
01368                 m_DC_deltE[indexDCAll]=(*iter_dtag)->deltaE();
01369                 m_DC_mass[indexDCAll]=(*iter_dtag)->mass();
01370                 indexDCAll++;
01371           }
01372         }
01373   }
01374   numDtagCand = pDtagCand_cms.size();
01375 
01376   dtagTool.clear();
01377   m_numAllMode = indexDAll;
01378   m_numD0Mode = indexD0All;
01379   m_numDCMode = indexDCAll;
01380   m_recNumD02KPi = indexD[0];
01381   m_recNumD02KPiPi0 = indexD[1]; 
01382   m_recNumD02KPiPiPi = indexD[2];
01383   m_recNumD02KPiPiPiPi0 = indexD[3];
01384   m_recNumD02KsPiPi = indexD[4];
01385   m_recNumD02KsPiPiPi0 = indexD[5];
01386   m_recNumDC2KPiPi = indexD[6];
01387   m_recNumDC2KPiPiPi0 = indexD[7];
01388   m_recNumDC2KsPiPi0 = indexD[8];
01389   m_recNumDC2KsPiPiPi = indexD[9];
01390   m_recNumDC2KKPi = indexD[10];
01391   //m_numD02KPiPi0Gam = indexD[1];
01392   //m_numD02KPi_Tof = indexTof[0];
01393   //m_numD02KPi_Dedx = indexDedx[0];
01394   //m_numD02KPiPi0_Tof = indexTof[1];
01395   //m_numD02KPiPi0_Dedx = indexDedx[1];
01396   Vint iGood;
01397   //  int nCharge = 0;
01398   for(int i = 0; i < evtRecEvent->totalCharged(); i++){
01399         EvtRecTrackIterator itTrk=evtRecTrkCol->begin() + i;
01400         if(!(*itTrk)->isMdcTrackValid()) continue;
01401         if (!(*itTrk)->isMdcKalTrackValid()) continue;
01402         RecMdcTrack *mdcTrk = (*itTrk)->mdcTrack();
01403         RecMdcKalTrack *mdcKalTrk = (*itTrk)->mdcKalTrack();    
01404 
01405         double pch =mdcTrk->p();
01406         double x0  =mdcTrk->x();
01407         double y0  =mdcTrk->y();
01408         double z0  =mdcTrk->z();
01409         double phi0=mdcTrk->helix(1);
01410         double xv=xorigin.x();
01411         double yv=xorigin.y();
01412         double Rxy=fabs((x0-xv)*cos(phi0)+(y0-yv)*sin(phi0));
01413         m_vtxTrkId = (*itTrk)->trackId();
01414         m_vx0 = x0;
01415         m_vy0 = y0;
01416         m_vz0 = z0-xorigin.z();
01417         m_vr0 = Rxy;
01418         m_Vct=cos(mdcTrk->theta());
01419         if (m_writeTuple && m_checkVtx)
01420         m_tuple2->write();
01421         if(fabs(m_vz0) >= m_vz0cut) continue;
01422         if(m_vr0 >= m_vr0cut) continue;
01423         if(fabs(m_Vct)>=m_cthcut) continue;
01424 
01425         iGood.push_back((*itTrk)->trackId());
01426         nCharge += mdcTrk->charge();
01427   }
01428   int nGood = iGood.size();
01429 
01433 
01434   if ( m_writeTuple && (m_checkDedx || m_checkTof))
01435   {
01436         Vint iGam;
01437         iGam.clear();
01438         for(int i = evtRecEvent->totalCharged(); i< evtRecEvent->totalTracks(); i++) {
01439           EvtRecTrackIterator itTrk=evtRecTrkCol->begin() + i;
01440           if(!(*itTrk)->isEmcShowerValid()) continue;
01441           RecEmcShower *emcTrk = (*itTrk)->emcShower();
01442           Hep3Vector emcpos(emcTrk->x(), emcTrk->y(), emcTrk->z());
01443           // find the nearest charged track
01444           double dthe = 200.;
01445           double dphi = 200.;
01446           double dang = 200.; 
01447           for(int j = 0; j < evtRecEvent->totalCharged(); j++) {
01448                 EvtRecTrackIterator jtTrk = evtRecTrkCol->begin() + j;
01449                 if(!(*jtTrk)->isExtTrackValid()) continue;
01450                 RecExtTrack *extTrk = (*jtTrk)->extTrack();
01451                 if(extTrk->emcVolumeNumber() == -1) continue;
01452                 Hep3Vector extpos = extTrk->emcPosition();
01453                 //      double ctht = extpos.cosTheta(emcpos);
01454                 double angd = extpos.angle(emcpos);
01455                 double thed = extpos.theta() - emcpos.theta();
01456                 double phid = extpos.deltaPhi(emcpos);
01457                 thed = fmod(thed+CLHEP::twopi+CLHEP::twopi+pi, CLHEP::twopi) - CLHEP::pi;
01458                 phid = fmod(phid+CLHEP::twopi+CLHEP::twopi+pi, CLHEP::twopi) - CLHEP::pi;
01459 
01460                 if(fabs(thed) < fabs(dthe)) dthe = thed;
01461                 if(fabs(phid) < fabs(dphi)) dphi = phid;
01462                 if(angd < dang) dang = angd;
01463           }
01464           if(dang>=200) continue;
01465           double eraw = emcTrk->energy();
01466           dthe = dthe * 180 / (CLHEP::pi);
01467           dphi = dphi * 180 / (CLHEP::pi);
01468           dang = dang * 180 / (CLHEP::pi);
01469           double m_dthe = dthe;
01470           double m_dphi = dphi;
01471           double m_dang = dang;
01472           double m_eraw = eraw;
01473           double theta=cos(emcpos.theta());
01474           if ((fabs(theta) > m_Barrel_th && fabs(theta) < m_Endcap_th_1) || fabs(theta) > m_Endcap_th_2 ) continue;
01475           //put different cut on barrel and endcap!
01476           if((emcTrk->module()==0 || emcTrk->module()==2)&&eraw < m_energyThreshold_e) continue;
01477           else if( emcTrk->module()==1 && eraw < m_energyThreshold_b) continue;
01478           //if(eraw < m_energyThreshold) continue;
01479           if(dang < m_gammaAngCut) continue;
01480           //
01481           // good photon cut will be set here
01482           //
01483           iGam.push_back((*itTrk)->trackId());
01484         }
01485 
01486         //
01487         // Finish Good Photon Selection
01488         //
01489         int nGam = iGam.size();
01490 
01491         //
01492         // Assign 4-momentum to each photon
01493         // 
01494 
01495         Vp4 pGam;
01496         pGam.clear();
01497         for(int i = 0; i < nGam; i++) {
01498           EvtRecTrackIterator itTrk = evtRecTrkCol->begin() + iGam[i]; 
01499           RecEmcShower* emcTrk = (*itTrk)->emcShower();
01500           double eraw = emcTrk->energy();
01501           double phi = emcTrk->phi();
01502           double the = emcTrk->theta();
01503           HepLorentzVector ptrk;
01504           ptrk.setPx(eraw*sin(the)*cos(phi));
01505           ptrk.setPy(eraw*sin(the)*sin(phi));
01506           ptrk.setPz(eraw*cos(the));
01507           ptrk.setE(eraw);
01508           //    ptrk = ptrk.boost(-0.011,0,0);// boost to cms
01509           pGam.push_back(ptrk);
01510         }
01511 
01512         //
01513         // check dedx infomation
01514         //
01515         if (m_checkDedx)
01516         {
01517           m_ngch = nGood;
01518           for(int ii = 0; ii < nGood; ii++) {
01519                 // dedx
01520                 m_ptrk[ii] = 999.0;
01521                 m_chie[ii] = 999.0;
01522                 m_chimu[ii] = 999.0;
01523                 m_chipi[ii] = 999.0;
01524                 m_chik[ii] = 999.0;
01525                 m_chip[ii] = 999.0;
01526                 m_ghit[ii] = 999.0;
01527                 m_thit[ii] = 999.0;
01528                 m_probPH[ii] = 1999.0;
01529                 m_normPH[ii] = 999.0;
01530           }
01531           int  indx0=0;
01532           for(int i = 0; i < nGood; i++) {
01533                 EvtRecTrackIterator  itTrk = evtRecTrkCol->begin() + iGood[i];
01534                 if(!(*itTrk)->isMdcTrackValid()) continue;
01535                 if(!(*itTrk)->isMdcDedxValid())continue;
01536                 RecMdcTrack* mdcTrk = (*itTrk)->mdcTrack();
01537                 RecMdcDedx* dedxTrk = (*itTrk)->mdcDedx();
01538                 m_trkiddedx[indx0] = (*itTrk)->trackId();
01539                 m_ptrk[indx0] = mdcTrk->p();
01540                 m_chie[indx0] = dedxTrk->chiE();
01541                 m_chimu[indx0] = dedxTrk->chiMu();
01542                 m_chipi[indx0] = dedxTrk->chiPi();
01543                 m_chik[indx0] = dedxTrk->chiK();
01544                 m_chip[indx0] = dedxTrk->chiP();
01545                 m_ghit[indx0] = dedxTrk->numGoodHits();
01546                 m_thit[indx0] = dedxTrk->numTotalHits();
01547                 m_probPH[indx0] = dedxTrk->probPH();
01548                 m_normPH[indx0] = dedxTrk->normPH();
01549                 indx0++;
01550           }
01551           m_tuple11->write();
01552         }  // end chechDedx
01553         //
01554         // check TOF infomation
01555         //
01556         if (m_checkTof)
01557         {
01558           m_ngch2 = nGood;
01559           for(int ii = 0; ii < nGood; ii++) {
01560                 //endtof
01561                 m_cntr_etof[ii] = 999.0;
01562                 m_ptot_etof[ii] = 999.0;
01563                 m_ph_etof[ii]   = 999.0;
01564                 m_rhit_etof[ii] = 999.0;
01565                 m_qual_etof[ii] = 999.0;
01566                 m_te_etof[ii]   = 999.0;
01567                 m_tmu_etof[ii]  = 999.0;
01568                 m_tpi_etof[ii]  = 999.0;
01569                 m_tk_etof[ii]   = 999.0;
01570                 m_tp_etof[ii]   = 999.0;
01571                 m_ec_tof[ii]    =  999.0;
01572                 m_ec_toff_e[ii]   = 999.0;
01573                 m_ec_toff_mu[ii]  = 999.0;
01574                 m_ec_toff_pi[ii]  = 999.0;
01575                 m_ec_toff_k[ii]   = 999.0;
01576                 m_ec_toff_p[ii]   = 999.0;
01577                 m_ec_tsig_e[ii]   = 999.0;
01578                 m_ec_tsig_mu[ii]  = 999.0;
01579                 m_ec_tsig_pi[ii]  = 999.0;
01580                 m_ec_tsig_k[ii]   = 999.0;
01581                 m_ec_tsig_p[ii]   = 999.0;
01582 
01583                 // barrel tof
01584                 m_cntr_btof1[ii]  = 999.0;
01585                 m_ptot_btof1[ii]  = 999.0;
01586                 m_ph_btof1[ii]    = 999.0;
01587                 m_zhit_btof1[ii]  = 999.0;
01588                 m_qual_btof1[ii]  = 999.0;
01589                 m_te_btof1[ii]    = 999.0;
01590                 m_tmu_btof1[ii]   = 999.0;
01591                 m_tpi_btof1[ii]   = 999.0;
01592                 m_tk_btof1[ii]    = 999.0;
01593                 m_tp_btof1[ii]    = 999.0;
01594                 m_b1_tof[ii]      = 999.0;
01595                 m_b1_toff_e[ii]   = 999.0;
01596                 m_b1_toff_mu[ii]  = 999.0;
01597                 m_b1_toff_pi[ii]  = 999.0;
01598                 m_b1_toff_k[ii]   = 999.0;
01599                 m_b1_toff_p[ii]   = 999.0;
01600                 m_b1_tsig_e[ii]   = 999.0;
01601                 m_b1_tsig_mu[ii]  = 999.0;
01602                 m_b1_tsig_pi[ii]  = 999.0;
01603                 m_b1_tsig_k[ii]   = 999.0;
01604                 m_b1_tsig_p[ii]   = 999.0;
01605                 //pid
01606                 //m_dedx_pid[ii] = 999.0;
01607                 //m_tof1_pid[ii] = 999.0;
01608                 //m_tof2_pid[ii] = 999.0;
01609                 //m_prob_pid[ii] = 999.0; 
01610                 //m_ptrk_pid[ii] = 999.0;  
01611                 //m_cost_pid[ii] = 999.0;
01612           }
01613           int indx1=0;
01614           for(int i = 0; i < nGood; i++) {
01615                 EvtRecTrackIterator  itTrk = evtRecTrkCol->begin() + iGood[i];
01616                 if(!(*itTrk)->isMdcTrackValid()) continue;
01617                 if(!(*itTrk)->isTofTrackValid()) continue;
01618 
01619                 RecMdcTrack * mdcTrk = (*itTrk)->mdcTrack();
01620                 SmartRefVector<RecTofTrack> tofTrkCol = (*itTrk)->tofTrack();
01621 
01622                 double ptrk = mdcTrk->p();
01623                 SmartRefVector<RecTofTrack>::iterator iter_tof = tofTrkCol.begin();
01624                 for(;iter_tof != tofTrkCol.end(); iter_tof++ ) { 
01625                   TofHitStatus *status = new TofHitStatus; 
01626                   status->setStatus((*iter_tof)->status());
01627                   if(!(status->is_barrel())){//endcap
01628                         if( !(status->is_counter()) ) continue; // ? 
01629                         if( status->layer()!=1 ) continue;//layer1
01630                         double path=(*iter_tof)->path(); // ? 
01631                         double tof  = (*iter_tof)->tof();
01632                         double ph   = (*iter_tof)->ph();
01633                         double rhit = (*iter_tof)->zrhit();
01634                         double qual = 0.0 + (*iter_tof)->quality();
01635                         double cntr = 0.0 + (*iter_tof)->tofID();
01636                         double texp[5];
01637                         double tsig[5];
01638                         for(int j = 0; j < 5; j++) {//0 e, 1 mu, 2 pi, 3 K, 4 p
01639                           texp[j] = (*iter_tof)->texp(j);
01640                           //           tsig[j] = (*iter_tof)->sigma(j);
01641                           //           toffset[j] = (*iter_tof)->offset(j);
01642                         }
01643                         m_trkidtof[indx1] = (*itTrk)->trackId();
01644                         m_cntr_etof[indx1]  = cntr;
01645                         m_ptot_etof[indx1] = ptrk;
01646                         m_ph_etof[indx1]   = ph;
01647                         m_rhit_etof[indx1]  = rhit;
01648                         m_qual_etof[indx1]  = qual;
01649                         m_te_etof[indx1]    = tof - texp[0];
01650                         m_tmu_etof[indx1]   = tof - texp[1];
01651                         m_tpi_etof[indx1]   = tof - texp[2];
01652                         m_tk_etof[indx1]    = tof - texp[3];
01653                         m_tp_etof[indx1]    = tof - texp[4];
01654 
01655                         m_ec_tof[indx1]      =  tof;
01656 
01657                         m_ec_toff_e[indx1]   = (*iter_tof)->toffset(0);
01658                         m_ec_toff_mu[indx1]  = (*iter_tof)->toffset(1);
01659                         m_ec_toff_pi[indx1]  = (*iter_tof)->toffset(2);
01660                         m_ec_toff_k[indx1]   = (*iter_tof)->toffset(3);
01661                         m_ec_toff_p[indx1]   = (*iter_tof)->toffset(4);
01662 
01663                         m_ec_tsig_e[indx1]   = (*iter_tof)->sigma(0);
01664                         m_ec_tsig_mu[indx1]  = (*iter_tof)->sigma(1);
01665                         m_ec_tsig_pi[indx1]  = (*iter_tof)->sigma(2);
01666                         m_ec_tsig_k[indx1]   = (*iter_tof)->sigma(3);
01667                         m_ec_tsig_p[indx1]   = (*iter_tof)->sigma(4);
01668 
01669                   }
01670                   else {//barrel
01671                         if( !(status->is_cluster()) ) continue; // ? 
01672                         double path=(*iter_tof)->path(); // ? 
01673                         double tof  = (*iter_tof)->tof();
01674                         double ph   = (*iter_tof)->ph();
01675                         double rhit = (*iter_tof)->zrhit();
01676                         double qual = 0.0 + (*iter_tof)->quality();
01677                         double cntr = 0.0 + (*iter_tof)->tofID();
01678                         double texp[5];
01679                         for(int j = 0; j < 5; j++) {
01680                           texp[j] = (*iter_tof)->texp(j);
01681                         }
01682                         m_cntr_btof1[indx1]  = cntr;
01683                         m_ptot_btof1[indx1] = ptrk;
01684                         m_ph_btof1[indx1]   = ph;
01685                         m_zhit_btof1[indx1]  = rhit;
01686                         m_qual_btof1[indx1]  = qual;
01687                         m_te_btof1[indx1]    = tof - texp[0];
01688                         m_tmu_btof1[indx1]   = tof - texp[1];
01689                         m_tpi_btof1[indx1]   = tof - texp[2];
01690                         m_tk_btof1[indx1]    = tof - texp[3];
01691                         m_tp_btof1[indx1]    = tof - texp[4];
01692 
01693                         m_b1_tof[indx1]      =  tof;
01694 
01695                         m_b1_toff_e[indx1]   = (*iter_tof)->toffset(0);
01696                         m_b1_toff_mu[indx1]  = (*iter_tof)->toffset(1);
01697                         m_b1_toff_pi[indx1]  = (*iter_tof)->toffset(2);
01698                         m_b1_toff_k[indx1]   = (*iter_tof)->toffset(3);
01699                         m_b1_toff_p[indx1]   = (*iter_tof)->toffset(4);
01700 
01701                         m_b1_tsig_e[indx1]   = (*iter_tof)->sigma(0);
01702                         m_b1_tsig_mu[indx1]  = (*iter_tof)->sigma(1);
01703                         m_b1_tsig_pi[indx1]  = (*iter_tof)->sigma(2);
01704                         m_b1_tsig_k[indx1]   = (*iter_tof)->sigma(3);
01705                         m_b1_tsig_p[indx1]   = (*iter_tof)->sigma(4);
01706 
01707                   }
01708                   delete status; 
01709                 } 
01710                 indx1++;
01711           } // loop all charged track
01712 
01713           m_tuple12->write();
01714           m_nggneu = nGam;
01715           int iphoton = 0;
01716           for (int i=0; i<nGam; i++)
01717           {
01718                 EvtRecTrackIterator itTrk=evtRecTrkCol->begin() + iGam[i];
01719                 if (!(*itTrk)->isEmcShowerValid()) continue; 
01720                 RecEmcShower *emcTrk = (*itTrk)->emcShower();
01721                 m_shwidtof[iphoton] = (*itTrk)->trackId();
01722                 m_numHits[iphoton] = emcTrk->numHits();
01723                 m_secondmoment[iphoton] = emcTrk->secondMoment();
01724                 m_x[iphoton] = emcTrk->x();
01725                 m_y[iphoton]= emcTrk->y();
01726                 m_z[iphoton]= emcTrk->z();
01727                 m_cosemc[iphoton]   = cos(emcTrk->theta());
01728                 m_phiemc[iphoton]   = emcTrk->phi();      
01729                 m_energy[iphoton]   = emcTrk->energy();
01730                 m_eSeed[iphoton]    = emcTrk->eSeed(); 
01731                 m_e3x3[iphoton]     = emcTrk->e3x3();  
01732                 m_e5x5[iphoton]     = emcTrk->e5x5(); 
01733                 m_lat[iphoton]      = emcTrk->latMoment();
01734                 m_a20[iphoton]      = emcTrk->a20Moment();
01735                 m_a42[iphoton]      = emcTrk->a42Moment();
01736                 iphoton++;
01737           }
01738           m_tuple13->write();
01739         }       //end checkTof
01740   }
01741 
01742   // DQAHist
01743   if (m_regHist)
01744   {
01745 
01746         string strdecaymode[11] = {"D02KPi","D02KPiPi0","D02KPiPiPi","D02KPiPiPiPi0","D02KsPiPi","D02KsPiPiPi0",\
01747           "DC2KPiPi","DC2KPiPiPi0","DC2KsPiPi0","DC2KsPiPiPi","DC2KKPi"};
01748         int pi0cut[11] = {0,1,0,1,0,1,0,1,1,0,0};
01749         int kscut[11] = {0,0,0,0,1,1,0,0,1,1,0};
01750         for (int i=0;i!=11;i++)
01751         {
01752           fillHistogram(m_thsvc,strdecaymode[i],indexD[i],3,m_decayModeHistName,m_tuple4,pi0cut[i],kscut[i]);
01753         }
01754         string strDdecaymode[2] = {"D0","DC"};
01755         fillHistogram(m_thsvc,strDdecaymode[0],indexD0All,3,m_decayModeHistName,m_tuple4);
01756         fillHistogram(m_thsvc,strDdecaymode[1],indexDCAll,3,m_decayModeHistName,m_tuple4);
01757         string strtrack[14] = {"D02KPi_kaon","D02KPi_pion","D02KPiPi0_kaon","D02KPiPi0_pion","D02KsPiPi_pion1",\
01758           "D02KsPiPi_pion2","D02KsPiPi_pion3","D02KsPiPi_pion4","DC2KPiPi_kaon","DC2KPiPi_pion1","DC2KPiPi_pion2",\
01759                 "DC2KsPiPi0_pion1","DC2KsPiPi0_pion2","DC2KsPiPi0_pion3" };
01760         int indexnum[14] = {indexD[0],indexD[0],indexD[1],indexD[1],indexD[4],indexD[4],indexD[4],indexD[4],indexD[6],\
01761           indexD[6],indexD[6],indexD[8],indexD[8],indexD[8]};
01762         for (int i=0;i!=14;i++)
01763         {
01764           fillHistogram(m_thsvc,strtrack[i],indexnum[i],10,m_trackHistName,m_tuple4);
01765         }
01766 
01767         string strtrackdedx[4] = {"D02KPi_kaon","D02KPi_pion","D02KPiPi0_kaon","D02KPiPi0_pion"};
01768         fillHistogram(m_thsvc,strtrackdedx[0],indexD[0],6,m_trackDedxHistName,m_tuple4);
01769         fillHistogram(m_thsvc,strtrackdedx[1],indexD[0],6,m_trackDedxHistName,m_tuple4);
01770         fillHistogram(m_thsvc,strtrackdedx[2],indexD[1],6,m_trackDedxHistName,m_tuple4);
01771         fillHistogram(m_thsvc,strtrackdedx[3],indexD[1],6,m_trackDedxHistName,m_tuple4);
01772         string strphoton[2] = {"D02KPiPi0_photon1","D02KPiPi0_photon2"};
01773         fillHistogram(m_thsvc,strphoton[0],indexD[1],5,m_photonHistName,m_tuple4);
01774         fillHistogram(m_thsvc,strphoton[1],indexD[1],5,m_photonHistName,m_tuple4);
01775         string strtracktof[4] = {"D02KPi_kaon","D02KPi_pion","D02KPiPi0_kaon","D02KPiPi0_pion"};
01776         fillTofHistogram(m_thsvc,strtrackdedx[0],indexD[0]);
01777         fillTofHistogram(m_thsvc,strtrackdedx[1],indexD[0]);
01778         fillTofHistogram(m_thsvc,strtrackdedx[2],indexD[1]);
01779         fillTofHistogram(m_thsvc,strtrackdedx[3],indexD[1]);
01780         string strpi0[2] = {"D02KPiPi0","DC2KsPiPi0"};
01781         fillHistogram(m_thsvc,strpi0[0],indexD[1],5,m_Pi0HistName,m_tuple4);
01782         fillHistogram(m_thsvc,strpi0[1],indexD[8],5,m_Pi0HistName,m_tuple4);
01783         string strks[2] = {"D02KsPiPi_ks","DC2KsPiPi0_ks"};
01784         fillHistogram(m_thsvc,strks[0],indexD[4],4,m_KsHistName,m_tuple4);
01785         fillHistogram(m_thsvc,strks[1],indexD[8],4,m_KsHistName,m_tuple4);
01786 
01787   }
01788   if (m_writeTuple)
01789   {
01790         m_tuple4->write();
01791   }
01792   if (m_writeTuple && m_checkEvt)
01793   {
01794         m_tuple5->write();
01795   }
01796 //if (m_checkDedx)
01797 //m_tuple16->write();
01798 //if (m_checkTof)
01799 //m_tuple17->write();
01800 
01801 setFilterPassed(true);
01802 
01803 return StatusCode::SUCCESS;
01804 }

00174                               {
00175   MsgStream log(msgSvc(), name());
00176         
00177   log << MSG::INFO << "in initialize()" << endmsg;
00178 
00179   m_irun=-100;
00180 
00181   StatusCode status;
00182 
00183   NTuplePtr nt2(ntupleSvc(), "DQAFILE/CheckVtx");
00184   if ( nt2 ) m_tuple2 = nt2;
00185   else {
00186         m_tuple2 = ntupleSvc()->book ("DQAFILE/CheckVtx", CLID_ColumnWiseTuple, "ks N-Tuple example");
00187         if ( m_tuple2 )    {
00188           status = m_tuple2->addItem ("vtxTrkId",m_vtxTrkId);
00189           status = m_tuple2->addItem ("vx0",   m_vx0);
00190           status = m_tuple2->addItem ("vy0",   m_vy0);
00191           status = m_tuple2->addItem ("vz0",   m_vz0);
00192           status = m_tuple2->addItem ("vr0",   m_vr0);
00193           status = m_tuple2->addItem ("vct",   m_Vct);
00194         }
00195         else    { 
00196           log << MSG::ERROR << "    Cannot book N-tuple:" << long(m_tuple2) << endmsg;
00197           return StatusCode::FAILURE;
00198         }
00199   }
00200 
00201   NTuplePtr nt4(ntupleSvc(), "DQAFILE/Dtag");
00202   if ( nt4 ) m_tuple4 = nt4;
00203   else {
00204         m_tuple4 = ntupleSvc()->book ("DQAFILE/Dtag", CLID_ColumnWiseTuple, "ks N-Tuple example");
00205         if ( m_tuple4 )    {
00206           status = m_tuple4->addItem ("run",    m_run);
00207           status = m_tuple4->addItem ("rec",    m_rec); 
00208           status = m_tuple4->addItem ("nch",    m_nch);
00209           status = m_tuple4->addItem ("nneu",   m_nneu);
00210           status = m_tuple4->addItem ("numD0Mode",              m_numD0Mode, 0 ,10000);
00211           status = m_tuple4->addIndexedItem ("D0_mBC",                  m_numD0Mode, m_D0_mBC);
00212           status = m_tuple4->addIndexedItem ("D0_deltE",                m_numD0Mode,m_D0_deltE);
00213           status = m_tuple4->addIndexedItem ("D0_mass",         m_numD0Mode,m_D0_mass);
00214           status = m_tuple4->addItem ("numDCMode",              m_numDCMode, 0 ,10000);
00215           status = m_tuple4->addIndexedItem ("DC_mBC",                  m_numDCMode, m_DC_mBC);
00216           status = m_tuple4->addIndexedItem ("DC_deltE",                m_numDCMode,m_DC_deltE);
00217           status = m_tuple4->addIndexedItem ("DC_mass",         m_numDCMode,m_DC_mass);
00218           status = m_tuple4->addItem ("numAllMode",             m_numAllMode, 0 ,10000);
00219           status = m_tuple4->addIndexedItem ("DtagbeamE",       m_numAllMode,m_dtagbeamE);
00220           status = m_tuple4->addIndexedItem ("mode",            m_numAllMode, m_mode);
00221           status = m_tuple4->addIndexedItem ("charge",          m_numAllMode,   m_charge);
00222           status = m_tuple4->addIndexedItem ("charm",           m_numAllMode,    m_charm);
00223           status = m_tuple4->addIndexedItem ("numchildren",     m_numAllMode, m_numofchildren);
00224           status = m_tuple4->addIndexedItem ("mass",            m_numAllMode, m_mass);
00225           status = m_tuple4->addIndexedItem ("mBC",                     m_numAllMode, m_mBC);
00226           status = m_tuple4->addIndexedItem ("deltaE",          m_numAllMode,m_deltae);
00227           status = m_tuple4->addIndexedItem ("E",                       m_numAllMode, m_e);
00228           status = m_tuple4->addIndexedItem ("notherTrk",       m_numAllMode, m_nothertrk);
00229           status = m_tuple4->addIndexedItem ("nTrk",            m_numAllMode,m_ntrk);
00230           status = m_tuple4->addIndexedItem ("nothershower",m_numAllMode, m_nothershower);
00231           status = m_tuple4->addIndexedItem ("nshower",         m_numAllMode,m_nshower);
00232           status = m_tuple4->addItem ( "recNumD02KPi",                  m_recNumD02KPi, 0,500);
00233           status = m_tuple4->addIndexedItem ("D02KPi_e",                m_recNumD02KPi, m_D02KPi_e);
00234           status = m_tuple4->addIndexedItem ("D02KPi_p",                m_recNumD02KPi, m_D02KPi_p);
00235           status = m_tuple4->addIndexedItem ("D02KPi_mBC",      m_recNumD02KPi, m_D02KPi_mBC);
00236           status = m_tuple4->addIndexedItem ("D02KPi_mass",     m_recNumD02KPi, m_D02KPi_mass);
00237           status = m_tuple4->addIndexedItem ("D02KPi_deltE",    m_recNumD02KPi, m_D02KPi_deltE);
00238           status = m_tuple4->addIndexedItem ("D02KPi_charm",    m_recNumD02KPi, m_D02KPi_charm);
00239           status = m_tuple4->addIndexedItem ("D02KPi_kaon_pp",          m_recNumD02KPi, m_D02KPi_kaon_pp);
00240           status = m_tuple4->addIndexedItem ("D02KPi_pion_pp",          m_recNumD02KPi, m_D02KPi_pion_pp);
00241           status = m_tuple4->addIndexedItem ("D02KPi_kaon_vx",          m_recNumD02KPi, m_D02KPi_kaon_vx);
00242           status = m_tuple4->addIndexedItem ("D02KPi_kaon_vy",          m_recNumD02KPi, m_D02KPi_kaon_vy);
00243           status = m_tuple4->addIndexedItem ("D02KPi_kaon_vz",          m_recNumD02KPi, m_D02KPi_kaon_vz);
00244           status = m_tuple4->addIndexedItem ("D02KPi_kaon_vr0",         m_recNumD02KPi, m_D02KPi_kaon_vr0);
00245           status = m_tuple4->addIndexedItem ("D02KPi_kaon_px",          m_recNumD02KPi, m_D02KPi_kaon_px);
00246           status = m_tuple4->addIndexedItem ("D02KPi_kaon_py",          m_recNumD02KPi, m_D02KPi_kaon_py);
00247           status = m_tuple4->addIndexedItem ("D02KPi_kaon_pz",          m_recNumD02KPi, m_D02KPi_kaon_pz);
00248           status = m_tuple4->addIndexedItem ("D02KPi_kaon_E",           m_recNumD02KPi, m_D02KPi_kaon_E);
00249           status = m_tuple4->addIndexedItem ("D02KPi_kaon_cos",         m_recNumD02KPi, m_D02KPi_kaon_cos);
00250           status = m_tuple4->addIndexedItem ("D02KPi_pion_vx",          m_recNumD02KPi, m_D02KPi_pion_vx);
00251           status = m_tuple4->addIndexedItem ("D02KPi_pion_vy",          m_recNumD02KPi, m_D02KPi_pion_vy);
00252           status = m_tuple4->addIndexedItem ("D02KPi_pion_vz",          m_recNumD02KPi, m_D02KPi_pion_vz);
00253           status = m_tuple4->addIndexedItem ("D02KPi_pion_vr0",         m_recNumD02KPi, m_D02KPi_pion_vr0);
00254           status = m_tuple4->addIndexedItem ("D02KPi_pion_px",          m_recNumD02KPi, m_D02KPi_pion_px);
00255           status = m_tuple4->addIndexedItem ("D02KPi_pion_py",          m_recNumD02KPi, m_D02KPi_pion_py);
00256           status = m_tuple4->addIndexedItem ("D02KPi_pion_pz",          m_recNumD02KPi, m_D02KPi_pion_pz);
00257           status = m_tuple4->addIndexedItem ("D02KPi_pion_E",           m_recNumD02KPi, m_D02KPi_pion_E);
00258           status = m_tuple4->addIndexedItem ("D02KPi_pion_cos",         m_recNumD02KPi, m_D02KPi_pion_cos);
00259           status = m_tuple4->addIndexedItem ("D02KPi_kaon_Id",          m_recNumD02KPi, m_D02KPi_kaon_Id);
00260           status = m_tuple4->addIndexedItem ("D02KPi_pion_Id",          m_recNumD02KPi, m_D02KPi_pion_Id);
00261           status = m_tuple4->addItem ( "recNumD02KPiPi0",               m_recNumD02KPiPi0, 0,500);
00262           status = m_tuple4->addIndexedItem ("D02KPiPi0_e",  m_recNumD02KPiPi0, m_D02KPiPi0_e);
00263           status = m_tuple4->addIndexedItem ("D02KPiPi0_p",     m_recNumD02KPiPi0, m_D02KPiPi0_p);
00264           status = m_tuple4->addIndexedItem ("D02KPiPi0_mBC",   m_recNumD02KPiPi0, m_D02KPiPi0_mBC);
00265           status = m_tuple4->addIndexedItem ("D02KPiPi0_mass",  m_recNumD02KPiPi0, m_D02KPiPi0_mass);
00266           status = m_tuple4->addIndexedItem ("D02KPiPi0_deltE", m_recNumD02KPiPi0, m_D02KPiPi0_deltE);
00267           status = m_tuple4->addIndexedItem ("D02KPiPi0_charm", m_recNumD02KPiPi0, m_D02KPiPi0_charm);
00268           status = m_tuple4->addIndexedItem ("D02KPiPi0_kaon_pp", m_recNumD02KPiPi0, m_D02KPiPi0_kaon_pp);
00269           status = m_tuple4->addIndexedItem ("D02KPiPi0_pion_pp", m_recNumD02KPiPi0, m_D02KPiPi0_pion_pp);
00270           status = m_tuple4->addIndexedItem ("D02KPiPi0_kaon_Id", m_recNumD02KPiPi0, m_D02KPiPi0_kaon_Id);
00271           status = m_tuple4->addIndexedItem ("D02KPiPi0_pion_Id", m_recNumD02KPiPi0, m_D02KPiPi0_pion_Id);
00272           status = m_tuple4->addIndexedItem ("D02KPiPi0_kaon_vx",       m_recNumD02KPiPi0,      m_D02KPiPi0_kaon_vx);
00273           status = m_tuple4->addIndexedItem ("D02KPiPi0_kaon_vy",       m_recNumD02KPiPi0,      m_D02KPiPi0_kaon_vy);
00274           status = m_tuple4->addIndexedItem ("D02KPiPi0_kaon_vz",       m_recNumD02KPiPi0,      m_D02KPiPi0_kaon_vz);
00275           status = m_tuple4->addIndexedItem ("D02KPiPi0_kaon_vr0",      m_recNumD02KPiPi0,      m_D02KPiPi0_kaon_vr0);
00276           status = m_tuple4->addIndexedItem ("D02KPiPi0_kaon_px",       m_recNumD02KPiPi0,      m_D02KPiPi0_kaon_px);
00277           status = m_tuple4->addIndexedItem ("D02KPiPi0_kaon_py",       m_recNumD02KPiPi0,      m_D02KPiPi0_kaon_py);
00278           status = m_tuple4->addIndexedItem ("D02KPiPi0_kaon_pz",       m_recNumD02KPiPi0,      m_D02KPiPi0_kaon_pz);
00279           status = m_tuple4->addIndexedItem ("D02KPiPi0_kaon_E",        m_recNumD02KPiPi0,      m_D02KPiPi0_kaon_E);
00280           status = m_tuple4->addIndexedItem ("D02KPiPi0_kaon_cos",      m_recNumD02KPiPi0,      m_D02KPiPi0_kaon_cos);
00281           status = m_tuple4->addIndexedItem ("D02KPiPi0_pion_vx",       m_recNumD02KPiPi0,      m_D02KPiPi0_pion_vx);
00282           status = m_tuple4->addIndexedItem ("D02KPiPi0_pion_vy",       m_recNumD02KPiPi0,      m_D02KPiPi0_pion_vy);
00283           status = m_tuple4->addIndexedItem ("D02KPiPi0_pion_vz",       m_recNumD02KPiPi0,      m_D02KPiPi0_pion_vz);
00284           status = m_tuple4->addIndexedItem ("D02KPiPi0_pion_vr0",      m_recNumD02KPiPi0,      m_D02KPiPi0_pion_vr0);
00285           status = m_tuple4->addIndexedItem ("D02KPiPi0_pion_px",       m_recNumD02KPiPi0,      m_D02KPiPi0_pion_px);
00286           status = m_tuple4->addIndexedItem ("D02KPiPi0_pion_py",       m_recNumD02KPiPi0,      m_D02KPiPi0_pion_py);
00287           status = m_tuple4->addIndexedItem ("D02KPiPi0_pion_pz",       m_recNumD02KPiPi0,      m_D02KPiPi0_pion_pz);
00288           status = m_tuple4->addIndexedItem ("D02KPiPi0_pion_E",        m_recNumD02KPiPi0,      m_D02KPiPi0_pion_E);
00289           status = m_tuple4->addIndexedItem ("D02KPiPi0_pion_cos",  m_recNumD02KPiPi0,  m_D02KPiPi0_pion_cos);
00290           status = m_tuple4->addIndexedItem ("D02KPiPi0_pi0mass", m_recNumD02KPiPi0, m_D02KPiPi0_pi0mass);
00291           status = m_tuple4->addIndexedItem ("D02KPiPi0_photon1_Id", m_recNumD02KPiPi0, m_D02KPiPi0_photon1_Id);
00292           status = m_tuple4->addIndexedItem ("D02KPiPi0_photon2_Id", m_recNumD02KPiPi0, m_D02KPiPi0_photon2_Id);
00293           status = m_tuple4->addIndexedItem ("D02KPiPi0_photon1_E", m_recNumD02KPiPi0, m_D02KPiPi0_photon1_E);
00294           status = m_tuple4->addIndexedItem ("D02KPiPi0_photon2_E", m_recNumD02KPiPi0, m_D02KPiPi0_photon2_E);
00295           status = m_tuple4->addIndexedItem ("D02KPiPi0_photon1_cos", m_recNumD02KPiPi0, m_D02KPiPi0_photon1_cos);
00296           status = m_tuple4->addIndexedItem ("D02KPiPi0_photon2_cos", m_recNumD02KPiPi0, m_D02KPiPi0_photon2_cos);
00297           status = m_tuple4->addItem ( "recNumD02KPiPiPi",                      m_recNumD02KPiPiPi, 0,500);
00298           status = m_tuple4->addIndexedItem ("D02KPiPiPi_e",            m_recNumD02KPiPiPi, m_D02KPiPiPi_e);
00299           status = m_tuple4->addIndexedItem ("D02KPiPiPi_p",            m_recNumD02KPiPiPi, m_D02KPiPiPi_p);
00300           status = m_tuple4->addIndexedItem ("D02KPiPiPi_mBC",          m_recNumD02KPiPiPi, m_D02KPiPiPi_mBC);
00301           status = m_tuple4->addIndexedItem ("D02KPiPiPi_mass",         m_recNumD02KPiPiPi, m_D02KPiPiPi_mass);
00302           status = m_tuple4->addIndexedItem ("D02KPiPiPi_deltE",        m_recNumD02KPiPiPi, m_D02KPiPiPi_deltE);
00303           status = m_tuple4->addIndexedItem ("D02KPiPiPi_charm",        m_recNumD02KPiPiPi, m_D02KPiPiPi_charm);
00304           status = m_tuple4->addIndexedItem ("D02KPiPiPi_kaon_pp",      m_recNumD02KPiPiPi, m_D02KPiPiPi_kaon_pp);
00305           status = m_tuple4->addIndexedItem ("D02KPiPiPi_pion1_pp",     m_recNumD02KPiPiPi, m_D02KPiPiPi_pion1_pp);
00306           status = m_tuple4->addIndexedItem ("D02KPiPiPi_pion2_pp",     m_recNumD02KPiPiPi, m_D02KPiPiPi_pion2_pp);
00307           status = m_tuple4->addIndexedItem ("D02KPiPiPi_pion3_pp",     m_recNumD02KPiPiPi, m_D02KPiPiPi_pion3_pp);
00308           status = m_tuple4->addIndexedItem ("D02KPiPiPi_kaon_Id",      m_recNumD02KPiPiPi, m_D02KPiPiPi_kaon_Id);
00309           status = m_tuple4->addIndexedItem ("D02KPiPiPi_pion1_Id",     m_recNumD02KPiPiPi, m_D02KPiPiPi_pion1_Id);
00310           status = m_tuple4->addIndexedItem ("D02KPiPiPi_pion2_Id",     m_recNumD02KPiPiPi, m_D02KPiPiPi_pion2_Id);
00311           status = m_tuple4->addIndexedItem ("D02KPiPiPi_pion3_Id",     m_recNumD02KPiPiPi, m_D02KPiPiPi_pion3_Id);
00312           status = m_tuple4->addItem ( "recNumD02KPiPiPiPi0",                   m_recNumD02KPiPiPiPi0,0,500);
00313           status = m_tuple4->addIndexedItem ("D02KPiPiPiPi0_e",                 m_recNumD02KPiPiPiPi0,m_D02KPiPiPiPi0_e);
00314           status = m_tuple4->addIndexedItem ("D02KPiPiPiPi0_p",                 m_recNumD02KPiPiPiPi0,m_D02KPiPiPiPi0_p);
00315           status = m_tuple4->addIndexedItem ("D02KPiPiPiPi0_mBC",       m_recNumD02KPiPiPiPi0,m_D02KPiPiPiPi0_mBC);
00316           status = m_tuple4->addIndexedItem ("D02KPiPiPiPi0_mass",      m_recNumD02KPiPiPiPi0,m_D02KPiPiPiPi0_mass);
00317           status = m_tuple4->addIndexedItem ("D02KPiPiPiPi0_deltE",     m_recNumD02KPiPiPiPi0,m_D02KPiPiPiPi0_deltE);
00318           status = m_tuple4->addIndexedItem ("D02KPiPiPiPi0_charm",             m_recNumD02KPiPiPiPi0,m_D02KPiPiPiPi0_charm);
00319           status = m_tuple4->addIndexedItem ("D02KPiPiPiPi0_kaon_pp",           m_recNumD02KPiPiPiPi0,m_D02KPiPiPiPi0_kaon_pp);
00320           status = m_tuple4->addIndexedItem ("D02KPiPiPiPi0_pion1_pp",  m_recNumD02KPiPiPiPi0,m_D02KPiPiPiPi0_pion1_pp);
00321           status = m_tuple4->addIndexedItem ("D02KPiPiPiPi0_pion2_pp",  m_recNumD02KPiPiPiPi0,m_D02KPiPiPiPi0_pion2_pp);
00322           status = m_tuple4->addIndexedItem ("D02KPiPiPiPi0_pion3_pp",  m_recNumD02KPiPiPiPi0,m_D02KPiPiPiPi0_pion3_pp);
00323           status = m_tuple4->addIndexedItem ("D02KPiPiPiPi0_kaon_Id",           m_recNumD02KPiPiPiPi0,m_D02KPiPiPiPi0_kaon_Id);
00324           status = m_tuple4->addIndexedItem ("D02KPiPiPiPi0_pion1_Id",  m_recNumD02KPiPiPiPi0,m_D02KPiPiPiPi0_pion1_Id);
00325           status = m_tuple4->addIndexedItem ("D02KPiPiPiPi0_pion2_Id",  m_recNumD02KPiPiPiPi0,m_D02KPiPiPiPi0_pion2_Id);
00326           status = m_tuple4->addIndexedItem ("D02KPiPiPiPi0_pion3_Id",  m_recNumD02KPiPiPiPi0,m_D02KPiPiPiPi0_pion3_Id);
00327           status = m_tuple4->addIndexedItem ("D02KPiPiPiPi0_pi0mass",   m_recNumD02KPiPiPiPi0,m_D02KPiPiPiPi0_pi0mass);
00328           status = m_tuple4->addIndexedItem ("D02KPiPiPiPi0_photon1_Id",        m_recNumD02KPiPiPiPi0,m_D02KPiPiPiPi0_photon1_Id);
00329           status = m_tuple4->addIndexedItem ("D02KPiPiPiPi0_photon2_Id",        m_recNumD02KPiPiPiPi0,m_D02KPiPiPiPi0_photon2_Id);
00330           status = m_tuple4->addIndexedItem ("D02KPiPiPiPi0_photon1_E", m_recNumD02KPiPiPiPi0,m_D02KPiPiPiPi0_photon1_E);
00331           status = m_tuple4->addIndexedItem ("D02KPiPiPiPi0_photon2_E", m_recNumD02KPiPiPiPi0,m_D02KPiPiPiPi0_photon2_E);
00332           status = m_tuple4->addIndexedItem ("D02KPiPiPiPi0_photon1_cos",       m_recNumD02KPiPiPiPi0,m_D02KPiPiPiPi0_photon1_cos);
00333           status = m_tuple4->addIndexedItem ("D02KPiPiPiPi0_photon2_cos",       m_recNumD02KPiPiPiPi0,m_D02KPiPiPiPi0_photon2_cos);
00334           status = m_tuple4->addItem ( "recNumD02KsPiPi",       m_recNumD02KsPiPi,0,500);
00335           status = m_tuple4->addIndexedItem ("D02KsPiPi_e",             m_recNumD02KsPiPi,m_D02KsPiPi_e);
00336           status = m_tuple4->addIndexedItem ("D02KsPiPi_p",             m_recNumD02KsPiPi,m_D02KsPiPi_p);
00337           status = m_tuple4->addIndexedItem ("D02KsPiPi_mBC",   m_recNumD02KsPiPi,m_D02KsPiPi_mBC);
00338           status = m_tuple4->addIndexedItem ("D02KsPiPi_mass",  m_recNumD02KsPiPi,m_D02KsPiPi_mass);
00339           status = m_tuple4->addIndexedItem ("D02KsPiPi_deltE", m_recNumD02KsPiPi,m_D02KsPiPi_deltE);
00340           status = m_tuple4->addIndexedItem ("D02KsPiPi_charm", m_recNumD02KsPiPi,m_D02KsPiPi_charm);
00341           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion1_pp", m_recNumD02KsPiPi,m_D02KsPiPi_pion1_pp);
00342           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion2_pp", m_recNumD02KsPiPi,m_D02KsPiPi_pion2_pp);
00343           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion3_pp", m_recNumD02KsPiPi,m_D02KsPiPi_pion3_pp);
00344           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion4_pp", m_recNumD02KsPiPi,m_D02KsPiPi_pion4_pp);
00345           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion1_Id", m_recNumD02KsPiPi,m_D02KsPiPi_pion1_Id);
00346           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion2_Id", m_recNumD02KsPiPi,m_D02KsPiPi_pion2_Id);
00347           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion3_Id", m_recNumD02KsPiPi,m_D02KsPiPi_pion3_Id);
00348           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion4_Id", m_recNumD02KsPiPi,m_D02KsPiPi_pion4_Id);
00349           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion1_vx",      m_recNumD02KsPiPi,      m_D02KsPiPi_pion1_vx);
00350           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion1_vy",      m_recNumD02KsPiPi,      m_D02KsPiPi_pion1_vy);
00351           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion1_vz",      m_recNumD02KsPiPi,      m_D02KsPiPi_pion1_vz);
00352           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion1_vr0",     m_recNumD02KsPiPi,      m_D02KsPiPi_pion1_vr0);
00353           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion1_px",      m_recNumD02KsPiPi,      m_D02KsPiPi_pion1_px);
00354           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion1_py",      m_recNumD02KsPiPi,      m_D02KsPiPi_pion1_py);
00355           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion1_pz",      m_recNumD02KsPiPi,      m_D02KsPiPi_pion1_pz);
00356           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion1_pp",      m_recNumD02KsPiPi,              m_D02KsPiPi_pion1_pp);
00357           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion1_E",       m_recNumD02KsPiPi,      m_D02KsPiPi_pion1_E);
00358           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion1_cos",     m_recNumD02KsPiPi,      m_D02KsPiPi_pion1_cos);
00359           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion2_vx",      m_recNumD02KsPiPi,      m_D02KsPiPi_pion2_vx);
00360           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion2_vy",      m_recNumD02KsPiPi,      m_D02KsPiPi_pion2_vy);
00361           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion2_vz",      m_recNumD02KsPiPi,      m_D02KsPiPi_pion2_vz);
00362           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion2_vr0",     m_recNumD02KsPiPi,      m_D02KsPiPi_pion2_vr0);
00363           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion2_px",      m_recNumD02KsPiPi,      m_D02KsPiPi_pion2_px);
00364           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion2_py",      m_recNumD02KsPiPi,      m_D02KsPiPi_pion2_py);
00365           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion2_pz",      m_recNumD02KsPiPi,      m_D02KsPiPi_pion2_pz);
00366           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion2_pp",      m_recNumD02KsPiPi,              m_D02KsPiPi_pion2_pp);
00367           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion2_E",       m_recNumD02KsPiPi,      m_D02KsPiPi_pion2_E);
00368           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion2_cos",     m_recNumD02KsPiPi,      m_D02KsPiPi_pion2_cos);
00369           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion3_vx",      m_recNumD02KsPiPi,      m_D02KsPiPi_pion3_vx);
00370           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion3_vy",      m_recNumD02KsPiPi,      m_D02KsPiPi_pion3_vy);
00371           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion3_vz",      m_recNumD02KsPiPi,      m_D02KsPiPi_pion3_vz);
00372           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion3_vr0",     m_recNumD02KsPiPi,      m_D02KsPiPi_pion3_vr0);
00373           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion3_px",      m_recNumD02KsPiPi,      m_D02KsPiPi_pion3_px);
00374           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion3_py",      m_recNumD02KsPiPi,      m_D02KsPiPi_pion3_py);
00375           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion3_pz",      m_recNumD02KsPiPi,      m_D02KsPiPi_pion3_pz);
00376           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion3_pp",      m_recNumD02KsPiPi,              m_D02KsPiPi_pion3_pp);
00377           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion3_E",       m_recNumD02KsPiPi,      m_D02KsPiPi_pion3_E);
00378           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion3_cos",     m_recNumD02KsPiPi,      m_D02KsPiPi_pion3_cos);
00379           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion4_vx",      m_recNumD02KsPiPi,      m_D02KsPiPi_pion4_vx);
00380           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion4_vy",      m_recNumD02KsPiPi,      m_D02KsPiPi_pion4_vy);
00381           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion4_vz",      m_recNumD02KsPiPi,      m_D02KsPiPi_pion4_vz);
00382           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion4_vr0",     m_recNumD02KsPiPi,      m_D02KsPiPi_pion4_vr0);
00383           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion4_px",      m_recNumD02KsPiPi,      m_D02KsPiPi_pion4_px);
00384           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion4_py",      m_recNumD02KsPiPi,      m_D02KsPiPi_pion4_py);
00385           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion4_pz",      m_recNumD02KsPiPi,      m_D02KsPiPi_pion4_pz);
00386           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion4_pp",      m_recNumD02KsPiPi,              m_D02KsPiPi_pion4_pp);
00387           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion4_E",       m_recNumD02KsPiPi,      m_D02KsPiPi_pion4_E);
00388           status = m_tuple4->addIndexedItem ("D02KsPiPi_pion4_cos",     m_recNumD02KsPiPi,      m_D02KsPiPi_pion4_cos);
00389           status = m_tuple4->addIndexedItem ("D02KsPiPi_ks_mass",       m_recNumD02KsPiPi,      m_D02KsPiPi_ks_mass);
00390           status = m_tuple4->addIndexedItem ("D02KsPiPi_ks_vfitmass",   m_recNumD02KsPiPi,m_D02KsPiPi_ks_vfitmass);
00391           status = m_tuple4->addIndexedItem ("D02KsPiPi_ks_pvchi2",     m_recNumD02KsPiPi,m_D02KsPiPi_ks_pvchi2);
00392           status = m_tuple4->addIndexedItem ("D02KsPiPi_ks_svchi2",     m_recNumD02KsPiPi,m_D02KsPiPi_ks_svchi2);
00393           status = m_tuple4->addIndexedItem ("D02KsPiPi_ks_dLength",    m_recNumD02KsPiPi,m_D02KsPiPi_ks_dLength);
00394           status = m_tuple4->addIndexedItem ("D02KsPiPi_ks_dLengthErrorRatio",          m_recNumD02KsPiPi,m_D02KsPiPi_ks_dLengthErrorRatio);
00395           status = m_tuple4->addItem ( "recNumD02KsPiPiPi0",            m_recNumD02KsPiPiPi0,0,500);
00396           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_e",          m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_e);
00397           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_p",          m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_p);
00398           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_mBC",        m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_mBC);
00399           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_mass",       m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_mass);
00400           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_deltE",      m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_deltE);
00401           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_charm",      m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_charm);
00402           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_pion1_pp",   m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_pion1_pp);
00403           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_pion2_pp",   m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_pion2_pp);
00404           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_pion3_pp",   m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_pion3_pp);
00405           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_pion4_pp",   m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_pion4_pp);
00406           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_pion1_Id",   m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_pion1_Id);
00407           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_pion2_Id",   m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_pion2_Id);
00408           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_pion3_Id",   m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_pion3_Id);
00409           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_pion4_Id",   m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_pion4_Id);
00410           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_pi0mass",    m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_pi0mass);
00411           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_photon1_Id", m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_photon1_Id);
00412           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_photon2_Id", m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_photon2_Id);
00413           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_photon1_E",  m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_photon1_E);
00414           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_photon2_E",  m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_photon2_E);
00415           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_photon1_cos",        m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_photon1_cos);
00416           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_photon2_cos",        m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_photon2_cos);
00417           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_ks_mass",    m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_ks_mass);
00418           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_ks_vfitmass",                m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_ks_vfitmass);
00419           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_ks_pvchi2",  m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_ks_pvchi2);
00420           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_ks_svchi2",  m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_ks_svchi2);
00421           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_ks_dLength",                 m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_ks_dLength);
00422           status = m_tuple4->addIndexedItem ("D02KsPiPiPi0_ks_dLengthErrorRatio",               m_recNumD02KsPiPiPi0,m_D02KsPiPiPi0_ks_dLengthErrorRatio);
00423           status = m_tuple4->addItem ( "recNumDC2KPiPi",                        m_recNumDC2KPiPi,0,500);
00424           status = m_tuple4->addIndexedItem ("DC2KPiPi_e",              m_recNumDC2KPiPi,m_DC2KPiPi_e);
00425           status = m_tuple4->addIndexedItem ("DC2KPiPi_p",              m_recNumDC2KPiPi,m_DC2KPiPi_p);
00426           status = m_tuple4->addIndexedItem ("DC2KPiPi_mBC",            m_recNumDC2KPiPi,m_DC2KPiPi_mBC);
00427           status = m_tuple4->addIndexedItem ("DC2KPiPi_mass",           m_recNumDC2KPiPi,m_DC2KPiPi_mass);
00428           status = m_tuple4->addIndexedItem ("DC2KPiPi_deltE",          m_recNumDC2KPiPi,m_DC2KPiPi_deltE);
00429           status = m_tuple4->addIndexedItem ("DC2KPiPi_charge",         m_recNumDC2KPiPi,m_DC2KPiPi_charge);
00430           status = m_tuple4->addIndexedItem ("DC2KPiPi_charm",          m_recNumDC2KPiPi,m_DC2KPiPi_charm);
00431           status = m_tuple4->addIndexedItem ("DC2KPiPi_kaon_pp",        m_recNumDC2KPiPi,m_DC2KPiPi_kaon_pp);
00432           status = m_tuple4->addIndexedItem ("DC2KPiPi_pion1_pp",       m_recNumDC2KPiPi,m_DC2KPiPi_pion1_pp);
00433           status = m_tuple4->addIndexedItem ("DC2KPiPi_pion2_pp",       m_recNumDC2KPiPi,m_DC2KPiPi_pion2_pp);
00434           status = m_tuple4->addIndexedItem ("DC2KPiPi_kaon_Id",        m_recNumDC2KPiPi,m_DC2KPiPi_kaon_Id);
00435           status = m_tuple4->addIndexedItem ("DC2KPiPi_pion1_Id",       m_recNumDC2KPiPi,m_DC2KPiPi_pion1_Id);
00436           status = m_tuple4->addIndexedItem ("DC2KPiPi_pion2_Id",       m_recNumDC2KPiPi,m_DC2KPiPi_pion2_Id);
00437           status = m_tuple4->addIndexedItem ("DC2KPiPi_kaon_vx",        m_recNumDC2KPiPi,       m_DC2KPiPi_kaon_vx);
00438           status = m_tuple4->addIndexedItem ("DC2KPiPi_kaon_vy",        m_recNumDC2KPiPi,       m_DC2KPiPi_kaon_vy);
00439           status = m_tuple4->addIndexedItem ("DC2KPiPi_kaon_vz",        m_recNumDC2KPiPi,       m_DC2KPiPi_kaon_vz);
00440           status = m_tuple4->addIndexedItem ("DC2KPiPi_kaon_vr0",       m_recNumDC2KPiPi,       m_DC2KPiPi_kaon_vr0);
00441           status = m_tuple4->addIndexedItem ("DC2KPiPi_kaon_px",        m_recNumDC2KPiPi,       m_DC2KPiPi_kaon_px);
00442           status = m_tuple4->addIndexedItem ("DC2KPiPi_kaon_py",        m_recNumDC2KPiPi,       m_DC2KPiPi_kaon_py);
00443           status = m_tuple4->addIndexedItem ("DC2KPiPi_kaon_pz",        m_recNumDC2KPiPi,       m_DC2KPiPi_kaon_pz);
00444           status = m_tuple4->addIndexedItem ("DC2KPiPi_kaon_E",         m_recNumDC2KPiPi,       m_DC2KPiPi_kaon_E);
00445           status = m_tuple4->addIndexedItem ("DC2KPiPi_kaon_cos",       m_recNumDC2KPiPi,       m_DC2KPiPi_kaon_cos);
00446           status = m_tuple4->addIndexedItem ("DC2KPiPi_pion1_vx",       m_recNumDC2KPiPi,       m_DC2KPiPi_pion1_vx);
00447           status = m_tuple4->addIndexedItem ("DC2KPiPi_pion1_vy",       m_recNumDC2KPiPi,       m_DC2KPiPi_pion1_vy);
00448           status = m_tuple4->addIndexedItem ("DC2KPiPi_pion1_vz",       m_recNumDC2KPiPi,       m_DC2KPiPi_pion1_vz);
00449           status = m_tuple4->addIndexedItem ("DC2KPiPi_pion1_vr0",      m_recNumDC2KPiPi,       m_DC2KPiPi_pion1_vr0);
00450           status = m_tuple4->addIndexedItem ("DC2KPiPi_pion1_px",       m_recNumDC2KPiPi,       m_DC2KPiPi_pion1_px);
00451           status = m_tuple4->addIndexedItem ("DC2KPiPi_pion1_py",       m_recNumDC2KPiPi,       m_DC2KPiPi_pion1_py);
00452           status = m_tuple4->addIndexedItem ("DC2KPiPi_pion1_pz",       m_recNumDC2KPiPi,       m_DC2KPiPi_pion1_pz);
00453           status = m_tuple4->addIndexedItem ("DC2KPiPi_pion1_E",        m_recNumDC2KPiPi,       m_DC2KPiPi_pion1_E);
00454           status = m_tuple4->addIndexedItem ("DC2KPiPi_pion1_cos",      m_recNumDC2KPiPi,       m_DC2KPiPi_pion1_cos);
00455           status = m_tuple4->addIndexedItem ("DC2KPiPi_pion2_vx",       m_recNumDC2KPiPi,       m_DC2KPiPi_pion2_vx);
00456           status = m_tuple4->addIndexedItem ("DC2KPiPi_pion2_vy",       m_recNumDC2KPiPi,       m_DC2KPiPi_pion2_vy);
00457           status = m_tuple4->addIndexedItem ("DC2KPiPi_pion2_vz",       m_recNumDC2KPiPi,       m_DC2KPiPi_pion2_vz);
00458           status = m_tuple4->addIndexedItem ("DC2KPiPi_pion2_vr0",      m_recNumDC2KPiPi,       m_DC2KPiPi_pion2_vr0);
00459           status = m_tuple4->addIndexedItem ("DC2KPiPi_pion2_px",       m_recNumDC2KPiPi,       m_DC2KPiPi_pion2_px);
00460           status = m_tuple4->addIndexedItem ("DC2KPiPi_pion2_py",       m_recNumDC2KPiPi,       m_DC2KPiPi_pion2_py);
00461           status = m_tuple4->addIndexedItem ("DC2KPiPi_pion2_pz",       m_recNumDC2KPiPi,       m_DC2KPiPi_pion2_pz);
00462           status = m_tuple4->addIndexedItem ("DC2KPiPi_pion2_E",        m_recNumDC2KPiPi,       m_DC2KPiPi_pion2_E);
00463           status = m_tuple4->addIndexedItem ("DC2KPiPi_pion2_cos",      m_recNumDC2KPiPi,       m_DC2KPiPi_pion2_cos);
00464           status = m_tuple4->addItem ( "recNumDC2KPiPiPi0",                     m_recNumDC2KPiPiPi0,0,500);
00465           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_e",           m_recNumDC2KPiPiPi0,m_DC2KPiPiPi0_e);
00466           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_p",           m_recNumDC2KPiPiPi0,m_DC2KPiPiPi0_p);
00467           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_mBC",         m_recNumDC2KPiPiPi0,m_DC2KPiPiPi0_mBC);
00468           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_mass",        m_recNumDC2KPiPiPi0,m_DC2KPiPiPi0_mass);
00469           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_deltE",       m_recNumDC2KPiPiPi0,m_DC2KPiPiPi0_deltE);
00470           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_charge",      m_recNumDC2KPiPiPi0,m_DC2KPiPiPi0_charge);
00471           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_charm",       m_recNumDC2KPiPiPi0,m_DC2KPiPiPi0_charm);
00472           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_kaon_pp",     m_recNumDC2KPiPiPi0,m_DC2KPiPiPi0_kaon_pp);
00473           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_pion1_pp",    m_recNumDC2KPiPiPi0,m_DC2KPiPiPi0_pion1_pp);
00474           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_pion2_pp",    m_recNumDC2KPiPiPi0,m_DC2KPiPiPi0_pion2_pp);
00475           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_kaon_Id",     m_recNumDC2KPiPiPi0,m_DC2KPiPiPi0_kaon_Id);
00476           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_pion1_Id",    m_recNumDC2KPiPiPi0,m_DC2KPiPiPi0_pion1_Id);
00477           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_pion2_Id",    m_recNumDC2KPiPiPi0,m_DC2KPiPiPi0_pion2_Id);
00478           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_kaon_vx",     m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_kaon_vx);
00479           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_kaon_vy",     m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_kaon_vy);
00480           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_kaon_vz",     m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_kaon_vz);
00481           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_kaon_vr0",    m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_kaon_vr0);
00482           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_kaon_px",     m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_kaon_px);
00483           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_kaon_py",     m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_kaon_py);
00484           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_kaon_pz",     m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_kaon_pz);
00485           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_kaon_E",      m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_kaon_E);
00486           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_kaon_cos",    m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_kaon_cos);
00487           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_pion1_vx",    m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_pion1_vx);
00488           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_pion1_vy",    m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_pion1_vy);
00489           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_pion1_vz",    m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_pion1_vz);
00490           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_pion1_vr0",   m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_pion1_vr0);
00491           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_pion1_px",    m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_pion1_px);
00492           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_pion1_py",    m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_pion1_py);
00493           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_pion1_pz",    m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_pion1_pz);
00494           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_pion1_E",     m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_pion1_E);
00495           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_pion1_cos",   m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_pion1_cos);
00496           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_pion2_vx",    m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_pion2_vx);
00497           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_pion2_vy",    m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_pion2_vy);
00498           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_pion2_vz",    m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_pion2_vz);
00499           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_pion2_vr0",   m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_pion2_vr0);
00500           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_pion2_px",    m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_pion2_px);
00501           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_pion2_py",    m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_pion2_py);
00502           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_pion2_pz",    m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_pion2_pz);
00503           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_pion2_E",     m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_pion2_E);
00504           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_pion2_cos",   m_recNumDC2KPiPiPi0,    m_DC2KPiPiPi0_pion2_cos);
00505           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_pi0mass",     m_recNumDC2KPiPiPi0,m_DC2KPiPiPi0_pi0mass);
00506           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_photon1_Id",  m_recNumDC2KPiPiPi0,m_DC2KPiPiPi0_photon1_Id);
00507           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_photon2_Id",  m_recNumDC2KPiPiPi0,m_DC2KPiPiPi0_photon2_Id);
00508           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_photon1_E",   m_recNumDC2KPiPiPi0,m_DC2KPiPiPi0_photon1_E);
00509           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_photon2_E",   m_recNumDC2KPiPiPi0,m_DC2KPiPiPi0_photon2_E);
00510           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_photon1_cos", m_recNumDC2KPiPiPi0,m_DC2KPiPiPi0_photon1_cos);
00511           status = m_tuple4->addIndexedItem ("DC2KPiPiPi0_photon2_cos", m_recNumDC2KPiPiPi0,m_DC2KPiPiPi0_photon2_cos);
00512           status = m_tuple4->addItem ( "recNumDC2KsPiPi0",              m_recNumDC2KsPiPi0,0,500);
00513           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_e",            m_recNumDC2KsPiPi0,m_DC2KsPiPi0_e);
00514           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_p",            m_recNumDC2KsPiPi0,m_DC2KsPiPi0_p);
00515           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_mBC",          m_recNumDC2KsPiPi0,m_DC2KsPiPi0_mBC);
00516           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_mass",         m_recNumDC2KsPiPi0,m_DC2KsPiPi0_mass);
00517           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_deltE",        m_recNumDC2KsPiPi0,m_DC2KsPiPi0_deltE);
00518           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_charge",       m_recNumDC2KsPiPi0,m_DC2KsPiPi0_charge);
00519           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_charm",        m_recNumDC2KsPiPi0,m_DC2KsPiPi0_charm);
00520           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion1_pp",     m_recNumDC2KsPiPi0,m_DC2KsPiPi0_pion1_pp);
00521           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion2_pp",     m_recNumDC2KsPiPi0,m_DC2KsPiPi0_pion2_pp);
00522           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion3_pp",     m_recNumDC2KsPiPi0,m_DC2KsPiPi0_pion3_pp);
00523           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion1_Id",     m_recNumDC2KsPiPi0,m_DC2KsPiPi0_pion1_Id);
00524           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion2_Id",     m_recNumDC2KsPiPi0,m_DC2KsPiPi0_pion2_Id);
00525           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion3_Id",     m_recNumDC2KsPiPi0,m_DC2KsPiPi0_pion3_Id);
00526           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion1_vx",     m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion1_vx);
00527           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion1_vy",     m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion1_vy);
00528           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion1_vz",     m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion1_vz);
00529           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion1_vr0",    m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion1_vr0);
00530           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion1_px",     m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion1_px);
00531           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion1_py",     m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion1_py);
00532           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion1_pz",     m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion1_pz);
00533           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion1_E",      m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion1_E);
00534           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion1_cos",    m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion1_cos);
00535           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion2_vx",     m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion2_vx);
00536           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion2_vy",     m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion2_vy);
00537           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion2_vz",     m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion2_vz);
00538           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion2_vr0",    m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion2_vr0);
00539           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion2_px",     m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion2_px);
00540           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion2_py",     m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion2_py);
00541           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion2_pz",     m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion2_pz);
00542           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion2_E",      m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion2_E);
00543           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion2_cos",    m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion2_cos);
00544           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion3_vx",     m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion3_vx);
00545           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion3_vy",     m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion3_vy);
00546           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion3_vz",     m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion3_vz);
00547           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion3_vr0",    m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion3_vr0);
00548           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion3_px",     m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion3_px);
00549           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion3_py",     m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion3_py);
00550           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion3_pz",     m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion3_pz);
00551           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion3_E",      m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion3_E);
00552           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pion3_cos",    m_recNumDC2KsPiPi0,     m_DC2KsPiPi0_pion3_cos);
00553           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_pi0mass",      m_recNumDC2KsPiPi0,m_DC2KsPiPi0_pi0mass);
00554           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_photon1_Id",   m_recNumDC2KsPiPi0,m_DC2KsPiPi0_photon1_Id);
00555           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_photon2_Id",   m_recNumDC2KsPiPi0,m_DC2KsPiPi0_photon2_Id);
00556           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_photon1_E",    m_recNumDC2KsPiPi0,m_DC2KsPiPi0_photon1_E);
00557           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_photon2_E",    m_recNumDC2KsPiPi0,m_DC2KsPiPi0_photon2_E);
00558           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_photon1_cos",  m_recNumDC2KsPiPi0,m_DC2KsPiPi0_photon1_cos);
00559           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_photon2_cos",  m_recNumDC2KsPiPi0,m_DC2KsPiPi0_photon2_cos);
00560           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_ks_mass",      m_recNumDC2KsPiPi0,m_DC2KsPiPi0_ks_mass);
00561           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_ks_vfitmass",  m_recNumDC2KsPiPi0,m_DC2KsPiPi0_ks_vfitmass);
00562           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_ks_pvchi2",    m_recNumDC2KsPiPi0,m_DC2KsPiPi0_ks_pvchi2);
00563           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_ks_svchi2",    m_recNumDC2KsPiPi0,m_DC2KsPiPi0_ks_svchi2);
00564           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_ks_dLength",   m_recNumDC2KsPiPi0,m_DC2KsPiPi0_ks_dLength);
00565           status = m_tuple4->addIndexedItem ("DC2KsPiPi0_ks_dLengthErrorRatio",         m_recNumDC2KsPiPi0,m_DC2KsPiPi0_ks_dLengthErrorRatio);
00566           status = m_tuple4->addItem ( "recNumDC2KsPiPiPi",                     m_recNumDC2KsPiPiPi,0,500);
00567           status = m_tuple4->addIndexedItem ("DC2KsPiPiPi_e",           m_recNumDC2KsPiPiPi,m_DC2KsPiPiPi_e);
00568           status = m_tuple4->addIndexedItem ("DC2KsPiPiPi_p",           m_recNumDC2KsPiPiPi,m_DC2KsPiPiPi_p);
00569           status = m_tuple4->addIndexedItem ("DC2KsPiPiPi_mBC",         m_recNumDC2KsPiPiPi,m_DC2KsPiPiPi_mBC);
00570           status = m_tuple4->addIndexedItem ("DC2KsPiPiPi_mass",        m_recNumDC2KsPiPiPi,m_DC2KsPiPiPi_mass);
00571           status = m_tuple4->addIndexedItem ("DC2KsPiPiPi_deltE",       m_recNumDC2KsPiPiPi,m_DC2KsPiPiPi_deltE);
00572           status = m_tuple4->addIndexedItem ("DC2KsPiPiPi_charge",      m_recNumDC2KsPiPiPi,m_DC2KsPiPiPi_charge);
00573           status = m_tuple4->addIndexedItem ("DC2KsPiPiPi_charm",       m_recNumDC2KsPiPiPi,m_DC2KsPiPiPi_charm);
00574           status = m_tuple4->addIndexedItem ("DC2KsPiPiPi_pion1_pp",    m_recNumDC2KsPiPiPi,m_DC2KsPiPiPi_pion1_pp);
00575           status = m_tuple4->addIndexedItem ("DC2KsPiPiPi_pion2_pp",    m_recNumDC2KsPiPiPi,m_DC2KsPiPiPi_pion2_pp);
00576           status = m_tuple4->addIndexedItem ("DC2KsPiPiPi_pion3_pp",    m_recNumDC2KsPiPiPi,m_DC2KsPiPiPi_pion3_pp);
00577           status = m_tuple4->addIndexedItem ("DC2KsPiPiPi_pion4_pp",    m_recNumDC2KsPiPiPi,m_DC2KsPiPiPi_pion4_pp);
00578           status = m_tuple4->addIndexedItem ("DC2KsPiPiPi_pion5_pp",    m_recNumDC2KsPiPiPi,m_DC2KsPiPiPi_pion5_pp);
00579           status = m_tuple4->addIndexedItem ("DC2KsPiPiPi_pion1_Id",    m_recNumDC2KsPiPiPi,m_DC2KsPiPiPi_pion1_Id);
00580           status = m_tuple4->addIndexedItem ("DC2KsPiPiPi_pion2_Id",    m_recNumDC2KsPiPiPi,m_DC2KsPiPiPi_pion2_Id);
00581           status = m_tuple4->addIndexedItem ("DC2KsPiPiPi_pion3_Id",    m_recNumDC2KsPiPiPi,m_DC2KsPiPiPi_pion3_Id);
00582           status = m_tuple4->addIndexedItem ("DC2KsPiPiPi_pion4_Id",    m_recNumDC2KsPiPiPi,m_DC2KsPiPiPi_pion4_Id);
00583           status = m_tuple4->addIndexedItem ("DC2KsPiPiPi_pion5_Id",    m_recNumDC2KsPiPiPi,m_DC2KsPiPiPi_pion5_Id);
00584           status = m_tuple4->addIndexedItem ("DC2KsPiPiPi_ks_mass",     m_recNumDC2KsPiPiPi,m_DC2KsPiPiPi_ks_mass);
00585           status = m_tuple4->addIndexedItem ("DC2KsPiPiPi_ks_vfitmass",         m_recNumDC2KsPiPiPi,m_DC2KsPiPiPi_ks_vfitmass);
00586           status = m_tuple4->addIndexedItem ("DC2KsPiPiPi_ks_pvchi2",m_recNumDC2KsPiPiPi,m_DC2KsPiPiPi_ks_pvchi2);
00587           status = m_tuple4->addIndexedItem ("DC2KsPiPiPi_ks_svchi2",m_recNumDC2KsPiPiPi,m_DC2KsPiPiPi_ks_svchi2);
00588           status = m_tuple4->addIndexedItem ("DC2KsPiPiPi_ks_dLength",  m_recNumDC2KsPiPiPi,m_DC2KsPiPiPi_ks_dLength);
00589           status = m_tuple4->addIndexedItem ("DC2KsPiPiPi_ks_dLengthErrorRatio",        m_recNumDC2KsPiPiPi,m_DC2KsPiPiPi_ks_dLengthErrorRatio);
00590           status = m_tuple4->addItem ( "recNumDC2KKPi",                         m_recNumDC2KKPi,0,500);
00591           status = m_tuple4->addIndexedItem ("DC2KKPi_e",               m_recNumDC2KKPi,m_DC2KKPi_e);
00592           status = m_tuple4->addIndexedItem ("DC2KKPi_p",               m_recNumDC2KKPi,m_DC2KKPi_p);
00593           status = m_tuple4->addIndexedItem ("DC2KKPi_mBC",             m_recNumDC2KKPi,m_DC2KKPi_mBC);
00594           status = m_tuple4->addIndexedItem ("DC2KKPi_mass",            m_recNumDC2KKPi,m_DC2KKPi_mass);
00595           status = m_tuple4->addIndexedItem ("DC2KKPi_deltE",           m_recNumDC2KKPi,m_DC2KKPi_deltE);
00596           status = m_tuple4->addIndexedItem ("DC2KKPi_charge",          m_recNumDC2KKPi,m_DC2KKPi_charge);
00597           status = m_tuple4->addIndexedItem ("DC2KKPi_charm",           m_recNumDC2KKPi,m_DC2KKPi_charm);
00598           status = m_tuple4->addIndexedItem ("DC2KKPi_kaon1_pp",        m_recNumDC2KKPi,m_DC2KKPi_kaon1_pp);
00599           status = m_tuple4->addIndexedItem ("DC2KKPi_kaon2_pp",        m_recNumDC2KKPi,m_DC2KKPi_kaon2_pp);
00600           status = m_tuple4->addIndexedItem ("DC2KKPi_pion_pp",         m_recNumDC2KKPi,m_DC2KKPi_pion_pp);
00601           status = m_tuple4->addIndexedItem ("DC2KKPi_kaon1_Id",        m_recNumDC2KKPi,m_DC2KKPi_kaon1_Id);
00602           status = m_tuple4->addIndexedItem ("DC2KKPi_kaon2_Id",        m_recNumDC2KKPi,m_DC2KKPi_kaon2_Id);
00603           status = m_tuple4->addIndexedItem ("DC2KKPi_pion_Id",         m_recNumDC2KKPi,m_DC2KKPi_pion_Id);
00604           status = m_tuple4->addIndexedItem ("D02KPi_kaon_Id_Dedx",                     m_recNumD02KPi, m_D02KPi_kaon_Id_Dedx); 
00605           status = m_tuple4->addIndexedItem ("D02KPi_pion_Id_Dedx",                     m_recNumD02KPi, m_D02KPi_pion_Id_Dedx); 
00606           status = m_tuple4->addIndexedItem ("D02KPi_kaon_chipi",                       m_recNumD02KPi, m_D02KPi_kaon_chipi); 
00607           status = m_tuple4->addIndexedItem ("D02KPi_kaon_chik",                        m_recNumD02KPi, m_D02KPi_kaon_chik); 
00608           status = m_tuple4->addIndexedItem ("D02KPi_kaon_ghit",                        m_recNumD02KPi, m_D02KPi_kaon_ghit);
00609           status = m_tuple4->addIndexedItem ("D02KPi_kaon_thit",                        m_recNumD02KPi, m_D02KPi_kaon_thit);
00610           status = m_tuple4->addIndexedItem ("D02KPi_kaon_probPH",                      m_recNumD02KPi, m_D02KPi_kaon_probPH);
00611           status = m_tuple4->addIndexedItem ("D02KPi_kaon_normPH",                      m_recNumD02KPi, m_D02KPi_kaon_normPH);
00612           status = m_tuple4->addIndexedItem ("D02KPi_pion_chipi",                       m_recNumD02KPi, m_D02KPi_pion_chipi);
00613           status = m_tuple4->addIndexedItem ("D02KPi_pion_chik",                        m_recNumD02KPi, m_D02KPi_pion_chik);
00614           status = m_tuple4->addIndexedItem ("D02KPi_pion_ghit",                        m_recNumD02KPi, m_D02KPi_pion_ghit);
00615           status = m_tuple4->addIndexedItem ("D02KPi_pion_thit",                        m_recNumD02KPi, m_D02KPi_pion_thit);
00616           status = m_tuple4->addIndexedItem ("D02KPi_pion_probPH",                      m_recNumD02KPi, m_D02KPi_pion_probPH);
00617           status = m_tuple4->addIndexedItem ("D02KPi_pion_normPH",                      m_recNumD02KPi, m_D02KPi_pion_normPH);
00618           status = m_tuple4->addIndexedItem ("D02KPiPi0_kaon_Id_Dedx",                  m_recNumD02KPiPi0,      m_D02KPiPi0_kaon_Id_Dedx); 
00619           status = m_tuple4->addIndexedItem ("D02KPiPi0_pion_Id_Dedx",                  m_recNumD02KPiPi0,      m_D02KPiPi0_pion_Id_Dedx); 
00620           status = m_tuple4->addIndexedItem ("D02KPiPi0_kaon_chipi",                    m_recNumD02KPiPi0,      m_D02KPiPi0_kaon_chipi); 
00621           status = m_tuple4->addIndexedItem ("D02KPiPi0_kaon_chik",                     m_recNumD02KPiPi0,      m_D02KPiPi0_kaon_chik); 
00622           status = m_tuple4->addIndexedItem ("D02KPiPi0_kaon_ghit",                     m_recNumD02KPiPi0,      m_D02KPiPi0_kaon_ghit);
00623           status = m_tuple4->addIndexedItem ("D02KPiPi0_kaon_thit",                     m_recNumD02KPiPi0,      m_D02KPiPi0_kaon_thit);
00624           status = m_tuple4->addIndexedItem ("D02KPiPi0_kaon_probPH",                   m_recNumD02KPiPi0,      m_D02KPiPi0_kaon_probPH);
00625           status = m_tuple4->addIndexedItem ("D02KPiPi0_kaon_normPH",                   m_recNumD02KPiPi0,      m_D02KPiPi0_kaon_normPH);
00626           status = m_tuple4->addIndexedItem ("D02KPiPi0_pion_chipi",                    m_recNumD02KPiPi0,      m_D02KPiPi0_pion_chipi);
00627           status = m_tuple4->addIndexedItem ("D02KPiPi0_pion_chik",                     m_recNumD02KPiPi0,      m_D02KPiPi0_pion_chik);
00628           status = m_tuple4->addIndexedItem ("D02KPiPi0_pion_ghit",                     m_recNumD02KPiPi0,      m_D02KPiPi0_pion_ghit);
00629           status = m_tuple4->addIndexedItem ("D02KPiPi0_pion_thit",                     m_recNumD02KPiPi0,      m_D02KPiPi0_pion_thit);
00630           status = m_tuple4->addIndexedItem ("D02KPiPi0_pion_probPH",                   m_recNumD02KPiPi0,      m_D02KPiPi0_pion_probPH);
00631           status = m_tuple4->addIndexedItem ("D02KPiPi0_pion_normPH",                   m_recNumD02KPiPi0,      m_D02KPiPi0_pion_normPH);
00632           status = m_tuple4->addIndexedItem ("D02KPi_kaon_Id_tof",              m_recNumD02KPi, m_D02KPi_kaon_Id_tof); 
00633           status = m_tuple4->addIndexedItem ("D02KPi_pion_Id_tof",              m_recNumD02KPi, m_D02KPi_pion_Id_tof); 
00634           status = m_tuple4->addIndexedItem ("D02KPi_kaon_ec_toff_mu",          m_recNumD02KPi, m_D02KPi_kaon_ec_toff_mu); 
00635           status = m_tuple4->addIndexedItem ("D02KPi_kaon_ec_toff_pi",          m_recNumD02KPi, m_D02KPi_kaon_ec_toff_pi); 
00636           status = m_tuple4->addIndexedItem ("D02KPi_kaon_ec_toff_k",           m_recNumD02KPi, m_D02KPi_kaon_ec_toff_k); 
00637           status = m_tuple4->addIndexedItem ("D02KPi_kaon_br_toff_mu",          m_recNumD02KPi, 7,m_D02KPi_kaon_br_toff_mu);
00638           status = m_tuple4->addIndexedItem ("D02KPi_kaon_br_toff_pi",          m_recNumD02KPi, 7,m_D02KPi_kaon_br_toff_pi);
00639           status = m_tuple4->addIndexedItem ("D02KPi_kaon_br_toff_k",           m_recNumD02KPi, 7,m_D02KPi_kaon_br_toff_k);
00640           status = m_tuple4->addIndexedItem ("D02KPi_pion_ec_toff_pi",          m_recNumD02KPi, m_D02KPi_pion_ec_toff_mu); 
00641           status = m_tuple4->addIndexedItem ("D02KPi_pion_ec_toff_mu",          m_recNumD02KPi, m_D02KPi_pion_ec_toff_pi); 
00642           status = m_tuple4->addIndexedItem ("D02KPi_pion_ec_toff_k",           m_recNumD02KPi, m_D02KPi_pion_ec_toff_k); 
00643           status = m_tuple4->addIndexedItem ("D02KPi_pion_br_toff_mu",          m_recNumD02KPi, 7,m_D02KPi_pion_br_toff_mu);
00644           status = m_tuple4->addIndexedItem ("D02KPi_pion_br_toff_pi",          m_recNumD02KPi, 7,m_D02KPi_pion_br_toff_pi);
00645           status = m_tuple4->addIndexedItem ("D02KPi_pion_br_toff_k",           m_recNumD02KPi, 7,m_D02KPi_pion_br_toff_k);
00646           status = m_tuple4->addIndexedItem ("D02KPiPi0_kaon_Id_tof",           m_recNumD02KPiPi0,      m_D02KPiPi0_kaon_Id_tof); 
00647           status = m_tuple4->addIndexedItem ("D02KPiPi0_pion_Id_tof",           m_recNumD02KPiPi0,      m_D02KPiPi0_pion_Id_tof); 
00648           status = m_tuple4->addIndexedItem ("D02KPiPi0_kaon_ec_toff_mu",               m_recNumD02KPiPi0,      m_D02KPiPi0_kaon_ec_toff_mu); 
00649           status = m_tuple4->addIndexedItem ("D02KPiPi0_kaon_ec_toff_pi",               m_recNumD02KPiPi0,      m_D02KPiPi0_kaon_ec_toff_pi); 
00650           status = m_tuple4->addIndexedItem ("D02KPiPi0_kaon_ec_toff_k",                m_recNumD02KPiPi0,      m_D02KPiPi0_kaon_ec_toff_k); 
00651           status = m_tuple4->addIndexedItem ("D02KPiPi0_kaon_br_toff_mu",               m_recNumD02KPiPi0,      7,m_D02KPiPi0_kaon_br_toff_mu);
00652           status = m_tuple4->addIndexedItem ("D02KPiPi0_kaon_br_toff_pi",               m_recNumD02KPiPi0,      7,m_D02KPiPi0_kaon_br_toff_pi);
00653           status = m_tuple4->addIndexedItem ("D02KPiPi0_kaon_br_toff_k",                m_recNumD02KPiPi0,      7,m_D02KPiPi0_kaon_br_toff_k);
00654           status = m_tuple4->addIndexedItem ("D02KPiPi0_pion_ec_toff_mu",               m_recNumD02KPiPi0,      m_D02KPiPi0_pion_ec_toff_mu); 
00655           status = m_tuple4->addIndexedItem ("D02KPiPi0_pion_ec_toff_pi",               m_recNumD02KPiPi0,      m_D02KPiPi0_pion_ec_toff_pi); 
00656           status = m_tuple4->addIndexedItem ("D02KPiPi0_pion_ec_toff_k",                m_recNumD02KPiPi0,      m_D02KPiPi0_pion_ec_toff_k); 
00657           status = m_tuple4->addIndexedItem ("D02KPiPi0_pion_br_toff_mu",               m_recNumD02KPiPi0,      7,m_D02KPiPi0_pion_br_toff_mu);
00658           status = m_tuple4->addIndexedItem ("D02KPiPi0_pion_br_toff_pi",               m_recNumD02KPiPi0,      7,m_D02KPiPi0_pion_br_toff_pi);
00659           status = m_tuple4->addIndexedItem ("D02KPiPi0_pion_br_toff_k",                m_recNumD02KPiPi0,      7,m_D02KPiPi0_pion_br_toff_k);
00660       status = m_tuple4->addIndexedItem ("D02KPiPi0_photon1_Id_Emc",            m_recNumD02KPiPi0,      m_D02KPiPi0_photon1_Id_Emc);
00661       status = m_tuple4->addIndexedItem ("D02KPiPi0_photon1_numHits",           m_recNumD02KPiPi0,      m_D02KPiPi0_photon1_numHits);
00662       status = m_tuple4->addIndexedItem ("D02KPiPi0_photon1_secondmoment",      m_recNumD02KPiPi0,      m_D02KPiPi0_photon1_secondmoment);
00663       status = m_tuple4->addIndexedItem ("D02KPiPi0_photon1_x",                 m_recNumD02KPiPi0,      m_D02KPiPi0_photon1_x);
00664       status = m_tuple4->addIndexedItem ("D02KPiPi0_photon1_y",                 m_recNumD02KPiPi0,      m_D02KPiPi0_photon1_y);
00665       status = m_tuple4->addIndexedItem ("D02KPiPi0_photon1_z",                 m_recNumD02KPiPi0,      m_D02KPiPi0_photon1_z);
00666       status = m_tuple4->addIndexedItem ("D02KPiPi0_photon2_Id_Emc",            m_recNumD02KPiPi0,      m_D02KPiPi0_photon2_Id_Emc);
00667       status = m_tuple4->addIndexedItem ("D02KPiPi0_photon2_numHits",           m_recNumD02KPiPi0,      m_D02KPiPi0_photon2_numHits);
00668       status = m_tuple4->addIndexedItem ("D02KPiPi0_photon2_secondmoment",      m_recNumD02KPiPi0,      m_D02KPiPi0_photon2_secondmoment);
00669       status = m_tuple4->addIndexedItem ("D02KPiPi0_photon2_x",                 m_recNumD02KPiPi0,      m_D02KPiPi0_photon2_x);
00670       status = m_tuple4->addIndexedItem ("D02KPiPi0_photon2_y",                 m_recNumD02KPiPi0,      m_D02KPiPi0_photon2_y);
00671       status = m_tuple4->addIndexedItem ("D02KPiPi0_photon2_z",                 m_recNumD02KPiPi0,      m_D02KPiPi0_photon2_z);
00672         }
00673         else    { 
00674           log << MSG::ERROR << "    Cannot book N-tuple:" << long(m_tuple4) << endmsg;
00675           return StatusCode::FAILURE;
00676         }
00677   }
00678 
00679   NTuplePtr nt5(ntupleSvc(), "DQAFILE/Events");
00680   if ( nt5 ) m_tuple5 = nt5;
00681   else {
00682         m_tuple5 = ntupleSvc()->book ("DQAFILE/Events", CLID_ColumnWiseTuple, "ks N-Tuple example");
00683         if ( m_tuple5 )    {
00684           status = m_tuple5->addItem ("evtrun",    m_evtrun);
00685           status = m_tuple5->addItem ("evtrec",    m_evtrec); 
00686           status = m_tuple5->addItem ("evtnch",    m_evtnch);
00687           status = m_tuple5->addItem ("evtnneu",   m_evtnneu);
00688           status = m_tuple5->addItem ("evtbeamE",   m_evtbeamE);
00689         }
00690         else    { 
00691           log << MSG::ERROR << "    Cannot book N-tuple:" << long(m_tuple5) << endmsg;
00692           return StatusCode::FAILURE;
00693         }
00694   }
00695   NTuplePtr nt11(ntupleSvc(), "DQAFILE/CheckDedxAll");
00696   if ( nt11 ) m_tuple11 = nt11;
00697   else {
00698         m_tuple11 = ntupleSvc()->book ("DQAFILE/CheckDedxAll", CLID_ColumnWiseTuple, "ks N-Tuple example");
00699         if ( m_tuple11 )    {
00700           status = m_tuple11->addItem ("ngch",  m_ngch, 0, 100);
00701           status = m_tuple11->addIndexedItem ("trkiddedx",  m_ngch, m_trkiddedx);
00702           status = m_tuple11->addIndexedItem ("mptrk",  m_ngch, m_ptrk);
00703           status = m_tuple11->addIndexedItem ("chie",   m_ngch, m_chie);
00704           status = m_tuple11->addIndexedItem ("chimu",  m_ngch,m_chimu);
00705           status = m_tuple11->addIndexedItem ("chipi",  m_ngch,m_chipi);
00706           status = m_tuple11->addIndexedItem ("chik",   m_ngch,m_chik);
00707           status = m_tuple11->addIndexedItem ("chip",   m_ngch,m_chip);
00708           status = m_tuple11->addIndexedItem ("probPH", m_ngch,m_probPH);
00709           status = m_tuple11->addIndexedItem ("normPH", m_ngch,m_normPH);
00710           status = m_tuple11->addIndexedItem ("ghit",   m_ngch,m_ghit);
00711           status = m_tuple11->addIndexedItem ("thit",   m_ngch,m_thit);
00712         }
00713         else    { 
00714           log << MSG::ERROR << "    Cannot book N-tuple:" << long(m_tuple11) << endmsg;
00715           return StatusCode::FAILURE;
00716         }
00717   }
00718 
00719   NTuplePtr nt12(ntupleSvc(), "DQAFILE/CheckTofChAll");
00720   if ( nt12 ) m_tuple12 = nt12;
00721   else {
00722         m_tuple12 = ntupleSvc()->book ("DQAFILE/CheckTofChAll", CLID_ColumnWiseTuple, "ks N-Tuple example");
00723         if ( m_tuple12 )    {
00724           status = m_tuple12->addItem ("ngch2",  m_ngch2, 0, 100);
00725           status = m_tuple12->addIndexedItem ("trkidtof",  m_ngch2,m_trkidtof);
00726           status = m_tuple12->addIndexedItem ("ptot_etof", m_ngch2,m_ptot_etof);
00727           status = m_tuple12->addIndexedItem ("cntr_etof", m_ngch2,m_cntr_etof);
00728           status = m_tuple12->addIndexedItem ("te_etof",   m_ngch2,m_te_etof);
00729           status = m_tuple12->addIndexedItem ("tmu_etof",  m_ngch2,m_tmu_etof);
00730           status = m_tuple12->addIndexedItem ("tpi_etof",  m_ngch2,m_tpi_etof);
00731           status = m_tuple12->addIndexedItem ("tk_etof",   m_ngch2,m_tk_etof);
00732           status = m_tuple12->addIndexedItem ("tp_etof",   m_ngch2,m_tp_etof);
00733           status = m_tuple12->addIndexedItem ("ph_etof",   m_ngch2,m_ph_etof);
00734           status = m_tuple12->addIndexedItem ("rhit_etof", m_ngch2,m_rhit_etof);
00735           status = m_tuple12->addIndexedItem ("qual_etof", m_ngch2,m_qual_etof);
00736           status = m_tuple12->addIndexedItem ("ec_toff_e", m_ngch2,m_ec_toff_e);
00737           status = m_tuple12->addIndexedItem ("ec_toff_mu",m_ngch2,m_ec_toff_mu);
00738           status = m_tuple12->addIndexedItem ("ec_toff_pi",m_ngch2,m_ec_toff_pi);
00739           status = m_tuple12->addIndexedItem ("ec_toff_k", m_ngch2,m_ec_toff_k);
00740           status = m_tuple12->addIndexedItem ("ec_toff_p", m_ngch2,m_ec_toff_p);
00741           status = m_tuple12->addIndexedItem ("ec_tsig_e", m_ngch2,m_ec_tsig_e);
00742           status = m_tuple12->addIndexedItem ("ec_tsig_mu",m_ngch2,m_ec_tsig_mu);
00743           status = m_tuple12->addIndexedItem ("ec_tsig_pi",m_ngch2,m_ec_tsig_pi);
00744           status = m_tuple12->addIndexedItem ("ec_tsig_k", m_ngch2,m_ec_tsig_k);
00745           status = m_tuple12->addIndexedItem ("ec_tsig_p", m_ngch2,m_ec_tsig_p);
00746           status = m_tuple12->addIndexedItem ("ec_tof",    m_ngch2,m_ec_tof);
00747           status = m_tuple12->addIndexedItem ("ptot_btof1",m_ngch2,m_ptot_btof1);
00748           status = m_tuple12->addIndexedItem ("cntr_btof1",m_ngch2,m_cntr_btof1);
00749           status = m_tuple12->addIndexedItem ("te_btof1",  m_ngch2,m_te_btof1);
00750           status = m_tuple12->addIndexedItem ("tmu_btof1", m_ngch2,m_tmu_btof1);
00751           status = m_tuple12->addIndexedItem ("tpi_btof1", m_ngch2,m_tpi_btof1);
00752           status = m_tuple12->addIndexedItem ("tk_btof1",  m_ngch2,m_tk_btof1);
00753           status = m_tuple12->addIndexedItem ("tp_btof1",  m_ngch2,m_tp_btof1);
00754           status = m_tuple12->addIndexedItem ("ph_btof1",  m_ngch2,m_ph_btof1);
00755           status = m_tuple12->addIndexedItem ("zhit_btof1",m_ngch2,m_zhit_btof1);
00756           status = m_tuple12->addIndexedItem ("qual_btof1",m_ngch2,m_qual_btof1);
00757           status = m_tuple12->addIndexedItem ("b1_toff_e", m_ngch2,m_b1_toff_e);
00758           status = m_tuple12->addIndexedItem ("b1_toff_mu",m_ngch2,m_b1_toff_mu);
00759           status = m_tuple12->addIndexedItem ("b1_toff_pi",m_ngch2,m_b1_toff_pi);
00760           status = m_tuple12->addIndexedItem ("b1_toff_k", m_ngch2,m_b1_toff_k);
00761           status = m_tuple12->addIndexedItem ("b1_toff_p", m_ngch2,m_b1_toff_p);
00762           status = m_tuple12->addIndexedItem ("b1_tsig_e", m_ngch2,m_b1_tsig_e);
00763           status = m_tuple12->addIndexedItem ("b1_tsig_mu",m_ngch2,m_b1_tsig_mu);
00764           status = m_tuple12->addIndexedItem ("b1_tsig_pi",m_ngch2,m_b1_tsig_pi);
00765           status = m_tuple12->addIndexedItem ("b1_tsig_k", m_ngch2,m_b1_tsig_k);
00766           status = m_tuple12->addIndexedItem ("b1_tsig_p", m_ngch2,m_b1_tsig_p);
00767           status = m_tuple12->addIndexedItem ("b1_tof",    m_ngch2,m_b1_tof);                                                    
00768 
00769           status = m_tuple12->addIndexedItem ("mdedx_pid", m_ngch2,m_dedx_pid);
00770           status = m_tuple12->addIndexedItem ("mtof1_pid", m_ngch2,m_tof1_pid);
00771           status = m_tuple12->addIndexedItem ("mtof2_pid", m_ngch2,m_tof2_pid);
00772           status = m_tuple12->addIndexedItem ("mprob_pid", m_ngch2,m_prob_pid);
00773           status = m_tuple12->addIndexedItem ("mptrk_pid", m_ngch2,m_ptrk_pid);
00774           status = m_tuple12->addIndexedItem ("mcost_pid", m_ngch2,m_cost_pid);
00775         }
00776         else    { 
00777           log << MSG::ERROR << "    Cannot book N-tuple:" << long(m_tuple12) << endmsg;
00778           return StatusCode::FAILURE;
00779         }
00780   }
00781 
00782   NTuplePtr nt13(ntupleSvc(), "DQAFILE/CheckTofNeAll");
00783   if ( nt13 ) m_tuple13 = nt13;
00784   else {
00785         m_tuple13 = ntupleSvc()->book ("DQAFILE/CheckTofNeAll", CLID_ColumnWiseTuple, "ks N-Tuple example");
00786         if ( m_tuple13 )    {
00787           status = m_tuple13->addItem ("ngneu",  m_nggneu, 0, 100);
00788           status = m_tuple13->addIndexedItem ("shwidtof",       m_nggneu,m_shwidtof);
00789           status = m_tuple13->addIndexedItem ("numHits",      m_nggneu,m_numHits);    // Total number of hits
00790           status = m_tuple13->addIndexedItem ("secondmoment", m_nggneu,m_secondmoment);
00791           status = m_tuple13->addIndexedItem ("mx",           m_nggneu,m_x);       //  Shower coordinates and errors
00792           status = m_tuple13->addIndexedItem ("my",           m_nggneu,m_y);
00793           status = m_tuple13->addIndexedItem ("mz",           m_nggneu,m_z);
00794           status = m_tuple13->addIndexedItem ("cosemc",       m_nggneu,m_cosemc);   // Shower Counter angles and errors
00795           status = m_tuple13->addIndexedItem ("phiemc",       m_nggneu,m_phiemc);
00796           status = m_tuple13->addIndexedItem ("energy",       m_nggneu,m_energy);  // Total energy observed in Emc
00797           status = m_tuple13->addIndexedItem ("eseed",        m_nggneu,m_eSeed);
00798           status = m_tuple13->addIndexedItem ("me9",          m_nggneu,m_e3x3); 
00799           status = m_tuple13->addIndexedItem ("me25",         m_nggneu,m_e5x5); 
00800           status = m_tuple13->addIndexedItem ("mlat",         m_nggneu,m_lat);
00801           status = m_tuple13->addIndexedItem ("ma20",         m_nggneu,m_a20);
00802           status = m_tuple13->addIndexedItem ("ma42",         m_nggneu,m_a42);  
00803         }
00804         else    { 
00805           log << MSG::ERROR << "    Cannot book N-tuple:" << long(m_tuple13) << endmsg;
00806           return StatusCode::FAILURE;
00807         }
00808   }
00809 
00810   if(service("THistSvc", m_thsvc).isFailure()) {
00811         log << MSG::ERROR << "Couldn't get THistSvc" << endreq;
00812         return StatusCode::FAILURE;
00813   }
00814 
00815   string strdecaymode[11] = {"D02KPi","D02KPiPi0","D02KPiPiPi","D02KPiPiPiPi0","D02KsPiPi","D02KsPiPiPi0",\
00816   "DC2KPiPi","DC2KPiPiPi0","DC2KsPiPi0","DC2KsPiPiPi","DC2KKPi"};
00817   for (int i=0;i!=11;i++)
00818   {
00819   regDecayModeHistogram(m_thsvc,strdecaymode[i]);
00820   }
00821   string strDdecaymode[2] = {"D0","DC"};
00822   regDecayModeHistogram(m_thsvc,strDdecaymode[0]);
00823   regDecayModeHistogram(m_thsvc,strDdecaymode[1]);
00824   string strtrack[14] = {"D02KPi_kaon","D02KPi_pion","D02KPiPi0_kaon","D02KPiPi0_pion","D02KsPiPi_pion1",\
00825         "D02KsPiPi_pion2","D02KsPiPi_pion3","D02KsPiPi_pion4","DC2KPiPi_kaon","DC2KPiPi_pion1","DC2KPiPi_pion2",\
00826           "DC2KsPiPi0_pion1","DC2KsPiPi0_pion2","DC2KsPiPi0_pion3" };
00827   for (int i=0;i!=14;i++)
00828   {
00829         regTrackHistogram(m_thsvc,strtrack[i]);
00830   }
00831   string strphoton[2] = {"D02KPiPi0_photon1","D02KPiPi0_photon2"};
00832   regPhotonHistogram(m_thsvc,strphoton[0]);
00833   regPhotonHistogram(m_thsvc,strphoton[1]);
00834   string strtrackdedx[4] = {"D02KPi_kaon","D02KPi_pion","D02KPiPi0_kaon","D02KPiPi0_pion"};
00835   for (int i=0;i!=4;i++)
00836   {
00837         regTrackDedxHistogram(m_thsvc,strtrackdedx[i]);
00838   }
00839   string strtracktof[4] = {"D02KPi_kaon","D02KPi_pion","D02KPiPi0_kaon","D02KPiPi0_pion"};
00840   for (int i=0;i!=4;i++)
00841   {
00842         regTrackTofHistogram(m_thsvc,strtrackdedx[i]);
00843   }
00844   string strpi0[2] = {"D02KPiPi0","DC2KsPiPi0"};
00845   for (int i=0;i!=2;i++)
00846   {
00847         regPi0Histogram(m_thsvc,strpi0[i]);
00848   }
00849   string strks[2] = {"D02KsPiPi","DC2KsPiPi0"};
00850   for (int i=0;i!=2;i++)
00851   {
00852         regKsHistogram(m_thsvc,strks[i]);
00853   }
00854   //
00855   //--------end of book--------
00856   //
00857 
00858   log << MSG::INFO << "successfully return from initialize()" <<endmsg;
00859   return StatusCode::SUCCESS;
00860 
00861 }

DQADtag Class Reference

Public Member Functions

Private Member Functions

Private Attributes

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation


Public Member Functions
	DQADtag (const std::string &name, ISvcLocator *pSvcLocator)
	DQADtag (const std::string &name, ISvcLocator *pSvcLocator)
StatusCode	execute ()
StatusCode	execute ()
StatusCode	finalize ()
StatusCode	finalize ()
StatusCode	initialize ()
StatusCode	initialize ()
Private Member Functions
void	fillDecayModeHistogram (ITHistSvc *&, std::string str, int size)
void	fillDecayModeHistogram (ITHistSvc *&, std::string str, int size)
void	fillDecayModeItem (DTagToolIterator iter_dtag, string str, int indexD)
void	fillDecayModeItem (DTagToolIterator iter_dtag, string str, int indexD)
void	fillHistogram (ITHistSvc &, string str, int size, int NUM, std::string fillHistName, int tuplenum, bool cut)
void	fillHistogram (ITHistSvc &, string str, int size, int NUM, std::string fillHistName, NTuple::Tuple *&m_tuple, int pi0cut, int kscut)
void	fillHistogram (ITHistSvc &, string str, int size, int NUM, std::string fillHistName, NTuple::Tuple *&m_tuple)
void	fillHistogram (ITHistSvc &, string str, int size, int NUM, std::string fillHistName, int tuplenum)
void	fillHistogram (ITHistSvc &, string str, int size, int NUM, std::string fillHistName, int tuplenum, bool cut)
void	fillHistogram (ITHistSvc &, string str, int size, int NUM, std::string fillHistName, NTuple::Tuple *&m_tuple, int pi0cut, int kscut)
void	fillHistogram (ITHistSvc &, string str, int size, int NUM, std::string fillHistName, NTuple::Tuple *&m_tuple)
void	fillHistogram (ITHistSvc &, string str, int size, int NUM, std::string fillHistName, int tuplenum)
void	fillKsHistogram (ITHistSvc *&, std::string str, int size)
void	fillKsHistogram (ITHistSvc *&, std::string str, int size)
void	fillKsItem (SmartRefVector< EvtRecTrack > tracks, string str, int indexD, DQAKsInfo *ksinfo)
void	fillKsItem (SmartRefVector< EvtRecTrack > tracks, string str, int indexD, DQAKsInfo *ksinfo)
void	fillPhotonItem (EvtRecTrack *shower, string str, int indexD)
void	fillPhotonItem (EvtRecTrack *shower, string str, int indexD)
void	fillPi0Histogram (ITHistSvc *&, std::string str, int size)
void	fillPi0Histogram (ITHistSvc *&, std::string str, int size)
void	fillPi0Item (SmartRefVector< EvtRecTrack > showers, string str, int indexD, DQAPi0Info *pi0info)
void	fillPi0Item (SmartRefVector< EvtRecTrack > showers, string str, int indexD, DQAPi0Info *pi0info)
void	fillTofHistogram (ITHistSvc *&, string str, int size)
void	fillTofHistogram (ITHistSvc *&, string str, int size)
void	fillTrackDedxHistogram (ITHistSvc *&, std::string str, int size)
void	fillTrackDedxHistogram (ITHistSvc *&, std::string str, int size)
int	fillTrackDedxItem (EvtRecTrack *track, string str, int type, int indexDedx)
int	fillTrackDedxItem (EvtRecTrack *track, string str, int type, int indexDedx)
void	fillTrackHistogram (ITHistSvc *&, std::string str, int size)