Home » Groups » STAR Protected

Run15 Jpsi RpA paper

Updated on Fri, 2022-01-07 09:07. Originally created by marr on 2019-07-25 15:59.
This page collects the supporting materials for Run15 J/psi RpA paper using MTD data.

1) Bug fix: in calculating event multiplicity at detector-level in simulations, the multiplicity for the first vertex is always used regardless which vertex is selected later
- Old code with bug

//_____________________________________________________________________________

Int_t StMyJpsiMcMaker::GetYRefMult(StMuDst *MuDst) const

{

        const Int_t nNodes = mMuDst->numberOfPrimaryTracks();

        Int_t nYRefMult = 0;

        for(Int_t i=0;i<nNodes;i++){

                StMuTrack* pMuTrack = mMuDst->primaryTracks(i);

                if(!pMuTrack) continue;

                //Float_t pt  = pMuTrack->pt();

                Float_t eta = pMuTrack->eta(); 

                Float_t dca = pMuTrack->dcaGlobal().mag();

 

                if(TMath::Abs(eta)>=1.0)          continue;

                if(pMuTrack->nHitsFit(kTpcId)<10) continue;

                if(dca>=1.0)                      continue;

                nYRefMult++;

        }

 

        return nYRefMult;

}


- New code

 

//_____________________________________________________________________________

Int_t StMyJpsiMcMaker::GetYRefMult(const int vtx_index) const

{

  mMuDst->setVertexIndex(vtx_index);

  const Int_t nNodes = mMuDst->numberOfPrimaryTracks();

  Int_t nYRefMult = 0;

  for(Int_t i=0;i<nNodes;i++){

    StMuTrack* pMuTrack = mMuDst->primaryTracks(i);

    if(!pMuTrack) continue;

    Float_t pt  = pMuTrack->pt();

    Float_t eta = pMuTrack->eta(); 

    Float_t dca = pMuTrack->dcaGlobal().mag();

    

    if(dca>=1.0)                      continue;

    if(pt<=0.2)                       continue;

    if(TMath::Abs(eta)>=1.0)          continue;

    if(pMuTrack->nHitsFit(kTpcId)<10) continue;

    nYRefMult++;

  }

  mMuDst->setVertexIndex(0);

  return nYRefMult;

}
Groups: