/[PAMELA software]/quicklook/QLflightS4_ND/S4_QL.cpp
ViewVC logotype

Diff of /quicklook/QLflightS4_ND/S4_QL.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.1 by pam-rm2, Thu Jun 15 14:00:30 2006 UTC revision 1.3 by pam-rm2, Wed Jul 12 09:06:35 2006 UTC
# Line 260  void S4_QL(TString base, TString outDir, Line 260  void S4_QL(TString base, TString outDir,
260    headBr->GetEntry(0);    headBr->GetEntry(0);
261    ph = eh->GetPscuHeader();    ph = eh->GetPscuHeader();
262    firstime = ph->GetOrbitalTime();    firstime = ph->GetOrbitalTime();
263    headBr->GetEntry(nevents-1);  
264    ph = eh->GetPscuHeader();    for (Int_t i = 0; i < nevents; i++){
265    lastime = ph->GetOrbitalTime();      headBr->GetEntry(i);
266        ph = eh->GetPscuHeader();
267        if((ph->GetOrbitalTime()) <= firstime) firstime=ph->GetOrbitalTime();
268        if((ph->GetOrbitalTime()) >= lastime) lastime=ph->GetOrbitalTime();
269      }
270    obmin=firstime;    obmin=firstime;
271    obmax=lastime;      obmax=lastime;  
272    const ULong_t nint=(((lastime-firstime)/(DeltaT*60000)));    const ULong_t nint=(((lastime-firstime)/(DeltaT*60000)));
# Line 279  void S4_QL(TString base, TString outDir, Line 283  void S4_QL(TString base, TString outDir,
283    title=filename+": S4 rate from Trigger Packet";    title=filename+": S4 rate from Trigger Packet";
284    
285    const ULong_t nint3=(lastime-firstime);    const ULong_t nint3=(lastime-firstime);
286    const Int_t size3 = (Int_t)((nint3)/3);    const Int_t size3 = (Int_t)((nint3)/100);
287    
288    TH1F *rate= new TH1F(title.c_str(), title.c_str(), nint3, obmin, obmax);    TH1F *rate= new TH1F(title.c_str(), title.c_str(), nint3, obmin, obmax);
289    TH1F *rateline= new TH1F(filename+": smoothed S4 rate", filename+": smoothed S4 rate", size3, obmin, obmax);    TH1F *rateline= new TH1F(filename+". S4 rate from Trigger Packet: mean value over 100 events", filename+". S4 rate from Trigger Packet: mean value over 100 events", size3, obmin, obmax);
290    
291    //------------------------------------------------------------------------------------------------------------------    //------------------------------------------------------------------------------------------------------------------
292    //------- fill histograms ---------//    //------- fill histograms ---------//
# Line 299  void S4_QL(TString base, TString outDir, Line 303  void S4_QL(TString base, TString outDir,
303      h1->Fill(s4e->S4_DATA);        h1->Fill(s4e->S4_DATA);  
304      Allev->Fill(obt,s4e->S4_DATA);      Allev->Fill(obt,s4e->S4_DATA);
305      Alltime->Fill(obt);      Alltime->Fill(obt);
306      if((i==n) && (s4rate >0)){      if ((trige->patterntrig[0] == 0) && (trige->patterntrig[1] != 0) &&(trige->patterntrig[2] == 0) && (trige->patterntrig[3] == 0) && (trige->patterntrig[4] == 0) && (trige->patterntrig[5] == 0)){
       rateline->Fill(obt,s4rate);  
     n=n+5;  
       }  
     if ((trige->patterntrig[0] == 0) && (trige->patterntrig[2] == 0) && (trige->patterntrig[3] == 0) && (trige->patterntrig[4] == 0) &&     (trige->patterntrig[5] == 0)){  
307        h2->Fill(s4e->S4_DATA);        h2->Fill(s4e->S4_DATA);
308        S4ev->Fill(obt,s4e->S4_DATA);        S4ev->Fill(obt,s4e->S4_DATA);
309        S4time->Fill(obt);        S4time->Fill(obt);
310        p=p+1;        p=p+1;
311      }      }
312    }    }
313      Int_t kk=0;
314      while (kk < nevents){
315        obt=0;
316        s4rate=0;
317        for(Int_t jj=kk; jj< (kk+100); jj++){
318          TriggerBr->GetEntry(jj);
319          headBr->GetEntry(jj);
320          obt = obt+(ph->GetOrbitalTime());  
321          s4rate= s4rate+(trige->s4calcount[0]);
322        }
323        rateline->Fill((obt/100),(s4rate/100));
324        kk=kk+100;
325      }
326    
327    //****************************** Canvases *******************************//    //****************************** Canvases *******************************//
328    //------------------- First Canvas --------------------------------//    //------------------- First Canvas --------------------------------//
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.23