JCalibrateNB.cc

Go to the documentation of this file.

#include <iostream>

#include "JDAQ/JDAQEventIO.hh"
#include "Jeep/JParser.hh"
#include "JSupport/JSupport.hh"
#include "JSupport/JMultipleFileScanner.hh"
#include "JSupport/JSupportToolkit.hh"
#include "JTrigger/JHitL0.hh"
#include "JTrigger/JBuildL0.hh"
#include "JDetector/JDetector.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JDetector/JModuleRouter.hh"
#include "JGizmo/JManager.hh"
#include "JPhysics/JDispersion.hh"

#include "TH1.h"
#include "TH2.h"
#include "TFile.h"

using namespace std;
using namespace JPP;
using namespace JDETECTOR;
using namespace KM3NETDAQ;

int main(int argc , char** argv){

  JMultipleFileScanner<JDAQEvent>  inputFiles;  
  JLimit_t&                        numberOfEvents = inputFiles.getLimit();
  string                           detectorFile;
  string                           outFile;
  JDetector                        detector;
  double                           depth;
  double                           wavelength;
    
  try {

    JParser<> zap;
    
    zap['f'] = make_field(inputFiles                           );
    zap['a'] = make_field(detectorFile                         );
    zap['n'] = make_field(numberOfEvents                       ) = JLimit::max();
    zap['o'] = make_field(outFile                              ) = "out.root";
    zap['h'] = make_field(depth         , "detector depth [m]" ) = 2500.0;
    zap['l'] = make_field(wavelength    , "LED wavelength [nm]") = 400.0;

    zap(argc,argv);
  }
  catch(const exception &error) {
    ERROR(error.what() << endl);
  }

  typedef JHitL0 JHit_t;
  typedef vector<JHit_t> JHitBuffer;

  JBuildL0<JHit_t> builder;
  
  JHitBuffer triggerBuffer;
  JHitBuffer snapshotBuffer;
  
  load(detectorFile , detector);

  JModuleRouter moduleRouter(detector);

  JDispersion jd (0.1 * depth);

  double cw = getSpeedOfLight() / jd.getIndexOfRefractionGroup(wavelength);
  
  JManager < string , TH1D >* timeDifferences;
  timeDifferences = new JManager <string , TH1D> (new TH1D("%" , "" , 501 , -250.5 , 250.5));

  while(inputFiles.hasNext()){
    
    JDAQEvent* event = inputFiles.next();
  
    triggerBuffer .clear();
    snapshotBuffer.clear();
    
    builder(*event, moduleRouter, true  , back_inserter(snapshotBuffer));
    builder(*event, moduleRouter, false , back_inserter( triggerBuffer)); 

    sort(triggerBuffer.begin(),triggerBuffer.end(),less<JHit>());
        
    JHitBuffer::const_iterator triggeredHit = triggerBuffer.begin();

    int triggeredString = moduleRouter.getModule( triggeredHit -> getModuleID() ).getString() ;
    int triggeredFloor  = moduleRouter.getModule( triggeredHit -> getModuleID() ).getFloor() ;

    for(JHitBuffer::const_iterator snapshotHit = snapshotBuffer.begin(); snapshotHit != snapshotBuffer.end(); ++snapshotHit){
                         
      int snapshotString = moduleRouter.getModule( snapshotHit -> getModuleID() ).getString() ;
      int snapshotFloor  = moduleRouter.getModule( snapshotHit -> getModuleID() ).getFloor() ;

      double ToF = triggeredHit->getDistance(*snapshotHit) / cw;
      double Dt  = snapshotHit->getT() - triggeredHit->getT();

      (*timeDifferences)[MAKE_STRING("S" + to_string(triggeredString) + "F" + to_string(triggeredFloor) + "_S" + to_string(snapshotString) + "F" + to_string(snapshotFloor)).c_str()]->Fill(Dt-ToF);
    }
  }

  TFile output(outFile.c_str() , "recreate") ;
  timeDifferences -> Write(output);
}