From 4c82df4c690fb5eb4c993ec2e78a57eb6a96ccf1 Mon Sep 17 00:00:00 2001
From: masera <masera@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Tue, 3 Jul 2007 15:57:59 +0000
Subject: [PATCH] New DA for SDD injectors (F.Prino)

---
 ITS/AliITSOnlineSDDInjectors.cxx |  10 +-
 ITS/AliITSOnlineSDDInjectors.h   |   2 +-
 ITS/ITSSDDINJda.cxx              | 205 +++++++++++++++++++++++++++++++
 3 files changed, 214 insertions(+), 3 deletions(-)
 create mode 100644 ITS/ITSSDDINJda.cxx

diff --git a/ITS/AliITSOnlineSDDInjectors.cxx b/ITS/AliITSOnlineSDDInjectors.cxx
index 1748abf2e40..e5e95aac128 100644
--- a/ITS/AliITSOnlineSDDInjectors.cxx
+++ b/ITS/AliITSOnlineSDDInjectors.cxx
@@ -17,6 +17,8 @@
 #include <TGraphErrors.h>
 #include <TMath.h>
 
+/* $Id$ */
+
 ///////////////////////////////////////////////////////////////////
 //                                                               //
 // Implementation of the class used for SDD injector analysis    //
@@ -364,6 +366,7 @@ void AliITSOnlineSDDInjectors::FindCentroids(){
 	fRMSCentroid[iii][ninj]=0.;
 	fGoodInj[iii][ninj]=0;
       }
+      if(fRMSCentroid[iii][ninj]==0) fGoodInj[iii][ninj]=0;
     }
   }
 }
@@ -382,11 +385,14 @@ void AliITSOnlineSDDInjectors::PrintCentroids(){
   }
 }
 //______________________________________________________________________
-void AliITSOnlineSDDInjectors::WriteToASCII(){
+void AliITSOnlineSDDInjectors::WriteToASCII(Int_t evNumb, UInt_t timeStamp, Int_t optAppend){
   //
   Char_t outfilnam[100];
   sprintf(outfilnam,"SDDinj_mod%03d_sid%d.data",fModuleId,fSide);
-  FILE* outf=fopen(outfilnam,"w");
+  FILE* outf;
+  if(optAppend==0) outf=fopen(outfilnam,"w");
+  else outf=fopen(outfilnam,"a");
+  fprintf(outf,"%d   %d   ",evNumb,timeStamp);
   for(Int_t ic=0;ic<fPolOrder+1;ic++){
     fprintf(outf,"%G ",fParam[ic]);
   }
diff --git a/ITS/AliITSOnlineSDDInjectors.h b/ITS/AliITSOnlineSDDInjectors.h
index c5cc8eff44f..298bc75b74e 100644
--- a/ITS/AliITSOnlineSDDInjectors.h
+++ b/ITS/AliITSOnlineSDDInjectors.h
@@ -45,7 +45,7 @@ class AliITSOnlineSDDInjectors : public AliITSOnlineSDD {
 
   void PrintInjMap();
   void PrintCentroids();
-  void WriteToASCII();
+  void WriteToASCII(Int_t evNumb, UInt_t timeStamp, Int_t optAppend=0);
 
   void Reset();      
   void AnalyzeEvent(TH2F* his);      
diff --git a/ITS/ITSSDDINJda.cxx b/ITS/ITSSDDINJda.cxx
new file mode 100644
index 00000000000..66c4ef2611f
--- /dev/null
+++ b/ITS/ITSSDDINJda.cxx
@@ -0,0 +1,205 @@
+//////////////////////////////////////////////////////////////////////////////
+// Detector Algorithm for analysis of SDD injector events.                  //
+//                                                                          //
+// Produces ASCII output files with:                                        //
+// 1. event number                                                          //
+// 2. event timestamp                                                       //
+// 3. Fit parameters of drift vel. vs. anode                                //
+// Tar Files are written to FXS                                             //
+//                                                                          //
+// Author: F. Prino (prino@to.infn.it)                                      //
+//                                                                          //
+//////////////////////////////////////////////////////////////////////////////
+
+
+
+extern "C" {
+#include "daqDA.h"
+}
+
+#include "event.h"
+#include "monitor.h"
+
+#include <stdio.h>
+#include <stdlib.h>
+
+// ROOT includes
+#include <TFile.h>
+#include <TH1F.h>
+#include <TH2F.h>
+#include <TSystem.h>
+
+// AliRoot includes
+#include "AliRawReaderDate.h"
+#include "AliITSOnlineSDDInjectors.h"
+#include "AliITSRawStreamSDD.h"
+/* Main routine
+      Arguments: list of DATE raw data files
+*/
+int main(int argc, char **argv) {
+
+  int status = 0;
+
+
+  /* log start of process */
+  printf("ITS SDD INJ algorithm program started\n");  
+
+
+  /* check that we got some arguments = list of files */
+  if (argc<2) {
+    printf("Wrong number of arguments\n");
+    return -1;
+  }
+
+
+  Int_t maxNEvents=10; // maximum number of events to be analyzed
+  const Int_t nSDDmodules=12;  // temp for test raw data
+  AliITSOnlineSDDInjectors **injan=new AliITSOnlineSDDInjectors*[2*nSDDmodules];
+  TH2F **histo=new TH2F*[2*nSDDmodules];
+  Int_t nWrittenEv[2*nSDDmodules];
+
+  Char_t hisnam[20];
+  for(Int_t imod=0; imod<nSDDmodules;imod++){
+    for(Int_t isid=0;isid<2;isid++){
+      Int_t index=2*imod+isid;
+      injan[index]=new AliITSOnlineSDDInjectors(imod,isid);
+      sprintf(hisnam,"his%03ds%d",imod,isid);
+      histo[index]=new TH2F(hisnam,"",256,-0.5,255.5,256,-0.5,255.5);
+      nWrittenEv[index]=0;
+    }
+  }
+  
+  /* report progress */
+  daqDA_progressReport(10);
+  Int_t iev=0;
+  /* read the data files */
+  int n;
+  for (n=1;n<argc;n++) {
+   
+    status=monitorSetDataSource( argv[n] );
+    if (status!=0) {
+      printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status));
+      return -1;
+    }
+
+    /* report progress */
+    /* in this example, indexed on the number of files */
+    // Progress report inside the event loop as well?
+    daqDA_progressReport(10+80*n/argc);
+
+    /* read the file */
+    for(;;) {
+      struct eventHeaderStruct *event;
+      eventTypeType eventT;
+
+      /* get next event */
+      status=monitorGetEventDynamic((void **)&event);
+      if (status==MON_ERR_EOF) break; /* end of monitoring file has been reached */
+      if (status!=0) {
+	printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));
+	return -1;
+      }
+
+      /* retry if got no event */
+      if (event==NULL) {
+	break;
+      }
+
+      iev++; 
+      if(iev>maxNEvents) break;
+      
+      /* use event - here, just write event id to result file */
+      eventT=event->eventType;
+      switch (event->eventType){
+	
+	/* START OF RUN */
+      case START_OF_RUN:
+	break;
+	
+	/* END OF RUN */
+      case END_OF_RUN:
+	break;
+	
+	//      case PHYSICS_EVENT:  // comment this line for test raw data
+	//	break;               // comment this line for test raw data
+
+      case CALIBRATION_EVENT:
+	break;  // uncomment this line for test raw data
+      case PHYSICS_EVENT: // uncomment this line for test raw data
+	printf(" event number = %i \n",iev);
+	AliRawReader *rawReader = new AliRawReaderDate((void*)event);
+	rawReader->RequireHeader(kFALSE);
+	rawReader->SelectEquipment(17,101,101);  // temp for test raw data
+
+	Int_t evtyp=0;
+	while(rawReader->ReadHeader()){
+	  const UInt_t *subev = rawReader->GetSubEventAttributes();
+	  if(subev[0]==0 && subev[1]==0 && subev[2]==0) evtyp=1; 
+	}
+
+	rawReader->Reset();
+	for(Int_t imod=0; imod<nSDDmodules;imod++){
+	  for(Int_t isid=0; isid<2;isid++){
+	    Int_t index=2*imod+isid;
+	    histo[index]->Reset();
+	  }
+	}
+	AliITSRawStreamSDD s(rawReader);
+	rawReader->SelectEquipment(17,101,101);  // temp for test raw data
+	
+	while(s.Next()){
+	  Int_t iddl=rawReader->GetDDLID();
+	  iddl=0; // temporary for test raw data
+	  Int_t isddmod=s.GetModuleNumber(iddl,s.GetCarlosId()); 
+	  isddmod-=240;  // to have SDD modules from 0 to 259
+	  isddmod=s.GetCarlosId(); // temporary for test raw data
+	  if(isddmod<nSDDmodules){ 
+	    Int_t index=2*isddmod+s.GetChannel(); 
+	    histo[index]->Fill(s.GetCoord2(),s.GetCoord1(),s.GetSignal());
+	  }
+	}
+	delete rawReader;
+	//	    UInt_t timeSt=rawReader->GetTimestamp();
+	UInt_t timeSt=iev*5000+12;  // fake timestamp for test
+	for(Int_t imod=0; imod<nSDDmodules;imod++){
+	  for(Int_t isid=0; isid<2;isid++){
+	    Int_t index=2*imod+isid;
+	    injan[index]->Reset();
+	    injan[index]->AnalyzeEvent(histo[index]);    
+	    injan[index]->WriteToASCII(iev,timeSt,nWrittenEv[index]);
+	    nWrittenEv[index]++;
+	  }
+	}
+	
+	/* free resources */
+	free(event);
+      }
+    }
+    
+  }
+    
+  Char_t filnam[100],command[120];
+  for(Int_t imod=0; imod<nSDDmodules;imod++){
+    for(Int_t isid=0; isid<2;isid++){
+      sprintf(filnam,"SDDinj_mod%03d_sid%d.data",imod,isid);
+      sprintf(command,"tar -rf SDDinj_LDC1.tar %s",filnam);
+      gSystem->Exec(command);
+    }  
+  }
+
+  /* write report */
+  printf("Run #%s, received %d injector events\n",getenv("DATE_RUN_NUMBER"),iev);
+
+  /* report progress */
+  daqDA_progressReport(90);
+
+
+  status=daqDA_FES_storeFile("./SDDinj_LDC1.tar","SDD_Calib");
+
+  /* report progress */
+  daqDA_progressReport(100);
+
+
+
+  return status;
+}
-- 
2.43.0