#include "TGeoPara.h"
#include "TGeoPhysicalNode.h"
#include "TGeoHalfSpace.h"
+#include "TTreeStream.h"
ClassImp(AliTPCv2)
// here I define a volume TPC
// retrive the medium name with "TPC_" as a leading string
//
- TGeoPcon *tpc = new TGeoPcon(0.,360.,20); //20 sections
+ TGeoPcon *tpc = new TGeoPcon(0.,360.,24); //24 sections
//
tpc->DefineSection(0,-291.,77.,278.);
tpc->DefineSection(1,-270,77.,278.);
tpc->DefineSection(9,-73.3,60.1,278.);
//
tpc->DefineSection(10,-73.3,56.9,278.);
- tpc->DefineSection(11,73.3,56.9,278.);
+ tpc->DefineSection(11,-68.5,56.9,278.);
//
- tpc->DefineSection(12,73.3,60.1,278.);
- tpc->DefineSection(13,74.0,60.1,278.);
+ tpc->DefineSection(12,-68.5,60.,278.);
+ tpc->DefineSection(13,-65.,60.,278.);
//
- tpc->DefineSection(14,74.0,60.8,278.);
- tpc->DefineSection(15,253.6,65.5,278.);
+ tpc->DefineSection(14,-65.,56.9,278.);
+ tpc->DefineSection(15,73.3,56.9,278.);
//
- tpc->DefineSection(16,253.6,65.6,278.);
- tpc->DefineSection(17,259.6,65.6,278.);
+ tpc->DefineSection(16,73.3,60.1,278.);
+ tpc->DefineSection(17,74.0,60.1,278.);
//
- tpc->DefineSection(18,259.6,70.0,278.);
- tpc->DefineSection(19,291.,77.,278.);
+ tpc->DefineSection(18,74.0,60.8,278.);
+ tpc->DefineSection(19,253.6,65.5,278.);
+ //
+ tpc->DefineSection(20,253.6,65.6,278.);
+ tpc->DefineSection(21,259.6,65.6,278.);
+ //
+ tpc->DefineSection(22,259.6,70.0,278.);
+ tpc->DefineSection(23,291.,77.,278.);
//
TGeoMedium *m1 = gGeoManager->GetMedium("TPC_Air");
TGeoVolume *v1 = new TGeoVolume("TPC_M",tpc,m1);
//
// segment outermost
//
- TGeoTubeSeg *t2 = new TGeoTubeSeg(76.6774,78.845,74.175,350.,109.6); // tedlar 38 microns
- TGeoTubeSeg *t3 = new TGeoTubeSeg(76.6812,78.8412,74.175,350.,109.6); // prepreg2 500 microns
- TGeoTubeSeg *t4 = new TGeoTubeSeg(76.7312,78.7912,74.175,350.,109.6); // prepreg3 300 microns
- TGeoTubeSeg *t5 = new TGeoTubeSeg(76.7612,78.7612,74.175,350.,109.6); // nomex 2 cm
- TGeoTubeSeg *tepox1 = new TGeoTubeSeg(76.6774,78.845,74.175,109.6,110.);//epoxy
+ TGeoTubeSeg *t2 = new TGeoTubeSeg(76.6774,78.845,74.175,350.,109.4); // tedlar 38 microns
+ TGeoTubeSeg *t3 = new TGeoTubeSeg(76.6812,78.8412,74.175,350.,109.4); // prepreg2 500 microns
+ TGeoTubeSeg *t4 = new TGeoTubeSeg(76.7312,78.7912,74.175,350.,109.4); // prepreg3 300 microns
+ TGeoTubeSeg *t5 = new TGeoTubeSeg(76.7612,78.7612,74.175,350.,109.4); // nomex 2 cm
+ TGeoTubeSeg *tepox1 = new TGeoTubeSeg(76.6774,78.845,74.175,109.4,110.);//epoxy
TGeoTubeSeg *tpr1 = new TGeoTubeSeg(78.845,78.885,74.175,109.,111.);
// volumes for the outer part
//
// segment middle
//
- TGeoTubeSeg *t6 = new TGeoTubeSeg(76.6774,78.795,5.,350.,109.6); // tedlar 38 microns
- TGeoTubeSeg *t7 = new TGeoTubeSeg(76.6812,78.7912,5.,350.,109.6); // prepreg2 250 microns
- TGeoTubeSeg *t8 = new TGeoTubeSeg(76.7062,78.7662,5.,350.,109.6); // prepreg3 300 microns
- TGeoTubeSeg *t9 = new TGeoTubeSeg(76.7362,78.7362,5.,350.,109.6); // nomex 2 cm
- TGeoTubeSeg *tepox2 = new TGeoTubeSeg(76.6774,78.795,5.,109.6,110.);//epoxy
+ TGeoTubeSeg *t6 = new TGeoTubeSeg(76.6774,78.795,5.,350.,109.4); // tedlar 38 microns
+ TGeoTubeSeg *t7 = new TGeoTubeSeg(76.6812,78.7912,5.,350.,109.4); // prepreg2 250 microns
+ TGeoTubeSeg *t8 = new TGeoTubeSeg(76.7062,78.7662,5.,350.,109.4); // prepreg3 300 microns
+ TGeoTubeSeg *t9 = new TGeoTubeSeg(76.7362,78.7362,5.,350.,109.4); // nomex 2 cm
+ TGeoTubeSeg *tepox2 = new TGeoTubeSeg(76.6774,78.795,5.,109.4,110.);//epoxy
TGeoTubeSeg *tpr2 = new TGeoTubeSeg(78.795,78.835,5.,109.,111.);
// volumes for the middle part
TGeoVolume *tv6 = new TGeoVolume("TPC_IFC7",t6,sm2);
//
// segment central part
//
- TGeoTubeSeg *t10 = new TGeoTubeSeg(76.6774,78.785,93.75,350.,109.6); // tedlar 38 microns
- TGeoTubeSeg *t11 = new TGeoTubeSeg(76.6812,78.7812,93.75,350.,109.6); // prepreg3 500 microns
- TGeoTubeSeg *t12 = new TGeoTubeSeg(76.7312,78.7312,93.75,350.,109.6); // nomex 2 cm
- TGeoTubeSeg *tepox3 = new TGeoTubeSeg(76.6774,78.785,93.75,109.6,110.);//epoxy
+ TGeoTubeSeg *t10 = new TGeoTubeSeg(76.6774,78.785,93.75,350.,109.4); // tedlar 38 microns
+ TGeoTubeSeg *t11 = new TGeoTubeSeg(76.6812,78.7812,93.75,350.,109.4); // prepreg3 500 microns
+ TGeoTubeSeg *t12 = new TGeoTubeSeg(76.7312,78.7312,93.75,350.,109.4); // nomex 2 cm
+ TGeoTubeSeg *tepox3 = new TGeoTubeSeg(76.6774,78.785,93.75,109.4,110.);//epoxy
TGeoTubeSeg *tpr3 = new TGeoTubeSeg(78.785,78.825,93.75,109.,111.);
// volumes for the central part
TGeoVolume *tv10 = new TGeoVolume("TPC_IFC11",t10,sm2);
//
TGeoVolume *cflv = new TGeoVolume("TPC_CDR",cfl,m3);
// sandwich
- TGeoTube *cd1 = new TGeoTube(60.6224,61.19,71.1);
- TGeoTube *cd2 = new TGeoTube(60.6262,61.1862,71.1);
- TGeoTube *cd3 = new TGeoTube(60.6462,61.1662,71.1);
- TGeoTube *cd4 = new TGeoTube(60.6562,61.1562,71.1);
+ TGeoTubeSeg *cd1 = new TGeoTubeSeg(60.6224,61.19,71.1,0.1,119.9);
+ TGeoTubeSeg *cd2 = new TGeoTubeSeg(60.6262,61.1862,71.1,0.1,119.9);
+ TGeoTubeSeg *cd3 = new TGeoTubeSeg(60.6462,61.1662,71.1,0.1,119.9);
+ TGeoTubeSeg *cd4 = new TGeoTubeSeg(60.6562,61.1562,71.1,0.1,119.9);
+ TGeoTubeSeg *tepox4 = new TGeoTubeSeg(60.6224,61.19,71.1,359.9,0.1);
+ // TGeoTube *cd1 = new TGeoTube(60.6224,61.19,71.1);
+ // TGeoTube *cd2 = new TGeoTube(60.6262,61.1862,71.1);
+ // TGeoTube *cd3 = new TGeoTube(60.6462,61.1662,71.1);
+ // TGeoTube *cd4 = new TGeoTube(60.6562,61.1562,71.1);
//
TGeoMedium *sm6 = gGeoManager->GetMedium("TPC_Prepreg1");
TGeoMedium *sm8 = gGeoManager->GetMedium("TPC_Epoxyfm");
TGeoVolume *cd2v = new TGeoVolume("TPC_CDR2",cd2,sm6);// prepreg1
TGeoVolume *cd3v = new TGeoVolume("TPC_CDR3",cd3,sm8); //epoxy film
TGeoVolume *cd4v = new TGeoVolume("TPC_CDR4",cd4,sm4); //nomex
+ TGeoVolume *tvep4 = new TGeoVolume("TPC_IFEPOX4",tepox4,sm1);
//
// seals for central drum 2 copies
//
TGeoRotation *ref = new TGeoRotation("ref",90.,0.,90.,90.,180.,0.);
//
- cd1v->AddNode(cd2v,1); cd2v->AddNode(cd3v,1); cd3v->AddNode(cd4v,1); cflv->AddNode(cd1v,1);
+ cd1v->AddNode(cd2v,1); cd2v->AddNode(cd3v,1); cd3v->AddNode(cd4v,1); //sandwich
+ // first segment
+ cflv->AddNode(cd1v,1); cflv->AddNode(tvep4,1);
+ // second segment
+ segrot = new TGeoRotation();
+ segrot->RotateZ(120.);
+ cflv->AddNode(cd1v,2,segrot); cflv->AddNode(tvep4,2,segrot);
+ // third segment
+ segrot = new TGeoRotation();
+ segrot->RotateZ(240.);
+ cflv->AddNode(cd1v,3,segrot); cflv->AddNode(tvep4,3,segrot);
//
v1->AddNode(siv,1,new TGeoTranslation(0.,0.,-69.9));
v1->AddNode(siv,2,new TGeoTranslation(0.,0.,69.9));
//
TGeoTube *ih = new TGeoTube(81.05,84.05,0.3);
TGeoTube *oh = new TGeoTube(250.,256.,0.5);
- TGeoTube *mem = new TGeoTube(84.05,250.,0.01);
+ TGeoTube *mem = new TGeoTube(84.05,250.,0.00115);
//
TGeoMedium *m4 = gGeoManager->GetMedium("TPC_G10");
gMC->SetMaxNStep(-30000); // max. number of steps increased
if (fPrimaryIonisation) {
+ // for FLUKA
gMC->Gstpar(idtmed[2],"PRIMIO_E", 20.77); // 1st ionisation potential
gMC->Gstpar(idtmed[2],"PRIMIO_N", 14.35);
gMC->Gstpar(idtmed[2],"LOSS", 14); // specific energy loss
gMC->Gstpar(idtmed[2],"STRA",4);
- } else {
- gMC->Gstpar(idtmed[2],"LOSS", 5); // specific energy loss
- }
+ }
+ // specific energy loss for geant3 is now defined in galice.cuts
AliDebug(1,"*** TPC version 2 initialized ***");
hits[4]=gMC->TrackTime();
AddHit(gAlice->GetMCApp()->GetCurrentTrackNumber(), vol,hits);
-
+ if (fDebugStreamer){
+ // You can dump here what you need
+ // function CreateDebugStremer() to be called in the Config.C macro
+ // if you want to enable it
+ // By default debug streaemer is OFF
+ Float_t edep = gMC->Edep();
+ Float_t tstep = gMC->TrackStep();
+ Int_t pid=gMC->TrackPid();
+ (*fDebugStreamer)<<"hit"<<
+ "x="<<hits[0]<< // hit position
+ "y="<<hits[1]<<
+ "z="<<hits[2]<<
+ "nel="<<hits[3]<< // number of electorns
+ "tof="<<hits[4]<< // hit TOF
+ "edep="<<edep<< // energy deposit
+ "pid="<<pid<< // pid
+ "step="<<tstep<<
+ "p.="<<&p<<
+ "\n";
+ }
+
} // step>0
} //within sector's limits
// Stemax calculation for the next step