]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - ZDC/AliZDCCalibData.cxx
Changing fabs into TMath::Abs
[u/mrichter/AliRoot.git] / ZDC / AliZDCCalibData.cxx
index 301b32f634efb97ee336b2ddce3df599f1e24899..f0dc6b457ac2be38d81000577c607c623d3d4788 100644 (file)
 ClassImp(AliZDCCalibData)
 
 //________________________________________________________________
-AliZDCCalibData::AliZDCCalibData()
+AliZDCCalibData::AliZDCCalibData():
+TNamed()
 {
-//  fHistMeanPed=0;
   Reset();
 }
 
 //________________________________________________________________
-AliZDCCalibData::AliZDCCalibData(const char* name)
+AliZDCCalibData::AliZDCCalibData(const char* name):
+TNamed(),
+fZEMEndValue(0),
+fZEMCutFraction(0),
+fDZEMSup(0),
+fDZEMInf(0),
+fEZN1MaxValue(0),
+fEZP1MaxValue(0),
+fEZDC1MaxValue(0),
+fEZN2MaxValue(0),
+fEZP2MaxValue(0),
+fEZDC2MaxValue(0)
 {
   // Constructor
   TString namst = "Calib_";
   namst += name;
   SetName(namst.Data());
   SetTitle(namst.Data());
-//  fHistMeanPed=0;
   Reset();
 }
 
@@ -48,15 +58,13 @@ AliZDCCalibData::AliZDCCalibData(const char* name)
 AliZDCCalibData::AliZDCCalibData(const AliZDCCalibData& calibda) :
   TNamed(calibda)
 {
-// copy constructor
+  // Copy constructor
   SetName(calibda.GetName());
   SetTitle(calibda.GetName());
   Reset();
-  for(int t=0; t<47; t++){
+  for(int t=0; t<48; t++){
      fMeanPedestal[t] = calibda.GetMeanPed(t);
      fMeanPedWidth[t] = calibda.GetMeanPedWidth(t);
-  }
-  for(int t=0; t<44; t++){
      fOOTPedestal[t]  = calibda.GetOOTPed(t);
      fOOTPedWidth[t]  = calibda.GetOOTPedWidth(t);
      fPedCorrCoeff[0][t] = calibda.GetPedCorrCoeff0(t);
@@ -64,8 +72,10 @@ AliZDCCalibData::AliZDCCalibData(const AliZDCCalibData& calibda) :
   }
   for(int t=0; t<6; t++)  fEnCalibration[t] = calibda.GetEnCalib(t);
   //
-  fEZEMEndValue    = calibda.GetEZEMEndValue();   
-  fEZEMCutFraction = calibda.GetEZEMCutFraction();
+  fZEMEndValue    = calibda.GetZEMEndValue();   
+  fZEMCutFraction = calibda.GetZEMCutFraction();
+  fDZEMSup       = calibda.GetDZEMSup();
+  fDZEMInf       = calibda.GetDZEMInf();
 }
 
 //________________________________________________________________
@@ -75,19 +85,17 @@ AliZDCCalibData &AliZDCCalibData::operator =(const AliZDCCalibData& calibda)
   SetName(calibda.GetName());
   SetTitle(calibda.GetName());
   Reset();
-  for(int t=0; t<47; t++){
+  for(int t=0; t<48; t++){
      fMeanPedestal[t] = calibda.GetMeanPed(t);
      fMeanPedWidth[t] = calibda.GetMeanPedWidth(t);
-  }
-  for(int t=0; t<44; t++){
      fOOTPedestal[t]  = calibda.GetOOTPed(t);
      fOOTPedWidth[t]  = calibda.GetOOTPedWidth(t);
      fPedCorrCoeff[0][t] = calibda.GetPedCorrCoeff0(t);
      fPedCorrCoeff[1][t] = calibda.GetPedCorrCoeff1(t);
   }
   for(int t=0; t<6; t++) fEnCalibration[t] = calibda.GetEnCalib(t);
-  fEZEMEndValue    = calibda.GetEZEMEndValue();
-  fEZEMCutFraction = calibda.GetEZEMCutFraction();
+  fZEMEndValue    = calibda.GetZEMEndValue();
+  fZEMCutFraction = calibda.GetZEMCutFraction();
 
   return *this;
 }
@@ -95,70 +103,37 @@ AliZDCCalibData &AliZDCCalibData::operator =(const AliZDCCalibData& calibda)
 //________________________________________________________________
 AliZDCCalibData::~AliZDCCalibData()
 {
-//  CleanHistos();
 }
 
 //________________________________________________________________
 void AliZDCCalibData::Reset()
 {
   // Reset
-  memset(fMeanPedestal,0,47*sizeof(Float_t));
-  memset(fMeanPedWidth,0,47*sizeof(Float_t));
-  memset(fOOTPedestal,0,44*sizeof(Float_t));
-  memset(fOOTPedWidth,0,44*sizeof(Float_t));
+  memset(fMeanPedestal,0,48*sizeof(Float_t));
+  memset(fMeanPedWidth,0,48*sizeof(Float_t));
+  memset(fOOTPedestal,0,48*sizeof(Float_t));
+  memset(fOOTPedWidth,0,48*sizeof(Float_t));
   memset(fEnCalibration,0,6*sizeof(Float_t));
+  memset(fZN1EqualCoeff,0,5*sizeof(Float_t));
+  memset(fZP1EqualCoeff,0,5*sizeof(Float_t));
+  memset(fZN2EqualCoeff,0,5*sizeof(Float_t));
+  memset(fZP2EqualCoeff,0,5*sizeof(Float_t));
 }                                                                                       
 
-/*
-//________________________________________________________________
-void AliZDCCalibData::CleanHistos()
-{
-  if (fHistMeanPed) delete fHistMeanPed; fHistMeanPed = 0;
-}
-
-//________________________________________________________________
-void AliZDCCalibData::PrepHistos()
-{
-//  CleanHistos(); // this gives a segm.viol!
-  Int_t   kNChannels = 47;
-  Float_t kMaxPedVal = 47.;
-  TString hname = GetName();  hname += "_Pedestals";
-  fHistMeanPed = new TH1F(hname.Data(),hname.Data(),kNChannels,0.,kMaxPedVal);
-  for(int i=0; i<47; i++)  fHistMeanPed->SetBinContent(i+1,GetMeanPed(i));
-}
-*/
 
 //________________________________________________________________
 void  AliZDCCalibData::Print(Option_t *) const
 {
    // Printing of calibration object
-   printf("\n #######  In-time pedestal values (mean value, sigma)     ####### \n");
-   for(int t=0; t<44; t++){
-     if(t==0 || t==24) printf("\n-------- ZN1 HighRes -------- \n");
-     else if(t==5 || t==29) printf("\n-------- ZN1 LowRes -------- \n");
-     else if(t==10 || t==34) printf("\n-------- ZP1 HighRes -------- \n");
-     else if(t==15 || t==39) printf("\n-------- ZP1 LowRes -------- \n");
-     else if(t==20) printf("\n-------- ZEM1 HighRes --------  \n");
-     else if(t==21) printf("\n-------- ZEM1 LowRes --------  \n");
-     else if(t==22) printf("\n-------- ZEM2 HighRes --------  \n");
-     else if(t==23) printf("\n-------- ZEM2 LowRes --------  \n");
-     printf("ADC%d (%.1f, %.1f)  ",t,fMeanPedestal[t],fMeanPedWidth[t]);
-   }
+   printf("\n ####### In-time pedestal values (mean value, sigma) ####### \n");
+   for(int t=0; t<48; t++) 
+      printf("\t ADC%d (%.1f, %.1f)\n",t,fMeanPedestal[t],fMeanPedWidth[t]);
    //
-   printf("\n\n\n #######      Out-of-time pedestal values (mean value, sigma) ####### \n");
-   for(int t=0; t<44; t++){
-     if(t==0 || t==24) printf("\n-------- ZN1 HighRes -------- \n");
-     else if(t==5 || t==29) printf("\n-------- ZN1 LowRes -------- \n");
-     else if(t==10 || t==34) printf("\n-------- ZP1 HighRes -------- \n");
-     else if(t==15 || t==39) printf("\n-------- ZP1 LowRes -------- \n");
-     else if(t==20) printf("\n-------- ZEM1 HighRes --------  \n");
-     else if(t==21) printf("\n-------- ZEM1 LowRes --------  \n");
-     else if(t==22) printf("\n-------- ZEM2 HighRes --------  \n");
-     else if(t==23) printf("\n-------- ZEM2 LowRes --------  \n");
-     printf("ADC%d (%.1f, %.1f)  ",t,fOOTPedestal[t],fOOTPedWidth[t]);
-   }
+   printf("\n\n ####### Out-of-time pedestal values (mean value, sigma) ####### \n");
+   for(int t=0; t<48; t++)
+      printf("\t ADC-OoT%d (%.1f, %.1f)\n",t,fOOTPedestal[t],fOOTPedWidth[t]);
  
-   printf("\n\n\n #######      Energy calibration coefficients ####### \n");
+   printf("\n\n ####### Energy calibration coefficients #######        \n");
    printf("  ZN1 = %.4f (E[TeV]/ADCch.) \n",fEnCalibration[0]);
    printf("  ZP1 = %.4f (E[TeV]/ADCch.) \n",fEnCalibration[1]);
    printf("  ZN2 = %.4f (E[TeV]/ADCch.) \n",fEnCalibration[2]);
@@ -166,42 +141,54 @@ void  AliZDCCalibData::Print(Option_t *) const
    printf("  ZEM1 = %.2f (E[TeV]/ADCch.) \n",fEnCalibration[4]);
    printf("  ZEM2 = %.2f (E[TeV]/ADCch.) \n",fEnCalibration[5]);
  
-   printf("\n\n\n #######      Equalization coefficients #######       \n");
+   printf("\n\n ####### Equalization coefficients ####### \n");
    printf("  ZN1 -> %1.2f %1.2f %1.2f %1.2f %1.2f  \n",
     fZN1EqualCoeff[0],fZN1EqualCoeff[1],fZN1EqualCoeff[2],fZN1EqualCoeff[3],fZN1EqualCoeff[4]);
    printf("  ZP1 -> %1.2f %1.2f %1.2f %1.2f %1.2f  \n",
     fZP1EqualCoeff[0],fZP1EqualCoeff[1],fZP1EqualCoeff[2],fZP1EqualCoeff[3],fZP1EqualCoeff[4]);
    printf("  ZN2 -> %1.2f %1.2f %1.2f %1.2f %1.2f  \n",
     fZN2EqualCoeff[0],fZN2EqualCoeff[1],fZN2EqualCoeff[2],fZN2EqualCoeff[3],fZN2EqualCoeff[4]);
-   printf("  ZP2 -> %1.2f %1.2f %1.2f %1.2f %1.2f  \n\n",
+   printf("  ZP2 -> %1.2f %1.2f %1.2f %1.2f %1.2f  \n",
     fZP2EqualCoeff[0],fZP2EqualCoeff[1],fZP2EqualCoeff[2],fZP2EqualCoeff[3],fZP2EqualCoeff[4]);
+   printf("\n\n ####### Parameters from EZDC vs. ZEM correlation #######       \n");
+   printf("  ZEMEndPoint = %1.2f, ZEMCutFraction = %1.2f \n"
+     "  DZEMInf = %1.2f, DZEMSup = %1.2f\n",
+     fZEMEndValue, fZEMCutFraction, fDZEMInf, fDZEMSup);
+   printf("\n\n ####### Parameters from EZDC vs. Nspec correlation #######     \n");
+   printf("  EZN1MaxValue = %1.2f, EZP1MaxValue = %1.2f, EZDC1MaxValue = %1.2f \n"
+     "  EZN2MaxValue = %1.2f, EZP2MaxValue = %1.2f, EZDC2MaxValue = %1.2f \n\n",
+     fEZN1MaxValue, fEZP1MaxValue, fEZDC1MaxValue,
+     fEZN2MaxValue, fEZP2MaxValue, fEZDC2MaxValue);
+
 } 
 
 //________________________________________________________________
 void AliZDCCalibData::SetMeanPed(Float_t* MeanPed)
 {
-  if(MeanPed) for(int t=0; t<47; t++) fMeanPedestal[t] = MeanPed[t];
-  else for(int t=0; t<47; t++) fMeanPedestal[t] = 0.;
+  if(MeanPed) for(int t=0; t<48; t++) fMeanPedestal[t] = MeanPed[t];
+  else for(int t=0; t<48; t++) fMeanPedestal[t] = 0.;
 }
 //________________________________________________________________
 void AliZDCCalibData::SetMeanPedWidth(Float_t* MeanPedWidth)
 {
-  if(MeanPedWidth) for(int t=0; t<47; t++) fMeanPedWidth[t] = MeanPedWidth[t];
-  else for(int t=0; t<47; t++) fMeanPedWidth[t] = 0.;
+  if(MeanPedWidth) for(int t=0; t<48; t++) fMeanPedWidth[t] = MeanPedWidth[t];
+  else for(int t=0; t<48; t++) fMeanPedWidth[t] = 0.;
 }
 
 //________________________________________________________________
 void AliZDCCalibData::SetOOTPed(Float_t* OOTPed)
 {
-  if(OOTPed) for(int t=0; t<44; t++) fOOTPedestal[t] = OOTPed[t];
-  else for(int t=0; t<44; t++) fOOTPedestal[t] = 0.;
+  if(OOTPed) for(int t=0; t<48; t++) fOOTPedestal[t] = OOTPed[t];
+  else for(int t=0; t<48; t++) fOOTPedestal[t] = 0.;
 }
 
 //________________________________________________________________
 void AliZDCCalibData::SetOOTPedWidth(Float_t* OOTPedWidth)
 {
-  if(OOTPedWidth) for(int t=0; t<44; t++) fOOTPedWidth[t] = OOTPedWidth[t];
-  else for(int t=0; t<44; t++) fOOTPedWidth[t] = 0.;
+  if(OOTPedWidth) for(int t=0; t<48; t++) fOOTPedWidth[t] = OOTPedWidth[t];
+  else for(int t=0; t<48; t++) fOOTPedWidth[t] = 0.;
 }
 
 //________________________________________________________________
@@ -210,13 +197,13 @@ void AliZDCCalibData:: SetPedCorrCoeff(Float_t* PedCorrCoeff)
   // Set coefficients for pedestal correlations
   if(PedCorrCoeff){
     for(Int_t j=0; j<2; j++){
-     for(int t=0; t<44; t++)
+     for(int t=0; t<48; t++)
        fPedCorrCoeff[j][t] = PedCorrCoeff[t];
     }
   }
   else{
     for(Int_t j=0; j<2; j++){
-     for(int t=0; t<44; t++)
+     for(int t=0; t<48; t++)
        fPedCorrCoeff[j][t] = 0.;
     }
   }
@@ -228,13 +215,13 @@ void AliZDCCalibData:: SetPedCorrCoeff(Float_t* PedCorrCoeff0, Float_t* PedCorrC
 {
   // Set coefficients for pedestal correlations
   if(PedCorrCoeff0 && PedCorrCoeff1){
-    for(int t=0; t<44; t++){
+    for(int t=0; t<48; t++){
        fPedCorrCoeff[0][t] = PedCorrCoeff0[t];
        fPedCorrCoeff[0][t] = PedCorrCoeff1[t];
     }
   }
   else{
-     for(int t=0; t<44; t++){
+     for(int t=0; t<48; t++){
        fPedCorrCoeff[0][t] = 0.;
        fPedCorrCoeff[1][t] = 0.;
     }
@@ -250,14 +237,6 @@ void AliZDCCalibData::SetEnCalib(Float_t* EnCalib)
   else for(int t=0; t<6; t++) fEnCalibration[t] = 0.;
 }
 
-//________________________________________________________________
-void AliZDCCalibData::SetPMTHVVal(Float_t* HVVal)
-{
-  // Set PMTs HV values
-  if(HVVal) for(int t=0; t<22; t++) fPMTHVVal[t] = HVVal[t];
-  else for(int t=0; t<22; t++) fPMTHVVal[t] = 0.;
-}
 //________________________________________________________________
 void AliZDCCalibData::SetZN1EqualCoeff(Float_t* EqualCoeff)
 {