AliEMCALHadronCorrectionv1 added.

diff --git a/EMCAL/AliEMCALHadronCorrectionv1.cxx b/EMCAL/AliEMCALHadronCorrectionv1.cxx
new file mode 100755 (executable)
index 0000000..e81b55d
--- /dev/null
+++ b/EMCAL/AliEMCALHadronCorrectionv1.cxx
@@ -0,0 +1,107 @@
+/**************************************************************************
+ * Copyright(c) 1998-2002, ALICE Experiment at CERN, All rights reserved. *
+ *                                                                        *
+ * Author: The ALICE Off-line Project.                                    *
+ * Contributors are mentioned in the code where appropriate.              *
+ *                                                                        *
+ * Permission to use, copy, modify and distribute this software and its   *
+ * documentation strictly for non-commercial purposes is hereby granted   *
+ * without fee, provided that the above copyright notice appears in all   *
+ * copies and that both the copyright notice and this permission notice   *
+ * appear in the supporting documentation. The authors make no claims     *
+ * about the suitability of this software for any purpose. It is          *
+ * provided "as is" without express or implied warranty.                  *
+ **************************************************************************/
+
+
+// Author : Mark Horner (LBL/UCT)
+
+#include "AliEMCALHadronCorrectionv1.h"
+#include "AliEMCALGeometry.h"
+static Double_t par_look_up[HCPARAMETERS][HCPARAMETERSETS] = 
+{  
+    {-1.68945e-04 , -1.68945e-04},   
+    { 2.09684e-02 ,  2.09684e-02},   
+    {-1.45683e-04 , -1.45683e-04},   
+    { 6.64803e-03 ,  6.64803e-03},   
+    { 2.00834e-02 ,  2.00834e-02},   
+    { 9.22074e-07 ,  9.22074e-07}
+};
+
+
+
+ClassImp(AliEMCALHadronCorrectionv1)
+
+AliEMCALHadronCorrectionv1* AliEMCALHadronCorrectionv1::fHadrCorr = 0;
+
+void AliEMCALHadronCorrectionv1::SetGeometry(AliEMCALGeometry *geometry)
+{
+    if (!geometry)
+    {
+       SetParameters();
+    }else
+    {
+        SetParameters(geometry->GetName());
+    }
+    return;    
+}      
+       
+void AliEMCALHadronCorrectionv1::SetGeometry(TString name)
+{
+  if ( name == ""              ||
+       name == "EMCAL_5655_21" ||
+       name == "EMCALArch1a"   ||
+       name == "EMCALArch1aN"  ||
+       name == "G56_2_55_19"   ||
+       name == "G56_2_55_19_104_14" )
+  { // set parameters to this hadron correction
+     for (Int_t i=0;i<6;i++)
+     {
+          fPar[i] = par_look_up[i][0];  
+     }
+  }else if( name == "EMCAL_6564_21" ||  
+           name == "G65_2_64_19" )
+  {      
+     for (Int_t i=0;i<6;i++)
+     {
+          fPar[i] = par_look_up[i][1];  
+     }
+  }else
+  {
+    printf("Geometry not defined in hadron correction\n"); 
+  }      
+       
+}      
+
+       
+AliEMCALHadronCorrectionv1::AliEMCALHadronCorrectionv1(const char *name,const char *title) 
+                           :AliEMCALHadronCorrection(name, title)
+{
+  fHadrCorr = this;
+}
+
+/*
+AliEMCALHadronCorrectionv1::AliEMCALHadronCorrectionv1(const char *name,const char *title,AliEMCALGeometry *geometry)
+{
+
+  fHadrCorr = this;
+  SetGeometry(geometry);  
+       
+}      
+*/
+
+AliEMCALHadronCorrectionv1*
+AliEMCALHadronCorrectionv1::Instance()
+{
+  if (! fHadrCorr) new AliEMCALHadronCorrectionv1();
+  return fHadrCorr;
+}
+
+Double_t 
+AliEMCALHadronCorrectionv1::GetEnergy(const Double_t pmom,const Double_t eta,const Int_t gid)
+{
+
+  Double_t value =  fPar[5]*pmom*pmom*pmom+ fPar[0]*pmom*pmom+fPar[1]*pmom +fPar[2]*pmom*eta +fPar[3]*eta + fPar[4];
+  return value;
+   
+}

diff --git a/EMCAL/AliEMCALHadronCorrectionv1.h b/EMCAL/AliEMCALHadronCorrectionv1.h
new file mode 100755 (executable)
index 0000000..e2d4549
--- /dev/null
+++ b/EMCAL/AliEMCALHadronCorrectionv1.h
@@ -0,0 +1,39 @@
+#ifndef ALIEMCALHADRONCORRECTIONV1_H
+#define ALIEMCALHADRONCORRECTIONV1_H
+/* Copyright(c) 1998-2002, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice     */
+
+//                  
+//*-- Author: Mark Horner (LBL/UCT)
+//
+#include "AliEMCALHadronCorrection.h"
+
+
+#define HCPARAMETERS    6
+#define HCPARAMETERSETS 2
+
+class AliEMCALGeometry;
+
+class AliEMCALHadronCorrectionv1: public AliEMCALHadronCorrection {
+ public:
+    static  AliEMCALHadronCorrectionv1* Instance();
+    virtual Double_t GetEnergy(const Double_t pmom,const Double_t eta,const Int_t gid); 
+    Double_t GetEnergy(const Double_t pmom, const Double_t eta) 
+       {return GetEnergy(pmom,eta,7);}
+    
+    void SetGeometry(TString name); 
+    void SetGeometry(AliEMCALGeometry *geometry); 
+    virtual ~AliEMCALHadronCorrectionv1() {}
+ protected:
+    AliEMCALHadronCorrectionv1(const char *name="HadronCorrectionv1", const char *title="Hadron Correction");
+    
+//    AliEMCALHadronCorrectionv1(const char *name="HadronCorrectionv1", const char *title="Hadron Correction",AliEMCALGeometry *geometry = NULL);
+ private:
+    void SetParameters(TString name = "");
+    static AliEMCALHadronCorrectionv1* fHadrCorr;
+    Double_t fPar[6];
+    
+    ClassDef(AliEMCALHadronCorrectionv1,1) // Hadron correction for EMC (version for MDC)
+};
+
+#endif // ALIEMCALHADRONCORRECTIONV1_H

diff --git a/EMCAL/EMCALLinkDef.h b/EMCAL/EMCALLinkDef.h
index 0e8a0dcc3cf06686623b48a4307b46b696263836..74a725fc1920da509a3da2b42cd0b4fb4c652113 100644 (file)
--- a/EMCAL/EMCALLinkDef.h
+++ b/EMCAL/EMCALLinkDef.h
@@ -23,7 +23,8 @@
 #pragma link C++ class AliEMCALJetMicroDst+;
 #pragma link C++ class AliEMCALFast+;
 #pragma link C++ class AliEMCALHadronCorrection+;               
-#pragma link C++ class AliEMCALHadronCorrectionv0+;             
+#pragma link C++ class AliEMCALHadronCorrectionv0+;
+#pragma link C++ class AliEMCALHadronCorrectionv1+;             
 #pragma link C++ class AliEMCALParton+;                
 #pragma link C++ class AliEMCALReconstructioner+;              
 #pragma link C++ class AliEMCALTrackSegment+;          

diff --git a/EMCAL/libEMCAL.pkg b/EMCAL/libEMCAL.pkg
index f0a126493ee07e99c82edb792a86642f9e2feef0..2bb6b54f7756d50882c5c95b10816a08e0394d81 100644 (file)
--- a/EMCAL/libEMCAL.pkg
+++ b/EMCAL/libEMCAL.pkg
@@ -20,6 +20,7 @@ AliEMCALClusterizer.cxx \
 AliEMCALClusterizerv1.cxx \
 AliEMCALHadronCorrection.cxx \
 AliEMCALHadronCorrectionv0.cxx \
+AliEMCALHadronCorrectionv1.cxx \
 AliEMCALJetMicroDst.cxx \
 AliEMCALParton.cxx \
 AliEMCALReconstructioner.cxx \