Macro and default object for the new Pb-Pb entry.

[u/mrichter/AliRoot.git] / EMCAL / AliEMCALRecParam.cxx

diff --git a/EMCAL/AliEMCALRecParam.cxx b/EMCAL/AliEMCALRecParam.cxx
index 20ec7cb8d5f86c7864026fcd9db95f2ed7d142e4..5c165c06cebb4baced25d4d50b842e378606b7ac 100644 (file)
--- a/EMCAL/AliEMCALRecParam.cxx
+++ b/EMCAL/AliEMCALRecParam.cxx
@@ -34,27 +34,27 @@
 
 ClassImp(AliEMCALRecParam)
   
-  TObjArray* AliEMCALRecParam::fgkMaps =0; //ALTRO mappings 
+TObjArray* AliEMCALRecParam::fgkMaps =0; //ALTRO mappings 
 
 AliEMCALRecParam::AliEMCALRecParam() :
   AliDetectorRecoParam(),
   fClusteringThreshold(0.5),
   fW0(4.5),
-  fMinECut(0.45), 
+  fMinECut(0.045), 
   fUnfold(kFALSE), 
   fLocMaxCut(0.03), 
   fTimeCut(1.),// high value, accept all
   fTimeMin(-1.),// small value, accept all
   fTimeMax(1.),// high value, accept all//clustering
-  fTrkCutX(6.0), 
-  fTrkCutY(6.0), 
-  fTrkCutZ(6.0),  
-  fTrkCutR(10.0),
-  fTrkCutAlphaMin(-50.0), 
-  fTrkCutAlphaMax(50.0), 
-  fTrkCutAngle(10000.0),
-  fTrkCutNITS(3.0),
-  fTrkCutNTPC(20.0), //track matching
+  fClusterizerFlag(AliEMCALRecParam::kClusterizerv1),
+  fNRowDiff(1),
+  fNColDiff(1),
+  fMthCutEta(0.025), 
+  fMthCutPhi(0.05),
+  fStep(50),
+  fTrkCutPt(0.0),
+  fTrkCutNITS(0.0),
+  fTrkCutNTPC(50.0), //track matching
   fHighLowGainFactor(16.0), 
   fOrderParameter(2), 
   fTau(2.35), 
@@ -86,7 +86,7 @@ AliEMCALRecParam::AliEMCALRecParam() :
   // Pb Pb
   
   // as a first step, all array elements are initialized to 0.0
-  Int_t i, j;
+  Int_t i=0, j=0;
   for (i = 0; i < 6; i++) {
     for (j = 0; j < 6; j++) {      
       fGamma[i][j] =  fPiZero[i][j] = fHadron[i][j] = 0.; 
@@ -220,7 +220,22 @@ AliEMCALRecParam::AliEMCALRecParam() :
   fHadronEnergyProb[3]=  2.030230e+00;
   fHadronEnergyProb[4]= -6.402242e-02;
   
-  
+  //unfolding  
+  fSSPars[0] = 0.9262;
+  fSSPars[1] = 3.365;
+  fSSPars[2] = 1.548;
+  fSSPars[3] = 0.1625;
+  fSSPars[4] = -0.4195;
+  fSSPars[5] = 0.;
+  fSSPars[6] = 0.;
+  fSSPars[7] = 2.332;
+  fPar5[0] = 12.31;
+  fPar5[1] = -0.007381;
+  fPar5[2] = -0.06936;
+  fPar6[0] = 0.05452; 
+  fPar6[1] = 0.0001228; 
+  fPar6[2] = 0.001361; 
+
 }
 
 //-----------------------------------------------------------------------------
@@ -234,13 +249,13 @@ AliEMCALRecParam::AliEMCALRecParam(const AliEMCALRecParam& rp) :
   fTimeCut(rp.fTimeCut), 
   fTimeMin(rp.fTimeMin),
   fTimeMax(rp.fTimeMax),//clustering
-  fTrkCutX(rp.fTrkCutX), 
-  fTrkCutY(rp.fTrkCutY), 
-  fTrkCutZ(rp.fTrkCutZ),  
-  fTrkCutR(rp.fTrkCutR),
-  fTrkCutAlphaMin(rp.fTrkCutAlphaMin), 
-  fTrkCutAlphaMax(rp.fTrkCutAlphaMax), 
-  fTrkCutAngle(rp.fTrkCutAngle), 
+  fClusterizerFlag(rp.fClusterizerFlag),
+  fNRowDiff(rp.fNRowDiff),
+  fNColDiff(rp.fNColDiff),
+  fMthCutEta(rp.fMthCutEta), 
+  fMthCutPhi(rp.fMthCutPhi),
+  fStep(rp.fStep),
+  fTrkCutPt(rp.fTrkCutPt),
   fTrkCutNITS(rp.fTrkCutNITS),
   fTrkCutNTPC(rp.fTrkCutNTPC), // track matching
   fHighLowGainFactor(rp.fHighLowGainFactor), 
@@ -256,7 +271,7 @@ AliEMCALRecParam::AliEMCALRecParam(const AliEMCALRecParam& rp) :
   //copy constructor
   
   //PID values
-  Int_t i, j;
+  Int_t i=0, j=0;
   for (i = 0; i < 6; i++) {
     for (j = 0; j < 6; j++) {
       fGamma[i][j]       = rp.fGamma[i][j];
@@ -271,6 +286,15 @@ AliEMCALRecParam::AliEMCALRecParam(const AliEMCALRecParam& rp) :
     
   }
   
+  //unfolding  
+  for (i = 0; i < 8; i++) {
+    fSSPars[i] = rp.fSSPars[i];
+  }
+  for (i = 0; i < 3; i++) {
+    fPar5[i] = rp.fPar5[i];
+    fPar6[i] = rp.fPar6[i];
+  }
+
 }
 
 //-----------------------------------------------------------------------------
@@ -287,15 +311,15 @@ AliEMCALRecParam& AliEMCALRecParam::operator = (const AliEMCALRecParam& rp)
     fTimeCut   = rp.fTimeCut;
     fTimeMax   = rp.fTimeMax;
     fTimeMin   = rp.fTimeMin;//clustering
-    fTrkCutX   = rp.fTrkCutX;
-    fTrkCutY   = rp.fTrkCutY;
-    fTrkCutZ   = rp.fTrkCutZ;
-    fTrkCutR   = rp.fTrkCutR;
-    fTrkCutAlphaMin = rp.fTrkCutAlphaMin;
-    fTrkCutAlphaMax = rp.fTrkCutAlphaMax;
-    fTrkCutAngle    = rp.fTrkCutAngle; 
-    fTrkCutNITS     = rp.fTrkCutNITS;
-    fTrkCutNTPC     = rp.fTrkCutNTPC; //track matching
+    fClusterizerFlag   = rp.fClusterizerFlag;
+    fNRowDiff  = rp.fNRowDiff;
+    fNColDiff  = rp.fNColDiff;
+    fMthCutEta         = rp.fMthCutEta;
+    fMthCutPhi         = rp.fMthCutPhi;
+    fStep              = rp.fStep;
+    fTrkCutPt          = rp.fTrkCutPt;
+    fTrkCutNITS        = rp.fTrkCutNITS;
+    fTrkCutNTPC        = rp.fTrkCutNTPC; //track matching
     fHighLowGainFactor = rp.fHighLowGainFactor; 
     fOrderParameter    = rp.fOrderParameter;
     fTau               = rp.fTau;
@@ -307,7 +331,7 @@ AliEMCALRecParam& AliEMCALRecParam::operator = (const AliEMCALRecParam& rp)
     fFitLEDEvents      = rp.fFitLEDEvents;//raw signal
          
     //PID values
-    Int_t i, j;
+    Int_t i=0, j=0;
     for (i = 0; i < 6; i++) {
       for (j = 0; j < 6; j++) {
        fGamma[i][j]       = rp.fGamma[i][j];
@@ -320,7 +344,15 @@ AliEMCALRecParam& AliEMCALRecParam::operator = (const AliEMCALRecParam& rp)
       fPiZeroEnergyProb[i] = rp.fPiZeroEnergyProb[i];
       fHadronEnergyProb[i] = rp.fHadronEnergyProb[i];
     }
-    
+    //unfolding  
+    for (i = 0; i < 8; i++) {
+      fSSPars[i] = rp.fSSPars[i];
+    }
+    for (i = 0; i < 3; i++) {
+      fPar5[i] = rp.fPar5[i];
+      fPar6[i] = rp.fPar6[i];
+    }
+
   }    
   
   return *this;
@@ -383,7 +415,7 @@ AliEMCALRecParam* AliEMCALRecParam::GetLowFluxParam()
   
   //PID parameters for pp  implemented 
   // as a first step, all array elements are initialized to 0.0
-  Int_t i, j;
+  Int_t i=0, j=0;
   for (i = 0; i < 6; i++) {
     for (j = 0; j < 6; j++) {
       params->SetGamma(i,j,0.);
@@ -542,6 +574,7 @@ AliEMCALRecParam* AliEMCALRecParam::GetHighFluxParam()
   params->SetName("High Flux - Pb+Pb");
   params->SetTitle("High Flux - Pb+Pb");
   params->SetEventSpecie(AliRecoParam::kHighMult);
+  params->SetTrkCutPt(0.15);//This value can be higher if necessary
   
   return params;
   
@@ -554,11 +587,26 @@ void AliEMCALRecParam::Print(Option_t * opt) const
   // if nothing is specified print all, if "reco" just clusterization/track matching
   // if "pid", just PID parameters, if "raw", just raw utils parameters.
   if(!strcmp("",opt) || !strcmp("reco",opt)){
+    AliInfo(Form("Clusterizer selected: %d", fClusterizerFlag));
     AliInfo(Form("Clusterization parameters :\n fClusteringThreshold=%.3f,\n fW0=%.3f,\n fMinECut=%.3f,\n fUnfold=%d,\n fLocMaxCut=%.3f,\n fTimeCut=%2.1f ns\n fTimeMin=%2.1f ns\n fTimeMax=%2.1f ns\n",
                 fClusteringThreshold,fW0,fMinECut,fUnfold,fLocMaxCut,fTimeCut*1.e9,fTimeMin*1e9,fTimeMax*1e9));
     
-    AliInfo(Form("Track-matching cuts :\n x %f, y %f, z %f, R %f \n alphaMin %f, alphaMax %f, Angle %f, NITS %f, NTPC %f\n", 
-                                fTrkCutX, fTrkCutY, fTrkCutZ, fTrkCutR,fTrkCutAlphaMin,fTrkCutAlphaMax, fTrkCutAngle,fTrkCutNITS,fTrkCutNTPC));
+    AliInfo(Form("Track-matching cuts :\n dEta<%f, dPhi<%f, step=%f[cm], pT>%f, NITS>%f, NTPC>%f\n", 
+                fMthCutEta, fMthCutPhi, fStep, fTrkCutPt, fTrkCutNITS,fTrkCutNTPC));
+
+    AliInfo(Form("Unfolding parameters, Shower shape function :\n")); 
+    for(Int_t i = 0; i < 8; i++){
+       printf(" %f, ", fSSPars[i]);
+    }
+    printf("\n Parameter 5 : ");
+    for(Int_t i = 0; i < 3; i++){
+      printf(" %f, ", fPar5[i]);
+    }
+    printf("\n Parameter 6 : ");
+    for(Int_t i = 0; i < 3; i++){
+      printf(" %f, ", fPar6[i]);
+    }
+    printf("\n");
   }
   
   if(!strcmp("",opt) || !strcmp("pid",opt)){