[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()
{
    gSystem->Load("libgeant321");
    new     TGeant3("C++ Interface to Geant3");
    TGeant3 *geant3 = (TGeant3*)gMC;

//=======================================================================
//  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(kPerTrack); 
     //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(kPerEvent); 
     gener->SetSigma(0,0,5.6);         // Sigma in (X,Y,Z) (cm) on IP position
     gener->SetProcess(kPyCharm); 
     gener->SetEnergyCMS(5500.);
     break;
     
 case param:
//*******************************************************
// Example for J/psi  Production from  Parameterisation *
//*******************************************************
     AliGenParam *gener = new AliGenParam(ntracks, AliGenMUONlib::kUpsilon);
     gener->SetMomentumRange(0,999);
     gener->SetPtRange(0,100.);
     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(kDiMuon);
     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(kPerTrack); 
     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, AliGenMUONlib::kPion);
     gener->AddGenerator(pions  , "Pions"  , 100);
     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	{
aee1d6a1	8	gSystem->Load("libgeant321");
	9	new TGeant3("C++ Interface to Geant3");
	10	TGeant3 geant3 = (TGeant3)gMC;
9efe4ffa	11
	12	//=======================================================================
	13	// Create the output file
	14
	15	TFile *rootfile = new TFile("galice.root","recreate");
	16	rootfile->SetCompressionLevel(2);
	17	TGeant3 geant3 = (TGeant3)gMC;
	18
	19
	20	//=======================================================================
	21	// ***** GEANT STEERING parameters FOR ALICE SIMULATION *****
	22	geant3->fGctrak->maxnst=1000000;
	23
	24	geant3->SetTRIG(1); //Number of events to be processed
	25	geant3->SetSWIT(4,100);
	26	geant3->SetDEBU(0,0,1);
	27	geant3->SetDCAY(1);
	28	geant3->SetPAIR(1);
	29	geant3->SetCOMP(1);
	30	geant3->SetPHOT(1);
	31	geant3->SetPFIS(0);
	32	geant3->SetDRAY(1);
	33	geant3->SetANNI(1);
	34	geant3->SetBREM(1);
	35	geant3->SetMUNU(1);
	36	geant3->SetCKOV(0);
	37	geant3->SetHADR(4); //Select pure GEANH (HADR 1) or GEANH/NUCRIN (HADR 3)
	38	geant3->SetLOSS(1);
	39	geant3->SetMULS(1);
	40	geant3->SetRAYL(0);
	41	geant3->SetAUTO(0); //Select automatic STMIN etc... calc. (AUTO 1) or manual (AUTO 0)
	42	geant3->SetABAN(0); //Restore 3.16 behaviour for abandoned tracks
	43	geant3->SetOPTI(2); //Select optimisation level for GEANT geometry searches (0,1,2)
	44	Float_t cut = 1.e-1; // 100MeV cut by default
	45	Float_t tofmax = 1.e10;
	46	// GAM ELEC NHAD CHAD MUON EBREM MUHAB EDEL MUDEL MUPA TOFMAX
	47	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);
	48
	49	gAlice->TrackingLimits( 700, 2200);
	50
	51	//
	52	//=======================================================================
	53	// *********** STEERING parameters FOR ALICE SIMULATION ************
	54	// --- Specify event type to be tracked through the ALICE setup
	55	// --- All positions are in cm, angles in degrees, and P and E in GeV
	56
	57	switch(gentype)
	58	{
	59	case gun:
	60	//*********************************************
	61	// Example for Fixed Particle Gun *
	62	//*********************************************
	63	AliGenFixed *gener = new AliGenFixed(ntracks);
	64	gener->SetMomentum(3);
	65	gener->SetPhiRange(0);
	66	gener->SetThetaRange(0);
	67	gener->SetOrigin(0,0,0); //vertex position
	68	gener->SetPart(kProton); //GEANT particle type
	69	break;
	70	case box:
	71	//*********************************************
	72	// Example for Moving Particle Gun *
	73	//*********************************************
	74	AliGenBox *gener = new AliGenBox(ntracks);
75	gener->SetMomentumRange(3,4);
76	gener->SetPhiRange(-360,360);
77	gener->SetThetaRange(2., 10.);
78	gener->SetOrigin(25,25,510.5);
aee1d6a1	79	gener->SetVertexSmear(kPerTrack);
9efe4ffa	80	//vertex position
	81	gener->SetSigma(1.8, 1.8,0); //Sigma in (X,Y,Z) (cm) on IP position
	82	gener->SetPart(kProton); //GEANT particle type
	83	break;
	84	case scan:
	85	//*********************************************
	86	// Scanning on a grid *
	87	//*********************************************
	88	AliGenScan *gener = new AliGenScan(-1);
	89	gener->SetMomentumRange(20,20);
	90	gener->SetPhiRange(90,90);
	91	gener->SetThetaRange(0,0);
	92	//vertex position
	93	gener->SetSigma(3,3,0); //Sigma in (X,Y,Z) (cm) on IP position
	94	gener->SetPart(kMuonMinus);
	95	gener->SetRange(20, -100, 100, 20, -100, 100, 1, 500, 500);
	96	break;
	97	case doublescan:
	98	//*********************************************
	99	// Scanning on a grid *
	100	//*********************************************
	101	AliGenDoubleScan *gener = new AliGenDoubleScan(-1);
	102	gener->SetMomentumRange(4,4);
	103	gener->SetPhiRange(0,360);
	104	gener->SetThetaRange(0,0);
	105	//vertex position
	106	gener->SetSigma(3,3,0); //Sigma in (X,Y,Z) (cm) on IP position
	107	gener->SetPart(8);
	108	gener->SetRange(20, -100, 100, 20, -100, 100, 1, 500, 500);
	109	gener->SetDistance(1);
	110
	111	break;
	112
	113	case hijing:
	114	AliGenHIJINGpara *gener = new AliGenHIJINGpara(ntracks);
	115	gener->SetMomentumRange(0,999);
	116	gener->SetPhiRange(0,360);
	117	gener->SetThetaRange(0.104,33.52);
	118	gener->SetOrigin(0,0,0); //vertex position
	119	gener->SetSigma(0,0,5.6); //Sigma in (X,Y,Z) (cm) on IP position
	120	break;
	121
	122	case pythia:
	123	//********************************************
	124	// Example for Charm Production with Pythia *
	125	//********************************************
	126
	127	AliGenPythia *gener = new AliGenPythia(ntracks);
	128	gener->SetMomentumRange(0,999);
	129	gener->SetPhiRange(0,360);
	130	gener->SetThetaRange(0., 180.);
	131	gener->SetYRange(-10,10);
	132	gener->SetPtRange(0,100);
	133	gener->SetOrigin(0,0,0); // vertex position
aee1d6a1	134	gener->SetVertexSmear(kPerEvent);
9efe4ffa	135	gener->SetSigma(0,0,5.6); // Sigma in (X,Y,Z) (cm) on IP position
aee1d6a1	136	gener->SetProcess(kPyCharm);
9efe4ffa	137	gener->SetEnergyCMS(5500.);
	138	break;
	139
	140	case param:
	141	//*******************************************************
	142	// Example for J/psi Production from Parameterisation *
	143	//*******************************************************
aee1d6a1	144	AliGenParam *gener = new AliGenParam(ntracks, AliGenMUONlib::kUpsilon);
9efe4ffa	145	gener->SetMomentumRange(0,999);
aee1d6a1	146	gener->SetPtRange(0,100.);
9efe4ffa	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
aee1d6a1	152	gener->SetForceDecay(kDiMuon);
9efe4ffa	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
aee1d6a1	185	gener->SetVertexSmear(kPerTrack);
9efe4ffa	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
aee1d6a1	206	pions = new AliGenParam(100, AliGenMUONlib::kPion);
9efe4ffa	207	gener->AddGenerator(pions , "Pions" , 100);
9efe4ffa	208	break;
	209	}
	210
	211	// Activate this line if you want the vertex smearing to happen
	212	// track by track
	213	//
	214
	215	gener->Init();
	216	gAlice->SetField(-2,1); //Specify maximum magnetic field in Tesla (neg. ==> default field)
	217	Int_t iMAG =0;
	218	Int_t iABSO =0;
	219	Int_t iDIPO =0;
	220	Int_t iSHIL =0;
	221	Int_t iPIPE =0;
	222	Int_t iMUON =1;
	223
	224	//=================== Alice BODY parameters =============================
	225	AliBODY *BODY = new AliBODY("BODY","Alice envelop");
	226
	227	if(iMAG) {
	228	//=================== MAG parameters ============================
	229	// --- Start with Magnet since detector layouts may be depending ---
	230	// --- on the selected Magnet dimensions ---
	231	AliMAG *MAG = new AliMAG("MAG","Magnet");
	232	}
	233
	234	if(iABSO) {
	235	//=================== ABSO parameters ============================
	236	AliABSO *ABSO = new AliABSOv1("ABSO","Muon Absorber");
	237	}
	238
	239	if(iDIPO) {
	240	//=================== DIPO parameters ============================
	241
	242	AliDIPO *DIPO = new AliDIPOv2("DIPO","Dipole version 2");
	243	}
	244
	245
	246	if(iSHIL) {
	247	//=================== SHIL parameters ============================
	248	AliSHIL *SHIL = new AliSHILv1("SHIL","Shielding");
	249	}
	250
	251
	252	if(iPIPE) {
	253	//=================== PIPE parameters ============================
	254	AliPIPE *PIPE = new AliPIPEv0("PIPE","Beam Pipe");
	255	}
	256
	257
	258	if(iMUON) {
	259	//=================== MUON parameters ===========================
	260
	261	AliMUON *MUON = new AliMUONv0("MUON","normal MUON");
	262	MUON->SetAcceptance(kTRUE, 0.1, 12.);
	263	MUON->SetIshunt(0);
	264	MUON->SetMaxStepGas(0.1);
	265	MUON->SetMaxStepAlu(0.1);
	266	}
	267	}