[u/mrichter/AliRoot.git] / AliFluka / macro / fluConfig.C

void Config()
{
  // ============================= 
  // Root file
  // ============================= 
 
  cout << "Processing config " << endl;

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

  // ============================= 
  // Alice modules
  // ============================= 

    // Select modules 
/*
  Int_t iABSO=1;
  Int_t iCASTOR=1;
  Int_t iDIPO=1;
  Int_t iFMD=1;
  Int_t iFRAME=1;
  Int_t iHALL=1;
  Int_t iITS=1;
  Int_t iMAG=1;
  Int_t iMUON=1;
  Int_t iPHOS=1;
  Int_t iPIPE=1;
  Int_t iPMD=1;
  Int_t iRICH=1;
  Int_t iSHIL=1;
  Int_t iSTART=1;
  Int_t iTOF=1;
  Int_t iTPC=1;
  Int_t iTRD=1;
  Int_t iZDC=1;
*/
  Int_t iABSO=0;
  Int_t iCASTOR=0;
  Int_t iDIPO=0;
  Int_t iFMD=0;
  Int_t iFRAME=0;
  Int_t iHALL=1;
  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;

    // Exclude detectors that do not work with Geant4
        iCASTOR=0;

// From G3 Config.C
// Without any modification 


  //=================== 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 AliFRAMEv1("FRAME","Space Frame");

  }

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

    AliSHIL *SHIL  = new AliSHILv0("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);
    //ITS->SetReadDet(kFALSE);
    //ITS->SetWriteDet("$ALICE_ROOT/ITS/ITSgeometry_vPPRasymm2.det");
    //ITS->SetThicknessDet1(300.);   // detector thickness on layer 1 must be in the range [100,300]
    //ITS->SetThicknessDet2(300.);   // detector thickness on layer 2 must be in the range [100,300]
    //ITS->SetThicknessChip1(300.);  // chip thickness on layer 1 must be in the range [150,300]
    //ITS->SetThicknessChip2(300.);  // chip thickness on layer 2 must be in the range [150,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);                                      
    //ITS->SetReadDet(kFALSE);
    //ITS->SetWriteDet("$ALICE_ROOT/ITS/ITSgeometry_vPPRsymm2.det");
    //ITS->SetThicknessDet1(300.);   // detector thickness on layer 1 must be in the range [100,300]
    //ITS->SetThicknessDet2(300.);   // detector thickness on layer 2 must be in the range [100,300]
    //ITS->SetThicknessChip1(300.);  // chip thickness on layer 1 must be in the range [150,300]
    //ITS->SetThicknessChip2(300.);  // chip thickness on layer 2 must be in the range [150,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 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 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
04c6e46a	1	void Config()
	2	{
	3	// =============================
	4	// Root file
	5	// =============================
	6
	7	cout << "Processing config " << endl;
	8
	9	// Create the output file
	10	TFile *rootfile = new TFile("galice.root","recreate");
	11	rootfile->SetCompressionLevel(2);
	12
	13	// =============================
	14	// Alice modules
	15	// =============================
	16
	17	// Select modules
	18	/*
	19	Int_t iABSO=1;
	20	Int_t iCASTOR=1;
	21	Int_t iDIPO=1;
	22	Int_t iFMD=1;
	23	Int_t iFRAME=1;
	24	Int_t iHALL=1;
	25	Int_t iITS=1;
	26	Int_t iMAG=1;
	27	Int_t iMUON=1;
	28	Int_t iPHOS=1;
	29	Int_t iPIPE=1;
	30	Int_t iPMD=1;
	31	Int_t iRICH=1;
	32	Int_t iSHIL=1;
	33	Int_t iSTART=1;
	34	Int_t iTOF=1;
	35	Int_t iTPC=1;
	36	Int_t iTRD=1;
	37	Int_t iZDC=1;
	38	*/
	39	Int_t iABSO=0;
	40	Int_t iCASTOR=0;
	41	Int_t iDIPO=0;
	42	Int_t iFMD=0;
	43	Int_t iFRAME=0;
	44	Int_t iHALL=1;
	45	Int_t iITS=0;
	46	Int_t iMAG=0;
	47	Int_t iMUON=0;
	48	Int_t iPHOS=0;
	49	Int_t iPIPE=0;
	50	Int_t iPMD=0;
	51	Int_t iRICH=0;
	52	Int_t iSHIL=0;
	53	Int_t iSTART=0;
	54	Int_t iTOF=0;
	55	Int_t iTPC=0;
	56	Int_t iTRD=0;
	57	Int_t iZDC=0;
	58
	59	// Exclude detectors that do not work with Geant4
	60	iCASTOR=0;
	61
	62	// From G3 Config.C
	63	// Without any modification
	64
65
66	//=================== Alice BODY parameters =============================
67	AliBODY *BODY = new AliBODY("BODY","Alice envelop");
68
69
70	if(iMAG) {
71	//=================== MAG parameters ============================
72	// --- Start with Magnet since detector layouts may be depending ---
73	// --- on the selected Magnet dimensions ---
74	AliMAG *MAG = new AliMAG("MAG","Magnet");
75	}
76
77
78	if(iABSO) {
79	//=================== ABSO parameters ============================
80	AliABSO *ABSO = new AliABSOv0("ABSO","Muon Absorber");
81	}
82
83	if(iDIPO) {
84	//=================== DIPO parameters ============================
85
86	AliDIPO *DIPO = new AliDIPOv2("DIPO","Dipole version 2");
87	}
88
89	if(iHALL) {
90	//=================== HALL parameters ============================
91
92	AliHALL *HALL = new AliHALL("HALL","Alice Hall");
93	}
94
95
96	if(iFRAME) {
97	//=================== FRAME parameters ============================
98
99	AliFRAME *FRAME = new AliFRAMEv1("FRAME","Space Frame");
100
101	}
102
103	if(iSHIL) {
104	//=================== SHIL parameters ============================
105
106	AliSHIL *SHIL = new AliSHILv0("SHIL","Shielding");
107	}
108
109
110	if(iPIPE) {
111	//=================== PIPE parameters ============================
112
113	AliPIPE *PIPE = new AliPIPEv0("PIPE","Beam Pipe");
114	}
115
116
117	if(iITS) {
118	//=================== ITS parameters ============================
119	//
120	// As the innermost detector in ALICE, the Inner Tracking System "impacts" on
121	// almost all other detectors. This involves the fact that the ITS geometry
122	// still has several options to be followed in parallel in order to determine
123	// the best set-up which minimizes the induced background. All the geometries
124	// available to date are described in the following. Read carefully the comments
125	// and use the default version (the only one uncommented) unless you are making
126	// comparisons and you know what you are doing. In this case just uncomment the
127	// ITS geometry you want to use and run Aliroot.
128	//
129	// Detailed geometries:
130	//
131	//
132	//
133	//AliITS *ITS = new AliITSv5symm("ITS","Updated ITS TDR detailed version with symmetric services");
134	//
135	AliITS *ITS = new AliITSv5asymm("ITS","Updates ITS TDR detailed version with asymmetric services");
136	//
137	//AliITSvPPRasymm *ITS = new AliITSvPPRasymm("ITS","New ITS PPR detailed version with asymmetric services");
138	//ITS->SetMinorVersion(2);
139	//ITS->SetReadDet(kFALSE);
140	//ITS->SetWriteDet("$ALICE_ROOT/ITS/ITSgeometry_vPPRasymm2.det");
141	//ITS->SetThicknessDet1(300.); // detector thickness on layer 1 must be in the range [100,300]
142	//ITS->SetThicknessDet2(300.); // detector thickness on layer 2 must be in the range [100,300]
143	//ITS->SetThicknessChip1(300.); // chip thickness on layer 1 must be in the range [150,300]
144	//ITS->SetThicknessChip2(300.); // chip thickness on layer 2 must be in the range [150,300]
145	//ITS->SetRails(1); // 1 --> rails in ; 0 --> rails out
146	//ITS->SetCoolingFluid(1); // 1 --> water ; 0 --> freon
147	//
148	//AliITSvPPRsymm *ITS = new AliITSvPPRsymm("ITS","New ITS PPR detailed version with symmetric services");
149	//ITS->SetMinorVersion(2);
150	//ITS->SetReadDet(kFALSE);
151	//ITS->SetWriteDet("$ALICE_ROOT/ITS/ITSgeometry_vPPRsymm2.det");
152	//ITS->SetThicknessDet1(300.); // detector thickness on layer 1 must be in the range [100,300]
153	//ITS->SetThicknessDet2(300.); // detector thickness on layer 2 must be in the range [100,300]
154	//ITS->SetThicknessChip1(300.); // chip thickness on layer 1 must be in the range [150,300]
155	//ITS->SetThicknessChip2(300.); // chip thickness on layer 2 must be in the range [150,300]
156	//ITS->SetRails(1); // 1 --> rails in ; 0 --> rails out
157	//ITS->SetCoolingFluid(1); // 1 --> water ; 0 --> freon
158	//
159	//
160	// Coarse geometries (warning: no hits are produced with these coarse geometries and they unuseful
161	// for reconstruction !):
162	//
163	//
164	//
165	//AliITSvPPRcoarseasymm *ITS = new AliITSvPPRcoarseasymm("ITS","New ITS coarse version with asymmetric services");
166	//ITS->SetRails(1); // 1 --> rails in ; 0 --> rails out
167	//ITS->SetSupportMaterial(0); // 0 --> Copper ; 1 --> Aluminum ; 2 --> Carbon
168	//
169	//AliITS *ITS = new AliITSvPPRcoarsesymm("ITS","New ITS coarse version with symmetric services");
170	//ITS->SetRails(1); // 1 --> rails in ; 0 --> rails out
171	//ITS->SetSupportMaterial(0); // 0 --> Copper ; 1 --> Aluminum ; 2 --> Carbon
172	//
173	//
174	//
175	// Geant3 <-> EUCLID conversion
176	// ============================
177	//
178	// SetEUCLID is a flag to output (=1) or not to output (=0) both geometry and
179	// media to two ASCII files (called by default ITSgeometry.euc and
180	// ITSgeometry.tme) in a format understandable to the CAD system EUCLID.
181	// The default (=0) means that you dont want to use this facility.
182	//
183	ITS->SetEUCLID(0);
184	}
185
186
187	if(iTPC) {
188	//============================ TPC parameters ================================
189	// --- This allows the user to specify sectors for the SLOW (TPC geometry 2)
190	// --- Simulator. SecAL (SecAU) <0 means that ALL lower (upper)
191	// --- sectors are specified, any value other than that requires at least one
192	// --- sector (lower or upper)to be specified!
193	// --- Reminder: sectors 1-24 are lower sectors (1-12 -> z>0, 13-24 -> z<0)
194	// --- sectors 25-72 are the upper ones (25-48 -> z>0, 49-72 -> z<0)
195	// --- SecLows - number of lower sectors specified (up to 6)
196	// --- SecUps - number of upper sectors specified (up to 12)
197	// --- Sens - sensitive strips for the Slow Simulator !!!
198	// --- This does NOT work if all S or L-sectors are specified, i.e.
199	// --- if SecAL or SecAU < 0
200	//
201	//
202	//-----------------------------------------------------------------------------
203
204	// gROOT->LoadMacro("SetTPCParam.C");
205	// AliTPCParam *param = SetTPCParam();
206	AliTPC *TPC = new AliTPCv2("TPC","Default");
207	// All sectors included
208	TPC->SetSecAL(-1);
209	TPC->SetSecAU(-1);
210
211	}
212
213	if(iTOF) {
214	//=================== TOF parameters ============================
215	AliTOF *TOF = new AliTOFv2("TOF","normal TOF");
216	}
217
218	if(iRICH) {
219	//=================== RICH parameters ===========================
220	AliRICH *RICH = new AliRICHv1("RICH","normal RICH");
221
222	}
223
224
225	if(iZDC) {
226	//=================== ZDC parameters ============================
227
228	AliZDC *ZDC = new AliZDCv1("ZDC","normal ZDC");
229	}
230
231	if(iCASTOR) {
232	//=================== CASTOR parameters ============================
233
234	AliCASTOR *CASTOR = new AliCASTORv1("CASTOR","normal CASTOR");
235	}
236
237	if(iTRD) {
238	//=================== TRD parameters ============================
239
240	AliTRD *TRD = new AliTRDv1("TRD","TRD slow simulator");
241
242	// Select the gas mixture (0: 97% Xe + 3% isobutane, 1: 90% Xe + 10% CO2)
243	TRD->SetGasMix(1);
244
245	// With hole in front of PHOS
246	TRD->SetPHOShole();
247	// With hole in front of RICH
248	TRD->SetRICHhole();
249	// Switch on TR
250	AliTRDsim *TRDsim = TRD->CreateTR();
251	}
252
253	if(iFMD) {
254	//=================== FMD parameters ============================
255
256	AliFMD *FMD = new AliFMDv0("FMD","normal FMD");
257	}
258
259	if(iMUON) {
260	//=================== MUON parameters ===========================
261
262	AliMUON *MUON = new AliMUONv1("MUON","default");
263	}
264
265	//=================== PHOS parameters ===========================
266
267	if(iPHOS) {
268	AliPHOS *PHOS = new AliPHOSv1("PHOS","GPS2");
269	}
270
271
272	if(iPMD) {
273	//=================== PMD parameters ============================
274
275	AliPMD *PMD = new AliPMDv1("PMD","normal PMD");
276	PMD->SetPAR(1., 1., 0.8, 0.02);
277	PMD->SetIN(6., 18., -580., 27., 27.);
278	PMD->SetGEO(0.0, 0.2, 4.);
279	PMD->SetPadSize(0.8, 1.0, 1.0, 1.5);
280
281	}
282
283	if(iSTART) {
284	//=================== START parameters ============================
285	AliSTART *START = new AliSTARTv1("START","START Detector");
286	}
287
288	}