start of database macro, adjustment of conservate defaults
authorloizides <cloizides@lbl.gov>
Wed, 22 Jan 2014 21:40:22 +0000 (22:40 +0100)
committerloizides <cloizides@lbl.gov>
Wed, 22 Jan 2014 21:40:22 +0000 (22:40 +0100)
PWG/EMCAL/macros/AddTaskEMCALTender.C
PWG/EMCAL/macros/AddTaskEMCALTenderUsingDatasetDef.C
PWG/EMCAL/macros/ConfigEmcalTenderSupply.C

index 53df1e1..732068b 100644 (file)
@@ -1,7 +1,6 @@
 // $Id$
 
 AliAnalysisTaskSE *AddTaskEMCALTender(
-  const char *perstr    = "LHC11h",//string defining period
   Bool_t distBC         = kTRUE,   //distance to bad channel
   Bool_t recalibClus    = kTRUE,   //recalibrate cluster energy
   Bool_t recalcClusPos  = kTRUE,   //recalculate cluster position
@@ -13,14 +12,14 @@ AliAnalysisTaskSE *AddTaskEMCALTender(
   Bool_t remBC          = kTRUE,   //remove bad channels
   UInt_t nonLinFunct    = AliEMCALRecoUtils::kBeamTestCorrected,
   Bool_t reclusterize   = kTRUE,   //reclusterize
-  Float_t seedthresh    = 0.3,     //seed threshold
-  Float_t cellthresh    = 0.05,    //cell threshold
+  Float_t seedthresh    = 0.100,   //seed threshold
+  Float_t cellthresh    = 0.050,   //cell threshold
   UInt_t clusterizer    = AliEMCALRecParam::kClusterizerv2,
-  Bool_t trackMatch     = kFALSE,  //track matching
+  Bool_t trackMatch     = kTRUE,   //track matching
   Bool_t updateCellOnly = kFALSE,  //only change if you run your own clusterizer task
   Float_t timeMin       = 100e-9,  //minimum time of physical signal in a cell/digit (s)
   Float_t timeMax       = 900e-9,  //maximum time of physical signal in a cell/digit (s)
-  Float_t timeCut       = 50e-9,   //maximum time difference between the digits inside EMC cluster (s)
+  Float_t timeCut       = 900e-9,  //maximum time difference between the digits inside EMC cluster (s)
   const char *pass      = 0        //string defining pass (use none if figured out from path)
 ) 
 {
@@ -40,12 +39,6 @@ AliAnalysisTaskSE *AddTaskEMCALTender(
   AliAnalysisTaskSE *ana = 0;
   AliAnalysisDataContainer *coutput1 = 0;
 
-  TString period(perstr);
-  period.ToLower();
-  if (period == "lhc12a15e")
-    nonLinFunct = AliEMCALRecoUtils::kPi0MCv3;
-  else if (period == "lhc12a15a")
-    nonLinFunct = AliEMCALRecoUtils::kPi0MCv2;
 
   gROOT->LoadMacro("$ALICE_ROOT/PWG/EMCAL/macros/ConfigEmcalTenderSupply.C");
 
@@ -53,11 +46,6 @@ AliAnalysisTaskSE *AddTaskEMCALTender(
                                                              fidRegion, calibEnergy, calibTime, remBC, nonLinFunct, reclusterize, seedthresh, 
                                                              cellthresh, clusterizer, trackMatch, updateCellOnly, timeMin, timeMax, timeCut);
 
-  if (!period.Contains("lhc10h") && !period.Contains("lhc11h")) {
-    ::Info("AddTaskEMCALTender.C","Switching on TrackMatching() since it is not PbPb");
-    EMCALSupply->SwitchOnTrackMatch();
-  }
-
   if (pass) 
     EMCALSupply->SetPass(pass);
 
@@ -68,13 +56,11 @@ AliAnalysisTaskSE *AddTaskEMCALTender(
     alitender->AddSupply(EMCALSupply);
     alitender->SetDefaultCDBStorage("raw://"); 
     ana = alitender;
-
     coutput1 = mgr->CreateContainer("emcal_tender_event", 
                                    AliESDEvent::Class(), 
                                    AliAnalysisManager::kExchangeContainer, 
                                    "default_tender");
-  }
-  else if (evhand->InheritsFrom("AliAODInputHandler")) {
+  } else if (evhand->InheritsFrom("AliAODInputHandler")) {
     AliEmcalTenderTask* emcaltender = mgr->GetTopTasks()->FindObject("AliEmcalTenderTask"); 
     if (!emcaltender)
       emcaltender = new  AliEmcalTenderTask("AliEmcalTenderTask");
@@ -84,8 +70,7 @@ AliAnalysisTaskSE *AddTaskEMCALTender(
                                    AliAODEvent::Class(), 
                                    AliAnalysisManager::kExchangeContainer, 
                                    "default_tender");
-  }
-  else {
+  } else {
     ::Error("AddTaskEMCALTender", "Input event handler not recognized, AOD/ESD expected. Returning...");
     return NULL;
   }
index 21ff735..fb86c9e 100644 (file)
@@ -1,7 +1,7 @@
 // $Id$
 
 AliAnalysisTaskSE *AddTaskEMCALTenderUsingDatasetDef(
-  const char *dstr    = "LHC11h",
+  const char *perstr  = "LHC11h",
   const char *pass    = 0 /*should not be needed*/
 ) 
 {
@@ -14,24 +14,31 @@ AliAnalysisTaskSE *AddTaskEMCALTenderUsingDatasetDef(
   Bool_t remExotic      = kTRUE;   //remove exotic cells
   Bool_t fidRegion      = kFALSE;  //apply fiducial cuts
   Bool_t calibEnergy    = kTRUE;   //calibrate energy
-  Bool_t calibTime      = kTRUE,   //calibrate timing
-  Bool_t remBC          = kTRUE,   //remove bad channels
+  Bool_t calibTime      = kTRUE;   //calibrate timing
+  Bool_t remBC          = kTRUE;   //remove bad channels
   UInt_t nonLinFunct    = AliEMCALRecoUtils::kBeamTestCorrected;
-  Bool_t reclusterize   = kFALSE,  //reclusterize
-  Float_t seedthresh    = 0.3,     //seed threshold
-  Float_t cellthresh    = 0.05,    //cell threshold
-  UInt_t clusterizer    = AliEMCALRecParam::kClusterizerv2,
-  Bool_t trackMatch     = kFALSE,  //track matching
-  Bool_t updateCellOnly = kFALSE,  //only change if you run your own clusterizer task
-  Float_t timeMin       = 100e-9,  //minimum time of physical signal in a cell/digit (s)
-  Float_t timeMax       = 900e-9,  //maximum time of physical signal in a cell/digit (s)
-  Float_t timeCut       = 50e-9,   //maximum time difference between the digits inside EMC cluster (s)
+  Bool_t reclusterize   = kFALSE;  //reclusterize
+  Float_t seedthresh    = 0.100;   //seed threshold
+  Float_t cellthresh    = 0.050;   //cell threshold
+  UInt_t clusterizer    = AliEMCALRecParam::kClusterizerv2;
+  Bool_t trackMatch     = kTRUE;   //track matching
+  Bool_t updateCellOnly = kFALSE;  //only change if you run your own clusterizer task
+  Float_t timeMin       = 100e-9;  //minimum time of physical signal in a cell/digit (s)
+  Float_t timeMax       = 900e-9;  //maximum time of physical signal in a cell/digit (s)
+  Float_t timeCut       = 900e-9;  //maximum time difference between the digits inside EMC cluster (s)
 
+  TString period(perstr);
+  period.ToLower();
+  if (period == "lhc12a15e")
+    nonLinFunct = AliEMCALRecoUtils::kPi0MCv3;
+  else if (period == "lhc12a15a")
+    nonLinFunct = AliEMCALRecoUtils::kPi0MCv2;
 
-AliAnalysisTaskSE *task = AddTaskEMCALTender(
-  "tothink",distBC,recalibClus,recalcClusPos,nonLinearCorr,remExotic,fidRegion,calibEnergy,calibTime,
-  remBC, nonLinFunct, reclusterize, seedthresh, cellthresh, clusterizer, trackMatch, updateCellOnly,
-  timeMin, timeMax, timeCut, 0);
+  AliAnalysisTaskSE *task = AddTaskEMCALTender(
+    distBC, recalibClus, recalcClusPos, nonLinearCorr,
+    remExotic, fidRegion, calibEnergy, calibTime, remBC, nonLinFunct, 
+    reclusterize, seedthresh, cellthresh, clusterizer, trackMatch, updateCellOnly,
+    timeMin, timeMax, timeCut, pass);
 
-  return et;
+  return task;
 }
index da0faa5..03820e4 100644 (file)
@@ -12,15 +12,15 @@ AliEMCALTenderSupply* ConfigEmcalTenderSupply(
   Bool_t remBC          = kTRUE,   //remove bad channels
   UInt_t nonLinFunct    = AliEMCALRecoUtils::kBeamTestCorrected,
   Bool_t reclusterize   = kTRUE,   //reclusterize
-  Float_t seedthresh    = 0.3,     //seed threshold
-  Float_t cellthresh    = 0.05,    //cell threshold
+  Float_t seedthresh    = 0.100,   //seed threshold
+  Float_t cellthresh    = 0.050,   //cell threshold
   UInt_t clusterizer    = AliEMCALRecParam::kClusterizerv2,
-  Bool_t trackMatch     = kFALSE,  //track matching
+  Bool_t trackMatch     = kTRUE,   //track matching
   Bool_t updateCellOnly = kFALSE,  //only change if you run your own clusterizer task
   Float_t timeMin       = 100e-9,  //minimum time of physical signal in a cell/digit (s)
   Float_t timeMax       = 900e-9,  //maximum time of physical signal in a cell/digit (s)
-  Float_t timeCut       = 50e-9    //maximum time difference between the digits inside EMC cluster (s)
-                                             )
+  Float_t timeCut       = 900e-9   //maximum time difference between the digits inside EMC cluster (s)
+)
 {
   AliEMCALTenderSupply *EMCALSupply = new AliEMCALTenderSupply("EMCALtender");  
   EMCALSupply->SetDebugLevel(2);