]> git.uio.no Git - u/mrichter/AliRoot.git/commitdiff
updates to comply with AliTOFGeometryV5 becoming AliTOFGeometry and Calibration objects
authorzampolli <zampolli@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 4 Oct 2007 13:32:16 +0000 (13:32 +0000)
committerzampolli <zampolli@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 4 Oct 2007 13:32:16 +0000 (13:32 +0000)
TOF/AliTOFDigitMap.cxx
TOF/AliTOFDigitMap.h
TOF/AliTOFDigitizer.cxx
TOF/AliTOFDigitizer.h
TOF/AliTOFHitMap.cxx
TOF/AliTOFHitMap.h
TOF/AliTOFPreprocessor.cxx
TOF/AliTOFPreprocessor.h

index 56c47786e9c11b9d9ec4104f132967878c42aef8..d120be1ba3e142321247ca523d74974ab4ee33da 100644 (file)
@@ -15,6 +15,9 @@
 
 /*
 $Log$
+Revision 1.12  2007/02/20 15:57:00  decaro
+Raw data update: to read the TOF raw data defined in UNPACKED mode
+
 */
 
 ////////////////////////////////////////////////////////////////////////
@@ -52,8 +55,7 @@ AliTOFDigitMap::AliTOFDigitMap():
   fNpx(-1),
   fNpz(-1),
   fMaxIndex(-1),
-  fDigitMap(0x0),
-  fTOFGeometry(new AliTOFGeometry())
+  fDigitMap(0x0)
 {
 //
 // Default ctor
@@ -61,7 +63,7 @@ AliTOFDigitMap::AliTOFDigitMap():
 
   fNSector = AliTOFGeometry::NSectors();
   fNplate = AliTOFGeometry::NPlates();
-  fNstrip = fTOFGeometry->NStripC();//fTOFGeometry->NMaxNstrip();
+  fNstrip = AliTOFGeometry::NStripC();
   fNpx  = AliTOFGeometry::NpadX();
   fNpz  = AliTOFGeometry::NpadZ();
   fMaxIndex=fNSector*fNplate*fNstrip*fNpx*fNpz;
@@ -80,8 +82,7 @@ AliTOFDigitMap::AliTOFDigitMap(const AliTOFDigitMap & /*digitMap*/):
   fNpx(-1),
   fNpz(-1),
   fMaxIndex(-1),
-  fDigitMap(0x0),
-  fTOFGeometry(new AliTOFGeometry())
+  fDigitMap(0x0)
 {
 //
 // Dummy copy constructor
@@ -101,7 +102,6 @@ AliTOFDigitMap::~AliTOFDigitMap()
     for (Int_t i=0; i<fMaxIndex; i++)  delete[] fDigitMap[i];
   }
 
-  fTOFGeometry = 0;
 
 }
 
index 64bfdf4b5e64fab0d3b7750a7acc0d3f6ecbc6d0..a5256606d31f5a1561bb80d5165ed61213432c76 100644 (file)
@@ -21,8 +21,6 @@
 
 #include "AliHitMap.h"
 
-class AliTOFGeometry;
-
 class AliTOFDigitMap : public TObject
 {
  public:
@@ -66,8 +64,6 @@ class AliTOFDigitMap : public TObject
     Int_t fMaxIndex;                      // maximum index in hit map
     Int_t **fDigitMap;                    // ! [fMaxIndex][kMaxDigitsPerPad]
 
-    AliTOFGeometry *fTOFGeometry;         // Pointer to the TOF geometry
-
-    ClassDef(AliTOFDigitMap,1) // Implements DigitMap as a 1-dim array
+    ClassDef(AliTOFDigitMap,2) // Implements DigitMap as a 1-dim array
 };
 #endif 
index 4dfb1bed856b3a183def5c2d6dac0d6eeb365697..be0928aca97ffeee9949ed6cf1a4ddbbd0c07cb9 100644 (file)
@@ -31,6 +31,7 @@
 #include "TH1F.h"
 #include "TTree.h"
 #include "TRandom.h"
+#include "TObjArray.h"
 
 #include "AliLoader.h"
 #include "AliLog.h"
 #include "AliRunLoader.h"
 #include "AliRun.h"
 
-#include "AliTOFCal.h"
 #include "AliTOFcalib.h"
-#include "AliTOFChannel.h"
+#include "AliTOFChannelOnline.h"
+#include "AliTOFChannelOffline.h"
 #include "AliTOFDigitizer.h"
 #include "AliTOFdigit.h"
 #include "AliTOFHitMap.h"
-#include "AliTOFGeometryV5.h"
+#include "AliTOFGeometry.h"
 #include "AliTOFSDigit.h"
 #include "AliTOF.h"
 
@@ -60,7 +61,6 @@ ClassImp(AliTOFDigitizer)
 //___________________________________________
   AliTOFDigitizer::AliTOFDigitizer()  :
     AliDigitizer(),
-    fGeom(0x0),
     fDigits(0x0),
     fSDigitsArray(0x0),
     fhitMap(0x0)
@@ -71,7 +71,6 @@ ClassImp(AliTOFDigitizer)
 //___________________________________________
 AliTOFDigitizer::AliTOFDigitizer(AliRunDigitizer* manager): 
   AliDigitizer(manager), 
-  fGeom(0x0),
   fDigits(0x0),
   fSDigitsArray(0x0),
   fhitMap(0x0)
@@ -82,7 +81,6 @@ AliTOFDigitizer::AliTOFDigitizer(AliRunDigitizer* manager):
 //------------------------------------------------------------------------
 AliTOFDigitizer::AliTOFDigitizer(const AliTOFDigitizer &source):
   AliDigitizer(source),
-  fGeom(0x0), 
   fDigits(0),
   fSDigitsArray(0),
   fhitMap(0)
@@ -91,7 +89,6 @@ AliTOFDigitizer::AliTOFDigitizer(const AliTOFDigitizer &source):
   this->fDigits=source.fDigits;
   this->fSDigitsArray=source.fSDigitsArray;
   this->fhitMap=source.fhitMap;
-  this->fGeom=source.fGeom; 
 
 }
 
@@ -102,7 +99,6 @@ AliTOFDigitizer::AliTOFDigitizer(const AliTOFDigitizer &source):
   this->fDigits=source.fDigits;
   this->fSDigitsArray=source.fSDigitsArray;
   this->fhitMap=source.fhitMap;
-  this->fGeom=source.fGeom; 
   return *this;
 
 }
@@ -146,21 +142,25 @@ void AliTOFDigitizer::Exec(Option_t* /*option*/)
      return;
    }
    
+  /*
   outrl->CdGAFile();
   TFile *in=(TFile*)gFile;
   TDirectory *savedir=gDirectory;
 
+   
+  //when fGeom was needed
+
   if (!in->IsOpen()) {
     AliWarning("Geometry file is not open default  TOF geometry will be used");
-    fGeom = new AliTOFGeometryV5();
+    fGeom = new AliTOFGeometry();
   }
   else {
     in->cd();
     fGeom = (AliTOFGeometry*)in->Get("TOFgeometry");
   }
-
+  
   savedir->cd();
-
+  */
   AliLoader* outgime = outrl->GetLoader("TOFLoader");
   if (outgime == 0x0)
    {
@@ -181,7 +181,7 @@ void AliTOFDigitizer::Exec(Option_t* /*option*/)
   fSDigitsArray=new TClonesArray("AliTOFSDigit",1000);
   
   // create hit map (to be created in Init())
-  fhitMap = new AliTOFHitMap(fSDigitsArray, fGeom);
+  fhitMap = new AliTOFHitMap(fSDigitsArray);
   
   // Loop over files to digitize
 
@@ -284,8 +284,8 @@ void AliTOFDigitizer::CreateDigits()
   } // end loop on sdigits - end digitizing all collected sdigits
 
   //Insert Decalibration 
-
-  AliTOFcalib * calib = new AliTOFcalib(fGeom);
+  AliInfo("in digitizer, create digits");
+  AliTOFcalib * calib = new AliTOFcalib();
   InitDecalibration(calib);
   DecalibrateTOFSignal(calib);
   delete calib;
@@ -406,16 +406,20 @@ void AliTOFDigitizer::CollectSDigit(AliTOFSDigit * sdigit)
 
 //_____________________________________________________________________________
 void AliTOFDigitizer::InitDecalibration( AliTOFcalib *calib) const {
-  calib->ReadSimParFromCDB("TOF/Calib", -1); // use AliCDBManager's number
+  calib->CreateSimCalArrays();
+  calib->ReadSimParOnlineFromCDB("TOF/Calib", -1); // use AliCDBManager's number
+  calib->ReadSimParOfflineFromCDB("TOF/Calib", -1); // use AliCDBManager's number
 }
 //---------------------------------------------------------------------
 void AliTOFDigitizer::DecalibrateTOFSignal( AliTOFcalib *calib){
 
   // Read Calibration parameters from the CDB
 
-  AliTOFCal * cal= calib->GetTOFCalSimArray();
+  TObjArray * calOnline= calib->GetTOFSimCalArrayOnline();
+  TObjArray * calOffline= calib->GetTOFSimCalArrayOffline();
 
-  AliDebug(2,Form("Size of AliTOFCal = %i",cal->NPads()));
+  AliDebug(2,Form("Size of array for Online Calibration = %i",calOnline->GetEntries()));
+  AliDebug(2,Form("Size of array for Offline Calibration = %i",calOffline->GetEntries()));
 
   // Initialize Quantities to Simulate ToT Spectra
 
@@ -454,7 +458,7 @@ void AliTOFDigitizer::DecalibrateTOFSignal( AliTOFcalib *calib){
   Bool_t isToTSimulated=kFALSE;
   Bool_t misCalibPars=kFALSE;
   if(hToT->GetEntries()>0)isToTSimulated=kTRUE;  
-  Int_t ndigits = fDigits->GetEntriesFast();    
+  Int_t ndigits = fDigits->GetEntriesFast();
   for (Int_t i=0;i<ndigits;i++){
     AliTOFdigit * dig = (AliTOFdigit*)fDigits->At(i);
     Int_t detId[5];
@@ -469,16 +473,17 @@ void AliTOFDigitizer::DecalibrateTOFSignal( AliTOFcalib *calib){
       //A realistic ToT Spectrum was found in input, 
       //decalibrated TOF Digits likely to be simulated....
  
-      Int_t index = calib->GetIndex(detId); // The channel index    
-      AliTOFChannel *calChannel = cal->GetChannel(index); //retrieve the info
-      Float_t timedelay = calChannel->GetDelay(); //The global channel delay
-      Float_t par[6];  // time slewing parameters
+      Int_t index = AliTOFGeometry::GetIndex(detId); // The channel index    
+      AliTOFChannelOnline *calChannelOnline = (AliTOFChannelOnline *)calOnline->At(index); //retrieve the info for time delay
+      AliTOFChannelOffline *calChannelOffline = (AliTOFChannelOffline *)calOffline->At(index); //retrieve the info for time slewing 
+      Double_t timedelay = (Double_t)calChannelOnline->GetDelay(); //The global channel delay, ns
+      Double_t par[6];  // time slewing parameters
   
       //check whether we actually ask for miscalibration
 
       if(timedelay!=0)misCalibPars=kTRUE;
       for (Int_t j = 0; j<6; j++){
-       par[j]=calChannel->GetSlewPar(j);
+       par[j]=(Double_t)calChannelOffline->GetSlewPar(j);
        if(par[j]!=0)misCalibPars=kTRUE;
       }
       AliDebug(2,Form(" Calib Pars = %f, %f, %f, %f, %f, %f ",par[0],par[1],par[2],par[3],par[4],par[5]));
@@ -502,16 +507,17 @@ void AliTOFDigitizer::DecalibrateTOFSignal( AliTOFcalib *calib){
       tToT= (Double_t) trix; // to apply slewing we start from ns..
       // transform TOF signal in ns
       AliDebug(2,Form(" The Initial Time (counts): %i: ",dig->GetTdc()));
-      AliDebug(2,Form(" Time before miscalibration (ps) %e: ",dig->GetTdc()*AliTOFGeometry::TdcBinWidth()));
+      AliDebug(2,Form(" Time before miscalibration (ps) %e: ",dig->GetTdc()*(Double_t)AliTOFGeometry::TdcBinWidth()));
       // add slewing effect
       timeCorr=par[0] + tToT*(par[1] +tToT*(par[2] +tToT*(par[3] +tToT*(par[4] +tToT*par[5])))); 
       AliDebug(2,Form(" The Time slewing (ns): %f: ",timeCorr));
+      AliDebug(2,Form(" The Time delay (ns): %f: ",timedelay));
       // add global time shift
       timeCorr = timeCorr + timedelay;
-      AliDebug(2,Form(" The Time Slewing+ delay (ns): %f: ",timeCorr));
+      AliDebug(2,Form(" The Time Slewing + delay (ns): %f: ",timeCorr));
       //convert to ps
       timeCorr*=1E3;
-      Double_t timeMis = (Double_t)(dig->GetTdc())*AliTOFGeometry::TdcBinWidth();
+      Double_t timeMis = (Double_t)(dig->GetTdc())*(Double_t)AliTOFGeometry::TdcBinWidth();
       timeMis = timeMis+timeCorr;
       AliDebug(2,Form(" The Miscalibrated time (ps): %e: ",timeMis));
 
@@ -525,7 +531,7 @@ void AliTOFDigitizer::DecalibrateTOFSignal( AliTOFcalib *calib){
       tToT*=1E3; //back to ps  
       Int_t tot=(Int_t)(tToT/AliTOFGeometry::ToTBinWidth());//(factor 1E3 as input ToT is in ns)
       dig->SetToT(tot); 
-      AliDebug(2,Form(" Final Time and ToT (counts): %i: ",dig->GetTdc(),dig->GetToT()));
+      AliDebug(2,Form(" Final Time and ToT (counts): %i: , %i:",dig->GetTdc(),dig->GetToT()));
       if(tdcCorr<0){
        AliWarning (Form(" The bad Slewed Time(TDC counts)= %i ", tdcCorr)); 
        AliWarning(Form(" The bad ToT (TDC counts)= %i ", tot)); 
index e0229deb225e0d645ed12752c2371407cf48be53..003992e7ec683c50f29cb05f193ca798d06bea44 100644 (file)
@@ -23,7 +23,6 @@ class AliRunDigitizer;
 class AliTOFHitMap;
 class AliTOFSDigit;
 class AliTOFcalib;
-class AliTOFGeometry;
 
 class AliTOFDigitizer : public AliDigitizer {
  public:
@@ -49,13 +48,12 @@ class AliTOFDigitizer : public AliDigitizer {
                                                     // due to the low
                                                     // noise expected
                                                     // level
-  AliTOFGeometry *fGeom;       // Pointer to the TOF geometry
   TClonesArray *fDigits;       //! array with digits
   TClonesArray *fSDigitsArray; //! List of summable digits; used as a
                               //container for all sdigits to be
                               //merged
   AliTOFHitMap *fhitMap ;      //! hit map used to perform the merging
   
-  ClassDef(AliTOFDigitizer,0)  // TOF/Merging/Digitization
+  ClassDef(AliTOFDigitizer,1)  // TOF/Merging/Digitization
 };    
 #endif
index f8b759d4cd3b4e14c5cfe9eb0097a92c8f108e9d..a52af1dc7c86e871c4dd6058fe3811e04e4c71a4 100644 (file)
@@ -50,8 +50,7 @@ AliTOFHitMap::AliTOFHitMap():
   fNpz(-1),
   fSDigits(0x0),
   fMaxIndex(-1),
-  fHitMap(0x0),
-  fTOFGeometry(0x0)
+  fHitMap(0x0)
 {
 //
 // Default ctor
@@ -59,7 +58,7 @@ AliTOFHitMap::AliTOFHitMap():
 }
 
 ////////////////////////////////////////////////////////////////////////
-AliTOFHitMap::AliTOFHitMap(TClonesArray *dig, AliTOFGeometry *tofGeom):
+AliTOFHitMap::AliTOFHitMap(TClonesArray *dig):
   fNSector(-1),
   fNplate(-1),
   fNstrip(-1),
@@ -67,8 +66,7 @@ AliTOFHitMap::AliTOFHitMap(TClonesArray *dig, AliTOFGeometry *tofGeom):
   fNpz(-1),
   fSDigits(dig),
   fMaxIndex(-1),
-  fHitMap(0x0),
-  fTOFGeometry(tofGeom)
+  fHitMap(0x0)
 {
 //
 // ctor
@@ -79,7 +77,7 @@ AliTOFHitMap::AliTOFHitMap(TClonesArray *dig, AliTOFGeometry *tofGeom):
 
   fNSector = AliTOFGeometry::NSectors();
   fNplate = AliTOFGeometry::NPlates();
-  fNstrip = fTOFGeometry->NStripC();//fTOFGeometry->NMaxNstrip();
+  fNstrip = AliTOFGeometry::NStripC();//fTOFGeometry->NMaxNstrip();
   fNpx  = AliTOFGeometry::NpadX();
   fNpz  = AliTOFGeometry::NpadZ();
   fMaxIndex=fNSector*fNplate*fNstrip*fNpx*fNpz;
@@ -97,8 +95,7 @@ AliTOFHitMap::AliTOFHitMap(const AliTOFHitMap & /*hitMap*/)
   fNpz(-1),
   fSDigits(0x0),
   fMaxIndex(-1),
-  fHitMap(0x0),
-  fTOFGeometry(0x0)
+  fHitMap(0x0)
 {
 //
 // Dummy copy constructor
@@ -115,8 +112,6 @@ AliTOFHitMap::~AliTOFHitMap()
 //
   delete[] fHitMap;
 
-  fTOFGeometry = 0x0;
-
 }
 
 ////////////////////////////////////////////////////////////////////////
index a16850f407b2a059033b11d0b18be6172cfff4e5..a917378bc0b0bd116714379d40e4a1825bb29a5f 100644 (file)
 #include "AliHitMap.h"
 
 class TClonesArray;
-class AliTOFGeometry;
 
 class AliTOFHitMap : public TObject
 {
  public:
     AliTOFHitMap();
-    AliTOFHitMap(TClonesArray *sdig, AliTOFGeometry *tofGeom);
+    AliTOFHitMap(TClonesArray *sdig);
     AliTOFHitMap(const AliTOFHitMap & hitMap);
     
     virtual ~AliTOFHitMap();
@@ -59,8 +58,6 @@ class AliTOFHitMap : public TObject
     Int_t fMaxIndex;                      // maximum index in hit map
     Int_t *fHitMap;                       // ! [fMaxIndex]         
 
-    AliTOFGeometry *fTOFGeometry;         // Pointer to the TOF geometry
-
-    ClassDef(AliTOFHitMap,0) // Implements HitMap as a 1-dim array
+    ClassDef(AliTOFHitMap,1) // Implements HitMap as a 1-dim array
 };
 #endif 
index e19251d665617f5489d6bee96cb67e985b045dd0..7bccb5d810b73dd5a79841ec276be188bec5874c 100644 (file)
 
 #include "AliCDBMetaData.h"
 #include "AliLog.h"
-#include "AliTOFCalOnline.h"
 #include "AliTOFChannelOnline.h"
 #include "AliTOFDataDCS.h"
-#include "AliTOFGeometryV5.h"
+#include "AliTOFGeometry.h"
 #include "AliTOFPreprocessor.h"
 
 class TF1;
@@ -61,10 +60,11 @@ AliTOFPreprocessor::AliTOFPreprocessor(AliShuttleInterface* shuttle) :
   fData(0),
   fh2(0),
   fCal(0),
-  fTOFGeometry(0),
+  fNChannels(0),
   fStoreRefData(kTRUE)
 {
   // constructor
+
 }
 
 //_____________________________________________________________________________
@@ -76,10 +76,9 @@ AliTOFPreprocessor::~AliTOFPreprocessor()
   fData = 0;
   delete fh2;
   fh2 = 0;
+  fCal->Clear();
   delete fCal;
   fCal = 0;
-  delete fTOFGeometry;
-  fTOFGeometry = 0;
 }
 
 //______________________________________________________________________________
@@ -96,9 +95,13 @@ void AliTOFPreprocessor::Initialize(Int_t run, UInt_t startTime,
 
        fData = new AliTOFDataDCS(fRun, fStartTime, fEndTime);
        fh2 = 0x0;
-       fTOFGeometry = new AliTOFGeometryV5();
-       fCal = new AliTOFCalOnline(fTOFGeometry);
-       fCal->CreateArray();
+       fNChannels = AliTOFGeometry::NSectors()*(2*(AliTOFGeometry::NStripC()+AliTOFGeometry::NStripB())+AliTOFGeometry::NStripA())*AliTOFGeometry::NpadZ()*AliTOFGeometry::NpadX();
+       fCal = new TObjArray(fNChannels);
+       fCal->SetOwner();
+       for (Int_t ich = 0; ich<fNChannels; ich ++){
+         AliTOFChannelOnline * calChOnline = new AliTOFChannelOnline();
+         fCal->AddAt(calChOnline,ich);
+       }
 }
 
 //_____________________________________________________________________________
@@ -226,8 +229,7 @@ UInt_t AliTOFPreprocessor::Process(TMap* dcsAliasMap)
            static const Int_t kSize=fh2->GetNbinsX();
            static const Int_t kNBins=fh2->GetNbinsY();
            static const Double_t kXBinmin=fh2->GetYaxis()->GetBinLowEdge(1);
-           Int_t npads = fCal->NPads();
-           if (kSize != npads){
+           if (kSize != fNChannels){
              Log(" number of bins along x different from number of pads, found only a subset of the histograms, TOF exiting from Shuttle");
              delete daqFile;
               return 7; //return error code for histograms not existing/junky
@@ -266,9 +268,9 @@ UInt_t AliTOFPreprocessor::Process(TMap* dcsAliasMap)
              sumw2=sumw2/nent; //<x^2>
              Double_t rmsmean= 0;
              rmsmean = TMath::Sqrt((sumw2-mean*mean)/nent);
-             if (ich<npads) {
-               AliTOFChannelOnline * ch = fCal->GetChannel(ich);
-               ch->SetDelay(mean);
+             if (ich<fNChannels) {
+               AliTOFChannelOnline * ch = (AliTOFChannelOnline *)fCal->At(ich);
+               ch->SetDelay((Double_t)mean*AliTOFGeometry::TdcBinWidth()*1.E-3);  // delay in ns
              }
            delete h1;
            h1=0x0;
index c0c4a1a072f1f8793cea9a0a633cf9b37ba279f1..4be1e830fe805a6d91c200d5325077ecd29a7a97 100644 (file)
@@ -14,9 +14,6 @@
 
 class AliTOFDataDCS;
 class TObjArray;
-class AliTOFCalOnline;
-class AliTOFGeometryV5;
-class AliTOFGeometry;
 class TH2S;
 
 class AliTOFPreprocessor : public AliPreprocessor
@@ -38,12 +35,13 @@ class AliTOFPreprocessor : public AliPreprocessor
     static const Int_t fgkBinRangeAve;   // number of bins where to 
                                          // calculate the mean
     static const Double_t fgkThrPar;     // parameter used to trigger the 
-                                         //calculation of the delay
-    AliTOFDataDCS *fData;    // CDB class that stores the data
-    TH2S *fh2;       // TH2S from DAQ for histograms for delays  
-    AliTOFCalOnline *fCal;         // TOF Calibration object
-    AliTOFGeometry *fTOFGeometry;  // TOF Geometry version
-    Bool_t fStoreRefData;  // Flag to decide storage of Ref Data
-    ClassDef(AliTOFPreprocessor, 1);
+                                         // calculation of the delay
+    AliTOFDataDCS *fData;                // CDB class that stores the data
+    TH2S *fh2;                           // TH2S from DAQ for histograms 
+                                         // for delays  
+    TObjArray *fCal;                     // TOF Calibration object
+    Int_t fNChannels;                    // number of TOF channels
+    Bool_t fStoreRefData;                // Flag to decide storage of Ref Data
+    ClassDef(AliTOFPreprocessor, 2);
 };
 #endif