[u/mrichter/AliRoot.git] / GRP / UpdateCDBGRPEntryMC.C

  void UpdateCDBGRPEntryMC() {
	// produce the GRP entry in CDB 
	// reading MC parameter from Config.C
	
	AliCDBManager* man = AliCDBManager::Instance();
  man->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
  man->SetRun(0);
	AliCDBId id(AliQA::GetGRPPath(), 0, AliCDBRunRange::Infinity());
  AliCDBMetaData *metadata= new AliCDBMetaData();

  // Get root version
  const char* rootv = gROOT->GetVersion();
	metadata->SetResponsible("Yves.Schutz@cern.ch");
	metadata->SetComment("GRP parameters for Monte Carlo");
 
		
	TMap *mappp = GetGRPList();
		
  //  Printf("Storing in CDB the default values for the GRP %d parameters produced with root %s and AliRoot version %s",list->GetEntries(),rootv,alirootv);
	
  man->Put(mappp,id,metadata);
}

TString ParseConfig(char * option)
{
	// Parse Config file to retrieve usefull information for the GRP
	
	ifstream in("Config.C", ios::in) ; 
  if (!in)
		AliError("Config.C file not found in current path") ; 	

	TString rv("") ; 
	
	TString soption(option) ;
	
	string line ; 
	TString sline ; 

	// search for the system p-p or Pb-Pb or ..... 
	if (soption.Contains("system") || soption.Contains("fAliceBeamType")){
		while (getline(in, line)) { 
			sline = line ;		
			if (sline.Contains("AliSimulation::Instance()->SetTriggerConfig")) 
				break ;
		}
		TString sarg1 = sline(sline.Index("[")+1, sline.Index("]") - sline.Index("[")-1) ; 
		sarg1 = sarg1.Strip(TString::kLeading) ; 
		sarg1 = sarg1.Strip(TString::kTrailing) ; 
		TString sarg2 = sline(sline.Index("(")+1, sline.Index("[") - sline.Index("(")-1) ; 
		sarg2 = sarg2.Strip(TString::kLeading) ; 
		sarg2 = sarg2.Strip(TString::kTrailing) ; 
		in.seekg(0) ;
		while (getline(in, line)) {
			sline = line ;
			if (sline.Contains(sarg1)) {
				sarg1 = sline(sline.Index("=")+1, sline.Index(";")-sline.Index("=")-1) ; 
				sarg1 = sarg1.Strip(TString::kLeading) ; 
				sarg1 = sarg1.Strip(TString::kTrailing) ; 
				break ; 
			}
		}
		in.seekg(0) ;
		Int_t index = -1 ; 
		while (getline(in, line)) {
			sline = line ;
			if (sline.Contains(sarg1)) {
				sline = sline.Strip(TString::kLeading) ; 				
				Ssiz_t deb = sline.Index(sarg1) ;
				Ssiz_t kom = sline.Index(",") ; 
				if ( deb == 0 || deb < kom)
					index = 1 ; 
				else 
					index = 2 ; 
				continue ; 
			}
			if (sline.Contains(sarg2)) {
				Ssiz_t deb = sline.First("\"") ; 
				if ( deb == -1 ) {
					getline(in, line) ; 
					sline = line ; 
				}
				sline.ReplaceAll("\"", "") ;
				sline = sline.Strip(TString::kLeading) ; 
				if (index == 1 ) 
					rv = sline(0, sline.Index(",")) ;
				else if (index == 2) 
					rv = sline(sline.Index(",")+1, sline.Length()-sline.Index(",")) ; 
				break ;
			}
		}
	} 
	// get the beam energy
	else if (soption.Contains("fAliceBeamEnergy")) {
		while (getline(in, line)) { 
			sline = line ;		
			if ( sline.Contains("->SetEnergyCMS(") ) {
				rv = sline(sline.Index("(")+1, sline.Index(")")-sline.Index("(")-1) ; 
				break ;		
			}
		}
	}
	// get the number of active detectors
	else if (soption.Contains("fNumberOfDetectors")) {
		UShort_t det = 0 ; 
		while (getline(in, line)) {
			sline = line ; 
			sline.ReplaceAll(" ", "") ; 
			if (sline.Contains("iITS=1") ||
					sline.Contains("iTPC=1") ||
					sline.Contains("iTRD=1") ||
					sline.Contains("iTOF=1") ||			
					sline.Contains("iPHOS=1") ||
					sline.Contains("iHMPID=1") ||
					sline.Contains("iEMCAL=1") ||
					sline.Contains("iMUON=1") ||
					sline.Contains("iFMD=1") ||
					sline.Contains("iPMD=1") ||
					sline.Contains("iT0=1") ||
					sline.Contains("iVZERO=1") ||
					sline.Contains("iZDC=1") ||
					sline.Contains("iACORDE=1"))
				det++ ; 
		}
		rv = Form("%d", det) ; 
	}
	else if (soption.Contains("fL3Current")) {
		TString sarg ;
		while (getline(in, line)) {
			sline = line ; 
			sline.ReplaceAll(" ", "") ; 
			if (sline.Contains("AliMagF*field=newAliMagF(")) {
				sarg = sline(sline.Last(',')+1, sline.Last(')')-sline.Last(',')-1) ; 
				break ; 
			}
		}
		in.seekg(0) ;
		while (getline(in, line)) {
			sline = line ; 
			sline.ReplaceAll(" ", "") ; 
			if (sline.Contains(Form("%s=",sarg.Data()))) {
				sarg = sline(sline.Index("=")+1, sline.Index(";")-sline.Index("=")-1) ;
				if (sarg == "k5kG") 
					rv = "30000" ;
				else if (sarg == "k2kG")
					rv = "12000" ;
				break ; 
			}
		}		
	}
		in.close() ;
		rv.Strip(TString::kLeading) ; 
		rv.Strip() ; 
	return rv ; 
}
	
//_______________________________________//
TMap *GetGRPList() {
 
  TString system = ParseConfig("fAliceBeamType") ;  
	TString fSystem = system;
  TMap *map = new TMap();
  map->SetName("MONTECARLO");

  //DAQ
	map->Add(new TObjString("fRunType"),new TObjString(AliQA::GetRunTypeName(AliQA::kPHYSICS)));
	map->Add(new TObjString("fAliceStartTime"),new TObjString("0"));
  map->Add(new TObjString("fAliceStopTime"),new TObjString("9999"));
	map->Add(new TObjString("fAliceBeamEnergy"),new TObjString(ParseConfig("fAliceBeamEnergy")));
  map->Add(new TObjString("fAliceBeamType"),new TObjString(system));
  map->Add(new TObjString("fNumberOfDetectors"),new TObjString(ParseConfig("fNumberOfDetectors")));
  map->Add(new TObjString("fDetectorMask"),new TObjString("1074790399"));
  map->Add(new TObjString("fLHCPeriod"),new TObjString("LHC08c"));

  //DCS
  map->Add(new TObjString("fLHCState"),new TObjString("0"));
  map->Add(new TObjString("fLHCCondition"),new TObjString("0"));
  map->Add(new TObjString("fLHCLuminosity"),new TObjString("0"));
  map->Add(new TObjString("fBeamIntensity"),new TObjString("0"));
  map->Add(new TObjString("fL3Current"),new TObjString(ParseConfig("fL3Current")));
  map->Add(new TObjString("fL3Polarity"),new TObjString("0"));
  map->Add(new TObjString("fDipoleCurrent"),new TObjString("6000"));
  map->Add(new TObjString("fDipolePolarity"),new TObjString("0"));
  map->Add(new TObjString("fCavernTemperature"),new TObjString("0"));
  map->Add(new TObjString("fCavernPressure"),new TObjString("0"));

  return map;
}
Commit	Line	Data
75958532	1	void UpdateCDBGRPEntryMC() {
	2	// produce the GRP entry in CDB
	3	// reading MC parameter from Config.C
	4
	5	AliCDBManager* man = AliCDBManager::Instance();
162637e4	6	man->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
75958532	7	man->SetRun(0);
	8	AliCDBId id(AliQA::GetGRPPath(), 0, AliCDBRunRange::Infinity());
	9	AliCDBMetaData *metadata= new AliCDBMetaData();
	10
	11	// Get root version
	12	const char* rootv = gROOT->GetVersion();
	13	metadata->SetResponsible("Yves.Schutz@cern.ch");
	14	metadata->SetComment("GRP parameters for Monte Carlo");
	15
	16
	17	TMap *mappp = GetGRPList();
	18
	19	// Printf("Storing in CDB the default values for the GRP %d parameters produced with root %s and AliRoot version %s",list->GetEntries(),rootv,alirootv);
	20
	21	man->Put(mappp,id,metadata);
	22	}
	23
	24	TString ParseConfig(char * option)
	25	{
	26	// Parse Config file to retrieve usefull information for the GRP
	27
	28	ifstream in("Config.C", ios::in) ;
	29	if (!in)
	30	AliError("Config.C file not found in current path") ;
	31
	32	TString rv("") ;
	33
	34	TString soption(option) ;
	35
	36	string line ;
	37	TString sline ;
	38
	39	// search for the system p-p or Pb-Pb or .....
	40	if (soption.Contains("system") \|\| soption.Contains("fAliceBeamType")){
	41	while (getline(in, line)) {
	42	sline = line ;
28da60d3	43	if (sline.Contains("AliSimulation::Instance()->SetTriggerConfig"))
75958532	44	break ;
	45	}
	46	TString sarg1 = sline(sline.Index("[")+1, sline.Index("]") - sline.Index("[")-1) ;
	47	sarg1 = sarg1.Strip(TString::kLeading) ;
	48	sarg1 = sarg1.Strip(TString::kTrailing) ;
	49	TString sarg2 = sline(sline.Index("(")+1, sline.Index("[") - sline.Index("(")-1) ;
	50	sarg2 = sarg2.Strip(TString::kLeading) ;
	51	sarg2 = sarg2.Strip(TString::kTrailing) ;
	52	in.seekg(0) ;
	53	while (getline(in, line)) {
	54	sline = line ;
	55	if (sline.Contains(sarg1)) {
	56	sarg1 = sline(sline.Index("=")+1, sline.Index(";")-sline.Index("=")-1) ;
	57	sarg1 = sarg1.Strip(TString::kLeading) ;
	58	sarg1 = sarg1.Strip(TString::kTrailing) ;
	59	break ;
	60	}
	61	}
	62	in.seekg(0) ;
	63	Int_t index = -1 ;
	64	while (getline(in, line)) {
	65	sline = line ;
	66	if (sline.Contains(sarg1)) {
	67	sline = sline.Strip(TString::kLeading) ;
	68	Ssiz_t deb = sline.Index(sarg1) ;
	69	Ssiz_t kom = sline.Index(",") ;
	70	if ( deb == 0 \|\| deb < kom)
	71	index = 1 ;
	72	else
	73	index = 2 ;
	74	continue ;
	75	}
	76	if (sline.Contains(sarg2)) {
	77	Ssiz_t deb = sline.First("\"") ;
	78	if ( deb == -1 ) {
	79	getline(in, line) ;
	80	sline = line ;
	81	}
	82	sline.ReplaceAll("\"", "") ;
	83	sline = sline.Strip(TString::kLeading) ;
	84	if (index == 1 )
	85	rv = sline(0, sline.Index(",")) ;
	86	else if (index == 2)
	87	rv = sline(sline.Index(",")+1, sline.Length()-sline.Index(",")) ;
	88	break ;
	89	}
	90	}
	91	}
	92	// get the beam energy
	93	else if (soption.Contains("fAliceBeamEnergy")) {
	94	while (getline(in, line)) {
	95	sline = line ;
	96	if ( sline.Contains("->SetEnergyCMS(") ) {
	97	rv = sline(sline.Index("(")+1, sline.Index(")")-sline.Index("(")-1) ;
	98	break ;
	99	}
	100	}
	101	}
	102	// get the number of active detectors
	103	else if (soption.Contains("fNumberOfDetectors")) {
	104	UShort_t det = 0 ;
	105	while (getline(in, line)) {
	106	sline = line ;
	107	sline.ReplaceAll(" ", "") ;
108	if (sline.Contains("iITS=1") \|\|
109	sline.Contains("iTPC=1") \|\|
110	sline.Contains("iTRD=1") \|\|
111	sline.Contains("iTOF=1") \|\|
112	sline.Contains("iPHOS=1") \|\|
113	sline.Contains("iHMPID=1") \|\|
114	sline.Contains("iEMCAL=1") \|\|
115	sline.Contains("iMUON=1") \|\|
116	sline.Contains("iFMD=1") \|\|
117	sline.Contains("iPMD=1") \|\|
118	sline.Contains("iT0=1") \|\|
119	sline.Contains("iVZERO=1") \|\|
120	sline.Contains("iZDC=1") \|\|
121	sline.Contains("iACORDE=1"))
122	det++ ;
123	}
124	rv = Form("%d", det) ;
125	}
126	else if (soption.Contains("fL3Current")) {
127	TString sarg ;
128	while (getline(in, line)) {
129	sline = line ;
130	sline.ReplaceAll(" ", "") ;
f7a1cc68	131	if (sline.Contains("AliMagF*field=newAliMagF(")) {
75958532	132	sarg = sline(sline.Last(',')+1, sline.Last(')')-sline.Last(',')-1) ;
	133	break ;
	134	}
	135	}
	136	in.seekg(0) ;
	137	while (getline(in, line)) {
	138	sline = line ;
	139	sline.ReplaceAll(" ", "") ;
	140	if (sline.Contains(Form("%s=",sarg.Data()))) {
	141	sarg = sline(sline.Index("=")+1, sline.Index(";")-sline.Index("=")-1) ;
	142	if (sarg == "k5kG")
	143	rv = "30000" ;
	144	else if (sarg == "k2kG")
	145	rv = "12000" ;
75958532	146	break ;
	147	}
	148	}
	149	}
	150	in.close() ;
	151	rv.Strip(TString::kLeading) ;
	152	rv.Strip() ;
	153	return rv ;
	154	}
	155
	156	//_______________________________________//
	157	TMap *GetGRPList() {
	158
	159	TString system = ParseConfig("fAliceBeamType") ;
	160	TString fSystem = system;
	161	TMap *map = new TMap();
	162	map->SetName("MONTECARLO");
	163
	164	//DAQ
	165	map->Add(new TObjString("fRunType"),new TObjString(AliQA::GetRunTypeName(AliQA::kPHYSICS)));
	166	map->Add(new TObjString("fAliceStartTime"),new TObjString("0"));
	167	map->Add(new TObjString("fAliceStopTime"),new TObjString("9999"));
	168	map->Add(new TObjString("fAliceBeamEnergy"),new TObjString(ParseConfig("fAliceBeamEnergy")));
	169	map->Add(new TObjString("fAliceBeamType"),new TObjString(system));
	170	map->Add(new TObjString("fNumberOfDetectors"),new TObjString(ParseConfig("fNumberOfDetectors")));
7e88424f	171	map->Add(new TObjString("fDetectorMask"),new TObjString("1074790399"));
7e88424f	172	map->Add(new TObjString("fLHCPeriod"),new TObjString("LHC08c"));
75958532	173
	174	//DCS
	175	map->Add(new TObjString("fLHCState"),new TObjString("0"));
	176	map->Add(new TObjString("fLHCCondition"),new TObjString("0"));
	177	map->Add(new TObjString("fLHCLuminosity"),new TObjString("0"));
	178	map->Add(new TObjString("fBeamIntensity"),new TObjString("0"));
	179	map->Add(new TObjString("fL3Current"),new TObjString(ParseConfig("fL3Current")));
	180	map->Add(new TObjString("fL3Polarity"),new TObjString("0"));
	181	map->Add(new TObjString("fDipoleCurrent"),new TObjString("6000"));
	182	map->Add(new TObjString("fDipolePolarity"),new TObjString("0"));
	183	map->Add(new TObjString("fCavernTemperature"),new TObjString("0"));
	184	map->Add(new TObjString("fCavernPressure"),new TObjString("0"));
	185
	186	return map;
	187	}