Particle definition and quality - SpiRIT-Collaboration/FLowAnalysis GitHub Wiki

Particle Definition

 void STKParticle::SetRecoTrack(STRecoTrack *atrack, STRecoTrack *atrackva)
{
 Clear();
 fRTrack  = atrack;  // trackArray
 fVATrack = atrackva; // trackVAArra
 fdEdx   = fVATrack -> GetdEdxWithCut(0, 0.7);
 forigP3 = fVATrack -> GetMomentumTargetPlane(); //v56
 fncl    = fVATrack -> GetClusterIDArray() -> size();
 frmom   = forigP3;
 fP    = forigP3.Mag();
 fChar = fVATrack->GetCharge();
 fGFChar = fVATrack->GetGenfitCharge();
 rVertexID      =  fRTrack -> GetVertexID();
 rHelixID       =  fRTrack -> GetHelixID();
 rdEdxPointSize =  fRTrack -> GetdEdxPointArray() -> size();
 rNDF           =  fRTrack -> GetNDF();
 rPOCAVertex    =  fRTrack -> GetPOCAVertex();
 rChi2          =  fRTrack -> GetChi2();
 rncl           =  fRTrack -> GetClusterIDArray() -> size();  //<- rClusterSize(STParticle.hh)
 // vertex pos
 Bool_t bgtrack =
   (abs(fRTrack->GetPOCAVertex().Z() + 14.85 ) <= 3.* 1.33) &&
   (abs(fRTrack->GetPOCAVertex().X()) <= 15.) &&
   (abs(fRTrack->GetPOCAVertex().Y() + 205.) <= 20. );
 if( !bgtrack )
   SetVertexAtTargetFlag(0);  //fVatTargetf
 SetGoodTrackFlag();
}

Setup Track quality flag

void STSpiRITTPCTask::SetupTrackQualityFlag(STKParticle *apart)
{
 if( bmA == 100 ) {
   cout << " vDistance " << apart->GetDistanceAtVertex() << " NDF " << apart->GetNDF() << endl;
   return ;
 }

if( apart->GetNCL() < 15 )  // since v55
  apart->SetNCLFlag(0); //// fnclf//

UInt_t pids = apart->GetPID_seq();

  if( pids < 7 ) {
    if( apart->GetP() < momRange[pids][0] ) {
      apart->SetMomentumFlag(0);
    }
  }
}
 void STKParticle::SetVertex(TVector3 value)
{
 fvertex = value;
 rDist = (rPOCAVertex - fvertex).Mag();

 if( rDist > 20 ) SetDistanceAtVertexFlag(0); //fdistanceatvertexf
}

Flag definition

void STKParticle::SetGoodTrackFlag()
{
 fgoodtrackf = fgoodtrackf != 0 ?
   1     * fdoublef    +
   10    * fmomentumf  +
   100   * fmassf      +
   1000  * fnclf       +
   10000 * fVatTargetf +
   100000* fdistanceatvertexf
   : 0;
}

Manual for Flow Analysis

Step 1

Reconstructed track files and BDC, NeuLAND, Katana, and Kyoto Array is combined into one files run by run. Track quality flag is set then reaction plane is obtained. If flattening correction database is prepared, corrected parameters are also stored in the output.

$ Edit run_analysis.sh

DBVERSION=25  --> Database version if it doesn't exist, no database is loaded.
VERSION=25    --> Output file version run####_BTt.v25.root
RUNNUMBER1=(${RNF132}) --> $RNF132 is defined in runList.sh

$ source run_analysis.sh

$ execb -→ Execute he first run defined by RUNNUMBER1

$ execc -→ Execute the whole runs from the second.

$ export MEVT= -→ if run for all events, otherwise, export MEVT=##

$ execa -→ test run with the first run

Second step

Calculate flattening correction factors then corrected. Databases of data/mlt_132Sn.v#..root and data/psi_132Sn.v.#.root are created.

$ Edit DoRPAnalysis.sh

export BEAMSN=$RNF108 <<-- or $RNF132 or $RNF112 or $RNF124
export DBVERSION=35
export VERSION=35.0         ##<- output file version
export MNDB=BTt             ##<- output file sufix

$ source DoAna.sh

$ docorrection ← Getting correction factor and correction is done sequentially.

Third step

$ Edit DoFlow.sh

##---->>> EDIT here
export MNMACRO=DoFlow_adv.C        ##<--- MACRO name
export MNRNF={$RNF132}             ##<---
export MNDB=BTt                    ##<---
export MNVERSION=39.0              ##   <------@@ input
export MNOVER=0
##<-----------

$ source DoFlow.sh

$ run -1 ← Psi angle correction and centrality dependence correction will be done.

The Last step

$ runbatch 2 ← 2: proton, 3: deuteron, 4: triton, 5: 3He, 6: 4He

Make plots for v1 and v2 as a function of rapidity and pt for each rapidity bin.