incllambda Class Reference

#include <incllambda.h>

List of all members.

Public Member Functions
	incllambda (const std::string &name, ISvcLocator *pSvcLocator)
StatusCode	initialize ()
StatusCode	execute ()
StatusCode	finalize ()
Private Attributes
int	m_pass [10]
double	m_vr0cut
double	m_vz0cut
int	m_checkDedx
int	m_checkTof
NTuple::Tuple *	m_tuple1
NTuple::Item< double >	m_np
NTuple::Item< double >	m_npi
NTuple::Item< double >	m_ctau
NTuple::Item< double >	m_lxyz
NTuple::Item< double >	m_elxyz
NTuple::Item< double >	m_chis
NTuple::Item< double >	m_lamRawMass
NTuple::Item< double >	m_plam
ITHistSvc *	m_thsvc

---

Detailed Description

Definition at line 15 of file incllambda.h.

---

Constructor & Destructor Documentation

incllambda::incllambda	(	const std::string &	name,
		ISvcLocator *	pSvcLocator
	)

Definition at line 61 of file incllambda.cxx.

References genRecEmupikp::i, m_checkDedx, m_checkTof, m_pass, m_tuple1, m_vr0cut, and m_vz0cut.

                                                                      :
  Algorithm(name, pSvcLocator) {
  
  m_tuple1 = 0;
  for(int i = 0; i < 10; i++) m_pass[i] = 0;

//Declare the properties  
  declareProperty("Vr0cut", m_vr0cut=20.0);
  declareProperty("Vz0cut", m_vz0cut=50.0);
  declareProperty("CheckDedx", m_checkDedx = 1);
  declareProperty("CheckTof",  m_checkTof = 1);

}

---

Member Function Documentation

StatusCode incllambda::execute

(

)

Definition at line 126 of file incllambda.cxx.

References TrackPool::AddTrack(), VertexFit::AddVertex(), VertexFit::BuildVirtualParticle(), DstMdcKalTrack::charge(), DstMdcTrack::charge(), SecondVertexFit::chisq(), SecondVertexFit::ctau(), SecondVertexFit::decayLength(), SecondVertexFit::decayLengthError(), DstMdcKalTrack::err(), calibUtil::ERROR, EventModel::EvtRec::EvtRecEvent, EventModel::EvtRec::EvtRecTrackCol, SecondVertexFit::Fit(), VertexFit::Fit(), DstMdcKalTrack::helix(), genRecEmupikp::i, Bes_Common::INFO, SecondVertexFit::init(), VertexFit::init(), SecondVertexFit::instance(), VertexFit::instance(), ParticleID::instance(), m_chis, m_ctau, m_elxyz, m_lamRawMass, m_lxyz, m_np, m_npi, m_pass, m_plam, m_thsvc, m_tuple1, m_vr0cut, m_vz0cut, mlam, mp, mpi, msgSvc(), SecondVertexFit::p4par(), pid, DstMdcKalTrack::pion, DstMdcKalTrack::proton, DstMdcKalTrack::px(), DstMdcKalTrack::py(), DstMdcKalTrack::pz(), DstMdcTrack::r(), VertexParameter::setEvx(), DstMdcKalTrack::setPidType(), SecondVertexFit::setVpar(), VertexParameter::setVx(), VertexFit::vpar(), VertexFit::wVirtualTrack(), and DstMdcTrack::z().

                               {
  StatusCode sc = StatusCode::SUCCESS;  

  MsgStream log(msgSvc(), name());
  log << MSG::INFO << "in execute()" << endreq;

  // DQA
  // Add the line below at the beginning of execute()
  //
  setFilterPassed(false);

  m_pass[0] += 1;

  SmartDataPtr<Event::EventHeader> eventHeader(eventSvc(),"/Event/EventHeader");
  SmartDataPtr<EvtRecEvent> evtRecEvent(eventSvc(), EventModel::EvtRec::EvtRecEvent);
  SmartDataPtr<EvtRecTrackCol> evtRecTrkCol(eventSvc(), EventModel::EvtRec::EvtRecTrackCol);

  Vint ilam, ip, ipi, iGood;
  iGood.clear();
  ilam.clear();
  ip.clear();
  ipi.clear();

  Vp4 ppi, pp;
  ppi.clear();
  pp.clear();

  int TotCharge = 0;
  for(int i = 0; i < evtRecEvent->totalCharged(); i++){
    EvtRecTrackIterator itTrk=evtRecTrkCol->begin() + i;
    if(!(*itTrk)->isMdcTrackValid()) continue;
    RecMdcTrack *mdcTrk = (*itTrk)->mdcTrack();
    if(fabs(mdcTrk->z()) >= m_vz0cut) continue;
    if(mdcTrk->r() >= m_vr0cut) continue;
    iGood.push_back(i);
    TotCharge += mdcTrk->charge();
  }
  //
  // Finish Good Charged Track Selection
  //
    int nGood = iGood.size();

  //
  // Charge track number cut
  //

  if((nGood < 2) || (TotCharge!=0)) return sc;

 m_pass[1] += 1;
  //
  // Assign 4-momentum to each charged track
  //
  ParticleID *pid = ParticleID::instance();
  for(int i = 0; i < nGood; i++) {
    EvtRecTrackIterator itTrk = evtRecTrkCol->begin() + iGood[i];
    //    if(pid) delete pid;
    pid->init();
    pid->setMethod(pid->methodProbability());
    pid->setChiMinCut(4);
    pid->setRecTrack(*itTrk);
    pid->usePidSys(pid->useDedx() | pid->useTof1() | pid->useTof2()); // use PID sub-system
    pid->identify(pid->onlyPion() | pid->onlyProton());  // seperater Pion/Kaon
    pid->calculate();
    if(!(pid->IsPidInfoValid())) continue;
//    RecMdcTrack* mdcTrk = (*itTrk)->mdcTrack();
    if(!(*itTrk)->isMdcKalTrackValid()) continue ;
    RecMdcKalTrack* mdcKalTrk = (*itTrk)->mdcKalTrack();

     if(mdcKalTrk->charge() == 0) continue; 

    if(pid->probPion() > 0.001 && (pid->probPion() > pid->probProton())){
      mdcKalTrk->setPidType(RecMdcKalTrack::pion);
      ipi.push_back(iGood[i]);
      HepLorentzVector ptrk;
      ptrk.setPx(mdcKalTrk->px());
      ptrk.setPy(mdcKalTrk->py());
      ptrk.setPz(mdcKalTrk->pz());
      double p3 = ptrk.mag();
      ptrk.setE(sqrt(p3*p3+mpi*mpi));
      ppi.push_back(ptrk);
    }
    if(pid->probProton() > 0.001 && (pid->probProton() > pid->probPion())){
      mdcKalTrk->setPidType(RecMdcKalTrack::proton);
      ip.push_back(iGood[i]);
      HepLorentzVector ptrk;
      ptrk.setPx(mdcKalTrk->px());
      ptrk.setPy(mdcKalTrk->py());
      ptrk.setPz(mdcKalTrk->pz());
      double p3 = ptrk.mag();
      ptrk.setE(sqrt(p3*p3+mp*mp));
      pp.push_back(ptrk);
    }
  }
  
  m_pass[2] += 1;

      int npi = ipi.size();
      int np = ip.size();
       m_npi = npi;
       m_np = np;
  if(npi < 1 || np <1) return sc;

  m_pass[3] += 1;
//
//****** Second Vertex Check************
//
  double chi=999.;
  double delm;

  HepPoint3D vx(0., 0., 0.);
  HepSymMatrix Evx(3, 0);
  double bx = 1E+6;
  double by = 1E+6;
  double bz = 1E+6;
  Evx[0][0] = bx*bx;
  Evx[1][1] = by*by;
  Evx[2][2] = bz*bz;
  VertexParameter vxpar;
  vxpar.setVx(vx);
  vxpar.setEvx(Evx);

   int ipion=-1, iproton=-1;

  VertexFit *vtxfit0 = VertexFit::instance();
  SecondVertexFit *vtxfit = SecondVertexFit::instance();

  for(unsigned int i1 = 0; i1 < ipi.size(); i1++) {
    RecMdcKalTrack *piTrk = (*(evtRecTrkCol->begin()+ipi[i1]))->mdcKalTrack();
     piTrk->setPidType(RecMdcKalTrack::pion);
    WTrackParameter wpitrk(mpi, piTrk->helix(), piTrk->err());

    for(unsigned int i2 = 0; i2 < ip.size(); i2++) {
      RecMdcKalTrack *protonTrk = (*(evtRecTrkCol->begin()+ip[i2]))->mdcKalTrack();
       protonTrk->setPidType(RecMdcKalTrack::proton);
      WTrackParameter wptrk(mp, protonTrk->helix(), protonTrk->err());

       int NCharge = 0;
       NCharge =piTrk->charge() + protonTrk->charge();

//      cout << "TotNcharge = " << NCharge << endl;

       if(NCharge != 0) continue;

      vtxfit0->init();
      vtxfit0->AddTrack(0, wpitrk);
      vtxfit0->AddTrack(1, wptrk);
      vtxfit0->AddVertex(0, vxpar, 0, 1);
      if(!(vtxfit0->Fit(0))) continue;
      vtxfit0->BuildVirtualParticle(0);

      vtxfit->init();
//      vtxfit->setPrimaryVertex(bs);
      vtxfit->AddTrack(0, vtxfit0->wVirtualTrack(0));
      vtxfit->setVpar(vtxfit0->vpar(0));
      if(!vtxfit->Fit()) continue;

      //double mass = vtxfit->p4par().m();
      //if(mass < m_massRangeLower) continue;
      //if(mass > m_massRangeHigher) continue;
      if(vtxfit->chisq() > chi) continue;
//       if(fabs((vtxfit->p4par()).m()-mlam) > delm) continue;
      delm = fabs((vtxfit->p4par()).m()-mlam);
      chi = vtxfit->chisq();
      ipion=ipi[i1];
      iproton=ip[i2];
    }
  }
//
//   Lammda check
//
      if(ipion==-1 || iproton==-1) return sc;
  m_pass[4] += 1;

    RecMdcKalTrack *pionTrk = (*(evtRecTrkCol->begin()+ipion))->mdcKalTrack();
     pionTrk->setPidType(RecMdcKalTrack::pion);
    WTrackParameter wpiontrk(mpi, pionTrk->helix(), pionTrk->err());

    RecMdcKalTrack *protonTrk = (*(evtRecTrkCol->begin()+iproton))->mdcKalTrack();
     protonTrk->setPidType(RecMdcKalTrack::proton);
    WTrackParameter wprotontrk(mp, protonTrk->helix(), protonTrk->err());
      vtxfit0->init();
      vtxfit0->AddTrack(0, wpiontrk);
      vtxfit0->AddTrack(1, wprotontrk);
      vtxfit0->AddVertex(0, vxpar, 0, 1);
      if(!(vtxfit0->Fit(0))) return sc;
      vtxfit0->BuildVirtualParticle(0);

      vtxfit->init();
//      vtxfit->setPrimaryVertex(bs);
      vtxfit->AddTrack(0, vtxfit0->wVirtualTrack(0));
      vtxfit->setVpar(vtxfit0->vpar(0));
      if(!vtxfit->Fit()) return sc;

      m_lamRawMass = vtxfit->p4par().m();
      m_ctau = vtxfit->ctau();
      m_lxyz = vtxfit->decayLength();
      m_elxyz = vtxfit->decayLengthError();
      m_chis = vtxfit->chisq();
      m_plam = vtxfit->p4par().rho();

      if(m_lxyz>3.5){

    TH1 *h(0);
    if (m_thsvc->getHist("/DQAHist/InclLam/InclLam_mass", h).isSuccess()) {
        h->Fill(m_lamRawMass);
    } else {
        log << MSG::ERROR << "Couldn't retrieve incllam_mass" << endreq;
    }

        }
    m_tuple1->write();
//DQA
// set tag and quality
  // Pid: 1 - electron, 2 - muon, 3 - pion, 4 - kaon, 5 - proton
//  (*(evtRecTrkCol->begin()+ipion))->setPartId(3);
//  (*(evtRecTrkCol->begin()+iproton))->setPartId(5);
  (*(evtRecTrkCol->begin()+ipion))->tagPion();
  (*(evtRecTrkCol->begin()+iproton))->tagProton();
   // Quality: defined by whether dE/dx or TOF is used to identify particle
  // 0 - no dE/dx, no TOF (can be used for dE/dx and TOF calibration)
  // 1 - only dE/dx (can be used for TOF calibration)
  // 2 - only TOF (can be used for dE/dx calibration)
  // 3 - Both dE/dx and TOF
  (*(evtRecTrkCol->begin()+ipion))->setQuality(1);
  (*(evtRecTrkCol->begin()+iproton))->setQuality(1);

//--------------------------------------------------
  // Add the line below at the end of execute(), (before return)
  //
  setFilterPassed(true);

  return StatusCode::SUCCESS;
}

StatusCode incllambda::finalize

(

)

Definition at line 362 of file incllambda.cxx.

References Bes_Common::INFO, m_pass, and msgSvc().

                                {

  MsgStream log(msgSvc(), name());
  log << MSG::INFO << "in finalize()" << endmsg;
  log << MSG::INFO << "Total Entries : " << m_pass[0] << endreq;
  log << MSG::INFO << "Qtot Cut : " << m_pass[1] << endreq;
  log << MSG::INFO << "Nks : " << m_pass[2] << endreq;
  log << MSG::INFO << "Good Ks cut : " << m_pass[3] << endreq;
  return StatusCode::SUCCESS;
}

StatusCode incllambda::initialize

(

)

Definition at line 76 of file incllambda.cxx.

References calibUtil::ERROR, Bes_Common::INFO, m_chis, m_ctau, m_elxyz, m_lamRawMass, m_lxyz, m_np, m_npi, m_plam, m_thsvc, m_tuple1, msgSvc(), and ntupleSvc().

                                 {
  MsgStream log(msgSvc(), name());

  log << MSG::INFO << "in initialize()" << endmsg;
  
  StatusCode status;

 if(service("THistSvc", m_thsvc).isFailure()) {
    log << MSG::ERROR << "Couldn't get THistSvc" << endreq;
    return StatusCode::FAILURE;
  }
  // "DQAHist" is fixed
  TH1F* incllambda_mass = new TH1F("InclLam_mass","INCLUSIVE_LAMBDA_MASS",150,1.11,1.125);
   incllambda_mass->GetXaxis()->SetTitle("M_{p#pi} (GeV)");
   incllambda_mass->GetYaxis()->SetTitle("Nentries/0.1MeV/c^{2}");
  if(m_thsvc->regHist("/DQAHist/InclLam/InclLam_mass", incllambda_mass).isFailure()) {
      log << MSG::ERROR << "Couldn't register incllambda_mass" << endreq;
  }

//*****************************************

    NTuplePtr nt1(ntupleSvc(), "DQAFILE/InclLam");
    if ( nt1 ) m_tuple1 = nt1;
    else {
     m_tuple1 = ntupleSvc()->book ("DQAFILE/InclLam", CLID_ColumnWiseTuple, "incllam Ntuple");
      if ( m_tuple1 )    {
        status = m_tuple1->addItem ("npi", m_npi);
        status = m_tuple1->addItem ("np", m_np);
        status = m_tuple1->addItem ("lxyz", m_lxyz);
        status = m_tuple1->addItem ("ctau", m_ctau);
        status = m_tuple1->addItem ("elxyz", m_elxyz);
        status = m_tuple1->addItem ("lamchi", m_chis);
        status = m_tuple1->addItem ("mlamraw", m_lamRawMass);
        status = m_tuple1->addItem ("plam", m_plam);
      }
      else    { 
        log << MSG::ERROR << "    Cannot book N-tuple:" << long(m_tuple1) << endmsg;
        return StatusCode::FAILURE;
      }
    }
  //
  //--------end of book--------
  //

  log << MSG::INFO << "successfully return from initialize()" <<endmsg;
  return StatusCode::SUCCESS;

}

---

Member Data Documentation

int incllambda::m_checkDedx [private]

Definition at line 32 of file incllambda.h.

Referenced by incllambda().

int incllambda::m_checkTof [private]

Definition at line 33 of file incllambda.h.

Referenced by incllambda().

NTuple::Item<double> incllambda::m_chis [private]

Definition at line 42 of file incllambda.h.

Referenced by execute(), and initialize().

NTuple::Item<double> incllambda::m_ctau [private]

Definition at line 39 of file incllambda.h.

Referenced by execute(), and initialize().

NTuple::Item<double> incllambda::m_elxyz [private]

Definition at line 41 of file incllambda.h.

Referenced by execute(), and initialize().

NTuple::Item<double> incllambda::m_lamRawMass [private]

Definition at line 43 of file incllambda.h.

Referenced by execute(), and initialize().

NTuple::Item<double> incllambda::m_lxyz [private]

Definition at line 40 of file incllambda.h.

Referenced by execute(), and initialize().

NTuple::Item<double> incllambda::m_np [private]

Definition at line 37 of file incllambda.h.

Referenced by execute(), and initialize().

NTuple::Item<double> incllambda::m_npi [private]

Definition at line 38 of file incllambda.h.

Referenced by execute(), and initialize().

int incllambda::m_pass[10] [private]

Definition at line 25 of file incllambda.h.

Referenced by execute(), finalize(), and incllambda().

NTuple::Item<double> incllambda::m_plam [private]

Definition at line 44 of file incllambda.h.

Referenced by execute(), and initialize().

ITHistSvc* incllambda::m_thsvc [private]

Definition at line 46 of file incllambda.h.

Referenced by execute(), and initialize().

NTuple::Tuple* incllambda::m_tuple1 [private]

Definition at line 36 of file incllambda.h.

Referenced by execute(), incllambda(), and initialize().

double incllambda::m_vr0cut [private]

Definition at line 28 of file incllambda.h.

Referenced by execute(), and incllambda().

double incllambda::m_vz0cut [private]

Definition at line 29 of file incllambda.h.

Referenced by execute(), and incllambda().

---