detector.cxx

Go to the documentation of this file.

#include "detector.h"
#include "Plex/PlexHandle.h"
#include "Plex/PlexSEIdAltL.h"

Detector::Detector()
{
  memset(channels,0,sizeof(channels));
  first_time=0;
  
  const char *dir = gSystem->Getenv("LISUMMARY_OUTPUT_DIR");
  if (!dir) {
    cerr << "Set LISUMMARY_OUTPUT_DIR first!\n";
    exit(1);
  }
  
  sprintf(filename,"%s/combined.root",dir);
  outfile = new TFile(filename,"NEW");
  if(!outfile->IsOpen()){
    cout << " ** file already exists, making different name ** " << endl;
    sprintf(filename,"%s/combined_1.root",dir);
    outfile = new TFile(filename,"NEW");
  }
  cout << " ** Created output file " << filename << " ** " << endl;
  InitialiseTree();
  cout << " ** Made detector ** " << endl;
}



void Detector::Update(Int_t new_led,
                      Int_t new_pulse_height,
                      Int_t new_pulse_width,
                      Int_t new_pulserbox,
                      Int_t new_calibpoint,
                      VldTimeStamp new_time, 
                      Int_t new_calibtype,
                      Int_t new_pulses,
                      Int_t new_period){
  

  
  if(first_time==0){
    first_time=1;
    calibpoint = new_calibpoint;
    pulse_height = new_pulse_height;
    pulse_width = new_pulse_width;
    pulserbox = new_pulserbox;
    led = new_led;
    time = new_time;
    period = new_period;
    calibtype = new_calibtype;
    pulses = new_pulses;
    
    cout << " ** Found first point " << new_pulserbox << " :: " 
         << new_led << " at time " << new_time << " ** " << endl;

    cout << " ** opening database ** " << endl;
    myplexhandle = new PlexHandle(validity);
    
  }
  
  if( (new_led != led) 
      || (new_pulse_height != pulse_height)
      || (new_pulse_width != pulse_width)
      || (new_pulserbox != pulserbox)
      || (new_calibpoint != calibpoint)){
    
    cout << " ** Found new point " << new_pulserbox << " :: " 
         << new_led << " at time " << new_time << " ** " << endl;
    
    // Write out data here...
    Int_t channel_address;
    
    for(Int_t crate_loop=0;crate_loop<10;crate_loop++){
      for(Int_t chad_loop=0;chad_loop<6000;chad_loop++){
        channel_address = (crate_loop*6000) + chad_loop;
        if(channels[channel_address] !=0){
          timestamp = new_time.GetSec();

          mean = channels[channel_address]->GetMean();
          rms = channels[channel_address]->GetRMS();
          no_entries = channels[channel_address]->GetEntries();
          
          crate = channels[channel_address]->GetChannelId().GetCrate();
          chAdd = channels[channel_address]->GetChannelId().GetChAdd();
          vaarc = channels[channel_address]->GetChannelId().GetVarcId();
          vmm = channels[channel_address]->GetChannelId().GetVmm();
          vaadc = channels[channel_address]->GetChannelId().GetVaAdcSel();
          
          chip = channels[channel_address]->GetChannelId().GetVaChip();
          channel = channels[channel_address]->GetChannelId().GetVaChannel();
          plane = myplexhandle->GetSEIdAltL(channels[channel_address]->GetChannelId()).GetPlane();

          outtree->Fill();
        }
      }
    }
    // Clean all channels..

    memset(channels,0,sizeof(channels));
    


    first_time=1;
    calibpoint = new_calibpoint;
    pulse_height = new_pulse_height;
    pulse_width = new_pulse_width;
    pulserbox = new_pulserbox;
    led = new_led;
    period = new_period;
    calibtype = new_calibtype;
    pulses = new_pulses;


  }


};

void Detector::AddEntry(RawChannelId channelid,
                        Double_t mean,
                        Double_t rms,
                        Int_t no_entries
                        ){
  Int_t channel_address=0;
  channel_address = (channelid.GetCrate()*6000) + channelid.GetChAdd();
  
  if(channels[ channel_address ] ==0){
    channels[ channel_address ] = new Channel(channelid);
  }
  
  channels[channel_address]->AddEntry(no_entries,mean,rms);
    
};
 


void Detector::InitialiseTree(){
  cout << " Making Tree " << endl;
  outfile = new TFile(filename,"RECREATE");
  outtree = new TTree("li_tree","li_tree");
  

  outtree->Branch("calibpoint",&calibpoint,"calibpoint/I",32000);
  outtree->Branch("pulserbox",&pulserbox,"pulserbox/I",32000); 
  outtree->Branch("led",&led,"led/I",32000);
  outtree->Branch("pulse_height",&pulse_height,"pulse_height/I",32000);
  outtree->Branch("pulse_width",&pulse_width,"pulse_width/I",32000);
  outtree->Branch("period",&period,"period/I",32000);
  outtree->Branch("calibtype",&calibtype,"calibtype/I",32000);
  outtree->Branch("pulses",&pulses,"pulses/I",32000);
  outtree->Branch("timestamp",&timestamp,"timestamp/I",32000);


  outtree->Branch("crate",&crate,"crate/I",32000);
  outtree->Branch("chAdd",&chAdd,"chAdd/I",32000);
  outtree->Branch("vaarc",&vaarc,"vaarc/I",32000);
  outtree->Branch("vmm",&vmm,"vmm/I",32000);
  outtree->Branch("vaadc",&vaadc,"vaadc/I",32000);


  outtree->Branch("plane",&plane,"plane/I",32000);
  outtree->Branch("chip",&chip,"chip/I",32000);
  outtree->Branch("channel",&channel,"channel/I",32000);

  outtree->Branch("no_entries",&no_entries,"no_entries/I",32000);
  outtree->Branch("mean",&mean,"mean/D",32000);
  outtree->Branch("rms",&rms,"rms/D",32000);

  cout <<" ** Made Tree ** " << endl;

};

---