}
if(rawReader->GetType() != 7) return; // skips non physical triggers
AliDebug(1,"entering MakeRaws\n");
- rawReader->SelectEquipment(17,fgkeqOffset,fgkeqOffset + fgkDDLidRange);
rawReader->Reset();
AliITSRawStream *stream;
static const Int_t fgkmodoffset = 240; // number of SPD modules
static const Int_t fgknAnode = 256; // anode per half-module
static const Int_t fgknSide =2; // side per module
- static const Int_t fgkeqOffset = 256; // DDL offset
- static const Int_t fgkDDLidRange = 24; // number of DDL:so DDL range is 257-280
static const Int_t fgkDDLIDshift = 0; // necessary option until RawStream Table is complete
static const Int_t fgkLADDonLAY3 = 14; // number of ladder on layer 3
static const Int_t fgkLADDonLAY4 = 22; // number of ladder on layer 4
Bool_t fHLTMode; // kTRUE mode C kFALSE mode A
// Used in online mode only
AliITSHLTforSDD *fHLTSDD; // used for offline QA as the HLT mode flag
- ClassDef(AliITSQASDDDataMakerRec,5) // description
+ ClassDef(AliITSQASDDDataMakerRec,6) // description
};
void AnalyzeSDDGainAllMod(Char_t *datafil, Int_t nDDL, Int_t firstEv=10, Int_t lastEv=16, Float_t pascalDAC=100){
- Int_t eqOffset = 256;
- Int_t DDLid_range = 24;
- // Int_t eqOffset = 100;
- // Int_t DDLid_range = 1;
const Int_t kTotDDL=24;
const Int_t kModPerDDL=12;
const Int_t kSides=2;
c0->Clear();
c0->Divide(4,6,0.001,0.001);
printf("Event # %d\n",iev);
- //rd->SelectEvents(7);
- rd->SelectEquipment(17,eqOffset,eqOffset+DDLid_range);
rd->Reset();
for(Int_t iddl=0; iddl<kTotDDL;iddl++){
for(Int_t imod=0; imod<kModPerDDL;imod++){
}
}
AliITSRawStreamSDD s(rd);
- rd->SelectEquipment(17,eqOffset,eqOffset+DDLid_range);
while(s.Next()){
Int_t iDDL=rd->GetDDLID();
Int_t iCarlos=s.GetCarlosId();
void AnalyzeSDDInjectorsAllMod(Char_t *datafil, Int_t nDDL, Int_t firstEv=10, Int_t lastEv=15){
- Int_t eqOffset = 256;
- Int_t DDLid_range = 24;
- // Int_t eqOffset = 100;
- // Int_t DDLid_range = 1;
+
const Int_t kTotDDL=24;
const Int_t kModPerDDL=12;
const Int_t kSides=2;
sprintf(hisnam,"h%02dc%02ds%d",iddl,imod,isid);
histo[index]=new TH2F(hisnam,"",256,-0.5,255.5,256,-0.5,255.5);
anal[index]=new AliITSOnlineSDDInjectors(iddl,imod,isid);
+/* Uncomment these lines for analysis of runs with 40 MHz sapling */
+// anal[index]->SetInjLineRange(0,20,50);
+// anal[index]->SetInjLineRange(1,90,160);
+// anal[index]->SetInjLineRange(2,170,240);
+// anal[index]->SetTimeStep(25.);
+/* END of lines to be uncommented */
nWrittenEv[index]=0;
}
}
c1->Divide(4,6,0.001,0.001);
printf("Event # %d\n",iev);
UInt_t timeSt=rd->GetTimestamp();
- //rd->SelectEvents(7);
- rd->SelectEquipment(17,eqOffset,eqOffset+DDLid_range);
rd->Reset();
for(Int_t iddl=0; iddl<kTotDDL;iddl++){
for(Int_t imod=0; imod<kModPerDDL;imod++){
}
AliITSRawStreamSDD s(rd);
- rd->SelectEquipment(17,eqOffset,eqOffset+DDLid_range);
while(s.Next()){
Int_t iDDL=rd->GetDDLID();
Int_t iCarlos=s.GetCarlosId();
// Origin: F. Prino (prino@to.infn.it)
void AnalyzeSDDNoiseAllMod(Char_t *datafil, Int_t nDDL, Int_t firstEv=10, Int_t lastEv=12){
- Int_t eqOffset = 256;
- Int_t DDLid_range = 24;
- // Int_t eqOffset = 100;
- // Int_t DDLid_range = 1;
+
const Int_t kTotDDL=24;
const Int_t kModPerDDL=12;
const Int_t kSides=2;
c0->Clear();
c0->Divide(4,6,0.001,0.001);
printf("Event # %d ",iev);
- //rd->SelectEvents(7);
- rd->SelectEquipment(17,eqOffset,eqOffset+DDLid_range);
rd->Reset();
for(Int_t iddl=0; iddl<kTotDDL;iddl++){
for(Int_t imod=0; imod<kModPerDDL;imod++){
}
}
AliITSRawStreamSDD s(rd);
- rd->SelectEquipment(17,eqOffset,eqOffset+DDLid_range);
while(s.Next()){
Int_t iDDL=rd->GetDDLID();
Int_t iCarlos=s.GetCarlosId();
c0->Clear();
c0->Divide(4,6,0.001,0.001);
printf("Event # %d ",iev);
- rd2->SelectEquipment(17,eqOffset,eqOffset+DDLid_range);
- // rd2->SelectEvents(7);
rd2->Reset();
for(Int_t iddl=0; iddl<kTotDDL;iddl++){
for(Int_t imod=0; imod<kModPerDDL;imod++){
}
AliITSRawStreamSDD s(rd2);
- rd2->SelectEquipment(17,eqOffset,eqOffset+DDLid_range);
while(s.Next()){
Int_t iDDL=rd2->GetDDLID();
Int_t iCarlos=s.GetCarlosId();
void DisplaySDDRawData(Char_t *datafil, Int_t firstEv=0, Int_t lastEv=5){
- Int_t eqOffset = 256;
- Int_t DDLid_range= 24;
AliITSDDLModuleMapSDD* ddlmap=new AliITSDDLModuleMapSDD();
ddlmap->SetJun08Map();
evtime->Start();
printf("Event # %d\n",iev);
- rd->SelectEquipment(17,eqOffset,eqOffset+DDLid_range);
rd->Reset();
hzphi3->Reset();
hzphi4->Reset();
AliITSRawStreamSDD s(rd);
- rd->SelectEquipment(17,eqOffset,eqOffset+DDLid_range);
Int_t iCountNext=0;
while(s.Next()){
iCountNext++;
}
- Int_t eqOffset = 256;
- Int_t DDLrange = 24;
Int_t maxNEvents=18; // maximum number of events to be analyzed
const Int_t kTotDDL=24;
const Int_t kModPerDDL=12;
case PHYSICS_EVENT: // uncomment this line for test raw data
printf(" event number = %i \n",iev);
AliRawReader *rawReader = new AliRawReaderDate((void*)event);
- rawReader->SelectEquipment(17,eqOffset+1,eqOffset+DDLrange);
Int_t evtyp=0;
}
- Int_t eqOffset = 256;
- Int_t DDLrange = 24;
Int_t maxNEvents=15; // maximum number of events to be analyzed
const Int_t kTotDDL=24;
const Int_t kModPerDDL=12;
case PHYSICS_EVENT: // uncomment this line for test raw data
printf(" event number = %i \n",iev);
AliRawReader *rawReader = new AliRawReaderDate((void*)event);
- rawReader->SelectEquipment(17,eqOffset+1,eqOffset+DDLrange);
Int_t evtyp=0;
while(rawReader->ReadHeader()){
}
- Int_t eqOffset = 256;
- Int_t DDLrange = 24;
Int_t maxNEvents=10; // maximum number of events to be analyzed
const Int_t kTotDDL=24;
const Int_t kModPerDDL=12;
case PHYSICS_EVENT: // uncomment this line for test raw data
printf(" event number = %i \n",iev);
AliRawReader *rawReader = new AliRawReaderDate((void*)event);
- rawReader->SelectEquipment(17,eqOffset+1,eqOffset+DDLrange);
UInt_t timeSt=rawReader->GetTimestamp();
//UInt_t timeSt=iev*5000+12; // fake timestamp for test
delete initgeom;
*/
- Int_t eqOffset = 256;
- Int_t DDLid_range = 24;
do{ // start loop on events
if(++evCounter > MaxEvts) { cout << MaxEvts << " events read, stop" << endl; evCounter--; break; }
cout << "Read Event: " << evCounter+FirstEvt-1 << endl;
rd->RequireHeader(kFALSE);
- rd->SelectEquipment(17,eqOffset+1,eqOffset+DDLid_range);
rd->Reset(); // reset the current position to the beginning of the event
AliITSRawStreamSDD s(rd); //This class provides access to ITS SDD digits in raw data.
strFile += "?EventType=7";
AliRawReader *rd = new AliRawReaderDate(strFile.Data(),FirstEvt); // open run
Int_t evCounter = 0;
- Int_t eqOffset = 256;
- Int_t DDLid_range = 24;
do{ // start loop on events
if(++evCounter > MaxEvts) { cout << MaxEvts << " events read, stop" << endl; evCounter--; break; }
cout << "Read Event: " << evCounter+FirstEvt-1 << endl;
rd->RequireHeader(kFALSE);
- rd->SelectEquipment(17,eqOffset+1,eqOffset+DDLid_range); //17 states for "DRorc acquisition"
rd->Reset(); // reset the current position to the beginning of the event
Int_t nSkip = 0; // number of skipped signals
void PlotSDDRawData(Char_t datafil[100], Int_t nDDL, Int_t firstEv=0, Int_t lastEv=5){
- Int_t eqOffset = 256;
- Int_t DDLid_range= 24;
const Int_t nHybrids=24;
TH2F** histo = new TH2F*[nHybrids];
evtime->Start();
printf("Event # %d\n",iev);
- rd->SelectEquipment(17,eqOffset,eqOffset+DDLid_range);
rd->Reset();
for(Int_t i=0;i<nHybrids;i++) histo[i]->Reset();
AliITSRawStreamSDD s(rd);
- rd->SelectEquipment(17,eqOffset,eqOffset+DDLid_range);
Int_t iCountNext=0;
while(s.Next()){
iCountNext++;