]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - EMCAL/AliEMCALRecParam.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff
? search:
summary | shortlog | log | commit | commitdiff | tree
raw | inline | side by side
increase array to cope with pbpb
[u/mrichter/AliRoot.git] / EMCAL / AliEMCALRecParam.cxx
diff --git a/EMCAL/AliEMCALRecParam.cxx b/EMCAL/AliEMCALRecParam.cxx
index a6ae3ef2878e9c5e0807bd08e85248df7d71701b..675911c4c2d99dd5e15a997188c68d630b04952a 100644 (file)
--- a/EMCAL/AliEMCALRecParam.cxx
+++ b/EMCAL/AliEMCALRecParam.cxx
@@ -34,7 +34,7 @@
 
 ClassImp(AliEMCALRecParam)
   
-  TObjArray* AliEMCALRecParam::fgkMaps =0; //ALTRO mappings 
+TObjArray* AliEMCALRecParam::fgkMaps =0; //ALTRO mappings 
 
 AliEMCALRecParam::AliEMCALRecParam() :
   AliDetectorRecoParam(),
@@ -46,6 +46,7 @@ AliEMCALRecParam::AliEMCALRecParam() :
   fTimeCut(1.),// high value, accept all
   fTimeMin(-1.),// small value, accept all
   fTimeMax(1.),// high value, accept all//clustering
+  fClusterizerFlag(AliEMCALRecParam::kClusterizerv1),
   fTrkCutX(6.0), 
   fTrkCutY(6.0), 
   fTrkCutZ(6.0),  
@@ -62,7 +63,8 @@ AliEMCALRecParam::AliEMCALRecParam() :
   fNPedSamples(5), 
   fRemoveBadChannels(kFALSE),
   fFittingAlgorithm(0), 
-  fUseFALTRO(kTRUE)//raw signal
+  fUseFALTRO(kTRUE), 
+  fFitLEDEvents(kFALSE)//raw signal
 {
   // default reco values
   
@@ -85,7 +87,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.; 
@@ -150,15 +152,12 @@ AliEMCALRecParam::AliEMCALRecParam() :
   fHadron[1][3] = -2.368500e-03; 
   fHadron[1][4] = 0.;
   
-  
   fHadron[2][0] = 0.;
   fHadron[2][1] =  -2.463152e-02;
   fHadron[2][2] = 1.349257e-01;
   fHadron[2][3] = -1.089440e-03;
   fHadron[2][4] = 0.;
   
-  
-  
   fHadron[3][0] = 0.;
   fHadron[3][1] = 5.101560e-01;
   fHadron[3][2] = 1.458679e-01;
@@ -222,10 +221,23 @@ 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; 
 
+}
 
 //-----------------------------------------------------------------------------
 AliEMCALRecParam::AliEMCALRecParam(const AliEMCALRecParam& rp) :
@@ -238,6 +250,7 @@ AliEMCALRecParam::AliEMCALRecParam(const AliEMCALRecParam& rp) :
   fTimeCut(rp.fTimeCut), 
   fTimeMin(rp.fTimeMin),
   fTimeMax(rp.fTimeMax),//clustering
+  fClusterizerFlag(rp.fClusterizerFlag),
   fTrkCutX(rp.fTrkCutX), 
   fTrkCutY(rp.fTrkCutY), 
   fTrkCutZ(rp.fTrkCutZ),  
@@ -254,12 +267,13 @@ AliEMCALRecParam::AliEMCALRecParam(const AliEMCALRecParam& rp) :
   fNPedSamples(rp.fNPedSamples),       
   fRemoveBadChannels(rp.fRemoveBadChannels),
   fFittingAlgorithm(rp.fFittingAlgorithm),  
-  fUseFALTRO(fUseFALTRO) //raw signal
+  fUseFALTRO(rp.fUseFALTRO),
+  fFitLEDEvents(rp.fFitLEDEvents) //raw signal
 {
   //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];
@@ -274,6 +288,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,9 +310,10 @@ AliEMCALRecParam& AliEMCALRecParam::operator = (const AliEMCALRecParam& rp)
     fMinECut   = rp.fMinECut;
     fUnfold    = rp.fUnfold;
     fLocMaxCut = rp.fLocMaxCut; 
-       fTimeCut   = rp.fTimeCut;
-       fTimeMax   = rp.fTimeMax;
-       fTimeMin   = rp.fTimeMin;//clustering
+    fTimeCut   = rp.fTimeCut;
+    fTimeMax   = rp.fTimeMax;
+    fTimeMin   = rp.fTimeMin;//clustering
+    fClusterizerFlag = rp.fClusterizerFlag;
     fTrkCutX   = rp.fTrkCutX;
     fTrkCutY   = rp.fTrkCutY;
     fTrkCutZ   = rp.fTrkCutZ;
@@ -306,10 +330,11 @@ AliEMCALRecParam& AliEMCALRecParam::operator = (const AliEMCALRecParam& rp)
     fNPedSamples       = rp.fNPedSamples; 
     fRemoveBadChannels = rp.fRemoveBadChannels;
     fFittingAlgorithm  = rp.fFittingAlgorithm;
-    fUseFALTRO         = rp.fUseFALTRO;//raw signal
+    fUseFALTRO         = rp.fUseFALTRO;
+    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];
@@ -322,7 +347,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;
@@ -385,7 +418,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.);
@@ -556,10 +589,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("Clusterization parameters :\n fClusteringThreshold=%.3f,\n fW0=%.3f,\n fMinECut=%.3f,\n fUnfold=%d,\n fLocMaxCut=%.3f,\n fTimeCut=%f \n fTimeMin=%f \n fTimeMax=%f \n",
-                fClusteringThreshold,fW0,fMinECut,fUnfold,fLocMaxCut,fTimeCut,fTimeMin,fTimeMax));
+    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 %\n", fTrkCutX, fTrkCutY, fTrkCutZ, fTrkCutR,fTrkCutAlphaMin,fTrkCutAlphaMax, fTrkCutAngle,fTrkCutNITS,fTrkCutNTPC));
+    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("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)){
@@ -597,8 +646,8 @@ void AliEMCALRecParam::Print(Option_t * opt) const
   if(!strcmp("",opt) || !strcmp("raw",opt)){
     AliInfo(Form("Raw signal parameters: \n gain factor=%f, order=%d, tau=%f, noise threshold=%d, nped samples=%d \n",
                 fHighLowGainFactor,fOrderParameter,fTau,fNoiseThreshold,fNPedSamples));
-    AliInfo(Form("Raw signal: with bad channels? %d, \n \t with fitting algorithm %d, \n \t Use FALTRO %d \n",
-                fRemoveBadChannels, fFittingAlgorithm, fUseFALTRO));
+    AliInfo(Form("Raw signal: remove bad channels? %d, \n \t with fitting algorithm %d, \n \t Use FALTRO %d, Fit LED events %d \n",
+                fRemoveBadChannels, fFittingAlgorithm, fUseFALTRO, fFitLEDEvents));
   }
 }