Updates by M. Horner and J. Klay
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 17 Mar 2003 14:26:06 +0000 (14:26 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 17 Mar 2003 14:26:06 +0000 (14:26 +0000)
EMCAL/AliEMCALHadronCorrectionv1.cxx
EMCAL/AliEMCALParton.cxx
EMCAL/AliEMCALParton.h

index e81b55d..5abd87f 100755 (executable)
 
 #include "AliEMCALHadronCorrectionv1.h"
 #include "AliEMCALGeometry.h"
+#include "iostream.h"
+#include "TMath.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}
+    {-2.82271e-4 , -2.39954e-4},   
+    {2.50796e-2  , 2.07172e-2},   
+    {1.02861e-3  , 1.48576e-3},   
+    {2.11539e-2  , -1.38473-2},   
+    {2.27003e-2  , 2.78252e-2},   
+    {1.65078e-6  , 1.51821e-6}
 };
 
 
@@ -39,15 +41,21 @@ void AliEMCALHadronCorrectionv1::SetGeometry(AliEMCALGeometry *geometry)
     if (!geometry)
     {
        SetParameters();
+       fSamplingFraction = geometry->GetSampling();
+       cout<<"Setting the sampling fraction to :"<<fSamplingFraction<<endl; 
     }else
     {
         SetParameters(geometry->GetName());
+       fSamplingFraction = geometry->GetSampling();
+       cout<<"Setting the sampling fraction to :"<<fSamplingFraction<<endl; 
     }
     return;    
 }      
        
-void AliEMCALHadronCorrectionv1::SetGeometry(TString name)
+void AliEMCALHadronCorrectionv1::SetGeometry(TString name,Double_t fs = 1.0)
 {
+  cout << "Setting sampling fraction to "<<fSamplingFraction<<endl;    
+   fSamplingFraction = fs;     
   if ( name == ""              ||
        name == "EMCAL_5655_21" ||
        name == "EMCALArch1a"   ||
@@ -55,16 +63,20 @@ void AliEMCALHadronCorrectionv1::SetGeometry(TString name)
        name == "G56_2_55_19"   ||
        name == "G56_2_55_19_104_14" )
   { // set parameters to this hadron correction
+     cout<<"HC parameters!"<<endl;
      for (Int_t i=0;i<6;i++)
      {
           fPar[i] = par_look_up[i][0];  
+          cout <<fPar[i]<<endl;
      }
   }else if( name == "EMCAL_6564_21" ||  
            name == "G65_2_64_19" )
   {      
+      cout<<"HC parameters!"<<endl;
      for (Int_t i=0;i<6;i++)
      {
           fPar[i] = par_look_up[i][1];  
+          cout <<fPar[i]<<endl;
      }
   }else
   {
@@ -100,8 +112,8 @@ AliEMCALHadronCorrectionv1::Instance()
 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;
+  Double_t etai = TMath::Abs(eta); 
+  Double_t value =  fPar[5]*pmom*pmom*pmom+ fPar[0]*pmom*pmom+fPar[1]*pmom +fPar[2]*pmom*etai +fPar[3]*etai + fPar[4];
+  return fSamplingFraction*value;
    
 }
index d542c6e..5cf41cc 100644 (file)
@@ -16,6 +16,9 @@
 /* $Id:  */
 /*
   $Log$
+  Revision 1.3  2003/02/05 17:24:25  morsch
+  Modifiactions by M. Horner.
+
   Revision 1.2  2002/10/14 14:55:35  hristov
   Merging the VirtualMC branch to the main development branch (HEAD)
 
 
 #include "AliEMCALParton.h"
 #include "Ecommon.h"
+#include <iostream.h>
 ClassImp(AliEMCALParton)   
     
 //____________________________________________________________________________
 AliEMCALParton::AliEMCALParton()
 {
   // Default constructor
-  fTrackEnergy = 0;
-  fTrackEta    = 0;
-  fTrackPhi    = 0;
-  fTrackPDG    = 0;
+//  fTrackEnergy = 0;
+//  fTrackEta    = 0;
+//  fTrackPhi    = 0;
+//  fTrackPDG    = 0;
   fNTracks     = 0;
 }
 
@@ -50,36 +54,24 @@ AliEMCALParton::AliEMCALParton(Float_t energy, Float_t phi, Float_t eta)
   fEnergy = energy;
   fPhi    = phi;
   fEta    = eta;
-  fTrackEnergy = 0;
-  fTrackEta    = 0;
-  fTrackPhi    = 0;
-  fTrackPDG    = 0;
+//  fTrackEnergy = 0;
+//  fTrackEta    = 0;
+//  fTrackPhi    = 0;
+//  fTrackPDG    = 0;
   fNTracks     = 0;
 }
 
 void AliEMCALParton::SetTrackList(Int_t NTracks, Float_t* Energy,  Float_t* Eta, Float_t* Phi, Int_t* PDG)
 {
 
-  if (fNTracks)
-  {
-    delete[] fTrackEnergy;
-    delete[] fTrackEta;
-    delete[] fTrackPhi;
-    delete[] fTrackPDG;
-  }
   fNTracks     = NTracks;
-  fTrackEnergy = new Float_t[NTracks];
-  fTrackEta    = new Float_t[NTracks];
-  fTrackPhi    = new Float_t[NTracks];
-  fTrackPDG    = new Int_t[NTracks];
-  
   for (Int_t i=0;i<NTracks;i++)
   {
     fTrackEnergy[i] = Energy[i] ;
     fTrackEta[i]    = Eta[i];
     fTrackPhi[i]    = Phi[i];
     fTrackPDG[i]    = PDG[i];
-  } 
+  }
 }
 
 void AliEMCALParton::GetTrackList(Float_t* Energy,  Float_t* Eta, Float_t* Phi, Int_t* PDG)
@@ -99,9 +91,5 @@ void AliEMCALParton::GetTrackList(Float_t* Energy,  Float_t* Eta, Float_t* Phi,
 AliEMCALParton::~AliEMCALParton()
 {
   // Destructor
-  delete[] fTrackEnergy;
-  delete[] fTrackEta;
-  delete[] fTrackPhi;
-  delete[] fTrackPDG;
   
 }
index 5345874..9d91221 100644 (file)
@@ -26,14 +26,14 @@ public:
   
 protected:
   Float_t   fEnergy;   // Jet Energy
-  Float_t*  fTrackEnergy;  // Jet Tracks Energy
-  Float_t*  fTrackEta;     // Jet Tracks Eta
-  Float_t*  fTrackPhi;     // Jet Tracks Phi
-  Int_t*    fTrackPDG;       // Jet Tracks PDG code
-  Int_t     fNTracks;
   Float_t   fEta;      // Jet Phi
   Float_t   fPhi;      // Jet Eta
-  ClassDef(AliEMCALParton,3) // Jet for EMCAL
+  Int_t     fNTracks;      // Number of tracks 
+  Float_t  fTrackEnergy[1000];  // Jet Tracks Energy
+  Float_t  fTrackEta[1000];     // Jet Tracks Eta
+  Float_t  fTrackPhi[1000];     // Jet Tracks Phi
+  Int_t    fTrackPDG[1000];     // Jet Tracks PDG code
+  ClassDef(AliEMCALParton,4) // Jet for EMCAL
     
 } ;