1st draft for pass0
[u/mrichter/AliRoot.git] / T0 / T0Physda.cxx
index 4e6921f..b643471 100644 (file)
@@ -125,6 +125,7 @@ int main(int argc, char **argv) {
    
   for(Int_t ic=0; ic<24; ic++) {
     hCFD1minCFD[ic] = new TH1F(Form("CFD1minCFD%d",ic+1),"CFD-CFD",kcbx,kclx,kcmx);
+    hCFD[ic] = new TH1F(Form("CFD1minCFD%d",ic+1),"CFD",kt0bx,kt0lx,kt0hx);
   }
   TH1F *hVertex = new TH1F("hVertex","T0 time",kt0bx,kt0lx,kt0hx);
   
@@ -183,13 +184,13 @@ int main(int argc, char **argv) {
       AliT0RawReader *start = new AliT0RawReader(reader, kTRUE);
       
       // Read raw data
-      Int_t allData[105][5];
-      for(Int_t i0=0;i0<105;i0++)
+      Int_t allData[110][5];
+      for(Int_t i0=0;i0<107;i0++)
        for(Int_t j0=0;j0<5;j0++)
          allData[i0][j0] = 0;
       
       if(start->Next()){
-       for (Int_t i=0; i<105; i++) {
+       for (Int_t i=0; i<107; i++) {
          for(Int_t iHit=0;iHit<5;iHit++){
            allData[i][iHit]= start->GetData(i,iHit);
          }
@@ -203,6 +204,7 @@ int main(int argc, char **argv) {
        Float_t meanShift[24];
        for (Int_t ik = 0; ik<24; ik++)
         { 
+          if(allData[ik+1][0]>0 ) hCFD[ik]->Fill(allData[ik+1][0]);
           if(ik<12 && allData[ik+1][0]>0 && allData[knpmtC][0]>0 ){
             hCFD1minCFD[ik]->Fill(allData[ik+1][0]-allData[knpmtC][0]);
           }
@@ -268,6 +270,8 @@ int main(int argc, char **argv) {
   for(Int_t j=0;j<24;j++){
     hCFD1minCFD[j]->SetDirectory(&hist);
     hCFD1minCFD[j]->Write();
+    hCFD[j]->Write();
+
   }
   hVertex->Write();
   hist.Close();