update some analysis parameters
authorgconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 1 Mar 2011 16:15:50 +0000 (16:15 +0000)
committergconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 1 Mar 2011 16:15:50 +0000 (16:15 +0000)
PWG4/macros/AddTaskPartCorr.C
PWG4/macros/AddTaskPartCorrM.C
PWG4/macros/AddTaskPi0.C

index 1ca8bab..61ba646 100644 (file)
@@ -32,13 +32,15 @@ AliAnalysisTaskParticleCorrelation *AddTaskPartCorr(TString inputDataType, TStri
   //if(!kSimulation) reader->SetFiredTriggerClassName("CINT1B-ABCE-NOPF-ALL");
   if(calorimeter == "EMCAL") {
     reader->SwitchOnEMCALCells();  
-    reader->SwitchOnEMCAL();
+    reader->SwitchOnEMCAL(); 
   }
   if(calorimeter == "PHOS") { 
     reader->SwitchOnPHOSCells();  
     reader->SwitchOnPHOS();
   }
   
+   reader->SwitchOnSuspiciousClustersRemoval();  //EMCAL
+
   // for case inputDataType="deltaAOD", no need to fill the EMCAL/PHOS cluster lists
   if(inputDataType.Contains("delta")){
     reader->SwitchOffEMCAL();
@@ -121,16 +123,16 @@ AliAnalysisTaskParticleCorrelation *AddTaskPartCorr(TString inputDataType, TStri
   anaphoton->SetMultiplicity(80, 120);
 
   if(calorimeter == "PHOS"){
-    anaphoton->SetNCellCut(0);// At least 2 cells
-    anaphoton->SetMinPt(0.);
+    anaphoton->SetNCellCut(2);// At least 3 cells
+    anaphoton->SetMinPt(0.3);
     anaphoton->SetMinDistanceToBadChannel(2, 4, 5);
   }
   else {//EMCAL
-    //anaphoton->SetNCellCut(0);// At least 2 cells
-    anaphoton->SetMinPt(0.1); // no effect minium EMCAL cut.
-    if(!kUseKinematics) anaphoton->SetTimeCut(400,900);// Time window of [400-900] ns
-    anaphoton->SetMinDistanceToBadChannel(6, 12, 18);
-    //anaphoton->SetMinDistanceToBadChannel(1, 2, 3);//For new releases.
+    anaphoton->SetNCellCut(1);// At least 2 cells
+    anaphoton->SetMinPt(0.3); 
+    //if(!kUseKinematics) anaphoton->SetTimeCut(400,900);// Time window of [400-900] ns
+    //anaphoton->SetMinDistanceToBadChannel(6, 12, 18);
+    anaphoton->SetMinDistanceToBadChannel(1, 2, 3);//For new releases.
   }
   anaphoton->SetCalorimeter(calorimeter);
   if(kUseKinematics) anaphoton->SwitchOnDataMC() ;//Access MC stack and fill more histograms
index a17a868..99b7428 100644 (file)
@@ -34,6 +34,8 @@ AliAnalysisTaskParticleCorrelationM *AddTaskPartCorrM(TString data, TString calo
     reader->SwitchOnPHOS();
   }
   
+   reader->SwitchOnSuspiciousClustersRemoval();  //EMCAL
+
   // for case data="deltaAOD", no need to fill the EMCAL/PHOS cluster lists
   if(data.Contains("delta")){
     reader->SwitchOffEMCAL();
@@ -148,10 +150,12 @@ AliAnalysisTaskParticleCorrelationM *AddTaskPartCorrM(TString data, TString calo
     
     if(calorimeter == "EMCAL"){
       anaphoton->SetNCellCut(1);
-      if(!kUseKinematics) anaphoton->SetTimeCut(525, 725);
-      anaphoton->SetMinDistanceToBadChannel(4, 5, 10);
+      anaphoton->SetMinPt(0.3); 
+      //if(!kUseKinematics) anaphoton->SetTimeCut(525, 725);
+      anaphoton->SetMinDistanceToBadChannel(1, 2, 3);
     }
     else{
+      anaphoton->SetMinPt(0.3);       
       anaphoton->SetMinDistanceToBadChannel(2, 4, 5);
       anaphoton->SetNCellCut(2);
     }
index eb65001..b82236b 100644 (file)
@@ -51,6 +51,8 @@ AliAnalysisTaskParticleCorrelation *AddTaskPi0(TString data, TString calorimeter
     reader->SwitchOnPHOSCells();  
     reader->SwitchOnPHOS();
   }
+
+   reader->SwitchOnSuspiciousClustersRemoval();  //EMCAL
   
   // for case data="deltaAOD", no need to fill the EMCAL/PHOS cluster lists
   if(data.Contains("delta")){
@@ -128,16 +130,16 @@ AliAnalysisTaskParticleCorrelation *AddTaskPi0(TString data, TString calorimeter
   AliAnaPhoton *anaphoton = new AliAnaPhoton();
   anaphoton->SetDebug(-1); //10 for lots of messages
   if(calorimeter == "PHOS"){
-    anaphoton->SetNCellCut(0);// At least 2 cells
-    anaphoton->SetMinPt(0.);
+    anaphoton->SetNCellCut(2);// At least 2 cells
+    anaphoton->SetMinPt(3.);
     anaphoton->SetMinDistanceToBadChannel(2, 4, 5);
   }
   else {//EMCAL
-    //anaphoton->SetNCellCut(0);// At least 2 cells
-    anaphoton->SetMinPt(0.1); // no effect minium EMCAL cut.
-    if(!kUseKinematics) anaphoton->SetTimeCut(400,900);// Time window of [400-900] ns
-    anaphoton->SetMinDistanceToBadChannel(6, 12, 18);//For officially produced ESDs/AODs
-    //anaphoton->SetMinDistanceToBadChannel(1, 2, 3);//For filtered AODs, new releases.
+    anaphoton->SetNCellCut(1);// At least 2 cells
+    anaphoton->SetMinPt(0.3); // no effect minium EMCAL cut.
+    //if(!kUseKinematics) anaphoton->SetTimeCut(400,900);// Time window of [400-900] ns
+    //anaphoton->SetMinDistanceToBadChannel(6, 12, 18);//For officially produced ESDs/AODs
+    anaphoton->SetMinDistanceToBadChannel(1, 2, 3);//For filtered AODs, new releases.
   }
   anaphoton->SetCalorimeter(calorimeter);
   if(kUseKinematics) anaphoton->SwitchOnDataMC() ;//Access MC stack and fill more histograms