[u/mrichter/AliRoot.git] / PWGCF / EBYE / PIDFluctuation / task / SteerAODAnalysisTaskPIDFluctuation.C

SteerAODAnalysisTaskPIDFluctuation(const Char_t *inputfilename, Int_t maxFiles = kMaxInt, Int_t maxEv = kMaxInt)
{

  /* include path for ACLic */
  gSystem->AddIncludePath("-I$ALICE_ROOT/include");
  gSystem->AddIncludePath("-I$ALICE_ROOT/TOF");
  /* load libraries */
  gSystem->Load("libANALYSIS");
  gSystem->Load("libANALYSISalice");
  /* build analysis task class */
  gROOT->LoadMacro("AliAnalysisTaskPIDFluctuation.cxx+g");

  /* setup input chain */
  TString str = inputfilename;
  const Char_t *filename;
  TChain *chain = new TChain("aodTree");
  if (str.EndsWith(".xml")) {
    TGrid::Connect("alien://");
    Info("", "reading data list from collection:");
    TAlienCollection coll(inputfilename, maxFiles);
    coll.Reset();
    while (coll.Next()) {
      filename = coll.GetTURL();
      Info("", Form("%s", filename));
      chain->Add(filename);
    }
  }
  else if (str.EndsWith(".txt")) {
    Info("", "reading data list from text file:");
    ifstream is(inputfilename);
    Char_t buf[4096];
    while(!is.eof()) {
      is.getline(buf, 4096);
      if (is.eof()) break;
      chain->Add(buf);
      Info("", Form("%s", buf));
    }
    is.close();
  }
  else {
    Info("", "single file:");
    filename = inputfilename;
    Info("", Form("%s", filename));
    chain->Add(filename);
  }
  Info("", Form("chain is ready: %d events", chain->GetEntries()));

  /* create analysis manager */
  AliAnalysisManager *mgr = new AliAnalysisManager("PIDFluctuation");

  /* define input event handler */
  AliAODInputHandler *aodh = new AliAODInputHandler();
  mgr->SetInputEventHandler(aodh);

  gROOT->LoadMacro("AddAnalysisTaskPIDFluctuation.C");
  AddAnalysisTaskPIDFluctuation(AliAODTrack::kTrkGlobal, 0., 10., -0.8, 0.8);
  AddAnalysisTaskPIDFluctuation(AliAODTrack::kTrkTPCOnly, 0., 10., -0.8, 0.8);
  AddAnalysisTaskPIDFluctuation(AliAODTrack::kTrkTPCOnly, 0.3, 1.5, -0.8, 0.8);
  AddAnalysisTaskPIDFluctuation(AliAODTrack::kTrkTPCOnly, 0.3, 1.5, -0.4, 0.4);

  /* start analysis */
  mgr->SetDebugLevel(0);
  if (!mgr->InitAnalysis()) return;
  mgr->PrintStatus();
  mgr->StartAnalysis("local", chain, maxEv);

  /* create dummy file to tell we are done */
  gSystem->Exec("touch done");

}
Commit	Line	Data
a6fd3cfe	1	SteerAODAnalysisTaskPIDFluctuation(const Char_t *inputfilename, Int_t maxFiles = kMaxInt, Int_t maxEv = kMaxInt)
	2	{
	3
	4	/* include path for ACLic */
	5	gSystem->AddIncludePath("-I$ALICE_ROOT/include");
	6	gSystem->AddIncludePath("-I$ALICE_ROOT/TOF");
	7	/* load libraries */
	8	gSystem->Load("libANALYSIS");
	9	gSystem->Load("libANALYSISalice");
	10	/* build analysis task class */
	11	gROOT->LoadMacro("AliAnalysisTaskPIDFluctuation.cxx+g");
	12
	13	/* setup input chain */
	14	TString str = inputfilename;
	15	const Char_t *filename;
	16	TChain *chain = new TChain("aodTree");
	17	if (str.EndsWith(".xml")) {
	18	TGrid::Connect("alien://");
	19	Info("", "reading data list from collection:");
	20	TAlienCollection coll(inputfilename, maxFiles);
	21	coll.Reset();
	22	while (coll.Next()) {
	23	filename = coll.GetTURL();
	24	Info("", Form("%s", filename));
	25	chain->Add(filename);
	26	}
	27	}
	28	else if (str.EndsWith(".txt")) {
	29	Info("", "reading data list from text file:");
	30	ifstream is(inputfilename);
	31	Char_t buf[4096];
	32	while(!is.eof()) {
	33	is.getline(buf, 4096);
	34	if (is.eof()) break;
	35	chain->Add(buf);
	36	Info("", Form("%s", buf));
	37	}
	38	is.close();
	39	}
	40	else {
	41	Info("", "single file:");
	42	filename = inputfilename;
	43	Info("", Form("%s", filename));
	44	chain->Add(filename);
	45	}
	46	Info("", Form("chain is ready: %d events", chain->GetEntries()));
	47
	48	/* create analysis manager */
	49	AliAnalysisManager *mgr = new AliAnalysisManager("PIDFluctuation");
	50
	51	/* define input event handler */
	52	AliAODInputHandler *aodh = new AliAODInputHandler();
	53	mgr->SetInputEventHandler(aodh);
	54
	55	gROOT->LoadMacro("AddAnalysisTaskPIDFluctuation.C");
	56	AddAnalysisTaskPIDFluctuation(AliAODTrack::kTrkGlobal, 0., 10., -0.8, 0.8);
	57	AddAnalysisTaskPIDFluctuation(AliAODTrack::kTrkTPCOnly, 0., 10., -0.8, 0.8);
	58	AddAnalysisTaskPIDFluctuation(AliAODTrack::kTrkTPCOnly, 0.3, 1.5, -0.8, 0.8);
	59	AddAnalysisTaskPIDFluctuation(AliAODTrack::kTrkTPCOnly, 0.3, 1.5, -0.4, 0.4);
	60
	61	/* start analysis */
	62	mgr->SetDebugLevel(0);
	63	if (!mgr->InitAnalysis()) return;
	64	mgr->PrintStatus();
65	mgr->StartAnalysis("local", chain, maxEv);
66
67	/* create dummy file to tell we are done */
68	gSystem->Exec("touch done");
69
70	}