[u/mrichter/AliRoot.git] / TOF / AliTOF.cxx

///////////////////////////////////////////////////////////////////////////////
//                                                                           //
//  Time Of Flight                                                           //
//  This class contains the basic functions for the Time Of Flight           //
//  detector. Functions specific to one particular geometry are              //
//  contained in the derived classes                                         //
//                                                                           //
//Begin_Html
/*
<img src="gif/AliTOFClass.gif">
*/
//End_Html
//                                                                           //
//                                                                           //
///////////////////////////////////////////////////////////////////////////////

#include "AliTOF.h"
#include <TNode.h>
#include <TTUBE.h>
#include <TBRIK.h>
#include "AliRun.h"
#include "AliMC.h"
#include "AliConst.h"
 
ClassImp(AliTOF)
 
//_____________________________________________________________________________
AliTOF::AliTOF()
{
  //
  // Default constructor
  //
  fIshunt   = 0;
}
 
//_____________________________________________________________________________
AliTOF::AliTOF(const char *name, const char *title)
       : AliDetector(name,title)
{
  //
  // AliTOF standard constructor
  // 
  fHits   = new TClonesArray("AliTOFhit",  405);
  //
  fIshunt     =  0;
  //
  SetMarkerColor(7);
  SetMarkerStyle(2);
  SetMarkerSize(0.4);
}

//_____________________________________________________________________________
void AliTOF::AddHit(Int_t track, Int_t *vol, Float_t *hits)
{
  //
  // Add a TOF hit
  //
  TClonesArray &lhits = *fHits;
  new(lhits[fNhits++]) AliTOFhit(fIshunt,track,vol,hits);
}
 
//_____________________________________________________________________________
void AliTOF::BuildGeometry()
{
  //
  // Build TOF ROOT geometry for the ALICE event viewver
  //
  TNode *Node, *Top;
  const int kColorTOF  = 27;
  //
  // Find top TNODE
  Top=gAlice->GetGeometry()->GetNode("alice");
  //
  // Define rotation matrixes
  new TRotMatrix("rot501","rot501",90,-18.94737,90,71.05263,0,0);
  new TRotMatrix("rot502","rot502",90,-37.89474,90,52.10526,0,0);
  new TRotMatrix("rot503","rot503",90,-56.84211,90,33.15789,0,0);
  new TRotMatrix("rot504","rot504",90,-75.78947,90,14.21053,0,0);
  new TRotMatrix("rot505","rot505",90,-94.73685,90,-4.736847,0,0);
  new TRotMatrix("rot506","rot506",90,-113.6842,90,-23.68421,0,0);
  new TRotMatrix("rot507","rot507",90,-132.6316,90,-42.63158,0,0);
  new TRotMatrix("rot508","rot508",90,-151.5789,90,-61.57895,0,0);
  new TRotMatrix("rot509","rot509",90,-170.5263,90,-80.52632,0,0);
  new TRotMatrix("rot510","rot510",90,-189.4737,90,-99.47369,0,0);
  new TRotMatrix("rot511","rot511",90,-208.4211,90,-118.4211,0,0);
  new TRotMatrix("rot512","rot512",90,-227.3684,90,-137.3684,0,0);
  new TRotMatrix("rot513","rot513",90,-246.3158,90,-156.3158,0,0);
  new TRotMatrix("rot514","rot514",90,-265.2632,90,-175.2632,0,0);
  new TRotMatrix("rot515","rot515",90,-284.2105,90,-194.2105,0,0);
  new TRotMatrix("rot516","rot516",90,-303.1579,90,-213.1579,0,0);
  new TRotMatrix("rot517","rot517",90,-322.1053,90,-232.1053,0,0);
  new TRotMatrix("rot518","rot518",90,-341.0526,90,-251.0526,0,0);
  new TRotMatrix("rot519","rot519",90,-360,90,-270,0,0);
  //
  // Define TOF basic volume
  new TBRIK("S_TOF1","TOF box","void",60.,6.,175.);
  //
  // Position the different copies
  Top->cd();
  Node = new TNode("FTO11","FTO11","S_TOF1",224.8018,288.8254,175,"rot502");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO12","FTO12","S_TOF1",306.4029,200.183,175,"rot503");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();  
  Node = new TNode("FTO13","FTO13","S_TOF1",354.8005,89.84769,175,"rot504");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node);
  //
  Top->cd();
  Node = new TNode("FTO14","FTO14","S_TOF1",364.7499,-30.22407,175,"rot505");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO15","FTO15","S_TOF1",335.173,-147.0205,175,"rot506");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();  
  Node = new TNode("FTO16","FTO16","S_TOF1",-335.1731,-147.0205,175,"rot513");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO17","FTO17","S_TOF1",-364.7499,-30.22407,175,"rot514");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO18","FTO18","S_TOF1",-354.8005,89.84778,175,"rot515");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO19","FTO19","S_TOF1",-306.4029,200.183,175,"rot516");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO110","FTO110","S_TOF1",-224.8019,288.8254,175,"rot517");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO111","FTO111","S_TOF1",224.8018,288.8254,-175,"rot502");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO112","FTO112","S_TOF1",306.4029,200.183,-175,"rot503");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO113","FTO113","S_TOF1",354.8005,89.84769,-175,"rot504");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO114","FTO114","S_TOF1",364.7499,-30.22407,-175,"rot505");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO115","FTO115","S_TOF1",335.173,-147.0205,-175,"rot506");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO116","FTO116","S_TOF1",-335.1731,-147.0205,-175,"rot513");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO117","FTO117","S_TOF1",-364.7499,-30.22407,-175,"rot514");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO118","FTO118","S_TOF1",-354.8005,89.84778,-175,"rot515");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO119","FTO119","S_TOF1",-306.4029,200.183,-175,"rot516");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO120","FTO120","S_TOF1",-224.8019,288.8254,-175,"rot517");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node);
  //
  // Define second TOF volume
  new TBRIK("S_TOF2","TOF box","void",60.,6.,100.);
  //
  // Position the volume
  Top->cd();
  Node = new TNode("FTO21","FTO21","S_TOF2",269.275,-247.885,250,"rot507");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO22","FTO22","S_TOF2",174.1967,-321.8874,250,"rot508");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO23","FTO23","S_TOF2",60.24158,-361.0082,250,"rot509");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO24","FTO24","S_TOF2",-60.24168,-361.0082,250,"rot510");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO25","FTO25","S_TOF2",-174.1967,-321.8874,250,"rot511");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO26","FTO26","S_TOF2",-269.275,-247.885,250,"rot512");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO27","FTO27","S_TOF2",269.275,-247.885,-250,"rot507");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO28","FTO28","S_TOF2",174.1967,-321.8874,-250,"rot508");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO29","FTO29","S_TOF2",60.24158,-361.0082,-250,"rot509");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO210","FTO210","S_TOF2",-60.24168,-361.0082,-250,"rot510");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO211","FTO211","S_TOF2",-174.1967,-321.8874,-250,"rot511");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO212","FTO212","S_TOF2",-269.275,-247.885,-250,"rot512");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  // Define third TOF volume
  new TBRIK("S_TOF3","TOF box","void",60.,6.,75.);
  //
  // Position it
  Top->cd();
  Node = new TNode("FTO31","FTO31","S_TOF3",118.84,346.1691,275,"rot501");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO32","FTO32","S_TOF3",-118.8399,346.1691,275,"rot518");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO33","FTO33","S_TOF3",0.,366.,275.,"rot519");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO34","FTO34","S_TOF3",118.84,346.1691,-275.,"rot501");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO35","FTO35","S_TOF3",-118.8399,346.1691,-275,"rot518");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
  //
  Top->cd();
  Node = new TNode("FTO36","FTO36","S_TOF3",0.,366.,-275.,"rot519");
  Node->SetLineColor(kColorTOF);
  fNodes->Add(Node); 
}

//_____________________________________________________________________________
void AliTOF::CreateGeometry()
{
  //
  // Common geometry code for version 2 and version 3 of the TOF
  //
  //Begin_Html
  /*
    <img src="gif/AliTOFv23.gif">
  */
  //End_Html
  //

  AliMC* pMC = AliMC::GetMC();
  
  const Double_t kPi=TMath::Pi();
  const Double_t kDegrad=kPi/180;
  const Double_t kRaddeg=180/kPi;
  //
  Float_t fil_rich;
  Int_t lmax;
  Float_t xtof, ytof, fil_step, phos_phi, phos_min, phos_max;
  Int_t lmax1;
  Float_t zcor1, zcor2, zcor3;
  Float_t ztof0, ztof1, ztof2, xcor2, ycor2;
  Int_t i;
  Float_t dx, dz, zl, xm, ym, phos_r;
  Int_t idrotm[101];
  Float_t phos_x;
  Float_t rp2, zm0, zm1, zm2;
  Float_t par[10], fil_min, fil_max;
  Float_t fil1;
  //
  Int_t *idtmed = gAlice->Idtmed();
  //
  phos_x = 214.6;
  phos_r = 467.;
  //phos_z = 260.;
  //rich_z = 472.5;
  xtof = 120.;
  ytof = 12.;
  ztof0 = 350.;
  ztof1 = 200.;
  ztof2 = 150.;
  //
  // frame thick along Z axis
  dz = 0.;
  //
  // frame thick along X axis
  dx = 0.;
  //
  // barrel length along Z axis
  zl = 720.;
  //
  // PHOS openings 
  fil_rich = 30.;
  phos_phi = TMath::ATan(phos_x / (phos_r * 2.));
  phos_min = (kPi - phos_phi * 4.) * kRaddeg;
  phos_max = (phos_phi * 4. + kPi) * kRaddeg;
  //
  // barrel radius in module contact point 
  par[0] = 355.;
  par[1] = 385.;
  par[2] = zl / 2.;
  pMC->Gsvolu("FBAR", "TUBE", idtmed[500], par, 3);
  //
  // --- Set module unseen ---
  pMC->Gspos("FBAR", 1, "ALIC", 0., 0., 0., 0, "ONLY");
  pMC->Gsatt("FBAR", "SEEN", 0);
  //
  // Number of TOF-block 
  lmax = 19;
  //
  // New size YTOF
  par[0] = xtof / 2.;
  par[1] = ytof / 2.;
  par[2] = ztof0 / 2.;
  pMC->Gsvolu("FTO1", "BOX ", idtmed[506], par, 3);
  pMC->Gsatt("FTO1", "SEEN", -2);
  par[2] = ztof1 / 2.;
  pMC->Gsvolu("FTO2", "BOX ", idtmed[506], par, 3);
  pMC->Gsatt("FTO2", "SEEN", -2);
  par[2] = ztof2 / 2.;
  pMC->Gsvolu("FTO3", "BOX ", idtmed[506], par, 3);
  pMC->Gsatt("FTO3", "SEEN", -2);
  //
  // Subtraction of TOF module boundaries 
  xm = xtof - dx * 2.;
  ym = ytof;
  zm0 = ztof0 - dz * 2.;
  zm1 = ztof1 - dz * 2.;
  zm2 = ztof2 - dz * 2.;
  //  
  // TOF module internal definitions
  //  
  TOFpc(xm, ym, zm0, zm1, zm2);
  //
  //rp1 = 382.;
  rp2 = 366.;
  fil_step = 360. / lmax;
  fil_min = phos_min - fil_step * .5;
  fil_max = phos_max + fil_step * .5;
  zcor1 = 175.;
  zcor2 = ztof0 - ztof1 / 2.;
  zcor3 = ztof0 - ztof2 / 2.;
  for (i = 1; i <= lmax; ++i) {
    fil1 = fil_step * i;
    //xcor1 = rp1 * TMath::Sin(fil1 * kDegrad) + dx * TMath::Cos(fil1 * kDegrad);
    //ycor1 = rp1 * TMath::Cos(fil1 * kDegrad) - dx * TMath::Sin(fil1 * kDegrad);
    xcor2 = rp2 * TMath::Sin(fil1 * kDegrad);
    ycor2 = rp2 * TMath::Cos(fil1 * kDegrad);
    lmax1 = i + lmax;
    AliMatrix(idrotm[i], 90., -fil1, 90., 90. - fil1, 0., 0.);
    if (fil1 >= fil_min  && fil1 <= fil_max) {
      pMC->Gspos("FTO2", i, "FBAR", xcor2, ycor2, zcor2, idrotm[i], "ONLY");
      pMC->Gspos("FTO2", lmax1, "FBAR", xcor2, ycor2, -zcor2, idrotm[i], "ONLY");
    } else if (fil1 <= fil_rich || fil1 >= 360. - fil_rich) {
      par[2] = ztof2 / 2.;
      pMC->Gspos("FTO3", i, "FBAR", xcor2, ycor2, zcor3, idrotm[i], "ONLY");
      pMC->Gspos("FTO3", lmax1, "FBAR", xcor2, ycor2, -zcor3, idrotm[i], "ONLY");
    } else {
      par[2] = ztof0 / 2.;
      pMC->Gspos("FTO1", i, "FBAR", xcor2, ycor2, zcor1, idrotm[i], "ONLY");
      pMC->Gspos("FTO1", lmax1, "FBAR", xcor2, ycor2, -zcor1, idrotm[i], "ONLY");
    }
  }
}

//_____________________________________________________________________________
void AliTOF::DrawDetector()
{
  //
  // Draw a shaded view of the common part of the TOF geometry
  // for versions 2 and 3
  //

  AliMC* pMC = AliMC::GetMC();
  
  // Set everything unseen
  pMC->Gsatt("*", "seen", -1);
  // 
  // Set ALIC mother transparent
  pMC->Gsatt("ALIC","SEEN",0);
  //
  // Set the volumes visible
  pMC->Gsatt("FBAR","SEEN",0);
  pMC->Gsatt("FTO1","SEEN",1);
  pMC->Gsatt("FTO2","SEEN",1);
  pMC->Gsatt("FTO3","SEEN",1);
  //
  pMC->Gdopt("hide", "on");
  pMC->Gdopt("shad", "on");
  pMC->Gsatt("*", "fill", 7);
  pMC->SetClipBox(".");
  pMC->SetClipBox("*", 0, 1000, -1000, 1000, -1000, 1000);
  pMC->DefaultRange();
  pMC->Gdraw("alic", 40, 30, 0, 12, 9.5, .02, .02);
  pMC->Gdhead(1111, "Time Of Flight");
  pMC->Gdman(18, 4, "MAN");
  pMC->Gdopt("hide","off");
}

//_____________________________________________________________________________
void AliTOF::CreateMaterials()
{
  //
  // Defines TOF materials for all versions
  // Authors :   Maxim Martemianov, Boris Zagreev (ITEP)   18/09/98 
  //
  Int_t   ISXFLD = gAlice->Field()->Integ();
  Float_t SXMGMX = gAlice->Field()->Max();
  //
  Float_t ac[2]   = { 12.,16. };
  Float_t zc[2]   = { 6.,8. };
  Float_t wc[2]   = { 1.,2. };
  Float_t ag10[4] = { 12.,1.,16.,28. };
  Float_t zg10[4] = { 6.,1.,8.,14. };
  Float_t wmatg10[4] = { .259,.288,.248,.205 };
  Float_t adme[5] = { 12.,1.,16.,19.,79. };
  Float_t zdme[5] = { 6.,1.,8.,9.,35. };
  Float_t wmatdme[5] = { .4056,.0961,.2562,.1014,.1407 };
  Float_t aal[2] = { 27.,16. };
  Float_t zal[2] = { 13.,8. };
  Float_t wmatal[2] = { 2.,3. };
  //
  Int_t nlmatdme;
  Float_t epsil, stmin, dc, densg10, densal, deemax, stemax;
  Float_t densdme;
  Int_t nlmatg10, nlmatal;
  //
  // --- Vacuum
  AliMaterial(0, "Vacuum$", 1e-16, 1e-16, 1e-16, 1e16, 1e16);
  // --- Air
  AliMaterial(1, "Air$",14.61,7.3,0.001205,30423.24,67500.);
  // --- CO2 
  dc = .001977;
  AliMixture(7, "CO2$", ac, zc, dc, -2, wc);
  // --- G10 
  densg10  = 1.7;
  nlmatg10 = -4;
  AliMixture(5, "G10$", ag10, zg10, densg10, nlmatg10, wmatg10);
  // --- DME 
  densdme  = .00205;
  nlmatdme = 5;
  AliMixture(6, "DME ", adme, zdme, densdme, nlmatdme, wmatdme);
  // ---- ALUMINA (AL203) 
  densal  = 2.3;
  nlmatal = -2;
  AliMixture(8, "ALUMINA$", aal, zal, densal, nlmatal, wmatal);
  //       Previous data 
  //       EPSIL  =  0.1   ! Tracking precision, 
  //       STEMAX = 0.1      ! Maximum displacement for multiple scattering
  //       DEEMAX = 0.1    ! Maximum fractional energy loss, DLS 
  //       STMIN  = 0.1 
  //       New data from 
  epsil  = .001;
  stemax = -1.;
  deemax = -.3;
  stmin  = -.8;
  AliMedium(500, "Vacuum  $", 0, 0, ISXFLD, SXMGMX, 10., stemax, deemax, epsil, stmin);
  AliMedium(501, "Air $", 0, 0, ISXFLD, SXMGMX, 10., stemax, deemax, epsil, stmin);
  AliMedium(505, "G10$", 5, 0, ISXFLD, SXMGMX, 10., stemax, deemax, epsil, stmin);
  AliMedium(506, "DME$", 6, 0, ISXFLD, SXMGMX, 10., stemax, deemax, epsil, stmin);
  AliMedium(507, "CO2$", 7, 0, ISXFLD, SXMGMX, 10., -.01, -.1, .01, -.01);
  AliMedium(508, "ALUMINA$", 8, 0, ISXFLD, SXMGMX, 10., stemax, deemax, epsil, stmin);
  AliMedium(510, "DME$", 6, 1, ISXFLD, SXMGMX, 10., stemax, deemax, epsil, stmin);
}

//_____________________________________________________________________________
Int_t AliTOF::DistancetoPrimitive(Int_t , Int_t )
{
  //
  // Returns distance from mouse pointer to detector, default version
  //
  return 9999;
}
 
//_____________________________________________________________________________
void AliTOF::Init()
{
  //
  // Initialise TOF detector after it has been built
  //
  Int_t i;
  AliMC *pMC=AliMC::GetMC();
  //
  printf("\n");
  for(i=0;i<35;i++) printf("*");
  printf(" TOF_INIT ");
  for(i=0;i<35;i++) printf("*");
  printf("\n");
  //
  // Set id of TOF sensitive volume
  fIdSens=pMC->VolId("FPG2");
  //
  for(i=0;i<80;i++) printf("*");
  printf("\n");
}

 
ClassImp(AliTOFhit)
 
//___________________________________________
AliTOFhit::AliTOFhit(Int_t shunt, Int_t track, Int_t *vol, Float_t *hits):
  AliHit(shunt, track)
{
  //
  // Store a TOF hit
  //
  Int_t i;
  for (i=0;i<3;i++) fVolume[i] = vol[i];
  //
  // Position
  fX=hits[0];
  fY=hits[1];
  fZ=hits[2];
  //
  // Momentum
  fPx=hits[3];
  fPy=hits[4];
  fPz=hits[5];
  fPmom=hits[6];
  //
  // Time Of Flight
  fTof=hits[7];
}
Commit	Line	Data
fe4da5cc	1	///////////////////////////////////////////////////////////////////////////////
	2	// //
	3	// Time Of Flight //
	4	// This class contains the basic functions for the Time Of Flight //
	5	// detector. Functions specific to one particular geometry are //
	6	// contained in the derived classes //
	7	// //
	8	//Begin_Html
	9	/*
	10	<img src="gif/AliTOFClass.gif">
	11	*/
	12	//End_Html
	13	// //
	14	// //
	15	///////////////////////////////////////////////////////////////////////////////
	16
	17	#include "AliTOF.h"
	18	#include <TNode.h>
	19	#include <TTUBE.h>
	20	#include <TBRIK.h>
	21	#include "AliRun.h"
	22	#include "AliMC.h"
	23	#include "AliConst.h"
	24
	25	ClassImp(AliTOF)
	26
	27	//_____________________________________________________________________________
	28	AliTOF::AliTOF()
	29	{
	30	//
	31	// Default constructor
	32	//
	33	fIshunt = 0;
	34	}
	35
	36	//_____________________________________________________________________________
	37	AliTOF::AliTOF(const char name, const char title)
	38	: AliDetector(name,title)
	39	{
	40	//
	41	// AliTOF standard constructor
	42	//
	43	fHits = new TClonesArray("AliTOFhit", 405);
	44	//
	45	fIshunt = 0;
	46	//
	47	SetMarkerColor(7);
	48	SetMarkerStyle(2);
	49	SetMarkerSize(0.4);
	50	}
	51
	52	//_____________________________________________________________________________
	53	void AliTOF::AddHit(Int_t track, Int_t vol, Float_t hits)
	54	{
	55	//
	56	// Add a TOF hit
	57	//
	58	TClonesArray &lhits = *fHits;
	59	new(lhits[fNhits++]) AliTOFhit(fIshunt,track,vol,hits);
	60	}
	61
	62	//_____________________________________________________________________________
	63	void AliTOF::BuildGeometry()
	64	{
65	//
66	// Build TOF ROOT geometry for the ALICE event viewver
67	//
68	TNode Node, Top;
69	const int kColorTOF = 27;
70	//
71	// Find top TNODE
72	Top=gAlice->GetGeometry()->GetNode("alice");
73	//
74	// Define rotation matrixes
75	new TRotMatrix("rot501","rot501",90,-18.94737,90,71.05263,0,0);
76	new TRotMatrix("rot502","rot502",90,-37.89474,90,52.10526,0,0);
77	new TRotMatrix("rot503","rot503",90,-56.84211,90,33.15789,0,0);
78	new TRotMatrix("rot504","rot504",90,-75.78947,90,14.21053,0,0);
79	new TRotMatrix("rot505","rot505",90,-94.73685,90,-4.736847,0,0);
80	new TRotMatrix("rot506","rot506",90,-113.6842,90,-23.68421,0,0);
81	new TRotMatrix("rot507","rot507",90,-132.6316,90,-42.63158,0,0);
82	new TRotMatrix("rot508","rot508",90,-151.5789,90,-61.57895,0,0);
83	new TRotMatrix("rot509","rot509",90,-170.5263,90,-80.52632,0,0);
84	new TRotMatrix("rot510","rot510",90,-189.4737,90,-99.47369,0,0);
85	new TRotMatrix("rot511","rot511",90,-208.4211,90,-118.4211,0,0);
86	new TRotMatrix("rot512","rot512",90,-227.3684,90,-137.3684,0,0);
87	new TRotMatrix("rot513","rot513",90,-246.3158,90,-156.3158,0,0);
88	new TRotMatrix("rot514","rot514",90,-265.2632,90,-175.2632,0,0);
89	new TRotMatrix("rot515","rot515",90,-284.2105,90,-194.2105,0,0);
90	new TRotMatrix("rot516","rot516",90,-303.1579,90,-213.1579,0,0);
91	new TRotMatrix("rot517","rot517",90,-322.1053,90,-232.1053,0,0);
92	new TRotMatrix("rot518","rot518",90,-341.0526,90,-251.0526,0,0);
93	new TRotMatrix("rot519","rot519",90,-360,90,-270,0,0);
94	//
95	// Define TOF basic volume
96	new TBRIK("S_TOF1","TOF box","void",60.,6.,175.);
97	//
98	// Position the different copies
99	Top->cd();
100	Node = new TNode("FTO11","FTO11","S_TOF1",224.8018,288.8254,175,"rot502");
101	Node->SetLineColor(kColorTOF);
102	fNodes->Add(Node);
103	//
104	Top->cd();
105	Node = new TNode("FTO12","FTO12","S_TOF1",306.4029,200.183,175,"rot503");
106	Node->SetLineColor(kColorTOF);
107	fNodes->Add(Node);
108	//
109	Top->cd();
110	Node = new TNode("FTO13","FTO13","S_TOF1",354.8005,89.84769,175,"rot504");
111	Node->SetLineColor(kColorTOF);
112	fNodes->Add(Node);
113	//
114	Top->cd();
115	Node = new TNode("FTO14","FTO14","S_TOF1",364.7499,-30.22407,175,"rot505");
116	Node->SetLineColor(kColorTOF);
117	fNodes->Add(Node);
118	//
119	Top->cd();
120	Node = new TNode("FTO15","FTO15","S_TOF1",335.173,-147.0205,175,"rot506");
121	Node->SetLineColor(kColorTOF);
122	fNodes->Add(Node);
123	//
124	Top->cd();
125	Node = new TNode("FTO16","FTO16","S_TOF1",-335.1731,-147.0205,175,"rot513");
126	Node->SetLineColor(kColorTOF);
127	fNodes->Add(Node);
128	//
129	Top->cd();
130	Node = new TNode("FTO17","FTO17","S_TOF1",-364.7499,-30.22407,175,"rot514");
131	Node->SetLineColor(kColorTOF);
132	fNodes->Add(Node);
133	//
134	Top->cd();
135	Node = new TNode("FTO18","FTO18","S_TOF1",-354.8005,89.84778,175,"rot515");
136	Node->SetLineColor(kColorTOF);
137	fNodes->Add(Node);
138	//
139	Top->cd();
140	Node = new TNode("FTO19","FTO19","S_TOF1",-306.4029,200.183,175,"rot516");
141	Node->SetLineColor(kColorTOF);
142	fNodes->Add(Node);
143	//
144	Top->cd();
145	Node = new TNode("FTO110","FTO110","S_TOF1",-224.8019,288.8254,175,"rot517");
146	Node->SetLineColor(kColorTOF);
147	fNodes->Add(Node);
148	//
149	Top->cd();
150	Node = new TNode("FTO111","FTO111","S_TOF1",224.8018,288.8254,-175,"rot502");
151	Node->SetLineColor(kColorTOF);
152	fNodes->Add(Node);
153	//
154	Top->cd();
155	Node = new TNode("FTO112","FTO112","S_TOF1",306.4029,200.183,-175,"rot503");
156	Node->SetLineColor(kColorTOF);
157	fNodes->Add(Node);
158	//
159	Top->cd();
160	Node = new TNode("FTO113","FTO113","S_TOF1",354.8005,89.84769,-175,"rot504");
161	Node->SetLineColor(kColorTOF);
162	fNodes->Add(Node);
163	//
164	Top->cd();
165	Node = new TNode("FTO114","FTO114","S_TOF1",364.7499,-30.22407,-175,"rot505");
166	Node->SetLineColor(kColorTOF);
167	fNodes->Add(Node);
168	//
169	Top->cd();
170	Node = new TNode("FTO115","FTO115","S_TOF1",335.173,-147.0205,-175,"rot506");
171	Node->SetLineColor(kColorTOF);
172	fNodes->Add(Node);
173	//
174	Top->cd();
175	Node = new TNode("FTO116","FTO116","S_TOF1",-335.1731,-147.0205,-175,"rot513");
176	Node->SetLineColor(kColorTOF);
177	fNodes->Add(Node);
178	//
179	Top->cd();
180	Node = new TNode("FTO117","FTO117","S_TOF1",-364.7499,-30.22407,-175,"rot514");
181	Node->SetLineColor(kColorTOF);
182	fNodes->Add(Node);
183	//
184	Top->cd();
185	Node = new TNode("FTO118","FTO118","S_TOF1",-354.8005,89.84778,-175,"rot515");
186	Node->SetLineColor(kColorTOF);
187	fNodes->Add(Node);
188	//
189	Top->cd();
190	Node = new TNode("FTO119","FTO119","S_TOF1",-306.4029,200.183,-175,"rot516");
191	Node->SetLineColor(kColorTOF);
192	fNodes->Add(Node);
193	//
194	Top->cd();
195	Node = new TNode("FTO120","FTO120","S_TOF1",-224.8019,288.8254,-175,"rot517");
196	Node->SetLineColor(kColorTOF);
197	fNodes->Add(Node);
198	//
199	// Define second TOF volume
200	new TBRIK("S_TOF2","TOF box","void",60.,6.,100.);
201	//
202	// Position the volume
203	Top->cd();
204	Node = new TNode("FTO21","FTO21","S_TOF2",269.275,-247.885,250,"rot507");
205	Node->SetLineColor(kColorTOF);
206	fNodes->Add(Node);
207	//
208	Top->cd();
209	Node = new TNode("FTO22","FTO22","S_TOF2",174.1967,-321.8874,250,"rot508");
210	Node->SetLineColor(kColorTOF);
211	fNodes->Add(Node);
212	//
213	Top->cd();
214	Node = new TNode("FTO23","FTO23","S_TOF2",60.24158,-361.0082,250,"rot509");
215	Node->SetLineColor(kColorTOF);
216	fNodes->Add(Node);
217	//
218	Top->cd();
219	Node = new TNode("FTO24","FTO24","S_TOF2",-60.24168,-361.0082,250,"rot510");
220	Node->SetLineColor(kColorTOF);
221	fNodes->Add(Node);
222	//
223	Top->cd();
224	Node = new TNode("FTO25","FTO25","S_TOF2",-174.1967,-321.8874,250,"rot511");
225	Node->SetLineColor(kColorTOF);
226	fNodes->Add(Node);
227	//
228	Top->cd();
229	Node = new TNode("FTO26","FTO26","S_TOF2",-269.275,-247.885,250,"rot512");
230	Node->SetLineColor(kColorTOF);
231	fNodes->Add(Node);
232	//
233	Top->cd();
234	Node = new TNode("FTO27","FTO27","S_TOF2",269.275,-247.885,-250,"rot507");
235	Node->SetLineColor(kColorTOF);
236	fNodes->Add(Node);
237	//
238	Top->cd();
239	Node = new TNode("FTO28","FTO28","S_TOF2",174.1967,-321.8874,-250,"rot508");
240	Node->SetLineColor(kColorTOF);
241	fNodes->Add(Node);
242	//
243	Top->cd();
244	Node = new TNode("FTO29","FTO29","S_TOF2",60.24158,-361.0082,-250,"rot509");
245	Node->SetLineColor(kColorTOF);
246	fNodes->Add(Node);
247	//
248	Top->cd();
249	Node = new TNode("FTO210","FTO210","S_TOF2",-60.24168,-361.0082,-250,"rot510");
250	Node->SetLineColor(kColorTOF);
251	fNodes->Add(Node);
252	//
253	Top->cd();
254	Node = new TNode("FTO211","FTO211","S_TOF2",-174.1967,-321.8874,-250,"rot511");
255	Node->SetLineColor(kColorTOF);
256	fNodes->Add(Node);
257	//
258	Top->cd();
259	Node = new TNode("FTO212","FTO212","S_TOF2",-269.275,-247.885,-250,"rot512");
260	Node->SetLineColor(kColorTOF);
261	fNodes->Add(Node);
262	//
263	// Define third TOF volume
264	new TBRIK("S_TOF3","TOF box","void",60.,6.,75.);
265	//
266	// Position it
267	Top->cd();
268	Node = new TNode("FTO31","FTO31","S_TOF3",118.84,346.1691,275,"rot501");
269	Node->SetLineColor(kColorTOF);
270	fNodes->Add(Node);
271	//
272	Top->cd();
273	Node = new TNode("FTO32","FTO32","S_TOF3",-118.8399,346.1691,275,"rot518");
274	Node->SetLineColor(kColorTOF);
275	fNodes->Add(Node);
276	//
277	Top->cd();
278	Node = new TNode("FTO33","FTO33","S_TOF3",0.,366.,275.,"rot519");
279	Node->SetLineColor(kColorTOF);
280	fNodes->Add(Node);
281	//
282	Top->cd();
283	Node = new TNode("FTO34","FTO34","S_TOF3",118.84,346.1691,-275.,"rot501");
284	Node->SetLineColor(kColorTOF);
285	fNodes->Add(Node);
286	//
287	Top->cd();
288	Node = new TNode("FTO35","FTO35","S_TOF3",-118.8399,346.1691,-275,"rot518");
289	Node->SetLineColor(kColorTOF);
290	fNodes->Add(Node);
291	//
292	Top->cd();
293	Node = new TNode("FTO36","FTO36","S_TOF3",0.,366.,-275.,"rot519");
294	Node->SetLineColor(kColorTOF);
295	fNodes->Add(Node);
296	}
297
298	//_____________________________________________________________________________
299	void AliTOF::CreateGeometry()
300	{
301	//
302	// Common geometry code for version 2 and version 3 of the TOF
303	//
304	//Begin_Html
305	/*
306	<img src="gif/AliTOFv23.gif">
307	*/
308	//End_Html
309	//
310
311	AliMC* pMC = AliMC::GetMC();
312
313	const Double_t kPi=TMath::Pi();
314	const Double_t kDegrad=kPi/180;
315	const Double_t kRaddeg=180/kPi;
316	//
317	Float_t fil_rich;
318	Int_t lmax;
319	Float_t xtof, ytof, fil_step, phos_phi, phos_min, phos_max;
320	Int_t lmax1;
321	Float_t zcor1, zcor2, zcor3;
322	Float_t ztof0, ztof1, ztof2, xcor2, ycor2;
323	Int_t i;
324	Float_t dx, dz, zl, xm, ym, phos_r;
325	Int_t idrotm[101];
326	Float_t phos_x;
327	Float_t rp2, zm0, zm1, zm2;
328	Float_t par[10], fil_min, fil_max;
329	Float_t fil1;
330	//
331	Int_t *idtmed = gAlice->Idtmed();
332	//
333	phos_x = 214.6;
334	phos_r = 467.;
335	//phos_z = 260.;
336	//rich_z = 472.5;
337	xtof = 120.;
338	ytof = 12.;
339	ztof0 = 350.;
340	ztof1 = 200.;
341	ztof2 = 150.;
342	//
343	// frame thick along Z axis
344	dz = 0.;
345	//
346	// frame thick along X axis
347	dx = 0.;
348	//
349	// barrel length along Z axis
350	zl = 720.;
351	//
352	// PHOS openings
353	fil_rich = 30.;
354	phos_phi = TMath::ATan(phos_x / (phos_r * 2.));
355	phos_min = (kPi - phos_phi * 4.) * kRaddeg;
356	phos_max = (phos_phi * 4. + kPi) * kRaddeg;
357	//
358	// barrel radius in module contact point
359	par[0] = 355.;
360	par[1] = 385.;
361	par[2] = zl / 2.;
362	pMC->Gsvolu("FBAR", "TUBE", idtmed[500], par, 3);
363	//
364	// --- Set module unseen ---
365	pMC->Gspos("FBAR", 1, "ALIC", 0., 0., 0., 0, "ONLY");
366	pMC->Gsatt("FBAR", "SEEN", 0);
367	//
368	// Number of TOF-block
369	lmax = 19;
370	//
371	// New size YTOF
372	par[0] = xtof / 2.;
373	par[1] = ytof / 2.;
374	par[2] = ztof0 / 2.;
375	pMC->Gsvolu("FTO1", "BOX ", idtmed[506], par, 3);
376	pMC->Gsatt("FTO1", "SEEN", -2);
377	par[2] = ztof1 / 2.;
378	pMC->Gsvolu("FTO2", "BOX ", idtmed[506], par, 3);
379	pMC->Gsatt("FTO2", "SEEN", -2);
380	par[2] = ztof2 / 2.;
381	pMC->Gsvolu("FTO3", "BOX ", idtmed[506], par, 3);
382	pMC->Gsatt("FTO3", "SEEN", -2);
383	//
384	// Subtraction of TOF module boundaries
385	xm = xtof - dx * 2.;
386	ym = ytof;
387	zm0 = ztof0 - dz * 2.;
388	zm1 = ztof1 - dz * 2.;
389	zm2 = ztof2 - dz * 2.;
390	//
391	// TOF module internal definitions
392	//
393	TOFpc(xm, ym, zm0, zm1, zm2);
394	//
395	//rp1 = 382.;
396	rp2 = 366.;
397	fil_step = 360. / lmax;
398	fil_min = phos_min - fil_step * .5;
399	fil_max = phos_max + fil_step * .5;
400	zcor1 = 175.;
401	zcor2 = ztof0 - ztof1 / 2.;
402	zcor3 = ztof0 - ztof2 / 2.;
403	for (i = 1; i <= lmax; ++i) {
404	fil1 = fil_step * i;
405	//xcor1 = rp1 * TMath::Sin(fil1 * kDegrad) + dx * TMath::Cos(fil1 * kDegrad);
406	//ycor1 = rp1 * TMath::Cos(fil1 * kDegrad) - dx * TMath::Sin(fil1 * kDegrad);
407	xcor2 = rp2 * TMath::Sin(fil1 * kDegrad);
408	ycor2 = rp2 * TMath::Cos(fil1 * kDegrad);
409	lmax1 = i + lmax;
410	AliMatrix(idrotm[i], 90., -fil1, 90., 90. - fil1, 0., 0.);
411	if (fil1 >= fil_min && fil1 <= fil_max) {
412	pMC->Gspos("FTO2", i, "FBAR", xcor2, ycor2, zcor2, idrotm[i], "ONLY");
413	pMC->Gspos("FTO2", lmax1, "FBAR", xcor2, ycor2, -zcor2, idrotm[i], "ONLY");
414	} else if (fil1 <= fil_rich \|\| fil1 >= 360. - fil_rich) {
415	par[2] = ztof2 / 2.;
416	pMC->Gspos("FTO3", i, "FBAR", xcor2, ycor2, zcor3, idrotm[i], "ONLY");
417	pMC->Gspos("FTO3", lmax1, "FBAR", xcor2, ycor2, -zcor3, idrotm[i], "ONLY");
418	} else {
419	par[2] = ztof0 / 2.;
420	pMC->Gspos("FTO1", i, "FBAR", xcor2, ycor2, zcor1, idrotm[i], "ONLY");
421	pMC->Gspos("FTO1", lmax1, "FBAR", xcor2, ycor2, -zcor1, idrotm[i], "ONLY");
422	}
423	}
424	}
425
426	//_____________________________________________________________________________
427	void AliTOF::DrawDetector()
428	{
429	//
430	// Draw a shaded view of the common part of the TOF geometry
431	// for versions 2 and 3
432	//
433
434	AliMC* pMC = AliMC::GetMC();
435
436	// Set everything unseen
437	pMC->Gsatt("*", "seen", -1);
438	//
439	// Set ALIC mother transparent
440	pMC->Gsatt("ALIC","SEEN",0);
441	//
442	// Set the volumes visible
443	pMC->Gsatt("FBAR","SEEN",0);
444	pMC->Gsatt("FTO1","SEEN",1);
445	pMC->Gsatt("FTO2","SEEN",1);
446	pMC->Gsatt("FTO3","SEEN",1);
447	//
448	pMC->Gdopt("hide", "on");
449	pMC->Gdopt("shad", "on");
450	pMC->Gsatt("*", "fill", 7);
451	pMC->SetClipBox(".");
452	pMC->SetClipBox("*", 0, 1000, -1000, 1000, -1000, 1000);
453	pMC->DefaultRange();
454	pMC->Gdraw("alic", 40, 30, 0, 12, 9.5, .02, .02);
455	pMC->Gdhead(1111, "Time Of Flight");
456	pMC->Gdman(18, 4, "MAN");
457	pMC->Gdopt("hide","off");
458	}
459
460	//_____________________________________________________________________________
461	void AliTOF::CreateMaterials()
462	{
463	//
464	// Defines TOF materials for all versions
465	// Authors : Maxim Martemianov, Boris Zagreev (ITEP) 18/09/98
466	//
467	Int_t ISXFLD = gAlice->Field()->Integ();
468	Float_t SXMGMX = gAlice->Field()->Max();
469	//
470	Float_t ac[2] = { 12.,16. };
471	Float_t zc[2] = { 6.,8. };
472	Float_t wc[2] = { 1.,2. };
473	Float_t ag10[4] = { 12.,1.,16.,28. };
474	Float_t zg10[4] = { 6.,1.,8.,14. };
475	Float_t wmatg10[4] = { .259,.288,.248,.205 };
476	Float_t adme[5] = { 12.,1.,16.,19.,79. };
477	Float_t zdme[5] = { 6.,1.,8.,9.,35. };
478	Float_t wmatdme[5] = { .4056,.0961,.2562,.1014,.1407 };
479	Float_t aal[2] = { 27.,16. };
480	Float_t zal[2] = { 13.,8. };
481	Float_t wmatal[2] = { 2.,3. };
482	//
483	Int_t nlmatdme;
484	Float_t epsil, stmin, dc, densg10, densal, deemax, stemax;
485	Float_t densdme;
486	Int_t nlmatg10, nlmatal;
487	//
488	// --- Vacuum
489	AliMaterial(0, "Vacuum$", 1e-16, 1e-16, 1e-16, 1e16, 1e16);
490	// --- Air
491	AliMaterial(1, "Air$",14.61,7.3,0.001205,30423.24,67500.);
492	// --- CO2
493	dc = .001977;
494	AliMixture(7, "CO2$", ac, zc, dc, -2, wc);
495	// --- G10
496	densg10 = 1.7;
497	nlmatg10 = -4;
498	AliMixture(5, "G10$", ag10, zg10, densg10, nlmatg10, wmatg10);
499	// --- DME
500	densdme = .00205;
501	nlmatdme = 5;
502	AliMixture(6, "DME ", adme, zdme, densdme, nlmatdme, wmatdme);
503	// ---- ALUMINA (AL203)
504	densal = 2.3;
505	nlmatal = -2;
506	AliMixture(8, "ALUMINA$", aal, zal, densal, nlmatal, wmatal);
507	// Previous data
508	// EPSIL = 0.1 ! Tracking precision,
509	// STEMAX = 0.1 ! Maximum displacement for multiple scattering
510	// DEEMAX = 0.1 ! Maximum fractional energy loss, DLS
511	// STMIN = 0.1
512	// New data from
513	epsil = .001;
514	stemax = -1.;
515	deemax = -.3;
516	stmin = -.8;
517	AliMedium(500, "Vacuum $", 0, 0, ISXFLD, SXMGMX, 10., stemax, deemax, epsil, stmin);
518	AliMedium(501, "Air $", 0, 0, ISXFLD, SXMGMX, 10., stemax, deemax, epsil, stmin);
519	AliMedium(505, "G10$", 5, 0, ISXFLD, SXMGMX, 10., stemax, deemax, epsil, stmin);
520	AliMedium(506, "DME$", 6, 0, ISXFLD, SXMGMX, 10., stemax, deemax, epsil, stmin);
521	AliMedium(507, "CO2$", 7, 0, ISXFLD, SXMGMX, 10., -.01, -.1, .01, -.01);
522	AliMedium(508, "ALUMINA$", 8, 0, ISXFLD, SXMGMX, 10., stemax, deemax, epsil, stmin);
523	AliMedium(510, "DME$", 6, 1, ISXFLD, SXMGMX, 10., stemax, deemax, epsil, stmin);
524	}
525
526	//_____________________________________________________________________________
527	Int_t AliTOF::DistancetoPrimitive(Int_t , Int_t )
528	{
529	//
530	// Returns distance from mouse pointer to detector, default version
531	//
532	return 9999;
533	}
534
535	//_____________________________________________________________________________
536	void AliTOF::Init()
537	{
538	//
539	// Initialise TOF detector after it has been built
540	//
541	Int_t i;
542	AliMC *pMC=AliMC::GetMC();
543	//
544	printf("\n");
545	for(i=0;i<35;i++) printf("*");
546	printf(" TOF_INIT ");
547	for(i=0;i<35;i++) printf("*");
548	printf("\n");
549	//
550	// Set id of TOF sensitive volume
551	fIdSens=pMC->VolId("FPG2");
552	//
553	for(i=0;i<80;i++) printf("*");
554	printf("\n");
555	}
556
557
558	ClassImp(AliTOFhit)
559
560	//___________________________________________
561	AliTOFhit::AliTOFhit(Int_t shunt, Int_t track, Int_t vol, Float_t hits):
562	AliHit(shunt, track)
563	{
564	//
565	// Store a TOF hit
566	//
567	Int_t i;
568	for (i=0;i<3;i++) fVolume[i] = vol[i];
569	//
570	// Position
571	fX=hits[0];
572	fY=hits[1];
573	fZ=hits[2];
574	//
575	// Momentum
576	fPx=hits[3];
577	fPy=hits[4];
578	fPz=hits[5];
579	fPmom=hits[6];
580	//
581	// Time Of Flight
582	fTof=hits[7];
583	}
584
585