+static Int_t eventsPerRun = 100;
void Config()
{
// 7-DEC-2000 09:00
// Switch on Transition Radiation simulation. 6/12/00 18:00
// iZDC=1 7/12/00 09:00
// ThetaRange is (0., 180.). It was (0.28,179.72) 7/12/00 09:00
+ // Theta range given through pseudorapidity limits 22/6/2001
// Set Random Number seed
// gRandom->SetSeed(12345);
gener->SetMomentumRange(0, 999);
gener->SetPhiRange(0, 360);
- // gener->SetThetaRange(0.28,179.72);
- gener->SetThetaRange(0., 180.);
+ // Set pseudorapidity range from -8 to 8.
+ Float_t thmin = EtaToTheta(8); // theta min. <---> eta max
+ Float_t thmax = EtaToTheta(-8); // theta max. <---> eta min
+ gener->SetThetaRange(thmin,thmax);
gener->SetOrigin(0, 0, 0); //vertex position
gener->SetSigma(0, 0, 0); //Sigma in (X,Y,Z) (cm) on IP position
gener->Init();
gAlice->SetField(-999, 2); //Specify maximum magnetic field in Tesla (neg. ==> default field)
Int_t iABSO = 1;
- Int_t iCASTOR = 1;
+ Int_t iCRT = 1;
Int_t iDIPO = 1;
Int_t iFMD = 1;
Int_t iFRAME = 1;
Int_t iTPC = 1;
Int_t iTRD = 1;
Int_t iZDC = 1;
+ Int_t iEMCAL = 0;
//=================== Alice BODY parameters =============================
AliBODY *BODY = new AliBODY("BODY", "Alice envelop");
{
//=================== SHIL parameters ============================
- AliSHIL *SHIL = new AliSHILv0("SHIL", "Shielding");
+ AliSHIL *SHIL = new AliSHILvF("SHIL", "Shielding");
}
{
//=================== ZDC parameters ============================
- AliZDC *ZDC = new AliZDCv1("ZDC", "normal ZDC");
+ AliZDC *ZDC = new AliZDCv2("ZDC", "normal ZDC");
}
- if (iCASTOR)
+ if (iCRT)
{
//=================== CASTOR parameters ============================
- AliCASTOR *CASTOR = new AliCASTORv1("CASTOR", "normal CASTOR");
+ AliCRT *CRT = new AliCRTv0("CRT", "normal CRT");
}
if (iTRD)
//=================== FMD parameters ============================
AliFMD *FMD = new AliFMDv1("FMD", "normal FMD");
- }
+ FMD->SetRingsSi1(256);
+ FMD->SetRingsSi2(64);
+ FMD->SetSectorsSi1(20);
+ FMD->SetSectorsSi2(24);
+ }
if (iMUON)
{
}
+ if (iEMCAL && !iRICH)
+ {
+ //=================== EMCAL parameters ============================
+ AliEMCAL *EMCAL = new AliEMCALv1("EMCAL", "EMCALArch1a");
+ }
+
if (iSTART)
{
//=================== START parameters ============================
}
+
+Float_t EtaToTheta(Float_t arg){
+ return (180./TMath::Pi())*2.*atan(exp(-arg));
+}