modifications in Dalitz task
authorfbock <fbock@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 23 Sep 2013 14:34:51 +0000 (14:34 +0000)
committerfbock <fbock@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 23 Sep 2013 14:34:51 +0000 (14:34 +0000)
PWGGA/GammaConv/AliAnalysisTaskGammaConvDalitzV1.cxx
PWGGA/GammaConv/AliAnalysisTaskGammaConvDalitzV1.h
PWGGA/GammaConv/AliDalitzElectronCuts.cxx
PWGGA/GammaConv/macros/AddTask_GammaConvDalitzV1_PbPb.C
PWGGA/GammaConv/macros/AddTask_GammaConvDalitzV1_pPb.C

index 25aea93..2f25e8c 100644 (file)
@@ -127,8 +127,11 @@ fV0Reader(NULL),
    hMCChiCPt(NULL),
    hMCChiCInAccPt(NULL),
    hESDEposEnegTruePi0DalitzInvMassPt(NULL),
+   hESDEposEnegTruePi0DalitzPsiPairDPhi(NULL),
    hESDEposEnegTrueEtaDalitzInvMassPt(NULL),
+   hESDEposEnegTrueEtaDalitzPsiPairDPhi(NULL),
    hESDEposEnegTruePhotonInvMassPt(NULL),
+   hESDEposEnegTruePhotonPsiPairDPhi(NULL),
    hESDEposEnegTrueJPsiInvMassPt(NULL),
    hESDTrueMotherChiCInvMassPt(NULL),
    hESDTrueMotherChiCDiffInvMassPt(NULL),
@@ -247,8 +250,11 @@ AliAnalysisTaskGammaConvDalitzV1::AliAnalysisTaskGammaConvDalitzV1( const char*
    hMCChiCPt(NULL),
    hMCChiCInAccPt(NULL),
    hESDEposEnegTruePi0DalitzInvMassPt(NULL),
+   hESDEposEnegTruePi0DalitzPsiPairDPhi(NULL),
    hESDEposEnegTrueEtaDalitzInvMassPt(NULL),
+   hESDEposEnegTrueEtaDalitzPsiPairDPhi(NULL),
    hESDEposEnegTruePhotonInvMassPt(NULL),
+   hESDEposEnegTruePhotonPsiPairDPhi(NULL),
    hESDEposEnegTrueJPsiInvMassPt(NULL),
    hESDTrueMotherChiCInvMassPt(NULL),
    hESDTrueMotherChiCDiffInvMassPt(NULL),
@@ -677,8 +683,11 @@ void AliAnalysisTaskGammaConvDalitzV1::UserCreateOutputObjects()
                        
       if ( fDoMesonQA ) {
        hESDEposEnegTruePi0DalitzInvMassPt           = new TH2F*[fnCuts];
+        hESDEposEnegTruePi0DalitzPsiPairDPhi         = new TH2F*[fnCuts];
        hESDEposEnegTrueEtaDalitzInvMassPt           = new TH2F*[fnCuts];
+        hESDEposEnegTrueEtaDalitzPsiPairDPhi         = new TH2F*[fnCuts];
        hESDEposEnegTruePhotonInvMassPt              = new TH2F*[fnCuts];
+        hESDEposEnegTruePhotonPsiPairDPhi            = new TH2F*[fnCuts];
        hESDEposEnegTrueJPsiInvMassPt                = new TH2F*[fnCuts];
       }
       
@@ -768,17 +777,34 @@ void AliAnalysisTaskGammaConvDalitzV1::UserCreateOutputObjects()
          fCutFolder[iCut]->Add(fTrueList[iCut]);
 
         if ( fDoMesonQA ) {
+
+
         hESDEposEnegTruePi0DalitzInvMassPt[iCut] = new TH2F("ESD_EposEneg_TruePi0Dalitz_InvMassPt","ESD_EposEneg_TruePi0Dalitz_InvMassPt",5000,0.,5.,100,0.,10.);
         fTrueList[iCut]->Add(hESDEposEnegTruePi0DalitzInvMassPt[iCut]);
+
+         hESDEposEnegTruePi0DalitzPsiPairDPhi[iCut] = new TH2F("ESD_EposEneg_TruePi0Dalitz_PsiPair_DPhi","ESD_EposEneg_TruePi0Dalitz_PsiPair_DPhi", 100, -1.0,1.0,100,-1.0,1.0 );
+         fTrueList[iCut]->Add(hESDEposEnegTruePi0DalitzPsiPairDPhi[iCut]);
+
         
         hESDEposEnegTrueEtaDalitzInvMassPt[iCut] = new TH2F("ESD_EposEneg_TrueEtaDalitz_InvMassPt","ESD_EposEneg_TrueEtaDalitz_InvMassPt",5000,0.,5.,100,0.,10.);
         fTrueList[iCut]->Add(hESDEposEnegTrueEtaDalitzInvMassPt[iCut]);
+
+         hESDEposEnegTrueEtaDalitzPsiPairDPhi[iCut] = new TH2F("ESD_EposEneg_TrueEtaDalitz_PsiPair_DPhi","ESD_EposEneg_TrueEtaDalitz_PsiPair_DPhi", 100, -1.0,1.0,100,-1.0,1.0 );
+         fTrueList[iCut]->Add(hESDEposEnegTrueEtaDalitzPsiPairDPhi[iCut]);
         
         hESDEposEnegTruePhotonInvMassPt[iCut] = new TH2F("ESD_EposEneg_TruePhoton_InvMassPt","ESD_EposEneg_TruePhoton_InvMassPt",5000,0.,5.,100,0.,10.);
         fTrueList[iCut]->Add(hESDEposEnegTruePhotonInvMassPt[iCut]);
-        
+        
+         hESDEposEnegTruePhotonPsiPairDPhi[iCut] = new TH2F("ESD_EposEneg_TruePhoton_PsiPair_DPhi","ESD_EposEneg_TruePhoton_PsiPair_DPhi", 100, -1.0,1.0,100,-1.0,1.0 );
+         fTrueList[iCut]->Add(hESDEposEnegTruePhotonPsiPairDPhi[iCut]);
         hESDEposEnegTrueJPsiInvMassPt[iCut] = new TH2F("ESD_EposEneg_TrueJPsi_InvMassPt","ESD_EposEneg_TrueJPsi_InvMassPt",5000,0.,5.,100,0.,10.);
         fTrueList[iCut]->Add(hESDEposEnegTrueJPsiInvMassPt[iCut]);
+
+
+
+
+
         }
 
 
@@ -1073,7 +1099,7 @@ Bool_t AliAnalysisTaskGammaConvDalitzV1::Notify()
 void AliAnalysisTaskGammaConvDalitzV1::Terminate(const Option_t *)
 {
 
-//Grid
+///Grid
 
 }
 //________________________________________________________________________
@@ -1423,10 +1449,10 @@ void AliAnalysisTaskGammaConvDalitzV1::ProcessElectronCandidates(){
         
          AliAODConversionPhoton *vphoton = new AliAODConversionPhoton(virtualPhoton); //To Apply PsiPairCut
            
-           if ( fDoMesonQA ) {
-             
+         if ( fDoMesonQA ) {
+
              if( fMCEvent ) {
-           
+
              TParticle *mcVgamma=virtualPhoton->GetMCParticle(fMCStack);
 
              if(mcVgamma){
@@ -1441,14 +1467,23 @@ void AliAnalysisTaskGammaConvDalitzV1::ProcessElectronCandidates(){
                else if( mcVgamma->GetPdgCode() == 221 ) isEtaDalitz = kTRUE;
              }
              }
-           
-           
-             if(isPhoton)    hESDEposEnegTruePhotonInvMassPt[fiCut]->Fill(vphoton->GetMass(),vphoton->Pt());
-             else if(isJPsi) hESDEposEnegTrueJPsiInvMassPt[fiCut]->Fill(vphoton->GetMass(),vphoton->Pt());
-             else if(isPi0Dalitz)hESDEposEnegTruePi0DalitzInvMassPt[fiCut]->Fill(vphoton->GetMass(),vphoton->Pt());
-             else if(isEtaDalitz)hESDEposEnegTrueEtaDalitzInvMassPt[fiCut]->Fill(vphoton->GetMass(),vphoton->Pt());
+
+             if(isPhoton){
+                    hESDEposEnegTruePhotonInvMassPt[fiCut]->Fill(vphoton->GetMass(),vphoton->Pt());
+                    hESDEposEnegTruePhotonPsiPairDPhi[fiCut]->Fill(deltaPhi,psiPair);
+              }
+             else if(isJPsi){
+                    hESDEposEnegTrueJPsiInvMassPt[fiCut]->Fill(vphoton->GetMass(),vphoton->Pt());
+              }
+             else if(isPi0Dalitz){
+                    hESDEposEnegTruePi0DalitzInvMassPt[fiCut]->Fill(vphoton->GetMass(),vphoton->Pt());
+                    hESDEposEnegTruePi0DalitzPsiPairDPhi[fiCut]->Fill(deltaPhi,psiPair);
+              }
+             else if(isEtaDalitz){
+                    hESDEposEnegTrueEtaDalitzInvMassPt[fiCut]->Fill(vphoton->GetMass(),vphoton->Pt());
+                    hESDEposEnegTrueEtaDalitzPsiPairDPhi[fiCut]->Fill(deltaPhi,psiPair);
+              }
            }
-           
          }
         
       
index 8815d0e..bb3e2c9 100644 (file)
@@ -152,8 +152,11 @@ class AliAnalysisTaskGammaConvDalitzV1: public AliAnalysisTaskSE
     TH1F **hMCChiCPt;
     TH1F **hMCChiCInAccPt;
     TH2F **hESDEposEnegTruePi0DalitzInvMassPt;
+    TH2F **hESDEposEnegTruePi0DalitzPsiPairDPhi;
     TH2F **hESDEposEnegTrueEtaDalitzInvMassPt;
+    TH2F **hESDEposEnegTrueEtaDalitzPsiPairDPhi;
     TH2F **hESDEposEnegTruePhotonInvMassPt;
+    TH2F **hESDEposEnegTruePhotonPsiPairDPhi;
     TH2F **hESDEposEnegTrueJPsiInvMassPt;
     TH2F **hESDTrueMotherChiCInvMassPt;
     TH2F **hESDTrueMotherChiCDiffInvMassPt;
index b94406b..f4a59c8 100644 (file)
@@ -1535,13 +1535,19 @@ Bool_t AliDalitzElectronCuts::SetPsiPairCut(Int_t psiCut) {
   case 1:\r
         fDoPsiPairCut = kTRUE;\r
         fPsiPairCut = 0.45; // Standard\r
-        fDeltaPhiCutMin = 0.;\r
+        fDeltaPhiCutMin = 0.0;\r
         fDeltaPhiCutMax = 0.12;\r
         break;\r
   case 2:\r
        fDoPsiPairCut = kTRUE;\r
         fPsiPairCut = 0.60; \r
-        fDeltaPhiCutMin = 0.;\r
+        fDeltaPhiCutMin = 0.0;\r
+        fDeltaPhiCutMax = 0.12;\r
+        break;\r
+  case 3:\r
+        fDoPsiPairCut = kTRUE;\r
+        fPsiPairCut = 0.52;\r
+        fDeltaPhiCutMin = 0.0;\r
         fDeltaPhiCutMax = 0.12;\r
         break;\r
     \r
index 94d35a1..d98c53d 100644 (file)
@@ -1,4 +1,5 @@
-void AddTask_GammaConvDalitzV1_PbPb(   Bool_t isMC   = kFALSE, //run MC 
+void AddTask_GammaConvDalitzV1_PbPb(   Int_t trainConfig = 1,
+                                       Bool_t isMC   = kFALSE, //run MC 
                                        Bool_t enableQAMesonTask = kTRUE, //enable QA in AliAnalysisTaskGammaConvDalitzV1
                                        Bool_t enableDoMesonChic = kFALSE, // enable additional Chic analysis
                                        TString fileNameInputForWeighting = "MCSpectraInput.root", // path to file for weigting input
@@ -35,7 +36,7 @@ void AddTask_GammaConvDalitzV1_PbPb(   Bool_t isMC   = kFALSE, //run MC
    // ================== GetAnalysisManager ===============================
    AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
    if (!mgr) {
-      Error("AddTask_GammaConvDalitzV1_PbPb", "No analysis manager found.");
+      Error(Form("AddTask_GammaConvDalitzV1_PbPb_%i",trainConfig), "No analysis manager found.");
       return ;
    }
 
@@ -100,40 +101,27 @@ void AddTask_GammaConvDalitzV1_PbPb(   Bool_t isMC   = kFALSE, //run MC
 
    if( !(AliDalitzElectronSelector*)mgr->GetTask("ElectronSelector") ){
 
-   AliDalitzElectronSelector *fElectronSelector = new AliDalitzElectronSelector("ElectronSelector");
+      AliDalitzElectronSelector *fElectronSelector = new AliDalitzElectronSelector("ElectronSelector");
 
-   // Set AnalysisCut Number
-
-   AliDalitzElectronCuts *fElecCuts=0;
-
-   TString ElecCuts = "900054000000020000";
-
-
-
-    if( ElecCuts!=""){
-
-       fElecCuts= new AliDalitzElectronCuts(ElecCuts.Data(),ElecCuts.Data());
-
-            if(fElecCuts->InitializeCutsFromCutString(ElecCuts.Data())){
-
-                fElectronSelector->SetDalitzElectronCuts(fElecCuts);
-
-                fElecCuts->SetFillCutHistograms("",kTRUE);
-
-            }
-
-    }
-
-    fElectronSelector->Init();
-    mgr->AddTask(fElectronSelector);
-
-    AliAnalysisDataContainer *cinput1  = mgr->GetCommonInputContainer();
-
-    //connect input V0Reader
+      // Set AnalysisCut Number
+      AliDalitzElectronCuts *fElecCuts=0;
+      TString ElecCuts = "900054000000020000";
+
+      if( ElecCuts!=""){
+         fElecCuts= new AliDalitzElectronCuts(ElecCuts.Data(),ElecCuts.Data());
+         if(fElecCuts->InitializeCutsFromCutString(ElecCuts.Data())){
+            fElectronSelector->SetDalitzElectronCuts(fElecCuts);
+            fElecCuts->SetFillCutHistograms("",kTRUE);
+         }
+      }
 
-    mgr->ConnectInput (fElectronSelector,0,cinput1);
+      fElectronSelector->Init();
+      mgr->AddTask(fElectronSelector);
 
-    }
+      AliAnalysisDataContainer *cinput1  = mgr->GetCommonInputContainer();
+      //connect input V0Reader
+      mgr->ConnectInput (fElectronSelector,0,cinput1);
+   }
 
 
 
@@ -146,35 +134,28 @@ void AddTask_GammaConvDalitzV1_PbPb(   Bool_t isMC   = kFALSE, //run MC
   
  
    AliAnalysisTaskGammaConvDalitzV1 *task=NULL;
-
-   task= new AliAnalysisTaskGammaConvDalitzV1("GammaConvDalitzV1");
-
-   task->SetIsHeavyIon(2);
+   task= new AliAnalysisTaskGammaConvDalitzV1(Form("GammaConvDalitzV1_%i",trainConfig));
+   task->SetIsHeavyIon(1);
    task->SetIsMC(isMC);
 
 
 
    // Cut Numbers to use in Analysis
-   Int_t numberOfCuts = 5;
+   Int_t numberOfCuts = 3;
 
    TString *ConvCutarray    = new TString[numberOfCuts];
-
    TString *ElecCutarray    = new TString[numberOfCuts];
-
    TString *MesonCutarray   = new TString[numberOfCuts];
 
-
-      
-
-
-   ConvCutarray[0]  = "312000104209297100322000000"; MesonCutarray[0] = "01522045009000"; ElecCutarray[0]  = "905162002552021217"; //PbPb  5-10% kAny
-   ConvCutarray[1]  = "112000104209297100322000000"; MesonCutarray[1] = "01522045009000"; ElecCutarray[1]  = "905162002552021217"; //PbPb 10-20% kAny
-   ConvCutarray[2]  = "124000104209297100322000000"; MesonCutarray[2] = "01522045009000"; ElecCutarray[2]  = "905162002552021217"; //PbPb 20-40% kAny
-   ConvCutarray[3]  = "146000104209297100322000000"; MesonCutarray[3] = "01522045009000"; ElecCutarray[3]  = "905162002552021217"; //PbPb 40-60% kAny
-   ConvCutarray[4]  = "168000104209297100322000000"; MesonCutarray[4] = "01522045009000"; ElecCutarray[4]  = "905162002552021217"; //PbPb 60-80% kAny
-
-
-
+   if( trainConfig == 1 ) {
+        ConvCutarray[0]  = "124000104209297100322000000"; MesonCutarray[0] = "01522045009000"; ElecCutarray[0]  = "905162002553025217"; //PbPb 20-40% kAny
+        ConvCutarray[1]  = "146000104209297100322000000"; MesonCutarray[1] = "01522045009000"; ElecCutarray[1]  = "905162002553025217"; //PbPb 40-60% kAny
+        ConvCutarray[2]  = "168000104209297100322000000"; MesonCutarray[2] = "01522045009000"; ElecCutarray[2]  = "905162002553025217"; //PbPb 60-80% kAny
+   } else if ( trainConfig == 2 ) {
+        ConvCutarray[0]  = "524000104209297100322000000"; MesonCutarray[0] = "01522045009000"; ElecCutarray[0]  = "905162002553025217"; //PbPb 20-40% kAny
+        ConvCutarray[1]  = "546000104209297100322000000"; MesonCutarray[1] = "01522045009000"; ElecCutarray[1]  = "905162002553025217"; //PbPb 40-60% kAny
+        ConvCutarray[2]  = "568000104209297100322000000"; MesonCutarray[2] = "01522045009000"; ElecCutarray[2]  = "905162002553025217"; //PbPb 60-80% kAny
+   }
 
    TList *ConvCutList  = new TList();
    TList *MesonCutList = new TList();
@@ -199,58 +180,43 @@ void AddTask_GammaConvDalitzV1_PbPb(   Bool_t isMC   = kFALSE, //run MC
 
 
    for(Int_t i = 0; i<numberOfCuts; i++){
-
-
       analysisCuts[i] = new AliConversionCuts();
       if( ! analysisCuts[i]->InitializeCutsFromCutString(ConvCutarray[i].Data()) ) {
             cout<<"ERROR: analysisCuts [" <<i<<"]"<<endl;
             return 0;
+      } else {
+         ConvCutList->Add(analysisCuts[i]);
+         analysisCuts[i]->SetFillCutHistograms("",kFALSE);
+         if( trainConfig == 1){
+            if (i == 0 && doWeighting) analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_2040V0M", "", "","Pi0_Fit_Data_PbPb_2760GeV_2040V0M");
+            if (i == 1 && doWeighting) analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_4060V0M", "", "","Pi0_Fit_Data_PbPb_2760GeV_4060V0M");
+            if (i == 2 && doWeighting) analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_6080V0M", "", "","Pi0_Fit_Data_PbPb_2760GeV_6080V0M");
+         } else if ( trainConfig == 2 ) {
+            if (i == 0 && doWeighting) analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_2040TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_2040V0M");
+            if (i == 1 && doWeighting) analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_4060TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_4060V0M");
+            if (i == 2 && doWeighting) analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_6080TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_6080V0M");
+         }
       }
-      else {
-      
-      ConvCutList->Add(analysisCuts[i]);
-      analysisCuts[i]->SetFillCutHistograms("",kFALSE);
-      if (i == 0 && doWeighting) analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_0510V0M", "", "","Pi0_Fit_Data_PbPb_2760GeV_0510V0M");
-      if (i == 1 && doWeighting) analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_1020V0M", "", "","Pi0_Fit_Data_PbPb_2760GeV_1020V0M");
-      if (i == 2 && doWeighting) analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_2040V0M", "", "","Pi0_Fit_Data_PbPb_2760GeV_2040V0M");
-      if (i == 3 && doWeighting) analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_4060V0M", "", "","Pi0_Fit_Data_PbPb_2760GeV_4060V0M");
-      if (i == 4 && doWeighting) analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_6080V0M", "", "","Pi0_Fit_Data_PbPb_2760GeV_6080V0M");
-
-
-
-      }
-
-
 
       analysisMesonCuts[i] = new AliConversionMesonCuts();
-    
       if( ! analysisMesonCuts[i]->InitializeCutsFromCutString(MesonCutarray[i].Data()) ) {
             cout<<"ERROR: analysisMesonCuts [ " <<i<<" ] "<<endl;
             return 0;
-      }
-      else {
+      } else {
         MesonCutList->Add(analysisMesonCuts[i]);
         analysisMesonCuts[i]->SetFillCutHistograms("");
       }
-
-
-       TString cutName( Form("%s_%s_%s",ConvCutarray[i].Data(),ElecCutarray[i].Data(),MesonCutarray[i].Data() ) );
-
-
-       analysisElecCuts[i] = new AliDalitzElectronCuts();
-       if( !analysisElecCuts[i]->InitializeCutsFromCutString(ElecCutarray[i].Data())) {
-
+      
+      TString cutName( Form("%s_%s_%s",ConvCutarray[i].Data(),ElecCutarray[i].Data(),MesonCutarray[i].Data() ) );
+      analysisElecCuts[i] = new AliDalitzElectronCuts();
+      if( !analysisElecCuts[i]->InitializeCutsFromCutString(ElecCutarray[i].Data())) {
             cout<< "ERROR:  analysisElecCuts [ " <<i<<" ] "<<endl;
             return 0;
-       }
-       else { 
-        ElecCutList->Add(analysisElecCuts[i]);
-        analysisElecCuts[i]->SetFillCutHistograms("",kFALSE,cutName); 
-       }
-
-
-        analysisCuts[i]->SetAcceptedHeader(HeaderList);
+      }  else { 
+         ElecCutList->Add(analysisElecCuts[i]);
+         analysisElecCuts[i]->SetFillCutHistograms("",kFALSE,cutName); 
+      }
+      analysisCuts[i]->SetAcceptedHeader(HeaderList);
 
    }
 
@@ -267,8 +233,8 @@ void AddTask_GammaConvDalitzV1_PbPb(   Bool_t isMC   = kFALSE, //run MC
 
    //connect containers
    AliAnalysisDataContainer *coutput =
-   mgr->CreateContainer("GammaConvDalitzV1", TList::Class(),
-                           AliAnalysisManager::kOutputContainer,"GammaConvV1.root");
+   mgr->CreateContainer(Form("GammaConvDalitzV1_%i",trainConfig), TList::Class(),
+                           AliAnalysisManager::kOutputContainer,Form("GammaConvV1Dalitz_%i.root",trainConfig));
 
    mgr->AddTask(task);
    mgr->ConnectInput(task,0,cinput);
index 95ebad9..c375982 100644 (file)
@@ -169,16 +169,12 @@ void AddTask_GammaConvDalitzV1_pPb(   Bool_t isMC       = kFALSE, //run MC
    stringShift = "pPb";
 
    ConvCutarray[0] = "800000008209360300220000000"; ElecCutarray[0] = "904754002582026217"; MesonCutarray[0] = "01039035009000"; //standard cut Pi0 PbPb 00-100
-   //ConvCutarray[1] = "800000008209360300220000000"; ElecCutarray[1] = "904754002582026217"; MesonCutarray[1] = "01039035000000"; //standard cut Pi0 PbPb 00-100 + No extra smearing
    ConvCutarray[1] = "800000008209360300220000000"; ElecCutarray[1] = "904754002582026117"; MesonCutarray[1] = "01039035009000"; //standard cut Pi0 PbPb 00-100 + Single Pt primary > 0.100 GeV
-   //ConvCutarray[3] = "800000008209360300220000000"; ElecCutarray[3] = "904754002582026317"; MesonCutarray[3] = "01039035009000"; //standard cut Pi0 PbPb 00-100 + Single Pt primary > 0.150 GeV
    ConvCutarray[2] = "800000008209460300220000000"; ElecCutarray[2] = "904754002582026217"; MesonCutarray[2] = "01039035009000"; //standard cut Pi0 PbPb 00-100 + dEdx electron gamma   -6 ,7 sigmas
    ConvCutarray[3] = "800000008209360300220300000"; ElecCutarray[3] = "904754002582026217"; MesonCutarray[3] = "01039035009000"; //standard cut Pi0 PbPb 00-100  do Aysemtri cut
-   ConvCutarray[4] = "800000008209460300220000000"; ElecCutarray[4] = "902754002582026217"; MesonCutarray[4] = "01039035009000"; //standard cut Pi0 PbPb 00-100  dEdx electrons primary and secondary -6 7 sigmas
+   ConvCutarray[4] = "800000008209360300220000000"; ElecCutarray[4] = "905154002582026217"; MesonCutarray[4] = "01039035009000"; //standard cut Pi0 PbPb 00-100
+  
    
-   //ConvCutarray[7] = "800000008249360300220000000"; ElecCutarray[7] = "904754002582026217"; MesonCutarray[7] = "01039035009000"; //standard cut Pi0 PbPb 00-100  Single Pt gamma > 0.075 GeV
-   //ConvCutarray[8] = "800000008219360300220000000"; ElecCutarray[8] = "904754002582026217"; MesonCutarray[8] = "01039035009000"; //standard cut Pi0 PbPb 00-100  Single Pt gamma > 0.100 GeV
-   //ConvCutarray[9] = "800000008209360300220000000"; ElecCutarray[9] = "904854002582026217"; MesonCutarray[9] = "01039035009000"; //standard cut Pi0 PbPb 00-100
  
 
    TList *ConvCutList  = new TList();
@@ -266,7 +262,7 @@ void AddTask_GammaConvDalitzV1_pPb(   Bool_t isMC       = kFALSE, //run MC
    //connect containers
    AliAnalysisDataContainer *coutput =
    mgr->CreateContainer("GammaConvDalitzV1", TList::Class(),
-                           AliAnalysisManager::kOutputContainer,"GammaConvV1.root");
+                           AliAnalysisManager::kOutputContainer,"GammaConvV1Dalitz.root");
 
    mgr->AddTask(task);
    mgr->ConnectInput(task,0,cinput);