[u/mrichter/AliRoot.git] / macros / ConfigBeautyPPR.C

void Config()
{
 
  //
  // Set Random Number seed
  TDatime dt;
  UInt_t curtime=dt.Get();
  UInt_t procid=gSystem->GetPid();
  UInt_t seed=curtime-procid;

  gRandom->SetSeed(seed);
  cerr<<"Seed for random number generation= "<<seed<<endl; 

  new AliGeant3("C++ Interface to Geant3");

  //=======================================================================
  //  Create the output file
   
  TFile *rootfile = new TFile("galice.root","recreate");
  rootfile->SetCompressionLevel(2);
  TGeant3 *geant3 = (TGeant3*)gMC;

  //
  // Set External decayer
  AliDecayer* decayer = new AliDecayerPythia();
  decayer->SetForceDecay(kAll);
  decayer->Init();
  gMC->SetExternalDecayer(decayer);

  //
  //
  //=======================================================================
  // ******* GEANT STEERING parameters FOR ALICE SIMULATION *******
  geant3->SetTRIG(1); //Number of events to be processed 
  geant3->SetSWIT(4,10);
  geant3->SetDEBU(0,0,1);
  //geant3->SetSWIT(2,2);
  geant3->SetDCAY(1);
  geant3->SetPAIR(1);
  geant3->SetCOMP(1);
  geant3->SetPHOT(1);
  geant3->SetPFIS(0);
  geant3->SetDRAY(0);
  geant3->SetANNI(1);
  geant3->SetBREM(1);
  geant3->SetMUNU(1);
  geant3->SetCKOV(1);
  geant3->SetHADR(1); //Select pure GEANH (HADR 1) or GEANH/NUCRIN (HADR 3)
  geant3->SetLOSS(2);
  geant3->SetMULS(1);
  geant3->SetRAYL(1);
  geant3->SetAUTO(1); //Select automatic STMIN etc... calc. (AUTO 1) or manual (AUTO 0)
  geant3->SetABAN(0); //Restore 3.16 behaviour for abandoned tracks
  geant3->SetOPTI(2); //Select optimisation level for GEANT geometry searches (0,1,2)
  geant3->SetERAN(5.e-7);

  Float_t cut    = 1.e-3; // 1MeV cut by default
  Float_t tofmax = 1.e10;
  //             GAM ELEC NHAD CHAD MUON EBREM MUHAB EDEL MUDEL MUPA TOFMAX
  geant3->SetCUTS(cut,cut, cut, cut, cut, cut,  cut,  cut, cut,  cut, tofmax);
  //
  //=======================================================================
  // ************* STEERING parameters FOR ALICE SIMULATION **************
  // --- Specify event type to be tracked through the ALICE setup
  // --- All positions are in cm, angles in degrees, and P and E in GeV
  // AliGenPythia *gener = new AliGenPythia(ntracks);
  AliGenPythia *gener = new AliGenPythia(-1);

  gener->SetMomentumRange(0,999);
  gener->SetPhiRange(-180.,180.);
  gener->SetThetaRange(0,180);
  gener->SetYRange(-999,999);
  //gener->SetPtRange(0,100);
  gener->SetOrigin(0,0,0);          // vertex position
  //gener->SetVertexSmear(kPerEvent); 
  gener->SetSigma(0,0,5.3);  // Sigma in (X,Y,Z) (cm) on IP position
  gener->SetTrackingFlag(0);
  // gener->SetForceDecay(kHadronicD);

  //
  // The following settings select the Pythia parameters tuned to agree
  // with beauty NLO calculation for Pb-Pb @ 5.5 TeV with MNR code.
  //
  gener->SetProcess(kPyBeautyPbMNR);
  gener->SetStrucFunc(kCTEQ4L);
  gener->SetPtHard(2.75,-1.0);
  gener->SetEnergyCMS(5500.);
  gener->SetNuclei(208,208); // Pb-Pb collisions
  // Force no decay heavy quark mesons
  gener->SetForceDecay(kNoDecayHeavy);

  gener->Init();
  // 
  // Activate this line if you want the vertex smearing to happen
  // track by track
  //
  //gener->SetVertexSmear(perTrack); 

  // Field (L3 0.4 T)
  AliMagFCM* field = new AliMagFCM(
	      "Map2","$(ALICE_ROOT)/data/field01.dat", 2, 1., 10.);
  field->SetSolenoidField(4.);
  gAlice->SetField(field);    

  Int_t iABSO=0;
  Int_t iCASTOR=0;
  Int_t iDIPO=0;
  Int_t iFMD=0;
  Int_t iFRAME=0;
  Int_t iHALL=0;
  Int_t iITS=0;
  Int_t iMAG=0;
  Int_t iMUON=0;
  Int_t iPHOS=0;
  Int_t iPIPE=0;
  Int_t iPMD=0;
  Int_t iRICH=0;
  Int_t iSHIL=0;
  Int_t iSTART=0;
  Int_t iTOF=0;
  Int_t iTPC=0;
  Int_t iTRD=0;
  Int_t iZDC=0;

  //=================== Alice BODY parameters =============================
  AliBODY *BODY = new AliBODY("BODY","Alice envelop");

  if(iMAG) {
    //=================== MAG parameters ============================
    // --- Start with Magnet since detector layouts may be depending ---
    // --- on the selected Magnet dimensions ---
    AliMAG *MAG  = new AliMAG("MAG","Magnet");
  }


  if(iABSO) {
    //=================== ABSO parameters ============================
    AliABSO *ABSO  = new AliABSOv0("ABSO","Muon Absorber");
  }

  if(iDIPO) {
    //=================== DIPO parameters ============================

    AliDIPO *DIPO  = new AliDIPOv2("DIPO","Dipole version 2");
  }

  if(iHALL) {
    //=================== HALL parameters ============================

    AliHALL *HALL  = new AliHALL("HALL","Alice Hall");
  }


  if(iFRAME) {
    //=================== FRAME parameters ============================

    AliFRAME *FRAME  = new AliFRAMEv2("FRAME","Space Frame");

  }

  if(iSHIL) {
    //=================== SHIL parameters ============================

    AliSHIL *SHIL  = new AliSHILv2("SHIL","Shielding");
  }


  if(iPIPE) {
    //=================== PIPE parameters ============================

    AliPIPE *PIPE  = new AliPIPEv0("PIPE","Beam Pipe");
  }


  if(iITS) {

//=================== ITS parameters ============================
    //
    // As the innermost detector in ALICE, the Inner Tracking System "impacts" on
    // almost all other detectors. This involves the fact that the ITS geometry
    // still has several options to be followed in parallel in order to determine
    // the best set-up which minimizes the induced background. All the geometries
    // available to date are described in the following. Read carefully the comments
    // and use the default version (the only one uncommented) unless you are making
    // comparisons and you know what you are doing. In this case just uncomment the
    // ITS geometry you want to use and run Aliroot.
    //
    // Detailed geometries:         
    //
    //
    //AliITS *ITS  = new AliITSv5symm("ITS","Updated ITS TDR detailed version with symmetric services");
    //
    //AliITS *ITS  = new AliITSv5asymm("ITS","Updates ITS TDR detailed version with asymmetric services");
    //
    AliITSvPPRasymm *ITS  = new AliITSvPPRasymm("ITS","New ITS PPR detailed version with asymmetric services");
    ITS->SetMinorVersion(2);                                         // don't touch this parameter if you're not an ITS developer
    ITS->SetReadDet(kFALSE);                                         // don't touch this parameter if you're not an ITS developer
    ITS->SetWriteDet("$ALICE_ROOT/ITS/ITSgeometry_vPPRasymm2.det");  // don't touch this parameter if you're not an ITS developer
    ITS->SetThicknessDet1(200.);   // detector thickness on layer 1 must be in the range [150,300]
    ITS->SetThicknessDet2(200.);   // detector thickness on layer 2 must be in the range [150,300]
    ITS->SetThicknessChip1(200.);  // chip thickness on layer 1 must be in the range [100,300]
    ITS->SetThicknessChip2(200.);  // chip thickness on layer 2 must be in the range [100,300]
    ITS->SetRails(1);	           // 1 --> rails in ; 0 --> rails out
    ITS->SetCoolingFluid(1);       // 1 --> water ; 0 --> freon
    //
    //AliITSvPPRsymm *ITS  = new AliITSvPPRsymm("ITS","New ITS PPR detailed version with symmetric services");
    //ITS->SetMinorVersion(2);                                       // don't touch this parameter if you're not an ITS developer
    //ITS->SetReadDet(kFALSE);                                       // don't touch this parameter if you're not an ITS developer
    //ITS->SetWriteDet("$ALICE_ROOT/ITS/ITSgeometry_vPPRsymm2.det"); // don't touch this parameter if you're not an ITS developer
    //ITS->SetThicknessDet1(300.);   // detector thickness on layer 1 must be in the range [150,300]
    //ITS->SetThicknessDet2(300.);   // detector thickness on layer 2 must be in the range [150,300]
    //ITS->SetThicknessChip1(300.);  // chip thickness on layer 1 must be in the range [100,300]
    //ITS->SetThicknessChip2(300.);  // chip thickness on layer 2 must be in the range [100,300]
    //ITS->SetRails(1);              // 1 --> rails in ; 0 --> rails out
    //ITS->SetCoolingFluid(1);       // 1 --> water ; 0 --> freon
    //
    //
    // Coarse geometries (warning: no hits are produced with these coarse geometries and they unuseful 
    // for reconstruction !):
    //                                                     
    //
    //AliITSvPPRcoarseasymm *ITS  = new AliITSvPPRcoarseasymm("ITS","New ITS PPR coarse version with asymmetric services");
    //ITS->SetRails(1);                // 1 --> rails in ; 0 --> rails out
    //ITS->SetSupportMaterial(0);      // 0 --> Copper ; 1 --> Aluminum ; 2 --> Carbon
    //
    //AliITS *ITS  = new AliITSvPPRcoarsesymm("ITS","New ITS PPR coarse version with symmetric services");
    //ITS->SetRails(1);                // 1 --> rails in ; 0 --> rails out
    //ITS->SetSupportMaterial(0);      // 0 --> Copper ; 1 --> Aluminum ; 2 --> Carbon
    //                      
    //
    //
    // Geant3 <-> EUCLID conversion
    // ============================
    //
    // SetEUCLID is a flag to output (=1) or not to output (=0) both geometry and
    // media to two ASCII files (called by default ITSgeometry.euc and
    // ITSgeometry.tme) in a format understandable to the CAD system EUCLID.
    // The default (=0) means that you dont want to use this facility.
    //
    ITS->SetEUCLID(0);  
  }
  

  if(iTPC) {
    //============================ TPC parameters ================================
    // --- This allows the user to specify sectors for the SLOW (TPC geometry 2)
    // --- Simulator. SecAL (SecAU) <0 means that ALL lower (upper)
    // --- sectors are specified, any value other than that requires at least one 
    // --- sector (lower or upper)to be specified!
    // --- Reminder: sectors 1-24 are lower sectors (1-12 -> z>0, 13-24 -> z<0)
    // ---           sectors 25-72 are the upper ones (25-48 -> z>0, 49-72 -> z<0)
    // --- SecLows - number of lower sectors specified (up to 6)
    // --- SecUps - number of upper sectors specified (up to 12)
    // --- Sens - sensitive strips for the Slow Simulator !!!
    // --- This does NOT work if all S or L-sectors are specified, i.e.
    // --- if SecAL or SecAU < 0
    //
    //
    //-----------------------------------------------------------------------------

    //  gROOT->LoadMacro("SetTPCParam.C");
    //  AliTPCParam *param = SetTPCParam();
    AliTPC *TPC  = new AliTPCv2("TPC","Default");
    // All sectors included 
    TPC->SetSecAL(-1);
    TPC->SetSecAU(-1);

  }


  if(iTOF) {
    //=================== TOF parameters ============================
    AliTOF *TOF  = new AliTOFv2("TOF","normal TOF");
  }

  if(iRICH) {
    //=================== RICH parameters ===========================
    AliRICH *RICH  = new AliRICHv1("RICH","normal RICH");    

  }


  if(iZDC) {
    //=================== ZDC parameters ============================

    AliZDC *ZDC  = new AliZDCv1("ZDC","normal ZDC");
  }

  if(iCASTOR) {
    //=================== CASTOR parameters ============================

    AliCASTOR *CASTOR  = new AliCASTORv1("CASTOR","normal CASTOR");
  }

  if(iTRD) {
    //=================== TRD parameters ============================
  
    AliTRD *TRD  = new AliTRDv1("TRD","TRD slow simulator");
  
    // Select the gas mixture (0: 97% Xe + 3% isobutane, 1: 90% Xe + 10% CO2)
    TRD->SetGasMix(1);
  
    // With hole in front of PHOS
    TRD->SetPHOShole();
    // With hole in front of RICH
    TRD->SetRICHhole();
    // Switch on TR
    AliTRDsim *TRDsim = TRD->CreateTR();
  }

  if(iFMD) {
    //=================== FMD parameters ============================

    AliFMD *FMD  = new AliFMDv0("FMD","normal FMD");
  }

  if(iMUON) {
    //=================== MUON parameters ===========================
    AliMUON *MUON  = new AliMUONv1("MUON","default");
  }
 
  //=================== PHOS parameters ===========================

  if(iPHOS) {
    AliPHOS *PHOS  = new AliPHOSv1("PHOS","GPS2");
  }


  if(iPMD) {
    //=================== PMD parameters ============================

    AliPMD *PMD  = new AliPMDv1("PMD","normal PMD");
    PMD->SetPAR(1., 1., 0.8, 0.02);
    PMD->SetIN(6., 18., -580., 27., 27.);
    PMD->SetGEO(0.0, 0.2, 4.);
    PMD->SetPadSize(0.8, 1.0, 1.0, 1.5);

  }

  if(iSTART) {
    //=================== START parameters ============================
    AliSTART *START  = new AliSTARTv1("START","START Detector");
  }

         
}
Commit	Line	Data
642b80f3	1	void Config()
	2	{
	3
	4	//
	5	// Set Random Number seed
	6	TDatime dt;
	7	UInt_t curtime=dt.Get();
	8	UInt_t procid=gSystem->GetPid();
	9	UInt_t seed=curtime-procid;
	10
	11	gRandom->SetSeed(seed);
	12	cerr<<"Seed for random number generation= "<<seed<<endl;
	13
	14	new AliGeant3("C++ Interface to Geant3");
	15
	16	//=======================================================================
	17	// Create the output file
	18
	19	TFile *rootfile = new TFile("galice.root","recreate");
	20	rootfile->SetCompressionLevel(2);
	21	TGeant3 geant3 = (TGeant3)gMC;
	22
	23	//
	24	// Set External decayer
	25	AliDecayer* decayer = new AliDecayerPythia();
	26	decayer->SetForceDecay(kAll);
	27	decayer->Init();
	28	gMC->SetExternalDecayer(decayer);
	29
	30	//
	31	//
	32	//=======================================================================
	33	// ***** GEANT STEERING parameters FOR ALICE SIMULATION *****
	34	geant3->SetTRIG(1); //Number of events to be processed
	35	geant3->SetSWIT(4,10);
	36	geant3->SetDEBU(0,0,1);
	37	//geant3->SetSWIT(2,2);
	38	geant3->SetDCAY(1);
	39	geant3->SetPAIR(1);
	40	geant3->SetCOMP(1);
	41	geant3->SetPHOT(1);
	42	geant3->SetPFIS(0);
	43	geant3->SetDRAY(0);
	44	geant3->SetANNI(1);
	45	geant3->SetBREM(1);
	46	geant3->SetMUNU(1);
	47	geant3->SetCKOV(1);
	48	geant3->SetHADR(1); //Select pure GEANH (HADR 1) or GEANH/NUCRIN (HADR 3)
	49	geant3->SetLOSS(2);
	50	geant3->SetMULS(1);
	51	geant3->SetRAYL(1);
	52	geant3->SetAUTO(1); //Select automatic STMIN etc... calc. (AUTO 1) or manual (AUTO 0)
	53	geant3->SetABAN(0); //Restore 3.16 behaviour for abandoned tracks
	54	geant3->SetOPTI(2); //Select optimisation level for GEANT geometry searches (0,1,2)
	55	geant3->SetERAN(5.e-7);
	56
	57	Float_t cut = 1.e-3; // 1MeV cut by default
	58	Float_t tofmax = 1.e10;
	59	// GAM ELEC NHAD CHAD MUON EBREM MUHAB EDEL MUDEL MUPA TOFMAX
	60	geant3->SetCUTS(cut,cut, cut, cut, cut, cut, cut, cut, cut, cut, tofmax);
	61	//
	62	//=======================================================================
	63	// *********** STEERING parameters FOR ALICE SIMULATION ************
	64	// --- Specify event type to be tracked through the ALICE setup
65	// --- All positions are in cm, angles in degrees, and P and E in GeV
66	// AliGenPythia *gener = new AliGenPythia(ntracks);
67	AliGenPythia *gener = new AliGenPythia(-1);
68
69	gener->SetMomentumRange(0,999);
70	gener->SetPhiRange(-180.,180.);
71	gener->SetThetaRange(0,180);
72	gener->SetYRange(-999,999);
73	//gener->SetPtRange(0,100);
74	gener->SetOrigin(0,0,0); // vertex position
75	//gener->SetVertexSmear(kPerEvent);
76	gener->SetSigma(0,0,5.3); // Sigma in (X,Y,Z) (cm) on IP position
77	gener->SetTrackingFlag(0);
78	// gener->SetForceDecay(kHadronicD);
79
80	//
81	// The following settings select the Pythia parameters tuned to agree
82	// with beauty NLO calculation for Pb-Pb @ 5.5 TeV with MNR code.
83	//
84	gener->SetProcess(kPyBeautyPbMNR);
85	gener->SetStrucFunc(kCTEQ4L);
86	gener->SetPtHard(2.75,-1.0);
87	gener->SetEnergyCMS(5500.);
88	gener->SetNuclei(208,208); // Pb-Pb collisions
89	// Force no decay heavy quark mesons
90	gener->SetForceDecay(kNoDecayHeavy);
91
92	gener->Init();
93	//
94	// Activate this line if you want the vertex smearing to happen
95	// track by track
96	//
97	//gener->SetVertexSmear(perTrack);
98
99	// Field (L3 0.4 T)
100	AliMagFCM* field = new AliMagFCM(
101	"Map2","$(ALICE_ROOT)/data/field01.dat", 2, 1., 10.);
102	field->SetSolenoidField(4.);
103	gAlice->SetField(field);
104
105	Int_t iABSO=0;
106	Int_t iCASTOR=0;
107	Int_t iDIPO=0;
108	Int_t iFMD=0;
109	Int_t iFRAME=0;
110	Int_t iHALL=0;
111	Int_t iITS=0;
112	Int_t iMAG=0;
113	Int_t iMUON=0;
114	Int_t iPHOS=0;
115	Int_t iPIPE=0;
116	Int_t iPMD=0;
117	Int_t iRICH=0;
118	Int_t iSHIL=0;
119	Int_t iSTART=0;
120	Int_t iTOF=0;
121	Int_t iTPC=0;
122	Int_t iTRD=0;
123	Int_t iZDC=0;
124
125	//=================== Alice BODY parameters =============================
126	AliBODY *BODY = new AliBODY("BODY","Alice envelop");
127
128	if(iMAG) {
129	//=================== MAG parameters ============================
130	// --- Start with Magnet since detector layouts may be depending ---
131	// --- on the selected Magnet dimensions ---
132	AliMAG *MAG = new AliMAG("MAG","Magnet");
133	}
134
135
136	if(iABSO) {
137	//=================== ABSO parameters ============================
138	AliABSO *ABSO = new AliABSOv0("ABSO","Muon Absorber");
139	}
140
141	if(iDIPO) {
142	//=================== DIPO parameters ============================
143
144	AliDIPO *DIPO = new AliDIPOv2("DIPO","Dipole version 2");
145	}
146
147	if(iHALL) {
148	//=================== HALL parameters ============================
149
150	AliHALL *HALL = new AliHALL("HALL","Alice Hall");
151	}
152
153
154	if(iFRAME) {
155	//=================== FRAME parameters ============================
156
157	AliFRAME *FRAME = new AliFRAMEv2("FRAME","Space Frame");
158
159	}
160
161	if(iSHIL) {
162	//=================== SHIL parameters ============================
163
e0559846	164	AliSHIL *SHIL = new AliSHILv2("SHIL","Shielding");
642b80f3	165	}
	166
	167
	168	if(iPIPE) {
	169	//=================== PIPE parameters ============================
	170
	171	AliPIPE *PIPE = new AliPIPEv0("PIPE","Beam Pipe");
	172	}
	173
	174
	175	if(iITS) {
	176
	177	//=================== ITS parameters ============================
	178	//
	179	// As the innermost detector in ALICE, the Inner Tracking System "impacts" on
	180	// almost all other detectors. This involves the fact that the ITS geometry
	181	// still has several options to be followed in parallel in order to determine
	182	// the best set-up which minimizes the induced background. All the geometries
	183	// available to date are described in the following. Read carefully the comments
	184	// and use the default version (the only one uncommented) unless you are making
	185	// comparisons and you know what you are doing. In this case just uncomment the
	186	// ITS geometry you want to use and run Aliroot.
	187	//
	188	// Detailed geometries:
	189	//
	190	//
	191	//AliITS *ITS = new AliITSv5symm("ITS","Updated ITS TDR detailed version with symmetric services");
	192	//
	193	//AliITS *ITS = new AliITSv5asymm("ITS","Updates ITS TDR detailed version with asymmetric services");
	194	//
	195	AliITSvPPRasymm *ITS = new AliITSvPPRasymm("ITS","New ITS PPR detailed version with asymmetric services");
	196	ITS->SetMinorVersion(2); // don't touch this parameter if you're not an ITS developer
	197	ITS->SetReadDet(kFALSE); // don't touch this parameter if you're not an ITS developer
	198	ITS->SetWriteDet("$ALICE_ROOT/ITS/ITSgeometry_vPPRasymm2.det"); // don't touch this parameter if you're not an ITS developer
	199	ITS->SetThicknessDet1(200.); // detector thickness on layer 1 must be in the range [150,300]
	200	ITS->SetThicknessDet2(200.); // detector thickness on layer 2 must be in the range [150,300]
	201	ITS->SetThicknessChip1(200.); // chip thickness on layer 1 must be in the range [100,300]
	202	ITS->SetThicknessChip2(200.); // chip thickness on layer 2 must be in the range [100,300]
	203	ITS->SetRails(1); // 1 --> rails in ; 0 --> rails out
	204	ITS->SetCoolingFluid(1); // 1 --> water ; 0 --> freon
	205	//
	206	//AliITSvPPRsymm *ITS = new AliITSvPPRsymm("ITS","New ITS PPR detailed version with symmetric services");
	207	//ITS->SetMinorVersion(2); // don't touch this parameter if you're not an ITS developer
	208	//ITS->SetReadDet(kFALSE); // don't touch this parameter if you're not an ITS developer
	209	//ITS->SetWriteDet("$ALICE_ROOT/ITS/ITSgeometry_vPPRsymm2.det"); // don't touch this parameter if you're not an ITS developer
	210	//ITS->SetThicknessDet1(300.); // detector thickness on layer 1 must be in the range [150,300]
	211	//ITS->SetThicknessDet2(300.); // detector thickness on layer 2 must be in the range [150,300]
	212	//ITS->SetThicknessChip1(300.); // chip thickness on layer 1 must be in the range [100,300]
	213	//ITS->SetThicknessChip2(300.); // chip thickness on layer 2 must be in the range [100,300]
	214	//ITS->SetRails(1); // 1 --> rails in ; 0 --> rails out
	215	//ITS->SetCoolingFluid(1); // 1 --> water ; 0 --> freon
	216	//
	217	//
	218	// Coarse geometries (warning: no hits are produced with these coarse geometries and they unuseful
	219	// for reconstruction !):
	220	//
	221	//
	222	//AliITSvPPRcoarseasymm *ITS = new AliITSvPPRcoarseasymm("ITS","New ITS PPR coarse version with asymmetric services");
	223	//ITS->SetRails(1); // 1 --> rails in ; 0 --> rails out
	224	//ITS->SetSupportMaterial(0); // 0 --> Copper ; 1 --> Aluminum ; 2 --> Carbon
	225	//
	226	//AliITS *ITS = new AliITSvPPRcoarsesymm("ITS","New ITS PPR coarse version with symmetric services");
	227	//ITS->SetRails(1); // 1 --> rails in ; 0 --> rails out
	228	//ITS->SetSupportMaterial(0); // 0 --> Copper ; 1 --> Aluminum ; 2 --> Carbon
229	//
230	//
231	//
232	// Geant3 <-> EUCLID conversion
233	// ============================
234	//
235	// SetEUCLID is a flag to output (=1) or not to output (=0) both geometry and
236	// media to two ASCII files (called by default ITSgeometry.euc and
237	// ITSgeometry.tme) in a format understandable to the CAD system EUCLID.
238	// The default (=0) means that you dont want to use this facility.
239	//
240	ITS->SetEUCLID(0);
241	}
242
243
244	if(iTPC) {
245	//============================ TPC parameters ================================
246	// --- This allows the user to specify sectors for the SLOW (TPC geometry 2)
247	// --- Simulator. SecAL (SecAU) <0 means that ALL lower (upper)
248	// --- sectors are specified, any value other than that requires at least one
249	// --- sector (lower or upper)to be specified!
250	// --- Reminder: sectors 1-24 are lower sectors (1-12 -> z>0, 13-24 -> z<0)
251	// --- sectors 25-72 are the upper ones (25-48 -> z>0, 49-72 -> z<0)
252	// --- SecLows - number of lower sectors specified (up to 6)
253	// --- SecUps - number of upper sectors specified (up to 12)
254	// --- Sens - sensitive strips for the Slow Simulator !!!
255	// --- This does NOT work if all S or L-sectors are specified, i.e.
256	// --- if SecAL or SecAU < 0
257	//
258	//
259	//-----------------------------------------------------------------------------
260
261	// gROOT->LoadMacro("SetTPCParam.C");
262	// AliTPCParam *param = SetTPCParam();
263	AliTPC *TPC = new AliTPCv2("TPC","Default");
264	// All sectors included
265	TPC->SetSecAL(-1);
266	TPC->SetSecAU(-1);
267
268	}
269
270
271	if(iTOF) {
272	//=================== TOF parameters ============================
273	AliTOF *TOF = new AliTOFv2("TOF","normal TOF");
274	}
275
276	if(iRICH) {
277	//=================== RICH parameters ===========================
278	AliRICH *RICH = new AliRICHv1("RICH","normal RICH");
279
280	}
281
282
283	if(iZDC) {
284	//=================== ZDC parameters ============================
285
286	AliZDC *ZDC = new AliZDCv1("ZDC","normal ZDC");
287	}
288
289	if(iCASTOR) {
290	//=================== CASTOR parameters ============================
291
292	AliCASTOR *CASTOR = new AliCASTORv1("CASTOR","normal CASTOR");
293	}
294
295	if(iTRD) {
296	//=================== TRD parameters ============================
297
298	AliTRD *TRD = new AliTRDv1("TRD","TRD slow simulator");
299
300	// Select the gas mixture (0: 97% Xe + 3% isobutane, 1: 90% Xe + 10% CO2)
301	TRD->SetGasMix(1);
302
303	// With hole in front of PHOS
304	TRD->SetPHOShole();
305	// With hole in front of RICH
306	TRD->SetRICHhole();
307	// Switch on TR
308	AliTRDsim *TRDsim = TRD->CreateTR();
309	}
310
311	if(iFMD) {
312	//=================== FMD parameters ============================
313
314	AliFMD *FMD = new AliFMDv0("FMD","normal FMD");
315	}
316
317	if(iMUON) {
318	//=================== MUON parameters ===========================
319	AliMUON *MUON = new AliMUONv1("MUON","default");
320	}
321
322	//=================== PHOS parameters ===========================
323
324	if(iPHOS) {
325	AliPHOS *PHOS = new AliPHOSv1("PHOS","GPS2");
326	}
327
328
329	if(iPMD) {
330	//=================== PMD parameters ============================
331
332	AliPMD *PMD = new AliPMDv1("PMD","normal PMD");
333	PMD->SetPAR(1., 1., 0.8, 0.02);
334	PMD->SetIN(6., 18., -580., 27., 27.);
335	PMD->SetGEO(0.0, 0.2, 4.);
336	PMD->SetPadSize(0.8, 1.0, 1.0, 1.5);
337
338	}
339
340	if(iSTART) {
341	//=================== START parameters ============================
342	AliSTART *START = new AliSTARTv1("START","START Detector");
343	}
344
345
346	}