submitted AddTasks for GG pPb
authorfbock <fbock@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 19 Sep 2013 14:12:48 +0000 (14:12 +0000)
committerfbock <fbock@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 19 Sep 2013 14:12:48 +0000 (14:12 +0000)
PWGGA/GammaConv/macros/AddTask_GammaConvV1_pPb.C

index 39bee85..a9e1266 100644 (file)
@@ -1,11 +1,11 @@
 void AddTask_GammaConvV1_pPb(  Int_t trainConfig = 1,  //change different set of cuts
-                              Bool_t isMC   = kFALSE, //run MC 
+                              Bool_t isMC   = kFALSE, //run MC
                               Bool_t enableQAMesonTask = kFALSE, //enable QA in AliAnalysisTaskGammaConvV1
                               Bool_t enableQAPhotonTask = kFALSE, // enable additional QA task
                               TString fileNameInputForWeighting = "MCSpectraInput.root", // path to file for weigting input
                               TString cutnumberAODBranch = "8000000060084000001500000" // cutnumber for AOD branch
                            ) {
-
    // ================= Load Librariers =================================
    gSystem->Load("libCore.so");  
    gSystem->Load("libTree.so");
@@ -24,14 +24,14 @@ void AddTask_GammaConvV1_pPb(  Int_t trainConfig = 1,  //change different set of
    gSystem->Load("libSTEERBase.so");
    gSystem->Load("libTENDER.so");
    gSystem->Load("libTENDERSupplies.so");
-      
+     
    // ================== GetAnalysisManager ===============================
    AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
    if (!mgr) {
       Error(Form("AddTask_GammaConvV1_%i",trainConfig), "No analysis manager found.");
       return ;
    }
-
    // ================== GetInputEventHandler =============================
    AliVEventHandler *inputHandler=mgr->GetInputEventHandler();
    
@@ -42,22 +42,22 @@ void AddTask_GammaConvV1_pPb(  Int_t trainConfig = 1,  //change different set of
    }
    
    //=========  Set Cutnumber for V0Reader ================================
-   TString cutnumber = "800000006008400100150000000"; 
+   TString cutnumber = "800000006008400100150000000";
    Bool_t doEtaShift = kFALSE;
    AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
    //========= Add V0 Reader to  ANALYSIS manager if not yet existent =====
    if( !(AliV0ReaderV1*)mgr->GetTask("V0ReaderV1") ){
       AliV0ReaderV1 *fV0ReaderV1 = new AliV0ReaderV1("V0ReaderV1");
-      
+     
       fV0ReaderV1->SetUseOwnXYZCalculation(kTRUE);
       fV0ReaderV1->SetCreateAODs(kFALSE);// AOD Output
       fV0ReaderV1->SetUseAODConversionPhoton(kTRUE);
-      
+     
       if (!mgr) {
          Error("AddTask_V0ReaderV1", "No analysis manager found.");
          return;
       }
-
       // Set AnalysisCut Number
       AliConversionCuts *fCuts=NULL;
       if(cutnumber!=""){
@@ -74,15 +74,15 @@ void AddTask_GammaConvV1_pPb(  Int_t trainConfig = 1,  //change different set of
          fV0ReaderV1->SetDeltaAODBranchName(Form("GammaConv_%s_gamma",cutnumberAODBranch.Data()));
       }
       fV0ReaderV1->Init();
-
       AliLog::SetGlobalLogLevel(AliLog::kInfo);
-
       //connect input V0Reader
       mgr->AddTask(fV0ReaderV1);
       mgr->ConnectInput(fV0ReaderV1,0,cinput);
-
    }
-
    //================================================
    //========= Add task to the ANALYSIS manager =====
    //================================================
@@ -93,65 +93,66 @@ void AddTask_GammaConvV1_pPb(  Int_t trainConfig = 1,  //change different set of
    task->SetIsMC(isMC);
    // Cut Numbers to use in Analysis
    Int_t numberOfCuts = 6;
-
    TString *cutarray = new TString[numberOfCuts];
    TString *mesonCutArray = new TString[numberOfCuts];
    Bool_t doEtaShiftIndCuts = kFALSE;
    TString stringShift = "";
-   if(trainConfig == 1){ 
+   if(trainConfig == 1){
       // Shifting in pPb direction
       doEtaShiftIndCuts = kTRUE;
       stringShift = "pPb";
-      cutarray[ 0] = "800000008009317200329000000"; mesonCutArray[ 0] = "01629045000000";  //standard cut Pi0 pPb RCut Changed
-      cutarray[ 1] = "800000008109317200329000000"; mesonCutArray[ 1] = "01629045000000";  //standard cut Pi0 pPb RCut Changed
-      cutarray[ 2] = "800000008309317200329000000"; mesonCutArray[ 2] = "01629045000000";  //standard cut Pi0 pPb RCut Changed
-      cutarray[ 3] = "800000008409317200329000000"; mesonCutArray[ 3] = "01629045000000";  //standard cut Pi0 pPb RCut Changed
-      cutarray[ 4] = "800000008509317200329000000"; mesonCutArray[ 4] = "01629045000000";  //standard cut Pi0 pPb RCut Changed
-      cutarray[ 5] = "800000008609317200329000000"; mesonCutArray[ 5] = "01629045000000";  //standard cut Pi0 pPb RCut Changed
-   } else if (trainConfig == 2) { 
-      doEtaShiftIndCuts = kTRUE;
-      stringShift = "pPb";
-      cutarray[ 0] = "800000008219317200329000000"; mesonCutArray[ 0] = "01629045000000";  //standard cut Pi0 pPb Single pT Cut changed
-      cutarray[ 1] = "800000008229317200329000000"; mesonCutArray[ 1] = "01629045000000";  //standard cut Pi0 pPb Single pT Cut changed
-      cutarray[ 2] = "800000008249317200329000000"; mesonCutArray[ 2] = "01629045000000";  //standard cut Pi0 pPb Single pT Cut changed
-      cutarray[ 3] = "800000008269317200329000000"; mesonCutArray[ 3] = "01629045000000";  //standard cut Pi0 pPb Single pT Cut changed
-      cutarray[ 4] = "800000008279317200329000000"; mesonCutArray[ 4] = "01629045000000";  //standard cut Pi0 pPb Single pT Cut changed
-      cutarray[ 5] = "800000008259317200329000000"; mesonCutArray[ 5] = "01629045000000";  //standard cut Pi0 pPb Single pT Cut changed
-   } else if (trainConfig == 3) { 
-      doEtaShiftIndCuts = kTRUE;
-      stringShift = "pPb";
-      cutarray[ 0] = "800000008209017200329000000"; mesonCutArray[ 0] = "01629045000000";  //standard cut Pi0 pPb dedx electron Cut changed
-      cutarray[ 1] = "800000008209117200329000000"; mesonCutArray[ 1] = "01629045000000";  //standard cut Pi0 pPb  dedx electron Cut changed
-      cutarray[ 2] = "800000008209217200329000000"; mesonCutArray[ 2] = "01629045000000";  //standard cut Pi0 pPb  dedx electron Cut changed
-      cutarray[ 3] = "800000008209417200329000000"; mesonCutArray[ 3] = "01629045000000";  //standard cut Pi0 pPb  dedx electron Cut changed
-      cutarray[ 4] = "800000008209517200329000000"; mesonCutArray[ 4] = "01629045000000";  //standard cut Pi0 pPb  dedx electron Cut changed
-      cutarray[ 5] = "800000008209617200329000000"; mesonCutArray[ 5] = "01629045000000";  //standard cut Pi0 pPb  dedx electron Cut changed
-   } else if (trainConfig == 4) { 
-      doEtaShiftIndCuts = kTRUE;
-      stringShift = "pPb";
-      cutarray[ 0] = "800000008209316200329000000"; mesonCutArray[ 0] = "01629045000000";  //standard cut Pi0 pPb dedx Pion Line changed
-      cutarray[ 1] = "800000008209310200329000000"; mesonCutArray[ 1] = "01629045000000";  //standard cut Pi0 pPb dedx Pion Line changed
-      cutarray[ 2] = "800000008209317400329000000"; mesonCutArray[ 2] = "01629045000000";  //standard cut Pi0 pPb dedx Pion Line changed
-      cutarray[ 3] = "800000008209317300329000000"; mesonCutArray[ 3] = "01629045000000";  //standard cut Pi0 pPb dedx Pion Line changed
-      cutarray[ 4] = "800000008209357200329000000"; mesonCutArray[ 4] = "01629045000000";  //standard cut Pi0 pPb dedx Pion Line changed
-      cutarray[ 5] = "800000008209307200329000000"; mesonCutArray[ 5] = "01629045000000";  //standard cut Pi0 pPb dedx Pion Line changed
-   } else if (trainConfig == 5) { 
-      doEtaShiftIndCuts = kFALSE;
-      stringShift = "pPb"; //GammaCuts
-      cutarray[ 0] = "800001100209326000380000000"; mesonCutArray[ 0] = "01631031009000"; 
-      cutarray[ 1] = "802001100209326000380000000"; mesonCutArray[ 1] = "01631031009000"; 
-      cutarray[ 2] = "804001100209326000380000000"; mesonCutArray[ 2] = "01631031009000"; 
-      cutarray[ 3] = "828001100209326000380000000"; mesonCutArray[ 3] = "01631031009000"; 
-      cutarray[ 4] = "848001100209326000380000000"; mesonCutArray[ 4] = "01631031009000"; 
-      cutarray[ 5] = "800001100209226000380000000"; mesonCutArray[ 5] = "01631031009000"; 
-   } else {
+      cutarray[ 0] = "800000008209317200329000000"; mesonCutArray[ 0] = "01629045009000";  //standard cut Pi0 pPb RCut Changed
+      cutarray[ 1] = "802000008209317200329000000"; mesonCutArray[ 1] = "01629045009000";  //standard cut Pi0 pPb RCut Changed
+      cutarray[ 2] = "824000008209317200329000000"; mesonCutArray[ 2] = "01629045009000";  //standard cut Pi0 pPb RCut Changed
+      cutarray[ 3] = "846000008209317200329000000"; mesonCutArray[ 3] = "01629045009000";  //standard cut Pi0 pPb RCut Changed
+      cutarray[ 4] = "868000008209317200329000000"; mesonCutArray[ 4] = "01629045009000";  //standard cut Pi0 pPb RCut Changed
+      cutarray[ 5] = "860000008209317200329000000"; mesonCutArray[ 5] = "01629045009000";  //standard cut Pi0 pPb RCut Changed
+   } else if (trainConfig == 2) {
+     // doEtaShiftIndCuts = kTRUE;
+     // stringShift = "pPb";
+      cutarray[ 0] = "800000007209317200329000000"; mesonCutArray[ 0] = "01627045009000";  //standard cut Pi0 pPb Single pT Cut changed
+      cutarray[ 1] = "802000007209317200329000000"; mesonCutArray[ 1] = "01627045009000";  //standard cut Pi0 pPb Single pT Cut changed
+      cutarray[ 2] = "824000007209317200329000000"; mesonCutArray[ 2] = "01627045009000";  //standard cut Pi0 pPb Single pT Cut changed
+      cutarray[ 3] = "846000007209317200329000000"; mesonCutArray[ 3] = "01627045009000";  //standard cut Pi0 pPb Single pT Cut changed
+      cutarray[ 4] = "868000007209317200329000000"; mesonCutArray[ 4] = "01627045009000";  //standard cut Pi0 pPb Single pT Cut changed
+      cutarray[ 5] = "860000007209317200329000000"; mesonCutArray[ 5] = "01627045009000";  //standard cut Pi0 pPb Single pT Cut changed
+   } else if (trainConfig == 3) {
+     // doEtaShiftIndCuts = kTRUE;
+     //  stringShift = "pPb";
+      cutarray[ 0] = "800000000209317200329000000"; mesonCutArray[ 0] = "01620045009000";  //standard cut Pi0 pPb Single pT Cut changed
+      cutarray[ 1] = "802000000209317200329000000"; mesonCutArray[ 1] = "01620045009000";  //standard cut Pi0 pPb Single pT Cut changed
+      cutarray[ 2] = "824000000209317200329000000"; mesonCutArray[ 2] = "01620045009000";  //standard cut Pi0 pPb Single pT Cut changed
+      cutarray[ 3] = "846000000209317200329000000"; mesonCutArray[ 3] = "01620045009000";  //standard cut Pi0 pPb Single pT Cut changed
+      cutarray[ 4] = "868000000209317200329000000"; mesonCutArray[ 4] = "01620045009000";  //standard cut Pi0 pPb Single pT Cut changed
+      cutarray[ 5] = "860000000209317200329000000"; mesonCutArray[ 5] = "01620045009000";  //standard cut Pi0 pPb Single pT Cut changed
+   }//  else if (trainConfig == 4) {
+   //    doEtaShiftIndCuts = kTRUE;
+   //    stringShift = "pPb";
+   //    cutarray[ 0] = "800000008209316200329000000"; mesonCutArray[ 0] = "01629045000000";  //standard cut Pi0 pPb dedx Pion Line changed
+   //    cutarray[ 1] = "800000008209310200329000000"; mesonCutArray[ 1] = "01629045000000";  //standard cut Pi0 pPb dedx Pion Line changed
+   //    cutarray[ 2] = "800000008209317400329000000"; mesonCutArray[ 2] = "01629045000000";  //standard cut Pi0 pPb dedx Pion Line changed
+   //    cutarray[ 3] = "800000008209317300329000000"; mesonCutArray[ 3] = "01629045000000";  //standard cut Pi0 pPb dedx Pion Line changed
+   //    cutarray[ 4] = "800000008209357200329000000"; mesonCutArray[ 4] = "01629045000000";  //standard cut Pi0 pPb dedx Pion Line changed
+   //    cutarray[ 5] = "800000008209307200329000000"; mesonCutArray[ 5] = "01629045000000";  //standard cut Pi0 pPb dedx Pion Line changed
+   // } else if (trainConfig == 5) {
+   //    doEtaShiftIndCuts = kFALSE;
+   //    stringShift = "pPb"; //GammaCuts
+   //    cutarray[ 0] = "800001100209326000380000000"; mesonCutArray[ 0] = "01631031009000";
+   //    cutarray[ 1] = "802001100209326000380000000"; mesonCutArray[ 1] = "01631031009000";
+   //    cutarray[ 2] = "804001100209326000380000000"; mesonCutArray[ 2] = "01631031009000";
+   //    cutarray[ 3] = "828001100209326000380000000"; mesonCutArray[ 3] = "01631031009000";
+   //    cutarray[ 4] = "848001100209326000380000000"; mesonCutArray[ 4] = "01631031009000";
+   //    cutarray[ 5] = "800001100209226000380000000"; mesonCutArray[ 5] = "01631031009000";
+   // }
+   else {
       Error(Form("GammaConvV1_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration");
       return;
    }
-
    TList *ConvCutList = new TList();
    TList *MesonCutList = new TList();
-
    TList *HeaderList = new TList();
    TObjString *Header1 = new TObjString("pi0_1");
    HeaderList->Add(Header1);
@@ -162,17 +163,17 @@ void AddTask_GammaConvV1_pPb(  Int_t trainConfig = 1,  //change different set of
    AliConversionCuts **analysisCuts = new AliConversionCuts*[numberOfCuts];
    MesonCutList->SetOwner(kTRUE);
    AliConversionMesonCuts **analysisMesonCuts = new AliConversionMesonCuts*[numberOfCuts];
-
-
    for(Int_t i = 0; i<numberOfCuts; i++){
-      
+     
       analysisCuts[i] = new AliConversionCuts();
       analysisCuts[i]->InitializeCutsFromCutString(cutarray[i].Data());
       if (doEtaShiftIndCuts) {
          analysisCuts[i]->DoEtaShift(doEtaShiftIndCuts);
          analysisCuts[i]->SetEtaShift(stringShift);
       }
-      
+     
       ConvCutList->Add(analysisCuts[i]);
          
       analysisCuts[i]->SetFillCutHistograms("",kFALSE);
@@ -182,23 +183,23 @@ void AddTask_GammaConvV1_pPb(  Int_t trainConfig = 1,  //change different set of
       analysisMesonCuts[i]->SetFillCutHistograms("");
       analysisCuts[i]->SetAcceptedHeader(HeaderList);
    }
-
    task->SetConversionCutList(numberOfCuts,ConvCutList);
    task->SetMesonCutList(numberOfCuts,MesonCutList);
    task->SetMoveParticleAccordingToVertex(kTRUE);
    task->SetDoMesonAnalysis(kTRUE);
    if (enableQAMesonTask) task->SetDoMesonQA(kTRUE); //Attention new switch for Pi0 QA
-   if (enableQAPhotonTask) task->SetDoPhotonQA(kTRUE);  //Attention new switch small for Photon QA
-
+   if (enableQAMesonTask) task->SetDoPhotonQA(kTRUE);  //Attention new switch small for Photon QA
    //connect containers
    AliAnalysisDataContainer *coutput =
       mgr->CreateContainer(Form("GammaConvV1_%i",trainConfig), TList::Class(),
                            AliAnalysisManager::kOutputContainer,Form("GammaConvV1_%i.root",trainConfig));
-
    mgr->AddTask(task);
    mgr->ConnectInput(task,0,cinput);
    mgr->ConnectOutput(task,1,coutput);
-
    return;
-
 }