Просмотр исходного кода

some new changes in StFemtoDstMaker_HFT and macro for pAu_2015

Nikita лет назад: 6
Родитель
Сommit
0c26b42291

+ 4 - 8
StRoot/StFemtoDstMaker/StFemtoDstMaker_HFT.cxx

@@ -49,10 +49,8 @@ StFemtoDstMaker_HFT::StFemtoDstMaker_HFT(StMuDstMaker *muMaker,
   // initialize single-particle cut variables
   mTrackMomentum[0] = 0.1;
   mTrackMomentum[1] = 2.;
-  mTrackDcaPerp[0] = 0.;
-  mTrackDcaPerp[1] = 2.;
-  mTrackDcaZ[0] = -3.;
-  mTrackDcaZ[1] = 3.;
+  mTrackDca[0] = 0.;
+  mTrackDca[1] = 3.;
   mTrackNHits[0] = 15;
   mTrackNHits[1] = 50;
   mTrackEta[0] = -1.1;
@@ -444,10 +442,8 @@ Bool_t StFemtoDstMaker_HFT::AcceptTrack(StMuTrack *trk, UShort_t vtxInd) {
                         trk->eta() >= mTrackEta[0] && trk->eta() <= mTrackEta[1] &&
                         trk->nHits() >= mTrackNHits[0] && trk->nHits() <= mTrackNHits[1] &&
                         trk->flag() >= mTrackFlag[0] && trk->flag() <= mTrackFlag[1] &&
-                        mDca.perp() >= mTrackDcaPerp[0] &&
-                        mDca.perp() <= mTrackDcaPerp[1] &&
-                        mDca.z() >= mTrackDcaZ[0] &&
-                        mDca.z() <= mTrackDcaZ[1];
+                        mDca.mag() >= mTrackDca[0] &&
+                        mDca.mag() <= mTrackDca[1];
 
   return mIsGoodTrack;
 }

+ 2 - 4
StRoot/StFemtoDstMaker/StFemtoDstMaker_HFT.h

@@ -78,8 +78,7 @@ class StFemtoDstMaker_HFT : public StMaker {
 
   // track
   void SetParticleMomentum(const Float_t& lo, const Float_t& hi) { mTrackMomentum[0]=lo; mTrackMomentum[1]=hi; }
-  void SetTrackDcaPerp(const Float_t& lo, const Float_t& hi)  { mTrackDcaPerp[0]=lo; mTrackDcaPerp[1]=hi; }
-  void SetTrackDcaZ(const Float_t& lo, const Float_t& hi)  { mTrackDcaZ[0]=lo; mTrackDcaZ[1]=hi; }
+  void SetTrackDca(const Float_t& lo, const Float_t& hi)  { mTrackDca[0]=lo; mTrackDca[1]=hi; }
   void SetTrackNHits(const Int_t& lo, const Int_t& hi)    { mTrackNHits[0]=lo; mTrackNHits[1]=hi; }
   void SetTrackEta(const Float_t& lo, const Float_t& hi)  { mTrackEta[0]=lo; mTrackEta[1]=hi; }
   void SetTrackFlag(const Short_t& lo, const Short_t& hi) { mTrackFlag[0]=lo; mTrackFlag[1]=hi;}
@@ -131,8 +130,7 @@ class StFemtoDstMaker_HFT : public StMaker {
   // Single-particle cuts
   //
   Float_t  mTrackMomentum[2];
-  Float_t  mTrackDcaPerp[2];
-  Float_t  mTrackDcaZ[2];
+  Float_t  mTrackDca[2];
   UShort_t mTrackNHits[2];
   Float_t  mTrackEta[2];
   Short_t  mTrackFlag[2];

+ 3 - 5
macros/femtoDst/FemtoDstMaker_pAu200_y2015.C

@@ -30,11 +30,10 @@ void FemtoDstMaker_pAu200_y2015(const Char_t* inFileList = defaultInFile,
   Float_t cVtxYShift = 0.;
   Float_t cVtxVpdVzDiff[2] = {-30., 30.};
   Float_t cParticleMom[2] = {0.15, 5.0};
-  Float_t cTrackDcaZ[2] = {-3., 3.};
-  Float_t cTrackDcaPerp[2] = {0., 2.};
+  Float_t cTrackDcaPerp[2] = {0., 3.};
   Int_t   cTrackNHits[2] = {10, 50};
   Float_t cTrackEta[2] = {-1., 1.};
-  Int_t   cTrackFlag[2] = {200, 700};
+  Int_t   cTrackFlag[2] = {1, 999};
 
   //Load libraries
   cout << "Loading libraries..." << endl;
@@ -88,8 +87,7 @@ void FemtoDstMaker_pAu200_y2015(const Char_t* inFileList = defaultInFile,
   femtoMaker->SetVtxShift(cVtxXShift, cVtxYShift);
   femtoMaker->SetVtxVpdVzDiffCut(cVtxVpdVzDiff[0], cVtxVpdVzDiff[1]);
   femtoMaker->SetParticleMomentum(cParticleMom[0], cParticleMom[1]);
-  femtoMaker->SetTrackDcaZ(cTrackDcaZ[0], cTrackDcaZ[1]);
-  femtoMaker->SetTrackDcaPerp(cTrackDcaPerp[0], cTrackDcaPerp[1]);
+  femtoMaker->SetTrackDca(cTrackDca[0], cTrackDca[1]);
   femtoMaker->SetTrackNHits(cTrackNHits[0], cTrackNHits[1]);
   femtoMaker->SetTrackEta(cTrackEta[0], cTrackEta[1]);
   femtoMaker->SetTrackFlag(cTrackFlag[0], cTrackFlag[1]);