Changes to run in central train
authordainese <dainese@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 15 Dec 2010 17:40:09 +0000 (17:40 +0000)
committerdainese <dainese@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 15 Dec 2010 17:40:09 +0000 (17:40 +0000)
PWG3/vertexingHF/AddD2HTrain.C
PWG3/vertexingHF/AliRDHFCutsDplustoKpipi.cxx
PWG3/vertexingHF/macros/AddTaskDStarSpectra.C
PWG3/vertexingHF/macros/AddTaskDplus.C

index 4016132..ef1f8c0 100644 (file)
@@ -2,14 +2,14 @@ Int_t AddD2HTrain(Bool_t readMC=kTRUE,
                  Int_t addD0Mass=1,
                  Int_t addD0MassLS=1,
                  Int_t addDplus=1,
-                 Int_t addLSD0=1,
-                 Int_t addLSJpsi=1,
-                 Int_t addCFD0=1,
+                 Int_t addLSD0=0,
+                 Int_t addLSJpsi=0,
+                 Int_t addCFD0=0,
                  Int_t addPromptD0=1,
-                 Int_t addDs=1,
+                 Int_t addDs=0,
                  Int_t addDStar=1,
-                 Int_t addDStarJets=1,
-                 Int_t addCFDStar=1) {
+                 Int_t addDStarJets=0,
+                 Int_t addCFDStar=0) {
   // 
   // Tasks of the D2H subgroup of PWG3 that can run in the Official Train
   //
index 9fcdeac..6b50a3d 100644 (file)
@@ -538,7 +538,6 @@ void AliRDHFCutsDplustoKpipi::SetStandardCutsPP2010() {
   anacutsval[6][10]=0.06;
   anacutsval[6][11]=0.06;
   anacutsval[6][12]=0.06;
-  anacutsval[6][13]=0.06;
   
   
   anacutsval[7][0]=0.08;
@@ -554,7 +553,6 @@ void AliRDHFCutsDplustoKpipi::SetStandardCutsPP2010() {
   anacutsval[7][10]=0.02;
   anacutsval[7][11]=0.02;
   anacutsval[7][12]=0.02;
-  anacutsval[7][13]=0.02;
   
   
   anacutsval[8][0]=0.5;
@@ -570,7 +568,6 @@ void AliRDHFCutsDplustoKpipi::SetStandardCutsPP2010() {
   anacutsval[8][10]=0.2;
   anacutsval[8][11]=0.2;
   anacutsval[8][12]=0.2;
-  anacutsval[8][13]=0.2;
     
     
   anacutsval[9][0]=0.95;
@@ -586,7 +583,6 @@ void AliRDHFCutsDplustoKpipi::SetStandardCutsPP2010() {
   anacutsval[9][10]=0.90;
   anacutsval[9][11]=0.90;
   anacutsval[9][12]=0.90;
-  anacutsval[9][13]=0.90;
   
   anacutsval[10][0]=0.0055;
   anacutsval[10][1]=0.0055;
@@ -601,9 +597,8 @@ void AliRDHFCutsDplustoKpipi::SetStandardCutsPP2010() {
   anacutsval[10][10]=0.0;
   anacutsval[10][11]=0.0;
   anacutsval[10][12]=0.0;
-  anacutsval[10][13]=0.0;
-  
   
+  SetGlobalIndex(nvars,nptbins);
   SetPtBins(nptbins+1,ptbins);
   SetCuts(nvars,nptbins,anacutsval);
   SetUsePID(kTRUE);
index a8ca832..f5efc53 100644 (file)
@@ -1,10 +1,10 @@
 //if like define a different number of signal for TPC PID\r
 //by default the task is anyway computing 1, 2 and 3 sigmas\r
-const Int_t    numberOfSigmasPID  = 3;\r
+const Int_t    numberOfSigmasPID  = 2;\r
 // option to switch on and off the TPC PID.\r
 const Bool_t usePID = kTRUE;\r
 // analysis type... TO BE REMOVED!!!\r
-const Bool_t anaType   = 1;//0 HD; 1 UU;\r
+const Bool_t anaType   = 0;//0 HD; 1 UU;\r
 //----------------------------------------------------\r
 \r
 AliAnalysisTaskSEDStarSpectra *AddTaskDStarSpectra(Bool_t theMCon=kFALSE)\r
@@ -18,14 +18,16 @@ AliAnalysisTaskSEDStarSpectra *AddTaskDStarSpectra(Bool_t theMCon=kFALSE)
   \r
   // cuts are stored in a TFile generated by makeTFile4CutsDStartoKpipi.C in ./macros/\r
   // set there the cuts!!!!!\r
+  Bool_t stdcuts=kFALSE;\r
   TFile* filecuts=new TFile("DStartoKpipiCuts.root");\r
   if(!filecuts->IsOpen()){\r
     cout<<"Input file not found: exit"<<endl;\r
-    return;\r
+    stdcuts=kTRUE;\r
   }\r
 \r
   AliRDHFCutsDStartoKpipi* RDHFDStartoKpipi=new AliRDHFCutsDStartoKpipi();\r
-  RDHFDStartoKpipi = (AliRDHFCutsDStartoKpipi*)filecuts->Get("DStartoKpipiCuts");\r
+  if(stdcuts) RDHFDStartoKpipi->SetStandardCutsPP2010();\r
+  else RDHFDStartoKpipi = (AliRDHFCutsDStartoKpipi*)filecuts->Get("DStartoKpipiCuts");\r
   RDHFDStartoKpipi->SetName("DStartoKpipiCuts");\r
 \r
   // mm let's see if everything is ok\r
index 1782fa0..ede268a 100644 (file)
@@ -1,6 +1,6 @@
 AliAnalysisTaskSEDplus *AddTaskDplus(Bool_t storeNtuple=kFALSE,
                                     Bool_t readMC=kFALSE,
-                                    TString filename="./DplustoKpipiCuts.root")
+                                    TString filename="DplustoKpipiCuts.root")
 {
   //                                                                                                                                    
   // Test macro for the AliAnalysisTaskSE for D+ candidates 
@@ -13,24 +13,26 @@ AliAnalysisTaskSEDplus *AddTaskDplus(Bool_t storeNtuple=kFALSE,
   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
   if (!mgr) {
     ::Error("AddTaskDplus", "No analysis manager to connect to.");
-    return NULL;
   }
 
+  Bool_t stdcuts=kFALSE;
   TFile* filecuts=new TFile(filename.Data());
   if(!filecuts->IsOpen()){
-    cout<<"Input file not found: exit"<<endl;
-    return;
+    cout<<"Input file not found: using standard cuts"<<endl;
+    stdcuts=kTRUE;
   }
   
   
   //Analysis Task
 
   
-   AliRDHFCutsDplustoKpipi* analysiscuts=new AliRDHFCutsDplustoKpipi();
-  analysiscuts = (AliRDHFCutsDplustoKpipi*)filecuts->Get("AnalysisCuts");
+  AliRDHFCutsDplustoKpipi* analysiscuts=new AliRDHFCutsDplustoKpipi();
+  if(stdcuts) analysiscuts->SetStandardCutsPP2010();
+  else analysiscuts = (AliRDHFCutsDplustoKpipi*)filecuts->Get("AnalysisCuts");
 
   AliRDHFCutsDplustoKpipi* prodcuts=new AliRDHFCutsDplustoKpipi();
-  prodcuts = (AliRDHFCutsDplustoKpipi*)filecuts->Get("ProdCuts");
+  if(stdcuts) prodcuts->SetStandardCutsPP2010();
+  else prodcuts = (AliRDHFCutsDplustoKpipi*)filecuts->Get("AnalysisCuts");
   
   //AliRDHFCutsDplustoKpipi *prodcuts = (AliRDHFCutsDplustoKpipi*)fileCuts->Get("ProdCuts");
   //AliRDHFCutsDplustoKpipi *analysiscuts = (AliRDHFCutsDplustoKpipi*)fileCuts->Get("AnalysisCuts");