updated for Non-HFE systemaitc study
authorssakai <ssakai@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 30 May 2013 13:28:49 +0000 (13:28 +0000)
committerssakai <ssakai@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 30 May 2013 13:28:49 +0000 (13:28 +0000)
PWGHF/hfe/AliAnalysisTaskHFECal.cxx
PWGHF/hfe/AliAnalysisTaskHFECal.h
PWGHF/hfe/macros/AddTaskHFECal.C
PWGHF/hfe/macros/configs/PbPb/ConfigHFECal.C

index 111c172..940e5e0 100644 (file)
@@ -107,6 +107,7 @@ AliAnalysisTaskHFECal::AliAnalysisTaskHFECal(const char *name)
   ,fOpeningAngleCut(0.1)
   ,fInvmassCut(0)       // no mass
   ,fSetMassConstraint(kTRUE)
+  ,fSetMassWidthCut(kFALSE)
   ,fNoEvents(0)
   ,fEMCAccE(0)
   ,hEMCAccE(0)
@@ -242,6 +243,7 @@ AliAnalysisTaskHFECal::AliAnalysisTaskHFECal()
   ,fOpeningAngleCut(0.1)
   ,fInvmassCut(0)       // no mass
   ,fSetMassConstraint(kTRUE)
+  ,fSetMassWidthCut(kFALSE)
   ,fNoEvents(0)
   ,fEMCAccE(0)
   ,hEMCAccE(0)
@@ -1530,6 +1532,8 @@ void AliAnalysisTaskHFECal::SelectPhotonicElectron(Int_t itrack, Double_t cent,
       }
     recg.GetMass(mass,width);
 
+    if(fSetMassWidthCut && width>1e10)continue;
+
         // angle   
     openingAngle = ge1.GetAngle(ge2);
     if(fFlagLS) fOpeningAngleLS->Fill(openingAngle);
index 7a8d5aa..d20a572 100644 (file)
@@ -53,6 +53,7 @@ class AliAnalysisTaskHFECal : public AliAnalysisTaskSE {
   void SetOpeningAngleCut (Double_t openingAngle) {fOpeningAngleCut = openingAngle;};
   void SetInvariantMassCut (Double_t invmass) {fInvmassCut = invmass;};
   void SetMassConstraint       (Bool_t MassConstraint)         { fSetMassConstraint    = MassConstraint; };
+  void SetMassWidthCut         (Bool_t MassWidthCut)           { fSetMassWidthCut      = MassWidthCut; };
   void SetQAHist (int qahist) {fqahist = qahist;};
   AliHFEpid *GetPID() const { return fPID; }
   void SetRejectKinkMother(Bool_t rejectKinkMother = kFALSE) { fRejectKinkMother = rejectKinkMother; };
@@ -89,6 +90,7 @@ class AliAnalysisTaskHFECal : public AliAnalysisTaskSE {
   Double_t             fOpeningAngleCut;       //openingAngle cut value
   Double_t             fInvmassCut;            //invariant mass cut value
   Bool_t                fSetMassConstraint;            // Set mass constraint
+  Bool_t                fSetMassWidthCut;              // Set mass constraint
  
   int ftriggers[48][60];//!
   int ftriggersCut[48][60];//!
index 7183fad..9d9ab6d 100644 (file)
@@ -1,4 +1,4 @@
-AliAnalysisTask *AddTaskHFECal(Bool_t MassConst)
+AliAnalysisTask *AddTaskHFECal(Bool_t MassConst, Bool_t MassWidthCut)
 {
   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
   if (!mgr) {
@@ -27,9 +27,9 @@ AliAnalysisTask *AddTaskHFECal(Bool_t MassConst)
   //analysis task 
   gROOT->LoadMacro("$ALICE_ROOT/PWGHF/hfe/macros/configs/PbPb/ConfigHFECal.C");
 
-  AliAnalysisTaskHFECal *hfetaskCent = ConfigHFECal(MCthere,MassConst);
-  AliAnalysisTaskHFECal *hfetaskTrig= ConfigHFECal(MCthere,MassConst);
-  AliAnalysisTaskHFECal *hfetaskTrig2= ConfigHFECal(MCthere,MassConst);
+  AliAnalysisTaskHFECal *hfetaskCent = ConfigHFECal(MCthere,MassConst,MassWidthCut);
+  AliAnalysisTaskHFECal *hfetaskTrig= ConfigHFECal(MCthere,MassConst,MassWidthCut);
+  AliAnalysisTaskHFECal *hfetaskTrig2= ConfigHFECal(MCthere,MassConst,MassWidthCut);
  
   mgr->AddTask(hfetaskCent);
   mgr->AddTask(hfetaskTrig);
@@ -74,7 +74,7 @@ AliAnalysisTask *AddTaskHFECal(Bool_t MassConst)
   //if(MCthere)
     //{
     //MB trigger
-    AliAnalysisTaskHFECal *hfetaskMB = ConfigHFECal(MCthere,MassConst);
+    AliAnalysisTaskHFECal *hfetaskMB = ConfigHFECal(MCthere,MassConst,MassWidthCut);
     mgr->AddTask(hfetaskMB);
     hfetaskMB->SelectCollisionCandidates(AliVEvent::kMB);
 
index 7ce0f3f..b5d5a22 100644 (file)
@@ -1,4 +1,4 @@
-AliAnalysisTaskHFECal* ConfigHFECal(Bool_t useMC,Bool_t MassConst){\r
+AliAnalysisTaskHFECal* ConfigHFECal(Bool_t useMC,Bool_t MassConst,Bool_t MassWidthCut){\r
   //\r
   // HFE standard task configuration\r
   //\r
@@ -22,6 +22,7 @@ AliAnalysisTaskHFECal* ConfigHFECal(Bool_t useMC,Bool_t MassConst){
   printf("task ------------------------ %p\n ", task);\r
   task->SetHFECuts(hfecuts);\r
   task->SetMassConstraint(MassConst);\r
+  task->SetMassWidthCut(MassWidthCut);\r
   Double_t masscut = 0.05;\r
   if(!MassConst)masscut = 0.1;\r
   task->SetInvariantMassCut(masscut);\r