Changeset 6486 for trunk/MagicSoft/Mars/mtemp/mifae/macros/computeCoeff.C

Timestamp:

02/15/05 15:36:45 (20 years ago)

Author:

rico

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mtemp/mifae/macros/computeCoeff.C (modified) (3 diffs)

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

@@ -10 +10 @@
   /* now it's just a cross check with the same MC */
   /************************************************/
-  const Int_t ebins = 10;              // number of bins to build spectrum
+  const Int_t ebins = 20;              // number of bins to build spectrum
   const Double_t emin=10;
-  const Double_t emax=600;
-  const Double_t logemin = TMath::Log10(emin);
-  const Double_t logemax = TMath::Log10(emax);
+  const Double_t emax=200;
+  const Int_t tbins = 2;
+  const Double_t thetab[tbins+1] = {0,10,20};
   calc.SetEbins(ebins);
-  //  calc.SetEsubbins(1);
-
+  calc.SetEmin(emin);
+  calc.SetEmax(emax);
+  calc.SetThetaBinning(tbins,thetab);
+  
   // define the funtion of the desired spectrum
-  calc.SetFunction("4.e9*pow(x,-2.6+1)",emin,emax);
+  calc.SetFunction("4.e9*pow(x,-2.6+1)");
   calc.ComputeAllFactors();
   
@@ -31 +33 @@
   ccut->SetAlias("theta","MMcEvt.fTelescopeTheta*180./3.14159");
 
+
+  const Double_t logemin = TMath::Log10(emin);
+  const Double_t logemax = TMath::Log10(emax);
   TH1D* hspec = new TH1D("hspec","Spectrum",ebins,logemin,logemax);
   hspec->Sumw2();
@@ -42 +47 @@
       Float_t corrval;
       if(effa)
-        corrval = uncorrval*unfold/effa;
+        corrval = uncorrval*unfold/effa*1e9;
       else
         corrval = 0;