gSystem->Load("libCORRFW.so");
gSystem->Load("libPHOSUtils.so");
gSystem->Load("libEMCALUtils.so");
- gSystem->Load("libPWG4CaloCalib.so");
+ gSystem->Load("libPWGGACaloTasks.so");
- /*
- // gSystem->Load("libPWG4omega3pi.so");
- // gSystem->Load("libCORRFW.so");
- // gSystem->Load("libPWGCFbase.so");
- // gSystem->Load("libPWGmuon.so");
- */
//--------------------------------------------------------
//If you want to use root and par files from aliroot
//--------------------------------------------------------
SetupPar("EMCALUtils");
//Create Geometry
TGeoManager::Import("geometry.root") ; //need file "geometry.root" in local dir!!!!
- SetupPar("PWG4CaloCalib");
+ SetupPar("PWGGACaloTasks");
*/
}
// gProof->ClearPackage("ESD");
// gProof->ClearPackage("AOD");
// gProof->ClearPackage("ANALYSIS");
- // gProof->ClearPackage("PWG4PartCorrBase");
- // gProof->ClearPackage("PWG4PartCorrDep");
-
+
// Enable the STEERBase Package
gProof->UploadPackage("STEERBase.par");
gProof->EnablePackage("STEERBase");
// Enable the PHOS geometry Package
//gProof->UploadPackage("PHOSUtils.par");
//gProof->EnablePackage("PHOSUtils");
- // Enable PartCorr analysis
- gProof->UploadPackage("PWG4PartCorrBase.par");
- gProof->EnablePackage("PWG4PartCorrBase");
- gProof->UploadPackage("PWG4PartCorrDep.par");
- gProof->EnablePackage("PWG4PartCorrDep");
gProof->ShowEnabledPackages();
}
if(kInputData=="ESD"){
printf("* Configure photon conversion analysis in macro \n");
TString arguments = "-run-on-train -use-own-xyz -force-aod -mc-off ";
- gROOT->LoadMacro("$ALICE_ROOT/PWG4/macros/ConfigGammaConversion.C");
+ gROOT->LoadMacro("$ALICE_ROOT/PWGGA/GammaConv/macros/ConfigGammaConversion.C");
AliAnalysisTaskGammaConversion * taskGammaConversion =
ConfigGammaConversion(arguments,mgr->GetCommonInputContainer());
taskGammaConversion->SelectCollisionCandidates();
Bool_t kPrint = kFALSE;
Bool_t deltaAOD = kFALSE;
- gROOT->LoadMacro("AddTaskCaloTrackCorr.C"); // $ALICE_ROOT/PWG4/macros
- gROOT->LoadMacro("AddTaskEMCALClusterize.C"); // $ALICE_ROOT/PWG4/CaloCalib/macros
+ gROOT->LoadMacro("AddTaskCaloTrackCorr.C"); // $ALICE_ROOT/PWGGA/CaloTrackCorrelations/macros
+ gROOT->LoadMacro("AddTaskEMCALClusterize.C"); // $ALICE_ROOT/PWGGA/EMCALTasks/macros
// ------
if (mode == mPROOF) {
//TProof::Mgr("ccalpmaster")->SetROOTVersion("ALICE_v5-27-06b");
gROOT->LoadMacro("/afs/in2p3.fr/group/alice/laf/EnableAliRootForLAF.C");
- TProof* proof = EnableAliRootForLAF("ccaplmaster",nPROOFWorkers.Data(),ccin2p3UserName.Data(),alienUserName.Data(),"",kFALSE,kTRUE,kTRUE,"OADB:ANALYSIS:ANALYSISalice:AOD:ESD:CORRFW:STEERBase:EMCALUtils:PHOSUtils:PWG4PartCorrBase:PWG4PartCorrDep:PWG4CaloCalib");
+ TProof* proof = EnableAliRootForLAF("ccaplmaster",nPROOFWorkers.Data(),ccin2p3UserName.Data(),alienUserName.Data(),"",kFALSE,kTRUE,kTRUE,"OADB:ANALYSIS:ANALYSISalice:AOD:ESD:CORRFW:STEERBase:EMCALUtils:PHOSUtils:PWGCaloTrackCorrBase:PWGGACaloTrackCorrelations:PWGGAEMCALTasks");
// TProof* proof = TProof::Open("ccaplmaster",Form("workers=%s",nPROOFWorkers.Data()));
// //proof->UploadPackage("JETAN");
// proof->UploadPackage("PHOSUtils");
// proof->UploadPackage("EMCALUtils");
- // proof->UploadPackage("PWG4PartCorrBase");
- // proof->UploadPackage("PWG4PartCorrDep");
- // proof->UploadPackage("PWG4CaloCalib");
+ // proof->UploadPackage("PWGCaloTrackBase");
+ // proof->UploadPackage("PWGGACaloTrackCorrelations");
+ // proof->UploadPackage("PWGGAEMCALTasks");
// proof->EnablePackage("STEERBase");
// proof->EnablePackage("ESD");
// //proof->EnablePackage("JETAN");
// proof->EnablePackage("PHOSUtils");
// proof->EnablePackage("EMCALUtils");
- // proof->EnablePackage("PWG4PartCorrBase");
- // proof->EnablePackage("PWG4PartCorrDep");
- // proof->EnablePackage("PWG4CaloCalib");
+ // proof->EnablePackage("PWGCaloTrackBase");
+ // proof->EnablePackage("PWGGACaloTrackCorrelations");
+ // proof->EnablePackage("PWGGAEMCALTasks");
return;
}
gSystem->Load("libPHOSUtils");
gSystem->Load("libEMCALUtils");
- gSystem->Load("libPWG4PartCorrBase");
- gSystem->Load("libPWG4PartCorrDep");
- gSystem->Load("libPWG4CaloCalib");
- //SetupPar("PWG4PartCorrBase");
- //SetupPar("PWG4PartCorrDep");
- //SetupPar("PWG4CaloCalib");
+ gSystem->Load("libPWGCaloTrackBase");
+ gSystem->Load("libPWGGACaloCorrelations");
+ gSystem->Load("libPWGGAEMCALTasks");
+ //SetupPar("PWGCaloTrackBase");
+ //SetupPar("PWGGACaloCorrelations");
+ //SetupPar("PWGGAEMCALTasks");
//gSystem->Load("libJETAN");
//gSystem->Load("FASTJETAN");
- //gSystem->Load("PWG4JetTasks");
+ //gSystem->Load("PWGJE");
gSystem->Load("libCORRFW.so");
- gSystem->Load("libPWG4GammaConv.so");
- //SetupPar("PWG4GammaConv");
+ gSystem->Load("libPWGGAGammaConv.so");
+ //SetupPar("PWGGAGammaConv");
// needed for plugin?
gSystem->AddIncludePath("-I$ALICE_ROOT");
//correlation analysis
gROOT->LoadMacro("AddTaskCaloTrackCorrM.C");
- AliAnalysisTaskParticleCorrelationM *taskEMCAL = AddTaskCaloTrackCorrM(kInputData,"EMCAL",kFALSE);
+ AliAnalysisTaskCaloTrackCorrelationM *taskEMCAL = AddTaskCaloTrackCorrM(kInputData,"EMCAL",kFALSE);
mgr->AddTask(taskEMCAL);
- AliAnalysisTaskParticleCorrelationM *taskPHOS = AddTaskCaloTrackCorrM(kInputData,"PHOS", kFALSE);
+ AliAnalysisTaskCaloTrackCorrelationM *taskPHOS = AddTaskCaloTrackCorrM(kInputData,"PHOS", kFALSE);
mgr->AddTask(taskPHOS);
//gROOT->LoadMacro("AddTaskChargeCorr.C");
- AliAnalysisTaskParticleCorrelationM *taskCharge = AddTaskCaloTrackCorrM(kInputData, "CTS",kFALSE);
+ AliAnalysisTaskCaloTrackCorrelationM *taskCharge = AddTaskCaloTrackCorrM(kInputData, "CTS",kFALSE);
// if(!kMC)
// taskCharge->SelectCollisionCandidates();
mgr->AddTask(taskCharge);
SetupPar("PHOSUtils");
SetupPar("EMCALUtils");
- SetupPar("PWG4CaloTrackCorrBase");
- SetupPar("PWG4CaloTrackCorrDep");
+ SetupPar("PWGCaloTrackCorrBase");
+ SetupPar("PWGGACaloTrackCorrelations");
}
else{
//--------------------------------------------------------
gSystem->Load("libANALYSISalice");
gSystem->Load("libPHOSUtils");
gSystem->Load("libEMCALUtils");
- gSystem->Load("libPWG4CaloTrackCorrBase");
- gSystem->Load("libPWG4CaloTrackCorrDep");
+ gSystem->Load("libPWGCaloTrackCorrBase");
+ gSystem->Load("libPWGGACaloTrackCorrelations");
}
}
TGeoHMatrix* matrix[10];
AliEMCALRecoUtils * reco = clusterize->GetRecoUtils();
- gROOT->LoadMacro("ConfigureEMCALRecoUtils.C"); // $ALICE_ROOT/PWG4/CaloCalib/macros
+ gROOT->LoadMacro("ConfigureEMCALRecoUtils.C"); // $ALICE_ROOT/PWGGA/EMCALTasks/macros
ConfigureEMCALRecoUtils(
reco,
bMC,
/* $Id: $ */
//--------------------------------------------------
-// Example macro to do analysis with the
-// analysis classes in PWG4PartCorr
+// Example macro to do EMCAL calibration analysis
// Can be executed with Root and AliRoot
//
// Pay attention to the options and definitions
// change whatever you need for your analysis case
// ------------------------------------------------------------------
LoadLibraries(mode) ;
- //gSystem->Unload("libPWG4CaloCalib.so");
+ //gSystem->Unload("libPWGGAEMCALTasks.so");
//Try to set the new library
- //gSystem->Load("./PWG4CaloCalib/libPWG4CaloCalib.so");
+ //gSystem->Load("./PWGGAEMCALTasks/libPWGGAEMCALTasks.so");
//gSystem->ListLibraries();
-
+
//-------------------------------------------------------------------------------------------------
//Create chain from ESD and from cross sections files, look below for options.
//-------------------------------------------------------------------------------------------------
cout<<"Wrong data type "<<kInputData<<endl;
break;
}
-
+
TChain *chain = new TChain(kTreeName) ;
CreateChain(mode, chain);
-
+
if(chain){
AliLog::SetGlobalLogLevel(AliLog::kError);//Minimum prints on screen
// Make the analysis manager
//-------------------------------------
AliAnalysisManager *mgr = new AliAnalysisManager("Manager", "Manager");
-
+
//input
if(kInputData == "ESD")
- {
- // ESD handler
- AliESDInputHandler *esdHandler = new AliESDInputHandler();
- mgr->SetInputEventHandler(esdHandler);
- esdHandler->SetReadFriends(kFALSE);
- cout<<"ESD handler "<<mgr->GetInputEventHandler()<<endl;
- }
+ {
+ // ESD handler
+ AliESDInputHandler *esdHandler = new AliESDInputHandler();
+ mgr->SetInputEventHandler(esdHandler);
+ esdHandler->SetReadFriends(kFALSE);
+ cout<<"ESD handler "<<mgr->GetInputEventHandler()<<endl;
+ }
if(kInputData == "AOD")
- {
- // AOD handler
- AliAODInputHandler *aodHandler = new AliAODInputHandler();
- mgr->SetInputEventHandler(aodHandler);
- cout<<"AOD handler "<<mgr->GetInputEventHandler()<<endl;
-
- }
+ {
+ // AOD handler
+ AliAODInputHandler *aodHandler = new AliAODInputHandler();
+ mgr->SetInputEventHandler(aodHandler);
+ cout<<"AOD handler "<<mgr->GetInputEventHandler()<<endl;
+
+ }
// mgr->SetDebugLevel(1);
//Define task, put here any other task that you want to use.
//-------------------------------------------------------------------------
- // ESD filter task
+/*
if(kInputData == "ESD"){
-
- //gROOT->LoadMacro("AddTaskPhysicsSelection.C");
- //AliPhysicsSelectionTask* physSelTask = AddTaskPhysicsSelection();
-
+
+ gROOT->LoadMacro("AddTaskPhysicsSelection.C");
+ AliPhysicsSelectionTask* physSelTask = AddTaskPhysicsSelection();
+
}
+*/
// Create containers for input/output
AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer();
-
+
AliAnalysisTaskEMCALPi0CalibSelection * pi0calib = new AliAnalysisTaskEMCALPi0CalibSelection ("EMCALPi0Calibration");
- if(kInputData == "ESD") pi0calib->SelectCollisionCandidates();
+ if(kInputData == "ESD") pi0calib->SelectCollisionCandidates();
//pi0calib->SetDebugLevel(10);
//pi0calib->UseFilteredEventAsInput();
pi0calib->SetClusterMinEnergy(0.5);
pi0calib->SetNCellsGroup(0);
pi0calib->SwitchOnSameSM();
//pi0calib->SwitchOnOldAODs();
-
+
+/*
TGeoHMatrix *matrix[4];
double rotationMatrix[4][9] = {-0.014585, -0.999892, -0.002031, 0.999892, -0.014589, 0.001950, -0.001979, -0.002003, 0.999996,
- -0.014585, 0.999892, 0.002031, 0.999892, 0.014589, -0.001950, -0.001979, 0.002003, -0.999996,
- -0.345861, -0.938280, -0.003412, 0.938281, -0.345869, 0.001950, -0.003010, -0.002527, 0.999992,
- -0.345861, 0.938280, 0.003412, 0.938281, 0.345869, -0.001950, -0.003010, 0.002527, -0.999992};
+ -0.014585, 0.999892, 0.002031, 0.999892, 0.014589, -0.001950, -0.001979, 0.002003, -0.999996,
+ -0.345861, -0.938280, -0.003412, 0.938281, -0.345869, 0.001950, -0.003010, -0.002527, 0.999992,
+ -0.345861, 0.938280, 0.003412, 0.938281, 0.345869, -0.001950, -0.003010, 0.002527, -0.999992};
double translationMatrix[4][3] = {0.367264, 446.508738, 175.97185+0.3,
- 1.078181, 445.826258, -174.026758+0.3,
- -153.843916, 418.304256, 175.956905+0.8,
- -152.649580, 417.621779, -174.040392+0.8};
+ 1.078181, 445.826258, -174.026758+0.3,
+ -153.843916, 418.304256, 175.956905+0.8,
+ -152.649580, 417.621779, -174.040392+0.8};
for(int j=0; j<4; j++)
- {
- matrix[j] = new TGeoHMatrix();
- matrix[j]->SetRotation(rotationMatrix[j]);
- matrix[j]->SetTranslation(translationMatrix[j]);
- matrix[j]->Print();
- pi0calib->SetGeometryMatrixInSM(matrix[j],j);
- }
+ {
+ matrix[j] = new TGeoHMatrix();
+ matrix[j]->SetRotation(rotationMatrix[j]);
+ matrix[j]->SetTranslation(translationMatrix[j]);
+ matrix[j]->Print();
+ pi0calib->SetGeometryMatrixInSM(matrix[j],j);
+ }
+
-
pi0calib->SwitchOnLoadOwnGeometryMatrices();
-
+*/
+
pi0calib->SwitchOnClusterCorrection();
AliEMCALRecoUtils * reco = pi0calib->GetEMCALRecoUtils();
- reco->SetParticleType(AliEMCALRecoUtils::kPhoton);
+ reco->SetParticleType(AliEMCALRecoUtils::kPhoton);
reco->SetW0(4.5);
+
+ //---------------------
+ // Geometry alignment
+ //---------------------
//reco->SetPositionAlgorithm(AliEMCALRecoUtils::kUnchanged);
-
reco->SetPositionAlgorithm(AliEMCALRecoUtils::kPosTowerGlobal);
+
//reco->SetMisalTransShift(0,1.134); reco->SetMisalTransShift(1,8.2); reco->SetMisalTransShift(2,1.197);
//reco->SetMisalTransShift(3,-3.093); reco->SetMisalTransShift(4,6.82);reco->SetMisalTransShift(5,1.635);
-
+
//reco->SetPositionAlgorithm(AliEMCALRecoUtils::kPosTowerIndex);
//reco->SetMisalTransShift(0,1.08); reco->SetMisalTransShift(1,8.35); reco->SetMisalTransShift(2,1.12);
//reco->SetMisalRotShift(3,-8.05); reco->SetMisalRotShift(4,8.05);
//reco->SetMisalTransShift(3,-0.42); reco->SetMisalTransShift(5,1.55);
-
+
+ //---------------------
+ // Non linearity
+ //---------------------
+
reco->SetNonLinearityFunction(AliEMCALRecoUtils::kNoCorrection);
//reco->SetNonLinearityFunction(AliEMCALRecoUtils::kPi0GammaGamma);
//reco->SetNonLinearityParam(0,1.04); reco->SetNonLinearityParam(1,-0.1445);
//reco->SetNonLinearityParam(2,1.046);
+
+ // reco->SetNonLinearityFunction(AliEMCALRecoUtils::kPi0GammaConversion);
+ // reco->SetNonLinearityParam(0,1.033); reco->SetNonLinearityParam(1,0.0566186);
+ // reco->SetNonLinearityParam(2,0.982133);
+
+
+ // reco->SetNonLinearityFunction(AliEMCALRecoUtils::kPi0MC);
+ // reco->SetNonLinearityParam(0,1.001); reco->SetNonLinearityParam(1,-0.01264);
+ // reco->SetNonLinearityParam(2,-0.03632);
+ // reco->SetNonLinearityParam(3,0.1798); reco->SetNonLinearityParam(4,-0.522);
+
+
+ //---------------------
+ // Recalibration
+ //---------------------
+
+ //
+ // reco->SwitchOnRecalibration();
+ // TFile * f = new TFile("RecalibrationFactors.root","read");
+ // TH2F * h0 = (TH2F*)f->Get("EMCALRecalFactors_SM0")->Clone();
+ // TH2F * h1 = (TH2F*)f->Get("EMCALRecalFactors_SM1")->Clone();
+ // TH2F * h2 = (TH2F*)f->Get("EMCALRecalFactors_SM2")->Clone();
+ // TH2F * h3 = (TH2F*)f->Get("EMCALRecalFactors_SM3")->Clone();
+
+ // reco->SetEMCALChannelRecalibrationFactors(0,h0);
+ // reco->SetEMCALChannelRecalibrationFactors(1,h1);
+ // reco->SetEMCALChannelRecalibrationFactors(2,h2);
+ // reco->SetEMCALChannelRecalibrationFactors(3,h3);
+
+ //---------------------
+ // Bad channels removal
+ //---------------------
-// reco->SetNonLinearityFunction(AliEMCALRecoUtils::kPi0GammaConversion);
-// reco->SetNonLinearityParam(0,1.033); reco->SetNonLinearityParam(1,0.0566186);
-// reco->SetNonLinearityParam(2,0.982133);
-
-
-// reco->SetNonLinearityFunction(AliEMCALRecoUtils::kPi0MC);
-// reco->SetNonLinearityParam(0,1.001); reco->SetNonLinearityParam(1,-0.01264);
-// reco->SetNonLinearityParam(2,-0.03632);
-// reco->SetNonLinearityParam(3,0.1798); reco->SetNonLinearityParam(4,-0.522);
-
-// reco->SwitchOnRecalibration();
-// TFile * f = new TFile("RecalibrationFactors.root","read");
-// TH2F * h0 = (TH2F*)f->Get("EMCALRecalFactors_SM0")->Clone();
-// TH2F * h1 = (TH2F*)f->Get("EMCALRecalFactors_SM1")->Clone();
-// TH2F * h2 = (TH2F*)f->Get("EMCALRecalFactors_SM2")->Clone();
-// TH2F * h3 = (TH2F*)f->Get("EMCALRecalFactors_SM3")->Clone();
-
-// reco->SetEMCALChannelRecalibrationFactors(0,h0);
-// reco->SetEMCALChannelRecalibrationFactors(1,h1);
-// reco->SetEMCALChannelRecalibrationFactors(2,h2);
-// reco->SetEMCALChannelRecalibrationFactors(3,h3);
-
+/*
reco->SwitchOnBadChannelsRemoval();
+
// SM0
reco->SetEMCALChannelStatus(0,3,13); reco->SetEMCALChannelStatus(0,44,1); reco->SetEMCALChannelStatus(0,3,13);
reco->SetEMCALChannelStatus(0,20,7); reco->SetEMCALChannelStatus(0,38,2);
reco->SetEMCALChannelStatus(2,19,22);
//SM3
reco->SetEMCALChannelStatus(3,4,7);
-
+*/
+
reco->SetNumberOfCellsFromEMCALBorder(1);
-
+
//reco->Print("");
-
+
pi0calib->PrintInfo();
mgr->AddTask(pi0calib);
-
+
AliAnalysisDataContainer *coutput2 =
mgr->CreateContainer("pi0calib", TList::Class(), AliAnalysisManager::kOutputContainer, "pi0calib.root");
AliAnalysisDataContainer *cout_cuts = mgr->CreateContainer("Cuts", TList::Class(),
- AliAnalysisManager::kOutputContainer, "pi0calib.root");
+ AliAnalysisManager::kOutputContainer, "pi0calib.root");
mgr->ConnectInput (pi0calib, 0, cinput1);
mgr->ConnectOutput (pi0calib, 1, coutput2 );
mgr->ConnectOutput (pi0calib, 2, cout_cuts);
-
+
//-----------------------
// Run the analysis
//-----------------------
mgr->InitAnalysis();
mgr->PrintStatus();
mgr->StartAnalysis(smode.Data(),chain);
-
-cout <<" Analysis ended sucessfully "<< endl ;
-
+
+ cout <<" Analysis ended sucessfully "<< endl ;
+
}
else cout << "Chain was not produced ! "<<endl;
// If you want to use already compiled libraries
// in the aliroot distribution
//--------------------------------------------------------
-
+
gSystem->Load("libSTEERBase.so");
gSystem->Load("libESD.so");
gSystem->Load("libAOD.so");
gSystem->Load("libANALYSISalice.so");
gSystem->Load("libANALYSISalice.so");
//TGeoManager::Import("geometry.root") ; //need file "geometry.root" in local dir!!!!
- gSystem->Load("libPWG4CaloCalib.so");
- //SetupPar("PWG4CaloCalib");
- /*
- // gSystem->Load("libPWG4omega3pi.so");
- // gSystem->Load("libCORRFW.so");
- // gSystem->Load("libPWGHFbase.so");
- // gSystem->Load("libPWGmuon.so");
- */
+ gSystem->Load("libPWGGAEMCALTasks.so");
+ //SetupPar("PWGGAEMCALTasks");
+
//--------------------------------------------------------
//If you want to use root and par files from aliroot
//--------------------------------------------------------
SetupPar("EMCALUtils");
//Create Geometry
TGeoManager::Import("geometry.root") ; //need file "geometry.root" in local dir!!!!
- SetupPar("PWG4CaloCalib");
-*/
+ SetupPar("PWGGAEMCALTasks");
+ */
}
-
+
//---------------------------------------------------------
// <<<<<<<<<< PROOF mode >>>>>>>>>>>>
//---------------------------------------------------------
// gProof->ClearPackage("ESD");
// gProof->ClearPackage("AOD");
// gProof->ClearPackage("ANALYSIS");
- // gProof->ClearPackage("PWG4PartCorrBase");
- // gProof->ClearPackage("PWG4PartCorrDep");
// Enable the STEERBase Package
gProof->UploadPackage("STEERBase.par");
// Enable the PHOS geometry Package
//gProof->UploadPackage("PHOSUtils.par");
//gProof->EnablePackage("PHOSUtils");
- // Enable PartCorr analysis
- gProof->UploadPackage("PWG4PartCorrBase.par");
- gProof->EnablePackage("PWG4PartCorrBase");
- gProof->UploadPackage("PWG4PartCorrDep.par");
- gProof->EnablePackage("PWG4PartCorrDep");
gProof->ShowEnabledPackages();
}
//Load par files, create analysis libraries
//For testing, if par file already decompressed and modified
//classes then do not decompress.
-
+
TString cdir(Form("%s", gSystem->WorkingDirectory() )) ;
TString parpar(Form("%s.par", pararchivename)) ;
if ( gSystem->AccessPathName(parpar.Data()) ) {
//If you want to add several ESD files sitting in a common directory INDIR
//Specify as environmental variables the directory (INDIR), the number of files
//to analyze (NFILES) and the pattern name of the directories with files (PATTERN)
-
+
if(gSystem->Getenv("INDIR"))
kInDir = gSystem->Getenv("INDIR") ;
else cout<<"INDIR not set, use default: "<<kInDir<<endl;
if ( kInDir && kFile) {
printf("Get %d files from directory %s\n",kFile,kInDir);
if ( ! gSystem->cd(kInDir) ) {//check if ESDs directory exist
- printf("%s does not exist\n", kInDir) ;
- return ;
+ printf("%s does not exist\n", kInDir) ;
+ return ;
}
-
-
+
+
cout<<"INDIR : "<<kInDir<<endl;
cout<<"NFILES : "<<kFile<<endl;
cout<<"PATTERN : " <<kPattern<<endl;
-
+
TString datafile="";
if(kInputData == "ESD") datafile = "AliESDs.root" ;
else if(kInputData == "AOD") datafile = "aod.root" ;
char file[120] ;
for (event = 0 ; event < kFile ; event++) {
- sprintf(file, "%s/%s%d/%s", kInDir,kPattern,event,datafile.Data()) ;
- TFile * fESD = 0 ;
- //Check if file exists and add it, if not skip it
- if ( fESD = TFile::Open(file)) {
- if ( fESD->Get(kTreeName) ) {
- printf("++++ Adding %s\n", file) ;
- chain->AddFile(file);
- }
- }
- else {
- printf("---- Skipping %s\n", file) ;
- skipped++ ;
- }
+ sprintf(file, "%s/%s%d/%s", kInDir,kPattern,event,datafile.Data()) ;
+ TFile * fESD = 0 ;
+ //Check if file exists and add it, if not skip it
+ if ( fESD = TFile::Open(file)) {
+ if ( fESD->Get(kTreeName) ) {
+ printf("++++ Adding %s\n", file) ;
+ chain->AddFile(file);
+ }
+ }
+ else {
+ printf("---- Skipping %s\n", file) ;
+ skipped++ ;
+ }
}
printf("number of entries # %lld, skipped %d\n", chain->GetEntries(), skipped*100) ;
}
else if(mode == mGRID){
//Get colection file. It is specified by the environmental
//variable XML
-
+
if(gSystem->Getenv("XML") )
kXML = gSystem->Getenv("XML");
else
return ;
}
else cout<<"XML file "<<kXML<<endl;
-
+
//Load necessary libraries and connect to the GRID
gSystem->Load("libNetx.so") ;
gSystem->Load("libRAliEn.so");
TGrid::Connect("alien://") ;
-
+
//Feed Grid with collection file
//TGridCollection * collection = (TGridCollection*)gROOT->ProcessLine(Form("TAlienCollection::Open(\"%s\", 0)", kXML));
TGridCollection * collection = (TGridCollection*) TAlienCollection::Open(kXML);
return kFALSE ;
}
TGridResult* result = collection->GetGridResult("",0 ,0);
-
+
// Makes the ESD chain
printf("*** Getting the Chain ***\n");
for (Int_t index = 0; index < result->GetEntries(); index++) {
gSystem->Load("libEMCALUtils.so");
//SetupPar("EMCALUtils");
- //SetupPar("PWG4CaloCalib");
+ //SetupPar("PWGGAEMCALTasks");
//TGeoManager::Import("geometry.root") ; //need file "geometry.root" in local dir!!!!
- gSystem->Load("libPWG4CaloCalib.so");
+ gSystem->Load("libPWGGAEMCALTasks.so");
/*
- // gSystem->Load("libPWG4omega3pi.so");
// gSystem->Load("libCORRFW.so");
// gSystem->Load("libPWGHFbase.so");
// gSystem->Load("libPWGmuon.so");
SetupPar("EMCALUtils");
//Create Geometry
TGeoManager::Import("geometry.root") ; //need file "geometry.root" in local dir!!!!
- SetupPar("PWG4CaloCalib");
+ SetupPar("PWGGAEMCALTasks");
*/
}
// gProof->ClearPackage("ESD");
// gProof->ClearPackage("AOD");
// gProof->ClearPackage("ANALYSIS");
- // gProof->ClearPackage("PWG4PartCorrBase");
- // gProof->ClearPackage("PWG4PartCorrDep");
// Enable the STEERBase Package
gProof->UploadPackage("STEERBase.par");
// Enable the PHOS geometry Package
//gProof->UploadPackage("PHOSUtils.par");
//gProof->EnablePackage("PHOSUtils");
- // Enable PartCorr analysis
- gProof->UploadPackage("PWG4PartCorrBase.par");
- gProof->EnablePackage("PWG4PartCorrBase");
- gProof->UploadPackage("PWG4PartCorrDep.par");
- gProof->EnablePackage("PWG4PartCorrDep");
gProof->ShowEnabledPackages();
}
mgr->ConnectInput(task, 0, cinput);
TString outputfile = AliAnalysisManager::GetCommonFileName();
- AliAnalysisDataContainer *coutput = mgr->CreateContainer("CaloConv", TList::Class(), AliAnalysisManager::kOutputContainer, Form("%s:PWG4_CaloConv",outputfile.Data()));
+ AliAnalysisDataContainer *coutput = mgr->CreateContainer("CaloConv", TList::Class(), AliAnalysisManager::kOutputContainer, Form("%s:PWGGA_CaloConv",outputfile.Data()));
mgr->ConnectOutput(task, 1, coutput);
- AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("CFCaloConv", TList::Class(), AliAnalysisManager::kOutputContainer, Form("%s:PWG4_CFCaloConv",outputfile.Data()));
+ AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("CFCaloConv", TList::Class(), AliAnalysisManager::kOutputContainer, Form("%s:PWGGA_CFCaloConv",outputfile.Data()));
mgr->ConnectOutput(task, 2, coutput2);
return task ;
// Connect I/O to the task
TString outputfile = AliAnalysisManager::GetCommonFileName();
- outputfile += ":PWG4_GammaConversion";
+ outputfile += ":PWGGA_GammaConversion";
outputfile += Form("_%s",kGCAnalysisCutSelectionId.Data());
AliAnalysisDataContainer *coutput3 = mgr->CreateContainer("histogramsAliGammaConversionDalitz", TList::Class(),AliAnalysisManager::kOutputContainer, outputfile);
mgr->ConnectInput ( dalitzTask, 0, mgr->GetCommonInputContainer() );
TString kGCoutputFileAppendix = "";
TString kGCdataList = "";
Bool_t kGCwriteNtuple = kFALSE;
-// WE DOO NOT NEED TO CHANGE THIS (kGCusePWG4PartCorr) ANYMORE SINCE IT IS TAKEN CARE OF AUTOMATICALLY NOW
-Bool_t kGCusePWG4PartCorr = kTRUE;
+// WE DOO NOT NEED TO CHANGE THIS (kGCusePWGGACaloTrackCorrelations) ANYMORE SINCE IT IS TAKEN CARE OF AUTOMATICALLY NOW
+Bool_t kGCusePWGGACaloTrackCorrelations = kTRUE;
/** Flag to enable running on train */
Bool_t kGCrunOnTrain = kFALSE;
}
void SetVersionLibrary(){
- // Check if the file $ALICE_ROOT/PWG4/GammaConv/AliAnalysisTaskGammaConversion.cxx exists.
- // If yes, we set kGCusePWG4PartCorr to false since we have a newer version
- // If no, kGCusePWG4PartCorr is true.
+ // Check if the file $ALICE_ROOT/PWGGA/GammaConv/AliAnalysisTaskGammaConversion.cxx exists.
+ // If yes, we set kGCusePWGGACaloTrackCorrelations to false since we have a newer version
+ // If no, kGCusePWGGACaloTrackCorrelations is true.
TString file = gSystem->Getenv("ALICE_ROOT");
- file+="/PWG4/PartCorr/AliAnalysisTaskGammaConversion.cxx";
+ file+="/PWGGA/PartCorr/AliAnalysisTaskGammaConversion.cxx";
ifstream stream;
stream.open(file.Data());
if(!stream){
- kGCusePWG4PartCorr=kFALSE;
+ kGCusePWGGACaloTrackCorrelations=kFALSE;
}
else{
- kGCusePWG4PartCorr=kTRUE;
+ kGCusePWGGACaloTrackCorrelations=kTRUE;
}
stream.close();
}
break;
}
- SetVersionLibrary(); // checks if PWG4GammaConv or PWG4PartCorr is used
+ SetVersionLibrary(); // checks if PWGGAGammaConv or PWGGACaloTrackCorrelations is used
if(cin_esd == NULL && kGCrunOnTrain == kTRUE){
cout<<"Error: kGCrunOnTrain flag is set to true but the input AliAnalysisDataContainer is NULL"<<endl;
// Declare Common Input Tchain
AliAnalysisDataContainer *cinput1 = NULL;
- if(kGCusePWG4PartCorr){
+ if(kGCusePWGGACaloTrackCorrelations){
if(kGCrunOnTrain == kFALSE){
cinput1 = mgr->CreateContainer("GammaConvChain",TChain::Class(),AliAnalysisManager::kInputContainer);
}
// Common Output Tree in common ??????default?????? output file
- // CKB kGCusePWG4PartCorr and writestandard are not mutually exclusive?
+ // CKB kGCusePWGGACaloTrackCorrelations and writestandard are not mutually exclusive?
AliAnalysisDataContainer *coutput1 = NULL;
if(kGCWriteAOD) {
TString outputfile = AliAnalysisManager::GetCommonFileName();
cout<<"Analyis cut selection ID is: "<<kGCAnalysisCutSelectionId.Data()<<endl;
- // outputfile += Form(":PWG4_GammaConversion_%s",kGCAnalysisCutSelectionId.Data());
- outputfile += Form(":PWG4_GammaConversion");
+ // outputfile += Form(":PWGGA_GammaConversion_%s",kGCAnalysisCutSelectionId.Data());
+ outputfile += Form(":PWGGA_GammaConversion");
// if(kGCrunNeutralMeson==kTRUE) outputfile +="1"; else outputfile +="0";
if( kGCrunDalitz ){
- gROOT->LoadMacro("$ALICE_ROOT/PWG4/macros/AddTaskGammaConvDalitz.C");
+ gROOT->LoadMacro("$ALICE_ROOT/PWGGA/macros/AddTaskGammaConvDalitz.C");
AddTaskGammaConvDalitz( v0Reader, kGCcalculateBackground, kGCRunStandalone );
}
gSystem->Load("libANALYSIS.so");
gSystem->Load("libANALYSISalice.so");
gSystem->Load("libCORRFW.so");
- gSystem->Load("libPWG4GammaConv.so");
+ gSystem->Load("libPWGGAGammaConv.so");
// gSystem->ChangeDirectory(pwd.Data());
gSystem->Load("libCORRFW.so");
////
- //Setting up PWG4GammaConv.par//
+ //Setting up PWGGAGammaConv.par//
////
- cout<<"compiling PWG4GammaConv"<<endl;
- setupPar("PWG4GammaConv");
- gSystem->Load("libPWG4GammaConv.so");
+ cout<<"compiling PWGGAGammaConv"<<endl;
+ setupPar("PWGGAGammaConv");
+ gSystem->Load("libPWGGAGammaConv.so");
gSystem->ChangeDirectory(pwd.Data());
}
/****************************************************************************************************************************
-****** provided by Gamma Conversion Group, PWG4, *****
+****** provided by Gamma Conversion Group, PWGGA, *****
****** Friederike Bock, friederike.bock@cern.ch *****
*****************************************************************************************************************************
*** This macro can be used to observe the effects of the different cuts used in the GammaConversion analysis, it can *****
char *StandardYAxis = "#gamma/ event scaled by multiplicity";
if(cutsel != ""){
- char *GammaDirectory = Form("PWG4_GammaConversion_%s", cutsel);
+ char *GammaDirectory = Form("PWGGA_GammaConversion_%s", cutsel);
cout << GammaDirectory << endl;
char *GammaList = Form("histogramsAliGammaConversion_%s", cutsel);
cout << GammaList << endl;
}else{
- char *GammaDirectory = "PWG4_GammaConversion";
+ char *GammaDirectory = "PWGGA_GammaConversion";
cout << GammaDirectory << endl;
char *GammaList = "histogramsAliGammaConversion";
cout << GammaList << endl;
Char_t filename_data[200] = (Form("%s%s",path,data));
TFile f(filename_data);
- TDirectory *fPWG4GammaConversion_data = new TDirectory(); // definition of first folder / list
+ TDirectory *fPWGGAGammaConversion_data = new TDirectory(); // definition of first folder / list
TList *fHistosGammaConversion_data = new TList(); // definition of first folder / list
TList *fESDContainer_data = new TList(); // definition of following folder / list
TList *fMappingContainer_data = new TList();
- if(!(fPWG4GammaConversion_data = (TDirectory*)f.Get(GammaDirectory))) cout <<"PWG4GammConversion TList NOT loaded correctly"<<endl;
- if(!(fHistosGammaConversion_data = (TList*)fPWG4GammaConversion_data->Get(GammaList))) cout<<"histogramsAliGammaConversion NOT loaded correctly!"<<endl;
+ if(!(fPWGGAGammaConversion_data = (TDirectory*)f.Get(GammaDirectory))) cout <<"PWGGAGammConversion TList NOT loaded correctly"<<endl;
+ if(!(fHistosGammaConversion_data = (TList*)fPWGGAGammaConversion_data->Get(GammaList))) cout<<"histogramsAliGammaConversion NOT loaded correctly!"<<endl;
if(!(fESDContainer_data = (TList*)fHistosGammaConversion_data->FindObject("ESD histograms"))) cout<<"ESD histograms NOT loaded correctly!"<<endl;
fMappingContainer_data = (TList*)fHistosGammaConversion_data->FindObject("Mapping histograms");
if(MC1){
// for new versions
- TDirectory *fPWG4GammaConversion_montecarlo = new TDirectory(); // definition of first folder / list
+ TDirectory *fPWGGAGammaConversion_montecarlo = new TDirectory(); // definition of first folder / list
TList *fHistosGammaConversion_montecarlo = new TList(); // definition of first folder / list
TList *fESDContainer_montecarlo = new TList(); // definition of following folder / list
TList *fMappingContainer_montecarlo = new TList(); // definition of following folder / list
- if(!(fPWG4GammaConversion_montecarlo = (TDirectory*)montecarlo->Get(GammaDirectory))) cout <<"PWG4GammConversion TList NOT loaded correctly"<<endl;
- if(!(fHistosGammaConversion_montecarlo = (TList*)fPWG4GammaConversion_montecarlo->Get(GammaList))) cout<<"histogramsAliGammaConversion NOT loaded correctly!"<<endl;
+ if(!(fPWGGAGammaConversion_montecarlo = (TDirectory*)montecarlo->Get(GammaDirectory))) cout <<"PWGGAGammConversion TList NOT loaded correctly"<<endl;
+ if(!(fHistosGammaConversion_montecarlo = (TList*)fPWGGAGammaConversion_montecarlo->Get(GammaList))) cout<<"histogramsAliGammaConversion NOT loaded correctly!"<<endl;
if(!(fESDContainer_montecarlo = (TList*)fHistosGammaConversion_montecarlo->FindObject("ESD histograms"))) cout<<"ESD histograms NOT loaded correctly!"<<endl;
fMappingContainer_montecarlo = (TList*)fHistosGammaConversion_montecarlo->FindObject("Mapping histograms");
-// provided by Gamma Conversion Group, $ALICE_ROOT/PWG4/GammaConv ;https://twiki.cern.ch/twiki/bin/view/ALICE/PWG4GammaConversion
+// provided by Gamma Conversion Group, $ALICE_ROOT/PWGGA/GammaConv ;https://twiki.cern.ch/twiki/bin/view/ALICE/PWG4GammaConversion
#include <stdlib.h>
#include <iostream>
TFile f(file.Data());
fCutSelection = cutSelection;
- TDirectory *pwg4dir =(TDirectory*)f.Get(Form("PWG4_GammaConversion_%s",fCutSelection.Data()));
+ TDirectory *pwg4dir =(TDirectory*)f.Get(Form("PWGGA_GammaConversion_%s",fCutSelection.Data()));
if(pwg4dir == NULL){
- cout<<Form("ERROR: PWG4_GammaConversion_%s",fCutSelection.Data())<<" is not found in the file"<<endl;
+ cout<<Form("ERROR: PWGGA_GammaConversion_%s",fCutSelection.Data())<<" is not found in the file"<<endl;
return;
}
TList *fHistosGammaConversion = (TList*)pwg4dir->Get(Form("histogramsAliGammaConversion_%s",fCutSelection.Data()));
-// provided by Gamma Conversion Group, PWG4, Kathrin Koch, kkoch@physi.uni-heidelberg.de
+// provided by Gamma Conversion Group, PWGGA, Kathrin Koch, kkoch@physi.uni-heidelberg.de
#include <fstream>
#include <Riostream.h>
c_mass->SetFillColor(0);
- TDirectory *pwg4dir =(TDirectory*)f.Get(Form("PWG4_GammaConversion_%s",cutSelection));
+ TDirectory *pwg4dir =(TDirectory*)f.Get(Form("PWGGA_GammaConversion_%s",cutSelection));
TList *fHistosGammaConversion = (TList*)pwg4dir->Get(Form("histogramsAliGammaConversion_%s",cutSelection));
TList *fESDContainer = (TList*)fHistosGammaConversion->FindObject("ESD histograms");
TList *fMCContainer = (TList*)fHistosGammaConversion->FindObject("MC histograms");
if(makeMappingPlots==kTRUE){
outputFile = new TFile(Form("%sPi0Characteristics.root",outputDir),"UPDATE");
}
- TDirectory *pwg4dir =(TDirectory*)f.Get(Form("PWG4_GammaConversion_%s",cutSelection));
+ TDirectory *pwg4dir =(TDirectory*)f.Get(Form("PWGGA_GammaConversion_%s",cutSelection));
TList *fHistosGammaConversion = (TList*)pwg4dir->Get(Form("histogramsAliGammaConversion_%s",cutSelection));
TList *fESDContainer = (TList*)fHistosGammaConversion->FindObject("ESD histograms");
TList *fMCContainer = (TList*)fHistosGammaConversion->FindObject("MC histograms");
-// provided by Gamma Conversion Group, PWG4, Kathrin Koch, kkoch@physi.uni-heidelberg.de
+// provided by Gamma Conversion Group, PWGGA, Kathrin Koch, kkoch@physi.uni-heidelberg.de
#include <Riostream>
#include <fstream>
TObject * o = f.Get(directories->At(entFile)->GetName()); // get the object in the base directory
- if(TString(o->IsA()->GetName())=="TDirectoryFile"){ // means that this is a directory (PWG4......)
+ if(TString(o->IsA()->GetName())=="TDirectoryFile"){ // means that this is a directory (PWGGA......)
TDirectory *pwg4dir =(TDirectory*)o;
TString baseDirName = pwg4dir->GetName();
- TString reconstructionFlagString = ""; // this is for new scheme where also the flags are coded in numbers in the PWG4.... name
+ TString reconstructionFlagString = ""; // this is for new scheme where also the flags are coded in numbers in the PWGGA.... name
if(baseDirName.Length()>31){
reconstructionFlagString = baseDirName(baseDirName.Index("GammaConversion_")+16,8);
/****************************************************************************************************************************
-****** provided by Gamma Conversion Group, PWG4, *****
+****** provided by Gamma Conversion Group, PWGGA, *****
****** Kathrin Koch, kkoch@physi.uni-heidelberg.de *****
****** Friederike Bock, friederike.bock@cern.ch *****
*****************************************************************************************************************************
Scaling->SetLineWidth(2);
if(cutsel != ""){
- char *GammaDirectory = Form("PWG4_GammaConversion_%s", cutsel);
+ char *GammaDirectory = Form("PWGGA_GammaConversion_%s", cutsel);
cout << GammaDirectory << endl;
char *GammaList = Form("histogramsAliGammaConversion_%s", cutsel);
cout << GammaList << endl;
}else{
- char *GammaDirectory = "PWG4_GammaConversion";
+ char *GammaDirectory = "PWGGA_GammaConversion";
cout << GammaDirectory << endl;
char *GammaList = "histogramsAliGammaConversion";
cout << GammaList << endl;
TFile f(filename_data);
//for new versions
- TDirectory *fPWG4GammaConversion_data = new TDirectory(); // definition of first folder / list
+ TDirectory *fPWGGAGammaConversion_data = new TDirectory(); // definition of first folder / list
TList *fHistosGammaConversion_data = new TList(); // definition of first folder / list
TList *fESDContainer_data = new TList(); // definition of following folder / list
TList *fMappingContainer_data = new TList();
- if(!(fPWG4GammaConversion_data = (TDirectory*)f.Get(GammaDirectory))) cout <<"PWG4GammConversion TList NOT loaded correctly"<<endl;
- if(!(fHistosGammaConversion_data = (TList*)fPWG4GammaConversion_data->Get(GammaList))) cout<<"histogramsAliGammaConversion NOT loaded correctly!"<<endl;
+ if(!(fPWGGAGammaConversion_data = (TDirectory*)f.Get(GammaDirectory))) cout <<"PWGGAGammConversion TList NOT loaded correctly"<<endl;
+ if(!(fHistosGammaConversion_data = (TList*)fPWGGAGammaConversion_data->Get(GammaList))) cout<<"histogramsAliGammaConversion NOT loaded correctly!"<<endl;
if(!(fESDContainer_data = (TList*)fHistosGammaConversion_data->FindObject("ESD histograms"))) cout<<"ESD histograms NOT loaded correctly!"<<endl;
fMappingContainer_data = (TList*)fHistosGammaConversion_data->FindObject("Mapping histograms");
// for new versions
- TDirectory *fPWG4GammaConversion_montecarlo = new TDirectory(); // definition of first folder / list
+ TDirectory *fPWGGAGammaConversion_montecarlo = new TDirectory(); // definition of first folder / list
TList *fHistosGammaConversion_montecarlo = new TList(); // definition of first folder / list
TList *fESDContainer_montecarlo = new TList(); // definition of following folder / list
TList *fMappingContainer_montecarlo = new TList(); // definition of following folder / list
TList *fMCtruth_montecarlo = new TList();
- if(!(fPWG4GammaConversion_montecarlo = (TDirectory*)montecarlo->Get(GammaDirectory))) cout <<"PWG4GammConversion TList NOT loaded correctly"<<endl;
- if(!(fHistosGammaConversion_montecarlo = (TList*)fPWG4GammaConversion_montecarlo->Get(GammaList))) cout<<"histogramsAliGammaConversion NOT loaded correctly!"<<endl;
+ if(!(fPWGGAGammaConversion_montecarlo = (TDirectory*)montecarlo->Get(GammaDirectory))) cout <<"PWGGAGammConversion TList NOT loaded correctly"<<endl;
+ if(!(fHistosGammaConversion_montecarlo = (TList*)fPWGGAGammaConversion_montecarlo->Get(GammaList))) cout<<"histogramsAliGammaConversion NOT loaded correctly!"<<endl;
if(!(fESDContainer_montecarlo = (TList*)fHistosGammaConversion_montecarlo->FindObject("ESD histograms"))) cout<<"ESD histograms NOT loaded correctly!"<<endl;
fMappingContainer_montecarlo = (TList*)fHistosGammaConversion_montecarlo->FindObject("Mapping histograms");
fMCtruth_montecarlo = (TList*)fHistosGammaConversion_montecarlo->FindObject("MC histograms");
/***********************************************************************************************
-*** provided by Gamma Conversion Group, PWG4, ******
+*** provided by Gamma Conversion Group, PWGGA, ******
*** Friederike Bock, fbock@physi.uni-heidelberg.de *** ******
************************************************************************************************
************************************************************************************************
Char_t filename_MC[100] = (Form("%s%s",path,MCfile));
if(cutsel != ""){
- char *GammaDirectory = Form("PWG4_GammaConversion_%s", cutsel);
+ char *GammaDirectory = Form("PWGGA_GammaConversion_%s", cutsel);
cout << GammaDirectory << endl;
char *GammaList = Form("histogramsAliGammaConversion_%s", cutsel);
cout << GammaList << endl;
}else{
- char *GammaDirectory = "PWG4_GammaConversion";
+ char *GammaDirectory = "PWGGA_GammaConversion";
cout << GammaDirectory << endl;
char *GammaList = "histogramsAliGammaConversion";
cout << GammaList << endl;
TFile f(filename_MC);
//************************** Container Loading ********************************************************************
- TDirectory *fPWG4GammaConversion_montecarlo = new TDirectory();
+ TDirectory *fPWGGAGammaConversion_montecarlo = new TDirectory();
TList *fHistosGammaConversion_montecarlo = new TList();
TList *fResolutionContainer_montecarlo = new TList();
TList *fESDContainer_montecarlo = new TList();
- if(!(fPWG4GammaConversion_montecarlo = (TDirectory*)f.Get(GammaDirectory))) cout <<"PWG4GammConversion TList NOT loaded correctly"<<endl;
- if(!(fHistosGammaConversion_montecarlo = (TList*)fPWG4GammaConversion_montecarlo->Get(GammaList))) cout<<"histogramsAliGammaConversion NOT loaded correctly!"<<endl;
+ if(!(fPWGGAGammaConversion_montecarlo = (TDirectory*)f.Get(GammaDirectory))) cout <<"PWGGAGammConversion TList NOT loaded correctly"<<endl;
+ if(!(fHistosGammaConversion_montecarlo = (TList*)fPWGGAGammaConversion_montecarlo->Get(GammaList))) cout<<"histogramsAliGammaConversion NOT loaded correctly!"<<endl;
if(fResolutionContainer_montecarlo = (TList*)fHistosGammaConversion_montecarlo->FindObject("Resolution histograms")) cout<<"list loaded correctly!"<<endl;
if(!(fESDContainer_montecarlo = (TList*)fHistosGammaConversion_montecarlo->FindObject("ESD histograms"))) cout<<"ESD histograms NOT loaded correctly!"<<endl;
TH1F * ESD_ConvGamma_Pt=fESDContainer_montecarlo->FindObject("ESD_ConvGamma_Pt");
-// provided by Gamma Conversion Group, PWG4, Kathrin Koch, kkoch@physi.uni-heidelberg.de and Friederike Bock, fbock@physi.uni-heidelberg.de
+// provided by Gamma Conversion Group, PWGGA, Kathrin Koch, kkoch@physi.uni-heidelberg.de and Friederike Bock, fbock@physi.uni-heidelberg.de
#include <Riostream.h>
#include <fstream>
// choice of dateset
if(cutsel != ""){
- char *GammaDirectory = Form("PWG4_GammaConversion_%s", cutsel);
+ char *GammaDirectory = Form("PWGGA_GammaConversion_%s", cutsel);
cout << GammaDirectory << endl;
char *GammaList = Form("histogramsAliGammaConversion_%s", cutsel);
cout << GammaList << endl;
}else{
- char *GammaDirectory = "PWG4_GammaConversion";
+ char *GammaDirectory = "PWGGA_GammaConversion";
cout << GammaDirectory << endl;
char *GammaList = "histogramsAliGammaConversion";
cout << GammaList << endl;
EtaRange->SetLineWidth(6);
//------------------------------- Reading FILES ----------------------------------------------------------------------
- TDirectory *fPWG4GammaConversion_input1 = new TDirectory(); // definition of first folder / list
+ TDirectory *fPWGGAGammaConversion_input1 = new TDirectory(); // definition of first folder / list
TList *fHistosGammaConversion_input1 = new TList(); // definition of first folder / list
TList *fESDContainer_input1 = new TList(); // definition of following folder / list
TList *fMappingContainer_input1 = new TList();
- fPWG4GammaConversion_input1 = (TDirectory*)f1.Get(GammaDirectory);
- fHistosGammaConversion_input1 = (TList*)fPWG4GammaConversion_input1->Get(GammaList);
+ fPWGGAGammaConversion_input1 = (TDirectory*)f1.Get(GammaDirectory);
+ fHistosGammaConversion_input1 = (TList*)fPWGGAGammaConversion_input1->Get(GammaList);
fMappingContainer_input1 = (TList*)fHistosGammaConversion_input1->FindObject("Mapping histograms");
fESDContainer_input1 = (TList*)fHistosGammaConversion_input1->FindObject("ESD histograms");
input2 = new TFile(Form("%s%s",path, secondinput));
- TDirectory *fPWG4GammaConversion_input2 = new TDirectory(); // definition of first folder / list
+ TDirectory *fPWGGAGammaConversion_input2 = new TDirectory(); // definition of first folder / list
TList *fHistosGammaConversion_input2 = new TList(); // definition of first folder / list
TList *fMappingContainer_input2 = new TList(); // definition of following folder / list
TList *fESDContainer_input2 = new TList(); // definition of following folder / list
- fPWG4GammaConversion_input2 = (TDirectory*)input2->Get(GammaDirectory);
- fHistosGammaConversion_input2 = (TList*)fPWG4GammaConversion_input2->Get(GammaList);
+ fPWGGAGammaConversion_input2 = (TDirectory*)input2->Get(GammaDirectory);
+ fHistosGammaConversion_input2 = (TList*)fPWGGAGammaConversion_input2->Get(GammaList);
fMappingContainer_input2 = (TList*)fHistosGammaConversion_input2->FindObject("Mapping histograms");
fESDContainer_input2 = (TList*)fHistosGammaConversion_input2->FindObject("ESD histograms");
/***********************************************************************************************
-*** provided by Gamma Conversion Group, PWG4, ******
+*** provided by Gamma Conversion Group, PWGGA, ******
*** Friederike Bock, fbock@physi.uni-heidelberg.de *** ******
************************************************************************************************/
/***********************************************************************************************
-*** provided by Gamma Conversion Group, PWG4, ******
+*** provided by Gamma Conversion Group, PWGGA, ******
*** Friederike Bock, fbock@physi.uni-heidelberg.de *** ******
************************************************************************************************/
fi
-root -b -q $ALICE_ROOT/PWG4/GammaConv/macros/MakeCutLog.C\(\"$RootFile\"\,\"$InputDirectory\"\,\"$OutputDirectory\"\)
+root -b -q $ALICE_ROOT/PWGGA/GammaConv/macros/MakeCutLog.C\(\"$RootFile\"\,\"$InputDirectory\"\,\"$OutputDirectory\"\)
# Read the different cuts form the Cut selection log file
exec<"$OutputDirectory/CutSelection.log"
while read cutSelection
do
echo CutSelection is $cutSelection;
- root -b -q $ALICE_ROOT/PWG4/GammaConv/macros/Extract_IntegratedPi0Yield.C\(\"$cutSelection\"\,\"$RootFile\"\,\"$InputDirectory\"\,\"$OutputDirectory\"\"\);
+ root -b -q $ALICE_ROOT/PWGGA/GammaConv/macros/Extract_IntegratedPi0Yield.C\(\"$cutSelection\"\,\"$RootFile\"\,\"$InputDirectory\"\,\"$OutputDirectory\"\"\);
- root -b -q $ALICE_ROOT/PWG4/GammaConv/macros/Extract_Pi0_Characteristics.C\(\"$cutSelection\"\,\"$RootFile\"\,\"$InputDirectory\"\,\"$OutputDirectory\"\);
+ root -b -q $ALICE_ROOT/PWGGA/GammaConv/macros/Extract_Pi0_Characteristics.C\(\"$cutSelection\"\,\"$RootFile\"\,\"$InputDirectory\"\,\"$OutputDirectory\"\);
done
- root -b -q $ALICE_ROOT/PWG4/GammaConv/macros/Plot_IntegratedPi0Yield.C\(\"$DatFilenameBase\"\,\"$OutputDirectory\"\,\"$Suffix\"\)
+ root -b -q $ALICE_ROOT/PWGGA/GammaConv/macros/Plot_IntegratedPi0Yield.C\(\"$DatFilenameBase\"\,\"$OutputDirectory\"\,\"$Suffix\"\)
- root -b -q $ALICE_ROOT/PWG4/GammaConv/macros/Plot_Pi0_Characteristics.C
\ No newline at end of file
+ root -b -q $ALICE_ROOT/PWGGA/GammaConv/macros/Plot_Pi0_Characteristics.C
\ No newline at end of file
//
// Usage example for EMCAL:
//
- // gSystem->Load("libPWG4UserTasks.so");
- // gROOT->LoadMacro("$ALICE_ROOT/PWG4/UserTasks/CaloCellQA/macros/AddTaskCaloCellsQA.C");
+ // gSystem->Load("libPWGGAPHOSTasks.so");
+ // gROOT->LoadMacro("$ALICE_ROOT/PWGGA/PHOSTasks/CaloCellQA/macros/AddTaskCaloCellsQA.C");
// AliAnalysisTaskCaloCellsQA *taskQA = AddTaskCaloCellsQA(10); // 10 supermodules
// taskQA->SelectCollisionCandidates(AliVEvent::kMB); // if necessary
// // taskQA->SetAvoidPileup(kFALSE); // some customization
//
// Usage example for PHOS:
//
- // gSystem->Load("libPWG4UserTasks.so");
- // gROOT->LoadMacro("$ALICE_ROOT/PWG4/UserTasks/CaloCellQA/macros/AddTaskCaloCellsQA.C");
+ // gSystem->Load("libPWGAPHOSTasks.so");
+ // gROOT->LoadMacro("$ALICE_ROOT/PWGGA/PHOSTasks/CaloCellQA/macros/AddTaskCaloCellsQA.C");
// AliAnalysisTaskCaloCellsQA *taskQA = AddTaskCaloCellsQA(4, 1);
// taskQA->SelectCollisionCandidates(AliVEvent::kMB); // if necessary
// taskQA->GetCaloCellsQA()->SetClusterEnergyCuts(0.3,0.1); // increase statistics
AliAnalysisTaskSE * AddTaskPHOSQA()
{
- gROOT->LoadMacro("$ALICE_ROOT/PWG4/UserTasks/CaloCellQA/macros/AddTaskCaloCellsQA.C");
+ gROOT->LoadMacro("$ALICE_ROOT/PWGGA/PHOSTasks/CaloCellQA/macros/AddTaskCaloCellsQA.C");
AliAnalysisTaskCaloCellsQA *taskPHOSCellQA = AddTaskCaloCellsQA(4, 1, "CaloCellsQA.root","PHOSCellsQA");
taskPHOSCellQA->GetCaloCellsQA()->SetClusterEnergyCuts(0.3,0.3,1.0);
// Example macro to run QA
// load relevant library
- gSystem->Load("libPWG4UserTasks.so");
+ gSystem->Load("libPWGGAPHOSTasks.so");
// change next line to a working code
TChain* chain = NULL; //CreateChain("wn.xml");
gSystem->Load("libANALYSIS");
gSystem->Load("libANALYSISalice"); //AliAnalysisTaskSE
- gSystem->Load("libPWG4UserTasks.so");
+ gSystem->Load("libPWGGAPHOSTasks.so");
gSystem->AddIncludePath("-I$ALICE_ROOT/include -I$ALICE_ROOT/PHOS");
TString outputfile = AliAnalysisManager::GetCommonFileName();
- AliAnalysisDataContainer *coutput = mgr->CreateContainer("omega3pi", TList::Class(), AliAnalysisManager::kOutputContainer, Form("%s:PWG4_omega3pi",outputfile.Data()));
+ AliAnalysisDataContainer *coutput = mgr->CreateContainer("omega3pi", TList::Class(), AliAnalysisManager::kOutputContainer, Form("%s:PWGGA_omega3pi",outputfile.Data()));
// Connect input/output
mgr->ConnectInput(omegaTask, 0, cinput);
//load analysis framework
gSystem->Load("libANALYSIS");
gSystem->Load("libANALYSISalice"); //AliAnalysisTaskSE
- gSystem->Load("libPWG4Omega3pi");
+ gSystem->Load("libPWGGAPHOSTasks");
gSystem->AddIncludePath("-I$ALICE_ROOT/include -I$ALICE_ROOT/PHOS");
// A task can be compiled dynamically with AClic
git://git.uio.no / u / mrichter / AliRoot.git / commitdiff