fill the fraction of cell energy in cluster in case of unfolding ON
authorgconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 19 Oct 2010 08:55:28 +0000 (08:55 +0000)
committergconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 19 Oct 2010 08:55:28 +0000 (08:55 +0000)
EMCAL/AliEMCALReconstructor.cxx
EMCAL/macros/RecParamDB/AliEMCALSetRecParamCDB.C
EMCAL/macros/RecParamDB/PrintEMCALRecParam.C

index 61c5955..aca000f 100644 (file)
@@ -479,11 +479,11 @@ void AliEMCALReconstructor::FillESD(TTree* digitsTree, TTree* clustersTree,
     for (Int_t iCell=0; iCell<cellMult; iCell++) {
       if (amplFloat[iCell] > 0) {
       absIdList[newCellMult] = (UShort_t)(digitInts[iCell]);
-      //Uncomment when unfolding is done
-      //if(emcCells.GetCellAmplitude(digitInts[iCell])>0)
-      //fracList[newCellMult] = amplFloat[iCell]/(emcCells.GetCellAmplitude(digitInts[iCell])*calibration);//get cell calibration value 
-      //else
-      fracList[newCellMult] = 0; 
+      //Calculate Fraction
+      if(emcCells.GetCellAmplitude(digitInts[iCell])>0 && GetRecParam()->GetUnfold())
+        fracList[newCellMult] = amplFloat[iCell]/(emcCells.GetCellAmplitude(digitInts[iCell]));//get cell calibration value 
+      else
+        fracList[newCellMult] = 0; 
       newCellMult++;
       }
     }
index d1f3b54..e63afff 100644 (file)
@@ -134,7 +134,7 @@ AliEMCALRecParam* GetHighMultiplicityParameters()
   params->SetClusteringThreshold(0.5);
   params->SetW0(4.5);
   params->SetMinECut(0.45);
-  params->SetUnfold(kFALSE);
+  params->SetUnfold(kTRUE);
   params->SetLocMaxCut(0.03);
   params->SetTimeCut(250e-9);//250 ns
   params->SetTimeMin(425e-9);//425 ns
@@ -293,7 +293,7 @@ AliEMCALRecParam* GetHighMultiplicityParameters()
   params->SetNoiseThreshold(3);
   params->SetNPedSamples(5);
   params->SetRemoveBadChannels(kFALSE);
-  params->SetFittingAlgorithm(AliEMCALRawUtils::kStandard);  
+  params->SetFittingAlgorithm(0);//(AliCaloConstants::kStandard);  
   params->SetFALTROUsage(kTRUE); 
   params->SetLEDFit(kFALSE);   
 
@@ -310,7 +310,7 @@ AliEMCALRecParam* GetLowMultiplicityParameters()
   AliEMCALRecParam* params =  AliEMCALRecParam::GetDefaultParameters();
   //params->SetClusterizerFlag(AliEMCALRecParam::kClusterizerNxN);
   params->SetClusterizerFlag(AliEMCALRecParam::kClusterizerv1);
-  
+  params->SetUnfold(kTRUE);
   params->SetClusteringThreshold(0.1); // 100 MeV                                             
   params->SetMinECut(0.01);  //10 MeV          
   params->SetTimeCut(250e-9);//250 ns
@@ -461,7 +461,7 @@ AliEMCALRecParam* GetLowMultiplicityParameters()
   params->SetNoiseThreshold(3);
   params->SetNPedSamples(5);
   params->SetRemoveBadChannels(kFALSE);
-  params->SetFittingAlgorithm(AliEMCALRawUtils::kStandard);  
+  params->SetFittingAlgorithm(0);//(AliCaloConstants::kStandard);  
   params->SetFALTROUsage(kTRUE); 
   params->SetLEDFit(kFALSE);   
 
index fe225fe..8ca3dd0 100644 (file)
@@ -4,7 +4,7 @@
 // Author: Gustavo Conesa (INFN-LNF)
 
 
-void PrintEMCALRecParam(char * file = "$ALICE_ROOT/OCDB/EMCAL/Calib/RecoParam/Run0_999999999_v0_s0.root"){
+void PrintEMCALRecParam(char * file = "$ALICE_ROOT/OCDB/EMCAL/Calib/RecoParam/Run0_999999999_v0_s2.root"){
 
 
 TFile * f = new TFile(file,"READ");
@@ -24,10 +24,10 @@ cout<<"============== "<<rparam->GetName()<<" ==============="<<endl;
 
 cout<<"================================================"<<endl;
 
-//rparam->Print("reco");//Print only clusterizer parameters
+rparam->Print("reco");//Print only clusterizer parameters
 //rparam->Print("pid");//Print only pid parameters
 //rparam->Print("raw");//Print only raw digitization parameters
-rparam->Print("");// Print all
+//rparam->Print("");// Print all
 
 
 }