Ignore:
Timestamp:
08/12/04 06:58:34 (20 years ago)
Author:
wittek
Message:
*** empty log message ***
File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/MagicSoft/Mars/macros/ONOFFAnalysis.C

    r3598 r4584  
    171171      //-----------------------------------------------
    172172      //TString tag = "040126";
    173       TString tag = "040127";
    174       //TString tag = "040215";
     173      //TString tag = "040127";
     174      TString tag = "040215";
    175175
    176176      //const char *offfile = "~magican/ct1test/wittek/offdata.preproc";
    … …  
    183183      //const char *offfile  = "*.OFF";
    184184      // 15 Feb 04
    185       const char *offfile  = "*.OFF";
     185      //const char *offfile  = "*.OFF";
     186      const char *offfile  = "174*.OFF";
    186187
    187188
    … …  
    191192      //const char *onfile  = "MCerPhot_output";
    192193      //const char *onfile  = "*.ON";
    193       const char *onfile  = "1216*.ON";
     194      //const char *onfile  = "1216*.ON";
    194195      // 27 Jan 04
    195196      //const char *onfile  = "12*.ON";
    … …  
    198199      // 15 Feb 04
    199200      //const char *onfile  = "*.ON";
    200 
    201 
    202      const char *mcfile  = "/data/MAGIC/mc_eth/magLQE_3/gh/0/0/G_M0_00_0_550*.root";
     201      const char *onfile  = "174*.ON";
     202
     203
     204      //const char *mcfile  = "/data/MAGIC/mc_eth/magLQE_3/gh/0/0/G_M0_00_0_550*.root";
    203205      //const char *mcfile  = "/data/MAGIC/mc_eth/magLQE_4/gh/0/0/G_M0_00_0_550*.root";
    204206      //const char *mcfile  = "/data/MAGIC/mc_eth/magLQE_5/gh/0/0/G_M0_00_0_550*.root";
     207      //const char *mcfile = "calibrated_gammas";
     208      const char *mcfile = "calibrated_data_david*";
     209
    205210      //-----------------------------------------------
    206211
    … …  
    224229
    225230     // 15 Feb 04, David
    226      if (tag == "040215")
    227        TString inPath = "/mnt/magicserv01/scratch/David/CalibratedData/Crab/2004_02_15/";
    228    
     231     //if (tag == "040215")
     232     //  TString inPath = "/mnt/magicserv01/scratch/David/CalibratedData/Crab/2004_02_15/";
     233     if (tag == "040215")   
     234       TString inPathON = "/.magic/magicserv01/scratch/Daniel/CalibratedData/MispointingTest/2004_02_15ForgottenData/finesam20/";
     235       TString inPathOFF = "/.magic/magicserv01/scratch/David/CalibratedData/Mkn421/2004_02_15/";
     236       TString inPathMC = "/.magic/magicserv01/scratch/David/MCData/MCApril2004/Data/gammas_highnoise/";
    229237
    230238      // path for output from Mars
    231239      //TString outPath = "~wittek/datacrab_feb04/";
    232240     //TString outPath = "~wittek/datacrab_01march04/";
    233       TString outPath = "~wittek/datacrab_19mar04/";
     241      TString outPath = "~wittek/datacrab_27april04/";
    234242
    235243      //-----------------------------------------------
    … …  
    248256
    249257    Bool_t JobA    = kTRUE; 
    250     Bool_t GPad    = kFALSE;    // generate padding histograms?
    251     Bool_t WPad    = kFALSE;   // write out padding histograms ?
    252     Bool_t RPad    = kFALSE;   // read in padding histograms ?
    253     Bool_t Wout    = kTRUE;   // write out root file ON1.root
    254                                // (or OFF1.root or MC1.root)?
     258    Bool_t GPadON  = kFALSE;    // \  generate padding histograms
     259    Bool_t GPadOFF = kFALSE;    //  | and write them onto disk
     260    Bool_t GPadMC  = kFALSE;    // /
     261    Bool_t Merge   = kFALSE;   // merge padding histograms
     262                               // and write them onto disk
     263    Bool_t Wout    = kTRUE;   // read in merged padding histograms and
     264                              // write out root file of padded data
     265                              // (ON1.root, OFF1.root or MC1.root)
    255266
    256267
    … …  
    310321    // Job E_XX : extended version of E_XX (including flux plots) 
    311322    //  - select g/h separation method XX
    312     //  - read MC root file 
    313     //  - calculate eff. collection area
     323    //  - read MC root file
     324    //  - calculate eff. collection area
    314325    //  - optimize energy estimator
    315326    //  - read ON root file
    … …  
    346357    gLog << "Macro ONOFFAnalysis : Start of Job A" << endl;
    347358    gLog << "" << endl;
    348     gLog << "Macro ONOFFAnalysis : JobA, WPad, RPad, Wout = "
    349          << (JobA ? "kTRUE" : "kFALSE")  << ",  "
    350          << (WPad ? "kTRUE" : "kFALSE")  << ",  "
    351          << (RPad ? "kTRUE" : "kFALSE")  << ",  "
    352          << (Wout ? "kTRUE" : "kFALSE")  << endl;
     359    gLog << "Macro ONOFFAnalysis : JobA, GPadON, GPadOFF, GPadMC, Merge, Wout = "
     360         << (JobA    ? "kTRUE" : "kFALSE")  << ",  "
     361         << (GPadON  ? "kTRUE" : "kFALSE")  << ",  "
     362         << (GPadOFF ? "kTRUE" : "kFALSE")  << ",  "
     363         << (GPadMC  ? "kTRUE" : "kFALSE")  << ",  "
     364         << (Merge   ? "kTRUE" : "kFALSE")  << ",  "
     365         << (Wout    ? "kTRUE" : "kFALSE")  << endl;
    353366   
    354367
    … …  
    358371
    359372
    360     TString fileON  = inPath;
     373    TString fileON  = inPathON;
    361374    fileON += onfile;
    362375    fileON += ".root";
    363376
    364     TString fileOFF = inPath;
     377    TString fileOFF = inPathOFF;
    365378    fileOFF += offfile;
    366379    fileOFF += ".root";
    367380
    368     TString fileMC  = mcfile;
     381    TString fileMC = inPathMC;
    369382    fileMC += mcfile;
    370383    fileMC += ".root";
    … …  
    373386         << fileOFF << ",  " << fileMC   << endl;
    374387
    375     // name of file to conatin the histograms for the padding
     388    // name of file to conatin the merged histograms for the padding
    376389    TString outNameSigTh = outPath;
    377390    outNameSigTh += "SigmaTheta";
    … …  
    379392
    380393    //--------------------------------------------------
     394    // name of files to contain the paddding histograms of ON, OFF and MC data
     395      TString NamePadON(outPath);
     396      NamePadON += "PadON";
     397      NamePadON += ".root";
     398
     399      TString NamePadOFF(outPath);
     400      NamePadOFF += "PadOFF";
     401      NamePadOFF += ".root";
     402
     403      TString NamePadMC(outPath);
     404      NamePadMC += "PadMC";
     405      NamePadMC += ".root";
     406
     407    //--------------------------------------------------
    381408    // type of data to be padded
    382     TString typeInput = "ON";
    383     //TString typeInput = "OFF";
    384     //TString typeInput = "MC";
     409      //TString typeInput = "ON";
     410      //TString typeInput = "OFF";
     411      TString typeInput = "MC";
    385412    gLog << "typeInput = " << typeInput << endl;
    386413
    … …  
    416443    //
    417444    // read ON, OFF and MC data files
    418     // generate (or read in) the padding histograms for ON and OFF data
    419     //                       and merge these histograms
     445    // generate (or read in) the padding histograms for ON, OFF and MC data
     446    //
    420447
    421448    MPad pad;
    … …  
    423450    pad.SetDataType(typeInput);
    424451
     452    if (GPadON || GPadOFF || GPadMC)
     453    {
    425454    // generate the padding histograms
    426     if (GPad)
    427     {
    428455      gLog << "=====================================================" << endl;
    429456      gLog << "Start generating the padding histograms" << endl;
     457    }
     458
    430459      //-----------------------------------------
    431460      // ON events
     461
     462      if (GPadON)
     463      {
    432464
    433465      gLog << "-----------" << endl;
    … …  
    438470      MParList pliston;
    439471
     472    MObservatory observon;
     473
    440474    MReadMarsFile  readON("Events", fileON);
    441475    readON.DisableAutoScheme();
    442     //MCT1ReadPreProc readON(fileON);
    443 
    444       //MFSelBasic selthetaon;
    445       //selthetaon.SetCuts(-100.0, 29.5, 35.5);
    446       //MContinue contthetaon(&selthetaon);
    447 
    448       MBlindPixelCalc blindon;
    449       blindon.SetUseBlindPixels();
     476
     477    MGeomApply        apply;
     478
     479    MSourcePosfromStarPos sourcefromstaron;
     480    sourcefromstaron.AddFile("~wittek/datacrab_26feb04/positions040215.txt", 0);
     481
     482      MBlindPixelsCalc2 blindon;
     483      //blindon.SetUseBlindPixels();
     484      blindon.SetCheckPedestalRms();
    450485
    451486      MFSelBasic selbasicon;
    … …  
    459494
    460495      MHSigmaTheta sigthON("SigmaThetaON");
    461       MFillH fillsigthetaON ("SigmaThetaON[MHSigmaTheta]", "MMcEvt");
     496      MFillH fillsigthetaON ("SigmaThetaON[MHSigmaTheta]", "MPointingPos");
    462497      fillsigthetaON.SetName("FillSigTheta");   
    463498 
    … …  
    467502      pliston.AddToList(&tliston);
    468503      InitBinnings(&pliston);
     504      pliston.AddToList(&observon);
    469505      pliston.AddToList(&blindON);
    470506      pliston.AddToList(&sigthON);
    … …  
    475511   
    476512      tliston.AddToList(&readON);
    477       //tliston.AddToList(&contthetaon);
     513      tliston.AddToList(&apply);
     514      tliston.AddToList(&sourcefromstaron);
    478515
    479516      tliston.AddToList(&blindon);
    … …  
    489526      evtloopon.SetProgressBar(&baron);
    490527
    491       Int_t maxeventson = -1;
    492       //Int_t maxeventson = 10000;
     528      //Int_t maxeventson = -1;
     529      Int_t maxeventson = 10000;
    493530      if ( !evtloopon.Eventloop(maxeventson) )
    494531          return;
    … …  
    507544      TH2D *blindnthon  = blindON.GetBlindN();
    508545
     546      gLog << "" << endl;
     547      gLog << "Macro ONOFFAnalysis : write padding plots for ON data from file "
     548           << NamePadON << endl;
     549
     550      TFile filewriteon(NamePadON, "RECREATE", "");
     551      sigthon->Write();
     552      sigpixthon->Write();
     553      diffpixthon->Write();
     554      blindidthon->Write();
     555      blindnthon->Write();
     556      filewriteon.Close();
     557
     558      gLog << "" << endl;
     559      gLog << "Macro ONOFFAnalysis : padding plots for ON data written onto file "
     560           << NamePadON << endl;
     561      }
     562      // end of GPadON
     563
     564
    509565      //-----------------------------------------
    510566      // OFF events
     567
     568      if (GPadOFF)
     569      {
    511570
    512571      gLog << "------------" << endl;
    … …  
    517576      MParList plistoff;
    518577
     578    MObservatory observoff;
     579
    519580    MReadMarsFile  readOFF("Events", fileOFF);
    520581    readOFF.DisableAutoScheme();
    521     //      MCT1ReadPreProc readOFF(fileOFF);
    522 
    523       MFSelBasic selthetaoff;
    524       selthetaoff.SetCuts(-100.0, 29.5, 35.5);
    525       MContinue contthetaoff(&selthetaoff);
    526 
    527       MBlindPixelCalc blindoff;
    528       blindoff.SetUseBlindPixels();
     582
     583    MSourcePosfromStarPos sourcefromstaroff;
     584    sourcefromstaroff.AddFile("~wittek/datacrab_26feb04/positions040215.txt", 0);
     585
     586      MBlindPixelsCalc2 blindoff;
     587      //blindoff.SetUseBlindPixels();
     588      blindoff.SetCheckPedestalRms();
    529589
    530590      MFSelBasic selbasicoff;
    … …  
    538598
    539599      MHSigmaTheta sigthOFF("SigmaThetaOFF");
    540       MFillH fillsigthetaOFF ("SigmaThetaOFF[MHSigmaTheta]", "MMcEvt");
     600      MFillH fillsigthetaOFF ("SigmaThetaOFF[MHSigmaTheta]", "MPointingPos");
    541601      fillsigthetaOFF.SetName("FillSigThetaOFF");     
    542602
    … …  
    546606      plistoff.AddToList(&tlistoff);
    547607      InitBinnings(&plistoff);
     608      plistoff.AddToList(&observoff);
    548609      plistoff.AddToList(&blindOFF);
    549610      plistoff.AddToList(&sigthOFF);
    … …  
    554615   
    555616      tlistoff.AddToList(&readOFF);
    556       //tlistoff.AddToList(&contthetaoff);
     617      tlistoff.AddToList(&sourcefromstaroff);
    557618
    558619      tlistoff.AddToList(&blindoff);
    … …  
    568629      evtloopoff.SetProgressBar(&baroff);
    569630
    570       Int_t maxeventsoff = -1;
    571       //Int_t maxeventsoff = 20000;
     631      //Int_t maxeventsoff = -1;
     632      Int_t maxeventsoff = 10000;
    572633      if ( !evtloopoff.Eventloop(maxeventsoff) )
    573634          return;
    … …  
    586647      TH2D *blindnthoff  = blindOFF.GetBlindN();
    587648
     649      gLog << "" << endl;
     650      gLog << "Macro ONOFFAnalysis : write padding plots for OFF data from file "
     651           << NamePadOFF << endl;
     652
     653      TFile filewriteoff(NamePadOFF, "RECREATE", "");
     654      sigthoff->Write();
     655      sigpixthoff->Write();
     656      diffpixthoff->Write();
     657      blindidthoff->Write();
     658      blindnthoff->Write();
     659      filewriteoff.Close();
     660
     661      gLog << "" << endl;
     662      gLog << "Macro ONOFFAnalysis : padding plots for OFF data written onto file "
     663           << NamePadOFF << endl;
     664      }
     665      // end of GPadOFF
     666
     667
    588668
    589669      //-----------------------------------------
    590670      // MC events
     671
     672      if (GPadMC)
     673      {
    591674
    592675      gLog << "------------" << endl;
    … …  
    597680      MParList plistmc;
    598681
     682    MObservatory observmc;
     683
    599684    MReadMarsFile  readMC("Events", fileMC);
    600685    readMC.DisableAutoScheme();
    601     //      MCT1ReadPreProc readMC(fileMC);
    602 
    603       MFSelBasic selthetamc;
    604       selthetamc.SetCuts(-100.0, 29.5, 35.5);
    605       MContinue contthetamc(&selthetamc);
    606 
    607       MBlindPixelCalc blindmc;
    608       blindmc.SetUseBlindPixels();
     686
     687    MSourcePosfromStarPos sourcefromstarmc;
     688
     689      MBlindPixelsCalc2 blindmc;
     690      //blindmc.SetUseBlindPixels();
     691      blindmc.SetCheckPedestalRms();
    609692
    610693      MFSelBasic selbasicmc;
    … …  
    618701
    619702      MHSigmaTheta sigthMC("SigmaThetaMC");
    620       MFillH fillsigthetaMC ("SigmaThetaMC[MHSigmaTheta]", "MMcEvt");
     703      MFillH fillsigthetaMC ("SigmaThetaMC[MHSigmaTheta]", "MPointingPos");
    621704      fillsigthetaMC.SetName("FillSigThetaMC");     
    622705
    … …  
    626709      plistmc.AddToList(&tlistmc);
    627710      InitBinnings(&plistmc);
     711      plistmc.AddToList(&observmc);
    628712      plistmc.AddToList(&blindMC);
    629713      plistmc.AddToList(&sigthMC);
    … …  
    634718   
    635719      tlistmc.AddToList(&readMC);
    636       //tlistmc.AddToList(&contthetamc);
     720      tlistmc.AddToList(&sourcefromstarmc);
    637721
    638722      tlistmc.AddToList(&blindmc);
    … …  
    648732      evtloopmc.SetProgressBar(&barmc);
    649733
    650       Int_t maxeventsmc = -1;
    651       //Int_t maxeventsmc = 20000;
     734      //Int_t maxeventsmc = -1;
     735      Int_t maxeventsmc = 10000;
    652736      if ( !evtloopmc.Eventloop(maxeventsmc) )
    653737          return;
    … …  
    666750      TH2D *blindnthmc  = blindMC.GetBlindN();
    667751
     752      gLog << "" << endl;
     753      gLog << "Macro ONOFFAnalysis : write padding plots for MC data from file "
     754           << NamePadMC << endl;
     755
     756      TFile filewritemc(NamePadMC, "RECREATE", "");
     757      sigthmc->Write();
     758      sigpixthmc->Write();
     759      diffpixthmc->Write();
     760      blindidthmc->Write();
     761      blindnthmc->Write();
     762      filewritemc.Close();
     763
     764      gLog << "" << endl;
     765      gLog << "Macro ONOFFAnalysis : padding plots for MC data written onto file "
     766           << NamePadMC << endl;
     767      }
     768      // end of GPadMC
    668769
    669770      //-----------------------------------------
    670771
     772    if (GPadON || GPadOFF || GPadMC)
     773    {
     774      gLog << "" << endl;
    671775      gLog << "End of generating the padding histograms" << endl;
    672776      gLog << "=====================================================" << endl;
    673 
    674       pad.MergeONOFFMC(sigthmc,      diffpixthmc,  sigpixthmc,
    675                        blindidthmc,  blindnthmc,
    676                        sigthon,      diffpixthon,  sigpixthon,
    677                        blindidthon,  blindnthon,
    678                        sigthoff,     diffpixthoff, sigpixthoff,
    679                        blindidthoff, blindnthoff);
    680 
    681 
    682       if (WPad)
    683       {
    684         // write the padding histograms onto a file  ---------
    685         pad.WritePaddingDist(outNameSigTh);     
    686       }
    687777    }
    688778
    689     // read the padding histograms ---------------------------
    690     if (RPad)
     779    if (Merge)
    691780    {
    692       pad.ReadPaddingDist(outNameSigTh);
     781
     782      //-----------------------------------
     783      TH2D sigon;
     784      TH3D sigpixon;
     785      TH3D diffpixon;
     786      TH2D blindidon;
     787      TH2D blindnon;
     788
     789      TFile filereadon(NamePadON);
     790      filereadon.ls();
     791      sigon.Read("2D-ThetaSigmabar(Inner)");
     792      sigpixon.Read("3D-ThetaPixSigma");
     793      diffpixon.Read("3D-ThetaPixDiff");
     794      blindidon.Read("2D-IdBlindPixels");
     795      blindnon.Read("2D-NBlindPixels");
     796
     797      TH2D *sigthon     = new TH2D( (const TH2D&)sigon     );
     798      TH3D *sigpixthon  = new TH3D( (const TH3D&)sigpixon  );
     799      TH3D *diffpixthon = new TH3D( (const TH3D&)diffpixon );
     800      TH2D *blindidthon = new TH2D( (const TH2D&)blindidon );
     801      TH2D *blindnthon  = new TH2D( (const TH2D&)blindnon  );
     802
     803      /*
     804      TH2D *sigthon     = (TH2D*)sigon.Clone();
     805      TH3D *sigpixthon  = (TH3D*)sigpixon.Clone();
     806      TH3D *diffpixthon = (TH3D*)diffpixon.Clone();
     807      TH2D *blindidthon = (TH2D*)blindidon.Clone();
     808      TH2D *blindnthon  = (TH2D*)blindnon.Clone();
     809      */
     810
     811      //filereadon.Close();
     812      gLog << "" << endl;
     813      gLog << "Macro ONOFFAnalysis : padding plots for ON data read from file "
     814           << NamePadON << endl;
     815
     816
     817      //-----------------------------------
     818      TH2D sigoff;
     819      TH3D sigpixoff;
     820      TH3D diffpixoff;
     821      TH2D blindidoff;
     822      TH2D blindnoff;
     823
     824      TFile filereadoff(NamePadOFF);
     825      filereadoff.ls();
     826      sigoff.Read("2D-ThetaSigmabar(Inner)");
     827      sigpixoff.Read("3D-ThetaPixSigma");
     828      diffpixoff.Read("3D-ThetaPixDiff");
     829      blindidoff.Read("2D-IdBlindPixels");
     830      blindnoff.Read("2D-NBlindPixels");
     831 
     832      TH2D *sigthoff     = new TH2D( (const TH2D&)sigoff     );
     833      TH3D *sigpixthoff  = new TH3D( (const TH3D&)sigpixoff  );
     834      TH3D *diffpixthoff = new TH3D( (const TH3D&)diffpixoff );
     835      TH2D *blindidthoff = new TH2D( (const TH2D&)blindidoff );
     836      TH2D *blindnthoff  = new TH2D( (const TH2D&)blindnoff  );
     837
     838      /*
     839      TH2D *sigthoff     = (TH2D*)sigoff.Clone();
     840      TH3D *sigpixthoff  = (TH3D*)sigpixoff.Clone();
     841      TH3D *diffpixthoff = (TH3D*)diffpixoff.Clone();
     842      TH2D *blindidthoff = (TH2D*)blindidoff.Clone();
     843      TH2D *blindnthoff  = (TH2D*)blindnoff.Clone();
     844      */
     845
     846      //filereadoff.Close();
     847      gLog << "" << endl;
     848      gLog << "Macro ONOFFAnalysis : padding plots for OFF data read from file "
     849           << NamePadOFF << endl;
     850
     851
     852      //-----------------------------------
     853      TH2D sigmc;
     854      TH3D sigpixmc;
     855      TH3D diffpixmc;
     856      TH2D blindidmc;
     857      TH2D blindnmc;
     858
     859      TFile filereadmc(NamePadMC);
     860      filereadmc.ls();
     861      sigmc.Read("2D-ThetaSigmabar(Inner)");
     862      sigpixmc.Read("3D-ThetaPixSigma");
     863      diffpixmc.Read("3D-ThetaPixDiff");
     864      blindidmc.Read("2D-IdBlindPixels");
     865      blindnmc.Read("2D-NBlindPixels");
     866 
     867      TH2D *sigthmc     = new TH2D( (const TH2D&)sigmc     );
     868      TH3D *sigpixthmc  = new TH3D( (const TH3D&)sigpixmc  );
     869      TH3D *diffpixthmc = new TH3D( (const TH3D&)diffpixmc );
     870      TH2D *blindidthmc = new TH2D( (const TH2D&)blindidmc );
     871      TH2D *blindnthmc  = new TH2D( (const TH2D&)blindnmc  );
     872
     873      /*
     874      TH2D *sigthmc     = (TH2D*)sigmc.Clone();
     875      TH3D *sigpixthmc  = (TH3D*)sigpixmc.Clone();
     876      TH3D *diffpixthmc = (TH3D*)diffpixmc.Clone();
     877      TH2D *blindidthmc = (TH2D*)blindidmc.Clone();
     878      TH2D *blindnthmc  = (TH2D*)blindnmc.Clone();
     879      */
     880
     881      //filereadmc.Close();
     882      gLog << "" << endl;
     883      gLog << "Macro ONOFFAnalysis : padding plots for MC data read from file "
     884           << NamePadMC << endl;
     885
     886      pad.MergeONOFFMC(*sigthmc,     *diffpixthmc, *sigpixthmc,
     887                       *blindidthmc, *blindnthmc,
     888                       *sigthon,     *diffpixthon, *sigpixthon,
     889                       *blindidthon, *blindnthon,
     890                       *sigthoff,    *diffpixthoff,*sigpixthoff,
     891                       *blindidthoff,*blindnthoff);
     892
     893      // write the target padding histograms onto a file  ---------
     894      pad.WritePaddingDist(outNameSigTh);     
    693895    }
     896    // end of Merge
     897
    694898
    695899
    … …  
    698902  if (Wout)
    699903  {
     904    // read the target padding histograms ---------------------------
     905    pad.ReadPaddingDist(outNameSigTh);
     906
     907
    700908    gLog << "=====================================================" << endl;
    701909    gLog << "Start the padding" << endl;
    … …  
    723931    MGeomApply        apply;
    724932
    725     //MPedestalWorkaround waround;
    726 
    727933    // a way to find out whether one is dealing with MC :
    728     MFDataMember f1("MRawRunHeader.fRunType", '>', 255.5);  // MC
    729     f1.SetName("Select MC");
    730     MFDataMember f2("MRawRunHeader.fRunType", '<', 255.5);  // data
    731     f2.SetName("Select Data");
    732 
    733     //if (typeInput ==  "ON")
    734     //{
    735     //   MCT1PointingCorrCalc pointcorr(sourceName, "MCT1PointingCorrCalc",
    736     //                                             "MCT1PointingCorrCalc");
    737     //}
     934    //MFDataMember f1("MRawRunHeader.fRunType", '>', 255.5);  // MC
     935    //f1.SetName("Select MC");
     936    //MFDataMember f2("MRawRunHeader.fRunType", '<', 255.5);  // data
     937    //f2.SetName("Select Data");
     938
    738939    MSourcePosfromStarPos sourcefromstar;
    739940    if (typeInput == "ON")
    740941    {
    741       sourcefromstar.SetSourceAndStarPosition(
    742                                "Crab", 22,  0, 52, 5, 34, 32.0,
    743                            "Zeta-Tau", 21,  8, 33, 5, 37, 38.7);
    744       sourcefromstar.AddStar("Tau114", 21, 56, 13, 5, 27, 38.1);
    745       //sourcefromstar.AddFile("~wittek/datacrab_26feb04/positionsOn.4.txt",0);
    746 
    747       if(inPath == "/.magic/magicserv01/scratch/calibrated/")
    748       {
    749         sourcefromstar.AddFile("~wittek/datacrab_26feb04/positions2stars.txt", 0);
    750         //sourcefromstar.AddFile("~wittek/datacrab_26feb04/positionsNoStar.txt", 0);
    751       }
    752       else if(inPath == "/.magic/magicserv01/scratch/calibrated26/")
    753         sourcefromstar.AddFile("~wittek/datacrab_26feb04/stars_2004_01_26", 0);
    754 
    755       else if(inPath == "/mnt/magicserv01/scratch/David/CalibratedData/Crab/2004_02_15/")
    756         sourcefromstar.AddFile("~wittek/datacrab_26feb04/positions040215.txt", 0);
     942      //sourcefromstar.SetSourceAndStarPosition(
     943      //                         "Crab", 22,  0, 52, 5, 34, 32.0,
     944      //                           "Zeta-Tau", 21,  8, 33, 5, 37, 38.7);
     945      //sourcefromstar.AddStar("Tau114", 21, 56, 13, 5, 27, 38.1);
     946      ////sourcefromstar.AddFile("~wittek/datacrab_26feb04/positionsOn.4.txt",0);
     947
     948      //if(inPath == "/.magic/magicserv01/scratch/calibrated/")
     949      //{
     950        //sourcefromstar.AddFile("~wittek/datacrab_26feb04/positions2stars.txt", 0);
     951        ////sourcefromstar.AddFile("~wittek/datacrab_26feb04/positionsNoStar.txt", 0);
     952      //}
     953      //else if(inPath == "/.magic/magicserv01/scratch/calibrated26/")
     954      //  sourcefromstar.AddFile("~wittek/datacrab_26feb04/stars_2004_01_26", 0);
     955
     956      //else if(inPath == "/mnt/magicserv01/scratch/David/CalibratedData/Crab/2004_02_15/")
     957      sourcefromstar.AddFile("~wittek/datacrab_26feb04/positions040215.txt", 0);
    757958    }
    758959    else if (typeInput == "OFF")
    759960    {
    760       if(inPath == "/.magic/magicserv01/scratch/calibrated/")
    761         sourcefromstar.AddFile("~wittek/datacrab_26feb04/positionsOff.txt", 0);
    762       else if(inPath == "/.magic/magicserv01/scratch/calibrated26/")
    763         sourcefromstar.AddFile("~wittek/datacrab_26feb04/positions_2004_01_26", 0);
    764       else if(inPath == "/mnt/magicserv01/scratch/David/CalibratedData/Crab/2004_02_15")
    765         sourcefromstar.AddFile("~wittek/datacrab_26feb04/positions040215.txt", 0);
     961      //if(inPath == "/.magic/magicserv01/scratch/calibrated/")
     962      //  sourcefromstar.AddFile("~wittek/datacrab_26feb04/positionsOff.txt", 0);
     963      //else if(inPath == "/.magic/magicserv01/scratch/calibrated26/")
     964      //  sourcefromstar.AddFile("~wittek/datacrab_26feb04/positions_2004_01_26", 0);
     965      //else if(inPath == "/mnt/magicserv01/scratch/David/CalibratedData/Crab/2004_02_15")
     966      sourcefromstar.AddFile("~wittek/datacrab_26feb04/positions040215.txt", 0);
    766967    }
    767 
    768 
    769     //MBlindPixelCalc blindbeforepad;
     968    else if (typeInput == "MC")
     969    {
     970    }
     971
     972    MBlindPixelsCalc2 blindbeforepad;
    770973    //blindbeforepad.SetUseBlindPixels();
    771     //blindbeforepad.SetName("BlindBeforePadding");
    772 
    773     //MBlindPixelCalc blind;
    774     //blind.SetUseBlindPixels();
    775     //blind.SetUseInterpolation();
    776     //blind.SetName("BlindAfterPadding");
    777 
    778     MSigmabarCalc sigbar;
     974    blindbeforepad.SetCheckPedestalRms();
     975    blindbeforepad.SetName("BlindBeforePadding");
    779976
    780977    //MBadPixelCalcRms blind;
    … …  
    794991    MSigmabarCalc sigbarcalc;
    795992
    796     MFillH fillsigtheta ("SigmaTheta[MHSigmaTheta]", "MMcEvt");
     993    MFillH fillsigtheta ("SigmaTheta[MHSigmaTheta]", "MPointingPos");
    797994    fillsigtheta.SetName("HSigmaTheta");
    798995
    … …  
    8451042      //write.AddContainer("MMcRunHeader",  "RunHeaders", kFALSE);
    8461043      //write.AddContainer("MTime",         "Events");
    847       write.AddContainer("MMcEvt",        "Events");
    848       //write.AddContainer("ThetaOrig",     "Events");
     1044      write.AddContainer("MPointingPos", "Events");
    8491045      write.AddContainer("MSrcPosCam",    "Events");
    8501046      write.AddContainer("MSigmabar",     "Events");
    … …  
    8681064   
    8691065    tliston.AddToList(&read);
    870 
    8711066    //tliston.AddToList(&f1);
    8721067    //tliston.AddToList(&f2);
    8731068    tliston.AddToList(&apply);
    874     //tliston.AddToList(&waround);
    8751069    tliston.AddToList(&sourcefromstar);
    8761070
    877     //tliston.AddToList(&blindbeforepad);
    878     //  tliston.AddToList(&pad);
    879     //  if (typeInput ==  "ON")
    880     //  tliston.AddToList(&pointcorr);
    881 
    882     tliston.AddToList(&sigbar);
     1071    tliston.AddToList(&blindbeforepad);
     1072
     1073    tliston.AddToList(&contbasic);
     1074    tliston.AddToList(&pad);
     1075
    8831076    tliston.AddToList(&blind);
    884     tliston.AddToList(&contbasic);
    885 
    8861077    tliston.AddToList(&fillblind);
    887     //tliston.AddToList(&sigbarcalc);
     1078    tliston.AddToList(&sigbarcalc);
    8881079    tliston.AddToList(&fillsigtheta);
    8891080    tliston.AddToList(&clean);
    … …  
    9131104    //  evtloop.Write();
    9141105
    915     Int_t maxevents = -1;
    916     //Int_t maxevents = 1000;
     1106    //Int_t maxevents = -1;
     1107    Int_t maxevents = 2000;
    9171108    if ( !evtloop.Eventloop(maxevents) )
    9181109        return;
    … …  
    11091300    matrixg.EnableGraphicalOutput();
    11101301
    1111     matrixg.AddColumn("cos(MMcEvt.fTelescopeTheta)");
     1302    matrixg.AddColumn("cos((MPointingPos.fZd)/kRad2Deg)");
    11121303    matrixg.AddColumn("MSigmabar.fSigmabar");
    11131304    matrixg.AddColumn("log10(MHillas.fSize)");
    … …  
    11981389    bing.SetEdges(10, 0., 1.0);
    11991390
    1200     MH3 gref("cos(MMcEvt.fTelescopeTheta)");
     1391    MH3 gref("cos((MPointingPos.fZd)/kRad2Deg)");
    12011392    gref.SetName(mgname);
    12021393    MH::SetBinning(&gref.GetHist(), &bing);
    … …  
    12311422      writetraing.AddContainer("MRawRunHeader", "RunHeaders");
    12321423      writetraing.AddContainer("MTime",         "Events");
    1233       writetraing.AddContainer("MMcEvt",        "Events");
     1424      writetraing.AddContainer("MPointingPos",  "Events");
    12341425      writetraing.AddContainer("ThetaOrig",     "Events");
    12351426      writetraing.AddContainer("MSrcPosCam",    "Events");
    … …  
    12481439      writetestg.AddContainer("MRawRunHeader", "RunHeaders");
    12491440      writetestg.AddContainer("MTime",         "Events");
    1250       writetestg.AddContainer("MMcEvt",        "Events");
     1441      writetestg.AddContainer("MPointingPos",  "Events");
    12511442      writetestg.AddContainer("ThetaOrig",     "Events");
    12521443      writetestg.AddContainer("MSrcPosCam",    "Events");
    … …  
    13311522    binh.SetEdges(10, 0., 1.0);
    13321523
    1333     //MH3 href("cos(MMcEvt.fTelescopeTheta)");
     1524    //MH3 href("cos((MPointingPos.fZd)/kRad2Deg)");
    13341525    //href.SetName(mhname);
    13351526    //MH::SetBinning(&href.GetHist(), &binh);
    … …  
    13661557      writetrainh.AddContainer("MRawRunHeader", "RunHeaders");
    13671558      writetrainh.AddContainer("MTime",         "Events");
    1368       writetrainh.AddContainer("MMcEvt",        "Events");
     1559      writetrainh.AddContainer("MPointingPos",  "Events");
    13691560      writetrainh.AddContainer("ThetaOrig",     "Events");
    13701561      writetrainh.AddContainer("MSrcPosCam",    "Events");
    … …  
    13821573      writetesth.AddContainer("MRawRunHeader", "RunHeaders");
    13831574      writetesth.AddContainer("MTime",         "Events");
    1384       writetesth.AddContainer("MMcEvt",        "Events");
     1575      writetesth.AddContainer("MPointingPos",  "Events");
    13851576      writetesth.AddContainer("ThetaOrig",     "Events");
    13861577      writetesth.AddContainer("MSrcPosCam",    "Events");
    … …  
    16671858      write.AddContainer("MRawRunHeader", "RunHeaders");
    16681859      write.AddContainer("MTime",         "Events");
    1669       write.AddContainer("MMcEvt",        "Events");
     1860      write.AddContainer("MPointingPos",  "Events");
    16701861      write.AddContainer("ThetaOrig",     "Events");
    16711862      write.AddContainer("MSrcPosCam",    "Events");
    … …  
    19922183      bin.SetEdges(10, 0., 1.0);
    19932184
    1994       MH3 mh3("cos(MMcEvt.fTelescopeTheta)");
     2185      MH3 mh3("cos((MPointingPos.fZd)/kRad2Deg)");
    19952186      mh3.SetName(mname);
    19962187      MH::SetBinning(&mh3.GetHist(), &bin);
    … …  
    23552546
    23562547      write.AddContainer("MRawRunHeader", "RunHeaders");
    2357       //write.AddContainer("MTime",         "Events");
    2358       write.AddContainer("MMcEvt",        "Events");
    2359       //write.AddContainer("ThetaOrig",     "Events");
     2548      write.AddContainer("MTime",         "Events");
     2549      write.AddContainer("MPointingPos",  "Events");
    23602550      write.AddContainer("MSrcPosCam",    "Events");
    23612551      write.AddContainer("MSigmabar",     "Events");
    … …  
    31373327
    31383328
    3139     MFillH filler("MHMcCollectionArea", "MMcEvt");
     3329    MFillH filler("MHMcCollectionArea", "MPointingPos");
    31403330    filler.SetName("CollectionArea");
    31413331
    … …  
    33263516      write.AddContainer("MRawRunHeader", "RunHeaders");
    33273517      write.AddContainer("MTime",         "Events");
    3328       write.AddContainer("MMcEvt",        "Events");
     3518      write.AddContainer("MPointingPos",  "Events");
    33293519      write.AddContainer("ThetaOrig",     "Events");
    33303520      write.AddContainer("MSrcPosCam",    "Events");
    … …  
    34003590    // new from Robert
    34013591
    3402     MFillH hfill6("MHTimeDiffTheta", "MMcEvt");
     3592    MFillH hfill6("MHTimeDiffTheta", "MPointingPos");
    34033593    hfill6.SetName("HTimeDiffTheta");
    34043594
    3405     MFillH hfill6a("MHTimeDiffTime", "MMcEvt");
     3595    MFillH hfill6a("MHTimeDiffTime", "MPointingPos");
    34063596    hfill6a.SetName("HTimeDiffTime");
    34073597
    … …  
    36113801
    36123802
    3613 
    3614 
    3615 
    3616 
    3617 
    3618 
    3619 
    3620 
    3621 
    3622 
    3623 
    3624 
    3625 
    3626 
    3627 
    3628 
    3629 
    3630 
    3631 
    3632 
    3633 
    3634 
     3803   
     3804
     3805
     3806
     3807
     3808
     3809
     3810
     3811
     3812
     3813
     3814
     3815
     3816
     3817
     3818
     3819
     3820
     3821
     3822
     3823
     3824
Note: See TracChangeset for help on using the changeset viewer.