AliAnalysisTask *AddTask_ConversionAODProduction(Int_t dataset=1){
+ // Before doing anything, we load the needed library
+ gSystem->Load("libPWGGAGammaConv.so");
// dataset 0: pp
// dataset 1: PbPb
//get the current analysis manager
AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
if (!mgr) {
- Error("AddTask_V0ReaderV1", "No analysis manager found.");
- return 0;
+ Error("AddTask_V0ReaderV1", "No analysis manager found.");
+ return 0;
}
TString analysiscut;
if(IsHeavyIon){
- analysiscut="900177009350113211200001000000000";
+ // Old cut string, no longer compatible with AliConversionCuts
+ // analysiscut="900177009350113211200001000000000";
+ // New cut string as of April 2013
+ analysiscut="1000000032091071001000000";
}
else{
- analysiscut="900397209450304221200000002000000";
+ // analysiscut="900397209450304221200000002000000";
+ analysiscut="0000000002084000002200000";
}
//========= Add V0 Reader to ANALYSIS manager =====
// Set AnalysisCut Number
AliConversionCuts *fCuts= new AliConversionCuts(analysiscut.Data(),analysiscut.Data());
if(fCuts->InitializeCutsFromCutString(analysiscut.Data())){
- fV0Reader->SetConversionCuts(fCuts);
+ fV0Reader->SetConversionCuts(fCuts);
}
fV0Reader->Init();
mgr->AddTask(fV0Reader);
else if (trainConfig.Contains("pPb")) IsHeavyIon=2;
Bool_t doEtaShift = kFALSE;
- Int_t forceEtaShift = 0; // Carefull !!! Should be zero otherwise will shift eta cut for all periods
- // Use doEtaShift flag for pPb or Pbp instead (1: shift +0.465, 2: shift -0.465)
-
+
TString cutnumber = "";
if(IsHeavyIon == 1){
cutnumber = "1000000002084001001500000";
fCuts->SelectCollisionCandidates(AliVEvent::kINT7);
fCuts->DoEtaShift(doEtaShift);
}
- fCuts->ForceEtaShift(forceEtaShift);
fV0ReaderV1->SetConversionCuts(fCuts);
fCuts->SetFillCutHistograms("",kTRUE);
}
cutarray[ 0] = "8020000082093172023290000"; mesonCutArray[0] = "01629045000"; //standard cut Pi0 Pb 00-20 shifted Eta 0.4
cutarray[ 1] = "8240000082093172023290000"; mesonCutArray[1] = "01629045000"; //standard cut Pi0 Pb 20-40 shifted Eta 0.4
cutarray[ 2] = "8460000082093172023290000"; mesonCutArray[2] = "01629045000"; //standard cut Pi0 Pb 40-60 shifted Eta 0.4
- cutarray[ 3] = "8680000082093172023290000"; mesonCutArray[3] = "01629045000"; //standard cut Pi0 Pb 60-80 shifted Eta 0.4
+ cutarray[ 3] = "8600000082093172023290000"; mesonCutArray[3] = "01629045000"; //standard cut Pi0 Pb 60-100 shifted Eta 0.4
} else {
cutarray[ 0] = "0000012002093663003800000"; mesonCutArray[0] = "01631031009"; //standard cut Pi0 pp 2.76TeV without SDD , only boxes
cutarray[ 1] = "0001012002093663003800000"; mesonCutArray[1] = "01631031009"; //standard cut Pi0 pp 2.76TeV without SDD, V0AND , only boxes
else if (trainConfig.Contains("pPb")) IsHeavyIon=2;
Bool_t doEtaShift = kFALSE;
- Int_t forceEtaShift = 0; // Carefull !!! Should be zero otherwise will shift eta cut for all periods
- // Use doEtaShift flag for pPb or Pbp instead (1: shift +0.465, 2: shift -0.465)
-
+
TString cutnumber = "";
if(IsHeavyIon == 1){
cutnumber = "1000000002084001001500000";
fCuts->SelectCollisionCandidates(AliVEvent::kINT7);
fCuts->DoEtaShift(doEtaShift);
}
- fCuts->ForceEtaShift(forceEtaShift);
fV0ReaderV1->SetConversionCuts(fCuts);
fCuts->SetFillCutHistograms("",kTRUE);
}
else if (trainConfig.Contains("pPb")) IsHeavyIon=2;
Bool_t doEtaShift = kFALSE;
- Int_t forceEtaShift = 0; // Carefull !!! Should be zero otherwise will shift eta cut for all periods
- // Use doEtaShift flag for pPb or Pbp instead (1: shift +0.465, 2: shift -0.465)
TString cutnumber = "";
if(IsHeavyIon == 1){
fCuts->SelectCollisionCandidates(AliVEvent::kINT7);
fCuts->DoEtaShift(doEtaShift);
}
- fCuts->ForceEtaShift(forceEtaShift);
fV0ReaderV1->SetConversionCuts(fCuts);
fCuts->SetFillCutHistograms("",kTRUE);
}
cutarray[ 0] = "8020000072093172023290000"; mesonCutArray[0] = "01627045000"; //standard cut Pi0 Pb 00-20 wo shifted Eta 0.3
cutarray[ 1] = "8240000072093172023290000"; mesonCutArray[1] = "01627045000"; //standard cut Pi0 Pb 20-40 wo shifted Eta 0.3
cutarray[ 2] = "8460000072093172023290000"; mesonCutArray[2] = "01627045000"; //standard cut Pi0 Pb 40-60 wo shifted Eta 0.3
- cutarray[ 3] = "8680000072093172023290000"; mesonCutArray[3] = "01627045000"; //standard cut Pi0 Pb 60-80 wo shifted Eta 0.3
+ cutarray[ 3] = "8600000072093172023290000"; mesonCutArray[3] = "01627045000"; //standard cut Pi0 Pb 60-100 wo shifted Eta 0.3
} else {
cutarray[ 0] = "0002011002093663003800000"; mesonCutArray[0] = "01631031009"; //standard cut Pi0 pp 2.76TeV with SDD , only Minbias MC
cutarray[ 1] = "0003011002093663003800000"; mesonCutArray[1] = "01631031009"; //standard cut Pi0 pp 2.76TeV with SDD, V0AND , only Minbias MC