MDQ_Thresholds
[Low Level Data Quality Checks]

Record thresholds from the ODB. More...

Functions
INT	MDQ_Thresholds_init (void)
INT	MDQ_Thresholds (EVENT_HEADER *, void *)
INT	MDQ_Thresholds_eor (INT)
Variables
HNDLE	hDB
TGlobalData *	gData
	Object to hold data used and produced by modules throughout alcapana stage of analysis.
TSetupData *	gSetup
	Hardware information about digitizers and detectors to be used during alcapana stage of analysis.
map< std::string, TH1F * >	DQ_Thresholds_histograms_map
	Map of bank name to histogram. FADC histograms have 2 entries, one each for an upper and lower trigger threshold. CAENs have single entry.
ANA_MODULE	MDQ_Thresholds_module

---

Detailed Description

Record thresholds from the ODB.

Author:

Andrew Edmonds

Record thresholds from ODB as entries in a histogram per bank.

---

Function Documentation

INT MDQ_Thresholds	(	EVENT_HEADER *	pheader,
		void *	pevent
	)

This method fills the histograms

Definition at line 238 of file MDQ_Thresholds.cpp.

{
        // Get the event number
        int midas_event_number = pheader->serial_number;

        // Some typedefs
        typedef map<string, vector<TPulseIsland*> > TStringPulseIslandMap;
        typedef pair<string, vector<TPulseIsland*> > TStringPulseIslandPair;
        typedef map<string, vector<TPulseIsland*> >::iterator map_iterator;

        // Don't need anything here

        return SUCCESS;
}

INT MDQ_Thresholds_eor

(

INT

run_number

)

This method does any last minute things to the histograms at the end of the run

Definition at line 115 of file MDQ_Thresholds.cpp.

References DQ_Thresholds_histograms_map, TSetupData::fBankToDetectorMap, hDB, TSetupData::Instance(), and TSetupData::IsFADC().

                                       {

  // Get the run duration to scale the histogram
  HNDLE hDB, hKey;
  char keyName[200];

  if(cm_get_experiment_database(&hDB, NULL) != CM_SUCCESS){
    printf("Warning: Could not connect to ODB database!\n");
    return false;
  }

  // Loop through the channels and get the thresholds
  std::map<std::string, std::string> bank_to_detector_map = gSetup->fBankToDetectorMap;
  for(std::map<std::string, std::string>::iterator mapIter = bank_to_detector_map.begin(); 
      mapIter != bank_to_detector_map.end(); mapIter++) { 

    std::string bankname = mapIter->first;

    if (TSetupData::Instance()->IsFADC(bankname)) {
      // get the FADC thresholds (both upper and lower)

      // first get the channel and address from the bankname
      int iChn = (int)(bankname[1] - 97);
      std::string iAddr = bankname.substr(2, 2);

      // Get the lower threshold
      sprintf(keyName, "/Equipment/Crate 9/Settings/NFADC %s/Channel %d/Lower threshold", iAddr.c_str(), iChn);
      if(db_find_key(hDB,0,keyName, &hKey) != SUCCESS){
        printf("Warning: Could not find key %s\n", keyName);
        return false;
      }
      KEY lower_threshold_key;
      if(db_get_key(hDB, hKey, &lower_threshold_key) != DB_SUCCESS){
        printf("Warning: Could not find key %s\n", keyName);
        return false;
      }
      INT LowerThresholds[lower_threshold_key.num_values];
      int size = sizeof(LowerThresholds);
      if(db_get_value(hDB, 0, keyName, LowerThresholds, &size, TID_INT, 0) != DB_SUCCESS){
        printf("Warning: Could not retrieve values for key %s\n", keyName);
        return false;
      }


      // Upper threshold
      sprintf(keyName, "/Equipment/Crate 9/Settings/NFADC %s/Channel %d/Upper threshold", iAddr.c_str(), iChn);
      if(db_find_key(hDB,0,keyName, &hKey) != SUCCESS){
        printf("Warning: Could not find key %s\n", keyName);
        return false;
      }
      KEY upper_threshold_key;
      if(db_get_key(hDB, hKey, &upper_threshold_key) != DB_SUCCESS){
        printf("Warning: Could not find key %s\n", keyName);
        return false;
      }
      INT UpperThresholds[upper_threshold_key.num_values];
      size = sizeof(UpperThresholds);
      if(db_get_value(hDB, 0, keyName, UpperThresholds, &size, TID_INT, 0) != DB_SUCCESS){
        printf("Warning: Could not retrieve values for key %s\n", keyName);
        return false;
      }

      // Print the results
      DQ_Thresholds_histograms_map[bankname]->Fill("lower threshold", LowerThresholds[0]);
      DQ_Thresholds_histograms_map[bankname]->Fill("upper threshold", UpperThresholds[0]);
    }
    else if (TSetupData::Instance()->IsHoustonCAEN(bankname)) {
      // get the UH CAEN thresholds

      // first get the channel and address from the bankname
      int iChn = (int)(bankname[1] - 97);

      // Get the threshold
      sprintf(keyName, "/Equipment/Crate 4/Settings/CAEN0/Ch%d/trigger threshhold", iChn);
      if(db_find_key(hDB,0,keyName, &hKey) != SUCCESS){
        printf("Warning: Could not find key %s\n", keyName);
        return false;
      }
      KEY threshold_key;
      if(db_get_key(hDB, hKey, &threshold_key) != DB_SUCCESS){
        printf("Warning: Could not find key %s\n", keyName);
        return false;
      }
      DWORD Thresholds[threshold_key.num_values];
      int size = sizeof(Thresholds);
      if(db_get_value(hDB, 0, keyName, Thresholds, &size, TID_DWORD, 0) != DB_SUCCESS){
        printf("Warning: Could not retrieve values for key %s\n", keyName);
        return false;
      }
      DQ_Thresholds_histograms_map[bankname]->Fill("trigger threshold", Thresholds[0]);
    }
    else if (TSetupData::Instance()->IsBostonCAEN(bankname)) {
      // get the BU CAEN thresholds

      // first get the channel and address from the bankname
      int iChn = (int)(bankname[1] - 97);

      // Get the threshold
      sprintf(keyName, "/Equipment/Crate 5/Settings/CAEN/Ch0%d/self_trigger_threshhold", iChn);
      if(db_find_key(hDB,0,keyName, &hKey) != SUCCESS){
        printf("Warning: Could not find key %s\n", keyName);
        return false;
      }
      KEY threshold_key;
      if(db_get_key(hDB, hKey, &threshold_key) != DB_SUCCESS){
        printf("Warning: Could not find key %s\n", keyName);
        return false;
      }
      float Thresholds[threshold_key.num_values];
      int size = sizeof(Thresholds);
      if(db_get_value(hDB, 0, keyName, Thresholds, &size, TID_FLOAT, 0) != DB_SUCCESS){
        printf("Warning: Could not retrieve values for key %s\n", keyName);
        return false;
      }
      DQ_Thresholds_histograms_map[bankname]->Fill("self_trigger_threshold", Thresholds[0]);
    }
  }

  return SUCCESS;
}

INT MDQ_Thresholds_init

(

)

This method initializes histograms.

Definition at line 81 of file MDQ_Thresholds.cpp.

References DQ_Thresholds_histograms_map, TSetupData::fBankToDetectorMap, and TSetupData::GetDetectorName().

{
    // See if the DataQuality_LowLevel/ directory already exists
  if (!gDirectory->Cd("DataQuality_LowLevel")) {
    
    std::string dir_name("DataQuality_LowLevel/");
    gDirectory->mkdir(dir_name.c_str());
    gDirectory->Cd(dir_name.c_str());
  }

  // Create a histogram for each channel
  std::map<std::string, std::string> bank_to_detector_map = gSetup->fBankToDetectorMap;
  for(std::map<std::string, std::string>::iterator mapIter = bank_to_detector_map.begin(); 
      mapIter != bank_to_detector_map.end(); mapIter++) { 

    std::string bankname = mapIter->first;
    std::string detname = gSetup->GetDetectorName(bankname);

    // hDQ_Thresholds_[DetName]_[BankName]
    std::string histname = "hDQ_Thresholds_" + detname + "_" + bankname;
    std::string histtitle = "ODB Thresholds in " + bankname;
    TH1F* hDQ_Histogram = new TH1F(histname.c_str(), histtitle.c_str(), 
                                   2,0,2); // will want an upper and lower threshold for FADC channels
    hDQ_Histogram->GetXaxis()->SetTitle("");
    hDQ_Histogram->GetYaxis()->SetTitle("Threshold");
    DQ_Thresholds_histograms_map[bankname] = hDQ_Histogram;
  }

  gDirectory->Cd("/MidasHists/");
  return SUCCESS;
}

---

Variable Documentation

DQ_Thresholds_histograms_map

Map of bank name to histogram. FADC histograms have 2 entries, one each for an upper and lower trigger threshold. CAENs have single entry.

Definition at line 63 of file MDQ_Thresholds.cpp.

Referenced by MDQ_Thresholds_eor(), and MDQ_Thresholds_init().

TGlobalData* gData

Object to hold data used and produced by modules throughout alcapana stage of analysis.

Definition at line 76 of file analyzer.cpp.

TSetupData* gSetup

Hardware information about digitizers and detectors to be used during alcapana stage of analysis.

Definition at line 80 of file analyzer.cpp.

HNDLE hDB

ANA_MODULE MDQ_Thresholds_module

Initial value:

{
        "MDQ_Thresholds",                    
        "Andrew Edmonds",              
        MDQ_Thresholds,                      
        NULL,                          
        MDQ_Thresholds_eor,                          
        MDQ_Thresholds_init,                 
        NULL,                          
        NULL,                          
        0,                             
        NULL,                          
}

Definition at line 65 of file MDQ_Thresholds.cpp.