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

Timestamp:

05/21/04 17:28:02 (21 years ago)

Author:

rico

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mtemp/mifae/macros/alpha.C (modified) (4 diffs)

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

@@ -13 +13 @@
 {    
     // constants
-    const Int_t nbins = 18;
+    const Int_t nbins = 9;
     TString psname("alphaMrk.gif");
     // general settings
@@ -36 +36 @@
     
     // ON data
-    ton->Add("/mnt/users/jrico/magic/mars/mars/mtemp/mifae/programs/srcPosPrueba.root");
+    ton->Add("/mnt/users/jrico/magic/mars/mars/mtemp/mifae/hillas/crab20040215OnRotateCalA-D.root");
 
     TChain* toff= new TChain("Parameters");
 
     // OFF data
-    toff->Add("/mnt/users/jrico/magic/mars/mars/mtemp/mifae/programs/srcPosOffPrueba.root");
+    toff->Add("/mnt/users/jrico/magic/mars/mars/mtemp/mifae/hillas/crab20040215OffRotateCalA-H.root");
     
     TCanvas *c1 = new TCanvas("c1","c1",800,500);
@@ -54 +54 @@
     ton->SetAlias("event","MRawEvtHeader.fDAQEvtNumber");
     ton->SetAlias("alpha","abs(MHillasSrc.fAlpha)");
+    ton->SetAlias("leakage","MNewImagePar.fInnerLeakage1");
     
     toff->SetAlias("length","MHillas.fLength*0.6/189");
@@ -62 +63 @@
     toff->SetAlias("event","MRawEvtHeader.fDAQEvtNumber");
     toff->SetAlias("alpha","abs(MHillasSrc.fAlpha)");
-        
+    toff->SetAlias("leakage","MNewImagePar.fInnerLeakage1");
+
     // define  cut(s)
-//    const char* eventcut="event%3==0 &&";
+//    const char* eventcut="event%4==0 &&";
     const char* eventcut="";
     const char varcut[512];
     //    sprintf(varcut,"size>%f && size<%f && length>%f  &&length<%f && width>%f &&  width<%f",cutsize,upcutsize,lengthlowcut,lengthcut,widthlowcut,widthcut);
-    sprintf(varcut,"size>%f && size<%f  && dist>0.2 && dist<1.1 && length>%f  &&length<%f && width>%f &&  width<%f",cutsize,upcutsize,lengthlowcut,lengthcut,widthlowcut,widthcut);
+    sprintf(varcut,"size>%f && size<%f  && dist>0.2 && dist<1.1 && leakage==0 && length>%f  &&length<%f && width>%f &&  width<%f",cutsize,upcutsize,lengthlowcut,lengthcut,widthlowcut,widthcut);
 //    sprintf(varcut,"size>%f && dist>0.2 && dist<0.8 && length<0.21 && width<0.125",cutsize,upcutsize,lengthcut,widthcut);
 //    sprintf(varcut,"size>%f && size<%f  && dist>0.2 && dist<0.8 && length<0.21 && width<0.125",cutsize,upcutsize);