LArTimePhysPrediction.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
#include "LArCalibUtils/LArTimePhysPrediction.h"
#include "CaloIdentifier/CaloGain.h"
#include "LArRawConditions/LArCaliWaveContainer.h"
#include "LArRawConditions/LArPhysCaliTdiffComplete.h"
#include "LArRawConditions/LArWaveHelper.h"
#include "CaloDetDescr/CaloDepthTool.h"
#include "GaudiKernel/ToolHandle.h"
#include "CaloIdentifier/CaloCell_ID.h"
#include "CaloDetDescr/CaloDetDescrManager.h"
#include "CaloDetDescr/CaloDetDescriptor.h"
#include "CaloDetDescr/CaloDetDescrElement.h"
#include "CaloIdentifier/LArEM_ID.h"
#include "CaloIdentifier/LArHEC_ID.h"
#include "CaloIdentifier/LArFCAL_ID.h"
#include "LArIdentifier/LArOnlineID.h"
#include "GaudiKernel/INTupleSvc.h"
#include "GaudiKernel/NTuple.h"
#include "GaudiKernel/IToolSvc.h"
#include <TMath.h>
 
using CaliCellIt = LArCaliWaveContainer::ConstConditionsMapIterator;
using CaliWaveIt = LArCaliWaveContainer::LArCaliWaves::const_iterator;
 
LArTimePhysPrediction::LArTimePhysPrediction (const std::string& name, ISvcLocator* pSvcLocator) : 
  AthAlgorithm(name, pSvcLocator),
  m_CaloDepthTool(nullptr),
  m_nt(nullptr)  
{
  declareProperty("KeyInput", m_keyinput="LArCaliWave");
  declareProperty("KeyOutput", m_keyoutput="LArPhysCaliTdiff");
  declareProperty("vLCalib_EMB",m_vLCalib_EMB);
  declareProperty("vEtaMin_EMB",m_vEtaMin_EMB);
  declareProperty("vEtaMax_EMB",m_vEtaMax_EMB);
  declareProperty("vLCalib_EMEC",m_vLCalib_EMEC = 0);
  declareProperty("vLCalib_HEC",m_vLCalib_HEC = 0);
  declareProperty("vLSignal_EMB",m_vLSignal_EMB);
  declareProperty("vLSignal_EMEC",m_vLSignal_EMEC = 0);
  declareProperty("vLSignal_HEC",m_vLSignal_HEC = 0);
  declareProperty("vLSignal_FCAL",m_vLSignal_FCAL);
  declareProperty("vDeltaTTC_EMB",m_vDeltaTTC_EMB);
  declareProperty("vDeltaTTC_EC",m_vDeltaTTC_EC);
  declareProperty("vDeltaTTC_ECC_SPEC",m_vDeltaTTC_ECC_SPEC);
  declareProperty("vDeltaTTC_ECA_SPEC",m_vDeltaTTC_ECA_SPEC);
  declareProperty("sCalib",m_sCalib=5.5);
  declareProperty("sSignal",m_sSignal=5.5);
  declareProperty("sPig",m_sPig=5.5);
  declareProperty("sLTP",m_sLTP=5.5);
  declareProperty("nchannels_max",m_nchannels_max=1000000);
  declareProperty("GroupingType",m_groupingType="SubDetector");
}
 
LArTimePhysPrediction::~LArTimePhysPrediction() 
= default;
 
StatusCode LArTimePhysPrediction::initialize() 
{
  ATH_MSG_INFO ( "LArTimePhysPrediction in initialize()" );
  ATH_CHECK( m_thistSvc.retrieve() );
  
  ATH_CHECK( m_cablingKey.initialize() );
  ATH_CHECK( m_calibMapKey.initialize() );
  ATH_CHECK( m_caloMgrKey.initialize() );
 
  //Initialize ntuples
  NTupleFilePtr file1(ntupleSvc(),"/NTUPLES/FILE1");
  if (!file1){
    ATH_MSG_ERROR ( "Could not get NTupleFilePtr: failed" );
    return StatusCode::FAILURE;
  }
  NTuplePtr nt(ntupleSvc(),"/NTUPLES/FILE1/MyNtuple");
  if (!nt) {
    nt=ntupleSvc()->book("/NTUPLES/FILE1/MyNtuple",CLID_ColumnWiseTuple,"Timing ntuple");
  }
  if (!nt){
    ATH_MSG_ERROR ( "Booking of NTuple failed" );
    return StatusCode::FAILURE; 
  }
  
  m_nt=nt;
    
  //Book leaves 
  ATH_CHECK( nt->addItem("Chid",m_Chid) );
  ATH_CHECK( nt->addItem("Channel",m_Channel) );
  ATH_CHECK( nt->addItem("CalibLine",m_CalibLine) );
  ATH_CHECK( nt->addItem("is_lar_em",m_is_lar_em) );
  ATH_CHECK( nt->addItem("is_lar_hec",m_is_lar_hec) );
  ATH_CHECK( nt->addItem("is_lar_fcal",m_is_lar_fcal) );
  ATH_CHECK( nt->addItem("pos_neg",m_pos_neg) );
  ATH_CHECK( nt->addItem("barrel_ec",m_barrel_ec) );
  ATH_CHECK( nt->addItem("FT",m_FT) );
  ATH_CHECK( nt->addItem("FEBid",m_FEBid) );
  ATH_CHECK( nt->addItem("slot",m_slot) );
  ATH_CHECK( nt->addItem("eta",m_eta) );
  ATH_CHECK( nt->addItem("phi",m_phi) );
  ATH_CHECK( nt->addItem("layer",m_layer) );
  ATH_CHECK( nt->addItem("real_eta",m_real_eta) );
  ATH_CHECK( nt->addItem("real_phi",m_real_phi) );
  ATH_CHECK( nt->addItem("t0",m_t0) );
  ATH_CHECK( nt->addItem("tcali",m_tcali) );
  ATH_CHECK( nt->addItem("tCalibPredicted",m_tCalibPredicted) );
  ATH_CHECK( nt->addItem("CalibCables",m_CalibCables) );
  ATH_CHECK( nt->addItem("SignalCables",m_SignalCables) );
  ATH_CHECK( nt->addItem("TOF",m_TOF) );
  ATH_CHECK( nt->addItem("DeltaTTC",m_DeltaTTC) );
  ATH_CHECK( nt->addItem("tPhysPredicted",m_tPhysPredicted) );
  
  return StatusCode::SUCCESS;
}
 
StatusCode LArTimePhysPrediction::stop() 
{
  ATH_MSG_INFO ( "LArTimePhysPrediction in stop()" );
  //Intermediate variables declaration (should be removed in an updated version)
  const double meter2ns = 3.33564095;
  int Channel;
  int CalibLine;
  double tphys;
  double tcali;
  double dt;
  double TOF;
  double CalibCables=0;
  double LSignalFCAL=0;
  double DeltaTTC;
  double t0;
  int barrel_ec;
  int layer;
  int eta;
  int phi;  
  int FT;  
  int slot;
  int pos_neg;
  double real_eta;
  double real_phi;
  double radius=0;
  
  //Retrieve the LArCaliwaveContainer
  const LArCaliWaveContainer* larCaliWaveContainer = nullptr;
  ATH_CHECK( detStore()->retrieve(larCaliWaveContainer, m_keyinput) );
  ATH_MSG_INFO ( "Loaded LArCaliWaveContainer with key = " << m_keyinput );
 
  //Create the LArPhysCaliTdiffComplete object
  LArPhysCaliTdiffComplete *larPhysCaliTdiffComplete = new LArPhysCaliTdiffComplete();
  ATH_CHECK( larPhysCaliTdiffComplete->setGroupingType(m_groupingType,msg()) );
  ATH_CHECK( larPhysCaliTdiffComplete->initialize() );
 
  IAlgTool* algTool = nullptr;
  ATH_CHECK( toolSvc()->retrieveTool("CaloDepthTool", algTool, this) );
  m_CaloDepthTool=dynamic_cast<CaloDepthTool*>(algTool);
  ATH_MSG_INFO ( "CaloDepthTool retrieved with name " << m_CaloDepth );
  
  SG::ReadCondHandle<LArCalibLineMapping> clHdl{m_calibMapKey};
  ATH_CHECK(clHdl.isValid());
  const LArCalibLineMapping *clCont{*clHdl};
 
  SG::ReadCondHandle<LArOnOffIdMapping> cablingHdl{m_cablingKey};
  ATH_CHECK(cablingHdl.isValid());
  const LArOnOffIdMapping* cabling{*cablingHdl};
 
  SG::ReadCondHandle<CaloDetDescrManager> caloMgrHandle{m_caloMgrKey};
  ATH_CHECK(caloMgrHandle.isValid());
  const CaloDetDescrManager* caloDDM{*caloMgrHandle};
  
  //Define helpers
  LArWaveHelper larWaveHelper;
  
  const LArOnlineID* onlineHelper = nullptr;
  ATH_CHECK( detStore()->retrieve(onlineHelper, "LArOnlineID") );
  
  const CaloCell_ID* caloCID = nullptr;
  ATH_CHECK( detStore()->retrieve(caloCID, "CaloCell_ID") );
 
  //Get identifiers
  const LArEM_ID* emId = caloCID->em_idHelper();
  const LArHEC_ID* hecId = caloCID->hec_idHelper();
  const LArFCAL_ID* fcalId = caloCID->fcal_idHelper();
  
  //------------------------------------------------------------------------------------------------------------------------------
  //--------------Start to loop on the LArCaliWaveContainer------------------------------------------------------------------------
    for ( unsigned gain_it = CaloGain::LARHIGHGAIN ; gain_it < CaloGain::LARNGAIN ; ++gain_it ) { // Gains
      
      CaliCellIt cell_it = larCaliWaveContainer->begin(gain_it) ;
      CaliCellIt cell_it_e = larCaliWaveContainer->end(gain_it) ;
      
      if ( cell_it == cell_it_e ) {
    ATH_MSG_DEBUG ( "LArCaliWaveContainer (key = " << m_keyinput << ") has no wave with gain = " << gain_it );
    continue;
      } else {
    ATH_MSG_INFO ( "Processing LArCaliWaveContainer (key = " << m_keyinput << ") in gain = " << gain_it );
      }
      
      //counters for channels and waves
      int nchannels = 0;
      
      for ( ; cell_it != cell_it_e; ++cell_it) { // Channels
    
    if(nchannels==m_nchannels_max){continue;}
    std::cout << "nchannels " << nchannels << std::endl;
    nchannels++;
 
    CaliWaveIt wave_it = cell_it->begin();
    CaliWaveIt wave_it_e = cell_it->end();
    if ( wave_it == wave_it_e ) {
      ATH_MSG_DEBUG ( "Empty channel found..." );
      continue; // skip empty channels
    } 
      
    //use HWIdentifier
    HWIdentifier chid = cell_it.channelId();
 
    Identifier id;
    
    try {
      id = cabling->cnvToIdentifier(chid);   
    } catch ( const LArID_Exception& ) {
      ATH_MSG_ERROR ( "LArCabling exception caught for channel " << MSG::hex << chid << MSG::dec );
      continue;
    }
 
    Channel = onlineHelper->channel(chid);
 
    const std::vector<HWIdentifier>& calibLineV = clCont->calibSlotLine(chid);
        std::vector<HWIdentifier>::const_iterator calibLineIt = calibLineV.begin();   
    CalibLine = onlineHelper->channel(*calibLineIt) ;
    
    //indexes eta/phi/layer
    if (emId->is_lar_em(id)) {
      eta=emId->eta(id); 
      phi=emId->phi(id);
      layer=emId->sampling(id);}
    else if (hecId->is_lar_hec(id)) {
      eta=hecId->eta(id); 
      phi=hecId->phi(id);
      layer=hecId->sampling(id);}
    else if (fcalId->is_lar_fcal(id)) {
      eta=fcalId->eta(id); 
      phi=fcalId->phi(id);
      layer=fcalId->module(id);}
    else {
          ATH_MSG_INFO ( "cell not in the calorimeters " );
          continue;
        }
      
    //identification using the online helper (common to all LAr calorimeters)
    pos_neg = onlineHelper->pos_neg(chid);
    barrel_ec = onlineHelper->barrel_ec(chid);
    FT   = onlineHelper->feedthrough(chid);
    slot = onlineHelper->slot(chid);
 
    //get the FT online identifier
    HWIdentifier febid = onlineHelper->feb_Id(chid);
    m_FEBid = febid.get_identifier32().get_compact();
    
    //real eta and phi: need the hash identifier
    IdentifierHash theHash = caloCID->calo_cell_hash(id) ;
    const CaloDetDescrElement* theDDE = caloDDM->get_element(theHash) ;
    
    if(theDDE==nullptr) {
      ATH_MSG_INFO ( "CellIndex =  " << theHash << " has a DDE pointer NULL " );
      continue;
    }
    
    real_eta = theDDE->eta();
    real_phi = theDDE->phi();
    CaloCell_ID::CaloSample sample = theDDE->getSampling();//ok for EMB and EMEC
    //retrieve/compute the shower depth
    //use the "best" available parametrisation of the shower depth at the time of development
    //WARNING: use the CaloDepthTool's convention radius=r(barrel), radius=z(end-cap)
    //for HEC and FCAL: lengths could be moved in the job options
    if(emId->is_lar_em(id) && m_CaloDepthTool){
      radius = CaloDepthTool::cscopt2_parametrized(sample,real_eta,real_phi,caloDDM);
    }
    else if(hecId->is_lar_hec(id)){//assumption: "arrival point" = middle of the compartment 
      if(layer==0) radius=4398.;
      if(layer==1) radius=4806.;
      if(layer==2) radius=5359.;
      if(layer==3) radius=5840.;
    }
    else if(fcalId->is_lar_fcal(id)){//assumption: "arrival point" = middle of the compartment 
      if(layer==1) {radius=4916.;LSignalFCAL=m_vLSignal_FCAL[0];}
      if(layer==2) {radius=5366.;LSignalFCAL=m_vLSignal_FCAL[1];}
      if(layer==3) {radius=5816.;LSignalFCAL=m_vLSignal_FCAL[2];}
    }
    
    for ( ; wave_it != wave_it_e; ++wave_it) { // DACs <==> iterator = the caliwave
      
      //initialize tphys
      tphys=-999.;
 
      //----step 1: get the calibration time
      dt = wave_it->getDt();
      tcali  = dt * larWaveHelper.getMax(*wave_it);
      //modified version of getT0
      const unsigned Nbase=5;
      double rT0=0.;
      double asamp1=0.;
      double asamp2=0.;
      double base = larWaveHelper.getBaseline(*wave_it,Nbase);
      double amax = larWaveHelper.getMaxAmp(*wave_it);
      unsigned int imax=larWaveHelper.getMax(*wave_it);
      
      for (unsigned int i=1; i<imax ; i++) {
        asamp1=wave_it->getSample(i);
        asamp2=wave_it->getSample(i-1);
        if((asamp1-base)>amax*0.10 && (asamp2-base)<amax*0.10)
          rT0 = i-1+(0.01*amax-(asamp2-base))/(asamp1-asamp2);
      }
      t0=rT0;
          
      //----step 2: compute the TOF
      if(emId->is_lar_em(id) && barrel_ec==0) TOF = fabs(radius)*TMath::CosH(real_eta)* meter2ns/1000;//EMB
      else TOF = fabs(radius)/TMath::TanH(fabs(real_eta))* meter2ns/1000;//EC
      
      //----step 3: deduction of the propagation time due to the propagation of the calibration pulse
      //from the pulse until the calibration board
      if(emId->is_lar_em(id) && barrel_ec==0){//EMB
        for(int ieta=0;ieta<16;ieta++){
          if(fabs(real_eta)>ieta*0.1 && fabs(real_eta)<(ieta+1)*0.1) CalibCables=0.7*m_sPig+m_vLCalib_EMB[layer][ieta]*m_sCalib; 
        }
      }
      else if(emId->is_lar_em(id) && abs(barrel_ec)==1) CalibCables=0.9*m_sPig+m_vLCalib_EMEC*m_sCalib;//EMEC
      else if(hecId->is_lar_hec(id)) CalibCables=m_vLCalib_HEC*m_sCalib;//HEC
      else if(fcalId->is_lar_fcal(id)) CalibCables=-LSignalFCAL*m_sSignal;//FCAL
              
      //----step 4: deduction of the propagation times due to the optical fibers from USA15 to FEC
      if(emId->is_lar_em(id) && barrel_ec==0){
        if(pos_neg==1) DeltaTTC=m_vDeltaTTC_EMB[1][FT];//EMBA
        else DeltaTTC=m_vDeltaTTC_EMB[0][FT];//EMBC 
      }
      else {
        if(real_eta>0) DeltaTTC=m_vDeltaTTC_EC[1][FT];//ECA
        else DeltaTTC=m_vDeltaTTC_EC[0][FT];//ECC
        //correction for special crates
        //A FTs 2/3
        if(emId->is_lar_em(id) && abs(barrel_ec)==1 && FT==2 && real_eta>0 && slot<11) DeltaTTC=m_vDeltaTTC_ECA_SPEC[0][0];
        if((emId->is_lar_em(id) && abs(barrel_ec)==1 && FT==2 && real_eta>0 && slot>10) || (emId->is_lar_em(id) && abs(barrel_ec)==1 && FT==3 && real_eta>0 && slot<3)) DeltaTTC=m_vDeltaTTC_ECA_SPEC[0][1];
        if(hecId->is_lar_hec(id) && abs(barrel_ec)==1 && FT==3 && real_eta>0) DeltaTTC=m_vDeltaTTC_ECA_SPEC[0][2];
        //A FTs 9/10
        if(emId->is_lar_em(id) && abs(barrel_ec)==1 && FT==9 && real_eta>0 && slot<11) DeltaTTC=m_vDeltaTTC_ECA_SPEC[1][0];
        if((emId->is_lar_em(id) && abs(barrel_ec)==1 && FT==9 && real_eta>0 && slot>10) || (emId->is_lar_em(id) && abs(barrel_ec)==1 && FT==10 && real_eta>0 && slot<3)) DeltaTTC=m_vDeltaTTC_ECA_SPEC[1][1];
        if(hecId->is_lar_hec(id) && abs(barrel_ec)==1 && FT==10 && real_eta>0) DeltaTTC=m_vDeltaTTC_ECA_SPEC[1][2];
        //A FTs 15/16
        if(emId->is_lar_em(id) && abs(barrel_ec)==1 && FT==15 && real_eta>0 && slot<11) DeltaTTC=m_vDeltaTTC_ECA_SPEC[2][0];
        if((emId->is_lar_em(id) && abs(barrel_ec)==1 && FT==15 && real_eta>0 && slot>10) || (emId->is_lar_em(id) && abs(barrel_ec)==1 && FT==16 && real_eta>0 && slot<3)) DeltaTTC=m_vDeltaTTC_ECA_SPEC[2][1];
        if(hecId->is_lar_hec(id) && abs(barrel_ec)==1 && FT==16 && real_eta>0) DeltaTTC=m_vDeltaTTC_ECA_SPEC[2][2];
        //A FTs 21/22
        if(emId->is_lar_em(id) && abs(barrel_ec)==1 && FT==21 && real_eta>0 && slot<11) DeltaTTC=m_vDeltaTTC_ECA_SPEC[3][0];
        if((emId->is_lar_em(id) && abs(barrel_ec)==1 && FT==21 && real_eta>0 && slot>10) || (emId->is_lar_em(id) && abs(barrel_ec)==1 && FT==21 && real_eta>0 && slot<3)) DeltaTTC=m_vDeltaTTC_ECA_SPEC[3][1];
        if(hecId->is_lar_hec(id) && abs(barrel_ec)==1 && FT==22 && real_eta>0) DeltaTTC=m_vDeltaTTC_ECA_SPEC[3][2];
        //C FTs 2/3
        if(emId->is_lar_em(id) && abs(barrel_ec)==1 && FT==2 && real_eta<0 && slot<11) DeltaTTC=m_vDeltaTTC_ECC_SPEC[0][0];
        if((emId->is_lar_em(id) && abs(barrel_ec)==1 && FT==2 && real_eta<0 && slot>10) || (emId->is_lar_em(id) && abs(barrel_ec)==1 && FT==3 && real_eta<0 && slot<3)) DeltaTTC=m_vDeltaTTC_ECC_SPEC[0][1];
        if(hecId->is_lar_hec(id) && abs(barrel_ec)==1 && FT==3 && real_eta>0) DeltaTTC=m_vDeltaTTC_ECC_SPEC[0][2];
        //C FTs 9/10
        if(emId->is_lar_em(id) && abs(barrel_ec)==1 && FT==9 && real_eta<0 && slot<11) DeltaTTC=m_vDeltaTTC_ECC_SPEC[1][0];
        if((emId->is_lar_em(id) && abs(barrel_ec)==1 && FT==9 && real_eta<0 && slot>10) || (emId->is_lar_em(id) && abs(barrel_ec)==1 && FT==10 && real_eta<0 && slot<3)) DeltaTTC=m_vDeltaTTC_ECC_SPEC[1][1];
        if(hecId->is_lar_hec(id) && abs(barrel_ec)==1 && FT==10 && real_eta<0) DeltaTTC=m_vDeltaTTC_ECC_SPEC[1][2];
        //C FTs 15/16
        if(emId->is_lar_em(id) && abs(barrel_ec)==1 && FT==15 && real_eta<0 && slot<11) DeltaTTC=m_vDeltaTTC_ECC_SPEC[2][0];
        if((emId->is_lar_em(id) && abs(barrel_ec)==1 && FT==15 && real_eta<0 && slot>10) || (emId->is_lar_em(id) && abs(barrel_ec)==1 && FT==16 && real_eta<0 && slot<3)) DeltaTTC=m_vDeltaTTC_ECC_SPEC[2][1];
        if(hecId->is_lar_hec(id) && abs(barrel_ec)==1 && FT==16 && real_eta<0) DeltaTTC=m_vDeltaTTC_ECC_SPEC[2][2];
        //C FTs 21/22
        if(emId->is_lar_em(id) && abs(barrel_ec)==1 && FT==21 && real_eta<0 && slot<11) DeltaTTC=m_vDeltaTTC_ECC_SPEC[3][0];
        if((emId->is_lar_em(id) && abs(barrel_ec)==1 && FT==21 && real_eta<0 && slot>10) || (emId->is_lar_em(id) && abs(barrel_ec)==1 && FT==21 && real_eta<0 && slot<3)) DeltaTTC=m_vDeltaTTC_ECC_SPEC[3][1];
        if(hecId->is_lar_hec(id) && abs(barrel_ec)==1 && FT==22 && real_eta<0) DeltaTTC=m_vDeltaTTC_ECC_SPEC[3][2];
      }                  
      
      //----step 5: computation of tphys
      //tphys=tcali-CalibCables+TOF-DeltaTTC;
      //add the effect from LTP cables: LTPI->LTP and LTP->LTP (hard coded)
      //for A-C sides difference: a 8 ns is set (derived from September 2008 splash events)
      if(pos_neg==1 && emId->is_lar_em(id) && abs(barrel_ec)==0) tphys=tcali-CalibCables+TOF+DeltaTTC+8+0.3*m_sLTP;//EMBA
      else if(pos_neg==0 && emId->is_lar_em(id) && abs(barrel_ec)==0) tphys=tcali-CalibCables+TOF+DeltaTTC+0.3*m_sLTP;//EMBC
      else if(pos_neg==1 && emId->is_lar_em(id) && abs(barrel_ec)==1) tphys=tcali-CalibCables+TOF+DeltaTTC+8;//EMECA
      else if(real_eta>0 && (emId->is_lar_hec(id)||emId->is_lar_fcal(id))) tphys=tcali-CalibCables+TOF+DeltaTTC+8+0.3*m_sLTP;//HEC,FCAL A
      else if(real_eta<0 && (emId->is_lar_hec(id)||emId->is_lar_fcal(id))) tphys=tcali-CalibCables+TOF+DeltaTTC+0.3*m_sLTP;//HEC,FCAL C 
      else tphys=tcali-CalibCables+TOF+DeltaTTC;//EMECC
      
      //prediction of the expected calibration time
      //not finalized (EMEC signal cables lengths not implemented)
      if(emId->is_lar_em(id) && abs(barrel_ec)==0){//EMB
        for(int ieta=0;ieta<16;ieta++){
          if(fabs(real_eta)>ieta*0.1 && fabs(real_eta)<(ieta+1)*0.1 && (layer==0 || (layer==1 && fabs(real_eta)<0.6))) m_SignalCables=0.9*m_sPig+m_vLSignal_EMB[layer][ieta]*m_sSignal;
          else if(fabs(real_eta)>ieta*0.1 && fabs(real_eta)<(ieta+1)*0.1) m_SignalCables=0.7*m_sPig+m_vLSignal_EMB[layer][ieta]*m_sSignal;
        }
      }
      else if(emId->is_lar_em(id) && abs(barrel_ec)==1) m_SignalCables=0.9*m_sPig+m_vLSignal_EMEC*m_sSignal;//EMEC
      else if(hecId->is_lar_hec(id)) m_SignalCables=m_vLSignal_HEC*m_sSignal;//HEC
      else if(fcalId->is_lar_fcal(id)) m_SignalCables=0.;//FCAL
      
      m_tCalibPredicted=m_SignalCables+CalibCables+tcali-t0;
      if(fcalId->is_lar_fcal(id))m_tCalibPredicted=tcali-t0;
      
      //fill larPhysCaliTdiffComplete (needed to fill the DB)
      //larPhysCaliTdiffComplete->set(chid,gain_it,tphys);
      
      //for debugging purposes: check the computed times
      //write the ntuple
      if(emId->is_lar_em(id)) m_is_lar_em=1;
      else m_is_lar_em=0;
      if(emId->is_lar_hec(id)) m_is_lar_hec=1;
      else m_is_lar_hec=0;
      if(emId->is_lar_fcal(id)) m_is_lar_fcal=1;
      else m_is_lar_fcal=0;
      m_Chid=chid.get_identifier32().get_compact();
      m_Channel=Channel;
      m_CalibLine=CalibLine;
      m_pos_neg=pos_neg;
      m_barrel_ec=barrel_ec;
      m_FT=FT;
      m_slot=slot;
      m_eta=eta;
      m_phi=phi;
      m_layer=layer;
      m_real_eta=real_eta;
      m_real_phi=real_phi;
      m_t0=t0;
      m_tcali=tcali;
      m_CalibCables=CalibCables;
      m_TOF=TOF;
      m_DeltaTTC=DeltaTTC;
      m_tPhysPredicted=tphys;
      //m_SignalCables: already filled
      //m_tCalibPredicted: already filled (computation to be finalized) 
      
      ATH_MSG_VERBOSE ( "Try to write to ntuple " );
      ATH_CHECK( ntupleSvc()->writeRecord(m_nt) );
          ATH_MSG_VERBOSE ( "wave written to ntuple" );
    } // end of loop over DACs  
    
      }// end of loop over Channels
      
    } // end of loop over Gains
    
    // Record LArPhysCaliTdiffComplete (needed to fill the DB)
    //-----------------------------------end of the main loop-------------------------------------------------------------------------
    //--------------------------------------------------------------------------------------------------------------------------------
    
    ATH_MSG_INFO ( "end of stop" );
    return StatusCode::SUCCESS ;
}