Changeset 4387 for trunk/MagicSoft/Mars/mtemp/mifae/macros

Timestamp:

07/15/04 11:23:59 (20 years ago)

Author:

jlopez

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mtemp/mifae/macros/lightcurve.C (modified) (10 diffs)

trunk/MagicSoft/Mars/mtemp/mifae/macros/lightcurve.C

@@ -1 @@
-Bool_t debug = kFALSE;
-
-
 Double_t ChiSquareNDof(TH1D *h1, TH1D *h2)
 {
@@ -61 +58 @@
 {
 
+  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 +118 @@
   //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 +190 @@
   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 +532 @@
                           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 +622 @@
   
   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 +885 @@
 //     }
   
+
+  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 +927 @@
 
   
-  TString psname = "./prueba.ps";
+
   TString openpsname = psname + "(";
   TString closepsname = psname + ")";
@@ -917 +933 @@
   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 +939 @@
   lightcurvegraph->SetMarkerSize(0.03);
   lightcurvegraph->Draw("AP");
+  lightcurvegraph->GetXaxis()->SetTimeDisplay(1);
   c1->Print(openpsname);