Beispiel #1
0
void FitterUtils::prepare_PDFs(string trigStr, string BDTVar, double BDTcut,
      string signalfile, string partrecofile, string combfile, string JpsiLeakfile,
      double minBMass, double maxBMass,
      string signaltree, string partrecotree, string combtree, string JpsiLeaktree)
{


   //***********Get the datasets
   TFile* fSignal = new TFile(signalfile.c_str());
   TTree* tSignal = (TTree*)fSignal->Get(signaltree.c_str());
   TFile* fPartReco = new TFile(partrecofile.c_str());
   TTree* tPartReco = (TTree*)fPartReco->Get(partrecotree.c_str());
   TFile* fComb = new TFile(combfile.c_str());
   TTree* tComb = (TTree*)fComb->Get(combtree.c_str()); 
   TFile* fJpsiLeak = new TFile(JpsiLeakfile.c_str());
   TTree* tJpsiLeak = (TTree*)fJpsiLeak->Get(JpsiLeaktree.c_str()); 


   //**********Define variables
   RooRealVar trigVar(trigStr.c_str(), trigStr.c_str(), -10, 10);
   RooRealVar BDTRooRealVar(BDTVar.c_str(), BDTVar.c_str(), -1,1);
   RooRealVar B_plus_M("B_plus_M", "M_{visible}", minBMass, maxBMass, "MeV");
   RooRealVar misPT("misPT", "p_{#perp}", 0, 5000, "MeV");
   RooRealVar B_plus_DTFM_M_zero("B_plus_DTFM_M_zero", "M_{constr}", 0, 20000, "MeV"); 
   RooRealVar e_plus_BremMultiplicity("e_plus_BremMultiplicity","e_plus_BremMultiplicity", -1,2);
   RooRealVar e_minus_BremMultiplicity("e_minus_BremMultiplicity","e_minus_BremMultiplicity", -1,2);

   RooRealVar weightPartReco("weightPartReco", "weightPartReco", 0, 10);
   RooRealVar weightLeakage("weightLeakage", "weightLeakage", 0, 10);
   RooRealVar dataMCWeightee("DataMCWeightee", "DataMCWeightee",0, 30);


   //***********Set only variables needed

   tSignal->SetBranchStatus("*", 0); tSignal->SetBranchStatus("B_plus_M", 1); tSignal->SetBranchStatus("misPT", 1); tSignal->SetBranchStatus("B_plus_DTFM_M_zero", 1); tSignal->SetBranchStatus(BDTVar.c_str(),1);
   tSignal->SetBranchStatus("e_plus_BremMultiplicity", 1); tSignal->SetBranchStatus("e_minus_BremMultiplicity", 1); tSignal->SetBranchStatus(trigStr.c_str()); tSignal->SetBranchStatus("DataMCWeightee",1);

   tPartReco->SetBranchStatus("*", 0); tPartReco->SetBranchStatus("B_plus_M", 1); tPartReco->SetBranchStatus("misPT", 1); tPartReco->SetBranchStatus("B_plus_DTFM_M_zero", 1);tPartReco->SetBranchStatus(BDTVar.c_str(),1);
   tPartReco->SetBranchStatus("e_plus_BremMultiplicity", 1); tPartReco->SetBranchStatus("e_minus_BremMultiplicity", 1); tPartReco->SetBranchStatus(trigStr.c_str()); tPartReco->SetBranchStatus("weightPartReco",1);

   tComb->SetBranchStatus("*", 0); tComb->SetBranchStatus("B_plus_M", 1); tComb->SetBranchStatus("misPT", 1); tComb->SetBranchStatus("B_plus_DTFM_M_zero", 1);tComb->SetBranchStatus(BDTVar.c_str(),1);
   tComb->SetBranchStatus("e_plus_BremMultiplicity", 1); tComb->SetBranchStatus("e_minus_BremMultiplicity", 1); tComb->SetBranchStatus(trigStr.c_str());

   tJpsiLeak->SetBranchStatus("*", 0); tJpsiLeak->SetBranchStatus("B_plus_M", 1); tJpsiLeak->SetBranchStatus("misPT", 1); 
   tJpsiLeak->SetBranchStatus("B_plus_DTFM_M_zero", 1);tJpsiLeak->SetBranchStatus(BDTVar.c_str(),1);
   tJpsiLeak->SetBranchStatus("e_plus_BremMultiplicity", 1); tJpsiLeak->SetBranchStatus("e_minus_BremMultiplicity", 1); tJpsiLeak->SetBranchStatus(trigStr.c_str()); 
   tJpsiLeak->SetBranchStatus("weightLeakage",1);


   //***********Set Binning


   RooBinning defaultMBins(floor((maxBMass-minBMass)/(40.)), B_plus_M.getMin(), B_plus_M.getMax() ); 
   RooBinning defaultMisPTBins(floor(40), misPT.getMin(), misPT.getMax()); 
   RooBinning broaderMBins(floor((maxBMass-minBMass)/(80.)), B_plus_M.getMin(), B_plus_M.getMax()); 
   RooBinning broaderMisPTBins(floor(40), misPT.getMin(), misPT.getMax()); 

   B_plus_M.setBinning( defaultMBins);
   misPT.setBinning( defaultMisPTBins );
   B_plus_M.setBinning( broaderMBins, "broaderBins");
   misPT.setBinning( broaderMisPTBins, "broaderBins" );

   B_plus_DTFM_M_zero.setBins(100);

   RooArgSet argset(BDTRooRealVar, B_plus_DTFM_M_zero, misPT,  B_plus_M, trigVar, e_plus_BremMultiplicity, e_minus_BremMultiplicity);
   RooArgSet argsetPartReco(BDTRooRealVar, B_plus_DTFM_M_zero, misPT,  B_plus_M, trigVar, e_plus_BremMultiplicity, e_minus_BremMultiplicity, weightPartReco);
   RooArgSet argsetLeakage(BDTRooRealVar, B_plus_DTFM_M_zero, misPT,  B_plus_M, trigVar, e_plus_BremMultiplicity, e_minus_BremMultiplicity, weightLeakage);
   RooArgSet argsetSignal(BDTRooRealVar, B_plus_DTFM_M_zero, misPT,  B_plus_M, trigVar, e_plus_BremMultiplicity, e_minus_BremMultiplicity, dataMCWeightee);

   cout<<"getting the datasets:"<<endl;

   RooDataSet* dataSetSignalZeroGamma;
   RooDataSet* dataSetSignalOneGamma;
   RooDataSet* dataSetSignalTwoGamma;
   RooDataSet* dataSetPartReco;
   RooDataSet* dataSetJpsiLeak;
   RooDataSet* dataSetComb;

   TFile* fw(NULL);

   string BDTCutString( ("("+BDTVar+">"+d2s(BDTcut)+")").c_str()  );

   dataSetSignalZeroGamma = new RooDataSet("dataSetSignalZeroGamma", "dataSetSignalZeroGamma", argsetSignal, Import(*tSignal), Cut(( " ("+trigStr+"  > 0.9) && "+BDTCutString+" && ((e_plus_BremMultiplicity+e_minus_BremMultiplicity) > -0.5) && ((e_plus_BremMultiplicity+e_minus_BremMultiplicity) < 0.5) && B_plus_M > "+d2s(minBMass)+" && B_plus_M < "+d2s(maxBMass)).c_str()), WeightVar("DataMCWeightee")  );

   dataSetSignalOneGamma = new RooDataSet("dataSetSignalOneGamma", "dataSetSignalOneGamma", argsetSignal, Import(*tSignal), Cut(( " ("+trigStr+"  > 0.9) && "+BDTCutString+" && ((e_plus_BremMultiplicity+e_minus_BremMultiplicity) > 0.5) && ((e_plus_BremMultiplicity+e_minus_BremMultiplicity) < 1.5) && B_plus_M > "+d2s(minBMass)+" && B_plus_M < "+d2s(maxBMass)).c_str()), WeightVar("DataMCWeightee")  );

   dataSetSignalTwoGamma = new RooDataSet("dataSetSignalTwoGamma", "dataSetSignalTwoGamma", argsetSignal, Import(*tSignal), Cut(( " ("+trigStr+"  > 0.9) && "+BDTCutString+" && ((e_plus_BremMultiplicity+e_minus_BremMultiplicity) > 1.5) && ((e_plus_BremMultiplicity+e_minus_BremMultiplicity) < 2.5) && B_plus_M > "+d2s(minBMass)+" && B_plus_M < "+d2s(maxBMass)).c_str()), WeightVar("DataMCWeightee")  );


   dataSetPartReco = new RooDataSet("dataSetPartReco", "dataSetPartReco",  argsetPartReco, Import(*tPartReco),Cut(("("+trigStr+"  > 0.9) && "+BDTCutString+ " && B_plus_M > "+d2s(minBMass)+" && B_plus_M < "+d2s(maxBMass)).c_str()), WeightVar("weightPartReco"));

   dataSetJpsiLeak = new RooDataSet("dataSetJpsiLeak", "dataSetJpsiLeak",  argsetLeakage, Import(*tJpsiLeak),Cut(("B_plus_M > "+d2s(minBMass)+" && B_plus_M < "+d2s(maxBMass)).c_str()), WeightVar("weightLeakage"));

    // dataSetComb = new RooDataSet("dataSetComb", "dataSetComb", tComb, argset, ("("+trigStr+"  > 0.9) && (UBDT3R > "+d2s(BDTcut-0.03)+")  && B_plus_M > "+d2s(minBMass)+" && B_plus_M < "+d2s(maxBMass)).c_str());
   dataSetComb = new RooDataSet("dataSetComb", "dataSetComb", tComb, argset, ("("+trigStr+"  > 0.9) && "+BDTCutString+"  && B_plus_M > "+d2s(minBMass)+" && B_plus_M < "+d2s(maxBMass)).c_str());

   cout<<"Number of zero: "<< dataSetSignalZeroGamma->sumEntries()<<endl;
   cout<<"Number of one: "<< dataSetSignalOneGamma->sumEntries()<<endl;
   cout<<"Number of two: "<< dataSetSignalTwoGamma->sumEntries()<<endl;
   cout<<"Number of PartReco: "<< dataSetPartReco->sumEntries()<<endl;
   cout<<"Number of Jpsi leaking:"<< dataSetJpsiLeak->sumEntries()<<endl;
   cout<<"Number of combinatorial events:"<< dataSetComb->sumEntries()<<endl;


   cout<<"binning the datasets:"<<endl;

   RooArgSet argset2(B_plus_M);
   if (fit2D) argset2.add(misPT);

   RooDataHist dataHistSignalZeroGamma("dataHistSignalZeroGamma", "dataHistSignalZeroGamma", argset2, *dataSetSignalZeroGamma); 
   RooDataHist dataHistSignalOneGamma("dataHistSignalOneGamma", "dataHistSignalOneGamma", argset2, *dataSetSignalOneGamma); 
   RooDataHist dataHistSignalTwoGamma("dataHistSignalTwoGamma", "dataHistSignalTwoGamma", argset2, *dataSetSignalTwoGamma); 
   RooDataHist dataHistComb("dataHistComb", "dataHistComb", argset2, *dataSetComb); 
   RooDataHist dataHistPartReco("dataHistPartReco", "dataHistPartReco", argset2, *dataSetPartReco); 
   RooDataHist dataHistJpsiLeak("dataHistJpsiLeak", "dataHistJpsiLeak", argset2, "broaderBins");
   dataHistJpsiLeak.add(*dataSetJpsiLeak); 

   //*************** Compute Error on J/psi leak

   double ErrorJpsi(0);
   if(dataSetJpsiLeak->sumEntries(("("+trigStr+"  > 0.9) && "+BDTCutString+" && B_plus_M > "+d2s(minBMass)+" && B_plus_M < "+d2s(maxBMass)).c_str()) > 0) ErrorJpsi = 1./sqrt(dataSetJpsiLeak->sumEntries(("("+trigStr+"  > 0.9) && "+BDTCutString+" && B_plus_M > "+d2s(minBMass)+" && B_plus_M < "+d2s(maxBMass)).c_str()));
   RooRealVar fractionalErrorJpsiLeak("fractionalErrorJpsiLeak", "fractionalErrorJpsiLeak", ErrorJpsi);
   cout<<"JPSI LEAK: "<<dataSetJpsiLeak->sumEntries(("("+trigStr+"  > 0.9) && "+BDTCutString+" && B_plus_M > "+d2s(minBMass)+" && B_plus_M < "+d2s(maxBMass)).c_str());
   cout<<"JPSI LEAK fractional Error: "<<ErrorJpsi<<endl;


   //***************Create 2D histogram estimates from data


   cout<<"Preparing the 3 2D histPdf: 1";
   //   RooArgSet argset2(B_plus_M);
   RooHistPdf histPdfSignalZeroGamma("histPdfSignalZeroGamma", "histPdfSignalZeroGamma", argset2, dataHistSignalZeroGamma,2); cout<<" 2";
   RooHistPdf histPdfSignalOneGamma("histPdfSignalOneGamma", "histPdfSignalOneGamma", argset2, dataHistSignalOneGamma,2); cout<<" 3";
   RooHistPdf histPdfSignalTwoGamma("histPdfSignalTwoGamma", "histPdfSignalTwoGamma", argset2, dataHistSignalTwoGamma,2); cout<<" 4";
   RooHistPdf histPdfPartReco("histPdfPartReco", "histPdfPartReco", argset2, dataHistPartReco,2); cout<<" 5";
   RooHistPdf histPdfJpsiLeak("histPdfJpsiLeak", "histPdfJpsiLeak", argset2, dataHistJpsiLeak,2); cout<<" 6";


   //***************Create combinatorial from fit to data

   RooRealVar expoConst("expoConst", "expoConst", -1e-3, -1, 1);
   RooRealVar T("T", "T", 97, 0, 200);
   RooRealVar n("n", "n", 3.5, 1., 5.5);
   RooAbsPdf *combPDF;

   if (fit2D)
   {  
      combPDF =  new RooPTMVis("PTMVis", "PTMVis", misPT, B_plus_M, T, n, expoConst);
   }
   else
   {
      combPDF =  new RooExponential("histPdfComb", "histPdfComb", B_plus_M, expoConst);
   }



   combPDF->fitTo(*dataSetComb); // 


   if (fit2D)
   {
      T.setConstant(true);
      n.setConstant(true);
      std::cout<<"T generated is: "<<T.getVal()<<std::endl;
   }

   RooRealVar trueExp("trueExp","trueExp", expoConst.getVal(), -1000, 1000);
   trueExp.setError(expoConst.getError());

   RooRealVar trueT("trueT","trueT", T.getVal(), -1000, 1000);
   trueT.setError(T.getError());
   RooRealVar trueN("trueN","trueN", n.getVal(), -1000, 1000);
   trueN.setError(n.getError());


   //***************Save everything on a workspace
   RooWorkspace workspace("workspace", "workspace");
   workspace.import(B_plus_DTFM_M_zero);
   workspace.import(B_plus_M);
   workspace.import(misPT);
   workspace.import(expoConst);
   workspace.import(trueExp);
   workspace.import(T);
   workspace.import(n);
   workspace.import(*dataSetSignalZeroGamma);
   workspace.import(*dataSetSignalOneGamma);
   workspace.import(*dataSetSignalTwoGamma);
   workspace.import(*dataSetPartReco);
   workspace.import(*dataSetComb);
   workspace.import(*dataSetJpsiLeak);
   workspace.import(histPdfSignalZeroGamma);
   workspace.import(histPdfSignalOneGamma);
   workspace.import(histPdfSignalTwoGamma);
   workspace.import(histPdfPartReco);
   workspace.import(histPdfJpsiLeak);
   workspace.import(fractionalErrorJpsiLeak);
   workspace.import(trueT);
   workspace.import(trueN);
   workspace.writeToFile(workspacename.c_str());



   delete fComb;
   delete fSignal;
   delete fPartReco;
   if(fw!=0 && fw != NULL) delete fw;
   delete combPDF;

   delete dataSetSignalZeroGamma; 
   delete dataSetSignalOneGamma;
   delete dataSetSignalTwoGamma;
   delete dataSetPartReco;
   delete dataSetJpsiLeak;
   delete dataSetComb;
}
void CEneH0102NightSpider::SeqCrash()
{
	//-----------------------------------------------------
	//破壊された
	//-----------------------------------------------------
	Sint32 ax,ay;

	if( m_bDead )
	{
		m_bDisp = gxFalse;
	}

	if( m_sCrashSeq == 0 )
	{
		m_pEffFire = new CEffFire( m_Pos.x , m_Pos.y-8*100 , PRIO_BG2+1);
		viiMus::PlaySound( enSoundExplosionMiddle );
		viiMus::PlaySound( enSoundExplosionLarge );
	}

	if( m_pEffFire )
	{
		m_pEffFire->m_Pos.x = m_Pos.x;
		m_pEffFire->m_Pos.y = m_Pos.y;
		m_pEffFire->m_Pos.z = PRIO_ENEMY;
	}

	if( m_sCrashSeq <= 96 )
	{
		new CEffectBlackFire ( m_Pos.x + (-64+viiSub::Rand()%128)*100, m_Pos.y  + (-32+viiSub::Rand()%64)*100, NULL );

		m_pEnemyAnime->PosCalc();
		Gimmick();

		m_stObj.mov.x += (0-m_stObj.mov.x)/30;

		m_stObj.pos.x += m_stObj.mov.x;
		m_stObj.pos.y += m_stObj.mov.y;

		m_Pos.x = m_stObj.pos.x;
		m_Pos.y = m_stObj.pos.y;

		m_HitShield[3].SetHantei( ID_ENEMY_DEF , &m_Pos );
		m_HitKurai.SetHantei( ID_ENEMY_DEF , &m_Pos );

		if(m_sCrashSeq%8==0)
		{
			ax = m_Pos.x        +(-8+viiSub::Rand()%16)*1400;
			ay = m_Pos.y-128*100 +(-8+viiSub::Rand()%16)*800;

			for(Sint32 n=0; n<4; n++ )
			{
				viiEff::SetBombDeluxe( ax , ay , 1 , 1+viiSub::Rand()%2 );
				//EffectMakeBigBomb( ax+(-32+viiSub::Rand()%64)*100, ay+(-32+viiSub::Rand()%64)*100 );
			}
			viiMus::PlaySound( enSoundExplosionMiddle );
		}
	}

	if( m_sCrashSeq == 96 )
	{
		//どかーーーん

		if( !m_bDead )
		{
			ax = m_Pos.x;
			ay = m_Pos.y-128*100;

			viiEff::Set(EFF_BIGBOMB , m_Pos.x, m_Pos.y , argset(5));

			for(Sint32 ii=0; ii<enMax ;ii++)
			{
				StPartsInfo *p = m_pEnemyAnime->GetParts( ii );

				_opt_[0] = 0;//(Sint32)&SprCEneH0102NightSpider[ p->sDispSprite ];	//64bitポインタ対応
				_opt_[1] = (m_sDirection==DIR_LEFT)? ATR_FLIP_X : ATR_DFLT;
				_opt_[2] = ARGB_DFLT;
				_opt_[3] = p->dr;
				_opt_[4] = 100;
				_opt_[5] = 1;

				viiEff::Set( EFF_DAMAGE_PARTS,m_Pos.x+p->dx*100*m_sDirection , m_Pos.y + p->dy*100 +m_Offset.y , _opt_ ,&SprCEneH0102NightSpider[ p->sDispSprite ] );
				viiEff::SetBombSimple( m_Pos.x+p->dx*100*m_sDirection , m_Pos.y + p->dy*100 +m_Offset.y , 1 , 9 );
			}

			viiEff::SetBombDeluxe( ax , ay , 1 , 1+viiSub::Rand()%2 );

			CEffectDmgprt *p;
			p = new CEffectDmgprt( ax, ay -1600, 180+15 , 800 );
			p = new CEffectDmgprt( ax, ay -1600, 180+65 , 600 );
			p = new CEffectDmgprt( ax, ay -1600, -60    , 800 );
			p = new CEffectDmgprt( ax, ay -1600, -20    , 700 );

			p = new CEffectDmgprt( ax, ay -1600, 180-15 , 800/2 );
			p = new CEffectDmgprt( ax, ay -1600, 180-65 , 600/2 );
			p = new CEffectDmgprt( ax, ay -1600, 60    , 800/2 );
			p = new CEffectDmgprt( ax, ay -1600, 20    , 700/2 );

			delete m_pEffFire;
			m_pEffFire = NULL;
			pGame->QuakeOn( 32 );
			viiMus::PlaySound( enSoundCrashA );
			ScorePlus( enScore );
		}

		m_bDead = gxTrue;
	}

	if( m_bClear )
	{
		SetActionSeq( enActionSeqEnd );
	}

	m_sCrashSeq ++;
}
void CEneEx0006Tank::LogicAI()
{
	//-----------------------------------------------------
	//移動ロジック
	//-----------------------------------------------------

	if( viiSub::GameCounter()%4 == 0 )
	{
		viiEff::Set(EFF_ENERGYBALL , m_Pos.x , m_Pos.y , argset(300,120,0,200));
	}

	if( m_sAtackWait > 0 )
	{
		m_sAtackWait --;
	}
	else
	{
		Atack();
	}

	switch( m_sSeq ){
	case enPhase1:
		//通常
		m_bDamage = gxFalse;
		m_fRotation += m_fRotationAdd;
		break;

	case enPhase2:
		//コア
		{
			Float32 fRot = ZERO;
			if( m_fRotationAdd > 0.f )
			{
				m_fRotationAdd -= 0.01f;
				if( m_fRotationAdd < 0.5f )
				{
					m_fRotationAdd = 0.5f;
				}
			}
			fRot = m_fRotationAdd;

			CEffectBlackFire *pBlack;
			Sint32 ax,ay;
			ax = m_Pos.x+(-40+viiSub::Rand()%80)*100;
			ay = m_Pos.y+(-40+viiSub::Rand()%80)*100;

			//煙出しとく

			if( m_bDamage )//|| viiSub::Rand()%16 == 0 )
			{
				fRot = -1.f;
				pBlack = new CEffectBlackFire( ax , ay , NULL );
				pBlack->m_Pos.z = PRIO_BG2+1;
			}

//			m_bDamage = gxFalse;
			m_fRotation += fRot;
		}
		break;
	}

	Gimmick();
}
void CEneH0102NightSpider::GimmickLegOpen()
{
	//----------------------------------------------------
	// 脚部制御
	//----------------------------------------------------

	if( m_bLegOpen )
	{
		if( m_fLegKakudo < 100 )
		{
			m_fLegKakudo += (100-m_fLegKakudo)/10;
			if( m_fLegKakudo > 100 )
			{
				m_fLegKakudo = 100;
			}
		}
	}
	else
	{
		if( m_fLegKakudo > 0 )
		{
			m_fLegKakudo += (0-m_fLegKakudo)/10;
			if( m_fLegKakudo < 0 )
			{
				m_fLegKakudo = 0;
			}
		}
	}

	Sint32 sOpenKaku[3];

	//前足
	sOpenKaku[0] = s_cast<Sint32>(-10*m_fLegKakudo/100);
	//真ん中
	sOpenKaku[1] = s_cast<Sint32>(5+10*m_fLegKakudo/100);
	//後ろ
	sOpenKaku[2] = s_cast<Sint32>(10+10*m_fLegKakudo/100);

	m_pEnemyAnime->GetParts( enMAEASI    )->dr += sOpenKaku[0];
	m_pEnemyAnime->GetParts( enUSIROASI1 )->dr += sOpenKaku[1];
	m_pEnemyAnime->GetParts( enUSIROASI2 )->dr += sOpenKaku[2];

	Sint32 sRot;
	Sint32 sDst;
	Sint32 dx,dy;
	Sint32 ox,oy;
	Sint32 sObj;

	for(Sint32 jj=0;jj<3;jj++)
	{
		sRot = m_stTire[jj].rot+sOpenKaku[jj];
		sDst = m_stTire[jj].dist;
		sObj = enSTER11+jj;

		//親の表示座標
		ox = m_pEnemyAnime->GetParent( sObj )->dx;
		oy = m_pEnemyAnime->GetParent( sObj )->dy;

		//子の表示座標を計算
		dx = ( viiMath::Cos100( sRot )*sDst )/100;
		dy = ( viiMath::Sin100( sRot )*sDst )/100;

		m_pEnemyAnime->GetParts( sObj )->dx = ox+dx;
		m_pEnemyAnime->GetParts( sObj )->dy = oy+dy;

		m_PosTire[jj].x = m_Pos.x+m_Offset.x;
		m_PosTire[jj].y = m_Pos.y+m_Offset.y;

		m_pEnemyAnime->GetPos( &m_PosTire[jj] , enSTER11+jj , m_sDirection );
	}

	m_Offset.y = s_cast<Sint32>( (viiSub::Rand()%4+-128+(16*m_fLegKakudo/100))*100 );

	m_sTireRot += m_sTireAdd;
	m_sTireRot = (36000+m_sTireRot)%36000;

	if( ( m_Add.x > -220 && m_sControl&JOY_L) || ( m_Add.x < 220 && m_sControl&JOY_R))
	{
		for(Sint32 ii=0; ii<3; ii++ )
		{
			new CEffDashStorm( m_PosTire[ii].x+(-32+viiSub::Rand()%64)*100 , m_sGround, 180 , 60 );
		}
		viiEff::Set(EFF_FIREFLOWER , m_PosTire[0].x , m_PosTire[0].y , argset( viiSub::Rand()%360) );
		viiEff::Set(EFF_FIREFLOWER , m_PosTire[1].x , m_PosTire[1].y , argset( viiSub::Rand()%360) );
		viiEff::Set(EFF_FIREFLOWER , m_PosTire[2].x , m_PosTire[2].y , argset( viiSub::Rand()%360) );
	}

}