[u/mrichter/AliRoot.git] / ITS / AliITSsimulationFastPoints.cxx

/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                          *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

/*
$Log$
Revision 1.1.2.1  2000/06/11 20:16:05  barbera
New: Fast simulation class for the ITS, class as part of new ITS code
structure.

*/

#include <TParticle.h>
#include "AliITS.h"
#include "AliITSsimulationFastPoints.h"
#include "AliITSstatistics.h"

ClassImp(AliITSsimulationFastPoints)

AliITSsimulationFastPoints::AliITSsimulationFastPoints()
{
  //constructor
  fSx = new AliITSstatistics(2);
  fSz = new AliITSstatistics(2);
}

//----------------------------------------------------------
AliITSsimulationFastPoints::~AliITSsimulationFastPoints()
{
  //destructor
  delete fSx;
  delete fSz;

}

//-------------------------------------------------------------
void AliITSsimulationFastPoints::CreateFastRecPoints(AliITSmodule *mod){
  // Fast points simulator for all of the ITS.
  Int_t   nhit,h,trk,ifirst;
  Float_t x,y,z,t,e;// local coordinate (cm) and time of flight, and dedx.
  Float_t x1,y1,z1;
  AliITShit *hit;

  fSx->Reset(); // Start out with things clearly zeroed
  fSz->Reset(); // Start out with things clearly zeroed
  e = 0.; // Start out with things clearly zeroed
  Double_t weight=1.;
  nhit = mod->GetNhits();
  ifirst = 1;
  for(h=0;h<nhit;h++){
    hit = mod->GetHit(h);
    hit->GetPositionL(x,y,z,t);
    if(ifirst) {x1=x;y1=y;z1=z;}
    e += hit->GetIonization();
    trk = hit->GetTrack();
    fSx->AddValue((Double_t)x,weight);
    fSz->AddValue((Double_t)z,weight);
    ifirst = 0;
    if(hit->StatusExiting()||  // leaving volume
       hit->StatusDisappeared()|| // interacted/decayed...
       hit->StatusStop() // dropped below E cuts.
       ){ // exiting track, write out RecPoint.
      //      if(fSz->GetRMS()>1.E-1) {
      //	TParticle *part = hit->GetParticle();
      //	printf("idpart %d energy %f \n",part->GetPdgCode(),part->Energy());
      //	printf("diffx=%e diffy=%e diffz=%e\n",x-x1,y-y1,z-z1);
      //      }
      switch (mod->GetLayer()){
      case 1: case 2:  // SPDs
	AddSPD(e,mod,trk);
	break;
      case 3: case 4:  // SDDs
	AddSDD(e,mod,trk);
	break;
      case 5: case 6:  // SSDs
	AddSSD(e,mod,trk);
	break;
      } // end switch
      fSx->Reset();
      fSz->Reset();
      e = 0.;
      ifirst = 1;
      continue;
    }// end if
  } // end for h
}
//_______________________________________________________________________
void AliITSsimulationFastPoints::AddSPD(Float_t &e,
					 AliITSmodule *mod,Int_t trackNumber){
  const Float_t kcmTomicron = 1.0e4;
  //  const Float_t kdEdXtoQ = ;
  const Float_t kRMSx = 12.0; // microns ITS TDR Table 1.3
  const Float_t kRMSz = 70.0; // microns ITS TDR Table 1.3
  Float_t a1,a2; // general float.
  AliITSRecPoint rpSPD;
  Int_t *trk = rpSPD.GetTracks();

  trk[0] = trackNumber;
  trk[1] = 0; trk[2] = 0;
  rpSPD.SetX(kcmTomicron*fSx->GetMean());
  rpSPD.SetZ(kcmTomicron*fSz->GetMean());
  rpSPD.SetdEdX(0.0);
  rpSPD.SetQ(1.0);
  a1 = kcmTomicron*fSx->GetRMS(); a1 *= a1; a1 += kRMSx*kRMSx;
  //  if(a1>1.E5) printf("addSPD: layer=%d track #%d dedx=%e sigmaX2= %e ",
  //		    mod->GetLayer(),trackNumber,e,a1);
  rpSPD.SetSigmaX2(a1);
  a2 = kcmTomicron*fSz->GetRMS(); a2 *= a2; a2 += kRMSz*kRMSz;
  //  if(a1>1.E5) printf(" sigmaZ2= %e\n",a2);
  rpSPD.SetSigmaZ2(a2);
  rpSPD.SetProbability(1.0);

  (mod->GetITS())->AddRecPoint(rpSPD);
}
//_______________________________________________________________________
void AliITSsimulationFastPoints::AddSDD(Float_t &e,
					 AliITSmodule *mod,Int_t trackNumber){

  const Float_t kcmTomicron = 1.0e4;
  const Float_t kdEdXtoQ = 2.778e+8; // Boris Batyuna June 10 2000.
  const Float_t kRMSx = 38.0; // microns ITS TDR Table 1.3
  const Float_t kRMSz = 28.0; // microns ITS TDR Table 1.3
  Float_t a1,a2; // general float.
  AliITSRecPoint rpSDD;
  Int_t *trk = rpSDD.GetTracks();

  trk[0] = trackNumber;
  trk[1] = 0; trk[2] = 0;
  rpSDD.SetX(kcmTomicron*fSx->GetMean());
  rpSDD.SetZ(kcmTomicron*fSz->GetMean());
  rpSDD.SetdEdX(e);
  rpSDD.SetQ(kdEdXtoQ*e);
  a1 = kcmTomicron*fSx->GetRMS(); a1 *= a1; a1 += kRMSx*kRMSx;
  //  if(a1>1.E5) printf("addSDD: layer=%d track #%d dedx=%e sigmaX2= %e ",
  //		    mod->GetLayer(),trackNumber,e,a1);
  rpSDD.SetSigmaX2(a1);
  a2 = kcmTomicron*fSz->GetRMS(); a2 *= a2; a2 += kRMSz*kRMSz;
  //  if(a1>1.E5) printf(" sigmaZ2= %e\n",a2);
  rpSDD.SetSigmaZ2(a2);
  rpSDD.SetProbability(1.0);

  (mod->GetITS())->AddRecPoint(rpSDD);
}
//_______________________________________________________________________
void AliITSsimulationFastPoints::AddSSD(Float_t &e,
					 AliITSmodule *mod,Int_t trackNumber){

  const Float_t kcmTomicron = 1.0e4;
  const Float_t kdEdXtoQ = 2.778e+8; // Boris Batyuna June 10 2000.
  const Float_t kRMSx = 20.0; // microns ITS TDR Table 1.3
  const Float_t kRMSz = 830.0; // microns ITS TDR Table 1.3
  Float_t a1,a2; // general float.
  AliITSRecPoint rpSSD;
  Int_t *trk = rpSSD.GetTracks();

  trk[0] = trackNumber;
  trk[1] = 0; trk[2] = 0;
  rpSSD.SetX(kcmTomicron*fSx->GetMean());
  rpSSD.SetZ(kcmTomicron*fSz->GetMean());
  rpSSD.SetdEdX(e);
  rpSSD.SetQ(kdEdXtoQ*e);
  a1 = kcmTomicron*fSx->GetRMS(); a1 *= a1; a1 += kRMSx*kRMSx;
  //  if(a1>1.E5) printf("addSSD: layer=%d track #%d dedx=%e sigmaX2= %e ",
  //		    mod->GetLayer(),trackNumber,e,a1);
  rpSSD.SetSigmaX2(a1);
  a2 = kcmTomicron*fSz->GetRMS(); a2 *= a2; a2 += kRMSz*kRMSz;
  //  if(a1>1.E5) printf(" sigmaZ2= %e RMSx=%e RMSz=%e\n",a2,fSx->GetRMS(),fSz->GetRMS());
  rpSSD.SetSigmaZ2(a2);
  rpSSD.SetProbability(1.0);

  (mod->GetITS())->AddRecPoint(rpSSD);
}
//_______________________________________________________________________
Commit	Line	Data
b0f5e3fc	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
	16	/*
	17	$Log$
	18	Revision 1.1.2.1 2000/06/11 20:16:05 barbera
	19	New: Fast simulation class for the ITS, class as part of new ITS code
	20	structure.
	21
	22	*/
	23
	24	#include <TParticle.h>
	25	#include "AliITS.h"
	26	#include "AliITSsimulationFastPoints.h"
	27	#include "AliITSstatistics.h"
	28
	29	ClassImp(AliITSsimulationFastPoints)
	30
	31	AliITSsimulationFastPoints::AliITSsimulationFastPoints()
	32	{
	33	//constructor
	34	fSx = new AliITSstatistics(2);
	35	fSz = new AliITSstatistics(2);
	36	}
	37
	38	//----------------------------------------------------------
	39	AliITSsimulationFastPoints::~AliITSsimulationFastPoints()
	40	{
	41	//destructor
	42	delete fSx;
	43	delete fSz;
	44
	45	}
	46
	47	//-------------------------------------------------------------
	48	void AliITSsimulationFastPoints::CreateFastRecPoints(AliITSmodule *mod){
	49	// Fast points simulator for all of the ITS.
	50	Int_t nhit,h,trk,ifirst;
	51	Float_t x,y,z,t,e;// local coordinate (cm) and time of flight, and dedx.
	52	Float_t x1,y1,z1;
	53	AliITShit *hit;
	54
	55	fSx->Reset(); // Start out with things clearly zeroed
	56	fSz->Reset(); // Start out with things clearly zeroed
	57	e = 0.; // Start out with things clearly zeroed
	58	Double_t weight=1.;
	59	nhit = mod->GetNhits();
	60	ifirst = 1;
	61	for(h=0;h<nhit;h++){
	62	hit = mod->GetHit(h);
	63	hit->GetPositionL(x,y,z,t);
	64	if(ifirst) {x1=x;y1=y;z1=z;}
65	e += hit->GetIonization();
66	trk = hit->GetTrack();
67	fSx->AddValue((Double_t)x,weight);
68	fSz->AddValue((Double_t)z,weight);
69	ifirst = 0;
70	if(hit->StatusExiting()\|\| // leaving volume
71	hit->StatusDisappeared()\|\| // interacted/decayed...
72	hit->StatusStop() // dropped below E cuts.
73	){ // exiting track, write out RecPoint.
74	// if(fSz->GetRMS()>1.E-1) {
75	// TParticle *part = hit->GetParticle();
76	// printf("idpart %d energy %f \n",part->GetPdgCode(),part->Energy());
77	// printf("diffx=%e diffy=%e diffz=%e\n",x-x1,y-y1,z-z1);
78	// }
79	switch (mod->GetLayer()){
80	case 1: case 2: // SPDs
81	AddSPD(e,mod,trk);
82	break;
83	case 3: case 4: // SDDs
84	AddSDD(e,mod,trk);
85	break;
86	case 5: case 6: // SSDs
87	AddSSD(e,mod,trk);
88	break;
89	} // end switch
90	fSx->Reset();
91	fSz->Reset();
92	e = 0.;
93	ifirst = 1;
94	continue;
95	}// end if
96	} // end for h
97	}
98	//_______________________________________________________________________
99	void AliITSsimulationFastPoints::AddSPD(Float_t &e,
100	AliITSmodule *mod,Int_t trackNumber){
101	const Float_t kcmTomicron = 1.0e4;
102	// const Float_t kdEdXtoQ = ;
103	const Float_t kRMSx = 12.0; // microns ITS TDR Table 1.3
104	const Float_t kRMSz = 70.0; // microns ITS TDR Table 1.3
105	Float_t a1,a2; // general float.
106	AliITSRecPoint rpSPD;
107	Int_t *trk = rpSPD.GetTracks();
108
109	trk[0] = trackNumber;
110	trk[1] = 0; trk[2] = 0;
111	rpSPD.SetX(kcmTomicron*fSx->GetMean());
112	rpSPD.SetZ(kcmTomicron*fSz->GetMean());
113	rpSPD.SetdEdX(0.0);
114	rpSPD.SetQ(1.0);
115	a1 = kcmTomicronfSx->GetRMS(); a1 = a1; a1 += kRMSx*kRMSx;
116	// if(a1>1.E5) printf("addSPD: layer=%d track #%d dedx=%e sigmaX2= %e ",
117	// mod->GetLayer(),trackNumber,e,a1);
118	rpSPD.SetSigmaX2(a1);
119	a2 = kcmTomicronfSz->GetRMS(); a2 = a2; a2 += kRMSz*kRMSz;
120	// if(a1>1.E5) printf(" sigmaZ2= %e\n",a2);
121	rpSPD.SetSigmaZ2(a2);
122	rpSPD.SetProbability(1.0);
123
124	(mod->GetITS())->AddRecPoint(rpSPD);
125	}
126	//_______________________________________________________________________
127	void AliITSsimulationFastPoints::AddSDD(Float_t &e,
128	AliITSmodule *mod,Int_t trackNumber){
129
130	const Float_t kcmTomicron = 1.0e4;
131	const Float_t kdEdXtoQ = 2.778e+8; // Boris Batyuna June 10 2000.
132	const Float_t kRMSx = 38.0; // microns ITS TDR Table 1.3
133	const Float_t kRMSz = 28.0; // microns ITS TDR Table 1.3
134	Float_t a1,a2; // general float.
135	AliITSRecPoint rpSDD;
136	Int_t *trk = rpSDD.GetTracks();
137
138	trk[0] = trackNumber;
139	trk[1] = 0; trk[2] = 0;
140	rpSDD.SetX(kcmTomicron*fSx->GetMean());
141	rpSDD.SetZ(kcmTomicron*fSz->GetMean());
142	rpSDD.SetdEdX(e);
143	rpSDD.SetQ(kdEdXtoQ*e);
144	a1 = kcmTomicronfSx->GetRMS(); a1 = a1; a1 += kRMSx*kRMSx;
145	// if(a1>1.E5) printf("addSDD: layer=%d track #%d dedx=%e sigmaX2= %e ",
146	// mod->GetLayer(),trackNumber,e,a1);
147	rpSDD.SetSigmaX2(a1);
148	a2 = kcmTomicronfSz->GetRMS(); a2 = a2; a2 += kRMSz*kRMSz;
149	// if(a1>1.E5) printf(" sigmaZ2= %e\n",a2);
150	rpSDD.SetSigmaZ2(a2);
151	rpSDD.SetProbability(1.0);
152
153	(mod->GetITS())->AddRecPoint(rpSDD);
154	}
155	//_______________________________________________________________________
156	void AliITSsimulationFastPoints::AddSSD(Float_t &e,
157	AliITSmodule *mod,Int_t trackNumber){
158
159	const Float_t kcmTomicron = 1.0e4;
160	const Float_t kdEdXtoQ = 2.778e+8; // Boris Batyuna June 10 2000.
161	const Float_t kRMSx = 20.0; // microns ITS TDR Table 1.3
162	const Float_t kRMSz = 830.0; // microns ITS TDR Table 1.3
163	Float_t a1,a2; // general float.
164	AliITSRecPoint rpSSD;
165	Int_t *trk = rpSSD.GetTracks();
166
167	trk[0] = trackNumber;
168	trk[1] = 0; trk[2] = 0;
169	rpSSD.SetX(kcmTomicron*fSx->GetMean());
170	rpSSD.SetZ(kcmTomicron*fSz->GetMean());
171	rpSSD.SetdEdX(e);
172	rpSSD.SetQ(kdEdXtoQ*e);
173	a1 = kcmTomicronfSx->GetRMS(); a1 = a1; a1 += kRMSx*kRMSx;
174	// if(a1>1.E5) printf("addSSD: layer=%d track #%d dedx=%e sigmaX2= %e ",
175	// mod->GetLayer(),trackNumber,e,a1);
176	rpSSD.SetSigmaX2(a1);
177	a2 = kcmTomicronfSz->GetRMS(); a2 = a2; a2 += kRMSz*kRMSz;
178	// if(a1>1.E5) printf(" sigmaZ2= %e RMSx=%e RMSz=%e\n",a2,fSx->GetRMS(),fSz->GetRMS());
179	rpSSD.SetSigmaZ2(a2);
180	rpSSD.SetProbability(1.0);
181
182	(mod->GetITS())->AddRecPoint(rpSSD);
183	}
184	//_______________________________________________________________________