Star-tpc-l mailing list

[许月 <xuyue AT impcas.ac.cn>]

Hi all,

Please find my presentation slides using the link below:

https://drupal.star.bnl.gov/STAR/system/files/Acceptanceindffproduc_NewBichsel.pdf

I used the Bichsel you mentioned recalculate nSigma of triton and helium3:

dEdxexpected = StdEdxPull::EvalPred(poverm, 1, charge)) 
nsigma = StdEdxPull::Eval(dedx_measured,dedx_resolution,betagamma,fit,charge)

Although the  center value  of nsigma calculated by the new Bichsel function is closer to zero, the acceptance boundary in new production (P23ie) is still inconsistent.

Best regards,

Yue

-----原始邮件-----
发件人:"Fisyak, Yuri V" <fisyak AT bnl.gov>
发送时间:2023-12-08 05:27:55 (星期五)
收件人: "Ma, Rongrong" <marr AT bnl.gov>, "Richard Witt" <witt AT usna.edu>
抄送: "Star-tpc L" <Star-tpc-l AT lists.bnl.gov>, xuyue <xuyue AT impcas.ac.cn>
主题: Re: [Star-tpc-l] TPC Meeting Minutes

Hello Rongrong and Yue,

   With the new dE/dx model prediction has to be calculated via StdEdxPull::EvalPred(poverm, 1, charge)); 

See example in $STAR/StRoot/StBichsel/bichselG10.C where

…

//________________________________________________________________________________

Double_t bichselZM(Double_t *x,Double_t *par) {

  Double_t pove   = pOverM(x[0]);

  Double_t scale = 1;

  Double_t mass = par[0];

  if (mass < 0) {mass = - mass; scale = 2;}

  Double_t poverm = pove/mass;

  Double_t charge = 1.;

  Double_t dx2 = 1;

  if (par[1] > 1.0) {

    charge = par[1];

    poverm *= TMath::Abs(charge);

    dx2 = TMath::Log2(5.);

  }

  return TMath::Log10(1e6*scale*StdEdxPull::EvalPred(poverm, 1, charge));

}

…   

Or

#if defined (__TFG__VERSION__)

//_________________

Float_t StPicoTrack::dEdxPull(Float_t mass, UChar_t fit, Int_t charge) const {

  Float_t z = -999.;

  Float_t momentum  = gMom().Mag();

  Float_t betagamma = momentum * TMath::Abs(charge) / mass;

  Float_t dedx_measured, dedx_resolution = -1;

  if (! fit) { // I70

    dedx_measured = 1e-6*dEdx();

    dedx_resolution = dEdxError();

  }

  else if ( fit == 1) {     // Ifit

    dedx_measured = 1e-6*dEdx();

    dedx_resolution = dEdxError();

  }

else {     // dNdx

    dedx_measured = dNdx();

    dedx_resolution = dNdxError();

  }

  if (dedx_resolution <= 0) return z;

  z = StdEdxPull::Eval(dedx_measured,dedx_resolution,betagamma,fit,charge);

  return z;

}

…

Or in $STAR/StRoot/StPicoEvent/StPicoTrack.cxx

…

#if defined (__TFG__VERSION__)

//_________________

Float_t StPicoTrack::dEdxPull(Float_t mass, UChar_t fit, Int_t charge) const {

  Float_t z = -999.;

  Float_t momentum  = gMom().Mag();

  Float_t betagamma = momentum * TMath::Abs(charge) / mass;

  Float_t dedx_measured, dedx_resolution = -1;

  if (! fit) { // I70

    dedx_measured = 1e-6*dEdx();

    dedx_resolution = dEdxError();

  }

  else if ( fit == 1) {     // Ifit

    dedx_measured = 1e-6*dEdx();

    dedx_resolution = dEdxError();

  }

  else {     // dNdx

    dedx_measured = dNdx();

    dedx_resolution = dNdxError();

  }

  if (dedx_resolution <= 0) return z;

  z = StdEdxPull::Eval(dedx_measured,dedx_resolution,betagamma,fit,charge);

  return z;

}

….

 

                                                              Yuri Fisyak

STAR                                           Phone: +1 631 344 3913
Brookhaven National Laboratory  Fax:    +1 631 344 4206
510A/1-161
http://www.star.bnl.gov/~fisyak       E-mail: fisyak AT bnl.gov

 

 

 

From: Star-tpc-l <star-tpc-l-bounces AT lists.bnl.gov> on behalf of Ma, Rongrong via Star-tpc-l <star-tpc-l AT lists.bnl.gov>
Date: Thursday, December 7, 2023 at 1:39 PM
To: Richard Witt <witt AT usna.edu>
Cc: Star-tpc L <Star-tpc-l AT lists.bnl.gov>
Subject: Re: [Star-tpc-l] TPC Meeting Minutes

Hello All

 

You can find the slides from Yue Xu (https://drupal.star.bnl.gov/STAR/system/files/Acceptanceindffproduc_20231208.pdf) and Wei Zhang (https://drupal.star.bnl.gov/STAR/system/files/nsigmaE_Distribution.pdf) here. They posted the links in a small group discussion thread. 

 

Yue also mentioned:

And the expected de/dx (Bichsel fuction) is:

 

              

        1. For triton : dedx_expected = charge*charge*(TMath::Exp(Bichsel::Instance()->GetMostProbableZM(TMath::Log10(p*TMath::Abs(charge)/mass), 1)));

        2. For helium3 : 

            Float_t ParfHe3Mean[4] = {29.3642,-0.983422,0.0958638,0.169811};

            dedx_expected = ParfHe3Mean[0]*TMath::Power(p, ParfHe3Mean[1] + ParfHe3Mean[2]*log(p) + ParfHe3Mean[3]*log(p)*log(p));

 

Best

Rongrong

On Dec 7, 2023, at 1:03 PM, Richard Witt via Star-tpc-l <Star-tpc-l AT lists.bnl.gov> wrote:

 

 

Present: Alexei, Tommy, Gene, Jim, Flemming, Rongrong, Weizhang, Yuri, Yue, Sonoraj, and myself

 

== Hardware

— Monday Alexei fixed the bottle leak (using old style clamps)

— will watch over next 3 months

— Finished moving RDOs to make room for target for TPC survey

— now need to ensure all holes are clean enough for survey

— will request survey be done next week at 2 o’clock meeting today

— some beneficial adjustments made to scaffolding 

 

== Electronics

— nothing new reported this week

 

== Software

— Run 22 SC (pp500) work is ongoing

— some impact due to vertex finding recently due to pedestal correction

— Yuri just finished dE/dx calibration for Run 19 AuAu200

— checked and found no surprises

— Now working on revisiting Irakli’s alignment work with cosmics 

— all five year with iTPC

— lacked some statistics for Run 19 and 20

— recovered those and now making some progress

— see https://drupal.star.bnl.gov/STAR/blog/fisyak/AuAu200GeV2019-dEdx-recalibration-new-model-has-finished

 

== AOB

— Yue Xu present on light nuclei acceptance in Runs 21 and 23

— links to slides will be posted to TPC list

— requesting a 10k production, keeping event.root, with libraries SL21d and SL23e 

— Weizhang presented photonic electrons from  AuAu 14.6 GeV using SL21c and SL23d 

— comparison of nsignmaE (electron) between libraries

— links to slides will be posted to the TPC list 

 

Best Regards,

  Richard

 

Sent from my iPad

-- 
"Education is the most powerful weapon which you can use to change
the world."

--Nelson Mandela

~~~***~~~***~~~***~~~***~~~***~~~***~~~***~~~***~~~***~~~***~~~***~
 Richard Witt, Professor of Physics     Phone (USNA): 410-293-6675
    United States Naval Academy         Email: witt AT usna.edu
~~~***~~~***~~~***~~~***~~~***~~~***~~~***~~~***~~~***~~~***~~~***~

_______________________________________________
Star-tpc-l mailing list
Star-tpc-l AT lists.bnl.gov
https://lists.bnl.gov/mailman/listinfo/star-tpc-l