/*
$Log$
+Revision 1.5 1999/09/29 09:24:29 fca
+Introduction of the Copyright and cvs Log
+
*/
/////////////////////////////////////////////////////////////////////
Float_t x,y,z;
Float_t pstart[3]={4.5,10.7,5.3};
- // Float_t pscin[3]={0.,2.54/2.,1.5};
Float_t ppmt[3]={0.,1.3,3.5};
Float_t pdivider[3]={0.,1.2,1.75};
Float_t pdiv2[3]={0.,1.2,1.25};
gMC->Gspos("STRT",2,"ALIC",0.,0.,-zdet,idrotm[901],"ONLY");
//START interior
- // gMC->Gsvolu("SCIN","TUBE",idtmed[2102-1],pscin,3);
gMC->Gsvolu("PMT ","TUBE",idtmed[2103-1],ppmt,3);
gMC->Gsvolu("DIVI","TUBE",idtmed[2103-1],pdivider,3);
x=6.5*TMath::Sin(is*2*3.1415/13);
y=6.5*TMath::Cos(is*2*3.1415/13);
z=-pstart[2]+ppmt[2];
- printf(" is %d x %f y %f r%f\n",is,x,y,sqrt(x*x+y*y));
- // gMC->Gspos("SCIN",is,"STRT",x,y,z,0,"ONLY");
- // z=z+pscin[2]+ppmt[2];
gMC->Gspos("PMT ",is,"STRT",x,y,z,0,"ONLY");
- z=ppmt[2]+pdivider[2];
+ z+=ppmt[2]+pdiv2[2];
printf(" is %d, z Divider %f\n",is,z);
gMC->Gspos("DIVI",is,"STRT",x,y,z,0,"ONLY");
}
x=9.3*TMath::Sin(2.*3.1415/26+(is-13)*2*3.1415/20);
y=9.3*TMath::Cos(2.*3.1315/26+(is-13)*2*3.1415/20);
z=-pstart[2]+ppmt[2];
- // gMC->Gspos("SCIN",is,"STRT",x,y,z,0,"ONLY");
- //z=z+pscin[2]+ppmt[2];
gMC->Gspos("PMT ",is,"STRT",x,y,z,0,"ONLY");
- z=ppmt[2]+pdiv2[2];
+ z+=ppmt[2]+pdiv2[2];
gMC->Gspos("DIVI",is,"STRT",x,y,z,0,"ONLY");
}
// PMT
Int_t nbuf;
// Scintillator CH
- // Float_t ascin[2]={1.01,12.01};
- // Float_t zscin[2]={1,6};
- // Float_t wscin[2]={1,1};
- // Float_t denscin=1.03;
+ Float_t ascin[2]={1.01,12.01};
+ Float_t zscin[2]={1,6};
+ Float_t wscin[2]={1,1};
+ Float_t denscin=1.03;
// PMT glass SiO2
Float_t aglass[2]={28.0855,15.9994};
Float_t zglass[2]={14.,8.};
Float_t denribber=0.8;
-// AliMC* gMC = AliMC::GetMC();
-
- // Int_t *idtmed = gAlice->Idtmed();
- // Int_t imat;
-
//*** Definition Of avaible START materials ***
AliMaterial(0, "START Steel$", 55.850,26.,7.87,1.76,999);
AliMaterial(1, "START Vacuum$", 1.e-16,1.e-16,1.e-16,1.e16,999);
zcer[1]=z;
AliMixture( 9, "Ceramic $", acer, zcer, denscer, 2, wcer);
- // AliMixture( 5, "Scintillator$",ascin,zscin,denscin,-2,wscin);
+ AliMixture( 5, "Scintillator$",ascin,zscin,denscin,-2,wscin);
AliMixture( 6, "Brass $", abrass, zbrass, denbrass, 2, wbrass);
AliMixture( 7, "Ribber $",aribber,zribber,denribber,-3,wribber);
//**
AliMedium(1, "START Air$", 2, 0, ISXFLD, SXMGMX, 10., .1, 1., .003, .003);
- // AliMedium(2102, "Scintillator$", 5, 1, ISXFLD, SXMGMX, 10., .01, 1., .003, .003);
+ AliMedium(2, "Scintillator$", 5, 1, ISXFLD, SXMGMX, 10., .01, 1., .003, .003);
AliMedium(3, "Vacuum$", 1, 0, ISXFLD, SXMGMX, 10., .01, .1, .003, .003);
AliMedium(4, "Ceramic$", 9, 0, ISXFLD, SXMGMX, 10., .01, .1, .003, .003);
AliMedium(6, "Glass$", 4, 0, ISXFLD, SXMGMX, 10., .01, .1, .003, .003);
//
//Set volumes visible
gMC->Gsatt("STRT","SEEN",0);
-//gMC->Gsatt("SCIN","SEEN",1);
gMC->Gsatt("PMT ","SEEN",1);
gMC->Gsatt("DIVI","SEEN",1);
//
Float_t ttime=gMC->TrackTime();
hits[6]=ttime*1e9;
edep=0;
- // Float_t xV=geant3->Gckine()->vert[0];
- //Float_t yV=geant3->Gckine()->vert[1];
- //Float_t zV=geant3->Gckine()->vert[2];
- //Float_t tl=gMC -> TrackLength();
- // if(hits[6]<2.4){
- // for (i=0; i<=6; i++){
- // printf(" HITS on START entr %f\n",hits[i]);}
- //}
}
if(gMC->IsTrackInside())
{
}
//---------------------------------------------------------------------
}
- //}
//void DrawSTART()
{
- geant3->Gsatt("*", "seen", -1);
- geant3->Gsatt("alic", "seen", 0);
+ gMC->Gsatt("*", "seen", -1);
+ gMC->Gsatt("alic", "seen", 0);
gROOT->Macro("ViewSTART.C");
- geant3->Gdopt("hide", "on");
- geant3->Gdopt("shad", "on");
- geant3->Gsatt("*", "fill", 7);
- geant3->SetClipBox(".");
- geant3->SetClipBox("*", 0, 1000, -1000, 1000, -1000, 1000);
- geant3->DefaultRange();
- geant3->Gdraw("alic", 40, 30, 0, 6, 9, .08, .08);
- geant3->Gdhead(1111, "START Detector");
- geant3->Gdman(13, 9, "MAN");
+ gMC->Gdopt("hide", "on");
+ gMC->Gdopt("shad", "on");
+ gMC->Gsatt("*", "fill", 7);
+ gMC->SetClipBox(".");
+ // gMC->SetClipBox("*", 0, 1000, -1000, 1000, -1000, 1000);
+ gMC->DefaultRange();
+ gMC->Gdraw("alic", 40, 30, 0, 6, 9, .08, .08);
+ gMC->Gdhead(1111, "START Detector");
+ gMC->Gdman(13, 9, "MAN");
}