1.MC simulation

      pythia8->Set("HardQCD:gg2ccbar = on");

      pythia8->Set("HardQCD:qqbar2ccbar = on");

      pythia8->Set("HardQCD:gg2bbbar = on");

      pythia8->Set("HardQCD:qqbar2bbbar = on");

      pythia8->Set("421:onMode=0"); // switch off D+/- decays 
      pythia8->Set("-421:onMode=0"); // switch off D+/- decays 
      pythia8->Set("421:onIfMatch=211 -321"); // switch on for decays to pi+/- and pi0
      pythia8->Set("-421:onIfMatch=-211 321"); // switch on for decays to pi+/- and pi0

We also set the vertex Z in flat (-5,5)

  chain->EventLoop(n);

        vz = StarRandom::flat(-5,5);

        cout<<"vx:"<<vx<<" vy:"<<vy<<" vz:"<<vz<<endl;

        primary -> SetVertex( 0,0, vz);

        primary -> SetSigma( 0,0,0);

  primary->event()->Print();

 2.Filter

  while( ( p = (StarGenParticle*)Iterator.Next() ) )

  {

    if(p->GetId() == 421)

    {

      indexDmeson = p->GetIndex();//Find D meson

      DmesonDaughter1 = p->GetFirstDaughter();//D meson's first child's index

     }

  }

 

 

3.BFC, Reconstruction and Simulation

The options are:

set chain=tpcRS,y2014a,AgML,MakeEvent,tables,SimuDb,VFMinuit,pxlFastSim,istFastSim,Idst,BAna,l0,Tree,logger,  Sti,StiPxl,StiIst,StiHftC,-SstIT,-ssdDb,TpcHitMover,TpxClu,bbcSim,btofsim,mtdsim,-tags,emcY2,EEfs,evout,geantout,  IdTruth,geant,big,fzin,clearmem,HftMatTree,mtdMatch 

But the lib now is still a big problem. We need to use starver DEV. And copy the HFT packages to local directory. 

/star/u/amcw7777/pwg/hft/pp_bc/StRoot

And run makelinks.sh* to make the link to HFT package. By now only DEV lib works and STAR announced the newest version but it doesn't work anyway. 
This package I copied comes from lib in July 2014. 

3.1 The Geant Hits Simulation

The simulation of Geant Hits is not incorrect now. We can see the single track can have a lot of hits. And most of them are at wrong position.

4.Analysis

        mAssocMaker = (StAssociationMaker*) GetMaker("StAssociationMaker");

        if (mAssocMaker) {

                mRcTrackMap = mAssocMaker->rcTrackMap();

                mMcTrackMap = mAssocMaker->mcTrackMap();

        }

4.1 The vertex finder

The vertex in reconstruction, has problem. We assume that it's due to the vertex find mistakes D0 as primary vertex. 
The following plot shows that:

The first one shows the ∆R of D0 and primary vertex. The left one is minBias sample and right one is D0 sample. 

And we can enlarge the scale to 0.1 cm.


The conclustion is, the ∆R of Geant MC and simulation are much different. We suppose in some cases, the D0 is reconstructed
as primary vertex. 

4.2 Topological cuts

We have seen DCA_{prm}, but it may not be correct. Our expectation is that DCA is large with D0 sample than minBias sample.
So we cut on 100µm for it.


The next one is ∆R of K and π. Ideally it should be 0 for D0 decay, but we applied DCA_{Kπ}<150µm.


The last one is theta. It should be a very powerful cut. But it's not that good cause the vertex reconstruction. 

4.3 Efficiency after topological cut

Compared with Au-Au, the final topological cut is a little different. 

	p+p	Au+Au
DCAKπ	<150µm	<50µm
DCAprimary	(>50µm)?	>50µm
|η|	<1.0	<1.0
#pixel	2 sets	2
#TPC	>20	>15
#TPC/poss	>0.52	None
cos(θ)	>0.8	>0.98
Δm	None	<35MeV

And the efficiency of the topology cuts are in next plot.

4.4 S/B ratio

5.Submit JOB to Condor

<SandBox installer="ZIP">

             <Package name="mydemox">

                <File>file:/star/u/amcw7777/pwg/MCProduction/Produce/starsim.C</File>

                <File>file:/star/u/amcw7777/pwg/MCProduction/Produce/.sl64_gcc447/</File>

                <File>file:/star/u/amcw7777/pwg/MCProduction/Reconstruction/runReco.csh</File>

             </Package>

 </SandBox>