[u/mrichter/AliRoot.git] / test / vmctest / lego / commonConfig.C

// $Id$
//
// AliRoot Configuration for running aliroot with Monte Carlo.
// commonConfig() includes the common setting for all MCs
// which has to be called before MC is instantiated.
// Called from MC specific configs (g3Config.C, g4Config.C).
//
// Extracted from G3 specific Config.C 
// by I. Hrivnacova, IPN Orsay

enum PprTrigConf_t {
    kDefaultPPTrig, kDefaultPbPbTrig
};

const char * pprTrigConfName[] = {
    "p-p","Pb-Pb"
};

enum ConfigVersion_t {
    kConfigV0,  // default configuration  
    kConfigV1   // configuration for LHC 2010 production
};

// Options 
static AliMagF::BMap_t smag = AliMagF::k5kG;
static PprTrigConf_t strig = kDefaultPPTrig; // default PP trigger configuration
static TString comment;

// Functions
void  LoadPythia();

void commonConfig(const TString& det, 
                  ConfigVersion_t configVersion = kConfigV0)
{
  cout << "Running commonConfig.C ... " << endl;

    // Set Random Number seed
  gRandom->SetSeed(123456); // Set 0 to use the currecnt time
  AliLog::Message(AliLog::kInfo, Form("Seed for random number generation = %d",gRandom->GetSeed()), "Config.C", "Config.C", "Config()","Config.C", __LINE__);


  //=======================================================================
  // Load Pythia libraries
  //=======================================================================

  LoadPythia();

  //=======================================================================
  // ALICE steering object (AliRunLoader)
  //=======================================================================

  AliRunLoader* rl 
    = AliRunLoader::Open("galice.root",
			  AliConfig::GetDefaultEventFolderName(),
			  "recreate");
  if ( ! rl ) {
    gAlice->Fatal("Config.C","Can not instatiate the Run Loader");
    return;
  }
  rl->SetCompressionLevel(2);
  rl->SetNumberOfEventsPerFile(3);
  gAlice->SetRunLoader(rl);
  
  //======================================================================
  // Trigger configuration
  //=======================================================================

  AliSimulation::Instance()->SetTriggerConfig(pprTrigConfName[strig]);
  cout << "Trigger configuration is set to  " << pprTrigConfName[strig] << endl;

  // ============================= 
  // Magnetic field
  // ============================= 

  // Field (L3 0.5 T)
  AliMagF* field = new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG);
  TGeoGlobalMagField::Instance()->SetField(field);

  printf("\n \n Comment: %s \n \n", comment.Data());

  // ============================= 
  // Modules
  // ============================= 

  rl->CdGAFile();

  Int_t   iABSO  =  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   iHMPID =  0;
  Int_t   iSHIL  =  0;
  Int_t   iT0    =  0;
  Int_t   iTOF   =  0;
  Int_t   iTPC   =  0;
  Int_t   iTRD   =  0;
  Int_t   iZDC   =  0;
  Int_t   iEMCAL =  0;
  Int_t   iACORDE = 0;
  Int_t   iVZERO =  0;

  if ( det == "ABSO" )  iABSO  =  1;
  if ( det == "DIPO" )  iDIPO  =  1;
  if ( det == "FMD" )   iFMD   =  1;
  if ( det == "FRAME" )  iFRAME =  1;
  if ( det == "HALL" )  iHALL  =  1;
  if ( det == "ITS" )   iITS   =  1;
  if ( det == "MAG" )   iMAG   =  1;
  if ( det == "MUON" )  iMUON  =  1;
  if ( det == "PHOS" )  iPHOS  =  1;
  if ( det == "PIPE" )  iPIPE  =  1;
  if ( det == "PMD" )   iPMD   =  1;
  if ( det == "HMPID" )  iHMPID =  1;
  if ( det == "SHIL" )  iSHIL  =  1;
  if ( det == "T0" )    iT0    =  1;
  if ( det == "TOF" )   iTOF   =  1;
  if ( det == "TPC" )   iTPC   =  1;
  if ( det == "TRD" )   iTRD   =  1;
  if ( det == "ZDC" )   iZDC   =  1;
  if ( det == "EMCAL" )  iEMCAL =  1;
  if ( det == "ACORDE" ) iACORDE = 1;
  if ( det == "VZERO" )  iVZERO =  1;

  if ( det == "ALL" ) {
    iABSO  =  1;
    iDIPO  =  1;
    iFMD   =  1;
    iFRAME =  1;
    iHALL  =  1;
    iITS   =  1;
    iMAG   =  1;
    iMUON  =  1;
    iPHOS  =  1;
    iPIPE  =  1;
    iPMD   =  1;
    iHMPID =  1;
    iSHIL  =  1;
    iT0    =  1;
    iTOF   =  1;
    iTPC   =  1;
    iTRD   =  1;
    iZDC   =  1;
    iEMCAL =  1;
    iACORDE = 1;
    iVZERO =  1;
  }

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

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

      AliDIPO *DIPO = new AliDIPOv3("DIPO", "Dipole version 3");
  }

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

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


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

      AliFRAMEv2 *FRAME = new AliFRAMEv2("FRAME", "Space Frame");
      FRAME->SetHoles(1);
  }

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

      AliSHIL *SHIL = new AliSHILv3("SHIL", "Shielding Version 3");
  }


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

      AliPIPE *PIPE = new AliPIPEv3("PIPE", "Beam Pipe");
  }
 
  if (iITS)
  {
      //=================== ITS parameters ============================

      AliITS *ITS  = new AliITSv11("ITS","ITS v11");
  }

  if (iTPC)
  {
      //============================ TPC parameters ===================
      AliTPC *TPC = new AliTPCv2("TPC", "Default");
  }


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


  if (iHMPID)
  {
      //=================== HMPID parameters ===========================
      AliHMPID *HMPID = new AliHMPIDv3("HMPID", "normal HMPID");

  }


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

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

  if (iTRD)
  {
      //=================== TRD parameters ============================

      AliTRD *TRD = new AliTRDv1("TRD", "TRD slow simulator");
      if ( configVersion == kConfigV1 ) {
        AliTRDgeometry *geoTRD = TRD->GetGeometry();
        // Partial geometry: modules at 0,1,7,8,9,16,17
        // starting at 3h in positive direction
        geoTRD->SetSMstatus(2,0);
        geoTRD->SetSMstatus(3,0);
        geoTRD->SetSMstatus(4,0);
        geoTRD->SetSMstatus(5,0);
        geoTRD->SetSMstatus(6,0);
        geoTRD->SetSMstatus(11,0);
        geoTRD->SetSMstatus(12,0);
        geoTRD->SetSMstatus(13,0);
        geoTRD->SetSMstatus(14,0);
        geoTRD->SetSMstatus(15,0);
        geoTRD->SetSMstatus(16,0);
      }
  }

  if (iFMD)
  {
      //=================== FMD parameters ============================
      AliFMD *FMD = new AliFMDv1("FMD", "normal FMD");
  }

  if (iMUON)
  {
      //=================== MUON parameters ===========================
      // New MUONv1 version (geometry defined via builders)
      AliMUON *MUON = new AliMUONv1("MUON", "default");
  }
  //=================== PHOS parameters ===========================

  if (iPHOS)
  {
     if ( configVersion == kConfigV0 ) 
       AliPHOS *PHOS = new AliPHOSv1("PHOS", "IHEP");
     else if ( configVersion == kConfigV1 )  
       AliPHOS *PHOS = new AliPHOSv1("PHOS", "noCPV_Modules123"); 
  }


  if (iPMD)
  {
      //=================== PMD parameters ============================
      AliPMD *PMD = new AliPMDv1("PMD", "normal PMD");
  }

  if (iT0)
  {
      //=================== T0 parameters ============================
      AliT0 *T0 = new AliT0v1("T0", "T0 Detector");
  }

  if (iEMCAL)
  {
      //=================== EMCAL parameters ============================
    if ( configVersion == kConfigV0 ) 
      AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "EMCAL_COMPLETEV1");
    else if ( configVersion == kConfigV1 )  
      AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "EMCAL_FIRSTYEARV1");
  }

   if (iACORDE)
  {
      //=================== ACORDE parameters ============================
      AliACORDE *ACORDE = new AliACORDEv1("ACORDE", "normal ACORDE");
  }

   if (iVZERO)
  {
      //=================== VZERO parameters ============================
      AliVZERO *VZERO = new AliVZEROv7("VZERO", "normal VZERO");
  }

  AliLog::Message(AliLog::kInfo, "End of Config", "Config.C", "Config.C", "Config()"," Config.C", __LINE__);

  cout << "Running commonConfig.C finished ... " << endl;
}

void LoadPythia()
{
  // Load Pythia related libraries
  gSystem->Load("liblhapdf");      // Parton density functions
  gSystem->Load("libEGPythia6");   // TGenerator interface
  gSystem->Load("libpythia6");     // Pythia
  gSystem->Load("libAliPythia6");  // ALICE specific implementations
}
Commit	Line	Data
e2251cfe	1	// $Id$
	2	//
	3	// AliRoot Configuration for running aliroot with Monte Carlo.
	4	// commonConfig() includes the common setting for all MCs
	5	// which has to be called before MC is instantiated.
	6	// Called from MC specific configs (g3Config.C, g4Config.C).
	7	//
	8	// Extracted from G3 specific Config.C
	9	// by I. Hrivnacova, IPN Orsay
	10
	11	enum PprTrigConf_t {
	12	kDefaultPPTrig, kDefaultPbPbTrig
	13	};
	14
	15	const char * pprTrigConfName[] = {
	16	"p-p","Pb-Pb"
	17	};
	18
	19	enum ConfigVersion_t {
	20	kConfigV0, // default configuration
	21	kConfigV1 // configuration for LHC 2010 production
	22	};
	23
	24	// Options
	25	static AliMagF::BMap_t smag = AliMagF::k5kG;
	26	static PprTrigConf_t strig = kDefaultPPTrig; // default PP trigger configuration
	27	static TString comment;
	28
	29	// Functions
	30	void LoadPythia();
	31
	32	void commonConfig(const TString& det,
	33	ConfigVersion_t configVersion = kConfigV0)
	34	{
	35	cout << "Running commonConfig.C ... " << endl;
	36
	37	// Set Random Number seed
	38	gRandom->SetSeed(123456); // Set 0 to use the currecnt time
	39	AliLog::Message(AliLog::kInfo, Form("Seed for random number generation = %d",gRandom->GetSeed()), "Config.C", "Config.C", "Config()","Config.C", __LINE__);
	40
	41
	42	//=======================================================================
	43	// Load Pythia libraries
	44	//=======================================================================
	45
	46	LoadPythia();
	47
	48	//=======================================================================
	49	// ALICE steering object (AliRunLoader)
	50	//=======================================================================
	51
	52	AliRunLoader* rl
	53	= AliRunLoader::Open("galice.root",
	54	AliConfig::GetDefaultEventFolderName(),
	55	"recreate");
	56	if ( ! rl ) {
	57	gAlice->Fatal("Config.C","Can not instatiate the Run Loader");
	58	return;
	59	}
	60	rl->SetCompressionLevel(2);
	61	rl->SetNumberOfEventsPerFile(3);
	62	gAlice->SetRunLoader(rl);
	63
	64	//======================================================================
65	// Trigger configuration
66	//=======================================================================
67
68	AliSimulation::Instance()->SetTriggerConfig(pprTrigConfName[strig]);
69	cout << "Trigger configuration is set to " << pprTrigConfName[strig] << endl;
70
71	// =============================
72	// Magnetic field
73	// =============================
74
75	// Field (L3 0.5 T)
76	AliMagF* field = new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG);
77	TGeoGlobalMagField::Instance()->SetField(field);
78
79	printf("\n \n Comment: %s \n \n", comment.Data());
80
81	// =============================
82	// Modules
83	// =============================
84
85	rl->CdGAFile();
86
87	Int_t iABSO = 0;
88	Int_t iDIPO = 0;
89	Int_t iFMD = 0;
90	Int_t iFRAME = 0;
91	Int_t iHALL = 0;
92	Int_t iITS = 0;
93	Int_t iMAG = 0;
94	Int_t iMUON = 0;
95	Int_t iPHOS = 0;
96	Int_t iPIPE = 0;
97	Int_t iPMD = 0;
98	Int_t iHMPID = 0;
99	Int_t iSHIL = 0;
100	Int_t iT0 = 0;
101	Int_t iTOF = 0;
102	Int_t iTPC = 0;
103	Int_t iTRD = 0;
104	Int_t iZDC = 0;
105	Int_t iEMCAL = 0;
106	Int_t iACORDE = 0;
107	Int_t iVZERO = 0;
108
109	if ( det == "ABSO" ) iABSO = 1;
110	if ( det == "DIPO" ) iDIPO = 1;
111	if ( det == "FMD" ) iFMD = 1;
112	if ( det == "FRAME" ) iFRAME = 1;
113	if ( det == "HALL" ) iHALL = 1;
114	if ( det == "ITS" ) iITS = 1;
115	if ( det == "MAG" ) iMAG = 1;
116	if ( det == "MUON" ) iMUON = 1;
117	if ( det == "PHOS" ) iPHOS = 1;
118	if ( det == "PIPE" ) iPIPE = 1;
119	if ( det == "PMD" ) iPMD = 1;
120	if ( det == "HMPID" ) iHMPID = 1;
121	if ( det == "SHIL" ) iSHIL = 1;
122	if ( det == "T0" ) iT0 = 1;
123	if ( det == "TOF" ) iTOF = 1;
124	if ( det == "TPC" ) iTPC = 1;
125	if ( det == "TRD" ) iTRD = 1;
126	if ( det == "ZDC" ) iZDC = 1;
127	if ( det == "EMCAL" ) iEMCAL = 1;
128	if ( det == "ACORDE" ) iACORDE = 1;
129	if ( det == "VZERO" ) iVZERO = 1;
130
131	if ( det == "ALL" ) {
132	iABSO = 1;
133	iDIPO = 1;
134	iFMD = 1;
135	iFRAME = 1;
136	iHALL = 1;
137	iITS = 1;
138	iMAG = 1;
139	iMUON = 1;
140	iPHOS = 1;
141	iPIPE = 1;
142	iPMD = 1;
143	iHMPID = 1;
144	iSHIL = 1;
145	iT0 = 1;
146	iTOF = 1;
147	iTPC = 1;
148	iTRD = 1;
149	iZDC = 1;
150	iEMCAL = 1;
151	iACORDE = 1;
152	iVZERO = 1;
153	}
154
155	rl->CdGAFile();
156	//=================== Alice BODY parameters =============================
157	AliBODY *BODY = new AliBODY("BODY", "Alice envelop");
158
159	if (iMAG)
160	{
161	//=================== MAG parameters ============================
162	// --- Start with Magnet since detector layouts may be depending ---
163	// --- on the selected Magnet dimensions ---
164	AliMAG *MAG = new AliMAG("MAG", "Magnet");
165	}
166
167
168	if (iABSO)
169	{
170	//=================== ABSO parameters ============================
171	AliABSO *ABSO = new AliABSOv3("ABSO", "Muon Absorber");
172	}
173
174	if (iDIPO)
175	{
176	//=================== DIPO parameters ============================
177
178	AliDIPO *DIPO = new AliDIPOv3("DIPO", "Dipole version 3");
179	}
180
181	if (iHALL)
182	{
183	//=================== HALL parameters ============================
184
185	AliHALL *HALL = new AliHALLv3("HALL", "Alice Hall");
186	}
187
188
189	if (iFRAME)
190	{
191	//=================== FRAME parameters ============================
192
193	AliFRAMEv2 *FRAME = new AliFRAMEv2("FRAME", "Space Frame");
194	FRAME->SetHoles(1);
195	}
196
197	if (iSHIL)
198	{
199	//=================== SHIL parameters ============================
200
201	AliSHIL *SHIL = new AliSHILv3("SHIL", "Shielding Version 3");
202	}
203
204
205	if (iPIPE)
206	{
207	//=================== PIPE parameters ============================
208
209	AliPIPE *PIPE = new AliPIPEv3("PIPE", "Beam Pipe");
210	}
211
212	if (iITS)
213	{
214	//=================== ITS parameters ============================
215
07787d1d	216	AliITS *ITS = new AliITSv11("ITS","ITS v11");
e2251cfe	217	}
	218
	219	if (iTPC)
	220	{
	221	//============================ TPC parameters ===================
	222	AliTPC *TPC = new AliTPCv2("TPC", "Default");
	223	}
	224
	225
	226	if (iTOF) {
	227	//=================== TOF parameters ============================
	228	AliTOF *TOF = new AliTOFv6T0("TOF", "normal TOF");
	229	}
	230
	231
	232	if (iHMPID)
	233	{
	234	//=================== HMPID parameters ===========================
	235	AliHMPID *HMPID = new AliHMPIDv3("HMPID", "normal HMPID");
	236
	237	}
	238
	239
	240	if (iZDC)
	241	{
	242	//=================== ZDC parameters ============================
	243
	244	AliZDC *ZDC = new AliZDCv3("ZDC", "normal ZDC");
	245	}
	246
	247	if (iTRD)
	248	{
	249	//=================== TRD parameters ============================
	250
	251	AliTRD *TRD = new AliTRDv1("TRD", "TRD slow simulator");
	252	if ( configVersion == kConfigV1 ) {
	253	AliTRDgeometry *geoTRD = TRD->GetGeometry();
	254	// Partial geometry: modules at 0,1,7,8,9,16,17
	255	// starting at 3h in positive direction
	256	geoTRD->SetSMstatus(2,0);
	257	geoTRD->SetSMstatus(3,0);
	258	geoTRD->SetSMstatus(4,0);
	259	geoTRD->SetSMstatus(5,0);
	260	geoTRD->SetSMstatus(6,0);
	261	geoTRD->SetSMstatus(11,0);
	262	geoTRD->SetSMstatus(12,0);
	263	geoTRD->SetSMstatus(13,0);
	264	geoTRD->SetSMstatus(14,0);
	265	geoTRD->SetSMstatus(15,0);
	266	geoTRD->SetSMstatus(16,0);
	267	}
	268	}
	269
	270	if (iFMD)
	271	{
	272	//=================== FMD parameters ============================
	273	AliFMD *FMD = new AliFMDv1("FMD", "normal FMD");
	274	}
	275
	276	if (iMUON)
	277	{
	278	//=================== MUON parameters ===========================
	279	// New MUONv1 version (geometry defined via builders)
	280	AliMUON *MUON = new AliMUONv1("MUON", "default");
281	}
282	//=================== PHOS parameters ===========================
283
284	if (iPHOS)
285	{
286	if ( configVersion == kConfigV0 )
287	AliPHOS *PHOS = new AliPHOSv1("PHOS", "IHEP");
288	else if ( configVersion == kConfigV1 )
289	AliPHOS *PHOS = new AliPHOSv1("PHOS", "noCPV_Modules123");
290	}
291
292
293	if (iPMD)
294	{
295	//=================== PMD parameters ============================
296	AliPMD *PMD = new AliPMDv1("PMD", "normal PMD");
297	}
298
299	if (iT0)
300	{
301	//=================== T0 parameters ============================
302	AliT0 *T0 = new AliT0v1("T0", "T0 Detector");
303	}
304
305	if (iEMCAL)
306	{
307	//=================== EMCAL parameters ============================
308	if ( configVersion == kConfigV0 )
309	AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "EMCAL_COMPLETEV1");
310	else if ( configVersion == kConfigV1 )
311	AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "EMCAL_FIRSTYEARV1");
312	}
313
314	if (iACORDE)
315	{
316	//=================== ACORDE parameters ============================
317	AliACORDE *ACORDE = new AliACORDEv1("ACORDE", "normal ACORDE");
318	}
319
320	if (iVZERO)
321	{
322	//=================== VZERO parameters ============================
323	AliVZERO *VZERO = new AliVZEROv7("VZERO", "normal VZERO");
324	}
325
326	AliLog::Message(AliLog::kInfo, "End of Config", "Config.C", "Config.C", "Config()"," Config.C", __LINE__);
327
328	cout << "Running commonConfig.C finished ... " << endl;
329	}
330
331	void LoadPythia()
332	{
333	// Load Pythia related libraries
b0635849	334	gSystem->Load("liblhapdf"); // Parton density functions
	335	gSystem->Load("libEGPythia6"); // TGenerator interface
	336	gSystem->Load("libpythia6"); // Pythia
	337	gSystem->Load("libAliPythia6"); // ALICE specific implementations
e2251cfe	338	}