[u/mrichter/AliRoot.git] / MUON / Config_acc.C

enum gentype_t {hijing, gun, box, pythia, param, cocktail, fluka, halo, ntuple, scan, doublescan};

gentype_t gentype=param;
Int_t ntracks=1000;

void Config()
{

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

//=======================================================================
//  Create the output file
   
TFile *rootfile = new TFile("galice.root","recreate");
rootfile->SetCompressionLevel(2);
TGeant3 *geant3 = (TGeant3*)gMC;
 
 
//=======================================================================
// ******* GEANT STEERING parameters FOR ALICE SIMULATION *******
 geant3->fGctrak->maxnst=1000000;
 
geant3->SetTRIG(1);          //Number of events to be processed 
geant3->SetSWIT(4,100);
geant3->SetDEBU(0,0,1);
geant3->SetDCAY(1);
geant3->SetPAIR(1);
geant3->SetCOMP(1);
geant3->SetPHOT(1);
geant3->SetPFIS(0);
geant3->SetDRAY(1);
geant3->SetANNI(1);
geant3->SetBREM(1);
geant3->SetMUNU(1);
geant3->SetCKOV(0);
geant3->SetHADR(4); //Select pure GEANH (HADR 1) or GEANH/NUCRIN (HADR 3)
geant3->SetLOSS(1);
geant3->SetMULS(1);
geant3->SetRAYL(0);
geant3->SetAUTO(0); //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)
Float_t cut    = 1.e-1; // 100MeV cut by default
Float_t tofmax = 1.e10;
//              GAM    ELEC   NHAD   CHAD   MUON  EBREM  MUHAB EDEL MUDEL MUPA TOFMAX
geant3->SetCUTS(1.e-4, 1.e-4, 1.e-3, 1.e-4, 1.e-4, cut,  cut,  cut, cut,  cut, 1.e-5);

 gAlice->TrackingLimits( 700, 2200);
 
//
//=======================================================================
// ************* 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

 switch(gentype)
 {
 case gun:
//*********************************************
// Example for Fixed Particle Gun             *
//*********************************************
     AliGenFixed *gener = new AliGenFixed(ntracks);
     gener->SetMomentum(3);
     gener->SetPhiRange(0);
     gener->SetThetaRange(0);
     gener->SetOrigin(0,0,0);          //vertex position
     gener->SetPart(kProton);          //GEANT particle type
     break;
 case box:  
//*********************************************
// Example for Moving Particle Gun            *
//*********************************************
     AliGenBox *gener = new AliGenBox(ntracks);
     gener->SetMomentumRange(3,4);
     gener->SetPhiRange(-360,360);
     gener->SetThetaRange(2., 10.);
     gener->SetOrigin(25,25,510.5);   
     gener->SetVertexSmear(perTrack); 
     //vertex position
     gener->SetSigma(1.8, 1.8,0);           //Sigma in (X,Y,Z) (cm) on IP position
     gener->SetPart(kProton);                    //GEANT particle type
     break;
 case scan:  
//*********************************************
// Scanning on a grid                         *
//*********************************************
     AliGenScan *gener = new AliGenScan(-1);
     gener->SetMomentumRange(20,20);
     gener->SetPhiRange(90,90);
     gener->SetThetaRange(0,0);
     //vertex position
     gener->SetSigma(3,3,0);           //Sigma in (X,Y,Z) (cm) on IP position
     gener->SetPart(kMuonMinus); 
     gener->SetRange(20, -100, 100, 20, -100, 100, 1, 500, 500);
     break;
 case doublescan:  
//*********************************************
// Scanning on a grid                         *
//*********************************************
     AliGenDoubleScan *gener = new AliGenDoubleScan(-1);
     gener->SetMomentumRange(4,4);
     gener->SetPhiRange(0,360);
     gener->SetThetaRange(0,0);
     //vertex position
     gener->SetSigma(3,3,0);           //Sigma in (X,Y,Z) (cm) on IP position
     gener->SetPart(8); 
     gener->SetRange(20, -100, 100, 20, -100, 100, 1, 500, 500);
     gener->SetDistance(1);
     
     break;
     
 case hijing:
     AliGenHIJINGpara *gener = new AliGenHIJINGpara(ntracks);
     gener->SetMomentumRange(0,999);
     gener->SetPhiRange(0,360);
     gener->SetThetaRange(0.104,33.52);
     gener->SetOrigin(0,0,0);          //vertex position
     gener->SetSigma(0,0,5.6);         //Sigma in (X,Y,Z) (cm) on IP position
     break;

 case pythia:
//********************************************
// Example for Charm  Production with Pythia *
//********************************************

     AliGenPythia *gener = new AliGenPythia(ntracks);
     gener->SetMomentumRange(0,999);
     gener->SetPhiRange(0,360);
     gener->SetThetaRange(0., 180.);
     gener->SetYRange(-10,10);
     gener->SetPtRange(0,100);
     gener->SetOrigin(0,0,0);          // vertex position
     gener->SetVertexSmear(perEvent); 
     gener->SetSigma(0,0,5.6);         // Sigma in (X,Y,Z) (cm) on IP position
//     gener->SetStrucFunc(DO_Set_1);
     gener->SetProcess(mb); 
     gener->SetEnergyCMS(5500.);
     break;
     
 case param:
//*******************************************************
// Example for J/psi  Production from  Parameterisation *
//*******************************************************
     AliGenParam *gener = new AliGenParam(ntracks,upsilon_p);
//     AliGenParam *gener = new AliGenParam(ntracks, jpsi_p);
     gener->SetMomentumRange(0,999);
     gener->SetPhiRange(0,360);
     gener->SetYRange(2.5,4);
     gener->SetPtRange(0,50);
     gener->SetOrigin(0,0,0);          //vertex position
     gener->SetSigma(0,0,0);           //Sigma in (X,Y,Z) (cm) on IP position
     gener->SetForceDecay(dimuon);
     gener->SetCutOnChild(1);
     gener->SetChildThetaRange(2,9);
     gener->SetTrackingFlag(1);

     break;
     
 case fluka:
//*******************************************************
// Example for a FLUKA Boundary Source                  *
//*******************************************************
     AliGenFLUKAsource *gener = new AliGenFLUKAsource(-1);
     gener->AddFile("$(ALICE_ROOT)/data/all32.root"); 
     rootfile->cd();
     gener->SetPartFlag(7);
     gener->SetMomentumRange(0,999);
     gener->SetPhiRange(0,360);
     gener->SetThetaRange(0., 180.); 
     gener->SetAgeMax(1.e-5);
     
//  31.7 events     
     gener->SetFraction(0.0315);     
     break;

 case ntuple:
//*******************************************************
// Example for reading from a external file                  *
//*******************************************************
     AliGenExtFile *gener = new AliGenExtFile(-1); 
     gener->SetFileName("$(ALICE_ROOT)/data/pbpb.root");
     gener->SetThetaRange(0.104,33.52);
     gener->SetOrigin(0,0,0);          //vertex position
     gener->SetSigma(0,0,5.6);         //Sigma in (X,Y,Z) (cm) on IP position
     gener->SetVertexSmear(perTrack); 
     gener->SetTrackingFlag(1);
     break;

 case halo:
//*******************************************************
// Example for Tunnel Halo Source                       *
//*******************************************************
     AliGenHalo *gener = new AliGenHalo(ntracks); 
     gener->SetFileName("/h1/morsch/marsip/marsip5.mu");
     break;
     
 case cocktail:
//*******************************************************
// Example for a Cocktail                               *
//*******************************************************
     AliGenCocktail *gener = new AliGenCocktail();
     gener->SetMomentumRange(0,10);
     gener->SetPhiRange(0,360);
     gener->SetThetaRange(45.,135);
    
     pions   = new AliGenParam(100, pion_p);
//     kaons   = new AliGenParam(10 , kaon_p);
//     protons = new AliGenParam(10 , proton_p);
     gener->AddGenerator(pions  , "Pions"  , 100);
//     gener->AddGenerator(kaons  , "Kaons"  , 10);
//     gener->AddGenerator(protons, "Protons", 10);
	 
     break;
 }
 
// Activate this line if you want the vertex smearing to happen
// track by track
//

gener->Init();
gAlice->SetField(-2,1);    //Specify maximum magnetic field in Tesla (neg. ==> default field)
Int_t iMAG    =0;
Int_t iABSO   =0;
Int_t iDIPO   =0;
Int_t iSHIL   =0;
Int_t iPIPE   =0;
Int_t iMUON   =1;

//=================== 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 AliABSOv1("ABSO","Muon Absorber");
}

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

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


if(iSHIL) {
//=================== SHIL parameters ============================
    AliSHIL *SHIL  = new AliSHILv1("SHIL","Shielding");
}


if(iPIPE) {
//=================== PIPE parameters ============================
    AliPIPE *PIPE  = new AliPIPEv0("PIPE","Beam Pipe");
}


if(iMUON) {
//=================== MUON parameters ===========================

AliMUON *MUON  = new AliMUONv0("MUON","normal MUON");
 MUON->SetAcceptance(kTRUE, 0.1, 12.);
 MUON->SetIshunt(0);
 MUON->SetMaxStepGas(0.1);
 MUON->SetMaxStepAlu(0.1);
}
}
Commit	Line	Data
9efe4ffa	1	enum gentype_t {hijing, gun, box, pythia, param, cocktail, fluka, halo, ntuple, scan, doublescan};
	2
	3	gentype_t gentype=param;
	4	Int_t ntracks=1000;
	5
	6	void Config()
	7	{
	8
	9	new AliGeant3("C++ Interface to Geant3");
	10
	11	//=======================================================================
	12	// Create the output file
	13
	14	TFile *rootfile = new TFile("galice.root","recreate");
	15	rootfile->SetCompressionLevel(2);
	16	TGeant3 geant3 = (TGeant3)gMC;
	17
	18
	19	//=======================================================================
	20	// ***** GEANT STEERING parameters FOR ALICE SIMULATION *****
	21	geant3->fGctrak->maxnst=1000000;
	22
	23	geant3->SetTRIG(1); //Number of events to be processed
	24	geant3->SetSWIT(4,100);
	25	geant3->SetDEBU(0,0,1);
	26	geant3->SetDCAY(1);
	27	geant3->SetPAIR(1);
	28	geant3->SetCOMP(1);
	29	geant3->SetPHOT(1);
	30	geant3->SetPFIS(0);
	31	geant3->SetDRAY(1);
	32	geant3->SetANNI(1);
	33	geant3->SetBREM(1);
	34	geant3->SetMUNU(1);
	35	geant3->SetCKOV(0);
	36	geant3->SetHADR(4); //Select pure GEANH (HADR 1) or GEANH/NUCRIN (HADR 3)
	37	geant3->SetLOSS(1);
	38	geant3->SetMULS(1);
	39	geant3->SetRAYL(0);
	40	geant3->SetAUTO(0); //Select automatic STMIN etc... calc. (AUTO 1) or manual (AUTO 0)
	41	geant3->SetABAN(0); //Restore 3.16 behaviour for abandoned tracks
	42	geant3->SetOPTI(2); //Select optimisation level for GEANT geometry searches (0,1,2)
	43	Float_t cut = 1.e-1; // 100MeV cut by default
	44	Float_t tofmax = 1.e10;
	45	// GAM ELEC NHAD CHAD MUON EBREM MUHAB EDEL MUDEL MUPA TOFMAX
	46	geant3->SetCUTS(1.e-4, 1.e-4, 1.e-3, 1.e-4, 1.e-4, cut, cut, cut, cut, cut, 1.e-5);
	47
	48	gAlice->TrackingLimits( 700, 2200);
	49
	50	//
	51	//=======================================================================
	52	// *********** STEERING parameters FOR ALICE SIMULATION ************
	53	// --- Specify event type to be tracked through the ALICE setup
	54	// --- All positions are in cm, angles in degrees, and P and E in GeV
	55
	56	switch(gentype)
	57	{
	58	case gun:
	59	//*********************************************
	60	// Example for Fixed Particle Gun *
	61	//*********************************************
	62	AliGenFixed *gener = new AliGenFixed(ntracks);
	63	gener->SetMomentum(3);
	64	gener->SetPhiRange(0);
65	gener->SetThetaRange(0);
66	gener->SetOrigin(0,0,0); //vertex position
67	gener->SetPart(kProton); //GEANT particle type
68	break;
69	case box:
70	//*********************************************
71	// Example for Moving Particle Gun *
72	//*********************************************
73	AliGenBox *gener = new AliGenBox(ntracks);
74	gener->SetMomentumRange(3,4);
75	gener->SetPhiRange(-360,360);
76	gener->SetThetaRange(2., 10.);
77	gener->SetOrigin(25,25,510.5);
78	gener->SetVertexSmear(perTrack);
79	//vertex position
80	gener->SetSigma(1.8, 1.8,0); //Sigma in (X,Y,Z) (cm) on IP position
81	gener->SetPart(kProton); //GEANT particle type
82	break;
83	case scan:
84	//*********************************************
85	// Scanning on a grid *
86	//*********************************************
87	AliGenScan *gener = new AliGenScan(-1);
88	gener->SetMomentumRange(20,20);
89	gener->SetPhiRange(90,90);
90	gener->SetThetaRange(0,0);
91	//vertex position
92	gener->SetSigma(3,3,0); //Sigma in (X,Y,Z) (cm) on IP position
93	gener->SetPart(kMuonMinus);
94	gener->SetRange(20, -100, 100, 20, -100, 100, 1, 500, 500);
95	break;
96	case doublescan:
97	//*********************************************
98	// Scanning on a grid *
99	//*********************************************
100	AliGenDoubleScan *gener = new AliGenDoubleScan(-1);
101	gener->SetMomentumRange(4,4);
102	gener->SetPhiRange(0,360);
103	gener->SetThetaRange(0,0);
104	//vertex position
105	gener->SetSigma(3,3,0); //Sigma in (X,Y,Z) (cm) on IP position
106	gener->SetPart(8);
107	gener->SetRange(20, -100, 100, 20, -100, 100, 1, 500, 500);
108	gener->SetDistance(1);
109
110	break;
111
112	case hijing:
113	AliGenHIJINGpara *gener = new AliGenHIJINGpara(ntracks);
114	gener->SetMomentumRange(0,999);
115	gener->SetPhiRange(0,360);
116	gener->SetThetaRange(0.104,33.52);
117	gener->SetOrigin(0,0,0); //vertex position
118	gener->SetSigma(0,0,5.6); //Sigma in (X,Y,Z) (cm) on IP position
119	break;
120
121	case pythia:
122	//********************************************
123	// Example for Charm Production with Pythia *
124	//********************************************
125
126	AliGenPythia *gener = new AliGenPythia(ntracks);
127	gener->SetMomentumRange(0,999);
128	gener->SetPhiRange(0,360);
129	gener->SetThetaRange(0., 180.);
130	gener->SetYRange(-10,10);
131	gener->SetPtRange(0,100);
132	gener->SetOrigin(0,0,0); // vertex position
133	gener->SetVertexSmear(perEvent);
134	gener->SetSigma(0,0,5.6); // Sigma in (X,Y,Z) (cm) on IP position
135	// gener->SetStrucFunc(DO_Set_1);
136	gener->SetProcess(mb);
137	gener->SetEnergyCMS(5500.);
138	break;
139
140	case param:
141	//*******************************************************
142	// Example for J/psi Production from Parameterisation *
143	//*******************************************************
144	AliGenParam *gener = new AliGenParam(ntracks,upsilon_p);
145	// AliGenParam *gener = new AliGenParam(ntracks, jpsi_p);
146	gener->SetMomentumRange(0,999);
147	gener->SetPhiRange(0,360);
148	gener->SetYRange(2.5,4);
149	gener->SetPtRange(0,50);
150	gener->SetOrigin(0,0,0); //vertex position
151	gener->SetSigma(0,0,0); //Sigma in (X,Y,Z) (cm) on IP position
152	gener->SetForceDecay(dimuon);
153	gener->SetCutOnChild(1);
154	gener->SetChildThetaRange(2,9);
155	gener->SetTrackingFlag(1);
156
157	break;
158
159	case fluka:
160	//*******************************************************
161	// Example for a FLUKA Boundary Source *
162	//*******************************************************
163	AliGenFLUKAsource *gener = new AliGenFLUKAsource(-1);
164	gener->AddFile("$(ALICE_ROOT)/data/all32.root");
165	rootfile->cd();
166	gener->SetPartFlag(7);
167	gener->SetMomentumRange(0,999);
168	gener->SetPhiRange(0,360);
169	gener->SetThetaRange(0., 180.);
170	gener->SetAgeMax(1.e-5);
171
172	// 31.7 events
173	gener->SetFraction(0.0315);
174	break;
175
176	case ntuple:
177	//*******************************************************
178	// Example for reading from a external file *
179	//*******************************************************
180	AliGenExtFile *gener = new AliGenExtFile(-1);
181	gener->SetFileName("$(ALICE_ROOT)/data/pbpb.root");
182	gener->SetThetaRange(0.104,33.52);
183	gener->SetOrigin(0,0,0); //vertex position
184	gener->SetSigma(0,0,5.6); //Sigma in (X,Y,Z) (cm) on IP position
185	gener->SetVertexSmear(perTrack);
186	gener->SetTrackingFlag(1);
187	break;
188
189	case halo:
190	//*******************************************************
191	// Example for Tunnel Halo Source *
192	//*******************************************************
193	AliGenHalo *gener = new AliGenHalo(ntracks);
194	gener->SetFileName("/h1/morsch/marsip/marsip5.mu");
195	break;
196
197	case cocktail:
198	//*******************************************************
199	// Example for a Cocktail *
200	//*******************************************************
201	AliGenCocktail *gener = new AliGenCocktail();
202	gener->SetMomentumRange(0,10);
203	gener->SetPhiRange(0,360);
204	gener->SetThetaRange(45.,135);
205
206	pions = new AliGenParam(100, pion_p);
207	// kaons = new AliGenParam(10 , kaon_p);
208	// protons = new AliGenParam(10 , proton_p);
209	gener->AddGenerator(pions , "Pions" , 100);
210	// gener->AddGenerator(kaons , "Kaons" , 10);
211	// gener->AddGenerator(protons, "Protons", 10);
212
213	break;
214	}
215
216	// Activate this line if you want the vertex smearing to happen
217	// track by track
218	//
219
220	gener->Init();
221	gAlice->SetField(-2,1); //Specify maximum magnetic field in Tesla (neg. ==> default field)
222	Int_t iMAG =0;
223	Int_t iABSO =0;
224	Int_t iDIPO =0;
225	Int_t iSHIL =0;
226	Int_t iPIPE =0;
227	Int_t iMUON =1;
228
229	//=================== Alice BODY parameters =============================
230	AliBODY *BODY = new AliBODY("BODY","Alice envelop");
231
232	if(iMAG) {
233	//=================== MAG parameters ============================
234	// --- Start with Magnet since detector layouts may be depending ---
235	// --- on the selected Magnet dimensions ---
236	AliMAG *MAG = new AliMAG("MAG","Magnet");
237	}
238
239	if(iABSO) {
240	//=================== ABSO parameters ============================
241	AliABSO *ABSO = new AliABSOv1("ABSO","Muon Absorber");
242	}
243
244	if(iDIPO) {
245	//=================== DIPO parameters ============================
246
247	AliDIPO *DIPO = new AliDIPOv2("DIPO","Dipole version 2");
248	}
249
250
251	if(iSHIL) {
252	//=================== SHIL parameters ============================
253	AliSHIL *SHIL = new AliSHILv1("SHIL","Shielding");
254	}
255
256
257	if(iPIPE) {
258	//=================== PIPE parameters ============================
259	AliPIPE *PIPE = new AliPIPEv0("PIPE","Beam Pipe");
260	}
261
262
263	if(iMUON) {
264	//=================== MUON parameters ===========================
265
266	AliMUON *MUON = new AliMUONv0("MUON","normal MUON");
267	MUON->SetAcceptance(kTRUE, 0.1, 12.);
268	MUON->SetIshunt(0);
269	MUON->SetMaxStepGas(0.1);
270	MUON->SetMaxStepAlu(0.1);
271	}
272	}