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++;
}
}
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
params->SetNoiseThreshold(3);
params->SetNPedSamples(5);
params->SetRemoveBadChannels(kFALSE);
- params->SetFittingAlgorithm(AliEMCALRawUtils::kStandard);
+ params->SetFittingAlgorithm(0);//(AliCaloConstants::kStandard);
params->SetFALTROUsage(kTRUE);
params->SetLEDFit(kFALSE);
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
params->SetNoiseThreshold(3);
params->SetNPedSamples(5);
params->SetRemoveBadChannels(kFALSE);
- params->SetFittingAlgorithm(AliEMCALRawUtils::kStandard);
+ params->SetFittingAlgorithm(0);//(AliCaloConstants::kStandard);
params->SetFALTROUsage(kTRUE);
params->SetLEDFit(kFALSE);
// 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");
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
}