#include "AliZDCReconstructor.h"
#include "AliZDCPedestals.h"
#include "AliZDCEnCalib.h"
+#include "AliZDCSaturationCalib.h"
#include "AliZDCTowerCalib.h"
#include "AliZDCMBCalib.h"
#include "AliZDCTDCCalib.h"
AliZDCReconstructor:: AliZDCReconstructor() :
fPedData(GetPedestalData()),
fEnCalibData(GetEnergyCalibData()),
+ fSatCalibData(GetSaturationCalibData()),
fTowCalibData(GetTowerCalibData()),
fTDCCalibData(GetTDCCalibData()),
fRecoMode(0),
// if(fgRecoParam) delete fgRecoParam;
if(fPedData) delete fPedData;
if(fEnCalibData) delete fEnCalibData;
+ if(fSatCalibData) delete fSatCalibData;
if(fTowCalibData) delete fTowCalibData;
if(fgMBCalibData) delete fgMBCalibData;
if(fESDZDC) delete fESDZDC;
||((beamType.CompareTo("PP"))==0) || ((beamType.CompareTo("P-P"))==0)){
fRecoMode=1;
}
+ else if(((beamType.CompareTo("p-A"))==0) || ((beamType.CompareTo("A-p"))==0)
+ ||((beamType.CompareTo("P-A"))==0) || ((beamType.CompareTo("A-P"))==0)){
+ fRecoMode=1;
+ }
else if((beamType.CompareTo("A-A")) == 0 || (beamType.CompareTo("AA")) == 0){
fRecoMode=2;
if(!fgRecoParam) fgRecoParam = const_cast<AliZDCRecoParam*>(GetRecoParam());
fBeamEnergy = 1380.;
}
+ AliInfo(Form("\n ZDC reconstruction mode %d (1 -> p-p, 2-> A-A)\n\n",fRecoMode));
+
fESDZDC = new AliESDZDC();
}
||((beamType.CompareTo("PP"))==0) || ((beamType.CompareTo("P-P"))==0)){
fRecoMode=1;
}
+ else if(((beamType.CompareTo("p-A"))==0) || ((beamType.CompareTo("A-p"))==0)
+ ||((beamType.CompareTo("P-A"))==0) || ((beamType.CompareTo("A-P"))==0)){
+ fRecoMode=1;
+ }
else if((beamType.CompareTo("A-A")) == 0 || (beamType.CompareTo("AA")) == 0){
fRecoMode=2;
if(!fgRecoParam) fgRecoParam = const_cast<AliZDCRecoParam*>(GetRecoParam());
// **** Pb-Pb data taking 2010 -> subtracting some ch. from correlation ****
// Not interested in o.o.t. signals (ADC modules 2, 3)
//if(adcMod == 2 || adcMod == 3) continue;
- if(((det==1 && quad==0) || (det==3))){
- if(det == 1){
- if(adcMod==0 || adcMod==1){
- if(gain==0) adcZN1[quad] = rawData.GetADCValue();
- else adcZN1lg[quad] = rawData.GetADCValue();
- }
- else if(adcMod==2 || adcMod==3){
- if(gain==0) adcZN1oot[quad] = rawData.GetADCValue();
- else adcZN1ootlg[quad] = rawData.GetADCValue();
- }
+ // **** Pb-Pb data taking 2011 -> subtracting only ZEM from correlation ****
+ if(det==3){
+ if(adcMod==0 || adcMod==1){
+ if(gain==0) adcZEM[quad-1] = rawData.GetADCValue();
+ else adcZEMlg[quad-1] = rawData.GetADCValue();
}
- else if(det == 3){
- if(adcMod==0 || adcMod==1){
- if(gain==0) adcZEM[quad-1] = rawData.GetADCValue();
- else adcZEMlg[quad-1] = rawData.GetADCValue();
- }
- else if(adcMod==2 || adcMod==3){
- if(gain==0) adcZEMoot[quad-1] = rawData.GetADCValue();
- else adcZEMootlg[quad-1] = rawData.GetADCValue();
- }
+ else if(adcMod==2 || adcMod==3){
+ if(gain==0) adcZEMoot[quad-1] = rawData.GetADCValue();
+ else adcZEMootlg[quad-1] = rawData.GetADCValue();
}
}
// When oot values are read the ADC modules 2, 3 can be skipped!!!
// *************************************************************************
if(quad != 5){ // ZDCs (not reference PTMs)
- if(det==1 && quad!=0){
+ if(det==1){
pedindex = quad;
if(gain == 0) tZN1Corr[quad] += (Float_t) (rawData.GetADCValue()-meanPed[pedindex]);
else tZN1Corr[quad+5] += (Float_t) (rawData.GetADCValue()-meanPed[pedindex+kNch]);
if(itdc==iprevtdc) ihittdc++;
else ihittdc=0;
iprevtdc=itdc;
- tdcData[itdc][ihittdc] = rawData.GetZDCTDCDatum();
+ if(ihittdc<4) tdcData[itdc][ihittdc] = rawData.GetZDCTDCDatum();
// Ch. debug
//if(ihittdc==0) printf(" TDC%d %d ",itdc, tdcData[itdc][ihittdc]);
}// ZDC TDC DATA
sPMRef2[1] = pmReflg[0] - (corrCoeff1[23+kNch]*pmRefootlg[1]+corrCoeff0[23+kNch]);
}
if(fPedSubMode==0 && fRecoMode==2){
- // **** Pb-Pb data taking 2010 -> subtracting some ch. from correlation ****
- tZN1Corr[0] = adcZN1[0] - (corrCoeff1[0]*adcZN1oot[0]+corrCoeff0[0]);
- tZN1Corr[5] = adcZN1lg[0] - (corrCoeff1[kNch]*adcZN1ootlg[0]+corrCoeff0[kNch]);
+ // **** Pb-Pb data taking 2011 -> subtracting some ch. from correlation ****
+ //tZN1Corr[0] = adcZN1[0] - (corrCoeff1[0]*adcZN1oot[0]+corrCoeff0[0]);
+ //tZN1Corr[5] = adcZN1lg[0] - (corrCoeff1[kNch]*adcZN1ootlg[0]+corrCoeff0[kNch]);
// Ch. debug
//printf(" adcZN1 %d adcZN1oot %d tZN1Corr %1.2f \n", adcZN1[0],adcZN1oot[0],tZN1Corr[0]);
//printf(" adcZN1lg %d adcZN1ootlg %d tZN1Corrlg %1.2f \n", adcZN1lg[0],adcZN1ootlg[0],tZN1Corr[5]);
}
else if(fPedSubMode==0 && fRecoMode==1){
// **** p-p data taking 2011 -> temporary patch to overcome DA problem ****
- tZN1Corr[0] = adcZN1[0] - meanPed[0];
- tZN1Corr[5] = adcZN1lg[0] - meanPed[kNch];
//
dZEM1Corr[0] = adcZEM[0] - meanPed[10];
dZEM1Corr[1] = adcZEMlg[0] - meanPed[10+kNch];
equalCoeffZP2[ji] = fTowCalibData->GetZP2EqualCoeff(ji);
}
// --- Energy calibration factors ------------------------------------
- Float_t calibEne[6];
+ Float_t calibEne[6], calibSatZNA[4], calibSatZNC[4];
// **** Energy calibration coefficient set to 1
// **** (no trivial way to calibrate in p-p runs)
for(Int_t ij=0; ij<6; ij++) calibEne[ij] = fEnCalibData->GetEnCalib(ij);
+ for(Int_t ij=0; ij<4; ij++){
+ calibSatZNA[ij] = fSatCalibData->GetZNASatCalib(ij);
+ calibSatZNC[ij] = fSatCalibData->GetZNCSatCalib(ij);
+ }
// ****** Equalization of detector responses
Float_t equalTowZN1[10], equalTowZN2[10], equalTowZP1[10], equalTowZP2[10];
}
}
// Ch. debug
- /*printf("\n ------------- EQUALIZATION -------------\n");
+ printf("\n ------------- EQUALIZATION -------------\n");
printf(" ADCZN1 [%1.2f %1.2f %1.2f %1.2f %1.2f]\n",
equalTowZN1[0],equalTowZN1[1],equalTowZN1[2],equalTowZN1[3],equalTowZN1[4]);
printf(" ADCZP1 [%1.2f %1.2f %1.2f %1.2f %1.2f]\n",
equalTowZN2[0],equalTowZN2[1],equalTowZN2[2],equalTowZN2[3],equalTowZN2[4]);
printf(" ADCZP2 [%1.2f %1.2f %1.2f %1.2f %1.2f]\n",
equalTowZP2[0],equalTowZP2[1],equalTowZP2[2],equalTowZP2[3],equalTowZP2[4]);
- printf(" ----------------------------------------\n");*/
+ printf(" ----------------------------------------\n");
+
+ // *** p-A RUN 2013 -> new calibration object
+ // to take into account saturation in ZN PMC
+ // -> 5th order pol. fun. to be applied BEFORE en. calibration
+ equalTowZN1[0] = equalTowZN1[0] + calibSatZNC[0]*equalTowZN1[0]*equalTowZN1[0] +
+ calibSatZNC[1]*equalTowZN1[0]*equalTowZN1[0]*equalTowZN1[0] +
+ calibSatZNC[2]*equalTowZN1[0]*equalTowZN1[0]*equalTowZN1[0]*equalTowZN1[0] +
+ calibSatZNC[3]*equalTowZN1[0]*equalTowZN1[0]*equalTowZN1[0]*equalTowZN1[0]*equalTowZN1[0];
+ equalTowZN2[0] = equalTowZN2[0] + calibSatZNA[0]*equalTowZN2[0]*equalTowZN2[0] +
+ calibSatZNA[1]*equalTowZN2[0]*equalTowZN2[0]*equalTowZN2[0] +
+ calibSatZNA[2]*equalTowZN2[0]*equalTowZN2[0]*equalTowZN2[0]*equalTowZN2[0] +
+ calibSatZNA[3]*equalTowZN2[0]*equalTowZN2[0]*equalTowZN2[0]*equalTowZN2[0]*equalTowZN2[0];
+
+ // Ch. debug
+ printf("\n ------------- SATURATION CORRECTION -------------\n");
+ printf(" ZNC PMC %1.2f\n", equalTowZN1[0]);
+ printf(" ZNA PMC %1.2f\n", equalTowZN2[0]);
+ printf(" ----------------------------------------\n");
// ****** Summed response for hadronic calorimeter (SUMMED and then CALIBRATED!)
Float_t calibSumZN1[]={0,0}, calibSumZN2[]={0,0}, calibSumZP1[]={0,0}, calibSumZP2[]={0,0};
equalCoeffZP2[ji] = fTowCalibData->GetZP2EqualCoeff(ji);
}
// --- Energy calibration factors ------------------------------------
- Float_t calibEne[6];
- // The energy calibration object already takes into account of E_beam
- // -> the value from the OCDB can be directly used (Jul 2010)
+ Float_t calibEne[6], calibSatZNA[4], calibSatZNC[4];
+ // **** Energy calibration coefficient set to 1
+ // **** (no trivial way to calibrate in p-p runs)
for(Int_t ij=0; ij<6; ij++) calibEne[ij] = fEnCalibData->GetEnCalib(ij);
+ for(Int_t ij=0; ij<4; ij++){
+ calibSatZNA[ij] = fSatCalibData->GetZNASatCalib(ij);
+ calibSatZNC[ij] = fSatCalibData->GetZNCSatCalib(ij);
+ }
// ****** Equalization of detector responses
Float_t equalTowZN1[10], equalTowZN2[10], equalTowZP1[10], equalTowZP2[10];
printf(" ----------------------------------------\n");
*/
+ // *** p-A RUN 2013 -> new calibration object
+ // to take into account saturation in ZN PMC
+ // -> 5th order pol. fun. to be applied BEFORE en. calibration
+ equalTowZN1[0] = equalTowZN1[0] + calibSatZNC[0]*equalTowZN1[0]*equalTowZN1[0] +
+ calibSatZNC[1]*equalTowZN1[0]*equalTowZN1[0]*equalTowZN1[0] +
+ calibSatZNC[2]*equalTowZN1[0]*equalTowZN1[0]*equalTowZN1[0]*equalTowZN1[0] +
+ calibSatZNC[3]*equalTowZN1[0]*equalTowZN1[0]*equalTowZN1[0]*equalTowZN1[0]*equalTowZN1[0];
+ equalTowZN2[0] = equalTowZN2[0] + calibSatZNA[0]*equalTowZN2[0]*equalTowZN2[0] +
+ calibSatZNA[1]*equalTowZN2[0]*equalTowZN2[0]*equalTowZN2[0] +
+ calibSatZNA[2]*equalTowZN2[0]*equalTowZN2[0]*equalTowZN2[0]*equalTowZN2[0] +
+ calibSatZNA[3]*equalTowZN2[0]*equalTowZN2[0]*equalTowZN2[0]*equalTowZN2[0]*equalTowZN2[0];
+
// ****** Summed response for hadronic calorimeter (SUMMED and then CALIBRATED!)
Float_t calibSumZN1[]={0,0}, calibSumZN2[]={0,0}, calibSumZP1[]={0,0}, calibSumZP2[]={0,0};
for(Int_t gi=0; gi<5; gi++){
calibSumZP2[1] = calibSumZP2[1]*calibEne[3];
//
Float_t sumZEM[]={0,0}, calibZEM1[]={0,0}, calibZEM2[]={0,0};
- calibZEM1[0] = corrADCZEM1[0]*calibEne[4]*8.;
+ calibZEM1[0] = corrADCZEM1[0]*calibEne[4];
calibZEM1[1] = corrADCZEM1[1]*calibEne[4];
- calibZEM2[0] = corrADCZEM2[0]*calibEne[5]*8.;
+ calibZEM2[0] = corrADCZEM2[0]*calibEne[5];
calibZEM2[1] = corrADCZEM2[1]*calibEne[5];
for(Int_t k=0; k<2; k++) sumZEM[k] = calibZEM1[k] + calibZEM2[k];
Double_t b=0., bA=0., bC=0.;
if(fIsCalibrationMB == kFALSE){
- // ****** Reconstruction parameters ------------------
- if(!fgRecoParam) fgRecoParam = const_cast<AliZDCRecoParam*>(GetRecoParam());
- if(!fgRecoParam){
- AliError(" RecoParam object not retrieved correctly: not reconstructing event!!!");
- return;
- }
- TH1D* hNpartDist = fgRecoParam->GethNpartDist();
- TH1D* hbDist = fgRecoParam->GethbDist();
- Float_t fClkCenter = fgRecoParam->GetClkCenter();
- if(!hNpartDist || !hbDist){
- AliError("Something wrong in Glauber MC histos got from AliZDCREcoParamPbPb: NO EVENT RECO FOR ZDC DATA!!!\n\n");
- return;
- }
-
+ // ****** Reconstruction parameters ------------------
+ if(!fgRecoParam) fgRecoParam = const_cast<AliZDCRecoParam*>(GetRecoParam());
+ if(!fgRecoParam){
+ AliError(" RecoParam object not retrieved correctly: not reconstructing ZDC event!!!");
+ return;
+ }
+ TH1D* hNpartDist = fgRecoParam->GethNpartDist();
+ TH1D* hbDist = fgRecoParam->GethbDist();
+ Float_t fClkCenter = fgRecoParam->GetClkCenter();
+ if(!hNpartDist || !hbDist){
+ AliError("Something wrong in Glauber MC histos got from AliZDCREcoParamPbPb: NO EVENT RECO FOR ZDC DATA!!!\n\n");
+ //return;
+ }
+ else{
if(!fgMBCalibData) fgMBCalibData = const_cast<AliZDCMBCalib*>(GetMBCalibData());
TH2F *hZDCvsZEM = fgMBCalibData->GethZDCvsZEM();
TH2F *hZDCCvsZEM = fgMBCalibData->GethZDCCvsZEM();
delete line;
delete lineC; delete lineA;
-
+ }
} // ONLY IF fIsCalibrationMB==kFALSE
Bool_t energyFlag = kTRUE;
fESDZDC->SetZDCScaler(counts);
}
- Int_t tdcValues[32][4];
- Float_t tdcCorrected[32][4];
+ Int_t tdcValues[32][4] = {{0,}};
+ Float_t tdcCorrected[32][4] = {{0.,}};
for(Int_t jk=0; jk<32; jk++){
for(Int_t lk=0; lk<4; lk++){
tdcValues[jk][lk] = reco.GetZDCTDCData(jk, lk);
// we try to keep the TDC oscillations as low as possible!
for(Int_t jk=0; jk<32; jk++){
for(Int_t lk=0; lk<4; lk++){
- if(TMath::Abs(tdcValues[jk][lk])>1e-10){
+ if(tdcValues[jk][lk]!=0.){
tdcCorrected[jk][lk] = 0.025*(tdcValues[jk][lk]-tdcValues[15][0])+fMeanPhase;
// Sep 2011: TDC ch. from 8 to 13 centered around 0 using OCDB
if(jk>=8 && jk<=13) tdcCorrected[jk][lk] = tdcCorrected[jk][lk] - tdcOffset[jk-8];
//Ch. debug
- //if((jk>=8 && jk<=13) || jk==15) printf(" *** tdcOffset%d %f tdcCorr%d %f ",jk,tdcOffset[jk-8],tdcCorrected[jk][lk]);
+ //if(jk>=8 && jk<=13) printf(" *** tdcOffset%d %f tdcCorr%d %f \n",jk,tdcOffset[jk-8],tdcCorrected[jk][lk]);
}
}
return calibdata;
}
+//_____________________________________________________________________________
+AliZDCSaturationCalib* AliZDCReconstructor::GetSaturationCalibData() const
+{
+
+ // Getting energy and equalization calibration object for ZDC set
+
+ AliCDBEntry *entry = AliCDBManager::Instance()->Get("ZDC/Calib/SaturationCalib");
+ if(!entry) AliFatal("No calibration data loaded!");
+ entry->SetOwner(kFALSE);
+
+ AliZDCSaturationCalib *calibdata = dynamic_cast<AliZDCSaturationCalib*> (entry->GetObject());
+ if(!calibdata) AliFatal("Wrong calibration object in calibration file!");
+
+ return calibdata;
+}
+
//_____________________________________________________________________________
AliZDCTowerCalib* AliZDCReconstructor::GetTowerCalibData() const
{