float channel width
authoralla <alla@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 6 Nov 2007 12:40:27 +0000 (12:40 +0000)
committeralla <alla@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 6 Nov 2007 12:40:27 +0000 (12:40 +0000)
T0/AliT0Digitizer.cxx
T0/AliT0Parameters.h
T0/AliT0Reconstructor.cxx

index ebbc491..3caa768 100644 (file)
@@ -165,7 +165,7 @@ void AliT0Digitizer::Exec(Option_t* /*option*/)
 
   
   Int_t ph2Mip = fParam->GetPh2Mip();     
-  Int_t channelWidth = fParam->GetChannelWidth() ;  
+  Float_t channelWidth = fParam->GetChannelWidth() ;  
   Float_t delayVertex = fParam->GetTimeDelayTVD();
    
   zdetC = TMath::Abs(fParam->GetZPosition("T0/C/PMT1"));
@@ -287,7 +287,7 @@ void AliT0Digitizer::Exec(Option_t* /*option*/)
          sl = gr->Eval(qt);
 
          trLED = Int_t(( timeGaus[i] + 1000*sl )/channelWidth);
-         qtCh=Int_t (1000.*TMath::Log(qt)) / channelWidth;
+         qtCh=Int_t (1000.*TMath::Log(qt) / channelWidth);
          fADC0->AddAt(0,i);
          fADC->AddAt(qtCh,i);
          ftimeLED->AddAt(trLED,i); 
index 21cff17..f51a59f 100644 (file)
@@ -36,7 +36,7 @@ public:
   // Set various `Fixed' parameters 
   void SetPh2Mip(Int_t r=300)          { fPh2Mip = r; }
   void SetmV2Mip(Int_t r=50)          { fmV2Mip = r; }
-  void SetChannelWidth(Int_t s=24.4)   { fChannelWidth = s;}
+  void SetChannelWidth(Float_t s=24.4)   { fChannelWidth = s;}
   void SetmV2channel(Int_t size=320) { fmV2Channel = size; }
   void SetQTmin(Int_t qt=13) {fQTmin = qt;}
   void SetQTmax(Int_t qt=125) {fQTmax = qt;}
@@ -51,7 +51,7 @@ public:
   // Get `Fixed' various parameters
   Int_t GetPh2Mip()          const { return fPh2Mip; }
   Int_t GetmV2Mip()          const { return fmV2Mip; }
-  Int_t GetChannelWidth()     const { return fChannelWidth; }
+  Float_t GetChannelWidth()     const { return fChannelWidth; }
   Int_t GetmV2channel()     const { return fmV2Channel; }
   Int_t GetQTmin() const {return fQTmin;}
   Int_t GetQTmax() const {return fQTmax;}
@@ -89,7 +89,7 @@ public:
   Float_t   fT0zPosition[2] ;  // z-position of the two T0s
   Int_t     fPh2Mip;            // # photoelectrons per MIP in radiator
   Int_t     fmV2Mip;            // # mV per MIP in radiator
-  Int_t     fChannelWidth;          // channel width in ns   
+  Float_t     fChannelWidth;          // channel width in ns   
   Int_t     fmV2Channel;     // ADC mv  2  channel # (200000ps/(25*25).
   Int_t     fQTmin;                 //min  time for QTC
   Int_t     fQTmax;                 //max  time fro QTC 
index 9ae8e78..de8c0d1 100644 (file)
@@ -115,7 +115,7 @@ void AliT0Reconstructor::Reconstruct(TTree*digitsTree, TTree*clustersTree) const
 
   //  Int_t mV2Mip = param->GetmV2Mip();     
   //mV2Mip = param->GetmV2Mip();     
-  Int_t channelWidth = fParam->GetChannelWidth() ;  
+  Float_t channelWidth = fParam->GetChannelWidth() ;  
   Int_t meanT0 = fParam->GetMeanT0();
   
   AliDebug(1,Form("Start DIGITS reconstruction "));
@@ -243,7 +243,8 @@ void AliT0Reconstructor::Reconstruct(AliRawReader* rawReader, TTree*recTree) con
 
   //  Int_t mV2Mip = param->GetmV2Mip();     
   //mV2Mip = param->GetmV2Mip();     
-  Int_t channelWidth = fParam->GetChannelWidth() ;  
+  Float_t channelWidth = fParam->GetChannelWidth() ;  
+
   Int_t meanT0 = fParam->GetMeanT0();
  
   for (Int_t in=0; in<24; in++)
@@ -280,6 +281,8 @@ void AliT0Reconstructor::Reconstruct(AliRawReader* rawReader, TTree*recTree) con
       frecpoints->SetTime(ipmt,time[ipmt]);
       frecpoints->SetAmp(ipmt,adc[ipmt]);
       frecpoints->SetAmpLED(ipmt,qt);
+    AliDebug(1,Form(" QTC %f mv,  QTC  %f MIPS time in chann %f time %f ",adc[ipmt], adc[ipmt]/50.,time[ipmt], time[ipmt]*channelWidth));
+     
     }
     else {
       time[ipmt] = 0;
@@ -309,8 +312,7 @@ void AliT0Reconstructor::Reconstruct(AliRawReader* rawReader, TTree*recTree) con
   Float_t vertex = 0;
   if(besttimeA !=999999 && besttimeC != 999999 ){
     timeDiff = (besttimeC - besttimeA)*channelWidth;
-    //    meanTime = (besttimeA + besttimeC)/2.;
-    meanTime = (meanT0 - (besttimeA + besttimeC)/2) * channelWidth;
+    meanTime = (meanT0 - (besttimeA + besttimeC)/2) * channelWidth;   
     vertex = c*(timeDiff)/2. + (fdZ_A - fdZ_C)/2; //-(lenr-lenl))/2;
     AliDebug(1,Form("  timeDiff %f ps,  meanTime %f ps, vertex %f cm",timeDiff, meanTime,vertex ));
     frecpoints->SetVertex(vertex);
@@ -318,12 +320,14 @@ void AliT0Reconstructor::Reconstruct(AliRawReader* rawReader, TTree*recTree) con
     
   }
   //time in each channel as time[ipmt]-MeanTimeinThisChannel(with vertex=0)
+  /*
   for (Int_t ipmt=0; ipmt<24; ipmt++) {
     if(time[ipmt]>1) {
       time[ipmt] = (time[ipmt] - fTime0vertex[ipmt])*channelWidth;
       frecpoints->SetTime(ipmt,time[ipmt]);
     }
-  }
+    }
+  */
   recTree->Fill();