void AnalyzeSDDGainAllMod(Char_t *datafil,
Int_t adcfreq=20,
Int_t nDDL=0,
- Int_t firstEv=10,
- Int_t lastEv=16,
- Float_t pascalDAC=100,
- Int_t dataformat=1){
+ Int_t firstEv=18,
+ Int_t lastEv=22,
+ Float_t pascalDAC=100){
const Int_t kTotDDL=24;
const Int_t kModPerDDL=12;
const Int_t kSides=2;
+ Bool_t writtenoutput=kFALSE;
TH2F** histo = new TH2F*[kTotDDL*kModPerDDL*kSides];
}
}
}
- AliITSRawStream* s;
- if(dataformat==0){
- s=new AliITSRawStreamSDD(rd);
- }else{
- s=new AliITSRawStreamSDDCompressed(rd);
- if(dataformat==1) s->SetADCEncoded(kTRUE);
+
+ UChar_t cdhAttr=AliITSRawStreamSDD::ReadBlockAttributes(rd);
+ UInt_t amSamplFreq=AliITSRawStreamSDD::ReadAMSamplFreqFromCDH(cdhAttr);
+ AliITSRawStream* s=AliITSRawStreamSDD::CreateRawStreamSDD(rd,cdhAttr);
+ if(!writtenoutput){
+ printf("Use %s raw stream, sampling frequency %d MHz\n",s->ClassName(),amSamplFreq);
+ writtenoutput=kTRUE;
}
while(s->Next()){
Int_t iDDL=rd->GetDDLID();
}
-void AnalyzeSDDGainAllMod(Int_t nrun, Int_t n2, Char_t* dir="LHC08d_SDD",
+void AnalyzeSDDGainAllMod(Int_t nrun, Int_t n2, Int_t year=2009, Char_t* dir="LHC09b_SDD",
Int_t adcfreq=20,
Int_t nDDL=0,
- Int_t firstEv=15,
- Int_t lastEv=20,
- Float_t pascalDAC=100,
- Int_t dataformat=1){
+ Int_t firstEv=18,
+ Int_t lastEv=22,
+ Float_t pascalDAC=100){
+
TGrid::Connect("alien:",0,0,"t");
Char_t filnam[200];
- sprintf(filnam,"alien:///alice/data/2008/%s/%09d/raw/08%09d%03d.10.root",dir,nrun,nrun,n2);
+ sprintf(filnam,"alien:///alice/data/%d/%s/%09d/raw/%02d%09d%03d.10.root",year,dir,nrun,year-2000,nrun,n2);
printf("Open file %s\n",filnam);
- AnalyzeSDDGainAllMod(filnam,adcfreq,nDDL,firstEv,lastEv,pascalDAC,dataformat);
+ AnalyzeSDDGainAllMod(filnam,adcfreq,nDDL,firstEv,lastEv,pascalDAC);
}