Small additional update by Raphaelle

diff --git a/TRD/AliTRDCalibraFillHisto.cxx b/TRD/AliTRDCalibraFillHisto.cxx
index 24efd5bafe5f522ce1927fc2edb803714252fca6..e71cca374fb11eb2bef25540941884e6347f0d80 100644 (file)
--- a/TRD/AliTRDCalibraFillHisto.cxx
+++ b/TRD/AliTRDCalibraFillHisto.cxx
@@ -963,7 +963,7 @@ Int_t AliTRDCalibraFillHisto::ProcessEventDAQ(AliTRDRawStream *rawStream, Bool_t
       
       Int_t idetector = rawStream->GetDet();                            //  current detector
       if((fDetectorPreviousTrack != idetector) && (fDetectorPreviousTrack != -1)){
-       if(TMath::Mean(fTimeMax,phvalue)>0.0){
+       if(TMath::Mean(fTimeMax,phvalue)>13.0){
          withInput = 2;
          for(Int_t k = 0; k < fTimeMax; k++){
            UpdateDAQ(fDetectorPreviousTrack,0,0,k,phvalue[k],fTimeMax);
@@ -987,14 +987,14 @@ Int_t AliTRDCalibraFillHisto::ProcessEventDAQ(AliTRDRawStream *rawStream, Bool_t
       Int_t n       = 0;
       for(Int_t itime = iTimeBin; itime < fin; itime++){
        // should extract baseline here!
-       if(signal[n]>5.0) phvalue[itime] = signal[n];
+       if(signal[n]>13) phvalue[itime] = signal[n];
        n++;
       }
     }
   
     // fill the last one
     if(fDetectorPreviousTrack != -1){
-      if(TMath::Mean(fTimeMax,phvalue)>0.0){
+      if(TMath::Mean(fTimeMax,phvalue)>13.0){
        withInput = 2;
        for(Int_t k = 0; k < fTimeMax; k++){
          UpdateDAQ(fDetectorPreviousTrack,0,0,k,phvalue[k],fTimeMax);
@@ -1012,7 +1012,7 @@ Int_t AliTRDCalibraFillHisto::ProcessEventDAQ(AliTRDRawStream *rawStream, Bool_t
 
       Int_t idetector = rawStream->GetDet();                            //  current detector
       if((fDetectorPreviousTrack != idetector) && (fDetectorPreviousTrack != -1)){
-       if(TMath::Mean(nbtimebin,phvalue)>0.0){
+       if(TMath::Mean(nbtimebin,phvalue)>13.0){
          withInput = 2;
          for(Int_t k = 0; k < nbtimebin; k++){
            UpdateDAQ(fDetectorPreviousTrack,0,0,k,phvalue[k],nbtimebin);
@@ -1033,14 +1033,14 @@ Int_t AliTRDCalibraFillHisto::ProcessEventDAQ(AliTRDRawStream *rawStream, Bool_t
       Int_t n       = 0;
       for(Int_t itime = iTimeBin; itime < fin; itime++){
        // should extract baseline here!
-       if(signal[n]>5.0) phvalue[itime] = signal[n];
+       if(signal[n]>13) phvalue[itime] = signal[n];
        n++;
       }
     }
     
     // fill the last one
     if(fDetectorPreviousTrack != -1){
-      if(TMath::Mean(nbtimebin,phvalue)>0.0){
+      if(TMath::Mean(nbtimebin,phvalue)>13.0){
        withInput = 2;
        for(Int_t k = 0; k < nbtimebin; k++){
          UpdateDAQ(fDetectorPreviousTrack,0,0,k,phvalue[k],nbtimebin);

diff --git a/TRD/TRDda.cxx b/TRD/TRDda.cxx
index 51924c9c8ef1576c0b69b5a6f43ad74288dac430..7d5f8dbee1d10371395701be9e6669725a6af0d8 100644 (file)
--- a/TRD/TRDda.cxx
+++ b/TRD/TRDda.cxx
@@ -142,7 +142,7 @@ int main(int argc, char **argv) {
       continue;
     }
 
-    printf(" physic event number %d will be processed\n",(Int_t) nevents_physics);  
+    printf(" event number %d (physic event number %d) will be processed\n",(Int_t) nevents_total,(Int_t) nevents_physics);  
 
 
     /* use event - here, just write event id to result file */
@@ -154,7 +154,7 @@ int main(int argc, char **argv) {
       printf("pad status calibration\n");
       AliRawReader *rawReader = new AliRawReaderDate((void*)event);
       AliTRDRawStream *trdRawStream = new AliTRDRawStream((AliRawReader *) rawReader);
-      if(!calibpad.ProcessEvent(trdRawStream),(Bool_t)nevents_total)) passpadstatus = kFALSE;
+      if(!calibpad.ProcessEvent(trdRawStream,(Bool_t)nevents_total)) passpadstatus = kFALSE;
       
       delete trdRawStream;
       delete rawReader;
@@ -170,14 +170,17 @@ int main(int argc, char **argv) {
       Int_t result = calibra->ProcessEventDAQ(trdRawStream,(Bool_t)nevents_physics);
       if(!result) passvdrift = kFALSE;
       else nbvdrift += (Int_t) result/2;
-     
-          
-      nevents_physics++;
-
+             
+      
       delete trdRawStream;
       delete rawReader;
     } 
    
+    if(eventT==PHYSICS_EVENT){
+ 
+      nevents_physics++;
+     
+    }
     
     nevents_total++;
 
@@ -213,7 +216,7 @@ int main(int argc, char **argv) {
   if((nbvdrift > 0) && passvdrift){
     Double_t *stat = calibra->StatH((TH2 *)(calibra->GetPH2d()),1);
     // write only of enough statistics
-    if(stat[6] < 0.20) {
+    if((stat[6] < 0.20) && (stat[5] > 3000.0)) {
       calibra->GetPH2d()->Write();
       passwrite = kTRUE;
     }