Changeset 2474 for trunk/MagicSoft/Mars/macros

Timestamp:

11/05/03 14:47:52 (21 years ago)

Author:

marcos

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/macros/weights.C (modified) (4 diffs)

trunk/MagicSoft/Mars/macros/weights.C

@@ -23 +23 @@
 \* ======================================================================== */
 
+
 //////////////////////////////////////////////////////////////////////////////
 //                                                                          //
@@ -31 +32 @@
 //////////////////////////////////////////////////////////////////////////////
 
+
+// --------------------------------------------------------------------------
+//
+//
+Double_t myfunction(Double_t *x, Double_t *par)
+{
+  Double_t xx = x[0];
+
+  return pow(xx,-2)*exp(-xx/20);  
+}
+
+
+
+// --------------------------------------------------------------------------
+//
+//
 void weights(TString filename="/up1/data/Magic-MC/CameraAll/Gammas/zbin0/Gamma_zbin0_0_7_1000to1009_w0-4:4:2.root")
 {
@@ -63 +80 @@
     reader.EnableBranch("fImpact");
 
-    // ------------------
+
+    // -------------------------------------------------------------
     //
     // Option 1. Just change the slope of the MC power law spectrum
     //
-    //MMcWeightEnergySpecCalc wcalc(-2.0);
+    //MMcWeightEnergySpecCalc wcalc(-2.0);                //<-- Uncomment me
+
     //
-    // Option 2. A completely differente specturm
+    // Option 2. A completely differente specturm pass as a TF1 function
     //           e.g. spectrum with exponential cutoff
     //
-    TF1* spec = new TF1("spectrum","pow(x,[0])*exp(-x/[1])");
-    spec->SetParameter(0,-2.0); // Spectral index
-    spec->SetParameter(1,50); // Spectral index
-    MMcWeightEnergySpecCalc wcalc(spec);
+    //TF1 spec("spectrum","pow(x,[0])*exp(-x/[1])");      //<-- Uncomment me
+    //spec->SetParameter(0,-2.0);                         //<-- Uncomment me
+    //spec->SetParameter(1,50);                           //<-- Uncomment me
+    //MMcWeightEnergySpecCalc wcalc(spec);                //<-- Uncomment me
+ 
+    //
+    // Option 3. A completely differente specturm pass as a cahr*
+    //           
+    //char* func = "pow(x,-2)";                           //<-- Uncomment me
+    //MMcWeightEnergySpecCalc wcalc(func);                //<-- Uncomment me
+
+    //
+    // Option 4. A completely differente specturm pass as a c++ function
+    //     
+    MMcWeightEnergySpecCalc wcalc((void*)myfunction);   //<-- Uncomment me
+    //
+    //-------------------------------------------------------------
+
 
     MFillH hfill(&h1,"MMcEvt");
     MFillH hfill2(&h2,"MMcEvt");
     hfill2.SetWeight("MWeight");
-    //
-    //------------------
 
     tasklist.AddToList(&reader);
@@ -115 +146 @@
     hist1->DrawClone();
     hist2->DrawClone("same");    
+
 }