fSumMult(0),
fEff(0)
{
- // cout<<"AliSTARTDigitizer::AliSTARTDigitizer"<<endl;
// ctor which should be used
AliDebug(1,"processed");
AliLoader * pOutStartLoader = outRL->GetLoader("STARTLoader");
AliDebug(1,"start...");
- cout<<" AliSTARTDigitizer::Exec "<<endl;
//input loader
//
// From hits to digits
if (al>threshold) {
qt=Int_t (TMath::Log(al*ph2mV) * mV2channel);
qtAmp=Int_t (TMath::Log(al*10*ph2mV) * mV2channel);
- cout<<i<<" "<<qt<<" "<<qtAmp<<endl;
fADC->AddAt(qt,i);
ftimeTDC->AddAt(tr,i);
fADCAmp->AddAt(qtAmp,i);
// Float_t hc=197.326960*1.e6; //mev*nm
Double_t hc=1.973*1.e-6; //gev*nm
- // cout<<"AliSTARTDigitizer::RegisterPhotoE >> energy "<<energy<<endl;
Float_t lambda=hc/energy;
Int_t bin= fEff->GetXaxis()->FindBin(lambda);
Float_t eff=fEff->GetBinContent(bin);
#include "AliSTART.h"
#include "AliSTARTRawData.h"
#include "AliSTARTdigit.h"
-//#include "AliSTARTLoader.h"
-
-#include <AliLoader.h>
-#include <AliRunLoader.h>
-#include <TClonesArray.h>
#include <TTree.h>
#include "AliRawDataHeader.h"
//_____________________________________________________________________________
AliSTARTRawData::AliSTARTRawData():TObject()
{
-
+ /*
+- 48 channels (2 words each as in TOF DDL) for :
+word 1 :0-5bit number of PMT; word 2: 0-7 error sign, 8-31 TDC
+and the same but for amplified signal. Now I wrote the same time because
+CDF are not ready and differences didn't measured yet.
+
+- 96 channel for amplitude: very preliminary, QTC features are not
+known now, preliminary i put as T1 time signal for this PMT in first
+channel and T1+A in second, where A=Log(Amplitude);
+and the same for amplified but A=Log(10*Amplitude).
+
+- T0-A and T0-C 2 channels
+- T0A-T0C vertex information
+- Time Meaner where T0C TOF increase to the T0A TOF distance
+- 6 multiplicity signals the same way as amplitude and with the same
+uncertances
+ */
fIndex=-1;
fDigits = NULL;
- ftimeTDC = new TArrayI(24);
+ fTimeTDC = new TArrayI(24);
fADC = new TArrayI(24);
- ftimeTDCAmp = new TArrayI(24);
+ fTimeTDCAmp = new TArrayI(24);
fADCAmp = new TArrayI(24);
fSumMult = new TArrayI(6);
// this->Dump();
delete fDigits;
fDigits = NULL;
}
- delete ftimeTDC;
+ delete fTimeTDC;
delete fADC;
- delete ftimeTDCAmp;
+ delete fTimeTDCAmp;
delete fADCAmp;
delete fSumMult;
}
Int_t error=0;
// Get the digits array
- fDigits->GetTime(*ftimeTDC);
+ fDigits->GetTime(*fTimeTDC);
fDigits->GetADC(*fADC);
- fDigits->GetTimeAmp(*ftimeTDCAmp);
+ fDigits->GetTimeAmp(*fTimeTDCAmp);
fDigits->GetADCAmp(*fADCAmp);
fDigits->GetSumMult(*fSumMult);
// Loop through all PMT
for (Int_t det = 0; det < 24; det++) {
- Int_t time=ftimeTDC->At(det);
+ Int_t time=fTimeTDC->At(det);
Int_t qtc=fADC->At(det);
- Int_t timeAmp=ftimeTDCAmp->At(det);
+ Int_t timeAmp=fTimeTDCAmp->At(det);
Int_t qtcAmp=fADCAmp->At(det);
//conver ADC to time (preliminary algorithm)
}
//---------------------------------------------------------------------------------------
-Int_t AliSTARTRawData::RawDataSTART(AliSTARTdigit *fDigits){
-
- //This method creates the Raw data files for TOF detector
+Int_t AliSTARTRawData::RawDataSTART(AliSTARTdigit *fDigits)
+{
+ //This method creates the Raw data files for START detector
+
+
const Int_t kSize=512; //2*AliTOFGeometry::NpadXSector()
//max number of digits per DDL file times 2
UInt_t buf[kSize];
UInt_t word;
fIndex=-1;
- // TClonesArray*& digits = * (TClonesArray**) branch->GetAddress();
char fileName[15];
ofstream outfile; // logical name of the output file
#include "TObject.h"
class AliSTART;
class AliSTARTdigit;
-class TTRee;
-class TBranch;
+//class TTRee;
+//class TBranch;
//class AliRawDataHeader;
class AliSTARTRawData : public TObject {
//START digits arrays
- AliSTARTdigit *fDigits; //! The START digits manager
- TArrayI *timeTDC() {return ftimeTDC;}
+ TArrayI *TimeTDC() {return fTimeTDC;}
TArrayI *ADC() {return fADC;}
- TArrayI *timeTDCAmp() {return ftimeTDCAmp;}
+ TArrayI *TimeTDCAmp() {return fTimeTDCAmp;}
TArrayI *ADCAmp() {return fADCAmp;}
TArrayI *SumMult() {return fSumMult;}
Int_t fTimeDiff ; //time difference
Int_t fMeanTime ; // average time - ALICE start signal
- Int_t fBestTimeLeft;
- Int_t fBestTimeRight;
- TArrayI * ftimeTDC;
- TArrayI * fADC;
- TArrayI * ftimeTDCAmp;
- TArrayI * fADCAmp;
- TArrayI * fSumMult;
+ Int_t fBestTimeLeft; //first particle on the left
+ Int_t fBestTimeRight; //first particle on the right
+ TArrayI * fTimeTDC; //TDC on the each PMT
+ TArrayI * fADC; //QTC (ADC) on the each PMT
+ TArrayI * fTimeTDCAmp; // TDC with amplified signal
+ TArrayI * fADCAmp; //QTC amplified
+ TArrayI * fSumMult; //multiplicity: left, right, left&right + all amplified
+ AliSTARTdigit *fDigits; //! The START digits manager
ClassDef(AliSTARTRawData,1) // START raw data class