[u/mrichter/AliRoot.git] / JETAN / JetAnalysisManagerLoc.C

void JetAnalysisManagerLoc()
{
      gSystem->Load("libTree.so");
      gSystem->Load("libPhysics.so");
      gSystem->Load("libGeom.so");
      gSystem->Load("libVMC.so");
      gSystem->Load("libANALYSIS.so");
      gSystem->Load("libSTEERBase.so");
      gSystem->Load("libAOD.so");
      gSystem->Load("libESD.so");
      gSystem->Load("libANALYSISalice.so");
      gSystem->Load("libJETAN.so");

     //
    if (gApplication) gApplication->InitializeGraphics();
    // Create the chain
    //
    gROOT->LoadMacro("CreateESDChain.C");
    TChain* chain = new TChain("esdTree");
    chain->Add("~/alice/data/highpt/kPythia6Jets125_150/030/AliESDs.root");

    /////////////////////////////////////////////////////////////////////////////////// 
    // Create the analysis manager
    //
    // Input 
    AliESDInputHandler* inpHandler = new AliESDInputHandler();
    // Output
    AliAODHandler* aodHandler = new AliAODHandler();
    aodHandler->SetOutputFileName("aod.root");
    // MC Truth
    AliMCEventHandler* mcHandler = new AliMCEventHandler();
    

    AliAnalysisManager *mgr  = new AliAnalysisManager("Jet Manager", "Jet Manager");
    mgr->SetInputEventHandler  (inpHandler);
    mgr->SetOutputEventHandler (aodHandler);
    mgr->SetMCtruthEventHandler(mcHandler);
    mgr->SetDebugLevel(10);
    /////////////////////////////////////////////////////////////////////////////////// 
    
    
    //
    // Set of cuts
    // 
    // standard
    AliESDtrackCuts* esdTrackCutsL = new AliESDtrackCuts("AliESDtrackCuts", "Loose");
    esdTrackCutsL->SetMinNClustersTPC(50);
    esdTrackCutsL->SetMaxChi2PerClusterTPC(3.5);
    esdTrackCutsL->SetMaxCovDiagonalElements(2,2,0.5,0.5,2);
    esdTrackCutsL->SetRequireTPCRefit(kTRUE);
    esdTrackCutsL->SetMinNsigmaToVertex(3);
    esdTrackCutsL->SetRequireSigmaToVertex(kTRUE);
    esdTrackCutsL->SetAcceptKingDaughters(kFALSE);
    //
    // hard
    AliESDtrackCuts* esdTrackCutsH = new AliESDtrackCuts("AliESDtrackCuts", "Hard");
    esdTrackCutsH->SetMinNClustersTPC(100);
    esdTrackCutsH->SetMaxChi2PerClusterTPC(2.0);
    esdTrackCutsH->SetMaxCovDiagonalElements(2,2,0.5,0.5,2);
    esdTrackCutsH->SetRequireTPCRefit(kTRUE);
    esdTrackCutsH->SetMinNsigmaToVertex(2);
    esdTrackCutsH->SetRequireSigmaToVertex(kTRUE);
    esdTrackCutsH->SetAcceptKingDaughters(kFALSE);
    //
    AliAnalysisFilter* trackFilter = new AliAnalysisFilter("trackFilter");
    trackFilter->AddCuts(esdTrackCutsL);
    trackFilter->AddCuts(esdTrackCutsH);
    //
    AliAnalysisTaskESDfilter *esdfilter = new AliAnalysisTaskESDfilter("ESD Filter");
    esdfilter->SetTrackFilter(trackFilter);
    esdfilter->SetDebugLevel(10);
    mgr->AddTask(esdfilter);

    AliAnalysisTaskJets *jetana = new AliAnalysisTaskJets("JetAnalysis");
    jetana->SetDebugLevel(10);


    AliAnalysisTaskJets *jetanaMC = new AliAnalysisTaskJets("JetAnalysisMC");
    jetanaMC->SetDebugLevel(10);
    jetanaMC->SetConfigFile("ConfigJetAnalysisMC.C");
    jetanaMC->SetNonStdBranch("jetsMC");
    mgr->AddTask(jetanaMC);
    mgr->AddTask(jetana);

    //
    // Create containers for input/output
    AliAnalysisDataContainer *cinput1 = mgr->CreateContainer("cchain",TChain::Class(), 
							     AliAnalysisManager::kInputContainer);

    AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("tree", TTree::Class(),
							      AliAnalysisManager::kOutputContainer, "default");
    AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("histos", TList::Class(),
							      AliAnalysisManager::kOutputContainer, "histos.root");

    AliAnalysisDataContainer *coutputMC2 = mgr->CreateContainer("histosMC", TList::Class(),
							      AliAnalysisManager::kOutputContainer, "histosMC.root");

    mgr->ConnectInput  (esdfilter,  0, cinput1  );
    mgr->ConnectOutput (esdfilter,  0, coutput1 );

    mgr->ConnectInput  (jetana,     0, cinput1  );
    mgr->ConnectOutput (jetana,     0, coutput1 );
    mgr->ConnectOutput (jetana,     1, coutput2 );

    mgr->ConnectInput  (jetanaMC,     0, cinput1  );
    mgr->ConnectOutput (jetanaMC,     0, coutput1 );
    mgr->ConnectOutput (jetanaMC,     1, coutputMC2 );


    //
    // Run the analysis
    //    
    mgr->InitAnalysis();
    mgr->PrintStatus();
    mgr->StartAnalysis("local",chain);
}
Commit	Line	Data
	1	void JetAnalysisManagerLoc()
	2	{
	3	gSystem->Load("libTree.so");
	4	gSystem->Load("libPhysics.so");
	5	gSystem->Load("libGeom.so");
	6	gSystem->Load("libVMC.so");
	7	gSystem->Load("libANALYSIS.so");
	8	gSystem->Load("libSTEERBase.so");
	9	gSystem->Load("libAOD.so");
	10	gSystem->Load("libESD.so");
	11	gSystem->Load("libANALYSISalice.so");
	12	gSystem->Load("libJETAN.so");
	13
	14	//
	15	if (gApplication) gApplication->InitializeGraphics();
	16	// Create the chain
	17	//
	18	gROOT->LoadMacro("CreateESDChain.C");
	19	TChain* chain = new TChain("esdTree");
	20	chain->Add("~/alice/data/highpt/kPythia6Jets125_150/030/AliESDs.root");
	21
	22	///////////////////////////////////////////////////////////////////////////////////
	23	// Create the analysis manager
	24	//
	25	// Input
	26	AliESDInputHandler* inpHandler = new AliESDInputHandler();
	27	// Output
	28	AliAODHandler* aodHandler = new AliAODHandler();
	29	aodHandler->SetOutputFileName("aod.root");
	30	// MC Truth
	31	AliMCEventHandler* mcHandler = new AliMCEventHandler();
	32
	33
	34	AliAnalysisManager *mgr = new AliAnalysisManager("Jet Manager", "Jet Manager");
	35	mgr->SetInputEventHandler (inpHandler);
	36	mgr->SetOutputEventHandler (aodHandler);
	37	mgr->SetMCtruthEventHandler(mcHandler);
	38	mgr->SetDebugLevel(10);
	39	///////////////////////////////////////////////////////////////////////////////////
	40
	41
	42	//
	43	// Set of cuts
	44	//
	45	// standard
	46	AliESDtrackCuts* esdTrackCutsL = new AliESDtrackCuts("AliESDtrackCuts", "Loose");
	47	esdTrackCutsL->SetMinNClustersTPC(50);
	48	esdTrackCutsL->SetMaxChi2PerClusterTPC(3.5);
	49	esdTrackCutsL->SetMaxCovDiagonalElements(2,2,0.5,0.5,2);
	50	esdTrackCutsL->SetRequireTPCRefit(kTRUE);
	51	esdTrackCutsL->SetMinNsigmaToVertex(3);
	52	esdTrackCutsL->SetRequireSigmaToVertex(kTRUE);
	53	esdTrackCutsL->SetAcceptKingDaughters(kFALSE);
	54	//
	55	// hard
	56	AliESDtrackCuts* esdTrackCutsH = new AliESDtrackCuts("AliESDtrackCuts", "Hard");
	57	esdTrackCutsH->SetMinNClustersTPC(100);
	58	esdTrackCutsH->SetMaxChi2PerClusterTPC(2.0);
	59	esdTrackCutsH->SetMaxCovDiagonalElements(2,2,0.5,0.5,2);
	60	esdTrackCutsH->SetRequireTPCRefit(kTRUE);
	61	esdTrackCutsH->SetMinNsigmaToVertex(2);
	62	esdTrackCutsH->SetRequireSigmaToVertex(kTRUE);
	63	esdTrackCutsH->SetAcceptKingDaughters(kFALSE);
	64	//
	65	AliAnalysisFilter* trackFilter = new AliAnalysisFilter("trackFilter");
	66	trackFilter->AddCuts(esdTrackCutsL);
	67	trackFilter->AddCuts(esdTrackCutsH);
	68	//
	69	AliAnalysisTaskESDfilter *esdfilter = new AliAnalysisTaskESDfilter("ESD Filter");
	70	esdfilter->SetTrackFilter(trackFilter);
	71	esdfilter->SetDebugLevel(10);
	72	mgr->AddTask(esdfilter);
	73
	74	AliAnalysisTaskJets *jetana = new AliAnalysisTaskJets("JetAnalysis");
	75	jetana->SetDebugLevel(10);
	76
	77
	78
	79	AliAnalysisTaskJets *jetanaMC = new AliAnalysisTaskJets("JetAnalysisMC");
	80	jetanaMC->SetDebugLevel(10);
	81	jetanaMC->SetConfigFile("ConfigJetAnalysisMC.C");
	82	jetanaMC->SetNonStdBranch("jetsMC");
	83	mgr->AddTask(jetanaMC);
	84	mgr->AddTask(jetana);
	85
	86	//
	87	// Create containers for input/output
	88	AliAnalysisDataContainer *cinput1 = mgr->CreateContainer("cchain",TChain::Class(),
	89	AliAnalysisManager::kInputContainer);
	90
	91	AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("tree", TTree::Class(),
	92	AliAnalysisManager::kOutputContainer, "default");
	93	AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("histos", TList::Class(),
	94	AliAnalysisManager::kOutputContainer, "histos.root");
	95
	96	AliAnalysisDataContainer *coutputMC2 = mgr->CreateContainer("histosMC", TList::Class(),
	97	AliAnalysisManager::kOutputContainer, "histosMC.root");
	98
	99	mgr->ConnectInput (esdfilter, 0, cinput1 );
	100	mgr->ConnectOutput (esdfilter, 0, coutput1 );
	101
	102	mgr->ConnectInput (jetana, 0, cinput1 );
	103	mgr->ConnectOutput (jetana, 0, coutput1 );
	104	mgr->ConnectOutput (jetana, 1, coutput2 );
	105
	106	mgr->ConnectInput (jetanaMC, 0, cinput1 );
	107	mgr->ConnectOutput (jetanaMC, 0, coutput1 );
	108	mgr->ConnectOutput (jetanaMC, 1, coutputMC2 );
	109
	110
	111	//
	112	// Run the analysis
	113	//
	114	mgr->InitAnalysis();
	115	mgr->PrintStatus();
	116	mgr->StartAnalysis("local",chain);
	117	}