[u/mrichter/AliRoot.git] / PWGLF / FORWARD / analysis / SubmitFMDCorrections.C

void SubmitFMDCorrections(const Char_t* filename, Bool_t store, Float_t energy, Int_t trigger, Float_t mag, Int_t collsystem,Bool_t inelGtZero=kFALSE) {
  
  gSystem->Load("libANALYSIS");
  gSystem->Load("libANALYSISalice");
  gSystem->Load("libPWG0base");
  gSystem->Load("libPWGLFforward");
  
  AliFMDAnaParameters* pars = AliFMDAnaParameters::Instance();
  
  if(energy == 900)
    pars->SetEnergy(AliFMDAnaParameters::k900);
  else if(energy == 10000)
    pars->SetEnergy(AliFMDAnaParameters::k10000);
  else if(energy == 14000)
    pars->SetEnergy(AliFMDAnaParameters::k14000);
  else if(energy == 7000)
    pars->SetEnergy(AliFMDAnaParameters::k7000);
  else if(energy == 2400)
    pars->SetEnergy(AliFMDAnaParameters::k2400);
  else if(energy == 5500)
    pars->SetEnergy(AliFMDAnaParameters::k5500);
  else if(energy == 2750)
    pars->SetEnergy(AliFMDAnaParameters::k2750);
  
  if(trigger == 0)
    pars->SetTriggerDefinition(AliFMDAnaParameters::kMB1);
  else if(trigger == 1)
    pars->SetTriggerDefinition(AliFMDAnaParameters::kMB2);
  else if(trigger == 2)
    pars->SetTriggerDefinition(AliFMDAnaParameters::kSPDFASTOR);
  else if(trigger == 3)
    pars->SetTriggerDefinition(AliFMDAnaParameters::kNOCTP);
  else if(trigger == 4)
    pars->SetTriggerDefinition(AliFMDAnaParameters::kEMPTY);
  else if(trigger == 5)
    pars->SetTriggerDefinition(AliFMDAnaParameters::kNSD);

  if(mag==0)
    pars->SetMagField(AliFMDAnaParameters::k0G);
  else if(mag==1)
    pars->SetMagField(AliFMDAnaParameters::k5G);
  else if(mag==-1)
    pars->SetMagField(AliFMDAnaParameters::k5Gnegative);
  
  if(collsystem == 0)
    pars->SetCollisionSystem(AliFMDAnaParameters::kPP);
  else if(collsystem == 1)
    pars->SetCollisionSystem(AliFMDAnaParameters::kPbPb);
  
  pars->SetInelGtZero(inelGtZero);
  pars->PrintStatus();
  
  std::cout<<"creating background object"<<std::endl;
  AliFMDAnalysisTaskGenerateCorrection t;
  
  t.ReadFromFile(filename,store,0);
  std::cout<<"object created in backgroundFromFile.root "<<std::flush;
  if(store)
    std::cout<<" - and stored!"<<std::endl;
  else
    std::cout<<" - and not stored!"<<std::endl;
}


//
// EOF
//
Commit	Line	Data
507687cd	1	void SubmitFMDCorrections(const Char_t* filename, Bool_t store, Float_t energy, Int_t trigger, Float_t mag, Int_t collsystem,Bool_t inelGtZero=kFALSE) {
4b8bdb60	2
	3	gSystem->Load("libANALYSIS");
	4	gSystem->Load("libANALYSISalice");
	5	gSystem->Load("libPWG0base");
bd6f5206	6	gSystem->Load("libPWGLFforward");
4b8bdb60	7
	8	AliFMDAnaParameters* pars = AliFMDAnaParameters::Instance();
	9
	10	if(energy == 900)
	11	pars->SetEnergy(AliFMDAnaParameters::k900);
4b8bdb60	12	else if(energy == 10000)
	13	pars->SetEnergy(AliFMDAnaParameters::k10000);
	14	else if(energy == 14000)
	15	pars->SetEnergy(AliFMDAnaParameters::k14000);
efe6cf49	16	else if(energy == 7000)
	17	pars->SetEnergy(AliFMDAnaParameters::k7000);
	18	else if(energy == 2400)
	19	pars->SetEnergy(AliFMDAnaParameters::k2400);
4b8bdb60	20	else if(energy == 5500)
4b8bdb60	21	pars->SetEnergy(AliFMDAnaParameters::k5500);
efe6cf49	22	else if(energy == 2750)
efe6cf49	23	pars->SetEnergy(AliFMDAnaParameters::k2750);
4b8bdb60	24
	25	if(trigger == 0)
	26	pars->SetTriggerDefinition(AliFMDAnaParameters::kMB1);
	27	else if(trigger == 1)
	28	pars->SetTriggerDefinition(AliFMDAnaParameters::kMB2);
efe6cf49	29	else if(trigger == 2)
	30	pars->SetTriggerDefinition(AliFMDAnaParameters::kSPDFASTOR);
	31	else if(trigger == 3)
	32	pars->SetTriggerDefinition(AliFMDAnaParameters::kNOCTP);
	33	else if(trigger == 4)
	34	pars->SetTriggerDefinition(AliFMDAnaParameters::kEMPTY);
	35	else if(trigger == 5)
	36	pars->SetTriggerDefinition(AliFMDAnaParameters::kNSD);
4b8bdb60	37
	38	if(mag==0)
	39	pars->SetMagField(AliFMDAnaParameters::k0G);
	40	else if(mag==1)
	41	pars->SetMagField(AliFMDAnaParameters::k5G);
efe6cf49	42	else if(mag==-1)
efe6cf49	43	pars->SetMagField(AliFMDAnaParameters::k5Gnegative);
4b8bdb60	44
	45	if(collsystem == 0)
	46	pars->SetCollisionSystem(AliFMDAnaParameters::kPP);
	47	else if(collsystem == 1)
	48	pars->SetCollisionSystem(AliFMDAnaParameters::kPbPb);
	49
507687cd	50	pars->SetInelGtZero(inelGtZero);
4b8bdb60	51	pars->PrintStatus();
	52
	53	std::cout<<"creating background object"<<std::endl;
	54	AliFMDAnalysisTaskGenerateCorrection t;
	55
	56	t.ReadFromFile(filename,store,0);
	57	std::cout<<"object created in backgroundFromFile.root "<<std::flush;
	58	if(store)
	59	std::cout<<" - and stored!"<<std::endl;
	60	else
	61	std::cout<<" - and not stored!"<<std::endl;
	62	}
efe6cf49	63
	64
	65	//
	66	// EOF
	67	//
	68