Index: trunk/MagicSoft/Mars/mtemp/mifae/macros/lightcurve.C
===================================================================
--- trunk/MagicSoft/Mars/mtemp/mifae/macros/lightcurve.C	(revision 4386)
+++ trunk/MagicSoft/Mars/mtemp/mifae/macros/lightcurve.C	(revision 4387)
@@ -1,5 +1,2 @@
-Bool_t debug = kFALSE;
-
-
 Double_t ChiSquareNDof(TH1D *h1, TH1D *h2)
 {
@@ -61,9 +58,16 @@
 {
 
+  gROOT->Reset();
+  gStyle->SetTimeOffset(-3600);
+
+  // Configuration
+  const Bool_t debug = kFALSE;
+  const Double_t timebin = 480.; //[sec]
+  TString psname = "./20040422_Mrk421." + timebin + "s.ps";
+  cout << psname << endl;
+  return;
+
   //Constanst
-  const Double_t ConvDegToRad = TMath::ACos(-1)/180.;
-
-  // Time binning
-  const Double_t timebin = 500.; //[sec]
+  const Double_t kConvDegToRad = TMath::ACos(-1)/180.;
   const Double_t kSec = 1e6;  //[sec/microsec]
   const Double_t kDay = 24.*60.*60.;  //[Day/sec]
@@ -114,5 +118,5 @@
   //alpha plot integration for gammas
   Double_t sigexccmin = 0.;
-  Double_t sigexccmax = 30.;
+  Double_t sigexccmax = 20.;
   Double_t bkgnormmin = 40.;
   Double_t bkgnormmax = 80.;
@@ -186,5 +190,5 @@
   plist_on.AddToList(&hAlpha_on_abs);
   MBinning binsAlpha_on_abs("BinningAlpha_on_abs");
-  Int_t nbins_abs = 9;
+  Int_t nbins_abs = 18;
   Double_t minalpha_abs = 0.;
   Double_t maxalpha_abs =90.;
@@ -528,5 +532,5 @@
 			  hAlpha_on_abs_timebin->Fill(TMath::Abs(alpha));
 			  hSrcPos_on_timebin->Fill(srcposx,srcposy);
-			  hCosZenith_on_timebin->Fill(TMath::Cos(zenith*ConvDegToRad));
+			  hCosZenith_on_timebin->Fill(TMath::Cos(zenith*kConvDegToRad));
 			  
 			  if (zenith != 0 && zenith < zenithMinimumOn[numberTimeBins-1])
@@ -618,5 +622,5 @@
   
   for (UInt_t bin=0; bin<numberTimeBins; bin++)
-    cout << bin << " timeOn " << timeOn[bin] << " min-max zenithOn " << zenithMinimumOn[bin] << "-" <<  zenithMaximumOn[bin] << " min-max cos(zenithOn) " << TMath::Cos(zenithMinimumOn[bin]*ConvDegToRad) << "-" <<  TMath::Cos(zenithMaximumOn[bin]*ConvDegToRad) << " numberOnEvents " << numberOnEvents[bin] << endl;
+    cout << bin << " timeOn " << timeOn[bin] << " min-max zenithOn " << zenithMinimumOn[bin] << "-" <<  zenithMaximumOn[bin] << " min-max cos(zenithOn) " << TMath::Cos(zenithMinimumOn[bin]*kConvDegToRad) << "-" <<  TMath::Cos(zenithMaximumOn[bin]*kConvDegToRad) << " numberOnEvents " << numberOnEvents[bin] << endl;
 
 
@@ -881,6 +885,18 @@
 //     }
   
+
+  TArrayD meanTimeBinOnInSec(numberTimeBins);
+  TArrayD widthTimeBinOnInSec(numberTimeBins);
+  
   for (UInt_t bin=0; bin<numberTimeBins; bin++)
     {
+      
+      meanTimeBinOnInSec[bin]  = (meanTimeBinOn[bin]-(Int_t)meanTimeBinOn[bin]);
+      if (meanTimeBinOnInSec[bin] > 0.5)
+	meanTimeBinOnInSec[bin] = meanTimeBinOnInSec[bin] - 1 ;
+      meanTimeBinOnInSec[bin] *= kDay;
+
+      widthTimeBinOnInSec[bin] = widthTimeBinOn[bin]*kDay;
+
       numberOffEvents[bin] = (Double_t) ((TH1F*)alphaHistoOff[bin])->Integral((Int_t)sigexccmin*nbins_abs/90+1,(Int_t)sigexccmax*nbins_abs/90+1);
       numberBkgEventsToNormOff[bin] =  (Double_t) ((TH1F*)alphaHistoOff[bin])->Integral((Int_t)bkgnormmin*nbins_abs/90+1,(Int_t)bkgnormmax*nbins_abs/90+1);
@@ -911,5 +927,5 @@
 
   
-  TString psname = "./prueba.ps";
+
   TString openpsname = psname + "(";
   TString closepsname = psname + ")";
@@ -917,5 +933,5 @@
   TCanvas *c1 = new TCanvas;
   c1->cd(1);
-  TGraphErrors* lightcurvegraph = new TGraphErrors(numberTimeBins,meanTimeBinOn.GetArray(),numberExcessEventsPerMin.GetArray(),widthTimeBinOn.GetArray(),errorNumberExcessEventsPerMin.GetArray());
+  TGraphErrors* lightcurvegraph = new TGraphErrors(numberTimeBins,meanTimeBinOnInSec.GetArray(),numberExcessEventsPerMin.GetArray(),widthTimeBinOnInSec.GetArray(),errorNumberExcessEventsPerMin.GetArray());
   lightcurvegraph->SetTitle("LightCurve");
   lightcurvegraph->SetMinimum(0.);
@@ -923,4 +939,5 @@
   lightcurvegraph->SetMarkerSize(0.03);
   lightcurvegraph->Draw("AP");
+  lightcurvegraph->GetXaxis()->SetTimeDisplay(1);
   c1->Print(openpsname);